  Conventional full-waveform inversion (FWI) using the least-squares norm
($L^2$) as a misfit function is known to suffer from cycle skipping. This
increases the risk of computing a local rather than the global minimum of the
misfit. In our previous work, we proposed the quadratic Wasserstein metric
($W_2$) as a new misfit function for FWI. The $W_2$ metric has been proved to
have many ideal properties with regards to convexity and insensitivity to
noise. When the observed and predicted seismic data are regarded as two density
functions, the quadratic Wasserstein metric corresponds to the optimal cost of
rearranging one density into the other, where the transportation cost is
quadratic in distance. The difficulty of transforming seismic signals into
nonnegative density functions is discussed. Unlike the $L^2$ norm, $W_2$
measures not only amplitude differences, but also global phase shifts, which
helps to avoid cycle skipping issues. In this work, we build on our earlier
method to cover more realistic high-resolution applications by embedding the
$W_2$ technique into the framework of the adjoint-state method and applying it
to seismic relevant 2D examples: the Camembert, the Marmousi, and the 2004 BP
models. We propose a new way of using the $W_2$ metric trace-by-trace in FWI
and compare it to global $W_2$ via the solution of the Monge-Amp\`ere equation.
With corresponding adjoint source, the velocity model can be updated using the
l-BFGS method. Numerical results show the effectiveness of $W_2$ for
alleviating cycle skipping issues and sensitivity to noise. Both mathematical
theory and numerical examples demonstrate that the quadratic Wasserstein metric
is a good candidate for a misfit function in seismic inversion.
