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
Muhammad Anees, David Hindle, Ernesto Meneses Rioseco, Jonas Kley, Bernd Leiss, Mumtaz Muhammad Shah, and Javed Akhter Qureshi
EGUsphere, https://doi.org/10.5194/egusphere-2025-5252, https://doi.org/10.5194/egusphere-2025-5252, 2025
This preprint is open for discussion and under review for Solid Earth (SE).
Short summary
Short summary
We studied how heat is distributed underground in the western Himalaya to understand its potential for geothermal energy. Using computer models, we show that heat from radioactive rocks, uplift of mountains, and the shape of the land all influence temperatures in the crust. Deep valleys such as the Indus and Hunza may host accessible hot zones, meaning they could be good places to explore for geothermal energy.
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.
Muhammad Anees, David Hindle, Ernesto Meneses Rioseco, Jonas Kley, Bernd Leiss, Mumtaz Muhammad Shah, and Javed Akhter Qureshi
EGUsphere, https://doi.org/10.5194/egusphere-2025-5252, https://doi.org/10.5194/egusphere-2025-5252, 2025
This preprint is open for discussion and under review for Solid Earth (SE).
Short summary
Short summary
We studied how heat is distributed underground in the western Himalaya to understand its potential for geothermal energy. Using computer models, we show that heat from radioactive rocks, uplift of mountains, and the shape of the land all influence temperatures in the crust. Deep valleys such as the Indus and Hunza may host accessible hot zones, meaning they could be good places to explore for geothermal energy.
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.
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...