Geology Reports⌕ Search

Geology topics

Richard M. Iverson

Publications and source records attributed to Richard M. Iverson.

At least 37 records · Page 2Linked to original sources
USGS2017

New methodology for computing tsunami generation by subaerial landslides: Application to the 2015 Tyndall Glacier landslide, Alaska

Landslide-generated tsunamis pose significant hazards and involve complex, multiphase physics that are challenging to model. We present a new methodology in which our depth-averaged two-phase model D-Claw is used to seamlessly simulate all stages of landslide dynamics as well as tsunami generation, propagation, and inundation. Because the model describes the evolution of solid and fluid volume fractions, it treats both landslides and tsunamis as special cases of a more general class of phenomena. Therefore, the landslide and tsunami can be efficiently simulated as a single-layer continuum with evolving solid-grain concentrations, and with wave generation via direct longitudinal momentum transfer—a dominant physical mechanism that has not been previously addressed in this manner. To test our methodology, we used D-Claw to model a large subaerial landslide and resulting tsunami that occurred on 17 October 2015, in Taan Fjord near the terminus of Tyndall Glacier, Alaska. Modeled shoreline inundation patterns compare well with those observed in satellite imagery.

Alaska↗
USGS2016

Debris flow runup on vertical barriers and adverse slopes

Runup of debris flows against obstacles in their paths is a complex process that involves profound flow deceleration and redirection. We investigate the dynamics and predictability of runup by comparing results from large-scale laboratory experiments, four simple analytical models, and a depth-integrated numerical model (D-Claw). The experiments and numerical simulations reveal the important influence of unsteady, multidimensional flow on runup, and the analytical models highlight key aspects of the underlying physics. Runup against a vertical barrier normal to the flow path is dominated by rapid development of a shock, or jump in flow height, associated with abrupt deceleration of the flow front. By contrast, runup on sloping obstacles is initially dominated by a smooth flux of mass and momentum from the flow body to the flow front, which precedes shock development and commonly increases the runup height. D-Claw simulations that account for the emergence of shocks show that predicted runup heights vary systematically with the adverse slope angle and also with the Froude number and degree of liquefaction (or effective basal friction) of incoming flows. They additionally clarify the strengths and limitations of simplified analytical models. Numerical simulations based on a priori knowledge of the evolving dynamics of incoming flows yield quite accurate runup predictions. Less predictive accuracy is attained in ab initio simulations that compute runup based solely on knowledge of static debris properties in a distant debris flow source area. Nevertheless, the paucity of inputs required in ab initio simulations enhances their prospective value in runup forecasting.

Journal of Geophysical Research F: Earth Surface↗
USGS2016

Comment on “The reduction of friction in long-runout landslides as an emergent phenomenon” by Brandon C. Johnson et al.

Results from a highly idealized, 2-D computational model indicate that dynamic normal-stress rarefactions might cause friction reduction in long-runout landslides, but the physical relevance of the idealized dynamics has not been confirmed by experimental tests. More importantly, the model results provide no evidence that refutes alternative hypotheses about friction reduction mechanisms. One alternative hypothesis, which is strongly supported by field evidence, experimental data, and the predictions of a well-constrained computational model, involves development of high pore fluid pressures in deforming landslide material or overridden bed material. However, no scientific basis exists for concluding that a universal mechanism is responsible for friction reduction in all long-runout landslides.

Journal of Geophysical Research F: Earth Surface↗
USGS2016

Discussion of “The relation between dilatancy, effective stress and dispersive pressure in granular avalanches” by P. Bartelt and O. Buser (DOI: 10.1007/s11440-016-0463-7)

A paper recently published by Bartelt and Buser (hereafter identified as “the authors”) aims to clarify relationships between granular dilatancy and dispersive pressure and to question the effective stress principle and its application to shallow granular avalanches (Bartelt and Buser in Act Geotech 11:549–557, 2 ). The paper also criticizes our own recent work, which utilizes the concepts of evolving dilatancy and effective stress to model the initiation and dynamics of water-saturated landslides and debris flows. Here we first explain why we largely agree with the authors’ views of dilatancy and dispersive pressure as they apply to depth-integrated granular avalanche models, and why we disagree with their views of effective stress and pore-fluid pressure. We conclude by explaining why the authors’ characterization of our recently developed D-Claw model is inaccurate.

