Landslides, threshold slopes and the survival of relict terrain in the wake of the Mendocino Triple Junction

Establishing landscape response to uplift is critical for interpreting sediment fluxes, hazard potential, and topographic evolution. We assess how landslides shape terrain in response to a wave of uplift traversing the northern California Coast Ranges (United States) in the wake of the Mendocino Triple Junction. We extracted knickpoints, landslide erosion rates, and topographic metrics across the region modified by Mendocino Triple Junction migration. Landslide erosion rates mapped from aerial imagery are consistent with modeled uplift and exhumation, while hillslope gradient is invariant across the region, suggesting that landslides accommodate uplift, as predicted by the threshold slope model. Landslides are concentrated along steepened channel reaches downstream of knickpoints generated by base-level fall at channel outlets, and limit slope angles and relief. We find evidence that landslide-derived coarse sediment delivery may suppress catchment-wide channel incision and landscape denudation over the time required for the uplift wave to traverse the region. We conclude that a landslide cover effect may provide a mechanism for the survival of relict terrain and orogenic relief in the northern Californian Coast Ranges and elsewhere over millennial time scales

Beyond the angle of repose: A review and synthesis of landslide processes in response to rapid uplift, Eel River, Northern California

In mountainous settings, increases in rock uplift are often followed by a commensurate uptick in denudation as rivers incise and steepen hillslopes, making them increasingly prone to landsliding as slope angles approach a limiting value. For decades, the threshold slope model has been invoked to account for landslide-driven increases in sediment flux that limit topographic relief, but the manner by which slope failures organize themselves spatially and temporally in order for erosion to keep pace with rock uplift has not been well documented. Here, we review past work and present new findings from remote sensing, cosmogenic radionuclides, suspended sediment records, and airborne lidar data, to decipher patterns of landslide activity and geomorphic processes related to rapid uplift along the northward-migrating Mendocino Triple Junction in Northern California. From historical air photos and airborne lidar, we estimated the velocity and sediment flux associated with active, slow-moving landslides (or earthflows) in the mélange- and argillite-dominated Eel River watershed using the downslope displacement of surface markers such as trees and shrubs. Although active landslides that directly convey sediment into the channel network account for only 7% of the landscape surface, their sediment flux amounts to more than 50% of the suspended load recorded at downstream sediment gaging stations. These active slides tend to exhibit seasonal variations in velocity as satellite-based interferometry has demonstrated that rapid acceleration commences within 1 to 2 months of the onset of autumn rainfall events before slower deceleration ensues in the spring and summer months. Curiously, this seasonal velocity pattern does not appear to vary with landslide size, suggesting that complex hydrologic–mechanical feedbacks (rather than 1-D pore pressure diffusion) may govern slide dynamics. A new analysis of 14 yrs of discharge and sediment concentration data for the Eel River indicates that the characteristic mid-winter timing of earthflow acceleration corresponds with increased suspended concentration values, suggesting that the seasonal onset of landslide motion each year may be reflected in the export of sediments to the continental margin. The vast majority of active slides exhibit gullied surfaces and the gully networks, which are also seasonally active, may facilitate sediment export although the proportion of material produced by this pathway is poorly known. Along Kekawaka Creek, a prominent tributary to the Eel River, new analyses of catchment-averaged erosion rates derived from cosmogenic radionuclides reveal rapid erosion (0.76 mm/yr) below a prominent knickpoint and slower erosion (0.29 mm/yr) upstream. Such knickpoints are frequently observed in Eel tributaries and are usually comprised of massive (> 10 m) interlocking resistant boulders that likely persist in the landscape for long periods of time (> 105 yr). Upstream of these knickpoints, active landslides tend to be less frequent and average slope angles are slightly gentler than in downstream areas, which indicates that landslide density and average slope angle appear to increase with erosion rate. Lastly, we synthesize evidence for the role of large, catastrophic landslides in regulating sediment flux and landscape form. The emergence of resistant blocks within the mélange bedrock has promoted large catastrophic slides that have dammed the Eel River and perhaps generated outburst events in the past. The frequency and impact of these landslide dams likely depend on the spatial and size distributions of resistant blocks relative to the width and drainage area of adjacent valley networks. Overall, our findings demonstrate that landslides within the Eel River catchment do not occur randomly, but instead exhibit spatial and temporal patterns related to baselevel lowering, climate forcing, and lithologic variations. Combined with recent landscape evolution models that incorporate landslides, these results provide predictive capability for estimating erosion rates and managing hazards in mountainous regions

