Geology Reports⌕ Search

Geology topics

David L. George

Publications and source records attributed to David L. George.

33 records · Page 2Linked to original sources

Valid debris-flow models must avoid hot starts

Debris-flow experiments and models commonly use “hot-start” initial conditions in which downslope motion begins when a large force imbalance is abruptly imposed. By contrast, initiation of natural debris flows almost invariably results from small perturbations of static force balances that apply to debris masses poised in steep channels or on steep slopes. Models that neglect these static balances may violate physical law. Here we assess how the effects of hot starts are manifested in physical experiments, analytical dam-break models, and numerical models in which frictional resistance is too small to satisfy static force balances in debris-flow source areas. We then outline a numerical modeling framework that avoids use of hot starts. In this framework an initial static force balance is gradually perturbed by increasing pore-fluid pressure that may trigger the onset of debris motion. Subsequent increases in pore-fluid pressure, driven by debris motion, may then reduce the debris frictional strength, leading to high flow mobility.

Conference Paper↗

Seamless numerical simulation of a hazard cascade in which a landslide triggers a dam-breach flood and consequent debris flow

Numerical simulations of hazard cascades downstream from moraine-dammed lakes commonly must specify linkages between models of discrete processes such as wave overtopping, dam breaching, erosion, and downstream floods or debris flows. Such linkages can be rather arbitrary and can detract from the ability to accurately conserve mass and momentum during complex sequences of events. Here we describe an alternative methodology in which we use high-resolution lidar topography and 2-D, two-phase conservation laws to seamlessly simulate all stages of a hazard-cascade that culminates in a debris flow. Our simulations employ our depth-integrated numerical model D-Claw to evaluate hazards from prospective breaching of a moraine dam that impounds Carver Lake on the eastern flank of South Sister volcano in central Oregon, USA. We simulate a “worst-case scenario” sequence of events that begins with a hypothetical 1.6 million m3 landslide that originates near the summit of South Sister and enters Carver Lake. Wave generation and displacement of lake water then leads to dam overtopping, breach erosion, and a downstream debris flow that funnels into Whychus Creek and eventually reaches the community of Sisters, Oregon, about 20 km away. Notably, our simulations predict that much of the debris is directed away from Sisters as a result of natural avulsion and flow diversion that occurs near the head of a low-gradient alluvial fan upstream from Sisters. Consequently, predicted hazards to downtown Sisters are less severe than those predicted by 1-D shallow-water simulations of a Carver Lake dam breach that were performed in the 1980s.

Conference Paper↗

Overcoming barriers to progress in seismic monitoring and characterization of debris flows and lahars

Debris flows generate seismic signals that contain valuable information about events as they unfold. Though seismic waves have been used for along-channel debris-flow and lahar monitoring systems for decades, it has proven difficult to move beyond detection to more quantitative characterizations of flow parameters and event size. This is for two primary reasons: (1) our limited understanding of how the radiated wavefield relates to debris flow characteristics and dynamics, and (2) difficulties quantifying the effects of heterogeneous shallow earth structure on the observed wavefield. The latter issue, essentially our inability to sufficiently separate seismic path effects from source information, is a barrier to improving our understanding of the first issue. We review the progress that has been made toward establishing the theory, models and methods required to use seismic observations to make quantitative measurements of flows and summarize the practical, social, and scientific barriers to progress. We discuss some specific ongoing efforts to overcome some of these barriers, with a focus on how we are using large-scale seismic experiments at the U.S. Geological Survey debris-flow flume to develop methods for directly measuring path effects and to develop and validate theoretical debris flow seismicity models.

Conference Paper↗

Combining InSAR and GPS to determine transient movement and thickness of a seasonally active low-gradient translational landslide

The combined application of continuous Global Positioning System data (high temporal resolution) with spaceborne interferometric synthetic aperture radar data (high spatial resolution) can reveal much more about the complexity of large landslide movement than is possible with geodetic measurements tied to only a few specific measurement sites. This approach is applied to an ~4 km 2 reactivated translational landslide in the Columbia River Gorge (Washington State), which moves mainly during the winter rainy season. Results reveal the complex three-dimensional shape of the landslide mass, how onset of sliding relates to cumulative rainfall, how surface velocity during sliding varies with location on the topographically complex landslide surface, and how the ground surface subsides slightly in weeks prior to downslope sliding.

