  Waveform inversion is concerned with estimating a heterogeneous medium,
modeled by variable coefficients of wave equations, using sources that emit
probing signals and receivers that record the generated waves. It is an old and
intensively studied inverse problem with a wide range of applications, but the
existing inversion methodologies are still far from satisfactory. The typical
mathematical formulation is a nonlinear least squares data fit optimization and
the difficulty stems from the non-convexity of the objective function that
displays numerous local minima at which local optimization approaches stagnate.
This pathological behavior has at least three unavoidable causes: (1) The
mapping from the unknown coefficients to the wave field is nonlinear and
complicated. (2) The sources and receivers typically lie on a single side of
the medium, so only backscattered waves are measured. (3) The probing signals
are band limited and with high frequency content. There is a lot of activity in
the computational science and engineering communities that seeks to mitigate
the difficulty of estimating the medium by data fitting. In this paper we
present a different point of view, based on reduced order models (ROMs) of two
operators that control the wave propagation. The ROMs are called data driven
because they are computed directly from the measurements, without any knowledge
of the wave field inside the inaccessible medium. This computation is
non-iterative and uses standard numerical linear algebra methods. The resulting
ROMs capture features of the physics of wave propagation in a complementary way
and have surprisingly good approximation properties that facilitate waveform
inversion.
