Geology topics

Richard M. Iverson

Publications and source records attributed to Richard M. Iverson.

At least 55 records · Page 3Linked to original sources

USGS2012

Grain-size segregation and levee formation in geophysical mass flows

Data from large-scale debris-flow experiments are combined with modeling of particle-size segregation to explain the formation of lateral levees enriched in coarse grains. The experimental flows consisted of 10 m 3 of water-saturated sand and gravel, which traveled ∼80 m down a steeply inclined flume before forming an elongated leveed deposit 10 m long on a nearly horizontal runout surface. We measured the surface velocity field and observed the sequence of deposition by seeding tracers onto the flow surface and tracking them in video footage. Levees formed by progressive downslope accretion approximately 3.5 m behind the flow front, which advanced steadily at ∼2 m s −1 during most of the runout. Segregation was measured by placing ∼600 coarse tracer pebbles on the bed, which, when entrained into the flow, segregated upwards at ∼6–7.5 cm s −1 . When excavated from the deposit these were distributed in a horseshoe-shaped pattern that became increasingly elevated closer to the deposit termination. Although there was clear evidence for inverse grading during the flow, transect sampling revealed that the resulting leveed deposit was strongly graded laterally, with only weak vertical grading. We construct an empirical, three-dimensional velocity field resembling the experimental observations, and use this with a particle-size segregation model to predict the segregation and transport of material through the flow. We infer that coarse material segregates to the flow surface and is transported to the flow front by shear. Within the flow head, coarse material is overridden, then recirculates in spiral trajectories due to size-segregation, before being advected to the flow edges and deposited to form coarse-particle-enriched levees.

Journal of Geophysical Research↗

USGS2011

Modeling hazardous mass flows Geoflows09: Mathematical and computational aspects of modeling hazardous geophysical mass flows; Seattle, Washington, 9–11 March 2009

A recent workshop at the University of Washington focused on mathematical and computational aspects of modeling the dynamics of dense, gravity-driven mass movements such as rock avalanches and debris flows. About 30 participants came from seven countries and brought diverse backgrounds in geophysics; geology; physics; applied and computational mathematics; and civil, mechanical, and geotechnical engineering. The workshop was cosponsored by the U.S. Geological Survey Volcano Hazards Program, by the U.S. National Science Foundation through a Vertical Integration of Research and Education (VIGRE) in the Mathematical Sciences grant to the University of Washington, and by the Pacific Institute for the Mathematical Sciences. It began with a day of lectures open to the academic community at large and concluded with 2 days of focused discussions and collaborative work among the participants.

Eos, Earth and Space Science News↗

USGS2011

Mount St. Helens: A 30-year legacy of volcanism

The spectacular eruption of Mount St. Helens on 18 May 1980 electrified scientists and the public. Photodocumentation of the colossal landslide, directed blast, and ensuing eruption column—which reached as high as 25 kilometers in altitude and lasted for nearly 9 hours—made news worldwide. Reconnaissance of the devastation spurred efforts to understand the power and awe of those moments (Figure 1). The eruption remains a seminal historical event—studying it and its aftermath revolutionized the way scientists approach the field of volcanology. Not only was the eruption spectacular, but also it occurred in daytime, at an accessible volcano, in a country with the resources to transform disaster into scientific opportunity, amid a transformation in digital technology. Lives lost and the impact of the eruption on people and infrastructure downstream and downwind made it imperative for scientists to investigate events and work with communities to lessen losses from future eruptions.

Washington↗

USGS2011

Geomorphic and hydrologic dynamics of zero-order basins

The 1987 International Symposium on Erosion and Sedimentation in the Pacific Rim, held August 3–7, 1987, in Corvallis, Oreg., included a special session on the geomorphic and hydrologic dynamics of zero-order drainage basins. “Zero-order basin” is one of several terms used to describe unchanneled swales or hollows that may occupy considerable areas of higher-order drainage basins. These basins serve as important accumulation sites for water, sediment, and accompanying contaminants that subsequently may enter streams. Moreover, in many terrains along the Pacific rim, zero-order basins are significant geologic hazard sites in which rapid landslides and debris flows commonly begin. Widespread recognition o f the importance of zero-order basins has occurred only within the last 15 years, and this special session was aimed at collecting and assessing some of the pertinent knowledge that has been gained to date.

Oregon↗

USGS2011

Entrainment of bed sediment by debris flows: results from large-scale experiments

