Towards adjoint-based inversion of time-dependent mantle convection with nonlinear viscosity

We develop and study an adjoint-based inversion method for the simultaneous recovery of initial temperature conditions and viscosity parameters in time-dependent mantle convection from the current mantle temperature and historic plate motion. Based on a realistic rheological model with temperature-dependent and strain-rate-dependent viscosity, we formulate the inversion as a PDE-constrained optimization problem. The objective functional includes the misfit of surface velocity (plate motion) history, the misfit of the current mantle temperature, and a regularization for the uncertain initial condition. The gradient of this functional with respect to the initial temperature and the uncertain viscosity parameters is computed by solving the adjoint of the mantle convection equations. This gradient is used in a pre-conditioned quasi-Newton minimization algorithm. We study the prospects and limitations of the inversion, as well as the computational performance of the method using two synthetic problems, a sinking cylinder and a realistic subduction model. The subduction model is characterized by the migration of a ridge toward a trench whereby both plate motions and subduction evolve. The results demonstrate: (1) for known viscosity parameters, the initial temperature can be well recovered, as in previous initial condition-only inversions where the effective viscosity was given; (2) for known initial temperature, viscosity parameters can be recovered accurately, despite the existence of trade-offs due to ill-conditioning; (3) for the joint inversion of initial condition and viscosity parameters, initial condition and effective viscosity can be reasonably recovered, but the high dimension of the parameter space and the resulting ill-posedness may limit recovery of viscosity parameters.

[1]  M. Gurnis,et al.  Assimilating lithosphere and slab history in 4-D Earth models , 2019 .

[2]  M. Gurnis,et al.  Dynamic topography, gravity and the role of lateral viscosity variations from inversion of global mantle flow , 2016 .

[3]  Georg Stadler,et al.  Adjoint-based estimation of plate coupling in a non-linear mantle flow model: Theory and examples , 2015 .

[4]  B. Kaus,et al.  Geodynamic inversion to constrain the non-linear rheology of the lithosphere , 2015 .

[5]  Georg Stadler,et al.  Towards adjoint-based inversion for rheological parameters in nonlinear viscous mantle flow , 2014 .

[6]  H. Bunge,et al.  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 , 2014 .

[7]  S. Kramer,et al.  Interaction of subducted slabs with the mantle transition‐zone: A regime diagram from 2‐D thermo‐mechanical models with a mobile trench and an overriding plate , 2014 .

[8]  Jean-Luc Guermond,et al.  Stability analysis of explicit entropy viscosity methods for non-linear scalar conservation equations , 2013, Math. Comput..

[9]  Georg Stadler,et al.  Large-scale adaptive mantle convection simulation , 2013 .

[10]  Martin Kronbichler,et al.  High accuracy mantle convection simulation through modern numerical methods , 2012 .

[11]  Maria Seton,et al.  Global continental and ocean basin reconstructions since 200 Ma , 2012 .

[12]  Taras Gerya,et al.  A free plate surface and weak oceanic crust produce single‐sided subduction on Earth , 2012 .

[13]  Alfio Borzì,et al.  Computational Optimization of Systems Governed by Partial Differential Equations , 2012, Computational science and engineering.

[14]  T. Gerya Future directions in subduction modeling , 2011 .

[15]  R. Moucha,et al.  Supplementary information for “Changes in African topography driven by mantle convection , 2011 .

[16]  Georg Stadler,et al.  The Dynamics of Plate Tectonics and Mantle Flow: From Local to Global Scales , 2010, Science.

[17]  E. Burkett,et al.  Dynamics and implications of slab detachment due to ridge‐trench collision , 2009 .

[18]  M. Gurnis,et al.  Adjoint models of mantle convection with seismic, plate motion, and stratigraphic constraints: North America since the Late Cretaceous , 2009 .

[19]  Lijun Liu,et al.  Reconstructing Farallon Plate Subduction Beneath North America Back to the Late Cretaceous , 2008, Science.

[20]  Lijun Liu,et al.  Simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection , 2008 .

[21]  Magali I. Billen,et al.  Modeling the Dynamics of Subducting Slabs , 2008 .

[22]  Thomas R. Anderson,et al.  A dynamic model of oceanic sulfur (DMOS) applied to the Sargasso Sea: Simulating the dimethylsulfide (DMS) summer paradox , 2008 .

[23]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[24]  G. Hirth,et al.  Rheologic controls on slab dynamics , 2007 .

[25]  J. Mitrovica,et al.  Lateral variations in mantle rheology: implications for convection related surface observables and inferred viscosity models , 2007 .

[26]  G. Schubert,et al.  Inverse problem of thermal convection: numerical approach and application to mantle plume restoration , 2004 .

[27]  Michael Gurnis,et al.  Evolving force balance during incipient subduction , 2004 .

[28]  Richard G. Derwent,et al.  Radiative forcing in the 21st century due to ozone changes in the troposphere and the lower stratosphere , 2003 .

[29]  M. Gurnis,et al.  Seismic tomography, surface uplift, and the breakup of Gondwanaland: Integrating mantle convection backwards in time , 2003 .

[30]  H. Muhlhaus,et al.  Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics , 2002 .

[31]  C. R. Hagelberg,et al.  Mantle circulation models with variational data assimilation: inferring past mantle flow and structure from plate motion histories and seismic tomography , 2001 .

[32]  Niles A. Pierce,et al.  An Introduction to the Adjoint Approach to Design , 2000 .

[33]  Louis Moresi,et al.  Plate tectonics and convection in the Earth's mantle: toward a numerical simulation , 2000, Comput. Sci. Eng..

[34]  Ulrich R. Christensen,et al.  Dynamic Earth: Plates, Plumes and Mantle Convection , 2000 .

[35]  Richards,et al.  Time scales and heterogeneous structure in geodynamic earth models , 1998, Science.

[36]  Shijie Zhong,et al.  Towards a realistic simulation of plate margins in mantle convection , 1995 .

[37]  M. Gurnis,et al.  Mantle Convection with Plates and Mobile, Faulted Plate Margins , 1995, Science.

[38]  G. Schubert,et al.  A benchmark comparison of numerical methods for infinite Prandtl number thermal convection in two-dimensional Cartesian geometry , 1990 .

[39]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[40]  Manfred Koch,et al.  A benchmark comparison for mantle convection codes , 1989 .

[41]  N. Weiss,et al.  Convection in the earth's mantle: towards a numerical simulation , 1974, Journal of Fluid Mechanics.

[42]  T. Atwater Implications of Plate Tectonics for the Cenozoic Tectonic Evolution of Western North America , 1970 .

[43]  Carsten Burstedde,et al.  p 4 est : Scalable Algorithms for Parallel Adaptive Mesh Refinement , 2018 .

[44]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[45]  Max Gunzburger,et al.  Perspectives in flow control and optimization , 1987 .

[46]  Giorgio Ranalli,et al.  Rheology of the earth , 1987 .