More extensive land loss expected on coastal deltas due to rivers jumping course during sea-level rise.

River deltas are home to hundreds of millions of people worldwide and are in danger of sinking due to anthropogenic sea-level rise, land subsidence, and reduced sediment supply. Land loss is commonly forecast by averaging river sediment supply across the entire delta plain to assess whether deposition can keep pace with sea-level rise. However, land loss and deposition vary across the landscape because rivers periodically jump course, rerouting sediment to distinct subregions called delta lobes. Here, we developed a model to forecast land loss that resolves delta lobes and tested the model against a scaled laboratory experiment. Both the model and the experiment show that rivers build land on the active lobe, but the delta incurs gradual land loss on inactive lobes that are cut off from sediment after the river abandons course. The result is a band of terrain along the coast that is usually drowned but is nonetheless a sink for sediment when the lobe is active, leaving less of the total sediment supply available to maintain persistent dry land. Land loss is expected to be more extensive than predicted by classical delta-plain-averaged models. Estimates for eight large deltas worldwide suggest that roughly half of the riverine sediment supply is delivered to terrain that undergoes long periods of submergence. These results draw the sustainability of deltas further into question and provide a framework to plan engineered diversions at a pace that will mitigate land loss in the face of rising sea levels.

Orbital and In-Situ Investigation of Periodic Bedrock Ridges in Glen Torridon, Gale Crater, Mars.

Gale crater, the field site for NASA's Mars Science Laboratory Curiosity rover, contains a diverse and extensive record of aeolian deposition and erosion. This study focuses on a series of regularly spaced, curvilinear, and sometimes branching bedrock ridges that occur within the Glen Torridon region on the lower northwest flank of Aeolis Mons, the central mound within Gale crater. During Curiosity's exploration of Glen Torridon between sols ∼2300-3080, the rover drove through this field of ridges, providing the opportunity for in situ observation of these features. This study uses orbiter and rover data to characterize ridge morphology, spatial distribution, compositional and material properties, and association with other aeolian features in the area. Based on these observations, we find that the Glen Torridon ridges are consistent with an origin as wind-eroded bedrock ridges, carved during the exhumation of Mount Sharp. Erosional features like the Glen Torridon ridges observed elsewhere on Mars, termed periodic bedrock ridges (PBRs), have been interpreted to form transverse to the dominant wind direction. The size and morphology of the Glen Torridon PBRs are consistent with transverse formative winds, but the orientation of nearby aeolian bedforms and bedrock erosional features raise the possibility of PBR formation by a net northeasterly wind regime. Although several formation models for the Glen Torridon PBRs are still under consideration, and questions persist about the nature of PBR-forming paleowinds, the presence of PBRs at this site provides important constraints on the depositional and erosional history of Gale crater.

Amplification of downstream flood stage due to damming of fine-grained rivers.

River dams provide many benefits, including flood control. However, due to constantly evolving channel morphology, downstream conveyance of floodwaters following dam closure is difficult to predict. Here, we test the hypothesis that the incised, enlarged channel downstream of dams provides enhanced water conveyance, using a case study from the lower Yellow River, China. We find that, although flood stage is lowered for small floods, counterintuitively, flood stage downstream of a dam can be amplified for moderate and large floods. This arises because bed incision is accompanied by sediment coarsening, which facilitates development of large dunes that increase flow resistance and reduce velocity relative to pre-dam conditions. Our findings indicate the underlying mechanism for such flood amplification may occur in >80% of fine-grained rivers, and suggest the need to reconsider flood control strategies in such rivers worldwide.

Where rivers jump course.

