Sensitivity kernels for viscoelastic loading based on adjoint methods

• Published in: Geophysical Journal International
• Main Authors: Al-Attar, David, Tromp, Jeroen
• Format: Text
• Language: English
• Published: Oxford University Press 2014
• Subjects: Geodynamics and tectonics; Eta; Lagrange; Ice Sheet
• Online Access: http://gji.oxfordjournals.org/cgi/content/short/196/1/34; https://doi.org/10.1093/gji/ggt395

Observations of glacial isostatic adjustment (GIA) allow for inferences to be made about mantle viscosity, ice sheet history and other related parameters. Typically, this inverse problem can be formulated as minimizing the misfit between the given observations and a corresponding set of synthetic data. When the number of parameters is large, solution of such optimization problems can be computationally challenging. A practical, albeit non-ideal, solution is to use gradient-based optimization. Although the gradient of the misfit required in such methods could be calculated approximately using finite differences, the necessary computation time grows linearly with the number of model parameters, and so this is often infeasible. A far better approach is to apply the ‘adjoint method’, which allows the exact gradient to be calculated from a single solution of the forward problem, along with one solution of the associated adjoint problem. As a first step towards applying the adjoint method to the GIA inverse problem, we consider its application to a simpler viscoelastic loading problem in which gravitationally self-consistent ocean loading is neglected. The earth model considered is non-rotating, self-gravitating, compressible, hydrostatically pre-stressed, laterally heterogeneous and possesses a Maxwell solid rheology. We determine adjoint equations and Fréchet kernels for this problem based on a Lagrange multiplier method. Given an objective functional J defined in terms of the surface deformation fields, we show that its first-order perturbation can be written <f>$\delta J = \int _{M_{\mathrm{S}}}K_{\eta }\delta \ln \eta \,\mathrm{d}V +\int _{t_{0}}^{t_{1}}\int _{\partial M}K_{\dot{\sigma }} \delta \dot{\sigma } \,\mathrm{d}S \,\mathrm{d}t$</f>, where δ ln η = δη/η denotes relative viscosity variations in solid regions M S , d V is the volume element, <f>$\delta \dot{\sigma }$</f> is the perturbation to the time derivative of the surface load which is defined on the earth model's surface ∂ M ...