@article{WickertSchildgen2019, author = {Wickert, Andrew D. and Schildgen, Taylor F.}, title = {Long-profile evolution of transport-limited gravel-bed rivers}, series = {Earth surface dynamics}, volume = {7}, journal = {Earth surface dynamics}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-17-2019}, pages = {17 -- 43}, year = {2019}, abstract = {Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.}, language = {en} } @misc{WickertSchildgen2019, author = {Wickert, Andrew D. and Schildgen, Taylor F.}, title = {Long-profile evolution of transport-limited gravel-bed rivers}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {680}, issn = {1866-8372}, doi = {10.25932/publishup-42571}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-425718}, pages = {27}, year = {2019}, abstract = {Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift rate, and steepness increases with an increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: to match observed channel concavities, gravel-sized sediments must weather and fine - typically rapidly - and valleys typically should widen gradually. To match the empirical square-root width-discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.}, language = {en} } @misc{Wickert2016, author = {Wickert, Andrew D.}, title = {Open-source modular solutions for flexural isostasy}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {498}, issn = {1866-8372}, doi = {10.25932/publishup-40836}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-408366}, pages = {21}, year = {2016}, abstract = {Isostasy is one of the oldest and most widely applied concepts in the geosciences, but the geoscientific community lacks a coherent, easy-to-use tool to simulate flexure of a realistic (i.e., laterally heterogeneous) lithosphere under an arbitrary set of surface loads. Such a model is needed for studies of mountain building, sedimentary basin formation, glaciation, sea-level change, and other tectonic, geodynamic, and surface processes. Here I present gFlex (for GNU flexure), an open-source model that can produce analytical and finite difference solutions for lithospheric flexure in one (profile) and two (map view) dimensions. To simulate the flexural isostatic response to an imposed load, it can be used by itself or within GRASS GIS for better integration with field data. gFlex is also a component with the Community Surface Dynamics Modeling System (CSDMS) and Landlab modeling frameworks for coupling with a wide range of Earth-surface-related models, and can be coupled to additional models within Python scripts. As an example of this in-script coupling, I simulate the effects of spatially variable lithospheric thickness on a modeled Iceland ice cap. Finite difference solutions in gFlex can use any of five types of boundary conditions: 0-displacement, 0-slope (i.e., clamped); 0-slope, 0-shear; 0-moment, 0-shear (i.e., broken plate); mirror symmetry; and periodic. Typical calculations with gFlex require << 1 s to similar to 1 min on a personal laptop computer. These characteristics - multiple ways to run the model, multiple solution methods, multiple boundary conditions, and short compute time - make gFlex an effective tool for flexural isostatic modeling across the geosciences.}, language = {en} } @article{Wickert2016, author = {Wickert, Andrew D.}, title = {Open-source modular solutions for flexural isostasy: gFlex v1.0}, series = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, volume = {9}, journal = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1991-959X}, doi = {10.5194/gmd-9-997-2016}, pages = {997 -- 1017}, year = {2016}, abstract = {Isostasy is one of the oldest and most widely applied concepts in the geosciences, but the geoscientific community lacks a coherent, easy-to-use tool to simulate flexure of a realistic (i.e., laterally heterogeneous) lithosphere under an arbitrary set of surface loads. Such a model is needed for studies of mountain building, sedimentary basin formation, glaciation, sea-level change, and other tectonic, geodynamic, and surface processes. Here I present gFlex (for GNU flexure), an open-source model that can produce analytical and finite difference solutions for lithospheric flexure in one (profile) and two (map view) dimensions. To simulate the flexural isostatic response to an imposed load, it can be used by itself or within GRASS GIS for better integration with field data. gFlex is also a component with the Community Surface Dynamics Modeling System (CSDMS) and Landlab modeling frameworks for coupling with a wide range of Earth-surface-related models, and can be coupled to additional models within Python scripts. As an example of this in-script coupling, I simulate the effects of spatially variable lithospheric thickness on a modeled Iceland ice cap. Finite difference solutions in gFlex can use any of five types of boundary conditions: 0-displacement, 0-slope (i.e., clamped); 0-slope, 0-shear; 0-moment, 0-shear (i.e., broken plate); mirror symmetry; and periodic. Typical calculations with gFlex require << 1 s to similar to 1 min on a personal laptop computer. These characteristics - multiple ways to run the model, multiple solution methods, multiple boundary conditions, and short compute time - make gFlex an effective tool for flexural isostatic modeling across the geosciences.}, language = {en} } @article{TofeldeSchildgenSavietal.2017, author = {Tofelde, Stefanie and Schildgen, Taylor F. and Savi, Sara and Pingel, Heiko and Wickert, Andrew D. and Bookhagen, Bodo and Wittmann, Hella and Alonso, Ricardo N. and Cottle, John and Strecker, Manfred}, title = {100 kyr fluvial cut-and-fill terrace cycles since the Middle Pleistocene in the southern Central Andes, NW Argentina}, series = {Earth \& planetary science letters}, volume = {473}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2017.06.001}, pages = {141 -- 153}, year = {2017}, abstract = {Fluvial fill terraces in intermontane basins are valuable geomorphic archives that can record tectonically and/or climatically driven changes of the Earth-surface process system. However, often the preservation of fill terrace sequences is incomplete and/or they may form far away from their source areas, complicating the identification of causal links between forcing mechanisms and landscape response, especially over multi-millennial timescales. The intermontane Toro Basin in the southern Central Andes exhibits at least five generations of fluvial terraces that have been sculpted into several-hundred-meter-thick Quaternary valley-fill conglomerates. New surface-exposure dating using nine cosmogenic Be-10 depth profiles reveals the successive abandonment of these terraces with a 100 kyr cyclicity between 75 +/- 7 and 487 +/- 34 ka. Depositional ages of the conglomerates, determined by four Al-26/Be-10 burial samples and U-Pb zircon ages of three intercalated volcanic ash beds, range from 18 +/- 141 to 936 +/- 170 ka, indicating that there were multiple cut-and-fill episodes. Although the initial onset of aggradation at similar to 1 Ma and the overall net incision since ca. 500 ka can be linked to tectonic processes at the narrow basin outlet, the superimposed 100 kyr cycles of aggradation and incision are best explained by eccentricity-driven climate change. Within these cycles, the onset of river incision can be correlated with global cold periods and enhanced humid phases recorded in paleoclimate archives on the adjacent Bolivian Altiplano, whereas deposition occurred mainly during more arid phases on the Altiplano and global interglacial periods. We suggest that enhanced runoff during global cold phases - due to increased regional precipitation rates, reduced evapotranspiration, or both - resulted in an increased sediment-transport capacity in the Toro Basin, which outweighed any possible increases in upstream sediment supply and thus triggered incision. Compared with two nearby basins that record precessional (21-kyr) and long-eccentricity (400-kyr) forcing within sedimentary and geomorphic archives, the recorded cyclicity scales with the square of the drainage basin length. (C) 2017 Elsevier B.V. All rights reserved.}, language = {en} } @misc{TofeldeSaviWickertetal.2019, author = {Tofelde, Stefanie and Savi, Sara and Wickert, Andrew D. and Bufe, Aaron and Schildgen, Taylor F.}, title = {Alluvial channel response to environmental perturbations}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {762}, issn = {1866-8372}, doi = {10.25932/publishup-43718}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-437185}, pages = {609 -- 631}, year = {2019}, abstract = {The sensitivity of fluvial systems to tectonic and climatic boundary conditions allows us to use the geomorphic and stratigraphic records as quantitative archives of past climatic and tectonic conditions. Thus, fluvial terraces that form on alluvial fans and floodplains as well as the rate of sediment export to oceanic and continental basins are commonly used to reconstruct paleoenvironments. However, we currently lack a systematic and quantitative understanding of the transient evolution of fluvial systems and their associated sediment storage and release in response to changes in base level, water input, and sediment input. Such knowledge is necessary to quantify past environmental change from terrace records or sedimentary deposits and to disentangle the multiple possible causes for terrace formation and sediment deposition. Here, we use a set of seven physical experiments to explore terrace formation and sediment export from a single, braided channel that is perturbed by changes in upstream water discharge or sediment supply, or through downstream base-level fall. Each perturbation differently affects (1) the geometry of terraces and channels, (2) the timing of terrace cutting, and (3) the transient response of sediment export from the basin. In general, an increase in water discharge leads to near-instantaneous channel incision across the entire fluvial system and consequent local terrace cutting, thus preserving the initial channel slope on terrace surfaces, and it also produces a transient increase in sediment export from the system. In contrast, a decreased upstream sediment-supply rate may result in longer lag times before terrace cutting, leading to terrace slopes that differ from the initial channel slope, and also lagged responses in sediment export. Finally, downstream base-level fall triggers the upstream propagation of a diffuse knickzone, forming terraces with upstream-decreasing ages. The slope of terraces triggered by base-level fall mimics that of the newly adjusted active channel, whereas slopes of terraces triggered by a decrease in upstream sediment discharge or an increase in upstream water discharge are steeper compared to the new equilibrium channel. By combining fillterrace records with constraints on sediment export, we can distinguish among environmental perturbations that would otherwise remain unresolved when using just one of these records.}, language = {en} } @article{TofeldeSaviWickertetal.2019, author = {Tofelde, Stefanie and Savi, Sara and Wickert, Andrew D. and Bufe, Aaron and Schildgen, Taylor F.}, title = {Alluvial channel response to environmental perturbations}, series = {Earth Surface Dynamics}, volume = {7}, journal = {Earth Surface Dynamics}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-7-609-2019}, pages = {609 -- 631}, year = {2019}, abstract = {The sensitivity of fluvial systems to tectonic and climatic boundary conditions allows us to use the geomorphic and stratigraphic records as quantitative archives of past climatic and tectonic conditions. Thus, fluvial terraces that form on alluvial fans and floodplains as well as the rate of sediment export to oceanic and continental basins are commonly used to reconstruct paleoenvironments. However, we currently lack a systematic and quantitative understanding of the transient evolution of fluvial systems and their associated sediment storage and release in response to changes in base level, water input, and sediment input. Such knowledge is necessary to quantify past environmental change from terrace records or sedimentary deposits and to disentangle the multiple possible causes for terrace formation and sediment deposition. Here, we use a set of seven physical experiments to explore terrace formation and sediment export from a single, braided channel that is perturbed by changes in upstream water discharge or sediment supply, or through downstream base-level fall. Each perturbation differently affects (1) the geometry of terraces and channels, (2) the timing of terrace cutting, and (3) the transient response of sediment export from the basin. In general, an increase in water discharge leads to near-instantaneous channel incision across the entire fluvial system and consequent local terrace cutting, thus preserving the initial channel slope on terrace surfaces, and it also produces a transient increase in sediment export from the system. In contrast, a decreased upstream sediment-supply rate may result in longer lag times before terrace cutting, leading to terrace slopes that differ from the initial channel slope, and also lagged responses in sediment export. Finally, downstream base-level fall triggers the upstream propagation of a diffuse knickzone, forming terraces with upstream-decreasing ages. The slope of terraces triggered by base-level fall mimics that of the newly adjusted active channel, whereas slopes of terraces triggered by a decrease in upstream sediment discharge or an increase in upstream water discharge are steeper compared to the new equilibrium channel. By combining fillterrace records with constraints on sediment export, we can distinguish among environmental perturbations that would otherwise remain unresolved when using just one of these records.}, language = {en} } @article{TofeldeDuesingSchildgenetal.2018, author = {Tofelde, Stefanie and Duesing, Walter and Schildgen, Taylor F. and Wickert, Andrew D. and Wittmann, Hella and Alonso, Ricardo N. and Strecker, Manfred}, title = {Effects of deep-seated versus shallow hillslope processes on cosmogenic Be-10 concentrations in fluvial sand and gravel}, series = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, volume = {43}, journal = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, number = {15}, publisher = {Wiley}, address = {Hoboken}, issn = {0197-9337}, doi = {10.1002/esp.4471}, pages = {3086 -- 3098}, year = {2018}, abstract = {Terrestrial cosmogenic nuclide (TCN) concentrations in fluvial sediment, from which denudation rates are commonly inferred, can be affected by hillslope processes. TCN concentrations in gravel and sand may differ if localized, deep-excavation processes (e.g. landslides, debris flows) affect the contributing catchment, whereas the TCN concentrations of sand and gravel tend to be more similar when diffusional processes like soil creep and sheetwash are dominant. To date, however, no study has systematically compared TCN concentrations in different detrital grain-size fractions with a detailed inventory of hillslope processes from the entire catchment. Here we compare concentrations of the TCN Be-10 in 20 detrital sand samples from the Quebrada del Toro (southern Central Andes, Argentina) to a hillslope-process inventory from each contributing catchment. Our comparison reveals a shift from low-slope gullying and scree production in slowly denuding, low-slope areas to steep-slope gullying and landsliding in fast-denuding, steep areas. To investigate whether the nature of hillslope processes (locally excavating or more uniformly denuding) may be reflected in a comparison of the Be-10 concentrations of sand and gravel, we define the normalized sand-gravel index (NSGI) as the Be-10-concentration difference between sand and gravel divided by their summed concentrations. We find a positive, linear relationship between the NSGI and median slope, such that our NSGI values broadly reflect the shift in hillslope processes from low-slope gullying and scree production to steep-slope gullying and landsliding. Higher NSGI values characterize regions affected by steep-slope gullying or landsliding. We relate the large scatter in the relationship, which is exhibited particularly in low-slope areas, to reduced hillslope-channel connectivity and associated transient sediment storage within those catchments. While high NSGI values in well-connected catchments are a reliable signal of deep-excavation processes, hillslope excavation processes may not be reliably recorded by NSGI values where sediment experiences transient storage. (c) 2018 John Wiley \& Sons, Ltd.}, language = {en} } @article{MeyScherlerWickertetal.2016, author = {Mey, J{\"u}rgen and Scherler, Dirk and Wickert, Andrew D. and Egholm, David L. and Tesauro, Magdala and Schildgen, Taylor F. and Strecker, Manfred}, title = {Glacial isostatic uplift of the European Alps}, series = {Nature Communications}, volume = {7}, journal = {Nature Communications}, publisher = {Nature Publ. Group}, address = {London}, issn = {2041-1723}, doi = {10.1038/ncomms13382}, pages = {2357 -- 2371}, year = {2016}, language = {en} } @article{BufeTurowskiBurbanketal.2019, author = {Bufe, Aaron and Turowski, Jens M. and Burbank, Douglas W. and Paola, Chris and Wickert, Andrew D. and Tofelde, Stefanie}, title = {Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment}, series = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, volume = {44}, journal = {Earth surface processes and landforms : the journal of the British Geomorphological Research Group}, number = {14}, publisher = {Wiley}, address = {Hoboken}, issn = {0197-9337}, doi = {10.1002/esp.4710}, pages = {2823 -- 2836}, year = {2019}, abstract = {Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights.}, language = {en} } @article{BallatoCifelliHeidarzadehetal.2017, author = {Ballato, Paolo and Cifelli, Francesca and Heidarzadeh, Ghasem and Ghassemi, Mohammad R. and Wickert, Andrew D. and Hassanzadeh, Jamshid and Dupont-Nivet, Guillaume and Balling, Philipp and Sudo, Masafumi and Zeilinger, Gerold and Schmitt, Axel K. and Mattei, Massimo and Strecker, Manfred}, title = {Tectono-sedimentary evolution of the northern Iranian Plateau: insights from middle-late Miocene foreland-basin deposits}, series = {Basin research}, volume = {29}, journal = {Basin research}, publisher = {Wiley}, address = {Hoboken}, issn = {0950-091X}, doi = {10.1111/bre.12180}, pages = {417 -- 446}, year = {2017}, abstract = {Sedimentary basins in the interior of orogenic plateaus can provide unique insights into the early history of plateau evolution and related geodynamic processes. The northern sectors of the Iranian Plateau of the Arabia-Eurasia collision zone offer the unique possibility to study middle-late Miocene terrestrial clastic and volcaniclastic sediments that allow assessing the nascent stages of collisional plateau formation. In particular, these sedimentary archives allow investigating several debated and poorly understood issues associated with the long-term evolution of the Iranian Plateau, including the regional spatio-temporal characteristics of sedimentation and deformation and the mechanisms of plateau growth. We document that middle-late Miocene crustal shortening and thickening processes led to the growth of a basement-cored range (Takab Range Complex) in the interior of the plateau. This triggered the development of a foreland-basin (Great Pari Basin) to the east between 16.5 and 10.7Ma. By 10.7Ma, a fast progradation of conglomerates over the foreland strata occurred, most likely during a decrease in flexural subsidence triggered by rock uplift along an intraforeland basement-cored range (Mahneshan Range Complex). This was in turn followed by the final incorporation of the foreland deposits into the orogenic system and ensuing compartmentalization of the formerly contiguous foreland into several intermontane basins. Overall, our data suggest that shortening and thickening processes led to the outward and vertical growth of the northern sectors of the Iranian Plateau starting from the middle Miocene. This implies that mantle-flow processes may have had a limited contribution toward building the Iranian Plateau in NW Iran.}, language = {en} }