Articles | Volume 5, issue 2
https://doi.org/10.5194/se-5-1087-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
https://doi.org/10.5194/se-5-1087-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Method article
| 
18 Nov 2014
Method article |
| 18 Nov 2014
Using the level set method in geodynamical modeling of multi-material flows and Earth's free surface
B. Hillebrand
CORRESPONDING AUTHOR
Department of Earth Sciences, Utrecht University, the Netherlands
C. Thieulot
Department of Earth Sciences, Utrecht University, the Netherlands
Centre of Earth Evolution and Dynamics (CEED), University of Oslo, 0316 Oslo, Norway
T. Geenen
SurfSara, Amsterdam, the Netherlands
A. P. van den Berg
Department of Earth Sciences, Utrecht University, the Netherlands
W. Spakman
Department of Earth Sciences, Utrecht University, the Netherlands
Centre of Earth Evolution and Dynamics (CEED), University of Oslo, 0316 Oslo, Norway
Related authors
No articles found.
Erik van der Wiel, Cedric Thieulot, and Douwe J. J. van Hinsbergen
Solid Earth, 15, 861–875, https://doi.org/10.5194/se-15-861-2024, https://doi.org/10.5194/se-15-861-2024, 2024
Short summary
Short summary
Geodynamic models of mantle convection provide a powerful tool to study the structure and composition of the Earth's mantle. Comparing such models with other datasets is difficult. We explore the use of
configurational entropy, which allows us to quantify mixing in models. The entropy may be used to analyse the mixed state of the mantle as a whole and may also be useful to validate numerical models against anomalies in the mantle that are obtained from seismology and geochemistry.
Cedric Thieulot and Wolfgang Bangerth
EGUsphere, https://doi.org/10.5194/egusphere-2024-1668, https://doi.org/10.5194/egusphere-2024-1668, 2024
Short summary
Short summary
One of the main numerical methods in geodynamics is the finite-element method. Many types of elements have been used in the past decades in hundreds of publications. They usually fall under two categories: quadrilaterals and triangles. For the first time we compare results obtained with the most used elements of each type on a series of geodynamical benchmarks and draw conclusions as to which are the best ones and which are to be preferably avoided.
Rene Gassmöller, Juliane Dannberg, Wolfgang Bangerth, Elbridge Gerry Puckett, and Cedric Thieulot
Geosci. Model Dev., 17, 4115–4134, https://doi.org/10.5194/gmd-17-4115-2024, https://doi.org/10.5194/gmd-17-4115-2024, 2024
Short summary
Short summary
Numerical models that use simulated particles are a powerful tool for investigating flow in the interior of the Earth, but the accuracy of these models is not fully understood. Here we present two new benchmarks that allow measurement of model accuracy. We then document that better accuracy matters for applications like convection beneath an oceanic plate. Our benchmarks and methods are freely available to help the community develop better models.
Iris van Zelst, Cedric Thieulot, and Timothy J. Craig
Solid Earth, 14, 683–707, https://doi.org/10.5194/se-14-683-2023, https://doi.org/10.5194/se-14-683-2023, 2023
Short summary
Short summary
A common simplification in subduction zone models is the use of constant thermal parameters, while experiments have shown that they vary with temperature. We test various formulations of temperature-dependent thermal parameters and show that they change the thermal structure of the subducting slab. We recommend that modelling studies of the thermal structure of subduction zones take the temperature dependence of thermal parameters into account, especially when providing insights into seismicity.
Barend Cornelis Root, Josef Sebera, Wolfgang Szwillus, Cedric Thieulot, Zdeněk Martinec, and Javier Fullea
Solid Earth, 13, 849–873, https://doi.org/10.5194/se-13-849-2022, https://doi.org/10.5194/se-13-849-2022, 2022
Short summary
Short summary
Several alternative gravity modelling techniques and associated numerical codes with their own advantages and limitations are available for the solid Earth community. With upcoming state-of-the-art lithosphere density models and accurate global gravity field data sets, it is vital to understand the differences of the various approaches. In this paper, we discuss the four widely used techniques: spherical harmonics, tesseroid integration, triangle integration, and hexahedral integration.
Iris van Zelst, Fabio Crameri, Adina E. Pusok, Anne Glerum, Juliane Dannberg, and Cedric Thieulot
Solid Earth, 13, 583–637, https://doi.org/10.5194/se-13-583-2022, https://doi.org/10.5194/se-13-583-2022, 2022
Short summary
Short summary
Geodynamic modelling provides a powerful tool to investigate processes in the Earth’s crust, mantle, and core that are not directly observable. In this review, we present a comprehensive yet concise overview of the modelling process with an emphasis on best practices. We also highlight synergies with related fields, such as seismology and geology. Hence, this review is the perfect starting point for anyone wishing to (re)gain a solid understanding of geodynamic modelling as a whole.
Cedric Thieulot and Wolfgang Bangerth
Solid Earth, 13, 229–249, https://doi.org/10.5194/se-13-229-2022, https://doi.org/10.5194/se-13-229-2022, 2022
Short summary
Short summary
One of the main numerical methods to solve the mass, momentum, and energy conservation equations in geodynamics is the finite-element method. Four main types of elements have been used in the past decades in hundreds of publications. For the first time we compare results obtained with these four elements on a series of geodynamical benchmarks and applications and draw conclusions as to which are the best ones and which are to be preferably avoided.
Melchior Schuh-Senlis, Cedric Thieulot, Paul Cupillard, and Guillaume Caumon
Solid Earth, 11, 1909–1930, https://doi.org/10.5194/se-11-1909-2020, https://doi.org/10.5194/se-11-1909-2020, 2020
Short summary
Short summary
This paper presents a numerical method for restoring models of the subsurface to a previous state in their deformation history, acting as a numerical time machine for geological structures. The method relies on the assumption that rock layers can be modeled as highly viscous fluids. It shows promising results on simple setups, including models with faults and non-flat topography. While issues still remain, this could open a way to add more physics to reverse time structural modeling.
Menno Fraters, Cedric Thieulot, Arie van den Berg, and Wim Spakman
Solid Earth, 10, 1785–1807, https://doi.org/10.5194/se-10-1785-2019, https://doi.org/10.5194/se-10-1785-2019, 2019
Short summary
Short summary
Three-dimensional numerical modelling of geodynamic processes may benefit strongly from using realistic 3-D starting models that approximate, e.g. natural subduction settings in the geological past or at present. To this end, we developed the Geodynamic World Builder (GWB), which enables relatively straightforward parameterization of complex 3-D geometric structures associated with geodynamic processes. The GWB is an open-source community code designed to easily interface with geodynamic codes.
Cedric Thieulot
Solid Earth, 9, 1169–1177, https://doi.org/10.5194/se-9-1169-2018, https://doi.org/10.5194/se-9-1169-2018, 2018
Short summary
Short summary
I present the GHOST (Geoscientific Hollow Sphere Tessellation) software which allows for the fast generation of computational meshes in hollow sphere geometries counting up to a hundred million cells. Each mesh is composed of concentric spherical shells made of quadrilaterals or triangles. I focus here on three commonly used meshes used in the geodynamics/geophysics community and further benchmark the gravity and gravitational potential procedures in the simple case of a constant density.
Alexis Plunder, Cédric Thieulot, and Douwe J. J. van Hinsbergen
Solid Earth, 9, 759–776, https://doi.org/10.5194/se-9-759-2018, https://doi.org/10.5194/se-9-759-2018, 2018
Short summary
Short summary
The thermal state of the Earth's crust determines how it reacts to tectonic forces and to fluid flow responsible for ore formation. We hypothesize that the angle between plate motion and convergent boundaries determines the thermal regime of subduction zones (where a plate goes under another one). Computer models and a geological reconstruction of Turkey were used to validate this hypothesis.
This research was done to validate a hypothesis made on the basis of nonquantitative field data.
Anne Glerum, Cedric Thieulot, Menno Fraters, Constantijn Blom, and Wim Spakman
Solid Earth, 9, 267–294, https://doi.org/10.5194/se-9-267-2018, https://doi.org/10.5194/se-9-267-2018, 2018
Short summary
Short summary
A nonlinear viscoplastic rheology is implemented and benchmarked in the ASPECT software, allowing for the modeling of lithospheric deformation. We showcase the new functionality with a four-dimensional model of thermomechanically coupled subduction.
Cedric Thieulot
Solid Earth, 8, 1181–1191, https://doi.org/10.5194/se-8-1181-2017, https://doi.org/10.5194/se-8-1181-2017, 2017
Short summary
Short summary
I present a new family of analytical flow solutions to the incompressible Stokes equation in a spherical shell. The velocity is tangential to both inner and outer boundaries, the viscosity is radial, and the solution has been designed so that the expressions for velocity, pressure, and body force are simple to implement in (geodynamics) codes. This forms the basis of a numerical benchmark for convection codes, and I have implemented it in two finite-element codes.
C. Thieulot
Solid Earth Discuss., https://doi.org/10.5194/sed-6-1949-2014, https://doi.org/10.5194/sed-6-1949-2014, 2014
Revised manuscript has not been submitted
Related subject area
Geodynamics
Increased metamorphic conditions in the lower crust during oceanic transform fault evolution
ECOMAN: an open-source package for geodynamic and seismological modelling of mechanical anisotropy
How a volcanic arc influences back-arc extension: insight from 2D numerical models
Quantifying mantle mixing through configurational entropy
Various lithospheric deformation patterns derived from rheological contrasts between continental terranes: insights from 2-D numerical simulations
Magmatic underplating associated with Proterozoic basin formation: insights from gravity study over the southern margin of the Bundelkhand Craton, India
On the impact of true polar wander on heat flux patterns at the core–mantle boundary
The influence of viscous slab rheology on numerical models of subduction
Statistical appraisal of geothermal heat flow observations in the Arctic
Fast uplift in the southern Patagonian Andes due to long- and short-term deglaciation and the asthenospheric window underneath
Modeling liquid transport in the Earth's mantle as two-phase flow: effect of an enforced positive porosity on liquid flow and mass conservation
Thrusts control the thermal maturity of accreted sediments
The crustal structure of the Longmenshan fault zone and its implications for seismogenesis: new insight from aeromagnetic and gravity data
Earth's core variability from magnetic and gravity field observations
The role of continental lithospheric thermal structure in the evolution of orogenic systems: application to the Himalayan–Tibetan collision zone
Glacial-isostatic-adjustment strain rate–stress paradox in the Western Alps and impact on active faults and seismicity
The effect of temperature-dependent material properties on simple thermal models of subduction zones
Plume–ridge interactions: ridgeward versus plate-drag plume flow
Transport mechanisms of hydrothermal convection in faulted tight sandstones
A corrected finite-difference scheme for the flexure equation with abrupt changes in coefficient
Influence of heterogeneous thermal conductivity on the long-term evolution of the lower-mantle thermochemical structure: implications for primordial reservoirs
The role of edge-driven convection in the generation ofvolcanism – Part 2: Interaction with mantle plumes, applied to the Canary Islands
The effect of low-viscosity sediments on the dynamics and accretionary style of subduction margins
Thermal non-equilibrium of porous flow in a resting matrix applicable to melt migration: a parametric study
101 geodynamic modelling: how to design, interpret, and communicate numerical studies of the solid Earth
Crustal structure of the Volgo–Uralian subcraton revealed by inverse and forward gravity modelling
On the choice of finite element for applications in geodynamics
A new finite element approach to model microscale strain localization within olivine aggregates
Interpolation of magnetic anomalies over an oceanic ridge region using an equivalent source technique and crust age model constraint
Coupled dynamics and evolution of primordial and recycled heterogeneity in Earth's lower mantle
Buoyancy versus shear forces in building orogenic wedges
Comparing global seismic tomography models using varimax principal component analysis
Magma ascent mechanisms in the transition regime from solitary porosity waves to diapirism
Analytical solution for residual stress and strain preserved in anisotropic inclusion entrapped in an isotropic host
Gravity effect of Alpine slab segments based on geophysical and petrological modelling
The role of edge-driven convection in the generation of volcanism – Part 1: A 2D systematic study
Gravity modeling of the Alpine lithosphere affected by magmatism based on seismic tomography
Timescales of chemical equilibrium between the convecting solid mantle and over- and underlying magma oceans
The preserved plume of the Caribbean Large Igneous Plateau revealed by 3D data-integrative models
Impact of upper mantle convection on lithosphere hyperextension and subsequent horizontally forced subduction initiation
Pragmatic solvers for 3D Stokes and elasticity problems with heterogeneous coefficients: evaluating modern incomplete LDLT preconditioners
Combined numerical and experimental study of microstructure and permeability in porous granular media
Mapping undercover: integrated geoscientific interpretation and 3D modelling of a Proterozoic basin
Monitoring crustal CO2 flow: methods and their applications to the mofettes in West Bohemia
On the self-regulating effect of grain size evolution in mantle convection models: application to thermochemical piles
Deciphering the metamorphic evolution of the Pulo do Lobo metasedimentary domain (SW Iberian Variscides)
The impact of rheological uncertainty on dynamic topography predictions
The effect of effective rock viscosity on 2-D magmatic porosity waves
Density distribution across the Alpine lithosphere constrained by 3-D gravity modelling and relation to seismicity and deformation
Pore-scale permeability prediction for Newtonian and non-Newtonian fluids
Peter Haas, Myron F. H. Thomas, Christian Heine, Jörg Ebbing, Andrey Seregin, and Jimmy van Itterbeeck
Solid Earth, 15, 1419–1443, https://doi.org/10.5194/se-15-1419-2024, https://doi.org/10.5194/se-15-1419-2024, 2024
Short summary
Short summary
Transform faults are conservative plate boundaries where no material is added or destroyed. Oceanic fracture zones are their inactive remnants and record tectonic processes that formed oceanic crust. In this study, we combine high-resolution data sets along fracture zones in the Gulf of Guinea to demonstrate that their formation is characterized by increased metamorphic conditions. This is in line with previous studies that describe the non-conservative character of transform faults.
Manuele Faccenda, Brandon P. VanderBeek, Albert de Montserrat, Jianfeng Yang, Francesco Rappisi, and Neil Ribe
Solid Earth, 15, 1241–1264, https://doi.org/10.5194/se-15-1241-2024, https://doi.org/10.5194/se-15-1241-2024, 2024
Short summary
Short summary
The Earth's internal dynamics and structure can be well understood by combining seismological and geodynamic modelling with mineral physics, an approach that has been poorly adopted in the past. To this end we have developed ECOMAN, an open-source software package that is intended to overcome the computationally intensive nature of this multidisciplinary methodology and the lack of a dedicated and comprehensive computational framework.
Duo Zhang and J. Huw Davies
Solid Earth, 15, 1113–1132, https://doi.org/10.5194/se-15-1113-2024, https://doi.org/10.5194/se-15-1113-2024, 2024
Short summary
Short summary
We numerically model the influence of an arc on back-arc extension. The arc is simulated by placing a hot region on the overriding plate. We investigate how plate ages and properties of the hot region affect back-arc extension and present regime diagrams illustrating the nature of back-arc extension for these models. We find that back-arc extension occurs not only in the hot region but also, surprisingly, away from it, and a hot region facilitates extension on the overriding plate.
Erik van der Wiel, Cedric Thieulot, and Douwe J. J. van Hinsbergen
Solid Earth, 15, 861–875, https://doi.org/10.5194/se-15-861-2024, https://doi.org/10.5194/se-15-861-2024, 2024
Short summary
Short summary
Geodynamic models of mantle convection provide a powerful tool to study the structure and composition of the Earth's mantle. Comparing such models with other datasets is difficult. We explore the use of
configurational entropy, which allows us to quantify mixing in models. The entropy may be used to analyse the mixed state of the mantle as a whole and may also be useful to validate numerical models against anomalies in the mantle that are obtained from seismology and geochemistry.
Renxian Xie, Lin Chen, Jason P. Morgan, and Yongshun John Chen
Solid Earth, 15, 789–806, https://doi.org/10.5194/se-15-789-2024, https://doi.org/10.5194/se-15-789-2024, 2024
Short summary
Short summary
Continental terranes have various rheological strengths due to the differences in their ages, compositions, and structures. We applied four assumed rheological models to three terranes in a collisional model and obtained four styles of lithosphere deformation patterns of collision, subduction, thickening/delamination, and replacement. These simulation patterns are seen in observed lithosphere deformation patterns and structures in East Asia.
Ananya Parthapradip Mukherjee and Animesh Mandal
Solid Earth, 15, 711–729, https://doi.org/10.5194/se-15-711-2024, https://doi.org/10.5194/se-15-711-2024, 2024
Short summary
Short summary
Global gravity data are used to develop 2D models and a Moho depth map from 3D inversion, depicting the crustal structure below the region covered by Proterozoic sedimentary basins, south of the Bundelkhand Craton in central India. The observed thick mafic underplated layer above the Moho indicates Proterozoic plume activity. Thus, the study offers insights into the crustal configuration of this region, illustrating the geodynamic processes that led to the formation of the basins.
Thomas Frasson, Stéphane Labrosse, Henri-Claude Nataf, Nicolas Coltice, and Nicolas Flament
Solid Earth, 15, 617–637, https://doi.org/10.5194/se-15-617-2024, https://doi.org/10.5194/se-15-617-2024, 2024
Short summary
Short summary
Heat flux heterogeneities at the bottom of Earth's mantle play an important role in the dynamic of the underlying core. Here, we study how these heterogeneities are affected by the global rotation of the Earth, called true polar wander (TPW), which has to be considered to relate mantle dynamics with core dynamics. We find that TPW can greatly modify the large scales of heat flux heterogeneities, notably at short timescales. We provide representative maps of these heterogeneities.
Natalie Hummel, Susanne Buiter, and Zoltán Erdős
Solid Earth, 15, 567–587, https://doi.org/10.5194/se-15-567-2024, https://doi.org/10.5194/se-15-567-2024, 2024
Short summary
Short summary
Simulations of subducting tectonic plates often use material properties extrapolated from the behavior of small rock samples in a laboratory to conditions found in the Earth. We explore several typical approaches to simulating these extrapolated material properties and show that they produce very rigid subducting plates with unrealistic dynamics. Our findings imply that subducting plates deform by additional mechanisms that are less commonly implemented in simulations.
Judith Freienstein, Wolfgang Szwillus, Agnes Wansing, and Jörg Ebbing
Solid Earth, 15, 513–533, https://doi.org/10.5194/se-15-513-2024, https://doi.org/10.5194/se-15-513-2024, 2024
Short summary
Short summary
Geothermal heat flow influences ice sheet dynamics, making its investigation important for ice-covered regions. Here we evaluate the sparse measurements for their agreement with regional solid Earth models, as well as with a statistical approach. This shows that some points should be excluded from regional studies. In particular, the NGRIP point, which strongly influences heat flow maps and the distribution of high basal melts, should be statistically considered an outlier.
Veleda A. P. Muller, Pietro Sternai, and Christian Sue
Solid Earth, 15, 387–404, https://doi.org/10.5194/se-15-387-2024, https://doi.org/10.5194/se-15-387-2024, 2024
Short summary
Short summary
Uplift rates up to 40 mm yr−1 are measured by GNSS in the southern Patagonian Icefield, a remainder of the Patagonian Ice Sheet that covered the Andes in the Last Glacial Maximum (LGM) at 26 ka. Using numerical modelling, we estimate an increase of 150 to 200 °C of the asthenospheric temperature due to the slab window under southern Patagonia, and we show that post-glacial rebound, after the long-term LGM and the short-term Little Ice Age (400 a), contributed to the modern uplift rate budget.
Changyeol Lee, Nestor G. Cerpa, Dongwoo Han, and Ikuko Wada
Solid Earth, 15, 23–38, https://doi.org/10.5194/se-15-23-2024, https://doi.org/10.5194/se-15-23-2024, 2024
Short summary
Short summary
Fluids and melts in the mantle are key to the Earth’s evolution. The main driving force for their transport is the compaction of the porous mantle. Numerically, the compaction equations can yield unphysical negative liquid fractions (porosity), and it is necessary to enforce positive porosity. However, the effect of such a treatment on liquid flow and mass conservation has not been quantified. We found that although mass conservation is affected, the liquid pathways are well resolved.
Utsav Mannu, David Fernández-Blanco, Ayumu Miyakawa, Taras Gerya, and Masataka Kinoshita
Solid Earth, 15, 1–21, https://doi.org/10.5194/se-15-1-2024, https://doi.org/10.5194/se-15-1-2024, 2024
Short summary
Short summary
Accretion during subduction, in which one tectonic plate moves under another, forms a wedge where sediments can be transformed into hydrocarbons. We utilised realistic computer models to investigate this and, in particular, how accretion affects mobility in the wedge and found that the evolution of the wedge and the thrusts it develops fundamentally control the thermal maturity of sediments. This can help us better understand the history of subduction and the formation of hydrocarbons in wedges.
Hai Yang, Shengqing Xiong, Qiankun Liu, Fang Li, Zhiye Jia, Xue Yang, Haofei Yan, and Zhaoliang Li
Solid Earth, 14, 1289–1308, https://doi.org/10.5194/se-14-1289-2023, https://doi.org/10.5194/se-14-1289-2023, 2023
Short summary
Short summary
The Wenchuan (Ms 8.0) and Lushan (Ms 7.0) earthquakes show different geodynamic features and form a 40–60 km area void of aftershocks for both earthquakes. The inverse models suggest that the downward-subducted basement of the Sichuan Basin is irregular in shape and heterogeneous in magnetism and density. The different focal mechanisms of the two earthquakes and the genesis of the seismic gap may be closely related to the differential thrusting mechanism caused by basement heterogeneity.
Anita Thea Saraswati, Olivier de Viron, and Mioara Mandea
Solid Earth, 14, 1267–1287, https://doi.org/10.5194/se-14-1267-2023, https://doi.org/10.5194/se-14-1267-2023, 2023
Short summary
Short summary
To understand core dynamics, insight from several possible observables is needed. By applying several separation methods, we show spatiotemporal variabilities in the magnetic and gravity fields related to the core dynamics. A 7-year oscillation is found in all conducted analyses. The results in the magnetic field reflect the core processes and the variabilities in the gravity field exhibit new findings that might be an interesting input to build an enhanced model of the Earth’s core.
Mengxue Liu, Dinghui Yang, and Rui Qi
Solid Earth, 14, 1155–1168, https://doi.org/10.5194/se-14-1155-2023, https://doi.org/10.5194/se-14-1155-2023, 2023
Short summary
Short summary
The continuous subduction mainly occurs with a relatively cold overriding lithosphere (Tmoho ≤ 450 °C), while slab break-off dominates when the model has a relatively hot procontinental Moho temparature (Tmoho ≥ 500 °C). Hr is more prone to facilitating the deformation of the lithospheric upper part than altering the collision mode. The lithospheric thermal structure may have played a significant role in the development of Himalayan–Tibetan orogenic lateral heterogeneity.
Juliette Grosset, Stéphane Mazzotti, and Philippe Vernant
Solid Earth, 14, 1067–1081, https://doi.org/10.5194/se-14-1067-2023, https://doi.org/10.5194/se-14-1067-2023, 2023
Short summary
Short summary
In glaciated regions, induced lithosphere deformation is proposed as a key process contributing to fault activity and seismicity. We study the impact of this effect on fault activity in the Western Alps. We show that the response to the last glaciation explains a major part of the geodetic strain rates but does not drive or promote the observed seismicity. Thus, seismic hazard studies in the Western Alps require detailed modeling of the glacial isostatic adjustment (GIA) transient impact.
Iris van Zelst, Cedric Thieulot, and Timothy J. Craig
Solid Earth, 14, 683–707, https://doi.org/10.5194/se-14-683-2023, https://doi.org/10.5194/se-14-683-2023, 2023
Short summary
Short summary
A common simplification in subduction zone models is the use of constant thermal parameters, while experiments have shown that they vary with temperature. We test various formulations of temperature-dependent thermal parameters and show that they change the thermal structure of the subducting slab. We recommend that modelling studies of the thermal structure of subduction zones take the temperature dependence of thermal parameters into account, especially when providing insights into seismicity.
Fengping Pang, Jie Liao, Maxim D. Ballmer, and Lun Li
Solid Earth, 14, 353–368, https://doi.org/10.5194/se-14-353-2023, https://doi.org/10.5194/se-14-353-2023, 2023
Short summary
Short summary
Plume–ridge interaction is an intriguing geological process in plate tectonics. In this paper, we address the respective role of ridgeward vs. plate-drag plume flow in 2D thermomechanical models and compare the results with a compilation of observations on Earth. From a geophysical and geochemical analysis of Earth plumes and in combination with the model results, we propose that the absence of plumes interacting with ridges in the Pacific is largely caused by the presence of plate drag.
Guoqiang Yan, Benjamin Busch, Robert Egert, Morteza Esmaeilpour, Kai Stricker, and Thomas Kohl
Solid Earth, 14, 293–310, https://doi.org/10.5194/se-14-293-2023, https://doi.org/10.5194/se-14-293-2023, 2023
Short summary
Short summary
The physical processes leading to the kilometre-scale thermal anomaly in faulted tight sandstones are numerically investigated. The fluid-flow pathways, heat-transfer types and interactions among different convective and advective flow modes are systematically identified. The methodologies and results can be applied to interpret hydrothermal convection-related geological phenomena and to draw implications for future petroleum and geothermal exploration and exploitation in analogous settings.
David Hindle and Olivier Besson
Solid Earth, 14, 197–212, https://doi.org/10.5194/se-14-197-2023, https://doi.org/10.5194/se-14-197-2023, 2023
Short summary
Short summary
By making a change to the way we solve the flexure equation that describes how the Earth's outer layer bends when it is subjected to loading by ice sheets or mountains, we develop new ways of using an old method from geodynamics. This lets us study the Earth's outer layer by measuring a parameter called the elastic thickness, effectively how stiff and springy the outer layer is when it gets loaded and also how the Earth's outer layer gets broken around its edges and in its interior.
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.
Antonio Manjón-Cabeza Córdoba and Maxim D. Ballmer
Solid Earth, 13, 1585–1605, https://doi.org/10.5194/se-13-1585-2022, https://doi.org/10.5194/se-13-1585-2022, 2022
Short summary
Short summary
The origin of many volcanic archipelagos on the Earth remains uncertain. By using 3D modelling of mantle flow and melting, we investigate the interaction between the convective mantle near the continental–oceanic transition and rising hot plumes. We believe that this phenomenon is the origin behind some archipelagos, in particular the Canary Islands. Analysing our results, we reconcile observations that were previously enigmatic, such as the complex patterns of volcanism in the Canaries.
Adina E. Pusok, Dave R. Stegman, and Madeleine Kerr
Solid Earth, 13, 1455–1473, https://doi.org/10.5194/se-13-1455-2022, https://doi.org/10.5194/se-13-1455-2022, 2022
Short summary
Short summary
Sediments play an important role in global volatile and tectonic cycles, yet their effect on subduction dynamics is poorly resolved. In this study, we investigate how sediment properties influence subduction dynamics and obtain accretionary or erosive-style margins. Results show that even a thin layer of sediments can exert a profound influence on the emergent regional-scale subduction dynamics.
Laure Chevalier and Harro Schmeling
Solid Earth, 13, 1045–1063, https://doi.org/10.5194/se-13-1045-2022, https://doi.org/10.5194/se-13-1045-2022, 2022
Short summary
Short summary
Fluid flow through rock occurs in many geological settings on different scales, at different temperature conditions and with different flow velocities. Fluid is either in local thermal equilibrium with the host rock or not. We explore the parameters of porous flow and give scaling laws. These allow us to decide whether porous flows are in thermal equilibrium or not. Applied to magmatic systems, moving melts in channels or dikes moderately to strongly deviate from thermal equilibrium.
Iris van Zelst, Fabio Crameri, Adina E. Pusok, Anne Glerum, Juliane Dannberg, and Cedric Thieulot
Solid Earth, 13, 583–637, https://doi.org/10.5194/se-13-583-2022, https://doi.org/10.5194/se-13-583-2022, 2022
Short summary
Short summary
Geodynamic modelling provides a powerful tool to investigate processes in the Earth’s crust, mantle, and core that are not directly observable. In this review, we present a comprehensive yet concise overview of the modelling process with an emphasis on best practices. We also highlight synergies with related fields, such as seismology and geology. Hence, this review is the perfect starting point for anyone wishing to (re)gain a solid understanding of geodynamic modelling as a whole.
Igor Ognev, Jörg Ebbing, and Peter Haas
Solid Earth, 13, 431–448, https://doi.org/10.5194/se-13-431-2022, https://doi.org/10.5194/se-13-431-2022, 2022
Short summary
Short summary
We present a new 3D crustal model of Volgo–Uralia, an eastern segment of the East European craton. We built this model by processing the satellite gravity data and using prior crustal thickness estimation from regional seismic studies to constrain the results. The modelling revealed a high-density body on the top of the mantle and otherwise reflected the main known features of the Volgo–Uralian crustal architecture. We plan to use the obtained model for further geothermal analysis of the region.
Cedric Thieulot and Wolfgang Bangerth
Solid Earth, 13, 229–249, https://doi.org/10.5194/se-13-229-2022, https://doi.org/10.5194/se-13-229-2022, 2022
Short summary
Short summary
One of the main numerical methods to solve the mass, momentum, and energy conservation equations in geodynamics is the finite-element method. Four main types of elements have been used in the past decades in hundreds of publications. For the first time we compare results obtained with these four elements on a series of geodynamical benchmarks and applications and draw conclusions as to which are the best ones and which are to be preferably avoided.
Jean Furstoss, Carole Petit, Clément Ganino, Marc Bernacki, and Daniel Pino-Muñoz
Solid Earth, 12, 2369–2385, https://doi.org/10.5194/se-12-2369-2021, https://doi.org/10.5194/se-12-2369-2021, 2021
Short summary
Short summary
In the first part of this article, we present a new methodology that we have developed to model the deformation and the microstructural evolutions of olivine rocks, which make up the main part of the Earth upper mantle. In a second part, using this methodology we show that microstructural features such as small grain sizes and preferential grain orientations can localize strain at the same intensity and can act together to produce an even stronger strain localization.
Duan Li, Jinsong Du, Chao Chen, Qing Liang, and Shida Sun
Solid Earth Discuss., https://doi.org/10.5194/se-2021-117, https://doi.org/10.5194/se-2021-117, 2021
Revised manuscript not accepted
Short summary
Short summary
Oceanic magnetic anomalies are generally carried out using only few survey lines and thus there are many areas with data gaps. Traditional interpolation methods based on the morphological characteristics of data are not suitable for data with large gaps. The use of dual-layer equivalent-source techniques may improve the interpolation of magnetic anomaly fields in areas with sparse data which gives a good consideration to the extension of the magnetic lineation feature.
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.
Lorenzo G. Candioti, Thibault Duretz, Evangelos Moulas, and Stefan M. Schmalholz
Solid Earth, 12, 1749–1775, https://doi.org/10.5194/se-12-1749-2021, https://doi.org/10.5194/se-12-1749-2021, 2021
Short summary
Short summary
We quantify the relative importance of forces driving the dynamics of mountain building using two-dimensional computer simulations of long-term coupled lithosphere–upper-mantle deformation. Buoyancy forces can be as high as shear forces induced by far-field plate motion and should be considered when studying the formation of mountain ranges. The strength of rocks flooring the oceans and the density structure of the crust control deep rock cycling and the topographic elevation of orogens.
Olivier de Viron, Michel Van Camp, Alexia Grabkowiak, and Ana M. G. Ferreira
Solid Earth, 12, 1601–1634, https://doi.org/10.5194/se-12-1601-2021, https://doi.org/10.5194/se-12-1601-2021, 2021
Short summary
Short summary
As the travel time of seismic waves depends on the Earth's interior properties, seismic tomography uses it to infer the distribution of velocity anomalies, similarly to what is done in medical tomography. We propose analysing the outputs of those models using varimax principal component analysis, which results in a compressed objective representation of the model, helping analysis and comparison.
Janik Dohmen and Harro Schmeling
Solid Earth, 12, 1549–1561, https://doi.org/10.5194/se-12-1549-2021, https://doi.org/10.5194/se-12-1549-2021, 2021
Short summary
Short summary
In partially molten regions within the Earth, the melt is able to move separately from the surrounding rocks. This allows for the emergence of so-called solitary porosity waves, driven by compaction and decompaction due to the melt with higher buoyancy. Our numerical models can predict whether a partially molten region will ascend dominated by solitary waves or diapirism. Even in diapiris-dominated regions, solitary waves will build up and ascend as a swarm when the ascend time is long enough.
Xin Zhong, Marcin Dabrowski, and Bjørn Jamtveit
Solid Earth, 12, 817–833, https://doi.org/10.5194/se-12-817-2021, https://doi.org/10.5194/se-12-817-2021, 2021
Short summary
Short summary
Elastic thermobarometry is an useful tool to recover paleo-pressure and temperature. Here, we provide an analytical model based on the Eshelby solution to calculate the residual stress and strain preserved in a mineral inclusion exhumed from depth. The method applies to ellipsoidal, anisotropic inclusions in infinite isotropic hosts. A finite-element method is also used for a facet effect. Volumetrically averaged stress is shown to be a good proxy for the overall heterogeneous stress stage.
Maximilian Lowe, Jörg Ebbing, Amr El-Sharkawy, and Thomas Meier
Solid Earth, 12, 691–711, https://doi.org/10.5194/se-12-691-2021, https://doi.org/10.5194/se-12-691-2021, 2021
Short summary
Short summary
This study estimates the gravitational contribution from subcrustal density heterogeneities interpreted as subducting lithosphere beneath the Alps to the gravity field. We showed that those heterogeneities contribute up to 40 mGal of gravitational signal. Such density variations are often not accounted for in Alpine lithospheric models. We demonstrate that future studies should account for subcrustal density variations to provide a meaningful representation of the complex geodynamic Alpine area.
Antonio Manjón-Cabeza Córdoba and Maxim D. Ballmer
Solid Earth, 12, 613–632, https://doi.org/10.5194/se-12-613-2021, https://doi.org/10.5194/se-12-613-2021, 2021
Short summary
Short summary
The study of intraplate volcanism can inform us about underlying mantle dynamic processes and thermal and/or compositional anomalies. Here, we investigated numerical models of mantle flow and melting of edge-driven convection (EDC), a potential origin for intraplate volcanism. Our most important conclusion is that EDC can only produce moderate amounts of mantle melting. By itself, EDC is insufficient to support the formation of voluminous island-building volcanism over several millions of years.
Davide Tadiello and Carla Braitenberg
Solid Earth, 12, 539–561, https://doi.org/10.5194/se-12-539-2021, https://doi.org/10.5194/se-12-539-2021, 2021
Short summary
Short summary
We present an innovative approach to estimate a lithosphere density distribution model based on seismic tomography and gravity data. In the studied area, the model shows that magmatic events have increased density in the middle to lower crust, which explains the observed positive gravity anomaly. We interpret the densification through crustal intrusion and magmatic underplating. The proposed method has been tested in the Alps but can be applied to other geological contexts.
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.
Ángela María Gómez-García, Eline Le Breton, Magdalena Scheck-Wenderoth, Gaspar Monsalve, and Denis Anikiev
Solid Earth, 12, 275–298, https://doi.org/10.5194/se-12-275-2021, https://doi.org/10.5194/se-12-275-2021, 2021
Short summary
Short summary
The Earth’s crust beneath the Caribbean Sea formed at about 90 Ma due to large magmatic activity of a mantle plume, which brought molten material up from the deep Earth. By integrating diverse geophysical datasets, we image for the first time two fossil magmatic conduits beneath the Caribbean. The location of these conduits at 90 Ma does not correspond with the present-day Galápagos plume. Either this mantle plume migrated in time or these conduits were formed above another unknown plume.
Lorenzo G. Candioti, Stefan M. Schmalholz, and Thibault Duretz
Solid Earth, 11, 2327–2357, https://doi.org/10.5194/se-11-2327-2020, https://doi.org/10.5194/se-11-2327-2020, 2020
Short summary
Short summary
With computer simulations, we study the interplay between thermo-mechanical processes in the lithosphere and the underlying upper mantle during a long-term (> 100 Myr) tectonic cycle of extension–cooling–convergence. The intensity of mantle convection is important for (i) subduction initiation, (ii) the development of single- or double-slab subduction zones, and (iii) the forces necessary to initiate subduction. Our models are applicable to the opening and closure of the western Alpine Tethys.
Patrick Sanan, Dave A. May, Matthias Bollhöfer, and Olaf Schenk
Solid Earth, 11, 2031–2045, https://doi.org/10.5194/se-11-2031-2020, https://doi.org/10.5194/se-11-2031-2020, 2020
Short summary
Short summary
Mantle and lithospheric dynamics, elasticity, subsurface flow, and other fields involve solving indefinite linear systems. Tools include direct solvers (robust, easy to use, expensive) and advanced iterative solvers (complex, problem-sensitive). We show that a third option, ILDL preconditioners, requires less memory than direct solvers but is easy to use, as applied to 3D problems with parameter jumps. With included software, we hope to allow researchers to solve previously infeasible problems.
Philipp Eichheimer, Marcel Thielmann, Wakana Fujita, Gregor J. Golabek, Michihiko Nakamura, Satoshi Okumura, Takayuki Nakatani, and Maximilian O. Kottwitz
Solid Earth, 11, 1079–1095, https://doi.org/10.5194/se-11-1079-2020, https://doi.org/10.5194/se-11-1079-2020, 2020
Short summary
Short summary
To describe permeability, a key parameter controlling fluid flows in the Earth’s subsurface, an accurate determination of permeability on the pore scale is necessary. For this reason, we sinter artificial glass bead samples with various
porosities, determining the microstructure and permeability using both
experimental and numerical approaches. Based on this we provide
parameterizations of permeability, which can be used as input parameters for
large-scale numerical models.
Mark D. Lindsay, Sandra Occhipinti, Crystal Laflamme, Alan Aitken, and Lara Ramos
Solid Earth, 11, 1053–1077, https://doi.org/10.5194/se-11-1053-2020, https://doi.org/10.5194/se-11-1053-2020, 2020
Short summary
Short summary
Integrated interpretation of multiple datasets is a key skill required for better understanding the composition and configuration of the Earth's crust. Geophysical and 3D geological modelling are used here to aid the interpretation process in investigating anomalous and cryptic geophysical signatures which suggest a more complex structure and history of a Palaeoproterozoic basin in Western Australia.
Tomáš Fischer, Josef Vlček, and Martin Lanzendörfer
Solid Earth, 11, 983–998, https://doi.org/10.5194/se-11-983-2020, https://doi.org/10.5194/se-11-983-2020, 2020
Short summary
Short summary
Data on CO2 degassing help understanding the relations of the gas flow on geodynamic processes. Long-term gas flow measurements in rough field conditions present a challenge due to technical problems. We describe methods used for CO2 flow monitoring in West-Bohemia/Vogtland, which is typical for high CO2 flow, and present a new robust method based on pressure measurements in a water column. The results of 10 years of CO2 flow measurements and their relation to seismic activity are discussed.
Jana Schierjott, Antoine Rozel, and Paul Tackley
Solid Earth, 11, 959–982, https://doi.org/10.5194/se-11-959-2020, https://doi.org/10.5194/se-11-959-2020, 2020
Short summary
Short summary
We investigate the size of mineral grains of Earth's rocks in computer models of the whole Earth. This is relevant because grain size affects the stiffness (large grains are stiffer) and deformation of the Earth's mantle. We see that mineral grains grow inside stable non-deforming regions of the Earth. However, these regions are less stiff than expected. On the other hand, we find that grain size diminishes during deformation events such as when surface material comes down into the Earth.
Irene Pérez-Cáceres, David Jesús Martínez Poyatos, Olivier Vidal, Olivier Beyssac, Fernando Nieto, José Fernando Simancas, Antonio Azor, and Franck Bourdelle
Solid Earth, 11, 469–488, https://doi.org/10.5194/se-11-469-2020, https://doi.org/10.5194/se-11-469-2020, 2020
Short summary
Short summary
The metamorphism of the Pulo do Lobo unit (SW Iberian Massif) is described in this paper. To this end, three different and complementary methodologies have been applied. The new results reported here contribute to the knowledge of the metamorphic conditions of the Pulo do Lobo unit in relation to its deformation. Furthermore, the results are compared in order to assess the reliability of the different methods applied.
Ömer F. Bodur and Patrice F. Rey
Solid Earth, 10, 2167–2178, https://doi.org/10.5194/se-10-2167-2019, https://doi.org/10.5194/se-10-2167-2019, 2019
Short summary
Short summary
Convection in the deep Earth dynamically changes the elevation of plates. Amplitudes of those vertical motions predicted from numerical models are significantly higher than observations. We find that at small wavelengths (< 1000 km) this misfit can be due to the oversimplification in viscosity of rocks. By a suite of numerical experiments, we show that considering the non–Newtonian rheology of the mantle results in predictions in amplitude of dynamic topography consistent with observations.
Janik Dohmen, Harro Schmeling, and Jan Philipp Kruse
Solid Earth, 10, 2103–2113, https://doi.org/10.5194/se-10-2103-2019, https://doi.org/10.5194/se-10-2103-2019, 2019
Short summary
Short summary
In source regions of magmatic systems the temperature is above solidus and melt ascent is assumed to occur predominantly by two-phase flow. This two-phase flow allows for the emergence of solitary porosity waves. By now most solutions of these waves used strongly simplified viscosity laws, while in our laws the viscosity decreases rapidly for small melt fractions. The results show that for higher background porosities the phase velocities and the width of the wave are significantly decreased.
Cameron Spooner, Magdalena Scheck-Wenderoth, Hans-Jürgen Götze, Jörg Ebbing, György Hetényi, and the AlpArray Working Group
Solid Earth, 10, 2073–2088, https://doi.org/10.5194/se-10-2073-2019, https://doi.org/10.5194/se-10-2073-2019, 2019
Short summary
Short summary
By utilising both the observed gravity field of the Alps and their forelands and indications from deep seismic surveys, we were able to produce a 3-D structural model of the region that indicates the distribution of densities within the lithosphere. We found that the present-day Adriatic crust is both thinner and denser than the European crust and that the properties of Alpine crust are strongly linked to their provenance.
Philipp Eichheimer, Marcel Thielmann, Anton Popov, Gregor J. Golabek, Wakana Fujita, Maximilian O. Kottwitz, and Boris J. P. Kaus
Solid Earth, 10, 1717–1731, https://doi.org/10.5194/se-10-1717-2019, https://doi.org/10.5194/se-10-1717-2019, 2019
Short summary
Short summary
Prediction of rock permeability is of crucial importance for several research areas in geoscience. In this study, we enhance the finite difference code LaMEM to compute fluid flow on the pore scale using Newtonian and non-Newtonian rheologies. The accuracy of the code is demonstrated using several analytical solutions as well as experimental data. Our results show good agreement with analytical solutions and recent numerical studies.
Cited articles
Adalsteinsson, D., and Sethian, J.: The Fast Construction of Extension Velocities in Level Set Methods, J. Comput. Phys., 148, 2–22, 1999.
Andrews, E. and Billen, M.: Rheologic controls on the dynamics of slab detachment, Tectonophysics, 464, 60–69, 2009.
Androvičová, A., Čižková, H., and van den Berg, A.: The effects of rheological decoupling on slab deformation in the Earth's upper mantle, Stud. Geophys. Geod., 57, 460–481, 2013.
Bangerth, W. and Heister, T.: ASPECT: Advanced Solver for Problems in Earth's ConvecTion, Texas A&M University/Computational Infrastructure in Geodynamics, 2013.
Baumann, C., Gerya, T., and Connolly, J.: Numerical modelling of spontaneous slab breakoff dynamics during continental collision, Geological Society, London, Spec. Publicat., 332, 99–114, 2010.
Běhounková, M. and Čižková, H.: Long-wavelength character of subducted slabs in the lower mantle, Earth Planet. Sc. Lett., 275, 43–53, 2008.
Billen, M. and Hirth, G.: Rheologic controls on slab dynamics, Geochem. Geophy. Geosy., 8, Q08012, https://doi.org/10.1029/2007GC001597, 2007.
Bourgouin, L., Mühlhaus, H.-B., Hale, A., and Arsac, A.: Towards realistic simulations of lava dome growth using the level set method, Acta Geotech. Slov., 1, 225–236, 2006.
Bourgouin, L., Mühlhaus, H.-B., Hale, A., and Arsac, A.: Studying the influence of a solid shell on lava dome growth and evolution using the level set method, Geophys. J. Int., 170, 1431–1438, 2007.
Braun, J., Thieulot, C., Fullsack, P., DeKool, M., Beaumont, C., and Huismans, R.: DOUAR: a new three-dimensional creeping flow numerical model for the solution of geological problems, Phys. Earth Planet. In., 171, 76–91, 2008.
Brooks, A. and Hughes, T.: Stream-line upwind/Petrov–Galerkin formulation for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations, Comput. Methods Appl. Mech. Engrg., 32, 199–259, 1982.
Chertova, M. V., Geenen, T., van den Berg, A., and Spakman, W.: Using open sidewalls for modelling self-consistent lithosphere subduction dynamics, Solid Earth, 3, 313–326, https://doi.org/10.5194/se-3-313-2012, 2012.
Christensen, U.: The influence of trench migration on slab penetration into the lower mantle, Earth Planet. Sc. Lett., 140, 27–39, 1996.
Chopp, D.: Another look at velocity extensions in the level set method, SIAM J. Sci. Comput., 31, 3255–3273, 2009.
Čížková, H., van Hunen, J., and van den Berg, A.: Stress distribution within subducting slabs and their deformation in the transition zone, Phys. Earth Planet. In., 161, 202–214, 2007.
Crameri, F., Schmeling, H., Golabek, G., Duretz, T., Orendt, R., Buiter, S., May, D. A., Kaus, B., Gerya, T., and Tackley, P.: A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the "sticky air" method, Geophys. J. Int., 189, 38–54, 2012.
Duan, X., Ma, Y., and Zhang, R.: Optimal shape control of fluid flow using variational level set method, Phys. Lett. A, 372, 1374–1379, 2008.
Duretz, T., Gerya, T., and May, D.: Numerical modelling of spontaneous slab breakoff and subsequent topographic response, Tectonophysics, 502, 244–256, 2011.
Duretz, T., Schmalholz, S., and Gerya, T.: Dynamics of slab detachment, Geochem. Geophy. Geosy., 13, Q03020, https://doi.org/10.1029/2011GC004024, 2012.
Duretz, T., Gerya, T. V., and Spakman, W.: Slab detachment in laterally varying subduction zones: 3-D numerical modeling, Geophys. Res. Lett., 41, 1951–1956, 2014.
Enright, D., Fedkiw, R., Ferziger, J., and Mitchell, I.: A Hybrid Particle Level Set Method for Improved Interface Capturing, J. Comput. Phys., 183, 83–116, 2002.
Fullsack, P.: An arbitrary Lagrangian-Eulerian formulation for creeping flows and its application in tectonic models, Geophys. J. Int., 120, 1–23, 1995.
Fedkiw, R., Aslam, T., Merriman, B., and Osher, S.: A non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method), J. Comput. Phys., 152, 457–492, 1999.
Gerya, T., Fossati, D., Cantieni, C., and Seward, D.: Dynamic effects of aseismic ridge subduction: numerical modelling, Eur. J. Mineral., 21, 649–661, 2009.
Gottlieb, S. and Shu, C.: Total variation diminishing Runge–Kutta schemes, Math. Comput., 67, 73–85, 1998.
Gross, L., Bourgouin, L., Hale, A., and Mühlhaus, H.-B.: Interface modeling in incompressible media using level sets in Escript, Phys. Earth Planet. In., 163, 23–34, 2007.
Gurnis, M. and Hager, B.: Controls of the structure of subducted slabs, Nature, 335, 317–321, 1988.
Hale, A., Bourgouin, L., and Mühlhaus, H.: Using the level set method to model endogenous lava dome growth, J. Geophys. Res., 112, B03213, https://doi.org/10.1029/2006JB004445, 2007.
Hale, A., Gottschaldt, K., Rosenbaum, G., Bourgouin, L., Bauchy, M., and Mühlhaus, H.: Dynamics of slab tear faults: insights from numerical modelling, Tectonophysics, 483, 58–70, 2010.
Jiang, G. and Peng, D.: Weighted ENO schemes for Hamilton–Jacobi equations, SIAM J. Sci. Comput., 21, 2126–2143, 2000.
Kronbichler, M., Heister, T., and Bangerth, W.: High accuracy mantle convection simulation through modern numerical methods, Geophys. J. Int., 191, 12–29, 2012.
Lenardic, A. and Kaula, W.: A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces, J. Geophys. Res., 98, 8243–8260, 1993.
Magni, V., van Hunen, J., Funiciello, F., and Faccenna, C.: Numerical models of slab migration in continental collision zones, Solid Earth, 3, 293–306, https://doi.org/10.5194/se-3-293-2012, 2012.
Min, C.: On reinitializing level set functions, J. Comput. Phys., 229, 2764–2772, 2010.
Oka, H. and Ishii, K.: Numerical analysis on the motion of gas bubbles using level set method, J. Phys. Soc. Jpn., 68, 823–832, 1999.
Osher, S. and Fedkiw, R.: Level set methods: an overview and some recent results, J. Comput. Phys., 169, 463–502, 2001.
Osher, S. and Sethian, J.: Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations, J. Comput. Phys., 79, 12–49, 1988.
Osher, S. and Shu, C.-W.: High-order essentially non-oscillatory schemes for Hamilton–Jacobi equations, SIAM J. Numer. Anal., 28, 907–922, 1991.
Quinquis, M., Buiter, S., and Ellis, S.: The role of boundary conditions in numerical models of subduction zone dynamics, Tectonophysics, 497, 57–70, 2011.
Rao, R., Mondy, L., Noble, D., Moffat, H., D. B., A., and Notz, P.: A level set method to study foam processing: a validation study, Int. J. Numer. Meth. Fl., 68, 1362–1392, 2011.
Samuel, H. and Evonuk, M.: Modeling advection in geophysical flows with particle level sets, Geochem. Geophy. Geosy., 11, Q08020, https://doi.org/10.1029/2010GC003081, 2010.
Schmalholz, S.: A simple analytical solution for slab detachment, Earth Planet. Sc. Lett., 304, 45–54, 2011.
Schmeling, H., Babeyko, A., Enns, A., Faccenna, C., Funiciello, F., Gerya, T., Golabek, G., Grigull, S., Kaus, B., and Morra, G.: A benchmark comparison of spontaneous subduction models – towards a free surface, Phys. Earth Planet. In., 171, 198–223, 2008.
Segal, A. and Praagman, N.: The SEPRAN package. Technical report, Technical, Ingenieurs-Bureau Sepra, the Netherlands, available at: http://ta.twi.tudelft.nl/sepran/sepran.html, 2005.
Sethian, J.: Evolution, Implementation and Application of Level Set and Fast Marching Methods for Advancing Fronts, J. Comput. Phys., 169, 503–555, 2001.
Suckale, J., Hager, B., Elkins-Tanton, L., and Nave, J.: It takes three to tango: 2. Bubble dynamics in basaltic volcanoes and ramifications for modeling normal Strombolian activity, J. Geophys. Res., 115, B07410, https://doi.org/10.1029/2009JB006917, 2010.
Sussman, M., Fatemi, E., Smereka, P., and Osher, S.: A level set approach for computing solutions to incompressible two-phase flow II, in: Proceedings of the Sixth International Symposium on Computational Fluid Dynamics, Lake Tahoe, NV, 1995.
Tackley, P. and King, S.: Testing the tracer ratio method for modeling active compositional fields in mantle convection simulations, Geochem. Geophy. Geosy., 4, 8302, https://doi.org/10.1029/2001GC000214, 2003.
Thieulot, C.: FANTOM: two- and three-dimensional numerical modelling of creeping flows for the solution of geological problems, Phys. Earth Planet. In., 188, 47–68, 2011.
Thieulot, C.: ELEFANT: a user-friendly multipurpose geodynamics code, Solid Earth Discuss., 6, 1949–2096, https://doi.org/10.5194/sed-6-1949-2014, 2014.
Thieulot, C., Glerum, A., Hillebrand, B., Spakman, W., and Torsvik, T.: Multiphase geodynamical modelling using Aspect, in preparation, 2014a.
Thieulot, C., Schmalholz, S., Glerum, A., Hillebrand, B., and Spakman, W.: A two- and three-dimensional numerical comparison study of slab detachment, in preparation, 2014b.
van Hunen, J. and Allen, M.: Continental collision and slab break-off: a comparison of 3-D numerical models with observations, Earth Planet. Sc. Lett., 302, 21–37, 2011.
van Hunen, J. and van den Berg, A.: Plate tectonics on the early Earth: limitations imposed by strength and buoyancy of subducted lithosphere, Lithos, 103, 217–235, 2008.
Van Keken, P., King, S., Schmeling, H., Christensen, U., Neumeister, D., and Doin, M.: A comparison of methods for the modeling of thermochemical convection, J. Geophys. Res.-Sol. Ea., 102, 22477–22495, 1997.
van der Pijl, S., Segal, A., Vuik, C., and Wesseling, P.: Computing three-dimensional two-phase flows with a mass-conserving level set method, Comput. Visual, Sci., 11, 221–235, 2008.
Zlotnik, S., Fernandez, M., Diez, P., and Verges, J.: Modelling gravitational instabilities: slab break–off and Rayleigh–Taylor diapirism, Pure Appl. Geophys., 165, 1491–1510, 2008.
Short summary
Our paper demonstrates that the level set method is a viable method for material tracking in multi-material flow models. The different benchmarks illustate several advantages that the level set method provides over tracer-based methods. We therefore conclude that the level set method is well suited for geodynamical modeling.
Our paper demonstrates that the level set method is a viable method for material tracking in...
