skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Numerical modeling of groundwater‐driven stream network evolution in low‐relief post‐glacial landscapes
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
Award ID(s):
1656935 1331906
PAR ID:
10304024
Author(s) / Creator(s):
 ;  ;  ;  
Date Published:
Journal Name:
Earth Surface Processes and Landforms
ISSN:
0197-9337
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Landscape dissection by rivers is a common qualitative measure of surface maturity. Quantitative studies of fluvial development over time indicate that drainage development increases non-linearly and is influenced by lithology, however, these studies typically take place over short timescales (10s of years), cover small areas, and focus on steep landscapes. In this work we use the Central Lowlands physiographic province (CL) as a natural laboratory in which we investigate rates and controls on drainage development in a post-glacial lowland landscape. Portions of the CL have been glaciated repeatedly in the Quaternary, and its topography is dominated by a patchwork of glacial landforms that have been developing drainage for 10 thousand to more than 500 thousand years. We modify the National Hydrography Dataset to estimate pre-agriculture drainage density developed over different amounts of time to reveal rates of drainage development in the CL. We find that drainage density in the CL increases non-linearly, increasing rapidly following glaciation before slowly approaching a maximum value. Much of the development is accomplished by 50 ka, well within a typical interglacial period. The apparent maximum value, ~1.5 km/km2, is comparable to the median drainage density measured in regions in the CL that have not experienced Quaternary glaciation. Our study shows that this value is likely influenced by soil sand content and regional precipitation levels. We note that while drainage density increases to an apparent maximum within an interglacial, the fluvial network is unlikely to adjust to post-glacial base level conditions within that same length of time. Our results are most consistent with a model of drainage expansion driven by the connection of closed depressions, or ‘non-contributing area’ (NCA), the portion of a watershed that does not drain to a river. We find that NCA decreases in tandem with increasing drainage density, which implies that NCA could be a measure of landscape integration that is at least as sensitive as drainage density. 
    more » « less
  2. Abstract Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well‐established for high‐relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low‐relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low‐relief wetland‐rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low‐relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment‐scale studies for low‐relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low‐relief wetland‐rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes. 
    more » « less
  3. Abstract. Calculating flow routing across a landscape is a routine process in geomorphology, hydrology, planetary science, and soil and water conservation. Flow-routing calculations often require a preprocessing step to remove depressions from a DEM to create a “flow-routing surface” that can host a continuous, integrated drainage network. However, real landscapes contain natural depressions that trap water. These are an important part of the hydrologic system and should be represented in flow-routing surfaces. Historically, depressions (or “pits”) in DEMs have been viewed as data errors, but the rapid expansion of high-resolution, high-precision DEM coverage increases the likelihood that depressions are real-world features. To address this long-standing problem of emerging significance, we developed FlowFill, an algorithm that routes a prescribed amount of runoff across the surface in order to flood depressions if enough water is available. This mass-conserving approach typically floods smaller depressions and those in wet areas, integrating drainage across them, while permitting internal drainage and disruptions to hydrologic connectivity. We present results from two sample study areas to which we apply a range of uniform initial runoff depths and report the resulting filled and unfilled depressions, the drainage network structure, and the required compute time. For the reach- to watershed-scale examples that we ran, FlowFill compute times ranged from approximately 1 to 30 min, with compute times per cell of 0.0001 to 0.006 s. 
    more » « less
  4. ABSTRACT Floodplains along low‐gradient, meandering river systems contain diverse hydrogeomorphic features, ranging from isolated depressions to hydrologically‐connected channels. These ephemerally‐flooded features inundate prior to river water overtopping all banks, enhancing river‐floodplain connectivity during moderately high flow stages. Predicting when and where ecological functions occur in floodplains requires understanding the dynamic hydrologic processes of hydrogeomorphic features, including inundation and exchange. In this study, we examined storm event‐scale inundation and exchange dynamics along a lowland, meandering river system in central Illinois (USA). We monitored surface water presence/absence, surface water level, and groundwater level across floodplain hydrogeomorphic feature types (i.e., isolated depression, backwater channel, and flow‐through channel). Using these data, we evaluated inundation onset and recession characteristics, drivers of groundwater‐surface water interactions, and direction of hydrologic exchange with the river channel. Surface water presence/absence patterns suggested inundation onset timescales were primarily controlled by microtopography and recession timescales were correlated with floodplain elevation. Employing a novel hysteresis approach for characterising groundwater‐surface water interactions, we observed distinct patterns indicating differences in water sources across hydrogeomorphic units and event characteristics. Finally, differences in hydraulic head along floodplain channels revealed that channels with multiple inlets/outlets (i.e., flow‐through channels) conveyed down‐valley flow and channels with single inlets primarily functioned as sinks of river‐derived water to the floodplain with short source periods. These results highlight the heterogeneity of hydrologic processes that occur along lowland, meandering river‐floodplains, and more specifically, point to the important role hydrogeomorphic features play in controlling dynamic connectivity within the river corridor. 
    more » « less
  5. Genetic divergence along the central Baja California Peninsula, Mexico, has been hypothesized to reflect a Pliocene cross-peninsular seaway that previously isolated northern and southern populations of terrestrial plants and animals. One way to test this hypothesis is through quantitative analysis of relict channels preserved on low-relief paleo-surfaces. Recognition of tidal channels on relict landscapes offers a powerful tool for reconstructing past sea level in tectonically active arid coastal regions where crustal uplift results in relative sea-level fall and preservation of ancient channel networks. This method requires reliable criteria to distinguish fluvial versus tidal channels, which is challenging due to the overlap of standard metrics for the two channel types, and possible inheritance or overprinting of geometries. We improve the utility of existing metrics and explore the potential for identifying paleo-sea-level indicators by analyzing modern and ancient channels to identify unique patterns in planform geometry and to evaluate their applicability for classifying tidal versus fluvial origins. Preliminary measurements of geographically diverse modern systems reveal distinct, quantifiable differences between the two channel types in along-channel curvature, width, and wavelet spectra. Modern tidal channels display a pronounced and systematic down-channel increase in channel width and decrease in curvature. In contrast, modern fluvial channels do not display spatial patterns in channel width and curvature along their lengths. These patterns provide diagnostic criteria that can be paired with wavelet analysis of meander belts to classify the paleoenvironment of ancient channels based on their planform geometry. We apply this approach to evaluate the origin of channels preserved on relict landscapes in the San Ignacio trough in the central Baja California peninsula, a former low-relief embayment of the Pacific Ocean. Early results reveal the presence of ancient tidal channel networks at elevations of ~ 50-300 m above modern sea level on surfaces that are independently dated to be ca. 4-5 Ma. These findings provide evidence for post 4-Ma uplift in the mid-peninsular region and an ancient tidal environment that may have isolated northern and southern terrestrial populations. 
    more » « less