Regularization techniques for the trace and the traction of elastic waves potentials previously built for domains of the class C^2 are extended to the Lipschitz case. In particular, this yields an elementary way to establish the mapping properties of elastic wave potentials from those of the scalar Helmholtz equation without resorting to the more advanced theory for elliptic systems in the Lipschitz domains. Scalar Gunter derivatives of a function defined on the boundary of a three-dimensional domain are expressed as components (or their opposites) of the tangential vector rotational of this function in the canonical orthonormal basis of the ambient space. This, in particular, implies that these derivatives define bounded operators from H^s to H^{s-1} (0 ≤ s ≤ 1) on the boundary of the Lipschitz domain and can easily be implemented in boundary element codes. Representations of the Gunter operator and potentials of single and double layers of elastic waves in the two-dimensional case are provided.