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   
                    
                            
                            Drainage network reorganization in landscapes buried by glacial deposits
                        
                    
    
            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   
        
    
                            - Award ID(s):
- 2052938
- PAR ID:
- 10509062
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            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
- 
            Abstract Erosional perturbations from changes in climate or tectonics are recorded in the profiles of bedrock rivers, but these signals can be challenging to unravel in settings with non‐uniform lithology. In layered rocks, the surface lithology at a given location varies through time as erosion exposes different layers of rock. Recent modeling studies have used the Stream Power Model (SPM) to highlight complex variations in erosion rates that arise in bedrock rivers incising through layered rocks. However, these studies do not capture the effects of coarse sediment cover on channel evolution. We use the “Stream Power with Alluvium Conservation and Entrainment” (SPACE) model to explore how sediment cover influences landscape evolution and modulates the topographic expression of erodibility contrasts in horizontally layered rocks. We simulate river evolution through alternating layers of hard and soft rock over million‐year timescales with a constant and uniform uplift rate. Compared to the SPM, model runs with sediment cover have systematically higher channel steepness values in soft rock layers and lower channel steepness values in hard rock layers. As more sediment accumulates, the contrast in steepness between the two rock types decreases. Effective bedrock erodibilities back‐calculated assuming the SPM are strongly influenced by sediment cover. We also find that sediment cover can significantly increase total relief and timescales of adjustment toward landscape‐averaged steady‐state topography and erosion rates.more » « less
- 
            Abstract Contrasts in bedrock erodibility have been shown to drive landscape transience, but it is unclear whether horizontal tectonic displacements would enhance such effects. Furthermore, one might expect these factors to coexist, as tectonic convergence helps to create rock strength contrasts in settings like the Himalayas. How do landscapes respond when contacts separating units are raised vertically and shifted horizontally by tectonics? To evaluate such questions, we use landscape evolution models to simulate the exposure of a weak unit in a landscape equilibrated to a strong unit. We explore different simulations varying factors like weak unit erodibility, diffusivity, contact dip, and topographic advection rate. In these simulations, we assess the migration of the main drainage divide as well as changes in channel steepness and topographic relief within the strong unit. Our model results show that the horizontal movement of a contact does enhance drainage divide migration and increases in channel steepness, especially when the contact migrates along rivers with low drainage areas. Across all simulations, however, increases in topographic relief are minimal and temporary. Unexpected behaviors emerge in our simulations in which the mass balance of topography is influenced by horizontal tectonic displacements. For example, the exposure of the weak unit causes a gradual decline in the steepness of the strong unit. We interpret such behaviors to be artifacts related to the fixed boundaries of our domain and likely unrepresentative of natural landscapes. Instead, we focus on simulations where advection does not influence the mass balance of topography. These models show that the horizontal movement of contacts can enhance landscape transience, but this transience is marked by features one can use as diagnostic characteristics. Detecting such characteristics in natural landscapes featuring tectonic convergence would be difficult, however, due to the natural coincidence of factors such as faulting, folding, and landslides.more » « less
- 
            Ancient glaciated landscapes reveal interactions among ice dynamics, bed composition, and preglacial river networks. Subglacial landscapes, revealed in regions of recent ice-sheet retreat, provide a window into ice-sheet dynamics and interactions with evolving subglacial topography. Here, we document landscape evolution beneath the southern Laurentide Ice Sheet of North America since the end of the Pliocene, 2.6 million years (Ma) ago, by reconstructing the isostatically adjusted preglacial surface and modern bedrock topography at 250 m horizontal resolution. We use flow routing to reconstruct drainage networks and river longitudinal profiles, revealing the pattern and extent of their glacially forced reorganization. The overall mean Quaternary (2.6 Ma ago to present) erosion rate is 27 m/Ma, rising within ice-streaming corridors to 35 m/Ma (and locally reaching 400 m/Ma) and falling to 22 m/Ma in non–ice-streaming regions. Our results suggest that subglacial erosion was sufficient to lower the southern Laurentide Ice Sheet into warmer environments, thereby enhancing ablation and reducing ice-sheet extent over time.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    