Rivers can abruptly shift pathways in rare events called avulsions, which cause devastating floods. The controls on avulsion locations are poorly understood as a result of sparse data on such features. We analyzed nearly 50 years of satellite imagery and documented 113 avulsions across the globe that indicate three distinct controls on avulsion location. Avulsions on fans coincide with valley-confinement change, whereas avulsions on deltas are primarily clustered within the backwater zone, indicating a control by spatial flow deceleration or acceleration during floods. However, 38% of avulsions on deltas occurred upstream of backwater effects. These events occurred in steep, sediment-rich rivers in tropical and desert environments. Our results indicate that avulsion location on deltas is set by the upstream extent of flood-driven erosion, which is typically limited to the backwater zone but can extend far upstream in steep, sediment-laden rivers. Our findings elucidate how avulsion hazards might respond to land use and climate change.

Glacial isostatic adjustment directed incision of the Channeled Scabland by Ice Age megafloods.

During the last deglaciation, dozens of glacial outburst floods-among the largest known floods on Earth-scoured the Channeled Scabland landscape of eastern Washington. Over this same period, deformation of the Earth's crust in response to the growth and decay of ice sheets changed the topography by hundreds of meters. Here, we investigated whether glacial isostatic adjustment affected routing of the Missoula floods and incision of the Channeled Scabland from an impounded, glacial Lake Columbia. We used modern topography corrected for glacial isostatic adjustment as an input to flood models that solved the depth-averaged, shallow water equations and compared the results to erosion constraints. Results showed that floods could have traversed and eroded parts of two major tracts of the Channeled Scabland-Telford-Crab Creek and Cheney-Palouse-near 18 ka, whereas glacial isostatic adjustment limited flow into the Cheney-Palouse tract at 15.5 ka. Partitioning of flow between tracts was governed by tilting of the landscape, which affected the filling and overspill of glacial Lake Columbia directly upstream of the tracts. These results highlight the impact of glacial isostatic adjustment on megaflood routing and landscape evolution.

Impact of River Channel Lateral Migration on Microbial Communities across a Discontinuous Permafrost Floodplain.

Permafrost soils store approximately twice the amount of carbon currently present in Earth's atmosphere and are acutely impacted by climate change due to the polar amplification of increasing global temperature. Many organic-rich permafrost sediments are located on large river floodplains, where river channel migration periodically erodes and redeposits the upper tens of meters of sediment. Channel migration exerts a first-order control on the geographic distribution of permafrost and floodplain stratigraphy and thus may affect microbial habitats. To examine how river channel migration in discontinuous permafrost environments affects microbial community composition, we used amplicon sequencing of the 16S rRNA gene on sediment samples from floodplain cores and exposed riverbanks along the Koyukuk River, a large tributary of the Yukon River in west-central Alaska. Microbial communities are sensitive to permafrost thaw: communities found in deep samples thawed by the river closely resembled near-surface active-layer communities in nonmetric multidimensional scaling analyses but did not resemble floodplain permafrost communities at the same depth. Microbial communities also displayed lower diversity and evenness in permafrost than in both the active layer and permafrost-free point bars recently deposited by river channel migration. Taxonomic assignments based on 16S and quantitative PCR for the methyl coenzyme M reductase functional gene demonstrated that methanogens and methanotrophs are abundant in older permafrost-bearing deposits but not in younger, nonpermafrost point bar deposits. The results suggested that river migration, which regulates the distribution of permafrost, also modulates the distribution of microbes potentially capable of producing and consuming methane on the Koyukuk River floodplain. IMPORTANCE Arctic lowlands contain large quantities of soil organic carbon that is currently sequestered in permafrost. With rising temperatures, permafrost thaw may allow this carbon to be consumed by microbial communities and released to the atmosphere as carbon dioxide or methane. We used gene sequencing to determine the microbial communities present in the floodplain of a river running through discontinuous permafrost. We found that the river's lateral movement across its floodplain influences the occurrence of certain microbial communities-in particular, methane-cycling microbes were present on the older, permafrost-bearing eroding riverbank but absent on the newly deposited river bars. Riverbank sediment had microbial communities more similar to those of the floodplain active-layer samples than permafrost samples from the same depth. Therefore, spatial patterns of river migration influence the distribution of microbial taxa relevant to the warming Arctic climate.