Acta Geotechnica↗
USGS2016

Clawpack: Building an open source ecosystem for solving hyperbolic PDEs

Clawpack is a software package designed to solve nonlinear hyperbolic partial differential equations using high-resolution finite volume methods based on Riemann solvers and limiters. The package includes a number of variants aimed at different applications and user communities. Clawpack has been actively developed as an open source project for over 20 years. The latest major release, Clawpack 5, introduces a number of new features and changes to the code base and a new development model based on GitHub and Git submodules. This article provides a summary of the most significant changes, the rationale behind some of these changes, and a description of our current development model.

PeerJ Computer Science↗
USGS2016

Modelling landslide liquefaction, mobility bifurcation and the dynamics of the 2014 Oso disaster

Some landslides move slowly or intermittently downslope, but others liquefy during the early stages of motion, leading to runaway acceleration and high-speed runout across low-relief terrain. Mechanisms responsible for this disparate behaviour are represented in a two-phase, depth-integrated, landslide dynamics model that melds principles from soil mechanics, granular mechanics and fluid mechanics. The model assumes that gradually increasing pore-water pressure causes slope failure to nucleate at the weakest point on a basal slip surface in a statically balanced mass. Failure then spreads to adjacent regions as a result of momentum exchange. Liquefaction is contingent on pore-pressure feedback that depends on the initial soil state. The importance of this feedback is illustrated by using the model to study the dynamics of a disastrous landslide that occurred near Oso, Washington, USA, on 22 March 2014. Alternative simulations of the event reveal the pronounced effects of a landslide mobility bifurcation that occurs if the initial void ratio of water-saturated soil equals the lithostatic, critical-state void ratio. They also show that the tendency for bifurcation increases as the soil permeability decreases. The bifurcation implies that it can be difficult to discriminate conditions that favour slow landsliding from those that favour liquefaction and long runout.

Washington↗
USGS2015

Controls on the breach geometry and flood hydrograph during overtopping of non-cohesive earthen dams

Overtopping failure of non-cohesive earthen dams was investigated in 13 large-scale experiments with dams built of compacted, damp, fine-grained sand. Breaching was initiated by cutting a notch across the dam crest and allowing water escaping from a finite upstream reservoir to form its own channel. The channel developed a stepped profile, and upstream migration of the steps, which coalesced into a headcut, led to the establishment of hydraulic control (critical flow) at the channel head, or breach crest, an arcuate erosional feature that functions hydraulically as a weir. Novel photogrammetric methods, along with underwater videography, revealed that the retreating headcut maintained a slope near the angle of friction of the sand, while the cross section at the breach crest maintained a geometrically similar shape through time. That cross-sectional shape was nearly unaffected by slope failures, contrary to the assumption in many models of dam breaching. Flood hydrographs were quite reproducible--for sets of dams ranging in height from 0.55 m to 0.98 m--when the time datum was chosen as the time that the migrating headcut intersected the breach crest. Peak discharge increased almost linearly as a function of initial dam height. Early-time variability between flood hydrographs for nominally identical dams is probably a reflection of subtle experiment-to-experiment differences in groundwater hydrology and the interaction between surface water and groundwater.

Water Resources Research↗
USGS2015

Entrainment of bed material by Earth-surface mass flows: review and reformulation of depth-integrated theory

Earth-surface mass flows such as debris flows, rock avalanches, and dam-break floods can grow greatly in size and destructive potential by entraining bed material they encounter. Increasing use of depth-integrated mass- and momentum-conservation equations to model these erosive flows motivates a review of the underlying theory. Our review indicates that many existing models apply depth-integrated conservation principles incorrectly, leading to spurious inferences about the role of mass and momentum exchanges at flow-bed boundaries. Model discrepancies can be rectified by analyzing conservation of mass and momentum in a two-layer system consisting of a moving upper layer and static lower layer. Our analysis shows that erosion or deposition rates at the interface between layers must in general satisfy three jump conditions. These conditions impose constraints on valid erosion formulas, and they help determine the correct forms of depth-integrated conservation equations. Two of the three jump conditions are closely analogous to Rankine-Hugoniot conditions that describe the behavior of shocks in compressible gasses, and the third jump condition describes shear traction discontinuities that necessarily exist across eroding boundaries. Grain-fluid mixtures commonly behave as compressible materials as they undergo entrainment, because changes in bulk density occur as the mixtures mobilize and merge with an overriding flow. If no bulk density change occurs, then only the shear-traction jump condition applies. Even for this special case, however, accurate formulation of depth-integrated momentum equations requires a clear distinction between boundary shear tractions that exist in the presence or absence of bed erosion.

