skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 11:00 PM ET on Thursday, May 23 until 2:00 AM ET on Friday, May 24 due to maintenance. We apologize for the inconvenience.

Title: Controls on the size distributions of shallow landslides

Rainfall-triggered shallow landslides are destructive hazards and play an important role in landscape processes. A theory explaining the size distributions of such features remains elusive. Prior work connects size distributions to topography, but field-mapped inventories reveal pronounced similarities in the form, mode, and spread of distributions from diverse landscapes. We analyze nearly identical distributions occurring in the Oregon Coast Range and the English Lake District, two regions of strikingly different topography, lithology, and vegetation. Similarity in minimum sizes at these sites is partly explained by theory that accounts for the interplay of mechanical soil strength controls resisting failure. Maximum sizes, however, are not explained by current theory. We develop a generalized framework to account for the entire size distribution by unifying a mechanistic slope stability model with a flexible spatial-statistical description for the variability of hillslope strength. Using hillslope-scale numerical experiments, we find that landslides can occur not only in individual low strength areas but also across multiple smaller patches that coalesce. We show that reproducing observed size distributions requires spatial strength variations to be strongly localized, of large amplitude, and a consequence of multiple interacting factors. Such constraints can act together with the mechanical determinants of landslide initiation to produce size distributions of broadly similar character in widely different landscapes, as found in our examples. We propose that size distributions reflect the systematic scale dependence of the spatially averaged strength. Our results highlight the critical need to constrain the form, amplitude, and wavelength of spatial variability in material strength properties of hillslopes.

