An official website of the United States government
Abstract Mountain topography alters the phase, amount, and spatial distribution of precipitation. Past efforts focused on how orographic precipitation can alter spatial patterns in mean runoff, with less emphasis on how time‐varying runoff statistics may also vary with topography. Given the importance of the magnitude and frequency of runoff events to fluvial erosion, we evaluated whether orographic patterns in mean runoff and daily runoff variability can be constrained using the global WaterGAP3 water model data. Model runoff data are validated against observational data in the contiguous United States, showing agreement with mean runoff in all settings and daily runoff variability in settings where rainfall‐runoff predominates. In snowmelt‐influenced settings, runoff variability is overestimated by the water model data. Cognizant of these limitations, we use the water model data to develop relationships between mean runoff and daily runoff variability and how these are mediated by snowmelt fraction in mountain topography globally. A global analysis of topographic controls on hydroclimatic variables using a random forest model was ambiguous. Instead, relationships between topography and runoff parameters are better assessed at the mountain range scale. Rulesets linking topography to mean runoff and snowmelt fraction are developed for three mid‐latitude mountain landscapes—British Columbia, European Alps, and Greater Caucasus. Increasing topographic elevation and relief together leads to higher mean runoff and lower runoff variability due to the increasing contribution of snowmelt. The three sets of empirical relationships developed here serve as the basis for a suite of numerical experiments in our companion manuscript (Part 2, Forte & Rossi, 2024a,https://doi.org.10.1002/2023JF007327).
more »
« less
Stochastic in Space and Time: 2. Effects of Simulating Orographic Gradients in Daily Runoff Variability on Bedrock River Incision
Abstract The extent to which climate and tectonics can be coupled rests on the degree to which topography and erosion rates scales linearly. The stream power incision model (SPIM) is commonly used to interpret such relationships, but is limited in probing mechanisms. A promising modification to stream power models are stochastic‐threshold incision models (STIM) which incorporate both variability in discharge and a threshold to erosion. In this family of models, the form of the topography erosion rate relationship is largely controlled by runoff variability. Applications of STIM typically assume temporally variable, but spatially uniform and synchronous runoff generating events, an assumption that is likely broken in regions with complicated orography. To address this limitation, we develop a new 1D STIM model, which we refer to as spatial‐STIM. This modified version of STIM allows for stochasticity in both time and space and is driven by empirical relationships between topography and runoff statistics. Coupling between mean runoff and runoff variability via topography in spatial‐STIM generates highly nonlinear relationships between steady‐state topography and erosion rates. We find that whether the daily statistics of runoff are spatially linked or unlinked, which sets whether there is spatial synchronicity in the recurrence interval of runoff generating events, is a fundamental control on landscape evolution. Many empirical topography—erosion rate data sets are based on data that span across the endmember scenarios of linked versus unlinked behavior. It is thus questionable whether singular SPIM relationships fit to those data can be meaningfully related to their associated hydroclimatic conditions.
more »
« less
- PAR ID:
- 10499625
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 129
- Issue:
- 3
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present 17 new 10Be erosion rates from southern Peru sampled across an extreme orographic rainfall gradient. Using a rainfall-weighted variant of the normalized channel steepness index, ksnQ, we show that channel steepness values, and thus topography, are adjusted to spatially varying rainfall. Rocks with similar physical characteristics define distinct relationships between ksnQ and erosion rate (E), suggesting ksnQ is also resolving lithologic variations in erodibility. However, substantial uncertainty exists in parameters describing these relationships. By combining our new data with 38 published erosion rates from Peru and Bolivia, we collapse the range of compatible parameter values and resolve robust, nonlinear ksnQ–E relationships suggestive of important influences of erosional thresholds, rock properties, sediment characteristics, and temporal runoff variability. In contrast, neither climatic nor lithologic effects are clear using the traditional channel steepness metric, ksn. Our results highlight that accounting for spatial rainfall variations is essential for disentangling the multiple influences of climate, lithology, and tectonics common in mountain landscapes, which is a necessary first step toward greater understanding of how these landscapes evolve.more » « less
-
Abstract Generalizable relationships for how subdaily rainfall statistics imprint into runoff statistics are lacking. We use the Colorado Front Range, known for destructive rainfall‐triggered floods and landslides, to assess whether orographic patterns in runoff generation are a direct consequence of rainstorm climatology. Climatological analysis relies on a dense network of tipping‐bucket rain gauges and gridded precipitation frequency estimates from the National Oceanic and Atmospheric Administration to evaluate relationships among subdaily rainfall statistics, topography, and flood frequency throughout the South Platte River basin. We find that event‐scale rainfall statistics only weakly depend on elevation, suggesting that orographic gradients in runoff “extremes” are not simply a consequence of rainfall patterns. In contrast, bedrock exposure strongly varies with elevation in a way that plausibly explains enhanced runoff generation at lower elevations via reduced water storage capacity. These findings are suggestive of feedbacks between bedrock river evolution and hillslope hydrology not typically included in models of landscape evolution.more » « less
-
Abstract Prior numerical modeling work has suggested that incision into sub‐horizontal layered stratigraphy with variable erodibility induces non‐uniform erosion rates even if base‐level fall is steady and sustained. Erosion rates of cliff bands formed in the stronger rocks in a stratigraphic sequence can greatly exceed the rate of base‐level fall. Where quartz in downstream sediment is sourced primarily from the stronger, cliff‐forming units, erosion rates estimated from concentrations of cosmogenic beryllium‐10 (10Be) in detrital sediment will reflect the locally high erosion rates in retreating cliff bands. We derive theoretical relationships for threshold hillslopes and channels described by the stream‐power incision model as a quantitative guide to the potential magnitude of this amplification of10Be‐derived erosion rates above the rate of base‐level fall. Our analyses predict that the degree of erosion rate amplification is a function of bedding dip and either the ratio of rock erodibility in alternating strong and weak layers in the channel network, or the ratio of cliff to intervening‐slope gradient on threshold hillslopes. We test our predictions in the cliff‐and‐bench landscape of the Grand Staircase in southern Utah, USA. We show that detrital cosmogenic erosion rates in this landscape are significantly higher (median 300 m/Ma) than the base‐level fall rate (~75 m/Ma) determined from the incision rate of a trunk stream into a ~0.6 Ma basalt flow emplaced along a 16 km reach of the channel. We infer a 3–6‐fold range in rock strength from near‐surface P‐wave velocity measurements. The approximately four‐fold difference between the median10Be‐derived erosion rate and the long‐term rate of base‐level fall is consistent with our model and the observation that the stronger, cliff‐forming lithologies in this landscape are the primary source of quartz in detrital sediments. © 2020 John Wiley & Sons, Ltd.more » « less
-
We present a multimodel analysis for mechanistic hypothesis testing in landscape evolution theory. The study site is a watershed with well‐constrained initial and boundary conditions in which a river network locally incised 50 m over the last 13 ka. We calibrate and validate a set of 37 landscape evolution models designed to hierarchically test elements of complexity from four categories: hillslope processes, channel processes, surface hydrology, and representation of geologic materials. Comparison of each model to a base model, which uses stream power channel incision, uniform lithology, hillslope transport by linear diffusion, and surface water discharge proportional to drainage area, serves as a formal test of which elements of complexity improve model performance. Model fit is assessed using an objective function based on a direct difference between observed and simulated modern topography. A hybrid optimization scheme identifies optimal parameters and uncertainty. Multimodel analysis determines which elements of complexity improve simulation performance. Validation tests which model improvements persist when models are applied to an independent watershed. The three most important model elements are (1) spatial variation in lithology (differentiation between shale and glacial till), (2) a fluvial erosion threshold, and (3) a nonlinear relationship between slope and hillslope sediment flux. Due to nonlinear interactions between model elements, some process representations (e.g., nonlinear hillslopes) only become important when paired with the inclusion of other processes (e.g., erosion thresholds). This emphasizes the need for caution in identifying the minimally sufficient process set. Our approach provides a general framework for hypothesis testing in landscape evolution.more » « less
