Trust-region methods for time-domain full-waveform seismic tomography

Seismic tomography infers the material properties of the Earth based on the observation of seismograms. This can be formulated as a PDE-constrained optimization problem governed by the elastic wave equation. We analyze the inverse problem in a function space setting and establish results on the differentiability of the solution operator. Furthermore, we present a class of trust-region Newton-type methods that can handle additional constraints on the material parameters and work with inexactly computed derivatives.