Recovering Mantle Memories from River Profiles

Researchers use a closed-loop modeling strategy to validate regional uplift patterns recorded in river profiles across the African continent.

Short communication: Forward and inverse models relating river long profile to monotonic step-changes in tectonic rock uplift rate history: A theoretical perspective under a nonlinear slope-erosion dependency

The long profile of rivers is widely considered as a recorded of tectonic uplift rate. Knickpoints form in response to rate changes and faster rates produce steeper channel segments. However, when the exponent relating fluvial incision to river slope, n, is not unity, the links between tectonic rates and channel profile are complicated by channel dynamics that consume and form river segments. Here, we explore non-linear cases leading to channel segment consumption and develop a Lagrangian analytic model for knickpoint migration. We derive a criterion for knickpoint preservation and merging, and develop a forward analytic model that resolves knickpoint and long profile evolution before and after knickpoint merging. We further propose a linear inverse scheme to infer tectonic history from river profiles when all knickpoints are preserved. Our description provides a new framework to explore the links between tectonic uplift rates and river profile evolution when n is not unity.

Post-flood field investigation of the June 2020 flash flood in the upper Muráň River basin and the catastrophic flash flood scenario

After a dry spring, in June 2020 several intense storms occurred at the headwaters of the small basins of the Muráň and Zdychava rivers in the territory of the Muránska planina National Park (Slaná River basin, Slovakia). In the first part of the study– according to a hydrological survey made by the authors after the flash flood – the peak discharge was reconstructed at several Muráň River profiles. Next, the flash flood waves were reconstructed by the rainfall-runoff model NLC (non-linear cascade). The results of modelling based on field investigations show that, despite the extreme precipitation event (108 mm per 1 hour at the precipitation gauging station at Predná Hora), the peak flow rates were not exceptional in selected profiles on the Muráň River. The fact that extreme precipitation above 120 mm fell in a relatively very small area at the division of the Muráň and Zdychava rivers’ water contributed to this result. In the second part, a catastrophic 1000-year rainfall event scenario on the Zdychava River basin has been prepared. In analysing time series and identifying T-year daily rainfall depths, daily data was used from six precipitation stations in the vicinity of Muránska planina. Then, the 1000-year discharge of the Zdychava at Revúca was simulated by the calibrated NLC model. In such an extreme precipitation scenario, the peak flow rate would reach 105.15 m3 s−1, i.e. with a specific runoff of 1.78 m3 s−1 km–2. The total runoff in an 18-hour period would be 1.119 million m3, representing 21.11% of the rainfall (5.301 million m3).

Emergent simplicity despite local complexity in eroding fluvial landscapes

Much of our current understanding of continental topographic evolution is rooted in measuring and predicting the rates at which rivers erode the landscape. Flume tank and field observations indicate that stochasticity and local conditions play important roles in determining rates at small scales (e.g., <10 km, thousands of years). Obversely, preserved river profiles and common shapes of rivers atop uplifting topography indicate that erosion rates are predictable at larger scales. These observations indicate that the response of rivers to forcing can be scale dependent. I demonstrate that erosional thresholds can provide an explanation for why profile evolution can be very complicated and unique at small scales yet simple and predictable at large scales.

Emergent simplicity despite local complexity in eroding fluvial landscapes

Much understanding of continental topographic evolution is rooted in measuring and predicting rates at which rivers erode. Flume tank and field observations indicate that stochasticity and local conditions play important roles in determining rates at small scales (e.g. < 10 km, thousands of years). Obversely, preserved river profiles and common shapes of rivers atop uplifting topography indicate that erosion rates are predictable at larger scales. These observations indicate that the response of rivers to forcing can be scale dependent. Here I demonstrate that erosional thresholds can provide an explanation for why profile evolution can be very complicated and unique at small scales yet simple and predictable at large scales.

Quantifying Normal Fault Evolution from River Profile Analysis in the Northern Basin and Range Province, Southwest Montana, USA