When debris flows grow by entraining sediment, they can become especially hazardous owing to increased volume, speed, and runout. To investigate the entrainment process, we conducted eight largescale experiments in the USGS debris-flow flume. In each experiment, we released a 6 m3 water-saturated debris flow across a 47-m long, ~12-cm thick bed of partially saturated sediment lining the 31º flume. Prior to release, we used low-intensity overhead sprinkling and real-time monitoring to control the bed-sediment wetness. As each debris flow descended the flume, we measured the evolution of flow thickness, basal total normal stress, basal pore-fluid pressure, and sediment scour depth. When debris flows traveled over relatively dry sediment, net scour was minimal, but when debris flows traveled over wetter sediment (volumetric water content > 0.22), debris-flow volume grew rapidly and flow speed and runout were enhanced. Data from scour sensors showed that entrainment occurred by rapid (5-10 cm/s), progressive scour rather than by mass failure at depth. Overriding debris flows rapidly generated high basal pore-fluid pressures when they loaded and deformed bed sediment, and in wetter beds these pressures approached lithostatic levels. Reduction of intergranular friction within the bed sediment thereby enhanced scour efficiency, entrainment, and runout.

Conference Paper↗

USGS2011

A two-phase debris-flow model that includes coupled evolution of volume fractions, granular dilatancy, and pore-fluid pressure

Pore-fluid pressure plays a crucial role in debris flows because it counteracts normal stresses at grain contacts and thereby reduces intergranular friction. Pore-pressure feedback accompanying debris deformation is particularly important during the onset of debrisflow motion, when it can dramatically influence the balance of forces governing downslope acceleration. We consider further effects of this feedback by formulating a new, depth-averaged mathematical model that simulates coupled evolution of granular dilatancy, solid and fluid volume fractions, pore-fluid pressure, and flow depth and velocity during all stages of debris-flow motion. To illustrate implications of the model, we use a finite-volume method to compute one-dimensional motion of a debris flow descending a rigid, uniformly inclined slope, and we compare model predictions with data obtained in large-scale experiments at the USGS debris-flow flume. Predictions for the first 1 s of motion show that increasing pore pressures (due to debris contraction) cause liquefaction that enhances flow acceleration. As acceleration continues, however, debris dilation causes dissipation of pore pressures, and this dissipation helps stabilize debris-flow motion. Our numerical predictions of this process match experimental data reasonably well, but predictions might be improved by accounting for the effects of grain-size segregation.

Conference Paper↗

USGS2010

The perfect debris flow? Aggregated results from 28 large-scale experiments

Aggregation of data collected in 28 controlled experiments reveals reproducible debris-flow behavior that provides a clear target for model tests. In each experiment ∼10 m 3 of unsorted, water-saturated sediment composed mostly of sand and gravel discharged from behind a gate, descended a steep, 95-m flume, and formed a deposit on a nearly horizontal runout surface. Experiment subsets were distinguished by differing basal boundary conditions (1 versus 16 mm roughness heights) and sediment mud contents (1 versus 7 percent dry weight). Sensor measurements of evolving flow thicknesses, basal normal stresses, and basal pore fluid pressures demonstrate that debris flows in all subsets developed dilated, coarse-grained, high-friction snouts, followed by bodies of nearly liquefied, finer-grained debris. Mud enhanced flow mobility by maintaining high pore pressures in flow bodies, and bed roughness reduced flow speeds but not distances of flow runout. Roughness had these effects because it promoted debris agitation and grain-size segregation, and thereby aided growth of lateral levees that channelized flow. Grain-size segregation also contributed to development of ubiquitous roll waves, which had diverse amplitudes exhibiting fractal number-size distributions. Despite the influence of these waves and other sources of dispersion, the aggregated data have well-defined patterns that help constrain individual terms in a depth-averaged debris-flow model. The patterns imply that local flow resistance evolved together with global flow dynamics, contradicting the hypothesis that any consistent rheology applied. We infer that new evolution equations, not new rheologies, are needed to explain how characteristic debris-flow behavior emerges from the interactions of debris constituents.

Journal of Geophysical Research F: Earth Surface↗

USGS2010

Gravity-driven groundwater flow and slope failure potential: 1. Elastic effective-stress model

Hilly or mountainous topography influences gravity-driven groundwater flow and the consequent distribution of effective stress in shallow subsurface environments. Effective stress, in turn, influences the potential for slope failure. To evaluate these influences, we formulate a two-dimensional, steady state, poroelastic model. The governing equations incorporate groundwater effects as body forces, and they demonstrate that spatially uniform pore pressure changes do not influence effective stresses. We implement the model using two finite element codes. As an illustrative case, we calculate the groundwater flow field, total body force field, and effective stress field in a straight, homogeneous hillslope. The total body force and effective stress fields show that groundwater flow can influence shear stresses as well as effective normal stresses. In most parts of the hillslope, groundwater flow significantly increases the Coulomb failure potential Φ, which we define as the ratio of maximum shear stress to mean effective normal stress. Groundwater flow also shifts the locus of greatest failure potential toward the slope toe. However, the effects of groundwater flow on failure potential are less pronounced than might be anticipated on the basis of a simpler, one-dimensional, limit equilibrium analysis. This is a consequence of continuity, compatibility, and boundary constraints on the two-dimensional flow and stress fields, and it points to important differences between our elastic continuum model and limit equilibrium models commonly used to assess slope stability.

