Articles | Volume 17, issue 13
https://doi.org/10.5194/gmd-17-5057-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Special issue:
https://doi.org/10.5194/gmd-17-5057-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Automatic adjoint-based inversion schemes for geodynamics: reconstructing the evolution of Earth's mantle in space and time
Sia Ghelichkhan
CORRESPONDING AUTHOR
Research School of Earth Sciences, The Australian National University, Canberra, ACT, Australia
Institute for Water Futures, The Australian National University, Canberra, ACT, Australia
Angus Gibson
Research School of Earth Sciences, The Australian National University, Canberra, ACT, Australia
D. Rhodri Davies
Research School of Earth Sciences, The Australian National University, Canberra, ACT, Australia
Stephan C. Kramer
Department of Earth Science and Engineering, Imperial College London, London, UK
David A. Ham
Department of Mathematics, Imperial College London, London, UK
Related authors
William Scott, Mark Hoggard, Thomas Duvernay, Sia Ghelichkhan, Angus Gibson, Dale Roberts, Stephan C. Kramer, and D. Rhodri Davies
EGUsphere, https://doi.org/10.5194/egusphere-2025-4168, https://doi.org/10.5194/egusphere-2025-4168, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Melting ice sheets drive solid Earth deformation and sea-level change on timescales of decades to thousands of years. Here, we present G-ADOPT, which models movement of the solid Earth in response to surface loads. It has flexibility in domain geometry, deformation mechanism parameterisation, and is scalable on high performance computers. Automatic derivation of adjoint sensitivity kernels also provides a means to assimilate historical and modern observations into future sea-level forecasts.
Conor P. B. O'Malley, Gareth G. Roberts, James Panton, Fred D. Richards, J. Huw Davies, Victoria M. Fernandes, and Sia Ghelichkhan
Geosci. Model Dev., 17, 9023–9049, https://doi.org/10.5194/gmd-17-9023-2024, https://doi.org/10.5194/gmd-17-9023-2024, 2024
Short summary
Short summary
We wish to understand how the history of flowing rock within Earth's interior impacts deflection of its surface. Observations exist to address this problem, and mathematics and different computing tools can be used to predict histories of flow. We explore how modeling choices impact calculated vertical deflections. The sensitivity of vertical motions at Earth's surface to deep flow is assessed, demonstrating how surface observations can enlighten flow histories.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
William Scott, Mark Hoggard, Thomas Duvernay, Sia Ghelichkhan, Angus Gibson, Dale Roberts, Stephan C. Kramer, and D. Rhodri Davies
EGUsphere, https://doi.org/10.5194/egusphere-2025-4168, https://doi.org/10.5194/egusphere-2025-4168, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Melting ice sheets drive solid Earth deformation and sea-level change on timescales of decades to thousands of years. Here, we present G-ADOPT, which models movement of the solid Earth in response to surface loads. It has flexibility in domain geometry, deformation mechanism parameterisation, and is scalable on high performance computers. Automatic derivation of adjoint sensitivity kernels also provides a means to assimilate historical and modern observations into future sea-level forecasts.
Davor Dundovic, Joseph G. Wallwork, Stephan C. Kramer, Fabien Gillet-Chaulet, Regine Hock, and Matthew D. Piggott
Geosci. Model Dev., 18, 4023–4044, https://doi.org/10.5194/gmd-18-4023-2025, https://doi.org/10.5194/gmd-18-4023-2025, 2025
Short summary
Short summary
Accurate numerical studies of glaciers often require high-resolution simulations, which often prove too demanding even for modern computers. In this paper we develop a method that identifies whether different parts of a glacier require high or low resolution based on its physical features, such as its thickness and velocity. We show that by doing so we can achieve a more optimal simulation accuracy for the available computing resources compared to uniform-resolution simulations.
Conor P. B. O'Malley, Gareth G. Roberts, James Panton, Fred D. Richards, J. Huw Davies, Victoria M. Fernandes, and Sia Ghelichkhan
Geosci. Model Dev., 17, 9023–9049, https://doi.org/10.5194/gmd-17-9023-2024, https://doi.org/10.5194/gmd-17-9023-2024, 2024
Short summary
Short summary
We wish to understand how the history of flowing rock within Earth's interior impacts deflection of its surface. Observations exist to address this problem, and mathematics and different computing tools can be used to predict histories of flow. We explore how modeling choices impact calculated vertical deflections. The sensitivity of vertical motions at Earth's surface to deep flow is assessed, demonstrating how surface observations can enlighten flow histories.
Reuben W. Nixon-Hill, Daniel Shapero, Colin J. Cotter, and David A. Ham
Geosci. Model Dev., 17, 5369–5386, https://doi.org/10.5194/gmd-17-5369-2024, https://doi.org/10.5194/gmd-17-5369-2024, 2024
Short summary
Short summary
Scientists often use models to study complex processes, like the movement of ice sheets, and compare them to measurements for estimating quantities that are hard to measure. We highlight an approach that ensures accurate results from point data sources (e.g. height measurements) by evaluating the numerical solution at true point locations. This method improves accuracy, aids communication between scientists, and is well-suited for integration with specialised software that automates processes.
D. Rhodri Davies, Stephan C. Kramer, Sia Ghelichkhan, and Angus Gibson
Geosci. Model Dev., 15, 5127–5166, https://doi.org/10.5194/gmd-15-5127-2022, https://doi.org/10.5194/gmd-15-5127-2022, 2022
Short summary
Short summary
Firedrake is a state-of-the-art system that automatically generates highly optimised code for simulating finite-element (FE) problems in geophysical fluid dynamics. It creates a separation of concerns between employing the FE method and implementing it. Here, we demonstrate the applicability and benefits of Firedrake for simulating geodynamical flows, with a focus on the slow creeping motion of Earth's mantle over geological timescales, which is ultimately the engine driving our dynamic Earth.
Stephan C. Kramer, D. Rhodri Davies, and Cian R. Wilson
Geosci. Model Dev., 14, 1899–1919, https://doi.org/10.5194/gmd-14-1899-2021, https://doi.org/10.5194/gmd-14-1899-2021, 2021
Short summary
Short summary
Computational models of Earth's mantle require rigorous verification and validation. Analytical solutions of the underlying Stokes equations provide a method to verify that these equations are accurately solved for. However, their derivation in spherical and cylindrical shell domains with physically relevant boundary conditions is involved. This paper provides a number of solutions. They are provided in a Python package (Assess) and their use is demonstrated in a convergence study with Fluidity.
Lior Suchoy, Saskia Goes, Benjamin Maunder, Fanny Garel, and Rhodri Davies
Solid Earth, 12, 79–93, https://doi.org/10.5194/se-12-79-2021, https://doi.org/10.5194/se-12-79-2021, 2021
Short summary
Short summary
We use 2D numerical models to highlight the role of basal drag in subduction force balance. We show that basal drag can significantly affect velocities and evolution in our simulations and suggest an explanation as to why there are no trends in plate velocities with age in the Cenozoic subduction record (which we extracted from recent reconstruction using GPlates). The insights into the role of basal drag will help set up global models of plate dynamics or specific regional subduction models.
Cited articles
Abadi, M., Agarwal, A., Barham, P., Brevdo, E., Chen, Z., Citro, C., Corrado, G. S., Davis, A., Dean, J., Devin, M., Ghemawat, S., Goodfellow, I., Harp, A., Irving, G., Isard, M., Jia, Y., Jozefowicz, R., Kaiser, L., Kudlur, M., Levenberg, J., Mané, D., Monga, R., Moore, S., Murray, D., Olah, C., Schuster, M., Shlens, J., Steiner, B., Sutskever, I., Talwar, K., Tucker, P., Vanhoucke, V., Vasudevan, V., Viégas, F., Vinyals, O., Warden, P., Wattenberg, M., Wicke, M., Yu, Y., and Zheng, X.: TensorFlow: Large-Scale Machine Learning on Heterogeneous Systems, https://www.tensorflow.org/ (last access: 30 June 2024), 2015. a
Al-Attar, D. and Tromp, J.: Sensitivity kernels for viscoelastic loading based on adjoint methods, Geophys. J. Int., 196, 34–77, 2014. a
Alnæs, M. S., Logg, A., Ølgaard, K. B., Rognes, M. E., and Wells, G. N.: Unified form language: A domain-specific language for weak formulations of partial differential equations, ACM T. Math. Software, 40, 1–37, https://doi.org/10.1145/2566630, 2014. a, b
Balay, S., Gropp, W. D., McInnes, L. C., and Smith, B. F.: Efficient Management of Parallelism in Object Oriented Numerical Software Libraries, in: Modern Software Tools in Scientific Computing, edited by: Arge, E., Bruaset, A. M., and Langtangen, H. P., 163–202, Birkhäuser Press, 1997. a
Balay, S., Abhyankar, S., Adams, M. F., Benson, S., Brown, J., Brune, P., Buschelman, K., Constantinescu, E. M., Dalcin, L., Dener, A., Eijkhout, V., Faibussowitsch, J., Gropp, W. D., Hapla, V., Isaac, T., Jolivet, P., Karpeev, D., Kaushik, D., Knepley, M. G., Kong, F., Kruger, S., May, D. A., McInnes, L. C., Mills, R. T., Mitchell, L., Munson, T., Roman, J. E., Rupp, K., Sanan, P., Sarich, J., Smith, B. F., Zampini, S., Zhang, H., Zhang, H., and Zhang, J.: PETSc Web page, https://petsc.org/ (last access: 30 June 2024), 2023. a
Ball, P., Czarnota, K., White, N., Klöcking, M., and Davies, D. R.: Thermal structure of eastern Australia's upper mantle and its relationship to Cenozoic volcanic activity and dynamic topography, Geochem. Geophys. Geosyst., 22, e2021GC009717, https://doi.org/10.1029/2021gc009717, 2021. a
Bauer, S., Huber, M., Ghelichkhan, S., Mohr, M., Rüde, U., and Wohlmuth, B.: Large-scale simulation of mantle convection based on a new matrix-free approach, J. Comput. Sci., 31, 60–76, 2019. a
Bauer, S., Bunge, H.-P., Drzisga, D., Ghelichkhan, S., Huber, M., Kohl, N., Mohr, M., Rüde, U., Thönnes, D., and Wohlmuth, B. I.: TerraNeo–Mantle Convection Beyond a Trillion Degrees of Freedom, in: Software for Exascale Computing – SPPEXA 2016–2019, Lecture Notes in Computational Science and Engineering, edited by: Bungartz, H J., Reiz, S., Uekermann, B., Neumann, P., and Nagel, W., vol 136, Springer, Cham. https://doi.org/10.1007/978-3-030-47956-5_19, 2020. a
Bocher, M., Coltice, N., Fournier, A., and Tackley, P. J.: A sequential data assimilation approach for the joint reconstruction of mantle convection and surface tectonics, Geophys. J. Int., 204, 200–214, https://doi.org/10.1093/gji/ggv427, 2016. a
Bocher, M., Fournier, A., and Coltice, N.: Ensemble Kalman filter for the reconstruction of the Earth's mantle circulation, Nonlin. Processes Geophys., 25, 99–123, https://doi.org/10.5194/npg-25-99-2018, 2018. a, b
Bower, D. J., Gurnis, M., and Seton, M.: Lower mantle structure from paleogeographically constrained dynamic Earth models, Geochem. Geophys. Geosyst., 14, 44–63, https://doi.org/10.1029/2012GC004267, 2013. a
Bozdağ, E., Peter, D., Lefebvre, M., Komatitsch, D., Tromp, J., Hill, J., Podhorszki, N., and Pugmire, D.: Global adjoint tomography: first-generation model, Geophys. J. Int., 207, 1739–1766, https://doi.org/10.1093/gji/ggw356, 2016. a
Bunge, H.-P. and Grand, S. P.: Mesozoic plate-motion history below the northeast Pacific Ocean from seismic images of the subducted Farallon slab, Nature, 405, 337–340, 2000. a
Bunge, H.-P., Richards, M., and Baumgardner, J.: Mantle–circulation models with sequential data assimilation: inferring present–day mantle structure from plate–motion histories, Philos. T. Roy. Soc. Lond. A, 360, 2545–2567, https://doi.org/10.1098/rsta.2002.1080, 2002. a
Bunge, H.-P., Hagelberg, C. R., and Travis, B. J.: Mantle circulation models with variational data assimilation: inferring past mantle flow and structure from plate motion histories and seismic tomography, Geophys. J. Int., 152, 280–301, https://doi.org/10.1046/j.1365-246X.2003.01823.x, 2003. a, b, c, d, e, f
Carrassi, A., Bocquet, M., Bertino, L., and Evensen, G.: Data assimilation in the geosciences: An overview of methods, issues, and perspectives, WIREs Clim. Change, 9, e535, https://doi.org/10.1002/wcc.535, 2018. a
Chust, T. C., Steinle-Neumann, G., Dolejš, D., Schuberth, B. S. A., and Bunge, H.-P.: MMA-EoS: A Computational Framework for Mineralogical Thermodynamics, J. Geophys. Res.-Sol. Ea., 122, 9881–9920, https://doi.org/10.1002/2017JB014501, 2017. a
Cloetingh, S. and Haq, B. U.: Inherited landscapes and sea level change, Science, 347, 1258375, https://doi.org/10.1126/science.1258375, 2015. a
Colli, L., Stotz, I., Bunge, H.-P., Smethurst, M., Clark, S. R., Iaffaldano, G., Tassara, A., Guillocheau, F., and Bianchi, M. C.: Rapid South Atlantic spreading changes and coeval vertical motion in surrounding continents: Evidence for temporal changes of pressure-driven upper mantle flow, Tectonics, 33, 1304–1321, https://doi.org/10.1002/2014TC003612, 2014. a
Colli, L., Bunge, H.-P., and Schuberth, B. S. A.: On retrodictions of global mantle flow with assimilated surface velocities, Geophys. Res. Lett., 42, 8341–8348, https://doi.org/10.1002/2015GL066001, 2015. a
Colli, L., Ghelichkhan, S., Bunge, H.-P., and Oeser, J.: Retrodictions of Mid Paleogene mantle flow and dynamic topography in the Atlantic region from compressible high resolution adjoint mantle convection models: Sensitivity to deep mantle viscosity and tomographic input model, Gondwana Res., 53, 252–272, https://doi.org/10.1016/j.gr.2017.04.027, 2018. a, b, c, d, e
Colli, L., Bunge, H., and Oeser, J.: Impact of model inconsistencies on reconstructions of past mantle flow obtained using the adjoint method, Geophys. J. Int., 221, 617–639, https://doi.org/10.1093/gji/ggaa023, 2020. a, b
Coltice, N., Gérault, M., and Ulvrová, M.: A mantle convection perspective on global tectonics, Earth-Sci. Rev., 165, 120–150, https://doi.org/10.1016/j.earscirev.2016.11.006, 2017. a, b
Coltice, N., Blessing, S., Giering, R., and Tackley, P.: Sensitivity Analysis of Global Kinematics on Mantle Structure Using Automatically Generated Adjoint Thermochemical Convection Codes, Earth ArXiv, https://doi.org/10.31223/X5N09Q, 2023. a
Colton, D. L., Kress, R., and Kress, R.: Inverse acoustic and electromagnetic scattering theory, vol. 93, Springer, 1998. a
Crawford, O., Al-Attar, D., Tromp, J., Mitrovica, J. X., Austermann, J., and Lau, H. C.: Quantifying the sensitivity of post-glacial sea level change to laterally varying viscosity, Geophys. J. Int., 214, 1324–1363, https://doi.org/10.1093/gji/ggy184, 2018. a
Dalcin, L. D., Paz, R. R., Kler, P. A., and Cosimo, A.: Parallel distributed computing using Python, Adv. Water Resour., 34, 1124–1139, https://doi.org/10.1016/j.advwatres.2011.04.013, 2011. a
Dannberg, J. and Gassmöller, R.: Chemical trends in ocean islands explained by plume–slab interaction, P. Natl. Acad. Sci. USA, 115, 4351–4356, https://doi.org/10.1073/pnas.1714125115, 2018. a
Davies, D. R. and Davies, J. H.: Thermally-driven mantle plumes reconcile multiple hot-spot observations, Earth Planet. Sc. Lett., 278, 50–54, 2009. a
Davies, D. R., Goes, S., Davies, J. H., Schuberth, B. S. A., Bunge, H.-P., and Ritsema, J.: Reconciling dynamic and seismic models of Earth's lower mantle: The dominant role of thermal heterogeneity, Earth Planet. Sc. Lett., 353–354, 253–269, https://doi.org/10.1016/j.epsl.2012.08.016, 2012. a
Davies, D. R., Goes, S., and Lau, H. C. P.: Thermally Dominated Deep Mantle LLSVPs: A Review, pp. 441–477, Springer International Publishing, Cham, ISBN 978-3-319-15627-9, https://doi.org/10.1007/978-3-319-15627-9_14, 2015. a
Davies, D. R., Ghelichkhan, S., Hoggard, M., Valentine, A. P., and Richards, F. D.: Observations and models of dynamic topography: Current status and future directions, Dynamics of Plate Tectonics and Mantle Convection, 2023, 223–269, https://doi.org/10.1016/B978-0-323-85733-8.00017-2, 2023. a, b, c, d
Davies, G. F. and Richards, M. A.: Mantle Convection, J. Geol., 100, 151–206, https://doi.org/10.1086/629582, 1992. a
Dembo, R. S. and Steihaug, T.: Truncated-Newton algorithms for large-scale unconstrained optimization, Math. Program., 26, 190–212, https://doi.org/10.1007/BF02592055, 1983. a
Farrell, P. E., Ham, D. A., Funke, S. W., and Rognes, M. E.: Automated derivation of the adjoint of high-level transient finite element programs, SIAM J. Sci. Comput., 35, 1–27, https://doi.org/10.1137/120873558, 2013a. a, b, c
Farrell, P. E., Ham, D. A., Funke, S. W., and Rognes, M. E.: Automated derivation of the adjoint of high-level transient finite element programs, SIAM J. Sci. Comput., 35, C369–C393, 2013b. a
Fichtner, A., van Herwaarden, D.-P., Afanasiev, M., Simutė, S., Krischer, L., Çubuk-Sabuncu, Y., Taymaz, T., Colli, L., Saygin, E., Villaseñor, A., Trampert, J., Cupillard, P., Bunge, H. P., and Igel, H.: The collaborative seismic earth model: generation 1, Geophys. Res. Lett., 45, 4007–4016, https://doi.org/10.1029/2018GL077338, 2018. a
firedrake-zenodo: Software used in `Software for g-adopt v2.1.0' (Firedrake_20231027.0), Zenodo [code], https://doi.org/10.5281/zenodo.10047031, 2023. a
French, S. W. and Romanowicz, B. A.: Whole-mantle radially anisotropic shear velocity structure from spectral-element waveform tomography, Geophys. J. Int., 199, 1303–1327, https://doi.org/10.1093/gji/ggu334, 2014. a
Friedrich, A. M., Bunge, H.-P., Rieger, S. M., Colli, L., Ghelichkhan, S., and Nerlich, R.: Stratigraphic framework for the plume mode of mantle convection and the analysis of interregional unconformities on geological maps, Gondwana Res., 53, 159–188, https://doi.org/10.1016/j.gr.2017.06.003, 2018. a
Funke, S. W. and Farrell, P. E.: A framework for automated PDE-constrained optimisation, arXiv [preprint], arXiv:1302.3894, 2013. a
Ghelichkhan, S. and Bunge, H.-P.: The compressible adjoint equations in geodynamics: Derivation and numerical assessment, GEM-International J. Geomath., 7, 1–30, https://doi.org/10.1007/s13137-016-0080-5, 2016. a
Ghelichkhan, S. and Bunge, H.-P.: The adjoint equations for thermochemical compressible mantle convection: derivation and verification by twin experiments, P. Roy. Soc. A, 474, 20180329, https://doi.org/10.1098/rspa.2018.0329, 2018. a
Ghelichkhan, S., Murböck, M., Colli, L., Pail, R., and Bunge, H.-P.: On the observability of epeirogenic movement in current and future gravity missions, Gondwana Res., 53, 273–284, https://doi.org/10.1016/j.gr.2017.04.016, 2018. a
Ghelichkhan, S., Fuentes, J. J., Hoggard, M. J., Richards, F. D., and Mitrovica, J. X.: The precession constant and its long-term variation, Icarus, 358, 114172, https://doi.org/10.1016/j.icarus.2020.114172, 2020. a
Ghelichkhan, S., Bunge, H.-P., and Oeser, J.: Global mantle flow retrodictions for the early Cenozoic using an adjoint method: evolving dynamic topographies, deep mantle structures, flow trajectories and sublithospheric stresses, Geophys. J. Int., 226, 1432–1460, https://doi.org/10.1093/gji/ggab108, 2021. a, b, c, d, e, f, g
Gibson, A., Davies, R., Kramer, S., Ghelichkhan, S., Duvernay, T., and Turner, R.: g-adopt/g-adopt: v2.1.0 (v2.1.0), Zenodo [code], https://doi.org/10.5281/zenodo.10050733, 2023. a
Giering, R. and Kaminski, T.: Recipes for adjoint code construction, ACM T. Math. Software, 24, 437–474, https://doi.org/10.1145/293686.293695, 1998. a
Gunzburger, M.: Adjoint equation-based methods for control problems in incompressible, viscous flows, Flow Turbul. Combust., 65, 249–272, 2000. a
Gunzburger, M. D.: Perspectives in flow control and optimization, Society for Industrial and Applied Mathematics, https://doi.org/10.1137/1.9780898718720, 2002. a, b
Gurnis, M., Turner, M., Zahirovic, S., DiCaprio, L., Spasojevic, S., Müller, R. D., Boyden, J., Seton, M., Manea, V. C., and Bower, D. J.: Plate tectonic reconstructions with continuously closing plates, Comput. Geosci., 38, 35–42, https://doi.org/10.1016/j.cageo.2011.04.014, 2012. a, b, c
Ham, D. A., Kelly, P. H. J., Mitchell, L., Cotter, C. J., Kirby, R. C., Sagiyama, K., Bouziani, N., Vorderwuelbecke, S., Gregory, T. J., Betteridge, J., Shapero, D. R., Nixon-Hill, R. W., Ward, C. J., Farrell, P. E., Brubeck, P. D., Marsden, I., Gibson, T. H., Homolya, M., Sun, T., McRae, A. T. T., Luporini, F., Gregory, A., Lange, M., Funke, S. W., Rathgeber, F., Bercea, G.-T., and Markall, G. R.: Firedrake User Manual, Imperial College London and University of Oxford and Baylor University and University of Washington, 1st edn., https://doi.org/10.25561/104839, 2023. a, b
Hansen, P. C.: Analysis of discrete ill-posed problems by means of the L-curve, SIAM review, 34, 561–580, 1992. a
Hastie, T., Tibshirani, R., Friedman, J. H., and Friedman, J. H.: The elements of statistical learning: data mining, inference, and prediction, vol. 2, Springer, 2009. a
Hinze, M., Pinnau, R., Ulbrich, M., and Ulbrich, S.: Optimization with PDE constraints, vol. 23, Springer Science & Business Media, 2008. a
Hoggard, M., Austerman, J., Randel, C., and Stephenson, S.: Observing dynamic topography through space and time, AGU Geophys. Monogr. 263, Washington, DC, https://doi.org/10.1002/9781119528609.ch15, 2021. a
Hoggard, M. J., Czarnota, K., Richards, F. D., Huston, D. L., Jaques, A. L., and Ghelichkhan, S.: Global distribution of sediment-hosted metals controlled by craton edge stability, Nat. Geosci., 13, 504–510, https://doi.org/10.1038/s41561-020-0593-2, 2020. a
Horbach, A., Bunge, H.-P., and Oeser, J.: 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: Int. J. Geomath., 5, 163–194, https://doi.org/10.1007/s13137-014-0061-5, 2014. a, b
Iglesias, J. A., Sturm, K., and Wechsung, F.: Two-dimensional shape optimization with nearly conformal transformations, SIAM J, Sci, Comput,, 40, A3807–A3830, 2018. a
Jarvis, G. T. and McKenzie, D. P.: Convection in a compressible fluid with infinite Prandtl number, J. Fluid Mech., 96, 515–583, https://doi.org/10.1017/S002211208000225X, 1980. a
Koelemeijer, P., Ritsema, J., Deuss, A., and Van Heijst, H.-J.: SP12RTS: a degree-12 model of shear-and compressional-wave velocity for Earth's mantle, Geophys. J. Int., 204, 1024–1039, https://doi.org/10.1093/gji/ggv481, 2016. a
Kouri, D. P., Jakeman, J. D., Huerta, J. G., Smith, C. B., Walsh, T. F., Udell, M., and Uryasev, S.: Risk-Adaptive Experimental Design for High-Consequence Systems: LDRD Final Report, Tech. rep., Sandia National Lab.(SNL-NM), Albuquerque, NM (United States), 2021a. a
Kouri, D. P., Ridzal, D., and Tuminaro, R.: KKT preconditioners for PDE-constrained optimization with the Helmholtz equation, SIAM J. Sci. Comput., 43, S225–S248, 2021b. a
Kouri, D. P., Staudigl, M., and Surowiec, T. M.: A relaxation-based probabilistic approach for PDE-constrained optimization under uncertainty with pointwise state constraints, Computational Optimization and Applications, 85, 441–478, 2023. a
Kramer, S. C., Wilson, C. R., and Davies, D. R.: An implicit free surface algorithm for geodynamical simulations, Phys. Earth Planet. In., 194, 25–37, 2012. a
Li, D., Gurnis, M., and Stadler, G.: Towards adjoint-based inversion of time-dependent mantle convection with nonlinear viscosity, Geophys. J. Int., 209, 86–105, https://doi.org/10.1093/gji/ggw493, 2017. a, b, c
Lin, C.-J. and Moré, J. J.: Newton's method for large bound-constrained optimization problems, SIAM J. Optim., 9, 1100–1127, https://doi.org/10.1137/S1052623498345075, 1999. a
Liu, L. and Gurnis, M.: Dynamic subsidence and uplift of the Colorado Plateau, Geology, 38, 663–666, https://doi.org/10.1130/g30624.1, 2010. a, b, c
Liu, L., Spasojević, S., and Gurnis, M.: Reconstructing Farallon plate subduction beneath North America back to the Late Cretaceous, Science, 322, 934–938, https://doi.org/10.1126/science.1162921, 2008. a, b
Liu, L., Gurnis, M., Seton, M., Saleeby, J., Müller, R. D., and Jackson, J. M.: The role of oceanic plateau subduction in the Laramide orogeny, Nat. Geosci., 3, 353–357, https://doi.org/10.1038/ngeo829, 2010. a, b
Logg, A., Mardal, K.-A., and Wells, G.: Automated solution of differential equations by the finite element method: The FEniCS book, vol. 84, Springer Science & Business Media, https://doi.org/10.1007/978-3-642-23099-8, 2012. a
Martinec, Z., Sasgen, I., and Velímskỳ, J.: The forward sensitivity and adjoint-state methods of glacial isostatic adjustment, Geophys. J. Int., 200, 77–105, 2015. a
Merdith, A. S., Williams, S. E., Collins, A. S., Tetley, M. G., Mulder, J. A., Blades, M. L., Young, A., Armistead, S. E., Cannon, J., Zahirovic, S., and Müller, R. D.: Extending full-plate tectonic models into deep time: Linking the Neoproterozoic and the Phanerozoic, Earth-Sci. Rev., 214, 103477, https://doi.org/10.1016/j.earscirev.2020.103477, 2021. a, b
Mitusch, S. K., Funke, S. W., and Dokken, J. S.: dolfin-adjoint 2018.1: automated adjoints for FEniCS and Firedrake, J. Open Source Softw., 4, 1292, https://doi.org/10.21105/joss.01292, 2019. a, b, c
Mosca, I., Cobden, L., Deuss, A., Ritsema, J., and Trampert, J.: Seismic and mineralogical structures of the lower mantle from probabilistic tomography, J. Geophys. Res., 117, B06304, https://doi.org/10.1029/2011jb008851, 2012. a
Moses, W. S., Narayanan, S. H. K., Paehler, L., Churavy, V., Schanen, M., Hückelheim, J., Doerfert, J., and Hovland, P.: Scalable automatic differentiation of multiple parallel paradigms through compiler augmentation, in: Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis, SC '22, IEEE Press, ISBN 9784665454445, 2022. a
Moucha, R., Forte, A. M., Mitrovica, J. X., Rowley, D. B., Quéré, S., Simmons, N. A., and Grand, S. P.: Dynamic topography and long-term sea-level variations: There is no such thing as a stable continental platform, Earth Planet. Sc. Lett., 271, 101–108, 2008. a
Müller, R. D., Seton, M., Zahirovic, S., Williams, S. E., Matthews, K. J., Wright, N. M., Shephard, G. E., Maloney, K. T., Barnett-Moore, N., Hosseinpour, M., Bower, D. J., and Cannon, J.: Ocean Basin Evolution and Global-Scale Plate Reorganization Events Since Pangea Breakup, Annu. Rev. Earth Planet. Sci., 44, 107–138, https://doi.org/10.1146/annurev-earth-060115-012211, 2016. a
Müller, R. D., Cannon, J., Qin, X., Watson, R. J., Gurnis, M., Williams, S., Pfaffelmoser, T., Seton, M., Russell, S. H. J., and Zahirovic, S.: GPlates: Building a Virtual Earth Through Deep Time, Geochem. Geophys. Geosyst., 19, 2243–2261, https://doi.org/10.1029/2018GC007584, 2018. a
Müller, R. D., Zahirovic, S., Williams, S. E., Cannon, J., Seton, M., Bower, D. J., Tetley, M. G., Heine, C., Le Breton, E., Liu, S., Liu, S., Russell, S. H. J., Yang, T., Leonard, J., and Gurnis, M.: A global plate model including lithospheric deformation along major rifts and orogens since the Triassic, Tectonics, 38, 1884–1907, https://doi.org/10.1029/2018TC005462, 2019. a, b
Müller, R. D., Flament, N., Cannon, J., Tetley, M. G., Williams, S. E., Cao, X., Bodur, Ö. F., Zahirovic, S., and Merdith, A.: A tectonic-rules-based mantle reference frame since 1 billion years ago – implications for supercontinent cycles and plate–mantle system evolution, Solid Earth, 13, 1127–1159, https://doi.org/10.5194/se-13-1127-2022, 2022. a
Müller, R. D., Cannon, J., Qin, X., Watson, R. J., Gurnis, M., Williams, S., Pfaffelmoser, T., Seton, M., Russell, S. H. J., and Zahirovic, S.: GPlates: Building a Virtual Earth Through Deep Time, Geochem. Geophys. Geosyst., 19, 2243–2261, https://doi.org/10.1029/2018GC007584, 2018. a
Naumann, U.: The Art of Differentiating Computer Programs, Society for Industrial and Applied Mathematics, https://doi.org/10.1137/1.9781611972078, 2011. a, b, c
Nerlich, R., Colli, L., Ghelichkhan, S., Schuberth, B., and Bunge, H.-P.: Constraining central Neo-Tethys Ocean reconstructions with mantle convection models, Geophys. Res. Lett., 43, 9595–9603, https://doi.org/10.1002/2016gl070524, 2016. a
Nitsche, J.: Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind, in: Abhandlungen aus dem mathematischen Seminar der Universität Hamburg, vol. 36, 9–15, Springer, 1971. a
Nixon-Hill, R. W., Shapero, D., Cotter, C. J., and Ham, D. A.: Consistent Point Data Assimilation in Firedrake and Icepack, arXiv [preprint], https://doi.org/10.48550/arXiv.2304.06058, 2023. a
Nocedal, J. and Wright, S. J. (Eds.): Numerical optimization, Springer, Springer New York, NY, https://doi.org/10.1007/b98874, , 1999. a
Nocedal, J. and Wright, S. J.: Quadratic programming, Numerical optimization, 448–492, 2006. a
Panasyuk, S. V. and Hager, B. H.: Inversion for mantle viscosity profiles constrained by dynamic topography and the geoid, and their estimated errors, Geophys. J. Int., 143, 821–836, 2000. a
Panton, J., Davies, J. H., and Myhill, R.: The Stability of Dense Oceanic Crust Near the Core-Mantle Boundary, J. Geophys. Res.-Sol. Ea., 128, e2022JB025610, https://doi.org/10.1029/2022JB025610, 2023. a
Paszke, A., Gross, S., Massa, F., Lerer, A., Bradbury, J., Chanan, G., Killeen, T., Lin, Z., Gimelshein, N., Antiga, L., Desmaison, A., Kopf, A., Yang, E., DeVito, Z., Raison, M., Tejani, A., Chilamkurthy, S., Steiner, B., Fang, L., Bai, J., and Chintala, S.: PyTorch: An Imperative Style, High-Performance Deep Learning Library, in: Advances in Neural Information Processing Systems 32, 8024–8035, Curran Associates, Inc., http://papers.neurips.cc/paper/9015-pytorch-an-imperative-style-high-performance-deep-learning-library.pdf (last access: 30 June 2024), 2019. a
Plessix, R.-E.: A review of the adjoint-state method for computing the gradient of a functional with geophysical applications, Geophys. J. Int., 167, 495–503, 2006. a
Poore, H. R., Samworth, R., White, N. J., Jones, S. M., and McCave, I. N.: Neogene overflow of Northern Component Water at the Greenland-Scotland Ridge, Geochem. Geophys. Geosyst., 7, Q06010, https://doi.org/10.1029/2005GC001085, 2006. a
Price, M. G. and Davies, J. H.: Profiling the robustness, efficiency and limits of the forward-adjoint method for 3D mantle convection modelling, Geophys. J. Int., 212, 1450–1462, https://doi.org/10.1093/gji/ggx489, 2018. a
Rathgeber, F., Ham, D. A., Mitchell, L., Lange, M., Luporini, F., Mcrae, A. T. T., Bercea, G.-T., Markall, G. R., and Kelly, P. H. J.: Firedrake: Automating the Finite Element Method by Composing Abstractions, ACT Trans. Math. Softw., 43, 1–24, https://doi.org/10.1145/2998441, 2016. a, b, c
Ratnaswamy, V., Stadler, G., and Gurnis, M.: Adjoint-based estimation of plate coupling in a non-linear mantle flow model: theory and examples, Geophys. J. Int., 202, 768–786, 2015. a
Rawlinson, N., Pozgay, S., and Fishwick, S.: Seismic tomography: a window into deep Earth, Phys. Earth Planet. In., 178, 101–135, https://doi.org/10.1016/j.pepi.2009.10.002, 2010. a
Reuber, G. S., Holbach, L., Popov, A. A., Hanke, M., and Kaus, B. J.: Inferring rheology and geometry of subsurface structures by adjoint-based inversion of principal stress directions, Geophys. J. Int., 223, 851–861, https://doi.org/10.1093/gji/ggaa344, 2020. a, b
Ricard, Y.: Physics of mantle convection, Treatise on geophysics, 7, 31–87, 2007. a
Sabadini, R. and Vermeersen, L.: Influence of lithospheric and mantle stratification on global post-seismic deformation, Geophys. Res. Lett., 24, 2075–2078, 1997. a
Sambridge, M. and Mosegaard, K.: Monte Carlo methods in geophysical inverse problems, Rev. Geophys., 40, 3–1, https://doi.org/10.1029/2000rg000089, 2002. a
Schwedes, T., Funke, S. W., and Ham, D. A.: An iteration count estimate for a mesh-dependent steepest descent method based on finite elements and Riesz inner product representation, arXiv [preprint], arXiv:1606.08069, 2016. a
Schwedes, T., Ham, D. A., Funke, S. W., Piggott, M. D., Schwedes, T., Ham, D. A., Funke, S. W., and Piggott, M. D.: Mesh dependence in PDE-constrained optimisation, Springer, https://doi.org/10.1007/978-3-319-59483-5_2, 2017. a
Seton, M., Müller, R. D., Zahirovic, S., Gaina, C., Torsvik, T. H., Shephard, G., Talsma, A., Gurnis, M., Turner, M., Maus, S., and Chandler, M.: Global continental and ocean basin reconstructions since 200 Ma, Earth-Sci. Rev., 113, 212–270, https://doi.org/10.1016/j.earscirev.2012.03.002, 2012. a
Seton, M., Müller, R. D., Zahirovic, S., Williams, S., Wright, N. M., Cannon, J., Whittaker, J. M., Matthews, K. J., and McGirr, R.: A Global Data Set of Present-Day Oceanic Crustal Age and Seafloor Spreading Parameters, Geochem. Geophys. Geosyst., 21, e2020GC009214, https://doi.org/10.1029/2020GC009214, 2020. a
Shephard, G., Bunge, H.-P., Schuberth, B., Müller, R., Talsma, A., Moder, C., and Landgrebe, T.: Testing absolute plate reference frames and the implications for the generation of geodynamic mantle heterogeneity structure, Earth Planet. Sc. Lett., 317–318, 204–217, https://doi.org/10.1016/j.epsl.2011.11.027, 2012. a
Simmons, N. A., Myers, S. C., Johannesson, G., Matzel, E., and Grand, S. P.: Evidence for long-lived subduction of an ancient tectonic plate beneath the southern Indian Ocean, Geophys. Res. Lett., 42, 9270–9278, https://doi.org/10.1002/2015gl066237, 2015. a
Stadler, G., Gurnis, M., Burstedde, C., Wilcox, L. C., Alisic, L., and Ghattas, O.: The Dynamics of Plate Tectonics and Mantle Flow: From Local to Global Scales, Science, 329, 1033–1038, https://doi.org/10.1126/science.1191223, 2010. a
Stixrude, L. and Lithgow-Bertelloni, C.: Thermodynamics of mantle minerals - II. Phase equilibria, Geophys. J. Int., 184, 1180–1213, https://doi.org/10.1111/j.1365-246X.2010.04890.x, 2011. a
Stotz, I., Iaffaldano, G., and Davies, D. R.: Pressure-driven Poiseuille flow: A major component of the torque-balance governing Pacific plate motion, Geophys. Res. Lett., 45, 117–125, 2018. a
Stotz, I. L., Iaffaldano, G., and Davies, D. R.: Late Miocene Pacific plate kinematic change explained with coupled global models of mantle and lithosphere dynamics, Geophys. Res. Lett., 44, 7177–7186, https://doi.org/10.1002/2017GL073920, 2017. a
Styles, E., Davies, D. R., and Goes, S.: Mapping spherical seismic into physical structure: biases from 3-D phase-transition and thermal boundary-layer heterogeneity, Geophys. J. Int., 184, 1371–1378, https://doi.org/10.1111/j.1365-246x.2010.04914.x, 2011. a
Taiwo, A., Bunge, H.-P., Schuberth, B. S. A., Colli, L., and Vilacis, B.: Robust global mantle flow trajectories and their validation via dynamic topography histories, Geophys. J. Int., 234, 2160–2179, https://doi.org/10.1093/gji/ggad188, 2023. a
Talagrand, O.: Assimilation of observations, an introduction (gtspecial issueltdata assimilation in meteology and oceanography: Theory and practice), J. Meteorol. Soc. Jpn. Ser. II, 75, 191–209, https://doi.org/10.2151/jmsj1965.75.1B_191, 1997. a
Tetley, M. G., Williams, S. E., Gurnis, M., Flament, N., and Müller, R. D.: Constraining Absolute Plate Motions Since the Triassic, J. Geophys. Res.-Sol. Ea., 124, https://doi.org/10.1029/2019JB017442, 2019. a
Tijskens, E., Roose, D., Ramon, H., and De Baerdemaeker, J.: Automatic differentiation for solving nonlinear partial differential equations: an efficient operator overloading approach, Numer. Algorithms, 30, 259–301, 2002. a
Tikhonov, A. N.: Solution of incorrectly formulated problems and the regularization method., Sov. Dok., 4, 1035–1038, 1963. a
Virtanen, P., Gommers, R., Oliphant, T. E., Haberland, M., Reddy, T., Cournapeau, D., Burovski, E., Peterson, P., Weckesser, W., Bright, J., van der Walt, S. J., Brett, M., Wilson, J., Millman, K. J., Mayorov, N., Nelson, A. R. J., Jones, E., Kern, R., Larson, E., Carey, C. J., Polat, İ., Feng, Y., Moore, E. W., VanderPlas, J., Laxalde, D., Perktold, J., Cimrman, R., Henriksen, I., Quintero, E. A., Harris, C. R., Archibald, A. M., Ribeiro, A. H., Pedregosa, F., van Mulbregt, P., and SciPy 1.0 Contributors: SciPy 1.0: fundamental algorithms for scientific computing in Python, Nature Methods, 17, 261–272, https://doi.org/10.1038/s41592-019-0686-2, 2020. a
Vynnytska, L. and Bunge, H.-P.: Restoring past mantle convection structure through fluid dynamic inverse theory: Regularisation through surface velocity boundary conditions, GEM-Int. J. Geomath., 6, 83–100, https://doi.org/10.1007/s13137-014-0060-6, 2015. a
Wang, Z. R., Stotz, I. L., Bunge, H.-P., Vilacís, B., Hayek, J. N., Ghelichkhan, S., and Lebedev, S.: Cenozoic upper mantle flow history of the Atlantic realm based on Couette/Poiseuille models: Towards paleo-mantle-flowgraphy, Phys. Earth Planet. In., 340, 107045, https://doi.org/10.1016/j.pepi.2023.107045, 2023. a
Weismüller, J., Gmeiner, B., Ghelichkhan, S., Huber, M., John, L., Wohlmuth, B., Rüde, U., and Bunge, H.-P.: Fast asthenosphere motion in high-resolution global mantle flow models, Geophys. Res. Lett., 42, 7429–7435, https://doi.org/10.1002/2015GL063727, 2015. a
Wilkinson, M. D., Dumontier, M., Aalbersberg, I. J., Appleton, G., Axton, M., Baak, A., Blomberg, N., Boiten, J.-W., da Silva Santos, L. B., Bourne, P. E., Bouwman, J., Brookes, A. J., Clark, T., Crosas, M., Dillo, I., Dumon, O., Edmunds, S., Evelo, C. T., Finkers, R., Gonzalez-Beltran, A., Gray, A. J. G., Groth, P., Goble, C., Grethe, J. S., Heringa, J., ’t Hoen, P. A. C., Hooft, R., Kuhn, T., Kok, R., Kok, J., Lusher, S. J., Martone, M. E., Mons, A., Packer, A. L., Persson, B., Rocca-Serra, P., Roos, M., van Schaik, R., Sansone, A.-A., Schultes, E., Sengstag, T., Slater, T., Strawn, G., Swertz, M. A., Thompson, M., van der Lei, J., van Mulligen, E., Velterop, J., Waagmeester, A., Wittenburg, P., Wolstencroft, K., Zhao, J., and Mons, B.: The FAIR Guiding Principles for scientific data management and stewardship, Sci. Data, 3, 1–9, 2016. a
Williams, S., Flament, N., Müller, R. D., and Butterworth, N.: Absolute plate motions since 130 Ma constrained by subduction zone kinematics, Earth Planet. Sc. Lett., 418, 66–77, 2015. a
Wolstencroft, M., Davies, J. H., and Davies, D. R.: Nusselt–Rayleigh number scaling for spherical shell Earth mantle simulation up to a Rayleigh number of 109, Phys. Earth Planet. In., 176, 132–141, 2009. a
Wunsch, C.: The Ocean Circulation Inverse Problem, Cambridge University Press, Cambridge, ISBN 9780511629570, https://doi.org/10.1017/CBO9780511629570, 1996. a
Xu, K. and Darve, E.: Physics constrained learning for data-driven inverse modeling from sparse observations, J. Comput. Phys., 453, 110938, https://doi.org/10.48550/arXiv.2002.10521, 2022. a
Young, A., Flament, N., Williams, S. E., Merdith, A., Cao, X., and Müller, R. D.: Long-term Phanerozoic sea level change from solid Earth processes, Earth Planet. Sc. Lett., 584, 117451, https://doi.org/10.1016/j.epsl.2022.117451, 2022. a, b
Zaroli, C., Sambridge, M., Lévêque, J.-J., Debayle, E., and Nolet, G.: An objective rationale for the choice of regularisation parameter with application to global multiple-frequency S-wave tomography, Solid Earth, 4, 357–371, https://doi.org/10.5194/se-4-357-2013, 2013. a
Zhang, H., Constantinescu, E. M., and Smith, B. F.: PETSc TSAdjoint: a discrete adjoint ODE solver for first-order and second-order sensitivity analysis, SIAM J. Sci. Comput., 44, C1–C24, 2022. a
Zhong, S. and Rudolph, M. L.: On the temporal evolution of long-wavelength mantle structure of the Earth since the early Paleozoic, Geochem. Geophys. Geosyst., 16, 1599–1615, https://doi.org/10.1002/2015GC005782, 2015. a
Zhong, S., Yuen, D. A., Moresi, L. N., and Knepley, M.: Numerical methods for mantle convection, Treatise on geophysics, 7, 227–252, 2007. a
Short summary
We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse modelling of Earth system processes, with an initial focus on mantle dynamics. G-ADOPT is built upon Firedrake, Dolfin-Adjoint and the Rapid Optimisation Library, which work together to optimise models using an adjoint method, aligning them with seismic and geologic datasets. We demonstrate G-ADOPT's ability to reconstruct mantle evolution and thus be a powerful tool in geosciences.
We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse...
Special issue