Coal fly ash is a major carbon flux in the Chang Jiang (Yangtze River) basin.

Fly ash-the residuum of coal burning-contains a considerable amount of fossilized particulate organic carbon (FOCash) that remains after high-temperature combustion. Fly ash leaks into natural environments and participates in the contemporary carbon cycle, but its reactivity and flux remained poorly understood. We characterized FOCash in the Chang Jiang (Yangtze River) basin, China, and quantified the riverine FOCash fluxes. Using Raman spectral analysis, ramped pyrolysis oxidation, and chemical oxidation, we found that FOCash is highly recalcitrant and unreactive, whereas shale-derived FOC (FOCrock) was much more labile and easily oxidized. By combining mass balance calculations and other estimates of fly ash input to rivers, we estimated that the flux of FOCash carried by the Chang Jiang was 0.21 to 0.42 Mt C⋅y-1 in 2007 to 2008-an amount equivalent to 37 to 72% of the total riverine FOC export. We attributed such high flux to the combination of increasing coal combustion that enhances FOCash production and the massive construction of dams in the basin that reduces the flux of FOCrock eroded from upstream mountainous areas. Using global ash data, a first-order estimate suggests that FOCash makes up to 16% of the present-day global riverine FOC flux to the oceans. This reflects a substantial impact of anthropogenic activities on the fluxes and burial of fossil organic carbon that has been made less reactive than the rocks from which it was derived.

Early plant organics increased global terrestrial mud deposition through enhanced flocculation.

An irreversible increase in alluvial mudrock occurred with the Ordovician-Silurian evolution of bryophytes, challenging a paradigm that deep-rooted plants were responsible for this landscape shift. We tested the idea that increased primary production and plant organics promoted aggregation of clay into flocs in rivers and facilitated mud deposition on floodplains. In experiments, we observed that clay readily flocculated for organic and clay concentrations common to modern rivers, yielding settling velocities three orders of magnitude larger than those without organics. Using a transport model, we found that flocculation substantially increased mud deposition, resulting in muddier floodplains. Thus, organic-induced flocculation may have been more critical than deep-rooted plants in the proliferation of muddy floodplains.

Accelerated river avulsion frequency on lowland deltas due to sea-level rise.

Sea-level rise, subsidence, and reduced fluvial sediment supply are causing river deltas to drown worldwide, affecting ecosystems and billions of people. Abrupt changes in river course, called avulsions, naturally nourish sinking land with sediment; however, they also create catastrophic flood hazards. Existing observations and models conflict on whether the occurrence of avulsions will change due to relative sea-level rise, hampering the ability to forecast delta response to global climate change. Here, we combined theory, numerical modeling, and field observations to develop a mechanistic framework to predict avulsion frequency on deltas with multiple self-formed lobes that scale with backwater hydrodynamics. Results show that avulsion frequency is controlled by the competition between relative sea-level rise and sediment supply that drives lobe progradation. We find that most large deltas are experiencing sufficiently low progradation rates such that relative sea-level rise enhances aggradation rates-accelerating avulsion frequency and associated hazards compared to preindustrial conditions. Some deltas may face even greater risk; if relative sea-level rise significantly outpaces sediment supply, then avulsion frequency is maximized, delta plains drown, and avulsion locations shift inland, posing new hazards to upstream communities. Results indicate that managed deltas can support more frequent engineered avulsions to recover sinking land; however, there is a threshold beyond which coastal land will be lost, and mitigation efforts should shift upstream.

Universal relation with regime transition for sediment transport in fine-grained rivers.

