The adjoint method in geodynamics: derivation from a general operator formulation and application to the initial condition problem in a high resolution mantle circulation model
Horbach, André, HansPeter Bunge, and Jens Oeser (2014),
The adjoint method in geodynamics: derivation from a general operator formulation and application to the initial condition problem in a high resolution mantle circulation model,
GEM  International Journal on Geomathematics, 5(2), 163194, doi:10.1007/s1313701400615.
 Abstract
 The adjoint method is a computationally efficient way to compute the gradient of a physical observable or an associated objective function relative to its parameters. In geodynamics the observable can be thought of as a representation of the present day heterogeneity structure in the Earth’s mantle, inferred in some form through seismic imaging, while a crucial derivative of interest is that relative to an earlier convective system state. Since mantle convection is governed by coupled, nonlinear conservation equations for mass, momentum and energy, computation of the derivative consists of iterative solutions to the forward and the adjoint problem, rendering the approach superior to finite difference approximations, which become impractical at the resolution of modern geodynamic models. Moreover, similarities in the forward and adjoint equations allow one to apply existing numerical codes that solve the forward problem to the adjoint equations with little adaptation.
Bunge et al. (Geophys J Int 152(2):280–301 (2003)), have derived the adjoint equations for mantle convection using the concept of Lagrangian multipliers. Here we introduce a more general approach using an operator formulation in Hilbert spaces, in order to connect to recent work in seismology (Fichtner et al. Phys Earth Planet Int 157(1–2):86–104 (2006a)), where the approach was used to derive the adjoint equations for the scalar wave equation. We demonstrate the practicality of the method for use in a high resolution mantle circulation model with more than 80 million finite elements by restoring a representation of present day mantle heterogeneity derived from the global seismic shear wave study of Grand et al. (GSA Today 7(4):1–7 1997) back in time for the past 40 million years. An important result is our finding of a strong global minimum for the unknown initial condition, regardless of the assumed first guess for the initial heterogeneity structure, which we attribute to the uniqueness theorem by Serrin. Paleo mantle convection modelling will improve our ability to test assumptions about the internal structure and dynamics of the Earth’s mantle against the geologic record.
 Further information

 BibTeX

@article{id2012,
author = {Andr{\'e} Horbach and HansPeter Bunge and Jens Oeser},
journal = {GEM  International Journal on Geomathematics},
number = {2},
pages = {163194},
title = {{The adjoint method in geodynamics: derivation from a general operator formulation and application to the initial condition problem in a high resolution mantle circulation model}},
volume = {5},
year = {2014},
language = {en},
url = {http://dx.doi.org/10.1007/s1313701400615},
doi = {10.1007/s1313701400615},
}
 EndNote

%0 Journal Article
%A Horbach, André
%A Bunge, HansPeter
%A Oeser, Jens
%D 2014
%N 2
%V 5
%J GEM  International Journal on Geomathematics
%P 163194
%T The adjoint method in geodynamics: derivation from a general operator formulation and application to the initial condition problem in a high resolution mantle circulation model
%U http://dx.doi.org/10.1007/s1313701400615