Search USGS⌕ Search

Geology topics

Richard M. Iverson

Publications and source records attributed to Richard M. Iverson.

At least 73 records · Page 4Linked to original sources
USGS2008

Deciphering landslide behavior using large-scale flume experiments

Landslides can be triggered by a variety of hydrologic events and they can exhibit a wide range of movement dynamics. Effective prediction requires understanding these diverse behaviors. Precise evaluation in the field is difficult; as an alternative we performed a series of landslide initiation experiments in the large-scale, USGS debris-flow flume. We systematically investigated the effects of three different hydrologic triggering mechanisms, including groundwater exfiltration from bedrock, prolonged rainfall infiltration, and intense bursts of rain. We also examined the effects of initial soil porosity (loose or dense) relative to the soil’s critical-state porosity. Results show that all three hydrologic mechanisms can instigate landsliding, but water pathways, sensor response patterns, and times to failure differ. Initial soil porosity has a profound influence on landslide movement behavior. Experiments using loose soil show rapid soil contraction during failure, with elevated pore pressures liquefying the sediment and creating fast-moving debris flows. In contrast, dense soil dilated upon shearing, resulting in slow, gradual, and episodic motion. These results have fundamental implications for forecasting landslide behavior and developing effective warning systems.

Conference Paper↗
USGS2007

Video documentation of experiments at the USGS debris-flow flume 1992–2024

