An official website of the United States government
Abstract The hydrologic dynamics and geomorphic evolution of watersheds are intimately coupled—runoff generation and water storage are controlled by topography and properties of the surface and subsurface, while also affecting the evolution of those properties over geologic time. However, the large disparity between their timescales has made it difficult to examine interdependent controls on emergent hydrogeomorphic properties, such as hillslope length, drainage density, and extent of surface saturation. In this study, we develop a new model coupling hydrology and landscape evolution to explore how runoff generation affects long‐term catchment evolution, and analyze numerical results using a nondimensional scaling framework. We focus on hydrologic processes dominating in humid climates where storm runoff primarily arises from shallow subsurface flow and from precipitation on saturated areas. The model solves hydraulic groundwater equations to predict the water‐table elevation given prescribed, constant groundwater recharge. Water in excess of the subsurface capacity for transport becomes overland flow, which generates shear stress on the surface and may detach and transport sediment. This affects the landscape form that in turn affects runoff generation. We show that (a) four dimensionless parameters describe the possible steady state landscapes that coevolve under steady recharge; (b) hillslope length increases with increasing transmissivity relative to the recharge rate; (c) three topographic metrics—steepness index, Laplacian curvature, and topographic index—together provide a basis for interpreting landscapes that have coevolved with runoff generated via shallow subsurface flow. Finally we discuss the possibilities and limitations for quantitative comparisons between the model results and real landscapes.
more »
« less
Catchment Coevolution and the Geomorphic Origins of Variable Source Area Hydrology
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.
more »
« less
- PAR ID:
- 10577147
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 60
- Issue:
- 6
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
For landscapes to achieve a topographic steady state, they require steady tectonic uplift and climate, and a bedrock that is uniformly erodible in the vertical direction. Basic landscape evolution models predict that incising drainage networks will eventually reach a static geometric equilibrium – that is, the map-view channel pattern will remain constant. In contrast, natural rivers typically incise through heterogeneous bedrock, which can force reorganization of the drainage structure. To investigate how lithological variability can force landscape reorganization, we draw inspiration from formerly glaciated portions of the upper Mississippi Valley. In this region, depth-to-bedrock maps reveal buried dendritic river networks dissecting paleozoic sedimentary rock. During the Pleistocene, ice advance buried the bedrock topography with glacial till, resurfacing the landscape and resetting the landscape evolution clock. As newly formed drainage networks develop and incise into the till-covered surface, they exhume the buried bedrock topography. This then leads to a geomorphic "decision point": Will the rivers follow the course of the bedrock paleodrainage network? Or will they maintain their new pattern? Using a numerical landscape evolution model, we find that two parameters determine this decision: (1) the contrast between the rock erodibility of the glacial till (more erodible) and of the buried sedimentary rock (less erodible) and (2) the orientation of the surface drainage network with respect to the buried network. We find that as the erodibility contrast increases, the drainage pattern is more likely to reorganize to follow the buried bedrock valleys. Additionally, as the alignment of the two networks increases, the surface drainage network also tends to restructure itself to follow the paleodrainage network. However, when there is less contrast and/or alignment, the surface drainage pattern becomes superimposed on the bedrock topography, with streams cutting across buried bedrock ridges. Our results agree with field studies demonstrating that variability in erodibility exerts a first-order control on landscape evolution and morphology. Our findings can provide insight into how lithologic variation affects surface processes, drives drainage reorganization, and creates geopatterns.more » « less
-
Temporal and spatial variations of tectonic rock uplift are generally thought to be the main controls on long-term erosion rates in various landscapes. However, rivers continuously lengthen and capture drainages in strike-slip fault systems due to ongoing motion across the fault, which can induce changes in landscape forms, drainage networks, and local erosion rates. Located along the restraining bend of the San Andreas Fault, the San Bernardino Mountains provide a suitable location for assessing the influence of topographic disequilibrium from perturbations by tectonic forcing and channel reorganization on measured erosion rates. In this study, we measured 17 new basin-averaged erosion rates using cosmogenic 10Be in river sands (hereafter, 10Be-derived erosion rates) and compiled 31 10Be-derived erosion rates from previous work. We quantify the degree of topographic disequilibrium using topographic analysis by examining hillslope and channel decoupling, the areal extent of pre-uplift surface, and drainage divide asymmetry across various landscapes. Similar to previous work, we find that erosion rates generally increase from north to south across the San Bernardino Mountains, reflecting a southward increase in tectonic activity. However, a comparison between 10Be-derived erosion rates and various topographic metrics in the southern San Bernardino Mountains suggests that the presence of transient landscape features such as relict topography and drainage-divide migration may explain local variations in 10Be-derived erosion rates. Our work shows that coupled analysis of erosion rates and topographic metrics provides tools for assessing the influence of tectonic uplift and channel reorganization on landscape evolution and 10Be-derived erosion rates in an evolving strike-slip restraining bend.more » « less
-
Abstract The connection between topography and erosion rate is central to understanding landscape evolution and sediment hazards. However, investigation of this relationship in steep landscapes has been limited due to expectations of: (a) decoupling between erosion rate and “threshold” hillslope morphology; and (b) bias in detrital cosmogenic nuclide erosion rates due to deep‐seated landslides. Here we compile 120 new and published10Be erosion rates from catchments in the San Gabriel Mountains, California, and show that hillslope morphology and erosion rate are coupled for slopes approaching 50° due to progressive exposure of bare bedrock with increasing erosion rate. We find no evidence for drainage area dependence in10Be erosion rates in catchments as small as 0.09 km2, and we show that landslide deposits influence erosion rate estimates mainly by adding scatter. Our results highlight the potential and importance of sampling small catchments to better understand steep hillslope processes.more » « less
-
In the low-relief post-glacial landscapes of the Central Lowlands of the United States, fluvial networks formed and expanded following deglaciation despite the low slopes and large fraction of the land surface occupied by closed depressions. Low relief topography allows for subtle surface water divides and increases the likelihood that groundwater divides do not coincide with surface water divides. We investigate how groundwater transfer across subtle surface water divides facilitates channel network expansion using a numerical model built on the Landlab platform. Our model simulates surface and subsurface water routing and fluvial erosion. We consider two end-member scenarios for surface water routing, one in which surface water in closed depressions is forced to connect to basin outlets (routing) and one in which surface water in closed depressions is lost to evapotranspiration (no routing). Groundwater is modeled as fully saturated flow within a confined aquifer. Groundwater emerges as surface water where the landscape has eroded to a specified depth. We held the total water flux constant and varied the fraction of water introduced as groundwater versus precipitation. Channel growth is significantly faster in routing cases than no-routing cases given identical groundwater fractions. In both routing and no-routing cases, channel expansion is fastest when ~30% of the total water enters the system as groundwater. Groundwater contributions also produce distinctive morphology including steepened channel profiles below groundwater seeps. Groundwater head gradients evolve with topography and groundwater-fed channels can grow more quickly than channels with larger surface water catchments. We conclude that rates of channel network growth in low-relief post-glacial areas are sensitive to groundwater contributions. More broadly, our findings suggest that landscape evolution models may benefit from more detailed representation of hydrologic processes.more » « less
