Articles | Volume 18, issue 20
https://doi.org/10.5194/gmd-18-7389-2025
© Author(s) 2025. 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-18-7389-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
On stabilisation of compositional density jumps in compressible mantle convection simulations
Department of Earth and Planetary Sciences, ETH Zurich, 8092 Zurich, Switzerland
Related authors
Paul James Tackley
EGUsphere, https://doi.org/10.5194/egusphere-2025-1354, https://doi.org/10.5194/egusphere-2025-1354, 2025
Short summary
Short summary
Tracers are commonly used in geodynamical models to track composition, but a common problem is that over time, gaps in the tracer distribution can develop, as well as bunches. Here a method to correct such problems is presented and tested. The method perturbs or “nudges” the positions of tracers in such a way as to close gaps and eliminate bunching. Test results show that this tracer nudging method is highly effective. The computational cost is small.
Joshua Martin Guerrero, Frédéric Deschamps, Yang Li, Wen-Pin Hsieh, and Paul James Tackley
Solid Earth, 14, 119–135, https://doi.org/10.5194/se-14-119-2023, https://doi.org/10.5194/se-14-119-2023, 2023
Short summary
Short summary
The mantle thermal conductivity's dependencies on temperature, pressure, and composition are often suppressed in numerical models. We examine the effect of these dependencies on the long-term evolution of lower-mantle thermochemical structure. We propose that depth-dependent conductivities derived from mantle minerals, along with moderate temperature and compositional correction, emulate the Earth's mean lowermost-mantle conductivity values and produce a stable two-pile configuration.
Anna Johanna Pia Gülcher, Maxim Dionys Ballmer, and Paul James Tackley
Solid Earth, 12, 2087–2107, https://doi.org/10.5194/se-12-2087-2021, https://doi.org/10.5194/se-12-2087-2021, 2021
Short summary
Short summary
The lower mantle extends from 660–2890 km depth, making up > 50 % of the Earth’s volume. Its composition and structure, however, remain poorly understood. In this study, we investigate several hypotheses with computer simulations of mantle convection that include different materials: recycled, dense rocks and ancient, strong rocks. We propose a new integrated style of mantle convection including
piles,
blobs, and
streaksthat agrees with various observations of the deep Earth.
Daniela Paz Bolrão, Maxim D. Ballmer, Adrien Morison, Antoine B. Rozel, Patrick Sanan, Stéphane Labrosse, and Paul J. Tackley
Solid Earth, 12, 421–437, https://doi.org/10.5194/se-12-421-2021, https://doi.org/10.5194/se-12-421-2021, 2021
Short summary
Short summary
We use numerical models to investigate the thermo-chemical evolution of a solid mantle during a magma ocean stage. When applied to the Earth, our study shows that the solid mantle and a magma ocean tend toward chemical equilibration before crystallisation of this magma ocean. Our findings suggest that a very strong chemical stratification of the solid mantle is unlikely to occur (as predicted by previous studies), which may explain why the Earth’s mantle is rather homogeneous in composition.
Paul James Tackley
EGUsphere, https://doi.org/10.5194/egusphere-2025-1354, https://doi.org/10.5194/egusphere-2025-1354, 2025
Short summary
Short summary
Tracers are commonly used in geodynamical models to track composition, but a common problem is that over time, gaps in the tracer distribution can develop, as well as bunches. Here a method to correct such problems is presented and tested. The method perturbs or “nudges” the positions of tracers in such a way as to close gaps and eliminate bunching. Test results show that this tracer nudging method is highly effective. The computational cost is small.
Joshua Martin Guerrero, Frédéric Deschamps, Yang Li, Wen-Pin Hsieh, and Paul James Tackley
Solid Earth, 14, 119–135, https://doi.org/10.5194/se-14-119-2023, https://doi.org/10.5194/se-14-119-2023, 2023
Short summary
Short summary
The mantle thermal conductivity's dependencies on temperature, pressure, and composition are often suppressed in numerical models. We examine the effect of these dependencies on the long-term evolution of lower-mantle thermochemical structure. We propose that depth-dependent conductivities derived from mantle minerals, along with moderate temperature and compositional correction, emulate the Earth's mean lowermost-mantle conductivity values and produce a stable two-pile configuration.
Anna Johanna Pia Gülcher, Maxim Dionys Ballmer, and Paul James Tackley
Solid Earth, 12, 2087–2107, https://doi.org/10.5194/se-12-2087-2021, https://doi.org/10.5194/se-12-2087-2021, 2021
Short summary
Short summary
The lower mantle extends from 660–2890 km depth, making up > 50 % of the Earth’s volume. Its composition and structure, however, remain poorly understood. In this study, we investigate several hypotheses with computer simulations of mantle convection that include different materials: recycled, dense rocks and ancient, strong rocks. We propose a new integrated style of mantle convection including
piles,
blobs, and
streaksthat agrees with various observations of the deep Earth.
Daniela Paz Bolrão, Maxim D. Ballmer, Adrien Morison, Antoine B. Rozel, Patrick Sanan, Stéphane Labrosse, and Paul J. Tackley
Solid Earth, 12, 421–437, https://doi.org/10.5194/se-12-421-2021, https://doi.org/10.5194/se-12-421-2021, 2021
Short summary
Short summary
We use numerical models to investigate the thermo-chemical evolution of a solid mantle during a magma ocean stage. When applied to the Earth, our study shows that the solid mantle and a magma ocean tend toward chemical equilibration before crystallisation of this magma ocean. Our findings suggest that a very strong chemical stratification of the solid mantle is unlikely to occur (as predicted by previous studies), which may explain why the Earth’s mantle is rather homogeneous in composition.
Cited articles
Andrés-Martínez, M., Morgan, J. P., Pérez-Gussinyé, M., and Rüpke, L.: A new free-surface stabilization algorithm for geodynamical modelling: Theory and numerical tests, Phys. Earth Planet. Int., 246, 41–51, https://doi.org/10.1016/j.pepi.2015.07.003, 2015.
Bezanson, J., Edelman, A., Karpinski, S., and Shah, V. B.: Julia: A fresh approach to numerical computing, SIAM Review, 59, 65–98, https://doi.org/10.1137/141000671, 2017.
Christensen, U. R., and Hofmann, A. W.: Segregation of subducted oceanic crust in the convecting mantle, J. Geophys. Res., 99, 19867–19884, https://doi.org/10.1029/93JB03403, 1994.
Crameri, F., Schmeling, H., Golabek, G. J., Duretz, T., Orendt, R., Buiter, S. J. H., May, D. A., Kaus, B. J. P., Gerya, T. V., and Tackley, P. J.: A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the “sticky air” method, Geophysical Journal International, 189, 38–54, https://doi.org/10.1111/j.1365-246X.2012.05388.x, 2012.
Davaille, A.: Simultaneous generation of hotspots and superswells by convection in a heterogeneous planetary mantle, Nature, 402, 756–760, https://doi.org/10.1038/45461, 1999.
Deschamps, F., Kaminski, E., and Tackley, P. J.: A deep mantle origin for the primitive signature of ocean island basalt, Nature Geosci., 4, 879–882, https://doi.org/10.1038/ngeo1295, 2011.
Duretz, T., May, D. A., Gerya, T. V., and Tackley, P. J.: Discretisation errors and free surface stabilization in the finite difference and marker-in-cell method in geodynamic applications: A numerical study, Geochem. Geophys. Geosyst., 12, https://doi.org/10.1029/2011GC003567, 2011.
Dziewonski, A. M. and Anderson, D. L.: Preliminary reference Earth model, Phys. Earth Planet. Inter, 25, 297–356, https://doi.org/10.1016/0031-9201(81)90046-7, 1981.
Gassmöller, R., Dannberg, J., Bangerth, W., Heister, T., and Myhill, R.: On formulations of compressible mantle convection, Geophys. J. Int., 221, 1264–1280, https://doi.org/10.1093/gji/ggaa078, 2020.
Gerya, T. V. and Yuen, D. A.: Robust characteristics method for modelling multiphase visco-elasto-plastic thermo-mechanical problems, Physics of the Earth and Planetary Interiors, 163, 83–105, https://doi.org/10.1016/j.pepi.2007.04.015, 2007.
Gurnis, M.: The effects of chemical density differences on convective mixing in the Earth's mantle, J. Geophys. Res., 91, 1407–1419, https://doi.org/10.1029/JB091iB11p11407, 1986.
Harlow, F. H. and Welch, J. E.: Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface, Phys. Fluids A, 8, 2182–2189, https://doi.org/10.1063/1.1761178, 1965.
Kameyama, M., Kageyama, A., and Sato, T.: Multigrid-based simulation code for mantle convection in spherical shell using yin-yang grid, Phys. Earth Planet. Int., 171, 19–32, https://doi.org/10.1016/j.pepi.2008.06.025, 2008.
Kaus, B. J. P., Mühlhaus, H., and May, D. A.: A stabilization algorithm for geodynamic numerical simulations with a free surface, Phys. Earth Planet. Int., 181, 12–20, https://doi.org/10.1016/j.pepi.2010.04.007, 2010.
Kaus, B. J. P., Popov, A. A., Baumann, T. S., Püsök, A. E., Bauville, A., Fernandez, N., and Collignon, M.: Forward and inverse modelling of lithospheric deformation on geological timescales, NIC Series, 48, 299, ISBN 978-3-95806-109-5, 2016.
King, S. D., Lee, C., van Keken, P. E., Leng, W., Zhong, S., Tan, E., Tosi, N., and Kameyama, M. C.: A community benchmark for 2-D Cartesian compressible convection in the Earth's mantle, Geophysical Journal International, 180, 73–87, https://doi.org/10.1111/j.1365-246X.2009.04413.x, 2010.
Kramer, S. C., Wilson, C. R., and Davies, D. R.: An implicit free surface algorithm for geodynamical simulations, Physics of the Earth and Planetary Interiors, 194–195, 25–37, https://doi.org/10.1016/j.pepi.2012.01.001, 2012.
Lin, J., Gerya, T. V., Tackley, P. J., Yuen, D. A., and Golabek, G. J.: Numerical modeling of protocore destabilization during planetary accretion: Feedbacks from non-Newtonian rheology and energy dissipation, Icarus, 213, 24–42, https://doi.org/10.1016/j.icarus.2009.06.035, 2011.
Nakagawa, T. and Tackley, P. J.: Influence of plate tectonic mode on the coupled thermochemical evolution of Earth's mantle and core, Geochem. Geophys. Geosys., 16, 3400–3413, https://doi.org/10.1002/2015gc005996, 2015.
Ogawa, M.: Numerical models of magmatism in convecting mantle with temperature-dependent viscosity and their implications for Venus and Earth, J. Geophys. Res., 105, 6997–7012, https://doi.org/10.1029/1999JE001162, 2000.
Ogawa, M., Schubert, G., and Zebib, A.: Numerical simulations of 3-dimensional thermal convection in a fluid with strongly temperature-dependent viscosity, J. Fluid Mech., 233, 299–328, https://doi.org/10.1017/S0022112091000496, 1991.
Patankar, S. V.: Numerical Heat Transfer and Fluid Flow, CRC Press, https://doi.org/10.1201/9781482234213, 1980.
Popov, A. A. and Sobolev, S. V.: SLIM3D: A tool for three-dimensional thermomechanical modeling of lithospheric deformation with elasto-visco-plastic rheology, Physics of the Earth and Planetary Interiors, 171, 55–75, https://doi.org/10.1016/j.pepi.2008.03.007, 2008.
Rose, I., Buffett, B., and Heister, T.: Stability and accuracy of free surface time integration in viscous flows, Physics of the Earth and Planetary Interiors, 262, 90–100, https://doi.org/10.1016/j.pepi.2016.11.007, 2017.
Samuel, H. and Tackley, P. J.: Dynamics of core formation and equilibrium by negative diapirism, Geochem. Geophys. Geosyst., 9, https://doi.org/10.1029/2007GC001896, 2008.
Schubert, G., Turcotte, D. L., and Olson, P.: Mantle Convection in the Earth and Planets, Cambridge University Press, https://doi.org/10.1017/CBO9780511612879, 2001.
Tackley, P. J.: Effects of strongly temperature-dependent viscosity on time-dependent, 3-dimensional models of mantle convection, Geophys. Res. Lett., 20, 2187–2190, https://doi.org/10.1029/93GL02317, 1993.
Tackley, P. J.: Effects of strongly variable viscosity on three-dimensional compressible convection in planetary mantles, J. Geophys. Res., 101, 3311–3332, https://doi.org/10.1029/95JB03211, 1996.
Tackley, P. J.: Three-dimensional simulations of mantle convection with a thermochemical CMB boundary layer: D”?, in: The Core-Mantle Boundary Region, edited by: Gurnis, M., Wysession, M. E., Knittle, E., and Buffett, B. A., American Geophysical Union, 231–253, https://doi.org/10.1029/GD028, 1998.
Tackley, P. J.: Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid, Phys. Earth Planet. Int., https://doi.org/10.1016/j.pepi.2008.08.005, 2008.
Tackley, P. J. and ETH Zurich: Software for manuscript “On stabilisation of compositional density jumps in compressible mantle convection simulations” (1.1), Zenodo [code], https://doi.org/10.5281/zenodo.15115816, 2025.
Trompert, R. A. and Hansen, U.: The application of a finite-volume multigrid method to 3-dimensional flow problems in a highly viscous fluid with a variable viscosity, Geophys. Astrophys. Fluid Dyn., 83, 261–291, https://doi.org/10.1080/03091929608208968, 1996.
Short summary
Large density jumps in numerical simulations of solid Earth dynamics can cause numerical oscillations. An effective method to prevent these at a free surface already exists. Here this is tested for compositional layers deeper in the mantle. The stabilisation method works effectively if density gradients due purely to compositional gradients are used but produces severe artefacts if total density is used.
Large density jumps in numerical simulations of solid Earth dynamics can cause numerical...
Special issue