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: Unsaturated Flow Processes and the Onset of Seasonal Deformation in Slow‐Moving Landslides
Abstract Predicting rainfall‐induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow‐moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity‐duration threshold. Our model shows that the seasonal response time of slow‐moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.  more » « less
Award ID(s):
1658800
PAR ID:
10419422
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
126
Issue:
5
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The rate and timing of hydrologically forced landslides is a complex function of precipitation patterns, material properties, topography, and groundwater hydrology. In the simplest form, however, slopes fail when subsurface pore pressure grows large enough to exceed the Mohr‐Coulomb failure criterion. The capacity for pore pressure rise in a landslide is determined in part by the thickness of the unsaturated zone above the water table, which itself is set by weathering patterns that should have predictable patterns across different lithologies. To investigate how this structure affects landslide behavior, we exploit a multi‐year record of precipitation, pore pressure, and velocity from Oak Ridge earthflow, a slow‐moving landslide set in Franciscan mélange, northern California, USA. In conjunction with electrical resistivity tomography and hydraulic conductivity measurements, these data show that Oak Ridge has a thin weathered profile that is comparable in thickness to other mélange landslides in California. We propose that due to the inherently thin vadose zone, mélange landscapes experience an unusually high water table that frequently brings them close to movement; however, the capacity to increase stress is limited by the small amount of dynamic storage available. Instead, excess pore pressure is shed via springs and saturation overland flow once the water table reaches the surface. Linkages between weathering patterns, hydrology, and deformation can explain behavior patterns exhibited by Franciscan mélange earthflows across a large precipitation gradient. 
    more » « less
  2. Abstract Landslide motion is often simulated with interface‐like laws able to capture changes in frictional strength caused by the growth of the pore water pressure and the consequent reduction of the effective stress normal to the plane of sliding. Here it is argued that, although often neglected, the evolution of all the 3D stress components within the basal shear zone of landslides also contributes to changes in frictional strength and must be accounted for to predict changes in seasonal velocity. For this purpose, an augmented sliding‐consolidation model is proposed which allows for the computation of excess pore pressure development and downslope sliding with any constitutive law with 3D stress evolution. Simulations of idealised infinite slope models subjected to hydrologic forcing are used to study the role of in‐situ stress conditions and stress rate multiaxiality. Specifically, a Drucker‐Prager perfectly plastic model is used to replicate frictional failure and shear deformation at the base of landslides. The model reveals that conditions amenable to the shearing of a frictional interface are met only after numerous rainfall cycles, that is, when multiaxial stress rates are suppressed. In this case, the landslide is predicted to move through a seasonal ratcheting controlled only by the effective stress component normal to the plane of sliding. By contrast, in newly formed landslides, the multiaxial stress evolution is found to produce further regimes of motion, from plastic shakedown to cyclic failure, neither of which can be captured by interface‐like frictional laws. Notably, the model suggests that a transition across these regimes can emerge in response to an aggravation of the magnitude of forcing, implying that (i) fluctuations in climate may alter the seasonal trends of motion observed today; (ii) our ability to quantify landslide‐induced risks is impaired unless proper geomechanical models are used to examine their long‐term dynamics. 
    more » « less
  3. Whether Earth materials exhibit frictional creep or catastrophic failure is a crucial but unresolved problem in predicting landslide and earthquake hazards. Here, we show that field-scale observations of sliding velocity and pore water pressure at two creeping landslides are explained by velocity-strengthening friction, in close agreement with laboratory measurements on similar materials. This suggests that the rate-strengthening friction commonly measured in clay-rich materials may govern episodic slow slip in landslides, in addition to tectonic faults. Further, our results show more generally that transient slow slip can arise in velocity-strengthening materials from modulation of effective normal stress through pore pressure fluctuations. This challenges the idea that episodic slow slip requires a narrow range of transitional frictional properties near the stability threshold, or pore pressure feedbacks operating on initially unstable frictional slip. 
    more » « less
  4. Abstract Empirical thresholds for landslide warning systems have benefitted from the incorporation of soil‐hydrologic monitoring data, but the mechanistic basis for their predictive capabilities is limited. Although physically based hydrologic models can accurately simulate changes in soil moisture and pore pressure that promote landslides, their utility is restricted by high computational costs and nonunique parameterization issues. We construct a deep learning model using soil moisture, pore pressure, and rainfall monitoring data acquired from landslide‐prone hillslopes in Oregon, USA, to predict the timing and magnitude of hydrologic response at multiple soil depths for 36‐hr intervals. We find that observation records as short as 6 months are sufficient for accurate predictions, and our model captures hydrologic response for high‐intensity rainfall events even when those storm types are excluded from model training. We conclude that machine learning can provide an accurate and computationally efficient alternative to empirical methods or physical modeling for landslide hazard warning. 
    more » « less
  5. Abstract Catastrophic landslides are often preceded by slow, progressive, accelerating deformation that differs from the persistent motion of slow‐moving landslides. Here, we investigate the motion of a landslide that damaged 12 homes in Rolling Hills Estates (RHE), Los Angeles, California on 8 July 2023, using satellite‐based synthetic aperture radar interferometry (InSAR) and pixel tracking of satellite‐based optical images. To better understand the precursory motion of the RHE landslide, we compared its behavior with local precipitation and with several slow‐moving landslides nearby. Unlike the slow‐moving landslides, we found that RHE was a first‐time progressive failure that failed after one of the wettest years on record. We then applied a progressive failure model to interpret the failure mechanisms and further predict the failure time from the pre‐failure movement of RHE. Our work highlights the importance of monitoring incipient slow motion of landslides, particularly where no discernible historical displacement has been observed. 
    more » « less