Reviews of Geophysics↗
USGS2015

Lahars and their deposits

Lahars occur during volcanic eruptions--or, less predictably, through other processes on steep volcanic terrain--when large masses of water mixed with sediment sweep down and off volcano slopes and commonly incorporate additional sediment and water. Because lahars are water-saturated, both liquid and solid interactions influence their behavior and distinguish them from other related phenomena common to volcanoes, such as debris avalanches and floods. The rock fragments carried by lahars make them especially destructive; the abundant liquid contained in them allows them to flow over gentle gradients and inundate areas far away from their sources. People in such distal areas commonly neither expect the danger nor anticipate the destructive power of lahars.

Book chapter↗
USGS2015

Landslide mobility and hazards: implications of the 2014 Oso disaster

Landslides reflect landscape instability that evolves over meteorological and geological timescales, and they also pose threats to people, property, and the environment. The severity of these threats depends largely on landslide speed and travel distance, which are collectively described as landslide &ldquo;mobility&rdquo;. To investigate causes and effects of mobility, we focus on a disastrous landslide that occurred on 22 March 2014 near Oso, Washington, USA, following a long period of abnormally wet weather. The landslide's impacts were severe because its mobility exceeded that of prior historical landslides at the site, and also exceeded that of comparable landslides elsewhere. The &sim;8&times;10 6 m 3 landslide originated on a gently sloping (<20&deg;) riverside bluff only 180 m high, yet it traveled across the entire &sim;1 km breadth of the adjacent floodplain and spread laterally a similar distance. Seismological evidence indicates that high-speed, flowing motion of the landslide began after about 50 s of preliminary slope movement, and observational evidence supports the hypothesis that the high mobility of the landslide resulted from liquefaction of water-saturated sediment at its base. Numerical simulation of the event using a newly developed model indicates that liquefaction and high mobility can be attributed to compression- and/or shear-induced sediment contraction that was strongly dependent on initial conditions. An alternative numerical simulation indicates that the landslide would have been far less mobile if its initial porosity and water content had been only slightly lower. Sensitive dependence of landslide mobility on initial conditions has broad implications for assessment of landslide hazards.

Washington↗
USGS2015

Scaling and design of landslide and debris-flow experiments

Scaling plays a crucial role in designing experiments aimed at understanding the behavior of landslides, debris flows, and other geomorphic phenomena involving grain-fluid mixtures. Scaling can be addressed by using dimensional analysis or &ndash; more rigorously &ndash; by normalizing differential equations that describe the evolving dynamics of the system. Both of these approaches show that, relative to full-scale natural events, miniaturized landslides and debris flows exhibit disproportionately large effects of viscous shear resistance and cohesion as well as disproportionately small effects of excess pore-fluid pressure that is generated by debris dilation or contraction. This behavioral divergence grows in proportion to H 3 , where H is the thickness of a moving mass. Therefore, to maximize geomorphological relevance, experiments with wet landslides and debris flows must be conducted at the largest feasible scales. Another important consideration is that, unlike stream flows, landslides and debris flows accelerate from statically balanced initial states. Thus, no characteristic macroscopic velocity exists to guide experiment scaling and design. On the other hand, macroscopic gravity-driven motion of landslides and debris flows evolves over a characteristic time scale ( L / g ) 1/2 , where g is the magnitude of gravitational acceleration and L is the characteristic length of the moving mass. Grain-scale stress generation within the mass occurs on a shorter time scale, H /( gL ) 1/2 , which is inversely proportional to the depth-averaged material shear rate. A separation of these two time scales exists if the criterion H / L < < 1 is satisfied, as is commonly the case. This time scale separation indicates that steady-state experiments can be used to study some details of landslide and debris-flow behavior but cannot be used to study macroscopic landslide or debris-flow dynamics.

