A vector-based method for bank-material tracking in coupled models of meandering and landscape evolution

Sinuous channels commonly migrate laterally and interact with banks of different strengths—an interplay that links geomorphology and life and shapes diverse landscapes from the seafloor to planetary surfaces. To investigate feedbacks between meandering rivers and landscapes over geomorphic timescales, numerical models typically represent bank properties using grids; however, this approach produces results inherently dependent on grid resolution. Herein we assess existing techniques for tracking landscape and bank-strength evolution in numerical models of meandering channels and show that grid-based models implicitly include unintended thresholds for bank migration that can control simulated landscape evolution. Building on stratigraphic modeling techniques, we develop a vector-based method for land surface- and subsurface-material tracking that overcomes the resolution-dependence inherent in grid-based techniques by allowing high-fidelity representation of bank-material properties for curvilinear banks and low channel lateral migration rates. We illustrate four specific applications of the new technique: (1) the effect of resistant mud-rich deposits in abandoned meander cutoff loops on meander belt evolution; (2) the stratigraphic architecture of aggrading, alluvial meandering channels that interact with cohesive-bank and floodplain material; (3) the evolution of an incising, meandering river with mixed bedrock and alluvial banks within a confined bedrock valley; and (4) the effect of a bank-height dependent lateral-erosion rate for a meandering river in an aggrading floodplain. In all cases the vector-based approach overcomes numerical artifacts with the grid-based model. Because of its geometric flexibility, the vector-based material tracking approach provides new opportunities for exploring the coevolution of meandering rivers and surrounding landscapes over geologic timescales

Numerical simulations of bedrock valley evolution by meandering rivers with variable bank material

Bedrock river valleys are fundamental components of many landscapes, and their morphologies—from slot canyons with incised meanders to wide valleys with strath terraces—may record environmental history. Several formation mechanisms for particular valley types have been proposed that involve changes in climatic and tectonic forcing, but the uniqueness of valley evolution pathways and the long-term stability of valley morphology under constant forcing are unknown and are not predicted in existing numerical models for vertically incising rivers. Because rivers often migrate more rapidly through alluvium than through bedrock, we explore the hypothesis that the distribution of bank materials strongly influences river meandering kinematics and can explain the diversity of bedrock river valley morphology. Simulations using a numerical model of river meandering with vector-based bank-material tracking indicate that channel lateral erosion rate in sediment and bedrock, vertical erosion rate, and initial alluvial-belt width explain first-order differences in bedrock valley type; that bedrock-bound channels can evolve under steady forcing from alluvial states; and that weak bedrock and low vertical incision rates favor wide, shallow valleys, while resistant bedrock and high vertical incision rates favor narrow, deep valleys. During vertical incision, sustained planation of the valley floor is favored when bedrock boundaries restrict channel migration to a zone of thin sediment fill. The inherent unsteadiness of river meandering in space and time is enhanced by evolving spatial contrasts in bank strength between sediment and bedrock and can account for several valley features—including strath terraces and underfit valleys—commonly ascribed to external drivers

Numerical model predictions of autogenic fluvial terraces and comparison to climate change expectations

Terraces eroded into sediment (alluvial) and bedrock (strath) preserve an important history of river activity. River terraces are thought to form when a river switches from a period of slow vertical incision and valley widening to fast vertical incision and terrace abandonment. Consequently, terraces are often interpreted to reflect changing external drivers including tectonics, sea level, and climate. In contrast, the intrinsic unsteadiness of lateral migration in rivers may generate terraces even under constant rates of vertical incision without external forcing. To explore this mechanism, we simulate landscape evolution by a vertically incising, meandering river and isolate the age and geometry of autogenic river terraces. Modeled autogenic terraces form for a wide range of lateral and vertical incision rates and are often paired and longitudinally extensive for intermediate ratios of vertical-to-lateral erosion rate. Autogenic terraces have a characteristic reoccurrence time that scales with the time for relief generation. There is a preservation bias against older terraces due to reworking of previously visited parts of the valley. Evolving, spatial differences in bank strength between bedrock and sediment reduce terrace formation frequency and length, favor pairing, and can explain sublinear terrace margins at valley boundaries. Age differences and geometries for modeled autogenic terraces are consistent, in cases, with natural terraces and overlap with metrics commonly attributed to terrace formation due to climate change. We suggest a new phase space of terrace properties that may allow differentiation of autogenic terraces from terraces formed by external drivers

On the morphodynamics of a wide class of large-scale meandering rivers: Insights gained by coupling LES with sediment-dynamics