Fine-grained sediment (grain size under 2,000 µm) builds floodplains and deltas, and shapes the coastlines where much of humanity lives. However, a universal, physically based predictor of sediment flux for fine-grained rivers remains to be developed. Herein, a comprehensive sediment load database for fine-grained channels, ranging from small experimental flumes to megarivers, is used to find a predictive algorithm. Two distinct transport regimes emerge, separated by a discontinuous transition for median bed grain size within the very fine sand range (81 to 154 µm), whereby sediment flux decreases by up to 100-fold for coarser sand-bedded rivers compared to river with silt and very fine sand beds. Evidence suggests that the discontinuous change in sediment load originates from a transition of transport mode between mixed suspended bed load transport and suspension-dominated transport. Events that alter bed sediment size near the transition may significantly affect fluviocoastal morphology by drastically changing sediment flux, as shown by data from the Yellow River, China, which, over time, transitioned back and forth 3 times between states of high and low transport efficiency in response to anthropic activities.

Low-gradient, single-threaded rivers prior to greening of the continents.

The Silurian-age rise of land plants is hypothesized to have caused a global revolution in the mechanics of rivers. In the absence of vegetation-controlled bank stabilization effects, pre-Silurian rivers are thought to be characterized by shallow, multithreaded flows, and steep river gradients. This hypothesis, however, is at odds with the pancontinental scale of early Neoproterozoic river systems that would have necessitated extraordinarily high mountains if such river gradients were commonplace at continental scale, which is inconsistent with constraints on lithospheric thickness. To reconcile these observations, we generated estimates of paleogradients and morphologies of pre-Silurian rivers using a well-developed quantitative framework based on the formation of river bars and dunes. We combined data from previous work with original field measurements of the scale, texture, and structure of fluvial deposits in Proterozoic-age Torridonian Group, Scotland-a type-example of pancontinental, prevegetation fluvial systems. Results showed that these rivers were low sloping (gradients 10-5 to 10-4), relatively deep (4 to 15 m), and had morphology similar to modern, lowland rivers. Our results provide mechanistic evidence for the abundance of low gradient, single-threaded rivers in the Proterozoic eon, at a time well before the evolution and radiation of land plants-despite the absence of muddy and vegetated floodplains. Single-threaded rivers with stable floodplains appear to have been a persistent feature of our planet despite singular changes in its terrestrial biota.

Self-formed bedrock waterfalls.

Waterfalls are inspiring landforms that set the pace of landscape evolution as a result of bedrock incision1-3. They communicate changes in sea level or tectonic uplift throughout landscapes2,4 or stall river incision, disconnecting landscapes from downstream perturbations3,5. Here we use a flume experiment with constant water discharge and sediment feed to show that waterfalls can form from a planar, homogeneous bedrock bed in the absence of external perturbations. In our experiment, instabilities between flow hydraulics, sediment transport and bedrock erosion lead to undulating bedforms, which grow to become waterfalls. We propose that it is plausible that the origin of some waterfalls in natural systems can be attributed to this intrinsic formation process and we suggest that investigations to distinguish self-formed from externally forced waterfalls may help to improve the reconstruction of Earth history from landscapes.

Abrupt drainage basin reorganization following a Pleistocene river capture.

River capture is a dramatic natural process of internal competition through which mountainous landscapes evolve and respond to perturbations in tectonics and climate. River capture may occur when one river network grows at the expense of another, resulting in a victor that steals the neighboring headwaters. While river capture occurs regularly in numerical models, field observations are rare. Here we document a late Pleistocene river capture in the Yimeng Mountains, China that abruptly shifted 25 km2 of drainage area from one catchment to another. River terraces and imbricated cobbles indicate that the main channel incised 27 m into granitic bedrock within 80 kyr, following the capture event, and upstream propagating knickpoints and waterfalls reversed the flow direction of a major river. Topographic analysis shows that the capture shifted the river basins far from topographic equilibrium, and active divide migration is propagating the effects of the capture throughout the landscape.

Similarity of stream width distributions across headwater systems.

