skip to main content

Title: Unsaturated Flow Processes and the Onset of Seasonal Deformation in Slow‐Moving Landslides

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):
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
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. Abstract

    In water‐stressed regions of the world, the inundation of working landscapes to replenish aquifers—known as flood‐managed aquifer recharge (flood‐MAR)—has become a valuable tool for sustainable groundwater management. Due to their diverse land use histories, however, many potential recharge sites host nonpoint source contaminants (such as salts, pesticides, and fertilizers) within the vadose zone that may flush to groundwater during recharge operations. To identify the controls on contaminant migration, we perform stochastic simulations of flood‐MAR through a heterogeneous alluvial aquifer and apply transient particle tracking to evaluate conservative and reactive contaminant transport over 80 years of recharge operations. With semi‐annual recharge events, the water table begins to rise 0.13–1.84 years after the first inundation event while solutes take much longer (11 to 80 years) to transit the 45‐m thick unsaturated zone. We derive a parametric expression for the ratio of celerity (or rate of pressure transmission) to velocity of the flood‐MAR wetting front and show that this simplified expression agrees with values calculated from heterogeneous model simulations. Slow solute velocities (0.25–1.75 m year−1) allow for significant contaminant removal through denitrification, but the contaminant plume experiences minimal dispersion or dilution over this time, reaching the water table as a sharp front. Our results suggest that minimizing groundwater velocity and maximizing groundwater celerity during flood‐MAR should optimize increases in water supply while limiting water quality degradation.

    more » « less
  4. Abstract

    Predicting fluid biogeochemistry in the vadose zone is difficult because of time‐dependent variation in multiple controlling factors, such as temperature, moisture, and biological activity. Furthermore, soils are multicomponent, heterogeneous porous media where manifold reactions may be affecting solution chemistry. We postulated that ecosystem‐scale processes, such as carbon fixation and ecohydrologic partitioning, control subsurface biogeochemical reactions, including mineral weathering. To test this hypothesis, we applied a novel “instrumented pedon” research approach. Analysis of the data streams demonstrates the interactions between pulsed wetting events and biogeochemical processes in the soil profile, and along groundwater flow paths. Rapid wetting front propagation into dry soil resulted in a pulsed increase in CO2partial pressure in deeper soil layers, whereas wetting front propagation into a premoistened soil profile showed the opposite effect. The apparent respiratory quotient (ARQ), calculated from CO2and O2fluxes, deviated from expected oxidative ratios particularly during soil wetting events. These deviations were correlated in time with pore water geochemical responses, revealing that a fraction of the respired CO2was consumed locally in pulsed silicate weathering events that accompanied wetting‐front propagation. However, most of this CO2was dissolved in the soil pore water and transported downgradient, and along the soil‐bedrock interface, where a portion of it was further consumed in silicate weathering reactions, and another portion was degassed to the atmosphere. These results highlight the tight coupling that exists between physical, biological, and chemical processes, on event time scales, during incremental co‐evolution of the critical zone, particularly in water‐limited systems.

    more » « less
  5. Subsurface drainage has been widely accepted to mitigate the hazard of landslides in areas prone to flooding. Specifically, the use of drainage wells with pumping systems has been recognized as an effective short-term solution to lower the groundwater table. However, this method has not been well considered for long-term purposes due to potentially high labor costs. This study aims to investigate the idea of an autonomous pumping system for subsurface drainage by leveraging conventional geotechnical engineering solutions and a deep learning technique—Long-Short Term Memory (LSTM)—to establish a geotechnical cyber-physical system for rainfall-induced landslide prevention. For this purpose, a typical soil slope equipped with three pumps was considered in a computer simulation. Forty-eight cases of rainfall events with a wide range of varieties in duration, total rainfall depths, and different rainfall patterns were generated. For each rainfall event, transient seepage analysis was performed using newly proposed Python code to obtain the corresponding pump’s flow rate data. A policy of water pumping for maintaining groundwater at a desired level was assigned to the pumps to generate the data. The LSTM takes rainfall event data as the input and predicts the required pump’s flow rate. The results from the trained model were validated using evaluation metrics of root mean square error (RMSE), mean absolute error (MAE), and R2. The R2-scores of 0.958, 0.962, and 0.954 for the predicted flow rates of the three pumps exhibited high accuracy of the predictions using the trained LSTM model. This study is intended to make a pioneering step toward reaching an autonomous pumping system and lowering the operational costs in controlling geosystems. 
    more » « less