Topographic evolution of mountain belts controlled by rheology and surface process efficiency

Author(s):  
Sebastian G. Wolf ◽  
Ritske S. Huismans ◽  
Jean Braun ◽  
Xiaoping Yuan

<p>It has been a long-standing problem how mountain belts gain and loose topography during their tectonically active growth and inactive decay phase. It is widely recognized that mountain belt topography is generated by crustal shortening, and lowered by river bedrock erosion, linking climate to tectonics. However, it remains enigmatic how to reconcile high erosion rates in active orogens as observed in Taiwan or New Zealand, with long term survival of topography for 100s of Myrs as observed for example in the Uralides and Appalachians. Here we use for the first time a tight coupling between a landscape evolution model (FastScape) with an upper mantle scale tectonic (thermo-mechanical) model to investigate the different stages of mountain belt growth and decay. Using two end-member models, we demonstrate that growing orogens with high erosive power remain small (<200 km), reach steady state between tectonic in- and erosional material eff-flux, and are characterized by transverse valleys. Contrarily, mountain belts with medium to low erosive power will not reach growth steady state, grow wide, and are characterized by longitudinal rivers deflected by active thrusting. However, during growth both types of orogens reach the same height, controlled by rheology and independent of surface process efficiency. Erosional efficiency controls orogenic decay, which is counteracted by regional isostatic rebound. Rheological control of mountain height implies that there is a natural upper limit for the steepness index of rivers on Earth. To compare model results to various natural examples, we quantify the degree of longitudinal flow of modeled rivers with river “longitudinality” in several active or recently active orogens on Earth. Application of the river “longitudinality index” gives information whether (parts of) an orogen is or was at steady state during orogenic growth.</p>

2017 ◽  
Author(s):  
Jean Braun ◽  
Lorenzo Gemignani ◽  
Peter van der Beek

Abstract. The purpose of detrital thermochronology is to provide constraints on regional scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected in a system of river catchments, allows to extract information about the relative distribution of erosion rates in an eroding hinterland without relying on a steady-state assumption or the value of thermal parameters. The model is based on a relatively low number of parameters describing lithological variability among the various catchments and their sizes, and only uses the raw binned ages. In order to illustrate the method, we invert age distributions collected in the Eastern Himalaya, one of the most tectonically active places on Earth. From the inversion of the cooling age distributions we predict present day erosion rates of the catchments along the Siang-Tsangpo-Brahmaputra river system, as well as smaller tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i.e. from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The inversion additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in R.code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.


The Holocene ◽  
2017 ◽  
Vol 27 (12) ◽  
pp. 1914-1927 ◽  
Author(s):  
Reto Grischott ◽  
Florian Kober ◽  
Maarten Lupker ◽  
Juergen M Reitner ◽  
Ruth Drescher-Schneider ◽  
...  

Reconstructing paleo-denudation rates over Holocene timescales in an Alpine catchment provides a unique opportunity to isolate the climatic forcing of denudation from other tectonic or anthropogenic effects. Cosmogenic 10Be on two sediment cores from Lake Stappitz (Austrian Alps) were measured yielding a 15-kyr-long catchment-averaged denudation record of the upstream Seebach Valley. The persistence of a lake at the outlet of the valley fixed the baselevel, and the high mean elevation minimizes anthropogenic impacts. The 10Be record indicates a decrease in the proportion of paraglacial sediments from 15 to 7 kyr cal. BP after which the 10Be concentrations are considered to reflect hillslope erosion and thus can be converted to denudation rates. These ones significantly fluctuated over this time period: lower hillslope erosion rates of ca. 0.4 mm/year dated between 5 and 7 kyr cal. BP correlate with a stable climate, sparse flooding events and elevated temperatures that favoured the widespread growth of stabilizing soils and vegetation. Higher hillslope erosion rates of ca. 0.8 mm/year over the last ~4 kyr correlate with a variable, cooler climate where frequent flooding events enhance denudation of less protected hillslopes. Overall, our results suggest a tight coupling of climate and hillslope erosion in alpine landscapes as it has been observed in other parts of the Alps.


2019 ◽  
Vol 7 (1) ◽  
pp. 107-128 ◽  
Author(s):  
Odin Marc ◽  
Robert Behling ◽  
Christoff Andermann ◽  
Jens M. Turowski ◽  
Luc Illien ◽  
...  

Abstract. In active mountain belts with steep terrain, bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high-magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analysing earthquake- and monsoon-induced landslide inventories across different timescales. With time series of 5 m satellite images over four main valleys in central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured ∼5–8 times more landsliding than the pre-earthquake mean rate. The long-term size–frequency distribution of monsoon-induced landsliding (MIL) was derived from these inventories and from an inventory of landslides larger than ∼0.1 km2 that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size–frequency distribution for earthquake-induced landsliding (EQIL). These two distributions are dominated by infrequent, large and giant landslides but under-predict an estimated Holocene frequency of giant landslides (> 1 km3) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquakes of variable magnitude and when considering that a shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of ∼2±0.75 mm yr−1 in agreement with most thermo-chronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size–frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling timescale required to adequately capture the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic 10Be methods. This observation presents a strong caveat when interpreting spatial or temporal variability in erosion rates from this method. Thus, in areas where a very large, rare landslide contributes heavily to long-term erosion (as the Himalayas), we recommend 10Be sample in catchments with source areas > 10 000 km2 to reduce the method mean bias to below ∼20 % of the long-term erosion.