Geophysical Research Letters↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Parallelization of GeoClaw code for modeling geophysical flows with adaptive mesh refinement on many-core systems

We parallelized the GeoClaw code on one-level grid using OpenMP in March, 2011 to meet the urgent need of simulating tsunami waves at near-shore from Tohoku 2011 and achieved over 75% of the potential speed-up on an eight core Dell Precision T7500 workstation [1]. After submitting that work to SC11 the International Conference for High Performance Computing, we obtained an unreleased OpenMP version of GeoClaw from David George, who developed the GeoClaw code as part of his PH.D thesis. In this paper, we will show the complementary characteristics of the two approaches used in parallelizing GeoClaw and the speed-up obtained by combining the advantage of each of the two individual approaches with adaptive mesh refinement (AMR), demonstrating the capabilities of running GeoClaw efficiently on many-core systems. We will also show a novel simulation of the Tohoku 2011 Tsunami waves inundating the Sendai airport and Fukushima Nuclear Power Plants, over which the finest grid distance of 20 meters is achieved through a 4-level AMR. This simulation yields quite good predictions about the wave-heights and travel time of the tsunami waves.

Conference Paper↗

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↗

Reconnaissance of chemical and biological quality in the Owyhee River from the Oregon State line to the Owyhee Reservoir, Oregon, 2001&ndash;02

The Owyhee River drains an extremely rugged and sparsely populated landscape in northern Nevada, southwestern Idaho, and eastern Oregon. Most of the segment between the Oregon State line and Lake Owyhee is part of the National Wild and Scenic Rivers System, and few water-quality data exist for evaluating environmental impacts. As a result, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, assessed this river segment to characterize chemical and biological quality of the river, identify where designated beneficial uses are met and where changes in stream quality occur, and provide data needed to address activities related to environmental impact assessments and Total Maximum Daily Loads. Water-quality issues identified at one or more sites were water temperature, suspended sediment, dissolved oxygen, pH, nutrients, trace elements, fecal bacteria, benthic invertebrate communities, and periphyton communities. Generally, summer water temperatures routinely exceeded Oregon's maximum 7-day average criteria of 17.8 degrees Celsius. The presence of few coldwater taxa in benthic invertebrate communities supports this observation. Suspended-sediment concentrations during summer base flow were less than 10 milligrams per liter (mg/L). Dissolved solids concentrations ranged from 46 to 222 mg/L, were highest during base flow, and tended to increase in a downstream direction. Chemical compositions of water samples indicated that large proportions of upland-derived water extend to the lower reaches of the study area during spring runoff. Dissolved fluoride and arsenic concentrations were highest during base flow and may be a result of geothermal springs discharging to the river. No dissolved selenium was detected. Upstream from the Rome area, spring runoff concentrations of suspended sediment ranged from 0 to 52 mg/L, and all except at the Three Forks site were typically below 20 mg/L. Stream-bottom materials from the North Fork Owyhee River, an area with no mines, were enriched with nine trace elements, which indicates that this basin may be a natural source of these elements. Near Rome, the part of the study area not included in the National Wild and Scenic Rivers System, land-use impacts resulted in elevated populations of Escherichia coli bacteria ( E. coli ) during base flow and elevated concentrations of nitrogen and phosphorus during spring runoff. Sites in this area had the highest numbers of benthic invertebrates; the fewest Ephemeroptera, Plecoptera, and Trichoptera taxa; and the highest Hilsenhoff Biotic Index scores. These results suggest degraded stream quality. Periphyton communities at sites in this area approached nuisance levels and could cause significant dissolved oxygen depletions and pH values that exceed Oregon's recommended criteria. Stream-bottom materials from Jordan Creek were enriched with mercury and manganese, which probably were ultimately caused by past mining in that basin. Below Crooked Creek, elevated suspended sediment concentrations (142 mg/L), phosphorus concentrations (0.23 mg/L), and E. coli populations (370 most probable number per 100 milliliters) during the largest spring runoff event could be the result of inputs at the lower end of Jordan Valley and (or) inputs from Crooked Creek. The New Zealand Mud Snail, a highly competitive gastropod introduced to the Snake River in the 1980s, was collected just downstream from the Crooked Creek confluence.

Idaho;Nevada;Oregon↗