TY - JOUR A1 - Wickert, Andrew D. T1 - Open-source modular solutions for flexural isostasy: gFlex v1.0 JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2016 U6 - https://doi.org/10.5194/gmd-9-997-2016 SN - 1991-959X SN - 1991-9603 VL - 9 SP - 997 EP - 1017 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Mey, Jürgen A1 - Scherler, Dirk A1 - Wickert, Andrew D. A1 - Egholm, David L. A1 - Tesauro, Magdala A1 - Schildgen, Taylor F. A1 - Strecker, Manfred T1 - Glacial isostatic uplift of the European Alps JF - Nature Communications Y1 - 2016 U6 - https://doi.org/10.1038/ncomms13382 SN - 2041-1723 VL - 7 SP - 2357 EP - 2371 PB - Nature Publ. Group CY - London ER - TY - JOUR A1 - Ballato, Paolo A1 - Cifelli, Francesca A1 - Heidarzadeh, Ghasem A1 - Ghassemi, Mohammad R. A1 - Wickert, Andrew D. A1 - Hassanzadeh, Jamshid A1 - Dupont-Nivet, Guillaume A1 - Balling, Philipp A1 - Sudo, Masafumi A1 - Zeilinger, Gerold A1 - Schmitt, Axel K. A1 - Mattei, Massimo A1 - Strecker, Manfred T1 - Tectono-sedimentary evolution of the northern Iranian Plateau: insights from middle-late Miocene foreland-basin deposits JF - Basin research N2 - 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. Y1 - 2017 U6 - https://doi.org/10.1111/bre.12180 SN - 0950-091X SN - 1365-2117 VL - 29 SP - 417 EP - 446 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Tofelde, Stefanie A1 - Schildgen, Taylor F. A1 - Savi, Sara A1 - Pingel, Heiko A1 - Wickert, Andrew D. A1 - Bookhagen, Bodo A1 - Wittmann, Hella A1 - Alonso, Ricardo N. A1 - Cottle, John A1 - Strecker, Manfred T1 - 100 kyr fluvial cut-and-fill terrace cycles since the Middle Pleistocene in the southern Central Andes, NW Argentina JF - Earth & planetary science letters N2 - 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. KW - Be-10 depth-profiles KW - surface inflation KW - aggradation-incision cycles KW - glacial-interglacial cycles KW - landscape response to climate change KW - Eastern Cordillera Y1 - 2017 U6 - https://doi.org/10.1016/j.epsl.2017.06.001 SN - 0012-821X SN - 1385-013X VL - 473 SP - 141 EP - 153 PB - Elsevier CY - Amsterdam ER - TY - JOUR A1 - Tofelde, Stefanie A1 - Duesing, Walter A1 - Schildgen, Taylor F. A1 - Wickert, Andrew D. A1 - Wittmann, Hella A1 - Alonso, Ricardo N. A1 - Strecker, Manfred T1 - Effects of deep-seated versus shallow hillslope processes on cosmogenic Be-10 concentrations in fluvial sand and gravel JF - Earth surface processes and landforms : the journal of the British Geomorphological Research Group N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1002/esp.4471 SN - 0197-9337 SN - 1096-9837 VL - 43 IS - 15 SP - 3086 EP - 3098 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Wickert, Andrew D. A1 - Schildgen, Taylor F. T1 - Long-profile evolution of transport-limited gravel-bed rivers JF - Earth surface dynamics N2 - 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. Y1 - 2019 U6 - https://doi.org/10.5194/esurf-7-17-2019 SN - 2196-6311 SN - 2196-632X VL - 7 IS - 1 SP - 17 EP - 43 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Bufe, Aaron A1 - Turowski, Jens M. A1 - Burbank, Douglas W. A1 - Paola, Chris A1 - Wickert, Andrew D. A1 - Tofelde, Stefanie T1 - Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment JF - Earth surface processes and landforms : the journal of the British Geomorphological Research Group N2 - 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. KW - braided alluvial rivers KW - physical experiments KW - channel migration Y1 - 2019 U6 - https://doi.org/10.1002/esp.4710 SN - 0197-9337 SN - 1096-9837 VL - 44 IS - 14 SP - 2823 EP - 2836 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Tofelde, Stefanie A1 - Savi, Sara A1 - Wickert, Andrew D. A1 - Bufe, Aaron A1 - Schildgen, Taylor F. T1 - Alluvial channel response to environmental perturbations BT - fill-terrace formation and sediment-signal disruption JF - Earth Surface Dynamics N2 - 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. Y1 - 2019 U6 - https://doi.org/10.5194/esurf-7-609-2019 SN - 2196-6311 SN - 2196-632X VL - 7 SP - 609 EP - 631 PB - Copernicus Publ. CY - Göttingen ER -