Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

• Published in: Journal of Geophysical Research: Earth Surface Vol. 115; no. F3
• Main Authors: Baum, Rex L., Godt, Jonathan W., Savage, William Z.
• Format: Journal Article
• Language: English
• Published: Washington, DC Blackwell Publishing Ltd 01.09.2010; American Geophysical Union
• Subjects: Aeration zone; distributed model; Earth sciences; Earth, ocean, space; Exact sciences and technology; Geographic information systems; Hydrology; Infiltration; Landslides; Landslides & mudslides; Mathematical models; Moisture content; Physical properties; Position (location); Rainfall; rainfall infiltration; Remote sensing; Satellite navigation systems; Slope stability; Soil water; Spatial distribution; Unsaturated; Water table; Water tables; United States; Washington; geographic information systems; digital simulation; slopes; North America; Modeling; rupture; pore pressure; safety; Diffusion; soils; Analytical solution; landscapes; efficiency; rainfall; Initial condition; water table; infiltration; unsaturated zone; numerical models; landslides; transient phenomena; Timing; saturated zone; frame structure; slope stability
• ISSN: 0148-0227; 2169-9003; 2156-2202; 2169-9011
• DOI: 10.1029/2009JF001321

Shallow rainfall‐induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two‐layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure‐diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil‐water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one‐dimensional slope‐stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall‐induced pore‐pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.