@article{YuanBraunGueritetal.2019, author = {Yuan, Xiaoping P. and Braun, Jean and Guerit, Laure and Rouby, D. and Cordonnier, G.}, title = {A New Efficient Method to Solve the Stream Power Law Model Taking Into Account Sediment Deposition}, series = {Journal of geophysical research : Earth surface}, volume = {124}, journal = {Journal of geophysical research : Earth surface}, number = {6}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9003}, doi = {10.1029/2018JF004867}, pages = {1346 -- 1365}, year = {2019}, abstract = {The stream power law model has been widely used to represent erosion by rivers but does not take into account the role played by sediment in modulating erosion and deposition rates. Davy and Lague (2009, ) provide an approach to address this issue, but it is computationally demanding because the local balance between erosion and deposition depends on sediment flux resulting from net upstream erosion. Here, we propose an efficient (i.e., O(N) and implicit) method to solve their equation. This means that, unlike other methods used to study the complete dynamics of fluvial systems (e.g., including the transition from detachment-limited to transport-limited behavior), our method is unconditionally stable even when large time steps are used. We demonstrate its applicability by performing a range of simulations based on a simple setup composed of an uplifting region adjacent to a stable foreland basin. As uplift and erosion progress, the mean elevations of the uplifting relief and the foreland increase, together with the average slope in the foreland. Sediments aggrade in the foreland and prograde to reach the base level where sediments are allowed to leave the system. We show how the topography of the uplifting relief and the stratigraphy of the foreland basin are controlled by the efficiency of river erosion and the efficiency of sediment transport by rivers. We observe the formation of a steady-state geometry in the uplifting region, and a dynamic steady state (i.e., autocyclic aggradation and incision) in the foreland, with aggradation and incision thicknesses up to tens of meters.}, language = {en} } @article{CordonnierBovyBraun2019, author = {Cordonnier, Guillaume and Bovy, Benoit and Braun, Jean}, title = {A versatile, linear complexity algorithm for flow routing in topographies with depressions}, series = {Earth surface dynamics}, volume = {7}, journal = {Earth surface dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-549-2019}, pages = {549 -- 562}, year = {2019}, abstract = {We present a new algorithm for solving the common problem of flow trapped in closed depressions within digital elevation models, as encountered in many applications relying on flow routing. Unlike other approaches (e.g., the Priority-Flood depression filling algorithm), this solution is based on the explicit computation of the flow paths both within and across the depressions through the construction of a graph connecting together all adjacent drainage basins. Although this represents many operations, a linear time complexity can be reached for the whole computation, making it very efficient. Compared to the most optimized solutions proposed so far, we show that this algorithm of flow path enforcement yields the best performance when used in landscape evolution models. In addition to its efficiency, our proposed method also has the advantage of letting the user choose among different strategies of flow path enforcement within the depressions (i.e., filling vs. carving). Furthermore, the computed graph of basins is a generic structure that has the potential to be reused for solving other problems as well, such as the simulation of erosion. This sequential algorithm may be helpful for those who need to, e.g., process digital elevation models of moderate size on single computers or run batches of simulations as part of an inference study.}, language = {en} } @article{MargirierBraunGautheronetal.2019, author = {Margirier, Audrey and Braun, Jean and Gautheron, Cecile and Carcaillet, Julien and Schwartz, Stephane and Jamme, Rosella Pinna and Stanley, Jessica}, title = {Climate control on Early Cenozoic denudation of the Namibian margin as deduced from new thermochronological constraints}, series = {Earth \& planetary science letters}, volume = {527}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2019.115779}, pages = {11}, year = {2019}, abstract = {The processes that control long term landscape evolution in continental interiors and, in particular, along passive margins such as in southern Africa, are still the subject of much debate (e.g. Braun, 2018). Although today the Namibian margin is characterized by an arid climate, it has experienced climatic fluctuations during the Cenozoic and, yet, to date no study has documented the potential role of climate on its erosion history. In western Namibia, the Brandberg Massif, an erosional remnant or inselberg, provides a good opportunity to document the Cenozoic denudation history of the margin using the relationship between rock cooling or exhumation ages and their elevation. Here we provide new apatite (UThSm)/He dates on the Brandberg Inselberg that range from 151 +/- 12 to 30 +/- 2 Ma. Combined with existing apatite fission track data, they yield new constraints on the denudation history of the margin. These data document two main cooling phases since continental break-up 130 Myr ago, a rapid one (similar to 10 degrees C/Myr) following break-up and a slower one (similar to 12 degrees C/Myr) between 65 and 35 Ma. We interpret them respectively to be related to escarpment erosion following rifting and continental break-up and as a phase of enhanced denudation during the Early Eocene Climatic Optimum. We propose that during the Early Eocene Climatic Optimum chemical weathering was important and contributed significantly to the denudation of the Namibian margin and the formation of a pediplain around the Brandberg and enhanced valley incision within the massif. Additionally, aridification of the region since 35 Ma has resulted in negligible denudation rates since that time. (C) 2019 Elsevier B.V. All rights reserved.}, 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{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{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{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{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{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{PicoMitrovicaPerronetal.2019, author = {Pico, Tamara and Mitrovica, Jerry X. and Perron, J. Taylor and Ferrier, Ken L. and Braun, Jean}, title = {Influence of glacial isostatic adjustment on river evolution along the US mid-Atlantic coast}, series = {Earth \& planetary science letters}, volume = {522}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2019.06.026}, pages = {176 -- 185}, year = {2019}, abstract = {Long-term river evolution depends partly on crustal deformation, which shapes the topography crossed by rivers. On glacial timescales, ice-sheet growth and decay can produce crustal vertical motion of ∼10 mm/yr resulting from the solid Earth's adjustment to variations in ice and water loads, comparable to tectonically-driven rates in the most rapidly uplifting mountains on Earth. This process of glacial isostatic adjustment (GIA) can influence river courses and drainage basins substantially, particularly near former ice margins. We explore the extent to which GIA influenced the evolution of rivers along the United States east coast during the last glacial cycle. We compute gravitationally self-consistent GIA responses that incorporate recent constraints on the Laurentide Ice Sheet history through the last glacial build-up phase, and we connect the predicted variations in topography to abrupt changes in river dynamics recorded in the Hudson, Delaware, Susquehanna, and Potomac Rivers from 40 ka to present. To the extent that increases in sediment transport capacity imply increases in river incision rate, the GIA-driven changes in slope and drainage area are consistent with episodes of erosion and sedimentation observed in the Hudson, Delaware, and Potomac Rivers, but inconsistent with the observed accelerated river incision in the Susquehanna River at 30-14 ka. These analyses add to a growing body of evidence showing that GIA strongly influences river evolution over millennial timescales.}, language = {en} }