In meandering rivers, interactions between flow, sediment transport, and bed topography affect diverse processes, including bedform development and channel migration. Predicting how these interactions affect the spatial patterns and magnitudes of bed deformation in meandering rivers is essential for various river engineering and geoscience problems. Computational fluid dynamics simulations can predict river morphodynamics at fine temporal and spatial scales but have traditionally been challenged by the large scale of natural rivers. We conducted coupled large-eddy simulation (LES) and bed morphodynamics simulations to create a unique database of hydro-morphodynamic datasets for 42 meandering rivers with a variety of planform shapes and large-scale geometrical features that mimic natural meanders. For each simulated river, the database includes (i) bed morphology, (ii) three-dimensional mean velocity field, and (iii) bed shear stress distribution under bankfull flow conditions. The calculated morphodynamics results at dynamic equilibrium revealed the formation of scour and deposition patterns near the outer and inner banks, respectively, while the location of point bars and scour regions around the apexes of the meander bends is found to vary as a function of the radius of curvature of the bends to the width ratio. A new mechanism is proposed that explains this seemingly paradoxical finding. The high-fidelity simulation results generated in this work provide researchers and scientists with a rich numerical database for morphodynamics and bed shear stress distributions in large-scale meandering rivers to enable systematic investigation of the underlying phenomena and support a range of river engineering applications

Model predictions of long-lived storage of organic carbon in river deposits

The mass of carbon stored as organic matter in terrestrial systems is sufficiently large to play an important role in the global biogeochemical cycling of CO_2 and O_2. Field measurements of radiocarbon-depleted particulate organic carbon (POC) in rivers suggest that terrestrial organic matter persists in surface environments over millennial (or greater) timescales, but the exact mechanisms behind these long storage times remain poorly understood. To address this knowledge gap, we developed a numerical model for the radiocarbon content of riverine POC that accounts for both the duration of sediment storage in river deposits and the effects of POC cycling. We specifically target rivers because sediment transport influences the maximum amount of time organic matter can persist in the terrestrial realm and river catchment areas are large relative to the spatial scale of variability in biogeochemical processes. Our results show that rivers preferentially erode young deposits, which, at steady state, requires that the oldest river deposits are stored for longer than expected for a well-mixed sedimentary reservoir. This geometric relationship can be described by an exponentially tempered power-law distribution of sediment storage durations, which allows for significant aging of biospheric POC. While OC cycling partially limits the effects of sediment storage, the consistency between our model predictions and a compilation of field data highlights the important role of storage in setting the radiocarbon content of riverine POC. The results of this study imply that the controls on the terrestrial OC cycle are not limited to the factors that affect rates of primary productivity and respiration but also include the dynamics of terrestrial sedimentary systems

The Entropic Braiding Index (eBI): A Robust Metric to Account for the Diversity of Channel Scales in Multi-Thread Rivers

The Braiding Index (BI), defined as the average count of intercepted channels per cross-section, is a widely used metric for characterizing multi-thread river systems. However, it does not account for the diversity of channels (e.g., in terms of water discharge) within different cross-sections, omitting important information related to system complexity. Here we present a modification of BI, the Entropic Braiding Index (eBI), which augments the information content in BI by using Shannon Entropy to encode the diversity of channels in each cross section. eBI is interpreted as the number of “effective channels” per cross-section, allowing a direct comparison with the traditional BI. We demonstrate the potential of the ratio BI/eBI to quantify channel disparity, differentiate types of multi-thread systems (braided vs. anastomosed), and assess the effect of discharge variability, such as seasonal flooding, on river cross-section stability

Detailed stratigraphy and bed thickness of the Mars north and south polar layered deposits

The Mars polar layered deposits (PLD) likely hold an extensive record of recent climate during a period of high-amplitude orbit and obliquity cycles. Previous work has detected limited evidence for orbital signatures within PLD stratigraphy, but data from the High Resolution Imaging Science Experiment (HiRISE) permit renewed analysis of PLD stratigraphy at sub-meter scale. Topography derived from HiRISE images using stereogrammetry resolves beds previously detectable only as alternating light and dark bands in visible images. We utilize these data to measure the thickness of individual beds within the PLD, corrected for non-horizontal bed orientation. Stratigraphic columns and bed thickness profiles are presented for two sites within the NPLD, and show several sets of finely bedded units 1–2 m thick; isolated marker beds 3–4 m thick; and undifferentiated sections. Bed thickness measurements for three sites within the SPLD exhibit only one bed type based on albedo and morphology, and bed thicknesses have a larger mean and variance compared to measurements for the NPLD. Power spectra of brightness and slope derived along the measured stratigraphic sections confirm the regularity of NPLD fine bed thickness, and the lack of a dominant SPLD bed thickness. The regularity of fine bed thickness of the NPLD is consistent with quasiperiodic bed formation, albeit with unknown temporal period; the SPLD thickness measurements show no such regularity

Topography and flow model files for the Platte River, Nebraska, 2016-2017

The data set includes raw and processed data for a digital elevation model; files for executing an open-source flow model; and maps of flow depth produced by this model.Automated techniques for extracting channels from topography are well developed for convergent channel networks, and identify flow paths based on land-surface gradients. These techniques—even when they allow multiple flow paths—do not consistently capture channel networks with frequent bifurcations (e.g., in rivers, deltas, and alluvial fans). The project uses multithread rivers as a template to develop a new approach for channel extraction suitable for channel networks with divergences. This data set includes topography and flow model data used to demonstrate this new approach for automated identification of channels on planetary surfaces. The data set accompanies a 2017 publication in the Journal of Geophysical Research: Earth Surface.St. Anthony Falls Laboratory Industrial ConsortiumNational Center for Earth-surface Dynamics