skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: The Art of Landslides: How Stochastic Mass Wasting Shapes Topography and Influences Landscape Dynamics
Abstract 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.  more » « less
Award ID(s):
2148762 2104102
PAR ID:
10370327
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
127
Issue:
8
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics. 
    more » « less
  2. Rapid sediment accumulation rates (SAR) in a fan delta situated on the rapidly uplifting footwall of the Taormina normal fault in NE Sicily preserves a rare record of earthquakes and base level change for a tightly coupled source to sink system. We use this sedimentary archive to reconstruct the kinematics and slip history of the fault and further an understanding of how tectonic forcing across various scales are encoded in stratigraphy. A revised luminescence-based age model indicates that ~82 m of the Pagliara fan-delta foreset facies was deposited in ~11 ka at a mean SAR of ~0.74 cm/yr during MIS 7. Syn-depositional terrestrial cosmogenic nuclide (TCN) determined paleoerosion rates of 0.91±0.12 mm/yr and 1.31 ±0.61 mm/yr are similar to published modern erosion rates for the Pagliara basin of 0.97 ±0.11 mm/yr. At the stratigraphic scale, a time series of magnetic susceptibility (c) sampled at 1 m intervals in the foresets displays four ~2,800 yr / 20 m-thick cycles of growing c, bounded by sharp decreases that do not coincide with changes in sediment texture. The c of the low-grade metamorphic bedrock in the source is 20-100 times weaker than the c of rubified soils mantling the hillslopes, which is comparable to the c of the delta sediments. We propose that large, bedrock-cored landslides quasi-periodically deliver weak c sediment to the delta that dilutes a c signal otherwise dominated by the stripping of soil-mantled hillslopes. We propose that centennial-scale recurrence interval earthquakes are most capable at triggering a basin-scale landslide only after channel incision has increased relief of hillslopes to the threshold condition, which requires millennia to achieve. At the landscape scale of delta geometry and location, the Pagliara delta accumulated in a hanging wall basin that has since been inverted. We reconstruct the history of base level fall for the delta from an inversion of fluvial topography and apportion that record to its rock uplift, delta deposition, and eustatic components. We show that footwall uplift has been unsteady over the past 600 ka ranging from -1 to 3 mm/yr. The integration of our stratigraphic- and landscape scale observations furthers our understanding of the natural hazards related to normal fault earthquakes and their impact on sediment dynamics in this steep, active tectonic setting. 
    more » « less
  3. Abstract Landscape form represents the cumulative effects of de‐stabilizing events relative to recovery processes. Most geomorphic research has focused on the role of episodic rare events on landscape form with less attention paid to the role and persistence of chronic inputs. To better establish the interplay between chronic and episodic extreme events at regional scales, we used aerial photography and post‐flood sediment sampling to assess stream and hillslope response and recovery to a 100–300 yr. flood caused by Tropical Storm Irene in New England. Within a 14 000 km2study area, analysis of aerial photographs indicated that the storm initiated (n = 534) and reactivated (n = 460) a large number of landslides. These landslides dramatically increased overall estimates of regional erosion rates (from 0.0023 mm/yr. without Irene to 0.0072 mm/yr. with Irene). Similarly, Irene‐generated LWD inputs of 0.25–0.5 trees/km exceeded annual background rates in a single event, and these concentrated inputs (101–102of trees/landslide) are likely to result in large jams and snags that are particularly persistent and geomorphically effective. Finally, we found that landslide scars continue to provide elevated sediment inputs years after the event, as evidenced by sustained higher suspended sediment concentrations in streams with Irene‐generated landslides. Overall, our results indicate that infrequent, high‐magnitude events have a more important geomorphic role in tectonically stable, more moderate‐relief systems than has been previously recognized. Understanding the role of these events has particular relevance in regions such as New England, where the frequency and magnitude of extreme storms is expected to increase. Further, these effects may force reconsideration of conservation and restoration targets (for example in channel form and large wood loading and distribution) in fluvial systems. Copyright © 2016 John Wiley & Sons, Ltd. 
    more » « less
  4. Abstract Regional‐scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics‐based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three‐dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad‐scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km2) using 3D, limit‐equilibrium (LE)‐based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back‐calculated from inventoried landslides using 3D LE‐based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA. 
    more » « less
  5. Abstract Large earthquakes can construct mountainous topography by inducing rock uplift but also erode mountains by causing landslides. Observations following the 2008 Wenchuan earthquake show that landslide volumes in some cases match seismically induced uplift, raising questions about how the actions of individual earthquakes accumulate to build topography. Here we model the two‐dimensional surface displacement field generated over a full earthquake cycle accounting for coseismic deformation, postseismic relaxation, landslide erosion, and erosion‐induced isostatic compensation. We explore the related volume balance across different seismotectonic and topographic conditions and revisit the Wenchuan case in this context. The ratio (Ω) between landslide erosion and uplift is most sensitive to parameters determining landslide volumes (particularly earthquake magnitudeMw, seismic energy source depth, and failure susceptibility, as well as the seismological factor responsible for triggering landslides), and is moderately sensitive to the effective elastic thickness of lithosphere,Te. For a specified magnitude, more erosive events (higher Ω) tend to occur at shallower depth, in thicker‐Telithosphere, and in steeper, more landslide‐prone landscapes. For given landscape and seismotectonic conditions, the volumes of both landslides and uplift to first order positively scale withMwand seismic momentMo. However, higherMwearthquakes generate lower landslide and uplift volumes per unitMo, suggesting lower efficiency in the use of seismic energy to drive topographic change. With our model, we calculate the long‐term average seismic volume balance for the eastern Tibetan region and find that the net topographic effect of earthquakes in this region tends to be constructive rather than erosive. Overall, destructive events are rare when considering processes over the full earthquake cycle, although they are more likely if only considering the coseismic volume budget (as was the case for the 2008 Wenchuan earthquake where landsliding substantially offset coseismic uplift). Irrespective of the net budget, our results suggest that the erosive power of earthquakes plays an important role in mountain belt evolution, including by influencing structures and spatial patterns of deformation, for example affecting the wavelength of topography. 
    more » « less