Articles | Volume 16, issue 4
https://doi.org/10.5194/gmd-16-1315-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/gmd-16-1315-2023
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Modeling large‐scale landform evolution with a stream power law for glacial erosion (OpenLEM v37): benchmarking experiments against a more process-based description of ice flow (iSOSIA v3.4.3)
Department of Environment and Biodiversity, University of Salzburg, 5020 Salzburg, Austria
Jörg Robl
Department of Environment and Biodiversity, University of Salzburg, 5020 Salzburg, Austria
Stefan Hergarten
Institute of Earth and Environmental Sciences, University of Freiburg, 79104 Freiburg, Germany
David Lundbek Egholm
Department of Geoscience, Aarhus University, 8000 Aarhus, Denmark
Kurt Stüwe
Department of Earth Science, University of Graz, 8020 Graz, Austria
Related authors
No articles found.
Jörg Robl, Fabian Dremel, Kurt Stüwe, Stefan Hergarten, Christoph von Hagke, and Derek Fabel
Earth Surf. Dynam., 13, 745–770, https://doi.org/10.5194/esurf-13-745-2025, https://doi.org/10.5194/esurf-13-745-2025, 2025
Short summary
Short summary
The Bohemian Massif is one of several low mountain ranges in Europe that rises more than 1 km above the surrounding lowlands. Landscape characteristics indicate relief rejuvenation due to recent surface uplift. To constrain the pace of relief formation, we determined erosion rates of 20 catchments that range from 22 to 51 m Myr-1. Correlating these rates with topographic properties reveals that contrasts in bedrock erodibility represent a critical control of landscape evolution.
Stefan Hergarten
EGUsphere, https://doi.org/10.5194/egusphere-2025-2242, https://doi.org/10.5194/egusphere-2025-2242, 2025
This preprint is open for discussion and under review for Geoscientific Model Development (GMD).
Short summary
Short summary
Numerical glacier and ice-sheet models have been widely used in the context of climate change and landform evolution. While simulations of ice flow were numerically expensive for a long time, their performance has recently been boosted to an unprecedented level by machine learning techniques. This paper aims at keeping classical numerics competitive by introducing a novel numerical scheme, which allows for simulations at spatial resolutions of 25 m or even finer on standard desktop PCs.
Stefan Hergarten
Earth Surf. Dynam., 12, 1315–1327, https://doi.org/10.5194/esurf-12-1315-2024, https://doi.org/10.5194/esurf-12-1315-2024, 2024
Short summary
Short summary
Faceted topographies are impressive footprints of active tectonics in geomorphology. This paper investigates the evolution of faceted topographies at normal faults and their interaction with a river network theoretically and numerically. As a main result beyond several relations for the geometry of facets, the horizontal displacement associated with normal faults is crucial for the dissection of initially polygonal facets into triangular facets bounded by almost parallel rivers.
Stefan Hergarten
Earth Surf. Dynam., 12, 1193–1203, https://doi.org/10.5194/esurf-12-1193-2024, https://doi.org/10.5194/esurf-12-1193-2024, 2024
Short summary
Short summary
Toma hills are relatively isolated hills found in the deposits of rock avalanches, and their origin is still enigmatic. This paper presents the results of numerical simulations based on a modified version of a friction law that was originally introduced for snow avalanches. The model produces more or less isolated hills (which look much like toma hills) on the valley floor. The results provide, perhaps, the first explanation of the occurrence of toma hills based on a numerical model.
Stefan Hergarten
Geosci. Model Dev., 17, 781–794, https://doi.org/10.5194/gmd-17-781-2024, https://doi.org/10.5194/gmd-17-781-2024, 2024
Short summary
Short summary
The Voellmy rheology has been widely used for simulating snow and rock avalanches. Recently, a modified version of this rheology was proposed, which turned out to be able to predict the observed long runout of large rock avalanches theoretically. The software MinVoellmy presented here is the first numerical implementation of the modified rheology. It consists of MATLAB and Python classes, where simplicity and parsimony were the design goals.
Stefan Hergarten
Earth Surf. Dynam., 12, 219–229, https://doi.org/10.5194/esurf-12-219-2024, https://doi.org/10.5194/esurf-12-219-2024, 2024
Short summary
Short summary
Large landslides turn into an avalanche-like mode of flow at high velocities, which allows for a much longer runout than predicted for a sliding solid body. In this study, the Voellmy rheology widely used in models for hazard assessment is reinterpreted and extended. The new approach predicts the increase in runout length with volume observed in nature quite well and may thus be a major step towards a more consistent modeling of rock avalanches and improved hazard assessment.
Stefan Hergarten
Nat. Hazards Earth Syst. Sci., 23, 3051–3063, https://doi.org/10.5194/nhess-23-3051-2023, https://doi.org/10.5194/nhess-23-3051-2023, 2023
Short summary
Short summary
Rockslides are a major hazard in mountainous regions. In formerly glaciated regions, the disposition mainly arises from oversteepened topography and decreases through time. However, little is known about this decrease and thus about the present-day hazard of huge, potentially catastrophic rockslides. This paper presents a new theoretical framework that explains the decrease in maximum rockslide size through time and predicts the present-day frequency of large rockslides for the European Alps.
Jannick Strüven and Stefan Hergarten
Hydrol. Earth Syst. Sci., 27, 3041–3058, https://doi.org/10.5194/hess-27-3041-2023, https://doi.org/10.5194/hess-27-3041-2023, 2023
Short summary
Short summary
This study uses dendritic flow patterns to analyze the recession behavior of aquifer springs. The results show that the long-term recession becomes slower for large catchments. After a short recharge event, however, the short-term behavior differs strongly from the exponential recession that would be expected from a linear reservoir. The exponential component still accounts for more than 80 % of the total discharge, much more than typically assumed for karst aquifers.
Stefan Hergarten and Alexa Pietrek
Earth Surf. Dynam., 11, 741–755, https://doi.org/10.5194/esurf-11-741-2023, https://doi.org/10.5194/esurf-11-741-2023, 2023
Short summary
Short summary
The transition from hillslopes to channelized flow is typically attributed to a threshold catchment size in landform evolution models. Here we propose an alternative concept directly based on topography. Using this concept, channels and hillslopes self-organize, whereby the catchment size of the channel heads varies over some range. Our numerical results suggest that this concept works better than the established idea of a strict threshold catchment size.
Stefan Hergarten
Earth Surf. Dynam., 10, 671–686, https://doi.org/10.5194/esurf-10-671-2022, https://doi.org/10.5194/esurf-10-671-2022, 2022
Short summary
Short summary
Many studies on modeling landform evolution have focused on mountain ranges, while large parts of Earth's surface are quite flat and alluvial plains have been preferred locations for human settlements. Conducting large-scale simulations of fluvial erosion and sediment transport, this study reveals that rivers in a tectonically inactive foreland are much more dynamic than rivers in a mountain range; the local redistribution of deposits in the foreland is the main driver of the dynamics.
Stefan Hergarten and Jörg Robl
Geosci. Model Dev., 15, 2063–2084, https://doi.org/10.5194/gmd-15-2063-2022, https://doi.org/10.5194/gmd-15-2063-2022, 2022
Short summary
Short summary
The influence of climate on landform evolution has attracted great interest over the past decades. This paper presents a simple model for simulating the influence of topography on precipitation and the decrease in precipitation over large continental areas. The approach can be included in numerical models of large-scale landform evolution and causes only a moderate increase in the numerical complexity. It opens a door to investigating feedbacks between climate and landform evolution.
Stefan Hergarten
Earth Surf. Dynam., 9, 937–952, https://doi.org/10.5194/esurf-9-937-2021, https://doi.org/10.5194/esurf-9-937-2021, 2021
Short summary
Short summary
This paper presents a new approach to modeling glacial erosion on large scales. The formalism is similar to large-scale models of fluvial erosion, so glacial and fluvial processes can be easily combined. The model is simpler and numerically less demanding than established models based on a more detailed description of the ice flux. The numerical implementation almost achieves the efficiency of purely fluvial models, so that simulations over millions of years can be performed on standard PCs.
Anne-Laure Argentin, Jörg Robl, Günther Prasicek, Stefan Hergarten, Daniel Hölbling, Lorena Abad, and Zahra Dabiri
Nat. Hazards Earth Syst. Sci., 21, 1615–1637, https://doi.org/10.5194/nhess-21-1615-2021, https://doi.org/10.5194/nhess-21-1615-2021, 2021
Short summary
Short summary
This study relies on topography to simulate the origin and displacement of potentially river-blocking landslides. It highlights a continuous range of simulated landslide dams that go unnoticed in the field due to their small scale. The computation results show that landslide-dammed lake volume can be estimated from upstream drainage area and landslide volume, thus enabling an efficient hazard assessment of possible landslide-dammed lake volume – and flooding magnitude in case of dam failure.
Stefan Hergarten
Earth Surf. Dynam., 8, 841–854, https://doi.org/10.5194/esurf-8-841-2020, https://doi.org/10.5194/esurf-8-841-2020, 2020
Short summary
Short summary
Many contemporary models of large-scale fluvial erosion focus on the detachment-limited regime where all material entrained by the river is immediately excavated. This limitation facilitates the comparison with real river profiles and strongly reduces the numerical complexity. Here a simple formulation for the opposite case, transport-limited erosion, and a new numerical scheme that achieves almost the same numerical efficiency as detachment-limited models are presented.
Cited articles
Alley, R. B., Cuffey, K. M., and Zoet, L. K.: Glacial erosion: status and
outlook, Ann. Glaciol., 60, 1–13, https://doi.org/10.1017/aog.2019.38, 2019. a, b
Anderson, R. S., Molnar, P., and Kessler, M. A.: Features of glacial valley
profiles simply explained, J. Geophys. Res.-Earth,
111, F01004, https://doi.org/10.1029/2005jf000344, 2006. a
Bahr, D. B.: Width and length scaling of glaciers, J. Glaciol., 43,
557–562, https://doi.org/10.3189/S0022143000035164, 1997. a
Bernard, M., Steer, P., Gallagher, K., and Egholm, D. L.: The Impact of
Lithology on Fjord Morphology, Geophys. Res. Lett., 48,
e2021GL093101, https://doi.org/10.1029/2021GL093101, 2021. a, b, c
Braun, J., Zwartz, D., and Tomkin, J. H.: A new surface-processes model
combining glacial and fluvial erosion, Ann. Glaciol., 28, 282–290,
https://doi.org/10.3189/172756499781821797, 1999. a, b, c, d
Brocklehurst, S. H. and Whipple, K. X.: Hypsometry of glaciated landscapes,
Earth Surf. Proc. Land., 29, 907–926, https://doi.org/10.1002/esp.1083,
2004. a
Budd, W. F., Keage, P. L., and Blundy, N. A.: Empirical Studies of Ice Sliding,
J. Glaciol., 23, 157–170, https://doi.org/10.3189/S0022143000029804, 1979. a
Cook, S. J., Swift, D. A., Kirkbride, M. P., Knight, P. G., and Waller, R. I.:
The empirical basis for modelling glacial erosion rates, Nat. Commun., 11, 1–7, https://doi.org/10.1038/s41467-020-14583-8, 2020. a
Davy, P. and Lague, D.: Fluvial erosion/transport equation of landscape
evolution models revisited, J. Geophys. Res.-Earth, 114, F03007,
https://doi.org/10.1029/2008JF001146, 2009. a
Egholm, D. L., Nielsen, S. B., Pedersen, V. K., and Lesemann, J. E.: Glacial
effects limiting mountain height, Nature, 460, 884–887, https://doi.org/10.1038/nature08263, 2009. a
Egholm, D. L., Jansen, J. D., Braedstrup, C. F., Pedersen, V. K., Andersen,
J. L., Ugelvig, S. V., Larsen, N. K., and Knudsen, M. F.: Formation of
plateau landscapes on glaciated continental margins, Nat. Geosci., 10,
592–597, https://doi.org/10.1038/Ngeo2980, 2017. a, b, c, d
Fowler, A. C. and Larson, D. A.: On the flow of polythermal glaciers – I.
Model and preliminary analysis, Proc. R. Soc. Lond., 363, 217–242,
https://doi.org/10.1098/rspa.1978.0165, 1978. a
Hallet, B.: A theoretical model of glacial abrasion, J. Glaciol., 23,
39–50, https://doi.org/10.3189/S0022143000029725, 1979. a
Harbor, J. M.: Numerical modeling of the development of U-shaped valleys by
glacial erosion, Geol. Soc. Am. Bull., 104, 1364–1375,
https://doi.org/10.1130/0016-7606(1992)104<1364:Nmotdo>2.3.Co;2, 1992. a, b, c
Headley, R. M., Roe, G., and Hallet, B.: Glacier longitudinal profiles in
regions of active uplift, Earth Planet. Sci. Lett., 317, 354–362,
2012. a
Hergarten, S.: Transport-limited fluvial erosion – simple formulation and efficient numerical treatment, Earth Surf. Dynam., 8, 841–854, https://doi.org/10.5194/esurf-8-841-2020, 2020. a, b
Hergarten, S.: The influence of sediment transport on stationary and mobile
knickpoints in river profiles, J. Geophys. Res.-Earth, 126,
e2021JF006218, https://doi.org/10.1029/2021JF006218, 2021b. a, b, c
Hergarten, S.: Theoretical and numerical considerations of rivers in a tectonically inactive foreland, Earth Surf. Dynam., 10, 671–686, https://doi.org/10.5194/esurf-10-671-2022, 2022. a
Hergarten, S.: OpenLEM, OpenLEM [code], http://hergarten.at/openlem, last access: 20 February 2023. a
Hergarten, S. and Robl, J.: The linear feedback precipitation model (LFPM 1.0) – a simple and efficient model for orographic precipitation in the context of landform evolution modeling, Geosci. Model Dev., 15, 2063–2084, https://doi.org/10.5194/gmd-15-2063-2022, 2022. a, b
Hergarten, S., Wagner, T., and Stüwe, K.: Age and prematurity of the Alps
derived from topography, Earth Planet. Sci. Lett., 297, 453–460,
https://doi.org/10.1016/j.epsl.2010.06.048, 2010. a
Herman, F., Beyssac, O., Brughelli, M., Lane, S. N., Leprince, S., Adatte, T.,
Lin, J. Y. Y., Avouac, J. P., and Cox, S. C.: Erosion by an Alpine glacier,
Science, 350, 193–195, https://doi.org/10.1126/science.aab2386, 2015. a, b
Herman, F., Braun, J., Deal, E., and Prasicek, G.: The Response Time of Glacial
Erosion, J. Geophys. Res.-Earth, 123, 801–817,
https://doi.org/10.1002/2017JF004586, 2018. a
Howard, A. D.: A detachment-limited model for drainage basin evolution, Water
Resour. Res., 30, 2261–2285, https://doi.org/10.1029/94WR00757, 1994. a, b
Hutter, K.: Time-dependent surface elevation of an ice slope, J. Glaciology,
25, 247–266, https://doi.org/10.3189/S0022143000010479, 1980. a
Iverson, N. R.: A theory of glacial quarrying for landscape evolution models,
Geology, 40, 679–682, https://doi.org/10.1130/G33079.1, 2012. a
Kooi, H. and Beaumont, C.: Large-scale geomorphology: classical concepts
reconciled and integrated with contemporary ideas via a surface process
model, J. Geophys. Res., 101, 3361–3386, 1996. a
Koppes, M., Hallet, B., Rignot, E., Mouginot, J., Wellner, J. S., and Boldt,
K.: Observed latitudinal variations in erosion as a function of glacier
dynamics, Nature, 526, 100–103, https://doi.org/10.1038/nature15385, 2015. a, b
Kühni, A. and Pfiffner, O. A.: The relief of the Swiss Alps and adjacent
areas and its relation to lithology and structure: topographic analysis from
a 250-m DEM, Geomorphology, 41, 285–307, https://doi.org/10.1016/S0169-555x(01)00060-5, 2001. a
Lai, J. and Anders, A. M.: Climatic controls on mountain glacier basal thermal regimes dictate spatial patterns of glacial erosion, Earth Surf. Dynam., 9, 845–859, https://doi.org/10.5194/esurf-9-845-2021, 2021. a, b
Leith, K., Moore, J. R., Amann, F., and Loew, S.: Subglacial extensional
fracture development and implications for Alpine valley evolution, J. Geophys. Res.-Earth, 119, 62–81,
https://doi.org/10.1002/2012JF002691, 2014. a
Liebl, M.: C++ (OpenLEM, v37), C (iSOSIA, v3.4.3) and Matlab codes, Zenodo [code and data set], https://doi.org/10.5281/zenodo.6557805, 2022. a, b
Magrani, F., Valla, P. G., and Egholm, D.: Modelling alpine glacier geometry
and subglacial erosion patterns in response to contrasting climatic forcing,
Earth Surf. Proc. Land., 47, 1054–1072,
https://doi.org/10.1002/esp.5302, 2022. a, b, c
O'Callaghan, J. F. and Mark, D. M.: The extraction of drainage networks from
digital elevation data, Comput. Vision Graph., 28,
323–344, https://doi.org/10.1016/S0734-189X(84)80011-0, 1984. a
Pedersen, V. K. and Egholm, D. L.: Glaciations in response to climate
variations preconditioned by evolving topography, Nature, 493, 206–210,
https://doi.org/10.1038/nature11786, 2013. a, b, c
Pedersen, V. K., Egholm, D. L., and Nielsen, S. B.: Alpine glacial topography
and the rate of rock column uplift: a global perspective, Geomorphology, 122,
129–139, https://doi.org/10.1016/j.geomorph.2010.06.005, 2010. a, b
Penck, A.: Glacial features in the surface of the Alps, J. Geol., 13, 1–19, 1905. a
Prasicek, G., Larsen, I. J., and Montgomery, D. R.: Tectonic control on the
persistence of glacially sculpted topography, Nat. Commun., 6, 8028, https://doi.org/10.1038/ncomms9028, 2015. a
Prasicek, G., Hergarten, S., Deal, E., Herman, F., and Robl, J.: A glacial
buzzsaw effect generated by efficient erosion of temperate glaciers in a
steady state model, Earth Planet. Sci. Lett., 543, 116350,
https://doi.org/10.1016/j.epsl.2020.116350, 2020a. a
Prasicek, G., Hergarten, S., Deal, E., Herman, F., and Robl, J.: A glacial
buzzsaw effect generated by efficient erosion of temperate glaciers in a
steady state model, Earth Planet. Sci. Lett., 543, 116350,
https://doi.org/10.1016/j.epsl.2020.116350, 2020b. a, b
Robl, J., Prasicek, G., Hergarten, S., and Stüwe, K.: Alpine topography in
the light of tectonic uplift and glaciation, Glob. Planet. Change, 127,
34–49, https://doi.org/10.1016/j.gloplacha.2015.01.008, 2015. a, b, c, d
Robl, J., Hergarten, S., and Prasicek, G.: The topographic state of fluvially
conditioned mountain ranges, Earth Sci. Rev., 168, 290–317,
https://doi.org/10.1016/j.earscirev.2017.03.007, 2017. a
Robl, J., Hergarten, S., and Prasicek, G.: Glacial erosion promotes high
mountains on thin crust, Earth Planet. Sci. Lett., 538, 116196,
https://doi.org/10.1016/j.epsl.2020.116196, 2020. a
Roering, J. J., Kirchner, J. W., and Dietrich, W. E.: Evidence for nonlinear,
diffusive sediment transport on hillslopes and implications for landscape
morphology, Water Resour. Res., 35, 853–870,
https://doi.org/10.1029/1998wr900090, 1999. a
Royden, L. and Perron, J. T.: Solutions of the stream power equation and
application to the evolution of river longitudinal profiles, J. Geophys. Res.-Earth, 118, 497–518, https://doi.org/10.1002/jgrf.20031, 2013. a
Salcher, B., Prasicek, G., Baumann, S., and Kober, F.: Alpine relief limited
by glacial occupation time, Geology, 49, 1209–1213, https://doi.org/10.1130/G48639.1,
2021. a
Salcher, B. C., Kober, F., Kissling, E., and Willett, S. D.: Glacial impact on
short-wavelength topography and long-lasting effects on the denudation of a
deglaciated mountain range, Global Planet. Change, 115, 59–70,
https://doi.org/10.1016/j.gloplacha.2014.01.002, 2014. a
Schwanghart, W. and Kuhn, N. J.: TopoToolbox: A set of Matlab functions for
topographic analysis, Environ. Modell. Softw., 25, 770–781,
https://doi.org/10.1016/j.envsoft.2009.12.002, 2010. a
Schwanghart, W. and Scherler, D.: Short Communication: TopoToolbox 2 – MATLAB-based software for topographic analysis and modeling in Earth surface sciences, Earth Surf. Dynam., 2, 1–7, https://doi.org/10.5194/esurf-2-1-2014, 2014. a
Seguinot, J., Ivy-Ochs, S., Jouvet, G., Huss, M., Funk, M., and Preusser, F.: Modelling last glacial cycle ice dynamics in the Alps, The Cryosphere, 12, 3265–3285, https://doi.org/10.5194/tc-12-3265-2018, 2018. a, b
Sternai, P., Herman, F., Fox, M. R., and Castelltort, S.: Hypsometric analysis
to identify spatially variable glacial erosion, J. Geophys. Res.-Earth, 116, F03001, https://doi.org/10.1029/2010jf001823, 2011. a
Sternai, P., Herman, F., Valla, P. G., and Champagnac, J.-D.: Spatial and
temporal variations of glacial erosion in the Rhône valley (Swiss Alps):
Insights from numerical modeling, Earth Planet. Sci. Lett., 368,
119–131, https://doi.org/10.1016/j.epsl.2013.02.039, 2013. a
Weertman, J.: On the Sliding of Glaciers, J. Glaciol., 3, 33–38,
https://doi.org/10.3189/S0022143000024709, 1957. a
Whipple, K. X.: Bedrock rivers and the geomorphology of active orogens, Annu.
Rev. Earth Planet. Sci., 32, 151–185,
https://doi.org/10.1146/annurev.earth.32.101802.120356, 2004. a
Whipple, K. X. and Tucker, G. E.: Dynamics of the stream-power river incision
model: Implications for height limits of mountain ranges, landscape response
timescales, and research needs, J. Geophys. Res.-Sol. Ea.,
104, 17661–17674, https://doi.org/10.1029/1999JB900120, 1999. a
Whipple, K. X. and Tucker, G. E.: Implications of sediment-flux-dependent river
incision models for landscape evolution, J. Geophys. Res., 107, 2039,
https://doi.org/10.1029/2000JB000044, 2002. a
Wobus, C., Whipple, K. X., Kirby, E., Snyder, N., Johnson, J., Spyropolou, K.,
Crosby, B., and Sheehan, D.: Tectonics from topography: Procedures, promise,
and pitfalls, in: Tectonics, Climate, and Landscape Evolution, edited by:
Willett, S. D., Hovius, N., Brandon, M. T., and Fisher, D. M., vol. 398 of
GSA Special Papers, 55–74, Geological Society of America,
Boulder, Washington, D.C., https://doi.org/10.1130/2006.2398(04), 2006. a
Short summary
In this study, we benchmark a topography-based model for glacier erosion (OpenLEM) with a well-established process-based model (iSOSIA). Our experiments show that large-scale erosion patterns and particularly the transformation of valley length geometry from fluvial to glacial conditions are very similar in both models. This finding enables the application of OpenLEM to study the influence of climate and tectonics on glaciated mountains with reasonable computational effort on standard PCs.
In this study, we benchmark a topography-based model for glacier erosion (OpenLEM) with a...