The morphology and abundance of streams control the rates of hydraulic and biogeochemical exchange between streams, groundwater, and the atmosphere. In large river systems, the relationship between river width and abundance is fractal, such that narrow rivers are proportionally more common than wider rivers. However, in headwater systems, where many biogeochemical reactions are most rapid, the relationship between stream width and abundance is unknown. To constrain this uncertainty, we surveyed stream hydromorphology (wetted width and length) in several headwater stream networks across North America and New Zealand. Here, we find a strikingly consistent lognormal statistical distribution of stream width, including a characteristic most abundant stream width of 32 ± 7 cm independent of discharge or physiographic conditions. We propose a hydromorphic model that can be used to more accurately estimate the hydromorphology of streams, with significant impact on the understanding of the hydraulic, ecological, and biogeochemical functions of stream networks.

Time scale bias in erosion rates of glaciated landscapes.

Deciphering erosion rates over geologic time is fundamental for understanding the interplay between climate, tectonic, and erosional processes. Existing techniques integrate erosion over different time scales, and direct comparison of such rates is routinely done in earth science. On the basis of a global compilation, we show that erosion rate estimates in glaciated landscapes may be affected by a systematic averaging bias that produces higher estimated erosion rates toward the present, which do not reflect straightforward changes in erosion rates through time. This trend can result from a heavy-tailed distribution of erosional hiatuses (that is, time periods where no or relatively slow erosion occurs). We argue that such a distribution can result from the intermittency of erosional processes in glaciated landscapes that are tightly coupled to climate variability from decadal to millennial time scales. In contrast, we find no evidence for a time scale bias in spatially averaged erosion rates of landscapes dominated by river incision. We discuss the implications of our findings in the context of the proposed coupling between climate and tectonics, and interpreting erosion rate estimates with different averaging time scales through geologic time.

Progressive incision of the Channeled Scablands by outburst floods.

The surfaces of Earth and Mars contain large bedrock canyons that were carved by catastrophic outburst floods. Reconstructing the magnitude of these canyon-forming floods is essential for understanding the ways in which floods modify planetary surfaces, the hydrology of early Mars and abrupt changes in climate. Flood discharges are often estimated by assuming that the floods filled the canyons to their brims with water; however, an alternative hypothesis is that canyon morphology adjusts during incision such that bed shear stresses exceed the threshold for erosion by a small amount. Here we show that accounting for erosion thresholds during canyon incision results in near-constant discharges that are five- to ten-fold smaller than full-to-the-brim estimates for Moses Coulee, a canyon in the Channeled Scablands, which was carved during the Pleistocene by the catastrophic Missoula floods in eastern Washington, USA. The predicted discharges are consistent with flow-depth indicators from gravel bars within the canyon. In contrast, under the assumption that floods filled canyons to their brims, a large and monotonic increase in flood discharge is predicted as the canyon was progressively incised, which is at odds with the discharges expected for floods originating from glacial lake outbursts. These findings suggest that flood-carved landscapes in fractured rock might evolve to a threshold state for bedrock erosion, thus implying much lower flood discharges than previously thought.

Experimental river delta size set by multiple floods and backwater hydrodynamics.

River deltas worldwide are currently under threat of drowning and destruction by sea-level rise, subsidence, and oceanic storms, highlighting the need to quantify their growth processes. Deltas are built through construction of sediment lobes, and emerging theories suggest that the size of delta lobes scales with backwater hydrodynamics, but these ideas are difficult to test on natural deltas that evolve slowly. We show results of the first laboratory delta built through successive deposition of lobes that maintain a constant size. We show that the characteristic size of delta lobes emerges because of a preferential avulsion node-the location where the river course periodically and abruptly shifts-that remains fixed spatially relative to the prograding shoreline. The preferential avulsion node in our experiments is a consequence of multiple river floods and Froude-subcritical flows that produce persistent nonuniform flows and a peak in net channel deposition within the backwater zone of the coastal river. In contrast, experimental deltas without multiple floods produce flows with uniform velocities and delta lobes that lack a characteristic size. Results have broad applications to sustainable management of deltas and for decoding their stratigraphic record on Earth and Mars.

