skip to main content

Title: Vadose Zone Thickness Limits Pore‐Fluid Pressure Rise in a Large, Slow‐Moving Earthflow

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
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

    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
  2. Abstract

    Explanations for distinct adjacent ecosystems that extend across hilly landscapes typically point to differences in climate or land use. Here we document—within a similar climate—how contrasting regional plant communities correlate with distinct underlying lithology and reveal how differences in water storage capacity in the critical zone (CZ) explain this relationship. We present observations of subsurface CZ structure and groundwater dynamics from deep boreholes and quantify catchment‐wide dynamic water storage in two Franciscan rock types of the Northern California Coast Ranges. Our field sites have a Mediterranean climate, where rains are out of phase with solar energy, amplifying the importance of subsurface water storage for periods of peak ecosystem productivity in the dry season. In the deeply weathered (~30 m at ridge) Coastal Belt argillite and sandstone, ample, seasonally replenished rock moisture supports an evergreen forest and groundwater drainage sustains baseflow throughout the summer. In the Central Belt argillite‐matrix mélange, a thin CZ (~3 m at ridge) limits total dynamic water storage capacity (100–200 mm) and rapidly sheds winter rainfall via shallow storm and saturation overland flow, resulting in low plant‐available water (inferred from predawn tree water potential) and negligible groundwater storage that can drain to streams in summer. This storage limitation mechanism explains the presence of an oak savanna‐woodland bounded by seasonally ephemeral streams, despite >1,800 mm of average precipitation. Through hydrologic monitoring and subsurface characterization, we reveal a mechanism by which differences in rock type result in distinct regionally extensive plant communities under a similar climate.

    more » « less
  3. Abstract. To explore the sensitivity of rivers to blocking from landslidedebris, we exploit two similar geomorphic settings in California'sFranciscan mélange where slow-moving landslides, often referred to asearthflows, impinge on river channels with drainage areas that differ by afactor of 30. Analysis of valley widths and river long profiles over∼19 km of Alameda Creek (185 km2 drainage area) andArroyo Hondo (200 km2 drainage area) in central California shows avery consistent picture in which earthflows that intersect these channelsforce tens of meters of gravel aggradation for kilometers upstream, leadingto apparently long-lived sediment storage and channel burial at these sites.In contrast, over a ∼30 km section of the Eel River (5547 km2 drainage area), there are no knickpoints or aggradation upstreamof locations where earthflows impinge on its channel. Hydraulic andhydrologic data from United States Geological Survey (USGS) gages on Arroyo Hondo and the Eel River, combinedwith measured size distributions of boulders input by landslides for bothlocations, suggest that landslide derived boulders are not mobile at eithersite during the largest floods (>2-year recurrence) with field-measured flow depths. We therefore argue that boulder transport capacity isan unlikely explanation for the observed difference in sensitivity tolandslide inputs. At the same time, we find that earthflow fluxes per unitchannel width are nearly identical for Oak Ridge earthflow on Arroyo Hondo,where evidence for blocking is clear, and for the Boulder Creek earthflow onthe Eel River, where evidence for blocking is absent. These observationssuggest that boulder supply is also an unlikely explanation for the observedmorphological differences along the two rivers. Instead, we argue that thedramatically different sensitivity of the two locations to landslideblocking is related to differences in channel width relative to typicalseasonal displacements of earthflows. A synthesis of seasonal earthflowdisplacements in the Franciscan mélange shows that the channel width ofthe Eel River is ∼5 times larger than the largest annualseasonal displacement. In contrast, during wet winters, earthflows arecapable of crossing the entire channel width of Arroyo Hondo and AlamedaCreek. In support of this interpretation, satellite imagery shows thatimmobile earthflow-derived boulders are generally confined to the edges ofthe channel on the Eel River. By contrast, immobile earthflow-derivedboulders jam the entire channel on Arroyo Hondo. Our results imply that lower drainage area reaches of earthflow-dominated catchments may be particularly prone to blocking. By inhibiting the upstreampropagation of base-level signals, valley-blocking earthflows may thereforepromote the formation of so-called “relict topography”. 
    more » « less
  4. 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
  5. Abstract

    Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment production within a landscape. Here, we investigate saprolite and weathered bedrock properties across a ridge‐valley system in the Northern California Coast Ranges, USA, where topography varies with slope aspect such that north‐facing slopes have thicker soils and are more densely vegetated than south‐facing slopes. We use active source seismic refraction surveys to extend observations made in boreholes to the hillslope scale. Seismic velocity models across several ridges capture a high velocity gradient zone (from 1,000 to 2,500 m/s) located ∼4–13 m below ridgetops that coincides with transitions in material strength and chemical depletion observed in boreholes. Comparing this transition depth across multiple north‐ and south‐facing slopes, we find that the thickness of saprolite does not vary with slope aspects. Additionally, seismic survey lines perpendicular and parallel to bedding planes reveal weathering profiles that thicken upslope and taper downslope to channels. Using a rock physics model incorporating seismic velocity, we estimate the total porosity of the saprolite and find that inherited fractures contribute a substantial amount of pore space in the upper 6 m, and the lateral porosity structure varies strongly with hillslope position. The aspect‐independent weathering structure suggests that the contemporary critical zone structure at Rancho Venada is a legacy of past climate and vegetation conditions.

    more » « less