Articles | Volume 5, issue 1
https://doi.org/10.5194/se-5-461-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-461-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
the Creative Commons Attribution 3.0 License.
Research article
| 
10 Jun 2014
Research article |
| 10 Jun 2014
Practical analytical solutions for benchmarking of 2-D and 3-D geodynamic Stokes problems with variable viscosity
I. Yu. Popov
St. Petersburg National Research University of Information Technologies, Mechanics and Optics, 49 Kronverkskiy, St. Petersburg, 197101, Russia
I. S. Lobanov
St. Petersburg National Research University of Information Technologies, Mechanics and Optics, 49 Kronverkskiy, St. Petersburg, 197101, Russia
S. I. Popov
St. Petersburg National Research University of Information Technologies, Mechanics and Optics, 49 Kronverkskiy, St. Petersburg, 197101, Russia
A. I. Popov
St. Petersburg National Research University of Information Technologies, Mechanics and Optics, 49 Kronverkskiy, St. Petersburg, 197101, Russia
T. V. Gerya
Institute of Geophysics, Department of Earth Sciences, Swiss Federal Institute of Technology Zurich (ETH), 5 Sonnegstrasse, 8092 Zurich, Switzerland
Related authors
No articles found.
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.
Simon Preuss, Jean Paul Ampuero, Taras Gerya, and Ylona van Dinther
Solid Earth, 11, 1333–1360, https://doi.org/10.5194/se-11-1333-2020, https://doi.org/10.5194/se-11-1333-2020, 2020
Short summary
Short summary
In this paper, we present newly developed numerical models to simulate episodic growth of geological faults.
This growth of faults occurs during the seismic cycle, with spontaneously generated primary and secondary fault structures. With these models we are able to show the evolution of complex fault geometries. Additionally, we can quantify the impact of earthquakes on fault growth.
Related subject area
Geodynamics
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
Fast uplift in the Southern Patagonian Andes due to long and short term deglaciation and the asthenospheric window underneath
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
3-D crustal density model of the Sea of Marmara
Oblique rifting: the rule, not the exception
GHOST: Geoscientific Hollow Sphere Tessellation
A high-resolution lithospheric magnetic field model over southern Africa based on a joint inversion of CHAMP, Swarm, WDMAM, and ground magnetic field data
Mechanical models to estimate the paleostress state from igneous intrusions
Density structure and isostasy of the lithosphere in Egypt and their relation to seismicity
The effect of obliquity on temperature in subduction zones: insights from 3-D numerical modeling
Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography
Nonlinear viscoplasticity in ASPECT: benchmarking and applications to subduction
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.
Veleda Astarte Paiva Muller, Pietro Sternai, and Christian Sue
EGUsphere, https://doi.org/10.5194/egusphere-2023-1711, https://doi.org/10.5194/egusphere-2023-1711, 2023
Short summary
Short summary
High uplift rates measured in the Southern Patagonian Icefield reflect 10 % of ice loss since the Little Ice Age (LIA), a remain of the ice sheet that covered the Patagonian Andes in the Last Glacial Maximum (LGM). Using numerical modelling we estimate that the mantle potential temperature due to a slab window needs to be 150–200 °C higher than in normal subduction zones to generate the measured uplift rates, and that both LIA and LGM deglaciation contribute to the modern uplift rate budget.
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.
Ershad Gholamrezaie, Magdalena Scheck-Wenderoth, Judith Bott, Oliver Heidbach, and Manfred R. Strecker
Solid Earth, 10, 785–807, https://doi.org/10.5194/se-10-785-2019, https://doi.org/10.5194/se-10-785-2019, 2019
Short summary
Short summary
Based on geophysical data integration and 3-D gravity modeling, we show that significant density heterogeneities are expressed as two large high-density bodies in the crust below the Sea of Marmara. The location of these bodies correlates spatially with the bends of the main Marmara fault, indicating that rheological contrasts in the crust may influence the fault kinematics. Our findings may have implications for seismic hazard and risk assessments in the Marmara region.
Sascha Brune, Simon E. Williams, and R. Dietmar Müller
Solid Earth, 9, 1187–1206, https://doi.org/10.5194/se-9-1187-2018, https://doi.org/10.5194/se-9-1187-2018, 2018
Short summary
Short summary
Fragmentation of continents often involves obliquely rifting segments that feature a complex three-dimensional structural evolution. Here we show that more than ~ 70 % of Earth’s rifted margins exceeded an obliquity of 20° demonstrating that oblique rifting should be considered the rule, not the exception. This highlights the importance of three-dimensional approaches in modelling, surveying, and interpretation of those rift segments where oblique rifting is the dominant mode of deformation.
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.
Foteini Vervelidou, Erwan Thébault, and Monika Korte
Solid Earth, 9, 897–910, https://doi.org/10.5194/se-9-897-2018, https://doi.org/10.5194/se-9-897-2018, 2018
Tara L. Stephens, Richard J. Walker, David Healy, Alodie Bubeck, and Richard W. England
Solid Earth, 9, 847–858, https://doi.org/10.5194/se-9-847-2018, https://doi.org/10.5194/se-9-847-2018, 2018
Short summary
Short summary
We present mechanical models that use the attitude and opening angles of igneous sills to constrain stress axes, the stress ratio, and relative magma pressure during dilation. The models can be applied to any set of dilated structures, including dikes, sills, or veins. Comparison with paleostress analysis for coeval faults and deformation bands indicates that sills can be used to characterise the paleostress state in areas where other brittle deformation structures (e.g. faults) are not present.
Mikhail K. Kaban, Sami El Khrepy, and Nassir Al-Arifi
Solid Earth, 9, 833–846, https://doi.org/10.5194/se-9-833-2018, https://doi.org/10.5194/se-9-833-2018, 2018
Short summary
Short summary
We present an integrative model of the crust and upper mantle of Egypt based on an analysis of gravity, seismic, and geological data. These results are essential for deciphering the link between the dynamic processes in the Earth system and near-surface processes (particularly earthquakes) that influence human habitat. We identified the distinct fragmentation of the lithosphere of Egypt in several blocks. This division is closely related to the seismicity patterns in this region.
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.
Anthony Osei Tutu, Bernhard Steinberger, Stephan V. Sobolev, Irina Rogozhina, and Anton A. Popov
Solid Earth, 9, 649–668, https://doi.org/10.5194/se-9-649-2018, https://doi.org/10.5194/se-9-649-2018, 2018
Short summary
Short summary
The Earth's surface is characterized by numerous geological processes, formed throughout the Earth's history to present day. The interior (mantle), on which plates rest, undergoes convection motion, generating stresses in the lithosphere plate and also causing the plate motion. This study shows that shallow density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting topography, giving the importance depth sampling.
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.
Cited articles
Albers, M.: A local mesh refinement multigrid method for 3-D convection problems with strongly variable viscosity, J. Comput. Phys., 160, 126–150, 2000.
Blankenbach, B., Busse, F., Christensen, U. Cserepes, L., Gunkel, D., Hansen, U., Harder, H., Jarvis, G., Koch, M., Marquart, G., Moore, D., Olson, P., Schmeling H., and Schnaubelt, T.: A benchmark comparison for mantle convection codes, Geophys. J. Int., 98, 23–38, 1989.
Buiter, S. J. H., Babeyko, A. Y., Ellis, S., Gerya, T. V., Kaus, B. J. P., Kellner, A., Schreurs, G., and Yamada, Y.: The numerical sandbox: comparison of model results for a shortening and an extension experiment, in: Analogue and Numerical Modelling of Crustal-Scale Processes, Geol. Soc., London, Special Publications, 253, 29–64, 2006.
Busse, F. H., Christensen, U., Clever, R., Cserepes, L., Gable, C., Giannandrea, E., Guillon, L., Houseman, G., Nataf, H.-C., Ogawa, M., Parmentier, M., Sotin, C., and Travis, B.: 3-D convection at infinite Prandtl number in Cartesian geometry – a benchmark comparison, Geophys. Astro. Fluid, 75, 39–59, 1993.
Crameri, F., Schmeling, H., Golabek1, 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, Geophys. J. Int., 189, 38–54, 2012.
Dabrowski, M., Krotkiewski, M., and Schmid, D. W.: MILAMIN: MATLAB-based finite-element method solver for large problems, Geochem. Geophy. Geosy., 9, Q04030, https://doi.org/10.1029/2007GC001719, 2008.
Deubelbeiss, Y. and Kaus, B. J. P.: Comparison of Eulerian and Lagrangian numerical techniques for the Stokes equations in the presence of strongly varying viscosity, Phys. Earth Planet. In., 171, 92–111, 2008.
Duretz, T., May, D. A., Gerya, T. V., and Tackley, P. J.: Discretization errors and free surface stabilization in the finite difference and marker-in-cell method for applied geodynamics: a numerical study, Geochem. Geophy. Geosy., 12, Q07004, https://doi.org/10.1029/2011GC003567, 1–26, 2011.
Gerya, T.: Introduction to Numerical Geodynamic Modelling, Cambridge University Press, Cambridge, 2010.
Gerya, T. V. and Yuen, D. A.: Characteristic-based marker-in-cell method with conservative finite-difference schemes for modeling geological flows with strongly variable transport properties, Phys. Earth Planet. In., 140, 293–318, 2003.
Gerya, T. V. and Yuen, D. A.: Robust characteristics method for modeling multiphase visco-elasto-plastic thermo-mechanical problems, Phys. Earth Planet. In., 163, 83–105, 2007.
Gerya, T. V., May, D. A., and Duretz, T.: An adaptive staggered grid finite difference method for modeling geodynamic Stokes flows with strongly variable viscosity, Geochem. Geophy. Geosy., 14, 1200–1225, 2013.
Hager, B. H. and O'Connell, R. J.: A simple global model of plane dynamics and mantle convection, J. Geophys. Res., 86, 4843–4867, 1981.
Ismail-Zadeh, A. and Tackley, P.: Computational Methods in Geodynamics, Cambridge University Press, Cambridge, 2010.
Kaus, B. J. P.: Factors that control the angle of shear bands in geodynamic numerical models of brittle deformation, Tectonophysics, 484, 36–47, 2010.
Kaus, B. J. P. and Becker, T. W.: Effects of elasticity on the Rayleigh–Taylor instability: implications for large-scale geodynamics, Geophys. J. Int., 168, 843–862, 2007.
Kaus, B. J. P. and Schmalholz, S. M.: 3-D Finite amplitude folding: implications for stress evolution during crustal and lithospheric deformation, Geophys. Res. Lett., 33, L14309, https://doi.org/10.1029/2011GC003567, 2006.
Lemiale, V., Muhlhaus, H.-B., Moresi, L. and Stafford, J.: Shear banding analysis of plastic models formulated for incompressible viscous flows, Phys. Earth Planetary Interiors, 171, 177–186, 2008.
Lobanov, I. S., Popov, I. Yu., Popov, A. I., and Gerya, T. V.: Numerical approach to the Stokes problem with high contrasts in viscosity, Appl. Mathemat. Comput., 235, 17–25, 2014.
Moresi, L., Zhong, S., and Gurnis, M.: The accuracy of finite element solutions of Stokes' flow with strongly varying viscosity, Phys. Earth Planet. In., 97, 83–94, 1996.
Moresi, L., Dufour, F., and Mulhaus, H.-B.: A Lagrangian integration point finite element method for large deformation modelingof viscoelastic geomaterials, J. Comput. Phys., 184, 476–497, 2003.
Popov, A. and Sobolev, S.: SLIM3D: a tool for three-dimensional thermomechanical modeling of lithospheric deformation with elasto-visco-plastic rheology, Phys. Earth Planet. In., 171, 55–75, 2008.
Popov, I. Yu. and Makeev, I. V.: A benchmark solution for 2-D Stokes flow over cavity, Z. Angew. Math. Phys., 65, 339–348, 2014.
Ramberg, H.: Instability of layered system in the field of gravity, II, Phys. Earth Planet. In., 1, 448–474, 1968.
Revenaugh, J. and Parsons, B.: Dynamic topography and gravity anomalies for fluid layers whose viscosity varies exponentially with depth, Geophys. J. Int., 90, 349–368, 1987.
Schmeling, H., Babeyko, A. Y., Enns, A., Faccenna, C., Funiciello, F., Gerya, T. V., Golabek, G. J., Grigull, S., Kaus, B. J. P., Morra, G., Schmalholz, S. M., and van Hunen, J.: A benchmark comparison of spontaneous subduction models – towards a free surface, Phys. Earth Planet. In., 171, 198–223, 2008.
Schmid, D. V. and Podladchikov, Yu. Yu.: Analytical solutions for deformable elliptical inclusions in general shear, Geophys. J. Int., 155, 269–288, 2003.
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, 2008.
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. In., 171, 7–18, 2008.
Tackley, P. J. and King, S. D.: Testing the tracer ratio methodfor modeling active compositional fields in mantle convection simulations, Geochem. Geophy. Geosys., 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 Plane. Int., 188, 47–68, 2011.
Torrance, K. E. and Turcotte, D. L.: Thermal convection with large viscosity variations, J. Fluid Mech., 47, 113–125, 1971.
Turcotte, D. L. and Schubert, G.: Geodynamics, Cambridge University Press, Cambridge, 2002.
van Keken, P. E., King, S. D., Schmeling, H., Christensen, U. R., Neumeister, D., and Doin, M.-P.: A comparison of methods for the modeling of thermochemical convection, J. Geophys. Res.-Sol. Earth., 102, 22477–22495, 1997.
van Keken, P. E., Currie, C., King, S. D., Behn, M. D., Cangleoncle, A., He, J., Katz, R. F., Lin, S.-C., Parmentier, E. M., Spiegelman, M., and Wang, K.: A community benchmark for subduction zone modeling, Phys. Earth Planet. In., 171, 187–197, 2008.
Zhong, S.: Analytic solutions for Stokes flow with lateral variations in viscosity, Geophys. J. Int., 124, 18–28, 1996.
Zhong, S. and Gurnis, M.: The role of plates and temperature-dependent viscosity in phase change dynamics, J. Geophys. Res., 99, 15903–15917, 1994.
Zhong, S., McNamara, A., Tan, E., Moresi, L., and Gurnis, M.: A benchmark study on mantle convection in a 3-D spherical shell using CitcomS, Geochem. Geophy. Geosy., 9, Q10017, https://doi.org/10.1029/2008GC002048, 2008.