Quantitative bounds on morphodynamics and implications for reading the sedimentary record.

Sedimentary rocks are the archives of environmental conditions and ancient planetary surface processes that led to their formation. Reconstructions of Earth's past surface behaviour from the physical sedimentary record remain controversial, however, in part because we lack a quantitative framework to deconvolve internal dynamics of sediment-transport systems from environmental signal preservation. Internal dynamics of landscapes--a consequence of the coupling between bed topography, sediment transport and flow dynamics (morphodynamics)--result in regular and quasiperiodic landforms that abound on the Earth and other planets. Here, using theory and a data compilation of morphodynamic landforms that span a wide range of terrestrial, marine and planetary depositional systems, we show that the advection length for settling sediment sets bounds on the scales over which internal landscape dynamics operate. These bounds provide a universal palaeohydraulic reconstruction tool on planetary surfaces and allow for quantitative identification of depositional systems that may preserve tectonic, climatic and anthropogenic signals.

Amphitheater-headed canyons formed by megaflooding at Malad Gorge, Idaho.

Many bedrock canyons on Earth and Mars were eroded by upstream propagating headwalls, and a prominent goal in geomorphology and planetary science is to determine formation processes from canyon morphology. A diagnostic link between process and form remains highly controversial, however, and field investigations that isolate controls on canyon morphology are needed. Here we investigate the origin of Malad Gorge, Idaho, a canyon system cut into basalt with three remarkably distinct heads: two with amphitheater headwalls and the third housing the active Wood River and ending in a 7% grade knickzone. Scoured rims of the headwalls, relict plunge pools, sediment-transport constraints, and cosmogenic ((3)He) exposure ages indicate formation of the amphitheater-headed canyons by large-scale flooding ∼46 ka, coeval with formation of Box Canyon 18 km to the south as well as the eruption of McKinney Butte Basalt, suggesting widespread canyon formation following lava-flow diversion of the paleo-Wood River. Exposure ages within the knickzone-headed canyon indicate progressive upstream younging of strath terraces and a knickzone propagation rate of 2.5 cm/y over at least the past 33 ka. Results point to a potential diagnostic link between vertical amphitheater headwalls in basalt and rapid erosion during megaflooding due to the onset of block toppling, rather than previous interpretations of seepage erosion, with implications for quantifying the early hydrosphere of Mars.

Landslide-dammed paleolake perturbs marine sedimentation and drives genetic change in anadromous fish.

Large bedrock landslides have been shown to modulate rates and processes of river activity by forming dams, forcing upstream aggradation of water and sediment, and generating catastrophic outburst floods. Less apparent is the effect of large landslide dams on river ecosystems and marine sedimentation. Combining analyses of 1-m resolution topographic data (acquired via airborne laser mapping) and field investigation, we present evidence for a large, landslide-dammed paleolake along the Eel River, CA. The landslide mass initiated from a high-relief, resistant outcrop which failed catastrophically, blocking the Eel River with an approximately 130-m-tall dam. Support for the resulting 55-km-long, 1.3-km(3) lake includes subtle shorelines cut into bounding terrain, deltas, and lacustrine sediments radiocarbon dated to 22.5 ka. The landslide provides an explanation for the recent genetic divergence of local anadromous (ocean-run) steelhead trout (Oncorhynchus mykiss) by blocking their migration route and causing gene flow between summer run and winter run reproductive ecotypes. Further, the dam arrested the prodigious flux of sediment down the Eel River; this cessation is recorded in marine sedimentary deposits as a 10-fold reduction in deposition rates of Eel-derived sediment and constitutes a rare example of a terrestrial event transmitted through the dispersal system and recorded offshore.