We prove nonlinear stability of planar shock fronts for certain relaxation systems in two spatial dimensions. If the subcharacteristic condition is assumed and the initial perturbation is sufficiently small and the mass carried by the perturbations is not necessarily finite, then the solution converges to a shifted planar shock front solution as time t ↑ ∞. The asymptotic phase shift of shock fronts is, in general, non-zero and governed by a similarity solution to the heat equation. The asymptotic decay rate to the shock front is proved to be t−1/4 in L∞(R2) without imposing extra decay rates in space for the initial perturbations. The proofs are based on an elementary weighted energy analysis to the error equation.