2018 ◽  
Author(s):  
Odin Marc ◽  
Robert Behling ◽  
Christoff Andermann ◽  
Jens M. Turowski ◽  
Luc Illien ◽  
...  

Abstract. In active mountain belts with steep terrain bedrock landsliding is a major erosional agent. In the Himalayas, landsliding is driven by annual hydro-meteorological forcing due to the summer monsoon and by rarer, exceptional events, such as earthquakes. Independent methods yield erosion rate estimates that appear to increase with sampling time, suggesting that rare, high magnitude erosion events dominate the erosional budget. Nevertheless, until now, neither the contribution of monsoon and earthquakes to landslide erosion, nor the proportion of erosion due to rare, giant landslides have been quantified in the Himalayas. We address these challenges by combining and analyzing earthquake and monsoon induced landslide inventories across different timescales. With time-series of 5 m satellite images over four main valleys in Central Nepal, we comprehensively mapped landslides caused by the monsoon from 2010 to 2018. We found no clear correlation between monsoon properties and landsliding, and a similar mean landsliding rate for all valleys, except in 2015, where the valleys affected by the earthquake featured ~ 5–8 times more landsliding than the pre-earthquake mean rate. The long-term size-frequency distribution of monsoon induced landslides (MIL) was derived from these inventories and from an inventory of landslides larger than ~ 0.1 km2 that occurred between 1972 and 2014. Using a published landslide inventory for the Gorkha 2015 earthquake, we derive the size-frequency distribution for earthquake-induced landslides (EQIL). These two distributions are dominated by infrequent, large and giant landslides, but underpredict an estimated Holocene frequency of giant landslides (> 1 km3) which we derived from a literature compilation. This discrepancy can be resolved when modelling the effect of a full distribution of earthquake of variable magnitude and considering that shallower earthquake may cause larger landslides. In this case, EQIL and MIL contribute about equally to a total long-term erosion of ~ 2 ± 0.75 mm.yr−1 in agreement with most thermochronological data. Independently of the specific total and relative erosion rates, the heavy-tailed size-frequency distribution from MIL and EQIL and the very large maximal landslide size in the Himalayas indicate that mean landslide erosion rates increase with sampling time, as has been observed for independent erosion estimates. Further, we find that the sampling time scale required for adequately capturing the frequency of the largest landslides, which is necessary for deriving long-term mean erosion rates, is often much longer than the averaging time of cosmogenic 10Be methods. This observation presents a strong caveat when interpreting spatial or temporal variability of erosion rates from this method.


2020 ◽  
Vol 6 (42) ◽  
pp. eaaz3166 ◽  
Author(s):  
B. A. Adams ◽  
K. X. Whipple ◽  
A. M. Forte ◽  
A. M. Heimsath ◽  
K. V. Hodges

The ongoing debate about the nature of coupling between climate and tectonics in mountain ranges derives, in part, from an imperfect understanding of how topography, climate, erosion, and rock uplift are interrelated. Here, we demonstrate that erosion rate is nonlinearly related to fluvial relief with a proportionality set by mean annual rainfall. These relationships can be quantified for tectonically active landscapes, and calculations based on them enable estimation of erosion where observations are lacking. Tests of the predictive power of this relationship in the Himalaya, where erosion is well constrained, affirm the value of our approach. Our model allows estimation of erosion rates in fluvial landscapes using readily available datasets, and the underlying relationship between erosion and rainfall offers the promise of a deeper understanding of how climate and tectonic evolution affect erosion and topography in space and time and of the potential influence of climate on tectonics.


2019 ◽  
Vol 7 (1) ◽  
pp. 275-299 ◽  
Author(s):  
Lorenz Michel ◽  
Christoph Glotzbach ◽  
Sarah Falkowski ◽  
Byron A. Adams ◽  
Todd A. Ehlers

Abstract. The Olympic Mountains of Washington state (USA) represent the aerially exposed accretionary wedge of the Cascadia Subduction Zone and are thought to be in flux steady state, whereby the mass outflux (denudation) and influx (tectonic accretion) into the mountain range are balanced. We use a multi-method approach to investigate how temporal variations in the influx and outflux could affect previous interpretations of flux steady state. This includes the analysis of published and new thermochronometric ages for (U–Th) ∕ He dating of apatite and zircon (AHe and ZHe, respectively), fission-track dating of apatite and zircon (AFT and ZFT, respectively), 1-D thermo-kinematic modeling of thermochronometric data, and independent estimates of outflux and influx. In total, we present 61 new AHe, ZHe, AFT, and ZFT thermochronometric ages from 21 new samples. AHe ages are generally young (< 4 Ma), and, in some samples, AFT ages (5–8 Ma) overlap ZHe ages (7–9 Ma) within uncertainties. Thermo-kinematic modeling shows that exhumation rates are temporally variable, with rates decreasing from > 2 to < 0.3 km Myr−1 around 5–7 Ma. With the onset of Plio–Pleistocene glaciation, exhumation rates increased to values > 1 km Myr−1. This demonstrates that the material outflux varies through time, requiring a commensurate variation in influx to maintain flux steady state. Evaluation of the offshore and onshore sediment record shows that the material influx is also variable through time and that the amount of accreted sediment in the wedge is spatially variable. This qualitatively suggests that significant perturbations of steady state occur on shorter timescales (105–106 years), like those created by Plio–Pleistocene glaciation. Our quantitative assessment of influx and outflux indicates that the Olympic Mountains could be in flux steady state on long timescales (107 years).


Sign in / Sign up

Export Citation Format

Share Document