Water Resources Research↗

USGS2010

Gravity-driven groundwater flow and slope failure potential: 2. Effects of slope morphology, material properties, and hydraulic heterogeneity

Hillslope morphology, material properties, and hydraulic heterogeneities influence the role of groundwater flow in provoking slope instability. We evaluate these influences quantitatively by employing the elastic effective stress model and Coulomb failure potential concept described in our companion paper (Iverson and Reid, this issue). Sensitivity analyses show that of four dimensionless quantities that control model results (i.e., Poisson's ratio, porosity, topographic profile, and hydraulic conductivity contrast), slope profiles and hydraulic conductivity contrasts have the most pronounced and diverse effects on groundwater seepage forces, effective stresses, and slope failure potentials. Gravity-driven groundwater flow strongly influences the shape of equilibrium hillslopes, which we define as those with uniform near-surface failure potentials. For homogeneous slopes with no groundwater flow, equilibrium hillslope profiles are straight; but with gravity-driven flow, equilibrium profiles are concave or convex-concave, and the largest failure potentials exist near the bases of convex slopes. In heterogeneous slopes, relatively slight hydraulic conductivity contrasts of less than 1 order of magnitude markedly affect the seepage force field and slope failure potential. Maximum effects occur if conductivity contrasts are of four orders of magnitude or more, and large hydraulic gradients commonly result in particularly large failure potentials just upslope from where low-conductivity layers intersect the ground surface.

Water Resources Research↗

USGS2010

Groundwater Seepage Vectors and the Potential for Hillslope Failure and Debris Flow Mobilization

Insight for understanding the effect of groundwater flow on the potential for hillslope failure and liquefaction is provided by a novel limit‐equilibrium analysis of infinite slopes with steady, uniform Darcian seepage of arbitrary magnitude and direction. Normalization of the limit‐equilibrium solution shows that three dimensionless parameters govern completely the Coulomb failure potential of saturated, cohesionless, infinite homogeneous hillslopes: (1) the ratio of seepage force magnitude to gravitational body force magnitude; (2) the angle θ − Φ, where θ is the surface slope angle and Φ is the angle of internal friction of the soil; and (3) the angle λ + Φ, where λ is the angle of the seepage vector measured with respect to an outward‐directed surface‐normal vector. An additional dimensionless parameter affects the solution if soil cohesion is included in the analysis. Representation of the normalized solution as a single family of curves shows that minimum slope stability universally occurs when the seepage direction is given by λ = 90° − Φ. It also shows that for some upward seepage conditions, slope stability is limited by static liquefaction rather than by Coulomb failure. Close association between these liquefaction conditions and certain Coulomb failure conditions indicates that slope failure in such instances could be responsible for nearly spontaneous mobilization of destructive flowing soil masses on hillslopes.

Water Resources Research↗

USGS2010

Volcano collapse promoted by progressive strength reduction: New data from Mount St. Helens

Rock shear strength plays a fundamental role in volcano flank collapse, yet pertinent data from modern collapse surfaces are rare. Using samples collected from the inferred failure surface of the massive 1980 collapse of Mount St. Helens (MSH), we determined rock shear strength via laboratory tests designed to mimic conditions in the pre-collapse edifice. We observed that the 1980 failure shear surfaces formed primarily in pervasively shattered older dome rocks; failure was not localized in sloping volcanic strata or in weak, hydrothermally altered rocks. Our test results show that rock shear strength under large confining stresses is reduced &sim;20% as a result of large quasi-static shear strain, as preceded the 1980 collapse of MSH. Using quasi-3D slope-stability modeling, we demonstrate that this mechanical weakening could have provoked edifice collapse, even in the absence of transiently elevated pore-fluid pressures or earthquake ground shaking. Progressive strength reduction could promote collapses at other volcanic edifices.

Washington↗

USGS2010

Effects of soil aggregates on debris-flow mobilization: Results from ring-shear experiments

Rates and styles of landslide motion are sensitive to pore-water pressure changes caused by changes in soil porosity accompanying shear deformation. Soil may either contract or dilate upon shearing, depending upon whether its initial porosity is greater or less, respectively, than a critical-state porosity attained after sufficiently high strain. We observed complications in this behavior, however, during rate-controlled (0.02 m s −1 ) ring-shear experiments conducted on naturally aggregated dense loamy sand at low confining stresses (10.6 and 40 kPa). The aggregated soil first dilated and then contracted to porosities less than initial values, whereas the same soil with its aggregates destroyed monotonically dilated. We infer that aggregates persisted initially during shear and caused dilation before their eventual breakdown enabled net contraction. An implication of this contraction, demonstrated in experiments in which initial soil porosity was varied, is that the value of porosity distinguishing initially contractive from dilative behavior can be significantly larger than the critical-state porosity, which develops only after disaggregation ceases at high strains. In addition, post-dilative contraction may produce excess pore pressures, thereby reducing frictional strength and facilitating debris-flow mobilization. We infer that results of triaxial tests, which generally produce strains at least a factor of ∼ 4 smaller than those we observed at the inception of post-dilative contraction, do not allow soil contraction to be ruled out as a mechanism for debris-flow mobilization in dense soils containing aggregates.

Engineering Geology↗

USGS2008

An exact solution for ideal dam-break floods on steep slopes

The shallow‐water equations are used to model the flow resulting from the sudden release of a finite volume of frictionless, incompressible fluid down a uniform slope of arbitrary inclination. The hodograph transformation and Riemann's method make it possible to transform the governing equations into a linear system and then deduce an exact analytical solution expressed in terms of readily evaluated integrals. Although the solution treats an idealized case never strictly realized in nature, it is uniquely well‐suited for testing the robustness and accuracy of numerical models used to model shallow‐water flows on steep slopes.

Water Resources Research↗

USGS2008

Book review of Avalanche Dynamics by Shiva P. Pudasaini and Kolumban Hutter. Springer: Berlin-Heidelberg, 2007. 602 pages, 225 figures, 15 tables

This highly specialized book is interesting not only because of its important subject matter but also because of its egocentric perspective. The majority of the book provides a nearly exhaustive retrospective of the authors’ many contributions to the shallow‐flow theory of granular avalanches, and it also critiques contributions by others. Indeed, some readers (including this reviewer) might be dismayed to find their own contributions critiqued inaccurately! This brief review will not attempt to correct inaccuracies, however, but will emphasize the book’s main features.

Journal of Geology↗

USGS2008

Steady and intermittent slipping in a model of landslide motion regulated by pore-pressure feedback

This paper studies a parsimonious model of landslide motion, which consists of the one-dimensional diffusion equation (for pore pressure) coupled through a boundary condition to a first-order ODE (Newton's second law). Velocity weakening of sliding friction gives rise to nonlinearity in the model. Analysis shows that solutions of the model equations exhibit a subcritical Hopf bifurcation in which stable, steady sliding can transition to cyclical, stick-slip motion. Numerical computations confirm the analytical predictions of the parameter values at which bifurcation occurs. The existence of stick-slip behavior in part of the parameter space is particularly noteworthy because, unlike stick-slip behavior in classical models , here it arises in the absence of a reversible (elastic) driving force. Instead, the driving force is static (gravitational), mediated by the effects of pore-pressure diffusion on frictional resistance.

SIAM Journal on Applied Mathematics↗

USGS2008

Frictional properties of the Mount St. Helens gouge

Frictional properties of gouge bounding the solid dacite plug that extruded at Mount St. Helens during 2004 and 2005 may have caused stick-slip upward motion of the plug and associated seismicity. Laboratory experiments were performed with a ring-shear device to test the dependence of the peak and steady-state frictional strength of the gouge on shearing rate and hold time. A remolded gouge specimen (~0.012 m 3 ) was sheared under constant normal stresses ranging from 5 to 200 kPa and at rates ranging from 10 -6 to 10 -3 m/s. The gouge exhibited rate-weakening behavior at rates lower than 1×10 -4 m/s and rate-strengthening at rates above 5×10 -4 m/s. Peak strengths occurred during the onset of shearing, when displacements were generally less than 0.5 mm. In slide-holdslide tests, the peak strength of the gouge increased logarithmically as hold times increased from 3 s to almost 10 5 s. Rate-weakening friction is a requirement for stick-slip behavior that is satisfied by the Mount St. Helens gouge. Indeed, regular stick-slip oscillations were observed in two experiments performed at the highest normal stress and lowest rates of shear. The conditions under which this stick-slip motion occurred indicate that the gouge also satisfies a second criterion for stick-slip behavior of materials exhibiting rateand-state dependent friction-gouge stiffness exceeds that of the ascending magma that drives upward motion of the plug. The presence of highly compliant magma as a driving element may be crucial for generating stick-slip instabilities at the shallow earthquake focal depths observed during the eruption.

Washington↗

USGS2008

Dynamics of seismogenic volcanic extrusion resisted by a solid surface plug, Mount St. Helens, 2004-2005

The 2004-5 eruption of Mount St. Helens exhibited sustained, near-equilibrium behavior characterized by nearly steady extrusion of a solid dacite plug and nearly periodic occurrence of shallow earthquakes. Diverse data support the hypothesis that these earthquakes resulted from stick-slip motion along the margins of the plug as it was forced incrementally upward by ascending, solidifying, gas-poor magma. I formalize this hypothesis with a mathematical model derived by assuming that magma enters the base of the eruption conduit at a steady rate, invoking conservation of mass and momentum of the magma and plug, and postulating simple constitutive equations that describe magma and conduit compressibilities and friction along the plug margins. Reduction of the model equations reveals a strong mathematical analogy between the dynamics of the magma-plug system and those of a variably damped oscillator. Oscillations in extrusion velocity result from the interaction of plug inertia, a variable upward force due to magma pressure, and a downward force due to the plug weight. Damping of oscillations depends mostly on plug-boundary friction, and oscillations grow unstably if friction exhibits rate weakening similar to that observed in experiments. When growth of oscillations causes the extrusion rate to reach zero, however, gravity causes friction to reverse direction, and this reversal instigates a transition from unstable oscillations to self-regulating stick-slip cycles. The transition occurs irrespective of the details of rate-weakening behavior, and repetitive stick-slip cycles are, therefore, robust features of the system’s dynamics. The presence of a highly compressible elastic driving element (that is, magma containing bubbles) appears crucial for enabling seismogenic slip events to occur repeatedly at the shallow earthquake focal depths (<1 km) observed during the 2004-5 eruption. Computations show that fluctuations in magma pressure accompanying such slip events are <3 kPa, indicating that deviations from mechanical equilibrium are slight and that coseismic force drops are <10 8 N. These results imply that the system’s self-regulating behavior is not susceptible to dramatic change--provided that the rate of magma ascent remains similar to the rate of magma accretion at the base of the plug, that plug surface erosion more or less compensates for mass gain due to basal accretion, and that magma and rock properties do not change significantly. Even if disequilibrium initial conditions are imposed, the dynamics of the magma-plug system are strongly attracted to self-regulating stick-slip cycles, although this self-regulating behavior can be bypassed on the way to runaway behavior if the initial state is too far from equilibrium.

Washington↗

USGS2008

Mobility statistics and automated hazard mapping for debris flows and rock avalanches

Power-law equations that are physically motivated and statistically tested and calibrated provide a basis for forecasting areas likely to be inundated by debris flows, rock avalanches, and lahars with diverse volumes. The equations A=α 1 V 2/3 and B=α 2 V 2/3 are based on the postulate that the maximum valley cross-sectional area (A) and total valley planimetric area (B) likely to be inundated by a flow depend only on its volume (V) and the topography of the flow path. Testing of these equations involves determining whether or not they fit data for documented flows satisfactorily, and calibration entails determining best-fit values of the coefficients α 1 and α 2 for debris flows, rock avalanches, and lahars. This report describes statistical testing and calibration of the equations by using field data compiled from many sources, and it describes application of the equations to delineation of debris-flow hazard zones. Statistical results show that for each type of flow (debris flows, rock avalanches, and lahars), the dependence of A and B on V is described well by power laws with exponents equal to 2/3. This value of the exponent produces fits that are effectively indistinguishable from the best fits obtained by using adjustable power-law exponents. Statistically calibrated values of the coefficients α 1 and α 2 provide scale-invariant indices of the relative mobilities of rock avalanches (α 1 = 0.2, α 2 = 20), nonvolcanic debris flows (α 1 = 0.1, α 2 = 20), and lahars (α 1 = 0.05, α 2 = 200). These values show, for example, that a lahar of specified volume can be expected to inundate a planimetric area ten times larger than that inundated by a rock avalanche or nonvolcanic debris flow of the same volume. The utility of the calibrated debris-flow inundation equations A=0.1V 2/3 and B=20V 2/3 is demonstrated by using them within the GIS program LAHARZ to delineate nested hazard zones for future debris flows in an area bordering the Umpqua River in the south-central Oregon Coast Range. This application requires use of high-resolution topographic data derived form LIDAR surveys, knowledge of local geology to specify a suitable range of prospective debris-flow volumes, and development and use of a new algorithm for identification of prospective debris-flow source areas in finely dissected terrain.

Oregon↗