Given a strictly hyperbolic, genuinely nonlinear system of conservation laws, we prove the a priori bound u (t,) u (t,)= O (1)(1+ t) | ln | on the distance between an exact BV solution u and a viscous approximation u, letting the viscosity coefficient 0. In the proof, starting from u we construct an approximation of the viscous solution u by taking a mollification u* and inserting viscous shock profiles at the locations of finitely many large shocks for each fixed . Error estimates are then obtained by introducing new Lyapunov functionals that control interactions of shock waves in the same family and also interactions of waves in different families. 2004 Wiley Periodicals, Inc.