Articles | Volume 14, issue 2
https://doi.org/10.5194/se-14-197-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/se-14-197-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A corrected finite-difference scheme for the flexure equation with abrupt changes in coefficient
Department of Structural Geology and Geodynamics, Georg-August-Universität Göttingen, Goldschmidtstr. 3, 37077 Göttingen, Germany
Olivier Besson
Institut de mathématiques, Université de Neuchâtel, Rue Emile-Argand 11, 2000, Neuchâtel, Switzerland
Related authors
David Hindle and Jonas Kley
Solid Earth, 12, 2425–2438, https://doi.org/10.5194/se-12-2425-2021, https://doi.org/10.5194/se-12-2425-2021, 2021
Short summary
Short summary
Central western Europe underwent a strange episode of lithospheric deformation, resulting in a chain of small mountains that run almost west–east across the continent and that formed in the middle of a tectonic plate, not at its edges as is usually expected. Associated with these mountains, in particular the Harz in central Germany, are marine basins contemporaneous with the mountain growth. We explain how those basins came to be as a result of the mountains bending the adjacent plate.
David Hindle, Boris Sedov, Susanne Lindauer, and Kevin Mackey
Solid Earth, 10, 561–580, https://doi.org/10.5194/se-10-561-2019, https://doi.org/10.5194/se-10-561-2019, 2019
Short summary
Short summary
On one of the least studied boundaries between tectonic plates (North America–Okhotsk in northeastern Russia), which moves very similarly to the famous San Andreas fault in California, we have found the traces of earthquakes from the recent past, but before the time of historical records. This makes us a little more sure that the fault is still the place where movement between the plates takes place, and when it happens again, there could be dangerous earthquakes.
David Hindle and Jonas Kley
Solid Earth, 12, 2425–2438, https://doi.org/10.5194/se-12-2425-2021, https://doi.org/10.5194/se-12-2425-2021, 2021
Short summary
Short summary
Central western Europe underwent a strange episode of lithospheric deformation, resulting in a chain of small mountains that run almost west–east across the continent and that formed in the middle of a tectonic plate, not at its edges as is usually expected. Associated with these mountains, in particular the Harz in central Germany, are marine basins contemporaneous with the mountain growth. We explain how those basins came to be as a result of the mountains bending the adjacent plate.
David Hindle, Boris Sedov, Susanne Lindauer, and Kevin Mackey
Solid Earth, 10, 561–580, https://doi.org/10.5194/se-10-561-2019, https://doi.org/10.5194/se-10-561-2019, 2019
Short summary
Short summary
On one of the least studied boundaries between tectonic plates (North America–Okhotsk in northeastern Russia), which moves very similarly to the famous San Andreas fault in California, we have found the traces of earthquakes from the recent past, but before the time of historical records. This makes us a little more sure that the fault is still the place where movement between the plates takes place, and when it happens again, there could be dangerous earthquakes.
Related subject area
Subject area: Tectonic plate interactions, magma genesis, and lithosphere deformation at all scales | Editorial team: Geodynamics and quantitative modelling | Discipline: Geodynamics
How a volcanic arc influences back-arc extension: insight from 2D numerical models
Various lithospheric deformation patterns derived from rheological contrasts between continental terranes: insights from 2-D numerical simulations
The influence of viscous slab rheology on numerical models of subduction
Statistical appraisal of geothermal heat flow observations in the Arctic
Thrusts control the thermal maturity of accreted sediments
The role of continental lithospheric thermal structure in the evolution of orogenic systems: application to the Himalayan–Tibetan collision zone
The effect of temperature-dependent material properties on simple thermal models of subduction zones
Plume–ridge interactions: ridgeward versus plate-drag plume flow
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
A new finite element approach to model microscale strain localization within olivine aggregates
Buoyancy versus shear forces in building orogenic wedges
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.
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.
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.
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.
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.
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.
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.
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.
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.
Cited articles
Audet, P. and Mareschal, J.-C.: Variations in elastic thickness in the Canadian
Shield, Earth Planet. Sci. Lett., 226, 17–31,
https://doi.org/10.1016/j.epsl.2004.07.035, 2004. a
Axelsson, O.: Iterative Solution Methods, Cambridge University Press,
https://doi.org/10.1017/CBO9780511624100, 1994. a
Barrell, J.: The strength of the Earth's crust, J. Geol., 22,
655–683, https://doi.org/10.1086/622145, 1914. a, b
Beaumont, C.: Foreland basins, Geophys. J. Int., 65, 291–329,
https://doi.org/10.1111/j.1365-246X.1981.tb02715.x, 1981. a, b
Bleyer, J.: Numerical Tours of Computational Mechanics with FEniCS, Zenodo [code],
https://doi.org/10.5281/zenodo.1287832, 2018. a
Bodine, J., Steckler, M., and Watts, A.: Observations of flexure and the
rheology of the oceanic lithosphere, J. Geophys. Res.-Sol.
Ea., 86, 3695–3707, https://doi.org/10.1029/JB086iB05p03695, 1981. a
Braun, J., Deschamps, F., Rouby, D., and Dauteuil, O.: Flexure of the
lithosphere and the geodynamical evolution of non-cylindrical rifted passive
margins: Results from a numerical model incorporating variable elastic
thickness, surface processes and 3D thermal subsidence, Tectonophysics, 604,
72–82, https://doi.org/10.1016/j.tecto.2012.09.033, 2013. a
Burkhard, M. and Sommaruga, A.: Evolution of the western Swiss Molasse basin:
structural relations with the Alps and the Jura belt, Geological Society,
London, Special Publications, 134, 279–298,
https://doi.org/10.1144/GSL.SP.1998.134.01.13, 1998. a
Burov, E. and Watts, A.: The long-term strength of continental
lithosphere: “jelly sandwich” or “crème brûlée”?, GSA today, 16,
4, https://doi.org/10.1130/1052-5173(2006)016<4:TLTSOC>2.0.CO;2, 2006. a
DeCelles, P. G. and Giles, K. A.: Foreland basin systems, Basin Res., 8,
105–123, https://doi.org/10.1046/j.1365-2117.1996.01491.x, 1996. a
DLR – German Aerospace Center: TanDEM-X – Digital Elevation Model (DEM) –
Global, 12 m,
https://gdk.gdi-de.org/geonetwork/srv/api/records/5eecdf4c-de57-4624-99e9-60086b032aea,
(last access: 4 July 2021), 2018. a
Earth Resources Observation And Science (EROS) Center: Shuttle Radar
Topography Mission (SRTM) 1 Arc-Second Global, https://doi.org/10.5066/F7PR7TFT, 2017. a
Egan, S.: The flexural isostatic response of the lithosphere to extensional
tectonics, Tectonophysics, 202, 291–308, https://doi.org/10.1016/0040-1951(92)90115-M,
1992. a
Fisher, O.: Results of a Transcontinental series of gravity measurements,
Nature, 52, 433–435, https://doi.org/10.1038/052433a0, 1895. a
Gilbert, G.: The strength of the earth's crust, B. Geol.
Soc. Am., 1, 23–27, 1889. a
Gilbert, G. K.: Lake Bonneville, vol. 1, United States Geological Survey,
https://doi.org/10.3133/m1, 1890. a
Govers, R., Meijer, P., and Krijgsman, W.: Regional isostatic response to
Messinian Salinity Crisis events, Tectonophysics, 463, 109–129,
https://doi.org/10.1016/j.tecto.2008.09.026, 2009. a
Gunn, R.: A quantitative study of isobaric equilibrium and gravity anomalies in
the Hawaiian Islands, J. Frankl. Inst., 236, 373–390,
https://doi.org/10.1016/s0016-0032(43)90275-3, 1943a. a
Gunn, R.: Quantitative aspects of juxtaposed ocean deeps, mountain chains and
volcanic ranges, Geophysics, 12, 238–255, https://doi.org/10.1190/1.1437321,
1947. a
Hagen, E. S., Shuster, M. W., and Furlong, K. P.: Tectonic loading and
subsidence of intermontane basins: Wyoming foreland province, Geology, 13,
585–588, https://doi.org/10.1130/0091-7613(1985)13<585:TLASOI>2.0.CO;2, 1985. a, b
Hindle, D.: davidhindle /
flexure-1d-hs, GitHub repository [code],
https://github.com/davidhindle/flexure-1d-hs (last access: 30 March 2021), 2021. a
Hindle, D.: davidhindle/flexure-1d-hs: flexure-1d-hs_v1.0-beta (v1.1), Zenodo [code], https://doi.org/10.5281/zenodo.4643989, 2021. a
Hindle, D. and Kley, J.: The Subhercynian Basin: an example of an intraplate foreland basin due to a broken plate, Solid Earth, 12, 2425–2438, https://doi.org/10.5194/se-12-2425-2021, 2021. a
Horton, B. K. and DeCelles, P. G.: The modern foreland basin system adjacent to
the Central Andes, Geology, 25, 895–898,
https://doi.org/10.1130/0091-7613(1997)025<0895:TMFBSA>2.3.CO;2, 1997. a
Kaban, M., Chen, B., Tesauro, M., Petrunin, A., El Khrepy, S., and Al-Arifi,
N.: Reconsidering effective elastic thickness estimates by incorporating the
effect of sediments: A case study for Europe, Geophys. Res. Lett.,
45, 9523–9532, https://doi.org/10.1029/2018GL079732, 2018. a, b
Karner, G. and Watts, A.: Gravity anomalies and flexure of the lithosphere at
mountain ranges, J. Geophys. Res.-Sol. Ea., 88,
10449–10477, https://doi.org/10.1029/JB088iB12p10449, 1983. a
Langtangen, H. P.: Finite difference methods for wave motion,
http://hplgit.github.io/num-methods-for-PDEs/doc/pub/wave/pdf/wave-4print.pdf, last access:
3 November 2016. a
McKenzie, D. and Fairhead, D.: Estimates of the effective elastic thickness of
the continental lithosphere from Bouguer and free air gravity anomalies,
J. Geophys. Res.-Sol. Ea., 102, 27523–27552,
https://doi.org/10.1029/97JB02481, 1997. a
NASA/METI/AIST/Japan Spacesystems And U.S./Japan ASTER Science Team: ASTER
Global Digital Elevation Model, https://doi.org/10.5067/ASTER/ASTGTM.002, 2009. a
Pérez-Gussinyé, M. and Watts, A.: The long-term strength of Europe and
its implications for plate-forming processes, Nature, 436, 381–384,
https://doi.org/10.1038/nature03854, 2005. a
Sebben, S. and Baliga, B. R.: Some extensions of tridiagonal and pentadiagonal
matrix algorithms, Numerical Heat Transfer, Part B Fundamentals, 28,
323–351, https://doi.org/10.1080/10407799508928837, 1995. a
Thor, R.: Mapping the thickness of the Martian elastic lithosphere using
maximum likelihood estimation,
http://resolver.tudelft.nl/uuid:ca2a860f-a78e-4517-a829-eb03a7254f45, last access: 28 June 2016. a
Van der Meulen, M., Buiter, S., Meulenkamp, J., and Wortel, M.: An early
Pliocene uplift of the central Apenninic foredeep and its geodynamic
significance, Tectonics, 19, 300–313, https://doi.org/10.1029/1999TC900064, 2000. a, b
Walcott, R.: Flexural rigidity, thickness, and viscosity of the lithosphere,
J. Geophys. Res., 75, 3941–3954,
https://doi.org/10.1029/JB075i020p03941, 1970a. a, b, c
Walcott, R.: Flexure of the lithosphere at Hawaii, Tectonophysics, 9, 435–446,
https://doi.org/10.1016/0040-1951(70)90056-9, 1970b. a
Walcott, R.: Isostatic response to loading of the crust in Canada, Can.
J. Earth Sci., 7, 716–727, https://doi.org/10.1139/e70-070,
1970c. a
Walcott, R.: Late Quaternary vertical movements in eastern North America:
Quantitative evidence of glacio-isostatic rebound, Rev. Geophys., 10,
849–884, https://doi.org/10.1029/RG010i004p00849, 1972. a
Watts, A.: The effective elastic thickness of the lithosphere and the evolution
of foreland basins, Basin Res., 4, 169–178,
https://doi.org/10.1111/j.1365-2117.1992.tb00043.x, 1992. a, b
Wessel, P., Luis, J., Uieda, L., Scharroo, R., Wobbe, F., Smith, W., and Tian,
D.: The generic mapping tools version 6, Geochem. Geophys.
Geosyst., 20, 5556–5564, https://doi.org/10.1029/2019GC008515, 2019. a, b, c, d
Wickert, A. D.: Open-source modular solutions for flexural isostasy: gFlex v1.0, Geosci. Model Dev., 9, 997–1017, https://doi.org/10.5194/gmd-9-997-2016, 2016. a, b, c
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.
By making a change to the way we solve the flexure equation that describes how the Earth's outer...