Geomorphology↗
USGS2014

A depth-averaged debris-flow model that includes the effects of evolving dilatancy. I. Physical basis

To simulate debris-flow behaviour from initiation to deposition, we derive a depth-averaged, two-phase model that combines concepts of critical-state soil mechanics, grain-flow mechanics and fluid mechanics. The model's balance equations describe coupled evolution of the solid volume fraction, m , basal pore-fluid pressure, flow thickness and two components of flow velocity. Basal friction is evaluated using a generalized Coulomb rule, and fluid motion is evaluated in a frame of reference that translates with the velocity of the granular phase, v s . Source terms in each of the depth-averaged balance equations account for the influence of the granular dilation rate, defined as the depth integral of ∇⋅ v s . Calculation of the dilation rate involves the effects of an elastic compressibility and an inelastic dilatancy angle proportional to m − m eq , where m eq is the value of m in equilibrium with the ambient stress state and flow rate. Normalization of the model equations shows that predicted debris-flow behaviour depends principally on the initial value of m − m eq and on the ratio of two fundamental timescales. One of these timescales governs downslope debris-flow motion, and the other governs pore-pressure relaxation that modifies Coulomb friction and regulates evolution of m . A companion paper presents a suite of model predictions and tests.

Proceedings of the Royal Society A↗
USGS2014

A depth-averaged debris-flow model that includes the effects of evolving dilatancy: II. Numerical predictions and experimental tests.

We evaluate a new depth-averaged mathematical model that is designed to simulate all stages of debris-flow motion, from initiation to deposition. A companion paper shows how the model&rsquo;s five governing equations describe simultaneous evolution of flow thickness, solid volume fraction, basal pore-fluid pressure, and two components of flow momentum. Each equation contains a source term that represents the influence of state-dependent granular dilatancy. Here we recapitulate the equations and analyze their eigenstructure to show that they form a hyperbolic system with desirable stability properties. To solve the equations we use a shock-capturing numerical scheme with adaptive mesh refinement, implemented in an open-source software package we call D-Claw. As tests of D-Claw, we compare model output with results from two sets of large-scale debris-flow experiments. One set focuses on flow initiation from landslides triggered by rising pore-water pressures, and the other focuses on downstream flow dynamics, runout, and deposition. D-Claw performs well in predicting evolution of flow speeds, thicknesses, and basal pore-fluid pressures measured in each type of experiment. Computational results illustrate the critical role of dilatancy in linking coevolution of the solid volume fraction and pore-fluid pressure, which mediates basal Coulomb friction and thereby regulates debris-flow dynamics.

Proceedings of the Royal Society A↗
USGS2014

Debris flows: Behavior and hazard assessment

Debris flows are water-laden masses of soil and fragmented rock that rush down mountainsides, funnel into stream channels, entrain objects in their paths, and form lobate deposits when they spill onto valley floors. Because they have volumetric sediment concentrations that exceed 40 percent, maximum speeds that surpass 10 m/s, and sizes that can range up to ~109 m3, debris flows can denude slopes, bury floodplains, and devastate people and property. Computational models can accurately represent the physics of debris-flow initiation, motion and deposition by simulating evolution of flow mass and momentum while accounting for interactions of debris' solid and fluid constituents. The use of physically based models for hazard forecasting can be limited by imprecise knowledge of initial and boundary conditions and material properties, however. Therefore, empirical methods continue to play an important role in debris-flow hazard assessment.

Geology Today↗
USGS2013

Discussion: Numerical study on the entrainment of bed material into rapid landslides

A paper recently published in this journal (Pirulli & Pastor, 2012) uses numerical modelling to study the important problem of entrainment of bed material by landslides. Unfortunately, some of the basic equations employed in the study are flawed, because they violate the principle of linear momentum conservation. Similar errors exist in some other studies of entrainment, and the errors appear to stem from confusion about the role of bed-sediment inertia in differing frames of reference.

Geotechnique↗