Long-term patterns of hillslope erosion by earthquake-induced landslides shape mountain landscapes

Widespread triggering of landslides by large storms or earthquakes is a dominant mechanism of erosion in mountain landscapes. If landslides occur repeatedly in particular locations within a mountain range, then they will dominate the landscape evolution of that section and could leave a fingerprint in the topography. Here, we track erosion provenance using a novel combination of the isotopic and molecular composition of organic matter deposited in Lake Paringa, New Zealand. We find that the erosion provenance has shifted markedly after four large earthquakes over 1000 years. Postseismic periods eroded organic matter from a median elevation of 722 +329/−293 m and supplied 43% of the sediment in the core, while interseismic periods sourced from lower elevations (459 +256/−226 m). These results are the first demonstration that repeated large earthquakes can consistently focus erosion at high elevations, while interseismic periods appear less effective at modifying the highest parts of the topography

The Art of Landslides: How Stochastic Mass Wasting Shapes Topography and Influences Landscape Dynamics

Bedrock landslides shape topography and mobilize large volumes of sediment. Yet, interactions between landslide-produced sediment and fluvial systems that together govern large-scale landscape evolution are not well understood. To explain morphological patterns observed in steep, landslide-prone terrain, we explicitly model stochastic landsliding and associated sediment dynamics. The model accounts for several common landscape features such as slope frequency distributions, which include values in excess of regional stability limits, quasi-planar hillslopes decorated with straight, closely spaced channel-like features, and accumulation of sediment in valley networks rather than on hillslopes. Stochastic landsliding strongly affects the magnitude and timing of sediment supply to the fluvial system. We show that intermittent sediment supply is ultimately reflected in topography. At dynamic equilibrium, landslide-derived sediment pulses generate persistent landscape dynamism through the formation and breaching of landslide dams and epigenetic gorges as landslides force shifts in channel positions. Our work highlights the importance of interactions between landslides and sediment dynamics that ultimately control landscape-scale response to environmental change

Geomorphological and hydrological controls on sediment export in earthquake-affected catchments in the Nepal Himalaya

Large earthquakes can contribute to mountain growth by building topography but also contribute to mass removal from mountain ranges through widespread mass wasting. On annual to decadal or centennial timescales, large earthquakes also have the potential to significantly alter fluvial sediment dynamics if a significant volume of the sediment generated reaches the fluvial network. In this contribution, we focus on the Melamchi–Indrawati and Bhote Koshi rivers in central Nepal, which have both experienced widespread landsliding associated with the 2015 Gorkha (Nepal) earthquake. Using a time series of high-resolution satellite imagery, we have mapped exposed sediment along the rivers from 2012–2021 to identify zones of active channel deposition and document changes over time. Counter to expectations, we show negligible increases in coarse-sediment accumulation along both river corridors since the Gorkha earthquake. However, an extremely high-concentration flow event on 15 June 2021 caused an approximately 4-fold increase in exposed sediment along a 30 km reach of the channel with up to 12 m of channel aggradation in the Melamchi–Indrawati rivers; this event was localised and did not impact the neighbouring Bhote Koshi catchment. Based on published reports, new helicopter-based photography, and satellite data, we demonstrate that this event was sourced from a localised rainfall event between 4500 and 4800 m and that a significant fraction of the sediment was supplied from sources that were unrelated to the landslides generated by the Gorkha earthquake.</p

The fate of sediment after a large earthquake