Over the past few decades, tectonic geomorphology has been widely implemented to constrain spatial and temporal patterns of fault slip, especially where existing geologic or geodetic data are poor. We apply this practice along the eastern margin of Bull Mountain, Southwest Montana, where 15 transient channels are eroding into the flat, upstream relict landscape in response to an ongoing period of increased base level fall along the Western North Boulder fault. We aim to improve constraints on the spatial and temporal slip rates across the Western North Boulder fault zone by applying channel morphometrics, cosmogenic erosion rates, bedrock characteristics, and calibrated reproductions of the modern river profiles using a 1-dimensional stream power incision model that undergoes a change in the rate of base level fall. We perform over 104 base level fall simulations to explore a wide range of fault slip dynamics and stream power parameters. Our best fit simulations suggest that the Western North Boulder fault started as individual fault segments along the middle to southern regions of Bull Mountain that nucleated around 6.2 to 2.5 Ma, respectively. This was followed by the nucleation of fault segments in the northern region around 1.5 to 0.4 Ma. We recreate the evolution of the Western North Boulder fault to show that through time, these individual segments propagate at the fault tips and link together to span over 40 km, with a maximum slip of 462 m in the central portion of the fault. Fault slip rates range from 0.02 to 0.45 mm/yr along strike and are consistent with estimates for other active faults in the region. We find that the timing of fault initiation coincides well with the migration of the Yellowstone hotspot across the nearby Idaho-Montana border and thus attribute the initiation of extension to the crustal bulge from the migrating hotspot. Overall, we provide the first quantitative constraints on fault initiation and evolution of the Western North Boulder fault, perhaps the farthest north basin in the Northern Basin and Range province that such constraints exist. We show that river profiles are powerful tools for documenting the spatial and temporal patterns of normal fault evolution, especially where other geologic/geodetic methods are limited, proving to be a vital tool for accurate tectonic hazard assessments.

The ups and downs of the Missouri River from Pleistocene to present: Impact of climatic change and forebulge migration on river profiles, river course, and valley fill complexity

The Missouri River is a continent-scale river that has thus far escaped a rigorous reporting of valley fill trends within its trunk system. This study summarizes evolution of the lower Missouri River profile from the time of outwash in the Last Glacial Maximum (LGM) until establishment of the modern dominantly precipitation-fed river. This work relies on optically stimulated luminescence (OSL) dating, water-well data, and a collection of surficial geological maps of the valley compiled from U.S. Geological Survey EDMAP and National Science Foundation Research Experience for Undergrads projects. Mapping reveals five traceable surfaces within valley fill between Yankton, South Dakota, USA, and Columbia, Missouri, USA, that record two cycles of incision and aggradation between ca. 23 ka and ca. 8 ka. The river aggraded during the LGM to form the Malta Bend surface by ca. 26 ka. The Malta Bend surface is buried and fragmented but presumed to record a braided outwash plain. The Malta Bend surface was incised up to 18 m between ca. 23 ka and ca. 16 ka to form the Carrolton surface (ca. 16 ka to ca. 14 ka). The Carrollton surface ghosts a braided outwash morphology locally through overlying mud. Aggradation followed (ca. 14 ka to ca. 13.5 ka) to within 4 m of the modern floodplain surface and generated the Salix surface (ca. 13.5 to ca. 12 ka). By Salix time, the Missouri River was no longer an outwash river and formed a single-thread meandering pattern. Reincision at ca. 12 ka followed Salix deposition to form the short-lived Vermillion surface at approximately the grade of the earlier Carrolton surface. Rapid aggradation from ca. 10 ka to ca. 8 ka followed and formed the modern Omaha surface (ca. 8 ka to Present). The higher Malta Bend and Omaha profiles are at roughly the same grade, as are the lower Carrolton and Vermillion surfaces. The Salix surface is in between. All surfaces converge downstream as they enter the narrow and shallow bedrock valley just before reaching Columbia, Missouri. The maximum departure of the profiles is 18 m near Sioux City, Iowa, USA, at ∼100 km downstream from the James Lobe glacial input near Yankton, South Dakota. Incision and aggradation appear to be driven by relative changes in input of sediment and water related to glacial advance and retreat and then later by climatic changes near the Holocene transition. The incision from the Malta Bend to the Carrolton surface records the initial breakdown of the cryosphere at the end of the LGM, and this same incisional event is found in both the Ohio and Mississippi valleys. This incisional event records a “big wash” that resulted in the evacuation of sediment from each of the major outwash rivers of North America. The direction and magnitude of incision from the LGM to the modern does not fit with modeled glacioisostatic adjustment trends for the Missouri Valley. Glaciotectonics likely influenced the magnitude of incision and aggradation secondarily but does not appear to have controlled the overall timing or magnitude of either. Glaciotectonic valley tilting during the Holocene, however, did likely cause the Holocene channel to consistently migrate away from the glacial front, which argues for a forebulge axis south of the Missouri Valley during the Holocene and, by inference, earlier. This is at least 200 km south of where models predict the Holocene forebulge axis. The Missouri Valley thus appears to reside in the tectonic low between the ice front and the forebulge crest. The buffer valley component of incision caused by profile variation could explain as much as 25 m of the total ∼40 m of valley incision at Sioux City, Iowa. The Missouri Valley also hosted a glacial lobe as far south as Sioux City, Iowa, in pre-Wisconsinan time, which is also a factor in valley excavation.

