@article{WolfHuismansBraunetal.2022, author = {Wolf, Sebastian G. and Huismans, Ritske S. and Braun, Jean and Yuan, Xiaoping}, title = {Topography of mountain belts controlled by rheology and surface processes}, series = {Nature : the international weekly journal of science}, volume = {606}, journal = {Nature : the international weekly journal of science}, number = {7914}, publisher = {Nature portfolio}, address = {Berlin}, issn = {0028-0836}, doi = {10.1038/s41586-022-04700-6}, pages = {516 -- 521}, year = {2022}, abstract = {It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics(1-4). However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process(5,6) and mantle-scale tectonic model(7). End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state(8), strength controlled (Bm approximate to 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.}, language = {en} } @article{NeuharthBruneWronaetal.2022, author = {Neuharth, Derek and Brune, Sascha and Wrona, Thilo and Glerum, Anne and Braun, Jean and Yuan, Xiaoping}, title = {Evolution of rift systems and their fault networks in response to surface processes}, series = {Tectonics}, volume = {41}, journal = {Tectonics}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0278-7407}, doi = {10.1029/2021TC007166}, pages = {22}, year = {2022}, abstract = {Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes (SP) code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: (a) distributed deformation and coalescence, (b) fault system growth, (c) fault system decline and basinward localization, (d) rift migration, and (e) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that SP are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.}, language = {en} } @article{NeuharthBruneGlerumetal.2021, author = {Neuharth, Derek and Brune, Sascha and Glerum, Anne and Morley, Chris K. and Yuan, Xiaoping and Braun, Jean}, title = {Flexural strike-slip basins}, series = {Geology : a venture in earth science reporting / the Geological Society of America}, volume = {50}, journal = {Geology : a venture in earth science reporting / the Geological Society of America}, number = {3}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {0091-7613}, doi = {10.1130/G49351.1}, pages = {361 -- 365}, year = {2021}, abstract = {Strike-slip faults are classically associated with pull-apart basins where continental crust is thinned between two laterally offset fault segments. We propose a subsidence mechanism to explain the formation of a new type of basin where no substantial segment offset or synstrike-slip thinning is observed. Such "flexural strike-slip basins" form due to a sediment load creating accommodation space by bending the lithosphere. We use a two-way coupling between the geodynamic code ASPECT and surface-processes code FastScape to show that flexural strike-slip basins emerge if sediment is deposited on thin lithosphere close to a strike slip fault. These conditions were met at the Andaman Basin Central fault (Andaman Sea, Indian Ocean), where seismic reflection data provide evidence of a laterally extensive flexural basin with a depocenter located parallel to the strike-slip fault trace.}, language = {en} } @article{Braun2022, author = {Braun, Jean}, title = {Comparing the transport-limited and ξ-q models for sediment transport}, series = {Earth surface dynamics}, volume = {10}, journal = {Earth surface dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-10-301-2022}, pages = {301 -- 327}, year = {2022}, abstract = {Here I present a comparison between two of the most widely used reduced-complexity models for the representation of sediment transport and deposition processes, namely the transport-limited (or TL) model and the under-capacity (or xi-q) model more recently developed by Davy and Lague (2009). Using both models, I investigate the behavior of a sedimentary continental system of length L fed by a fixed sedimentary flux from a catchment of size A(0) in a nearby active orogen through which sediments transit to a fixed base level representing a large river, a lake or an ocean. This comparison shows that the two models share the same steady-state solution, for which I derive a simple 1D analytical expression that reproduces the major features of such sedimentary systems: a steep fan that connects to a shallower alluvial plain. The resulting fan geometry obeys basic observational constraints on fan size and slope with respect to the upstream drainage area, A(0). The solution is strongly dependent on the size of the system, L, in comparison to a distance L-0, which is determined by the size of A(0), and gives rise to two fundamentally different types of sedimentary systems: a constrained system where L < L-0 and open systems where L > L-0. I derive simple expressions that show the dependence of the system response time on the system characteristics, such as its length, the size of the upstream catchment area, the amplitude of the incoming sedimentary flux and the respective rate parameters (diffusivity or erodibility) for each of the two models. I show that the xi-q model predicts longer response times. I demonstrate that although the manner in which signals propagates through the sedimentary system differs greatly between the two models, they both predict that perturbations that last longer than the response time of the system can be recorded in the stratigraphy of the sedimentary system and in particular of the fan. Interestingly, the xi-q model predicts that all perturbations in the incoming sedimentary flux will be transmitted through the system, whereas the TL model predicts that rapid perturbations cannot. I finally discuss why and under which conditions these differences are important and propose observational ways to determine which of the two models is most appropriate to represent natural systems.}, language = {en} } @article{BraunGemignanivanderBeek2018, author = {Braun, Jean and Gemignani, Lorenzo and van der Beek, Peter}, title = {Extracting information on the spatial variability in erosion rate stored in detrital cooling age distributions in river sands}, series = {Earth surface dynamics}, volume = {6}, journal = {Earth surface dynamics}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-6-257-2018}, pages = {257 -- 270}, year = {2018}, abstract = {One of the main purposes of detrital thermochronology is to provide constraints on the 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. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its 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 method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.}, language = {en} } @article{HermanBraunDealetal.2018, author = {Herman, Frederic and Braun, Jean and Deal, Eric and Prasicek, Gunther}, title = {The response time of glacial erosion}, series = {Journal of geophysical research : Earth surface}, volume = {123}, journal = {Journal of geophysical research : Earth surface}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1002/2017JF004586}, pages = {801 -- 817}, year = {2018}, abstract = {There has been recent progress in the understanding of the evolution of Quaternary climate. Simultaneously, there have been improvements in the understanding of glacial erosion processes, with better parameter constraints. Despite this, there remains much debate about whether or not the observed cooling over the Quaternary has driven an increase in glacial erosion rates. Most studies agree that the erosional response to climate change must be transient; therefore, the time scale of the climatic change and the response time of glacial erosion must be accounted for. Here we analyze the equations governing glacial erosion in a steadily uplifting landscape with variable climatic forcing and derive expressions for two fundamental response time scales. The first time scale describes the response of the glacier and the second one the glacial erosion response. We find that glaciers have characteristic time scales of the order of 10 to 10,000 years, while the characteristic time scale for glacial erosion is of the order of a few tens of thousands to a few million years. We then use a numerical model to validate the approximations made to derive the analytical solutions. The solutions show that short period forcing is dampened by the glacier response time, and long period forcing (>1 Myr) may be dampened by erosional response of glaciers when the rock uplift rates are high. In most tectonic and climatic conditions, we expect to see the strongest response of glacial erosion to periodic climatic forcing corresponding to Plio-Pleistocene climatic cycles. Finally, we use the numerical model to predict the response of glacial systems to the observed climatic forcing of the Quaternary, including, but not limited to, the Milankovich periods and the long-term secular cooling trend. We conclude that an increase of glacial erosion in response to Quaternary cooling is physically plausible, and we show that the magnitude of the increase depends on rock uplift and ice accumulation rates.}, language = {en} } @article{DealBraunBotter2018, author = {Deal, Eric and Braun, Jean and Botter, Gianluca}, title = {Understanding the role of rainfall and hydrology in determining fluvial erosion efficiency}, series = {Journal of geophysical research : Earth surface}, volume = {123}, journal = {Journal of geophysical research : Earth surface}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1002/2017JF004393}, pages = {744 -- 778}, year = {2018}, abstract = {Due to the challenges in upscaling daily climatic forcing to geological time, physically realistic models describing how rainfall drives fluvial erosion are lacking. To bridge this gap between short-term hydrology and long-term geomorphology, we derive a theoretical framework for long-term fluvial erosion rates driven by realistic climate by integrating an established stochastic-mechanistic model of hydrology into a threshold-stochastic formulation of stream power. The hydrological theory provides equations for the daily streamflow probability distribution as a function of climatic boundary conditions. The new parameters introduced are rooted firmly in established climatic and hydrological theory. This allows us to account for how fluvial erosion rates respond to changes in rainfall intensity, frequency, evapotranspiration rates, and soil moisture dynamics in a way that is consistent with existing theories. We use this framework to demonstrate how hydroclimatic conditions and erosion threshold magnitude control the degree of nonlinearity between steepness index and erosion rate. We find that hydrological processes can have a significant influence on how erosive a particular climatic forcing will be. By accounting for the influence of hydrology on fluvial erosion, we conclude that climate is an important control on erosion rates and long-term landscape evolution.}, language = {en} } @article{PrasicekHermanRobletal.2018, author = {Prasicek, G{\"u}nther and Herman, Frederic and Robl, J{\"o}rg and Braun, Jean}, title = {Glacial steady state topography controlled by the coupled influence of tectonics and climate}, series = {Journal of geophysical research : Earth surface}, volume = {123}, journal = {Journal of geophysical research : Earth surface}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2017JF004559}, pages = {1344 -- 1362}, year = {2018}, abstract = {Glaciers and rivers are the main agents of mountain erosion. While in the fluvial realm empirical relationships and their mathematical description, such as the stream power law, improved the understanding of fundamental controls on landscape evolution, simple constraints on glacial topography and governing scaling relations are widely lacking. We present a steady state solution for longitudinal profiles along eroding glaciers in a coupled system that includes tectonics and climate. We combined the shallow ice approximation and a glacial erosion rule to calculate ice surface and bed topography from prescribed glacier mass balance gradient and rock uplift rate. Our approach is inspired by the classic application of the stream power law for describing a fluvial steady state but with the striking difference that, in the glacial realm, glacier mass balance is added as an altitude-dependent variable. From our analyses we find that ice surface slope and glacial relief scale with uplift rate with scaling exponents indicating that glacial relief is less sensitive to uplift rate than relief in most fluvial landscapes. Basic scaling relations controlled by either basal sliding or internal deformation follow a power law with the exponent depending on the exponents for the glacial erosion rule and Glen's flow law. In a mixed scenario of sliding and deformation, complicated scaling relations with variable exponents emerge. Furthermore, a cutoff in glacier mass balance or cold ice in high elevations can lead to substantially larger scaling exponents which may provide an explanation for high relief in high latitudes.}, language = {en} } @article{PicoMitrovicaBraunetal.2018, author = {Pico, T. and Mitrovica, J. X. and Braun, Jean and Ferrier, K. L.}, title = {Glacial isostatic adjustment deflects the path of the ancestral Hudson River}, series = {Geology}, volume = {46}, journal = {Geology}, number = {7}, publisher = {American Institute of Physics}, address = {Boulder}, issn = {0091-7613}, doi = {10.1130/G40221.1}, pages = {591 -- 594}, year = {2018}, abstract = {Quantifying the pace of ice-sheet growth is critical to understanding ice-age climate and dynamics. Here, we show that the diversion of the Hudson River (northeastern North America) late in the last glaciation phase (ca. 30 ka), which some previous studies have speculated was due to glacial isostatic adjustment (GIA), can be used to infer the timing of the Laurentide Ice Sheet's growth to its maximum extent. Landscapes in the vicinity of glaciated regions have likely responded to crustal deformation produced by ice-sheet growth and decay through river drainage reorganization, given that rates of uplift and subsidence are on the order of tens of meters per thousand years. We perform global, gravitationally self-consistent simulations of GIA and input the predicted crustal deformation field into a landscape evolution model. Our calculations indicate that the eastward diversion of the Hudson River at 30 ka is consistent with exceptionally rapid growth of the Laurentide Ice Sheet late in the glaciation phase, beginning at 50-35 ka.}, language = {en} } @article{BiswasHermanKingetal.2018, author = {Biswas, R. H. and Herman, F. and King, G. E. and Braun, Jean}, title = {Thermoluminescence of feldspar as a multi-thermochronometer to constrain the temporal variation of rock exhumation in the recent past}, series = {Earth \& planetary science letters}, volume = {495}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2018.04.030}, pages = {56 -- 68}, year = {2018}, abstract = {Natural thermoluminescence (TL) in rocks reflects a dynamic equilibrium between radiation-induced TL growth and decay via thermal and athermal pathways. When rocks exhume through Earth's crust and cool from high to low temperature, this equilibrium level increases as the temperature dependent thermal decay decreases. This phenomenon can be exploited to extract thermal histories of rocks. The main advantage of TL is that a single TL glow curve has a wide range of thermal stabilities (lifetime 100 °C/Ma, whereas deeper traps, i.e. with higher activation energies, provide constraints on thermal histories for higher cooling rates (>300 °C/Ma). Finally, we show how the path of rock exhumation (i.e., depth vs. time) can be constrained using an inverse approach. The newly developed methodology is applied to rapidly cooled samples from the Namche Barwa massif, eastern Himalaya to suggest a trend in exhumation rate with time that follows an inverse correlation with global temperature and glaciers equilibrium altitude line (ELA).}, language = {en} }