Large earthquakes rapidly denude hillslopes by triggering thousands of coseismic landslides. The sediment produced by these landslides is initially quickly mobilised from the landscape by an interconnected cascade of processes. This cascade can dramatically but briefly enhance local erosion rates. Hillslope and channel processes, such as landsliding and debris flows, interact to influence the total mass, caliber, and rate of sediment transport through catchments. Calculating the sediment budget of an earthquake lends insight into the nature of these interactions. Using satellite imagery derived landslide inventories, channel surveys and a literature review combined with a Monte Carlo simulation approach we present a constrained sediment budget of the first decade after the 2008 Mw7.9 Wenchuan earthquake. With this sediment budget we demonstrate that debris flows are dominant process for delivering sediment into channels and that large volumes of sediment remain in the landscape. In our study area over 88% (469.7 Mega tonnes) of the coseismically generated sediment remains on the hillslopes in 2018. Of the 12% of the sediment that was mobilised, 67% (45.2 ± 22 Mt) was mobilised by debris flows. Despite the large proportion of sediment remaining on the hillslope, the frequency of debris flows declined significantly over our observation period. The reduction in debris-flow frequency is not correlated to reductions in the frequency of triggering storms, suggesting changes in the mechanical properties of hillslope sediment may drive this observation. The stabilization of coseismically generated sediment greatly extends its residence time and may influence catchment sediment yields for centuries or millennia

The role of infrequently mobile boulders in modulating landscape evolution and geomorphic hazards

A landscape’s sediment grain size distribution is the product of, and an important influence on, earth surface processes and landscape evolution. Grains can be large enough that the motion of a single grain, infrequently mobile in size-selective transport systems, constitutes or triggers significant geomorphic change. We define these grains as boulders. Boulders affect landscape evolution; their dynamics and effects on landscape form have been the focus of substantial recent community effort. We review progress on five key questions related to how boulders influence the evolution of unglaciated, eroding landscapes: 1) What factors control boulder production on eroding hillslopes and the subsequent downslope evolution of the boulder size distribution? 2) How do boulders influence hillslope processes and long-term hillslope evolution? 3) How do boulders influence fluvial processes and river channel shape? 4) How do boulder-mantled channels and hillslopes interact to set the long-term form and evolution of boulder-influenced landscapes? 5) How do boulders contribute to geomorphic hazards, and how might improved understanding of boulder dynamics be used for geohazard mitigation? Boulders are produced on eroding hillslopes by landsliding, rockfall, and/or exhumation through the critical zone. On hillslopes dominated by local sediment transport, boulders affect hillslope soil production and transport processes such that the downslope boulder size distribution sets the form of steady-state hillslopes. Hillslopes dominated by nonlocal sediment transport are less likely to exhibit boulder controls on hillslope morphology as boulders are rapidly transported to the hillslope toe. Downslope transport delivers boulders to eroding rivers where the boulders act as large roughness elements that change flow hydraulics and the efficiency of erosion and sediment transport. Over longer timescales, river channels adjust their geometry to accommodate the boulders supplied from adjacent hillslopes such that rivers can erode at the baselevel fall rate given their boulder size distribution. The delivery of boulders from hillslopes to channels, paired with the channel response to boulder delivery, drives channel-hillslope feedbacks that affect the transient evolution and steady-state form of boulder-influenced landscapes. At the event scale, boulder dynamics in eroding landscapes represent a component of geomorphic hazards that can be mitigated with an improved understanding of the rates and processes associated with boulder production and mobility. Opportunities for future work primarily entail field-focused data collection across gradients in landscape boundary conditions (tectonics, climate, and lithology) with the goal of understanding boulder dynamics as one component of landscape self-organization

Boulders as a Lithologic Control on River and Landscape Response to Tectonic Forcing at the Mendocino Triple Junction

