Search for: All records

Gravel‐bed rivers that incise into bedrock are common worldwide. These systems have many similarities with other alluvial channels: they transport large amounts of sediment and adjust their forms in response to discharge and sediment supply. At the same time, the occurrence of bedrock incision implies behaviour that falls on a spectrum between fully detachment‐limited ‘bedrock channels’ and fully transport‐limited ‘alluvial channels’. Here, we present a mathematical model of river profile evolution that integrates bedrock erosion, gravel transport and the formation of channels whose hydraulic geometry is consistent with that of near‐threshold alluvial channels. We combine theory for five interrelated processes: bedload sediment transport in equilibrium gravel‐bed channels, channel width adjustment to flow and sediment characteristics, abrasion of bedrock by mobile sediment, plucking of bedrock and progressive loss of gravel‐sized sediment due to grain attrition. This model contributes to a growing class of models that seek to capture the dynamics of both bedrock incision and alluvial sediment transport. We demonstrate the model's ability to reproduce expected fluvial features such as inverse power law scaling between slope and area, and width and depth consistent with near‐threshold channel theory, and we discuss the role of sediment characteristics in influencing the mode of channel behaviour, erosional mechanism, channel steepness and profile concavity. 

Abstract Changes in the properties of rainfall distributions at sub‐daily scales are key to assessing soil erosion rates under climate transition. However, such changes are difficult to detect and model, especially over landscape evolution timescales. In this contribution, we validate a new catchment‐scale landscape evolution model against event‐scale runoff and sediment records. Through multi‐century numerical experiments, we also show that changes in the sub‐daily rainfall distribution, like those observed under modern climate change, can increase soil erosion rates by 40% but cannot be accurately inferred from changes in the average event properties and total rainfall. We quantify erosion and topographic trajectories associated with plausible changes in the sub‐daily rainfall distribution, highlighting scenarios in which shifting tail properties impact landscape evolution, at times, contrary to expectations based on changes in total rainfall. 

Abstract Large dendritic valley networks observed on Mars present a paleoclimate paradox. Geologic observations of Noachian units on Mars reveal a global extent of valley networks, which are believed to have been formed through incisions made by flowing water. However, most climate models predict global surface temperatures too far below the freezing point of water to support an active hydrological system. Conflicting observations and models have led to disparate theories for the climate of early Mars. In this work, we surveyed a large region of the cratered southern highlands to identify the location, elevation, and distribution of observed valley heads. These valley head locations were compared to landscape evolution simulations in which the spatial distribution of runoff was varied. The measured valley head distributions were compared to predictions from landscape evolution models for two end‐member hypotheses: (a) a warm wet climate that supported spatially distributed precipitation, and (b) surface runoff from ice cap margins, as envisioned by the Late Noachian Icy Highland model (LNIH). The observed elevation distribution in valley heads is consistent with the prediction of precipitation‐fed models, and inconsistent with models in which runoff derives exclusively from a single line‐source of high‐elevation ice‐melt. The results support the view that it is unlikely for ice caps to be the sole source of water and are consistent with the hypothesis that precipitation significantly contributed to valley network formation on ancient Mars. 

Abstract Features of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long‐term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro‐geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non‐dendritic, trellis‐like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro‐geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long‐term landscape evolution. 

Abstract Hillslope topographic change in response to climate and climate change is a key aspect of landscape evolution. The impact of short‐duration rainstorms on hillslope evolution in arid regions is persistently questioned but often not directly examined in landscape evolution studies, which are commonly based on mean climate proxies. This study focuses on hillslope surface processes responding to rainstorms in the driest regions of Earth. We present a numerical model for arid, rocky hillslopes with lithology of a softer rock layer capped by a cliff‐forming resistant layer. By representing the combined action of bedrock and clast weathering, cliff‐debris ravel, and runoff‐driven erosion, the model can reproduce commonly observed cliff‐profile morphology. Numerical experiments with a fixed base level were used to test hillslope response to cliff‐debris grain size, rainstorm intensities, and alternation between rainstorm patterns. The persistence of vertical cliffs and the pattern of sediment sorting depend on rainstorm intensities and the size of cliff debris. Numerical experiments confirm that these two variables could have driven the landscape in the Negev Desert (Israel) toward an observed spatial contrast in topographic form over the past 10⁵–10⁶ years. For a given total storm rain depth, short‐duration higher‐intensity rainstorms are more erosive, resulting in greater cliff retreat distances relative to longer, low‐intensity storms. Temporal alternation between rainstorm regimes produces hillslope profiles similar to those previously attributed to Quaternary oscillations in the mean climate. We suggest that arid hillslopes may undergo considerable geomorphic transitions solely by alternating intra‐storm patterns regardless of rainfall amounts. 

Abstract. Progress in better understanding and modeling Earth surface systems requires an ongoing integration of data and numerical models. Advances are currently hampered by technical barriers that inhibit finding, accessing, and executing modeling software with related datasets. We propose a design framework for Data Components, which are software packages that provide access to particular research datasets or types of data. Because they use a standard interface based on the Basic Model Interface (BMI), Data Components can function as plug-and-play components within modeling frameworks to facilitate seamless data–model integration. To illustrate the design and potential applications of Data Components and their advantages, we present several case studies in Earth surface processes analysis and modeling. The results demonstrate that the Data Component design provides a consistent and efficient way to access heterogeneous datasets from multiple sources and to seamlessly integrate them with various models. This design supports the creation of open data–model integration workflows that can be discovered, accessed, and reproduced through online data sharing platforms, which promotes data reuse and improves research transparency and reproducibility. 

Key Points Modeled ecosystem response to climate follows the “geo‐ecological law of distribution,” highlights the importance of ecohdyrologic refugia Woody Plant Encroachment is predicted as a three‐phase phenomenon: early establishment, rapid expansion, and woody plant equilibrium Regime shifts from grassland to shrubland are marked by vegetation cover thresholds