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: Drainage network signatures of landscape reorganization in post‐glacial landscapes
Earth's drainage networks encode clues that can be used to decipher geologic and geomorphic history. Dendritic drainage patterns, the most common, typically form on approximately homogeneous bedrock. Variations in rock properties or lithologic structure can disrupt dendritic patterns and form, e.g., trellis or rectangular networks. Although textbooks include such lithological–drainage links, the mechanisms driving drainage reorganization via lithologic variability remain poorly understood. To cast light on this mystery, we study drainage patterns in post-glacial landscapes of the Upper Mississippi River Valley (UMRV). Pleistocene glaciers deposited till across parts of this region, burying a landscape of fluvially dissected sedimentary rock whose buried valley network differs from modern-day drainage patterns. As the current river network erodes and exhumes the bedrock, it comes to a geomorphic "decision point": Does it reorganize to recreate the paleodrainage network, or does it maintain its new drainage pattern? To understand this decision-making, we conducted idealized landscape evolution modeling experiments. Modeled landscapes that reintegrated more of the paleodrainage network exhibited higher tortuosity, measured by dividing the real flowpath length by shortest path-length to the outlet, and obtuse tributary-junction angles. We apply this metric to two adjacent landscape types in the UMRV: (1) never glaciated (Driftless Area, DA) and (2) formerly till-mantled (Driftless-style Area, DSA), and measure the basin-averaged tortuosity for sub-basins draining streams of order 1 through 7. Across the UMRV, tortuosity increases as the maximum stream order of the sub-basin increases. For each order, tortuosity is statistically higher in areas that had been previously buried and re-exhumed (DSA) than the DA, indicating that the rivers in the DSA have reintegrated the paleodrainage network since deglaciation. For the 1st and 2nd order sub-basins, the mean basin-averaged tortuosity in the DSA is ~1-2% higher than the DA (p-value < 0.01) and ~10-14% higher (p-value < 0.01) in the 6th and 7th order sub-basins. Our analysis suggests that a drainage-based metric, tortuosity, can identify landscapes where lithological heterogeneity or structure plays a dominant role in drainage organization.  more » « less
Award ID(s):
2052938
PAR ID:
10509060
Author(s) / Creator(s):
; ;
Publisher / Repository:
American Geophysical Union
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Dendritic, i.e., tree-like, river networks are ubiquitous features on Earth’s landscapes; however, how and why river networks organize themselves into this form are incompletely understood. A branching pattern has been argued to be an optimal state. Therefore, we should expect models of river evolution to drastically reorganize (suboptimal) purely nondendritic networks into (more optimal) dendritic networks. To date, current physically based models of river basin evolution are incapable of achieving this result without substantial allogenic forcing. Here, we present a model that does indeed accomplish massive drainage reorganization. The key feature in our model is basin-wide lateral incision of bedrock channels. The addition of this submodel allows for channels to laterally migrate, which generates river capture events and drainage migration. An important factor in the model that dictates the rate and frequency of drainage network reorganization is the ratio of two parameters, the lateral and vertical rock erodibility constants. In addition, our model is unique from others because its simulations approach a dynamic steady state. At a dynamic steady state, drainage networks persistently reorganize instead of approaching a stable configuration. Our model results suggest that lateral bedrock incision processes can drive major drainage reorganization and explain apparent long-lived transience in landscapes on Earth. 
    more » « less
  3. 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
  4. Abstract Terrestrial cosmogenic nuclides (TCN) are widely employed to infer denudation rates in mountainous landscapes. The calculation of an inferred denudation rate (Dinf) from TCN concentrations is typically performed under the assumptions that denudation rates were steady during TCN accumulation and that soil chemical weathering negligibly impacted soil mineral abundances. In many landscapes, however, denudation rates were not steady and soil composition was significantly impacted by chemical weathering, which complicates interpretation of TCN concentrations. We present a landscape evolution model that computes transient changes in topography, soil thickness, soil mineralogy, and soil TCN concentrations. We used this model to investigate TCN responses in transient landscapes by imposing idealized perturbations in tectonically (rock uplift rate) and climatically sensitive parameters (soil production efficiency, hillslope transport efficiency, and mineral dissolution rate) on initially steady‐state landscapes. These experiments revealed key insights about TCN responses in transient landscapes. (a) Accounting for soil chemical erosion is necessary to accurately calculateDinf. (b) Responses ofDinfto tectonic perturbations differ from those to climatic perturbations, suggesting that spatial and temporal patterns inDinfare signatures of perturbation type and magnitude. (c) If soil chemical erosion is accounted for, basin‐averagedDinfinferred from TCN in stream sediment closely tracks actual basin‐averaged denudation rate, showing thatDinfis a reasonable proxy for actual denudation rate, even in many transient landscapes. (d) Response times ofDinfto perturbations increase with hillslope length, implying that response times should be sensitive to the climatic, biological, and lithologic processes that control hillslope length. 
    more » « less
  5. Taken together, lakes and drained lake basins may cover up to 80% of the lowland landscapes in permafrost regions of the Arctic. Lake formation, growth, and drainage in lowland permafrost regions create a terrestrial and aquatic landscape mosaic of importance to geomorphic and hydrologic processes, tundra vegetation communities, permafrost and ground-ice characteristics, biogeochemical cycling, wildlife habitat, and human land-use activities. Our project focuses on quantifying the role of thermokarst lake expansion, drainage, and drained lake basin evolution in the Arctic System. We did this through a combination of field studies, environmental sensor networks, remote sensing, and modeling. This dataset consists of environmental sensor records that record temperature and water level at three lakes in the Bugeye Lakes Complex in 2021 and 2022. Onset HOBO water level loggers (U2OL-04) were deployed in Bugeye Lakes 1, 2, and 4 in 2021 and 2022 to record temperature and pressure changes at 30-minute intervals. Water level was determined in Onset Hoboware Pro v. 3.7.23 using the barometric compensation assistant based on pressure transducer measurements below the water and from nearby atmospheric pressure measurements from a local pressure transducer mounted to pole on the tundra. The sensor data capture the partial drainage of Bugeye Lakes 2 and 4 into Bugeye Lake 1. Bugeye Lakes 1 and 2 also measure snow dam outburst flooding associated with the recently drained lake basins. 
    more » « less