Diese Seite ist aus Gründen der Barrierefreiheit optimiert für aktuelle Browser. Sollten Sie einen älteren Browser verwenden, kann es zu Einschränkungen der Darstellung und Benutzbarkeit der Website kommen!
Geophysics Homepage
Search:
Log in
print

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é, Hans-Peter 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), 163-194, doi:10.1007/s13137-014-0061-5.

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, non-linear 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 Hans-Peter Bunge and Jens Oeser},
  journal = {GEM - International Journal on Geomathematics},
  number = {2},
  pages = {163-194},
  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/s13137-014-0061-5},
  doi = {10.1007/s13137-014-0061-5},
}
EndNote
%0 Journal Article
%A Horbach, André
%A Bunge, Hans-Peter
%A Oeser, Jens
%D 2014
%N 2
%V 5
%J GEM - International Journal on Geomathematics
%P 163-194
%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/s13137-014-0061-5
ImprintPrivacy PolicyContact
Printed 17. Aug 2019 12:34