Mountain landscapes have dynamic climates that, together with tectonic processes, influence their topographic evolution. Spatial and temporal variations in rainfall are ubiquitous in these settings as orographic precipitation patterns evolve with climate change and topography. Despite important implications such changes have for river incision, their influence is understudied. Here, we investigate how changes in rainfall pattern should affect both the steady state form and transient evolution of river profiles at the catchment scale using the stream power model. We find that spatially varied rainfall patterns can complicate steady state relationships between mean rainfall, channel steepness and fluvial relief, depending on where rainfall is concentrated in catchments, and lead to unexpected transient behavior if they are neglected. Specifically, changes in rainfall pattern cause multi‐stage transient responses that differ from responses to uniform changes in rainfall. Disparate responses by rivers that experience different rainfall conditions, particularly trunk and tributary rivers, are also an important factor in understanding catchment‐wide responses to climate change. Accounting for such disparities in sampling strategies and topographic analyses may, therefore, be vital for detecting and quantifying climate's role in landscape evolution. Lastly, we show how explicitly accounting for rainfall patterns in channel steepness indices, and thus spatial variations in erosional efficiency, may advance understanding of landscape sensitivity to climate. These results have important implications for detecting transient responses to changes in rainfall pattern (and more broadly climate), interpretation of morphometrics in steady state and transient landscapes, and quantifying the sensitivity of landscapes and erosion rates to climate.
This content will become publicly available on March 25, 2025
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,
- NSF-PAR ID:
- 10499670
- 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
-
more » « less
Abstract -
more » « less
Forecasting the timing and magnitude of snowmelt and runoff is critical to managing mountain water resources. Warming temperatures are increasing the rain–snow transition elevation and are limiting the forecasting skill of statistical models relating historical snow water equivalent to streamflow. While physically based methods are available, they require accurate estimations of the spatial and temporal distribution of meteorological variables in complex terrain. Across many mountainous areas, measurements of precipitation and other meteorological variables are limited to a few reference stations and are not adequate to resolve the complex interactions between topography and atmospheric flow. In this paper, we evaluate the ability of the Weather Research and Forecasting (WRF) Model to approximate the inputs required for a physics-based snow model, iSnobal, instead of using meteorological measurements, for the Boise River Basin (BRB) in Idaho, United States. An iSnobal simulation using station data from 40 locations in and around the BRB resulted in an average root-mean-square error (RMSE) of 4.5 mm compared with 12 SNOTEL measurements. Applying WRF forcings alone was associated with an RMSE of 10.5 mm, while including a simple bias correction to the WRF outputs of temperature and precipitation reduced the RMSE to 6.5 mm. The results highlight the utility of using WRF outputs as input to snowmelt models, as all required input variables are spatiotemporally complete. This will have important benefits in areas with sparse measurement networks and will aid snowmelt and runoff forecasting in mountainous basins.
-
Abstract. Prediction of mean annual runoff is of great interest but still poses achallenge in ungauged basins. The present work diagnoses the prediction inmean annual runoff affected by the uncertainty in estimated distribution ofsoil water storage capacity. Based on a distribution function, a waterbalance model for estimating mean annual runoff is developed, in which theeffects of climate variability and the distribution of soil water storagecapacity are explicitly represented. As such, the two parameters in themodel have explicit physical meanings, and relationships between theparameters and controlling factors on mean annual runoff are established.The estimated parameters from the existing data of watershed characteristicsare applied to 35 watersheds. The results showed that the model couldcapture 88.2 % of the actual mean annual runoff on average across thestudy watersheds, indicating that the proposed new water balance model ispromising for estimating mean annual runoff in ungauged watersheds. Theunderestimation of mean annual runoff is mainly caused by theunderestimation of the area percentage of low soil water storage capacitydue to neglecting the effect of land surface and bedrock topography. Higherspatial variability of soil water storage capacity estimated through theheight above the nearest drainage (HAND) and topographic wetness index (TWI)indicated that topography plays a crucial role in determining the actualsoil water storage capacity. The performance of mean annual runoffprediction in ungauged basins can be improved by employing better estimationof soil water storage capacity including the effects of soil, topography,and bedrock. It leads to better diagnosis of the data requirement forpredicting mean annual runoff in ungauged basins based on a newly developedprocess-based model finally.more » « less
-
more » « less
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
Abstract Simulations of the Argentine Basin have large uncertainties associated with quantities such as air‐sea exchanges of heat and carbon in current generation climate models and ocean reanalysis products. This is due to the complex topography, profound undersampling until recent years, and strong currents and mixing of subpolar and subtropical water masses in the basin. Because mixing of water masses is important here, model resolution is hypothesized to play an important role in estimating ocean quantities and determining overall budgets. We construct three regional ocean models with biogeochemistry at 1/3°, 1/6°, and 1/12° resolutions for the year 2017 to investigate heat and carbon dynamics in the region and determine the effect of model resolution on these dynamics. Initial conditions and boundary forcing from BSOSE (the Biogeochemical Southern Ocean State Estimate (Verdy & Mazloff, 2017),
https://doi.org/10.1002/2016JC012650 ) and atmospheric forcing from ERA5 are used. The models are evaluated for accuracy by comparing output to Argo and BGC‐Argo float profiles, BSOSE, and other reanalyses and mapped products. We then quantify the effect of resolution on model upper ocean heat and carbon transport and the associated air‐sea exchanges. We determine that increasing the resolution from 1/3° to 1/12° enhances the upward vertical transport and surface exchanges of heat but causes no significant effect on surface carbon fluxes despite enhancing downward transport of anomalous DIC.