Constraining Earth’s sediment mass balance over geologic time requires a quantitative understanding of how landscapes respond to transient tectonic perturbations. However, the mechanisms by which bedrock lithology governs landscape response remain poorly understood. Rock type influences the size of sediment delivered to river channels, which controls how efficiently rivers respond to tectonic forcing. The Mendocino triple junction region of northern California, USA, is one landscape in which large boulders, delivered by hillslope failures to channels, may alter the pace of landscape response to a pulse of rock uplift. Boulders frequently delivered by earthflows in one lithology, the Franciscan mélange, have been hypothesized to steepen channels and slow river response to rock uplift, helping to preserve high-elevation, low-relief topography. Channels in other units (the Coastal Belt and the Franciscan schist) may experience little or no erosion inhibition due to boulder delivery. Here we investigate spatial patterns in channel steepness, an indicator of erosion resistance, and how it varies between mélange and non-mélange channels. We then ask whether lithologically controlled boulder delivery to rivers is a possible cause of steepness variations. We find that mélange channels are steeper than Coastal Belt channels but not steeper than schist channels. Though channels in all units steepen with increasing proximity to mapped hillslope failures, absolute steepness values near failures are much higher (∼2×) in the mélange and schist than in Coastal Belt units. This could reflect reduced rock erodibility or increased erosion rates in the mélange and schist, or disproportionate steepening due to enhanced boulder delivery by hillslope failures in those units. To investigate the possible influence of lithology-dependent boulder delivery, we map boulders at failure toes in the three units. We find that boulder size, frequency, and concentration are greatest in mélange channels and that Coastal Belt channels have the lowest concentrations. Using our field data to parameterize a mathematical model for channel slope response to boulder delivery, we find that the modeled influence of boulders in the mélange could be strong enough to account for some observed differences in channel steepness between lithologies. At the landscape scale, we lack the data to fully disentangle boulder-induced steepening from that due to spatially varying erosion rates and in situ rock erodibility. However, our boulder mapping and modeling results suggest that lithology-dependent boulder delivery to channels could retard landscape adjustment to tectonic forcing in the mélange and potentially also in the schist. Boulder delivery may modulate landscape response to tectonics and help preserve high-elevation, low-relief topography at the Mendocino triple junction and elsewhere

Experimental Study on Landslides of Loose Sediment Slope Induced by Stream Bed Incision

Uplift of the Qinghai-Tibetan Plateau has resulted in rapid incision of rivers along the margin of the plateau. Landslides occur frequently as a consequence of increasing bank slope and potential landslide energy due to stream bed incision or lateral bank erosion on the concave banks at bends. The Fencha Gully is on the eastern margin of the Qinghai-Tibetan Plateau and is developing on a huge landslide body. Flume experiments were conducted on the base of the field investigation to study the mechanism of landslides induced by stream bed incision. The experiments were designed with a length scale ratio of 1:20. Landslides and stream bed incision with loose sediment were observed and analyzed. The results show that landslides are induced as a result of stream bed incision. The potential landslide energy is defined, which increases quickly with an effective incision depth coupling vertical incision and lateral bank erosion. The occurrence of landslides can be attributed to increasing incision depth and potential landslide energy. Results indicate that the critical effective incision depth is 4.0–6.0 m. A critical value of the potential landslide energy is found from the experiments. Landslides occur if the potential energy exceeds the critical energy, which is 2.24×104 t·m/s2 for the Fencha Gully. The incision depth and potential energy of landslides from the Fencha Gully agree well with the results

The fate of landslide debris after large earthquakes, Wenchuan, China

Large earthquakes are significant geomorphic events in mountain ranges;they can cause uplift while also triggering huge volumes of bedrock landsliding. Constraining theerosion of coseismic landslide deposits is important forunderstandingthe long-term impact of earthquakes on landscape evolution and natural hazards. The 2008 Mw7.9 Wenchuan earthquake provides an excellent opportunity to constrain the processes and rates that erode landslide deposits within and out of mountain catchments. I constructed a large multitemporal landslide inventory of the epicentral regionof the earthquaketo map the erosion of landslide deposits in the 10 years after the earthquake. In these 10 years less than 13% of the sediment hadbeen deposited into the stores associated with the main trunk of the orogen draining river leaving up to 88% stored on the hillslope. Debris flows were the key remobilising process, transporting over 50% of all the sediment that is mobilised. Using slope stability modelling and analysis of the landslide inventory I identified most of the debris flows were triggered on the hill slope and their frequency decreased rapidly through time. Initially debris flows were triggered in small deposits with low drainage areasbut through time the landslides closer to channel network become more likely to be the source of debris flow activity, a pattern that highlights a change in the physical properties of the deposits. I evaluate several theories that may explain this pattern and highlight strategies to identify the stabilising processes. Finally, I present a 0-Dimensional model investigatingthe impact of storage of coseismic sediment on landscape evolution. This model suggests coseismic sediment can remain in the mountain range for timescales greater than the recurrence time of large earthquakes causing a decrease in exhumation ratesas sediment is remobilised. The long storage time of sediment in orogens reduces the impact of earthquakes on long term (>1000 years) exhumation records so that single earthquakes are only recorded in the most local or highest resolution records