Uplift and growth of the northwest Pamir

&lt;p&gt;The Pamir forms the northwestern tail of the Tibetan plateau and is a first-order tectonic feature of the Cenozoic Indo-Eurasian collision. The nature of the topographic uplift and orogenic growth of the entire northwestern margin of the Pamir is poorly constrained; however, this history can provide important constraints that are required to test geodynamic models of the tectonic evolution of the Pamir. Here we focus on the uplift history of the western and northwestern unglaciated margin of the Northern Pamir, the Darvaz and the Peter-the-First Ranges. These three ranges were formed by three major fault systems: the Main Pamir Thrust (MPT), the Darvaz and the Vakhsh fault zones (DFZ, VFZ). To assess the impact of tectonic uplift on the geomorphic evolution, we analyzed geomorphic characteristics of the topography, the longitudinal river profiles and the relief. To better constrain the regional crustal cooling history and uplift, we obtained thermochronologic cooling ages from the three regions.&lt;/p&gt;&lt;p&gt;We present 19 new zircon (U-Th-Sm)/He (ZHe) ages, 7 apatite fission track (AFT) ages, and 4 apatite (U-Th-Sm)/He (AHe) ages, ranging between &gt;200 and 4 Ma, 14 and 4 Ma, and 15 and 3 Ma, respectively. The three units are characterized by unique Neogene cooling pathways, suggesting that they exhumed independently.&lt;/p&gt;&lt;p&gt;We discovered extensive low-relief landscapes with Neogene sedimentary cover uplifted ~2 km in elevation above the present-day regional base level. Our analysis indicates that the Panj and Vakhsh rivers form the regional base levels for the river network draining the entire northern and western margin of the Pamir. In the hanging wall of DFZ, the Paleozoic bedrock is characterized by significant relief (&gt;1 km), the Neogene cover onlaps directly onto this Paleozoic bedrock. The tributary rivers crossing these landscapes are characterized by gentle, concave upstream longitudinal profiles at high elevation. These are interrupted by major knickpoint zones and steep downstream segments draining towards the deeply incised Panj and Vakhsh rivers. This indicates that the Darvaz Fault hanging wall had been uplifted and eroded prior to deposition of upper Neogene sediments, suggesting that the DFZ has a prolonged Neogene slip history. In contrast to the northeastern Pamir, here, the MPT-hanging-wall is characterized by reset late Oligocene-Early Miocene ZHe cooling ages ranging between 26 and 17 Ma. AFT and AHe-ages between 15 and 13 Ma suggest that exhumation suddenly terminated during the middle Miocene. In contrast, Jurassic sandstones exposed near the DFZ yield mostly un-reset Triassic-Jurassic ZHe ages (~250-170 Ma), a reset AFT age of ~5 Ma and a 2.5 Ma AHe age. Within the Peter-the-1st-Range, we obtained fully reset ~ 5 Ma ZHe ages, and ~4 Ma AFT ages. The rapid cooling trends since at least ~5 Ma suggest that deformation and a significant portion of crustal shortening propagated into the Tadjik foreland basin, causing enhanced uplift and erosion of the hanging wall of the VFZ and related faults. This deformation triggered ~2 km uplift of the entire northwest Pamir, recorded in uplifted paleo-landscapes and dissected tributaries of the Panj and Vakhsh rivers.&lt;/p&gt;