This set of videos presents about 19.1 hours of footage documenting the 183 experiments conducted at the USGS debris-flow flume from 1992 to 2024. Owing to improvements in video technology over the years, the quality of footage from recent experiments generally exceeds that from earlier experiments. Use the link below to access the individual videos, which are mostly grouped by date and subject matter. When a video is selected from the list, multiple video sequences are generally shown in succession, beginning with a far-field overview and proceeding to close-up views and post-experiment documentation [ https://pubs.usgs.gov/of/2007/1315/ ]. Interpretations and data from experiments at the USGS debris-flow flume are not provided here but can be found in published reports, many of which are available online at www.usgs.gov . A brief introduction to the flume facility is also available online in USGS Open-File Report 92–483 [ https://doi.org/10.3133/ofr92483 ] .

Open-File Report↗
USGS2006

Forecasting runout of rock and debris avalanches

Physically based mathematical models and statistically based empirical equations each may provide useful means of forecasting runout of rock and debris avalanches. This paper compares the foundations, strengths, and limitations of a physically based model and a statistically based forecasting method, both of which were developed to predict runout across three-dimensional topography. The chief advantage of the physically based model results from its ties to physical conservation laws and well-tested axioms of soil and rock mechanics, such as the Coulomb friction rule and effective-stress principle. The output of this model provides detailed information about the dynamics of avalanche runout, at the expense of high demands for accurate input data, numerical computation, and experimental testing. In comparison, the statistical method requires relatively modest computation and no input data except identification of prospective avalanche source areas and a range of postulated avalanche volumes. Like the physically based model, the statistical method yields maps of predicted runout, but it provides no information on runout dynamics. Although the two methods differ significantly in their structure and objectives, insights gained from one method can aid refinement of the other.

Conference Paper↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 2. Experimental tests

Scaling considerations indicate that miniature experiments can be used to test models of granular avalanches in which the effects of intergranular fluid and cohesion are negligible. To test predictions of a granular avalanche model described in a companion paper, we performed bench top experiments involving avalanches of dry sand across irregular basal topography that mimicked the complexity of natural terrain. The experiments employed a novel method of laser-assisted cartography to map the three-dimensional morphology of rapidly moving avalanches, thereby providing high-resolution data for comparison with model output. Model input consisted of two material properties (angles of internal and basal Coulomb friction of the sand), which were measured in independent tests, and of initial and boundary conditions that characterized the geometry of the experimental apparatus. Experimental results demonstrate that the model accurately predicts not only the gross behavior but also many details of avalanche motion, from initiation to deposition. We attribute this accuracy to a mathematical and computational formulation that conserves mass and momentum in three-dimensional physical space and satisfies the Coulomb equation in three-dimensional stress space. Our results support the hypothesis that a Coulomb proportionality between shear and normal stresses applies in moderately rapid granular flows and that complicated constitutive postulates are unnecessary if momentum conservation is strictly enforced in continuum avalanche models. Furthermore, predictions of our Coulomb continuum model contrast with those of a Coulomb point mass model, illustrating the importance of multidimensional modeling and model testing.

Journal of Geophysical Research Earth Surface↗
USGS2004

Granular avalanches across irregular three-dimensional terrain: 1. Theory and computation

To establish a theoretical basis for predicting and interpreting the behavior of rapid mass movements on Earth's surface, we develop and test a new computational model for gravity-driven motion of granular avalanches across irregular, three-dimensional (3-D) terrain. The principles embodied in the model are simple and few: continuum mass and momentum conservation and intergranular stress generation governed by Coulomb friction. However, significant challenges result from the necessity of satisfying these principles when deforming avalanches interact with steep and highly variable 3-D terrain. We address these challenges in four ways. (1) We formulate depth-averaged governing equations that are referenced to a rectangular Cartesian coordinate system (with z vertical) and that account explicitly for the effect of nonzero vertical accelerations on depth-averaged mass and momentum fluxes and stress states. (2) We compute fluxes of mass and momentum across vertical cell boundaries using a high-resolution finite volume method and Roe-type Riemann solver. Our algorithm incorporates flux difference splitting, an entropy correction for the flux, and eigenvector decomposition to embed the effects of driving and resisting forces in Riemann solutions. (3) We use a finite element method and avalanche displacements predicted by Riemann solutions to compute Coulomb stresses conjugate to the displacements in 3-D stress space. (4) We test the model output against analytical solutions, a sand cone conceptual experiment, and (in a companion paper) data from detailed laboratory experiments. Model results illustrate a complex interplay of basal traction and internal stress, and they successfully predict not only the gross behavior but also many details of avalanche motion from initiation to deposition.

Journal of Geophysical Research Earth Surface↗
USGS2004

Comment on “Piezometric response in shallow bedrock at CB1: Implications for runoff generation and landsliding” by David R. Montgomery, William E. Dietrich, and John T. Heffner

Piezometric responses to rainfall on hillslopes commonly dictate the timing of landsliding. Insight to this phenomenon can be gained by evaluating the timescales for pore pressure perturbations to propagate normal and parallel to the ground surface, and these timescales can be estimated using characteristic values of hydraulic diffusivity [ Iverson , 2000 ]. However, Montgomery et al. [2002] employed an erroneous definition of hydraulic diffusivity, leading to flawed assessment of the piezometric response timescales identified by Iverson [2000] . This comment aims to rectify the errors of Montgomery et al. [2002] and clarify the methods used by Iverson [2000] to estimate piezometric response timescales at the CB1 field site studied by Montgomery et al.

Water Resources Research↗
USGS2003

Gravity-driven mass flows

Gravity-driven mass flows, also known as sediment gravity flows, include a spectrum of phenomena in which more-or-less coherent mixtures of grains and intergranular fluid flow down slopes. At one end of this spectrum are dilute flows in which momentum is transferred mostly by fluid forces and sediment is largely a passive cargo that increases the effective fluid density. These dilute mass flows are part of a larger class of fluid dynamical phenomena known as gravity currents ( Simpson, 1987 ; Bonnecaze et al., 1993 ). At the other end of the spectrum are concentrated mass flows such as granular avalanches and debris flows. Forces exerted by interacting solid grains dominate momentum transfer in such flows, although fluid forces can mediate grain interactions ( Savage and Hutter, 1989 ; Iverson, 1997 ).

Book chapter↗
USGS2001

Mechanics of debris flows and debris-laden flash floods

A new mathematical model developed to predict behavior of debris flows and avalanches also holds promise for predicting behavior of debris-laden flash floods. The model assumes that debris flows behave as mixtures of interacting Newtonian fluids and Coulomb solids. Solid and fluid constituents obey three-dimensional mass and momentum balances, which are summed and depth-integrated to yield equations that describe shallow flows of the mixture as a whole. An important distinction between these mixture equations and standard shallow-water equations results from strong variation of flow resistance due to interacting solid and fluid forces. Partitioning of flow resistance between solid and fluid components depends on fluid pressure, which evolves as flow evolves. If fluid pressure supports the total weight of the flowing mass, all resistance results from hydrodynamic forces, and the equations reduce to the conventional shallow-water form. If fluid pressure supports none of the weight of the flowing mass, all flow resistance results from Coulomb friction between interacting solids, and the equations describe motion of granular avalanches. A combination of solid and fluid resistance typifies debris flows and debris-laden flash floods. In these flows solid resistance commonly is concentrated at the fronts of advancing bores that may be heavily freighted with rocks and woody debris. Riemann methods provide an effective tool for solving the shallow flow equations numerically and predicting unsteady motion of debris flows and flash floods along paths with arbitrary geometry and inclination

Conference Paper↗
USGS2000

Landslide triggering by rain infiltration

Landsliding in response to rainfall involves physical processes that operate on disparate timescales. Relationships between these timescales guide development of a mathematical model that uses reduced forms of Richards equation to evaluate effects of rainfall infiltration on landslide occurrence, timing, depth, and acceleration in diverse situations. The longest pertinent timescale is A / D 0 , where D 0 is the maximum hydraulic diffusivity of the soil and A is the catchment area that potentially affects groundwater pressures at a prospective landslide slip surface location with areal coordinates x , y and depth H . Times greater than A / D 0 are necessary for establishment of steady background water pressures that develop at ( x , y , H ) in response to rainfall averaged over periods that commonly range from days to many decades. These steady groundwater pressures influence the propensity for landsliding at ( x , y , H ), but they do not trigger slope failure. Failure results from rainfall over a typically shorter timescale H 2 / D 0 associated with transient pore pressure transmission during and following storms. Commonly, this timescale ranges from minutes to months. The shortest timescale affecting landslide responses to rainfall is √(H/g), where g is the magnitude of gravitational acceleration. Postfailure landslide motion occurs on this timescale, which indicates that the thinnest landslides accelerate most quickly if all other factors are constant. Effects of hydrologic processes on landslide processes across these diverse timescales are encapsulated by a response function, R(t*) = √(t*/π) exp (-1/t*) - erfc (1/√t*), which depends only on normalized time, t*. Use of R(t*) in conjunction with topographic data, rainfall intensity and duration information, an infinite-slope failure criterion, and Newton's second law predicts the timing, depth, and acceleration of rainfall-triggered landslides. Data from contrasting landslides that exhibit rapid, shallow motion and slow, deep-seated motion corroborate these predictions.

Water Resources Research↗
USGS1999

Experimental testing of flexible barriers for containment of debris flows

In June 1996, six experiments conducted at the U.S. Geological Survey Debris Flow Flume demonstrated that flexible, vertical barriers constructed of wire rope netting can stop small debris flows. All experimental debris flows consisted of water-saturated gravelly sand with less than two percent finer sediment by weight. All debris flows had volumes of about 10 cubic meters, masses of about 20 metre tons, and impact velocities of 5 to 9 meters per second. In four experiments, the debris flow impacted pristine, unreformed barriers of varying design; in the other two experiments, the debris flow impacted barriers already loaded with sediment from a previous flow. Differences in barrier design led to differences in barrier performance. Experiments were conducted with barriers constructed of square-mesh wire-rope netting with 30centimeter, 20centimeter, and 15 centimeter mesh openings as well as 30centimeter diameter interlocking steel rings. In all cases, sediment cascading downslope at the leading edge of the debris flows tended to spray through the nets. Nets fitted with finer-mesh chain link or chicken wire liners contained more sediment than did unlined nets, and a ring net fitted with a synthetic silt screen liner contained nearly 100 percent of the sediment. Irreversible net displacements of up to 2 meters and friction brake engagement on the support and anchor cables dissipated some of the impact energy. However, substantial forces developed in the steel support columns and the lateral and tie-back anchor cables attached to these columns. As predicted by elementary mechanics, the anchor cables experienced larger tensile forces when the support columns were hinged at the base rather than bolted rigidly to the foundation. Measured loads in the lateral anchor cables exceeded those in the tie-back anchor cables and the load cell capacity of 45 kilo-Newtons. Measurements also indicated that the peak loads in the tie- back anchors were highly transient and occurred at the points of maximum momentum impulse to the net.

Open-File Report↗
USGS1998

Objective delineation of lahar-inundation hazard zones

A new method of delineating lahar hazard zones in valleys that head on volcano flanks provides a rapid, objective, reproducible alternative to traditional methods. The rationale for the method derives from scaling analyses of generic lahar paths and statistical analyses of 27 lahar paths documented at nine volcanoes. Together these analyses yield semiempirical equations that predict inundated valley cross-sectional areas (A) and planimetric areas (B) as functions of lahar volume (V) . The predictive equations ( A = 0.05 V 2/3 and B = 200 V 2/3 ) provide all information necessary to calculate and plot inundation limits on topographic maps. By using a range of prospective lahar volumes to evaluate A and B , a range of inundation limits can be plotted for lahars of increasing volume and decreasing probability. Resulting hazard maps show graphically that lahar-inundation potentials are highest near volcanoes and along valley thalwegs, and diminish gradually as distances from volcanoes and elevations above valley floors increase. We automate hazard-zone delineation by embedding the predictive equations in a geographic information system (GIS) computer program that uses digital elevation models of topography. Lahar hazard zones computed for Mount Rainier, Washington, mimic those constructed on the basis of intensive field investigations. The computed hazard zones illustrate the potentially widespread impact of large lahars, which on average inundate planimetric areas 20 times larger than those inundated by rock avalanches of comparable volume.

Washington↗
USGS1997

Debris-flow mobilization from landslides

Field observations, laboratory experiments, and theoretical analyses indicate that landslides mobilize to form debris flows by three processes: ( a ) widespread Coulomb failure within a sloping soil, rock, or sediment mass, ( b ) partial or complete liquefaction of the mass by high pore-fluid pressures, and ( c ) conversion of landslide translational energy to internal vibrational energy (i.e. granular temperature). These processes can operate independently, but in many circumstances they appear to operate simultaneously and synergistically. Early work on debris-flow mobilization described a similar interplay of processes but relied on mechanical models in which debris behavior was assumed to be fixed and governed by a Bingham or Bagnold rheology. In contrast, this review emphasizes models in which debris behavior evolves in response to changing pore pressures and granular temperatures. One-dimensional infinite-slope models provide insight by quantifying how pore pressures and granular temperatures can influence the transition from Coulomb failure to liquefaction. Analyses of multidimensional experiments reveal complications ignored in one-dimensional models and demonstrate that debris-flow mobilization may occur by at least two distinct modes in the field.

Annual Review of Earth and Planetary Sciences↗
USGS1997

Debris-flow initiation experiments using diverse hydrologic triggers

Controlled debris-flow initiation experiments focused on three hydrologic conditions that can trigger slope failure: localized ground-water inflow; prolonged moderate-intensity rainfall; and high-intensity rainfall. Detailed monitoring of slope hydrology and deformation provided exceptionally complete data on conditions preceding and accompanying slope failure and debris-flow mobilization. Ground-water inflow and high-intensity sprinkling led to abrupt, complete failure whereas moderate-intensity sprinkling led to retrogressive, block-by-block failure. Failure during ground-water inflow and during moderate-intensity sprinkling occurred with a rising water table and positive pore pressures. Failure during high-intensity sprinkling occurred without widespread positive pore pressures. In all three cases, pore pressures in most locations increased dramatically (within 2-3 seconds) during failure. In some places, pressures in unsaturated materials rapidly 'flashed' from zero to elevated positive values. Transiently elevated pore pressures and partially liquefied soil enhanced debris-flow mobilization.

Conference Paper↗
USGS1997

Automated, reproducible delineation of zones at risk from inundation by large volcanic debris flows

Large debris flows can pose hazards to people and property downstream from volcanoes. We have developed a rapid, reproducible, objective, and inexpensive method to delineate distal debris-flow hazard zones. Our method employs the results of scaling and statistical analyses of the geometry of volcanic debris flows (lahars) to predict inundated valley cross-sectional areas (A) and planimetric areas (B) as functions of lahar volume. We use a range of specified lahar volumes to evaluate A and B. In a Geographic Information System (GIS) we employ the resulting range of predicted A and B to delineate gradations in inundation hazard, which is highest near the volcano and along valley thalwegs and diminishes as distances from the volcano and elevations above valley floors increase. Comparison of our computer-generated hazard maps with those constructed using traditional, field-based methods indicates that our method can provide an accurate means of delineating lahar hazard zones.

Conference Paper↗