more » « less
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range / eLocation ID:
Article No. e2021855118
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Shallow bedrock strength controls both landslide hazard and the rate and form of erosion, yet regional patterns in near‐surface mechanical properties are rarely known quantitatively due to the challenge in collectingin situmeasurements. Here we present seismic and geomechanical characterizations of the shallow subsurface across the central Himalayan Range in Nepal. By pairing widely distributed 1D shear wave velocity surveys and engineering outcrop descriptions per the Geological Strength Index classification system, we evaluate landscape‐scale patterns in near‐surface mechanical characteristics and their relation to environmental factors known to affect rock strength. We find that shallow bedrock strength is more dependent on the degree of chemical and physical weathering, rather than the mineral and textural differences between the metamorphic lithologies found in the central Himalaya. Furthermore, weathering varies systematically with topography. Bedrock ridge top sites are highly weathered and have S‐wave seismic velocities and shear strength characteristics that are more typical of soils, whereas sites near valley bottoms tend to be less weathered and characterized by high S‐wave velocities and shear strength estimates typical of rock. Weathering on hillslopes is significantly more variable, resulting in S‐wave velocities that range between the ridge and channel endmembers. We speculate that variability in the hillslope environment may be partly explained by the episodic nature of mass wasting, which clears away weathered material where landslide scars are recent. These results underscore the mechanical heterogeneity in the shallow subsurface and highlight the need to account for variable bedrock weathering when estimating strength parameters for regional landslide hazard analysis.

    more » « less
  2. Abstract Distributions of landslide size are hypothesized to reflect hillslope strength, and consequently weathering patterns. However, the association of weathering and critical zone architecture with mechanical strength properties of parent rock and soil are poorly-constrained. Here we use three-dimensional stability to analyze 7330 landslides in western Oregon to infer combinations of strength - friction angles and cohesion - through analysis of both failed and reconstructed landslide terrain. Under a range of conditions, our results demonstrate that the failure envelope that relates shear strength and normal stress in landslide terrain is nonlinear owing to an exchange in strength with landslide thickness. Despite the variability in material strength at large scales, the observed gradient in proportional cohesive strength with landslide thickness may serve as a proxy for subsurface weathering. We posit that the observed relationships between strength and landslide thickness are associated with the coalescence of zones of low shear strength driven by fractures and weathering, which constitutes a first-order control on the mechanical behavior of underlying soil and rock mass. 
    more » « less
  3. Abstract

    The potentially important influence of climate change on landscape evolution and on critical zone processes is not sufficiently understood. The relative contribution of hydro‐climatic factors on hillslope erosion rates may significantly vary with topography at the watershed scale. The objective of this study is to quantify the hydro‐geomorphic behavior of two contrasting landscapes in response to different climate change scenarios in the Luquillo Critical Zone Observatory, a site of particular geomorphological interest, in terms of hillslope erosion and rainfall‐triggered landslides. We investigate the extent to which hillslope erosion and landslide occurrence remain relatively invariant with future hydro‐climatic perturbations. The adjacent Mameyes and Icacos watersheds are studied, which are underlain by contrasting lithologies. A high resolution coupled hydro‐geomorphic model based on tRIBS (Triangulated Irregular Network‐based Real‐time Integrated Basin Simulator) is used. Observations of landslide activity and hillslope erosion are used to evaluate the model performance. The process‐based model quantifies feedbacks among different hydrologic processes, landslide occurrence, and topsoil erosion and deposition. Simulations suggest that the propensity for landslide occurrence in the Luquillo Mountains is controlled by tropical storms, subsurface water flow, and by non‐climatic factors, and is expected to remain significant through 2099. The Icacos watershed, which is underlain by quartz diorite, is dominated by relatively large landslides. The relative frequency of smaller landslides is higher at the Mameyes watershed, which is underlain by volcaniclastic rock. While projections of precipitation decrease at the study site may lead to moderate decline in hillslope erosion rates, the simulated erosional potential of the two diverse landscapes likely remains significant. © 2018 John Wiley & Sons, Ltd.

    more » « less
  4. Abstract

    The size distributions of sediment delivered from hillslopes to rivers profoundly influence river morphodynamics, including river incision into bedrock and the quality of aquatic habitat. Yet little is known about the factors that influence size distributions of sediment produced by weathering on hillslopes. We present results of a field study of hillslope sediment size distributions at Inyo Creek, a steep catchment in granitic bedrock of the Sierra Nevada, USA. Particles sampled near the base of hillslopes, adjacent to the trunk stream, show a pronounced decrease in sediment size with decreasing sample elevation across all but the coarsest size classes. Measured size distributions become increasingly bimodal with decreasing elevation, exhibiting a coarse, bouldery mode that does not change with elevation and a more abundant finer mode that shifts from cobbles at the highest elevations to gravel at mid elevations and finally to sand at low elevations. We interpret these altitudinal variations in hillslope sediment size to reflect changes in physical, chemical, and biological weathering that can be explained by the catchment's strong altitudinal gradients in topography, climate, and vegetation cover. Because elevation and travel distance to the outlet are closely coupled, the altitudinal trends in sediment size produce a systematic decrease in sediment size along hillslopes parallel to the trunk stream. We refer to this phenomenon as ‘downvalley fining.’ Forward modeling shows that downvalley fining of hillslope sediment is necessary for downstream fining of the long‐term average flux of coarse sediment in mountain landscapes where hillslopes and channels are coupled and long‐term net sediment deposition is negligible. The model also shows that abrasion plays a secondary role in downstream fining of coarse sediment flux but plays a dominant role in partitioning between the bedload and suspended load. Patterns observed at Inyo Creek may be widespread in mountain ranges around the world. © 2020 John Wiley & Sons Ltd

    more » « less
  5. Abstract. Carbonate rocks are highly reactive and can have higher ratios of chemical weathering to total denudation relative to most other rock types. Their chemical reactivity affects the first-order morphology of carbonate-dominated landscapes and their climate sensitivity to weathering.However, there have been few efforts to quantify the partitioning ofdenudation into mechanical erosion and chemical weathering in carbonatelandscapes such that their sensitivity to changing climatic and tectonicconditions remains elusive. Here, we compile bedrock and catchment-averagedcosmogenic calcite–36Cl denudation rates and compare them to weathering rates derived from stream water chemistry from the same regions. Local bedrock denudation and weathering rates are comparable, ∼20–40 mm ka−1, whereas catchment-averaged denudation rates are ∼2.7 times higher. The discrepancy between bedrock and catchment-averaged denudation is 5 times lower compared to silicate-rich rocks, illustrating that elevated weathering rates make denudation more spatially uniform in carbonate-dominated landscapes. Catchment-averaged denudation rates correlate well with topographic relief and hillslope gradients, and moderate correlations with runoff can be explained by concurrent increases in weathering rates. Comparing denudation rates with weathering rates shows that mechanical erosion processes contribute ∼50 % of denudation in southern France and ∼70 % in Greece and Israel. Our results indicate that the partitioning between largely slope-independent chemical weathering and slope-dependent mechanical erosion varies based on climate and tectonics and impacts the landscape morphology. This leads us to propose a conceptual model whereby in humid, slowly uplifting regions, carbonates are associated with low-lying, flat topography because slope-independent chemical weathering dominates denudation. In contrast, in arid climates with rapid rock uplift rates, carbonate rocks form steep mountains that facilitate rapid, slope-dependent mechanical erosion required to compensate for inefficient chemical weathering and runoff loss to groundwater systems. This result suggests that carbonates represent an end member for interactions between climate, tectonics, and lithology. 
    more » « less