Geology Reports⌕ Search

Geology topics

David L. George

Publications and source records attributed to David L. George.

At least 19 recordsLinked to original sources

Modeling floods, sediment entrainment, and downstream debris flows from hypothetical breaches of the blockage at Spirit Lake, Washington

This report describes a modeling investigation by the U.S. Geological Survey (USGS) of hazards in the Toutle and Cowlitz River valleys posed by hypothetical outburst floods from Spirit Lake, Washington. A massive debris avalanche resulting from the collapse of Mount St. Helens’ north flank during the May 18, 1980, eruption blocked Spirit Lake’s natural outlet into the North Fork Toutle River. Lacking a natural outlet, subsequent runoff in the Spirit Lake watershed contributed to a rising lake level, elevating the potential for debris-dam breaching or catastrophic failure. The influence of highly erodible bed sediment in the upper North Fork Toutle River on downstream flood and debris-flow dynamics and extent is assessed in this study. Simulations of clear-water (non-erosive) outburst floods were used as a baseline and compared to erosive flows that entrain large volumes of material and transition into debris flows along their flow path, revealing the influence of entrainment on hazard extent. Clear-water floods were modeled with the shallow water equations. Erosive flows were modeled with a two-phase granular fluid model that accommodates mobilization and incorporation of sediment from the bed into the overlying flow and resultant changes in flow rheology across a wide range of solid concentrations, from dilute suspensions to dense-granular debris flows. Entrainment of bed material was found to substantially increase the total flow volume (total volume of transported water and sediment is approximately 150 percent of the water volume for non-erosive flows). Erosive flows are shown to exhibit higher flow-front speeds and faster downstream arrival times than non-erosive flows, consistent with volume amplification effects near the actively mobilizing flow front. However, the larger total volume of transported material does not necessarily lead to an enhancement of total volume throughput (cumulative discharge) or inundation extent (total affected area) for all locations along the entire flow path; while entrainment leads to the displacement of a larger volume of material overall, much of this dislocated material (water and sediment) deposits upstream from the distal extent of the flows. These results are consistent with energetic considerations of initial potential energy and granular shear resistance.

Washington↗

Using the D-Claw software package to model lahars in the Middle Fork Nooksack River drainage and beyond, Mount Baker, Washington

Lahars, or volcanic mudflows, are the most hazardous eruption-related phenomena that will affect communities living along rivers that originate on Mount Baker. In the past 15,000 years, the largest lahars from Mount Baker have affected the Middle Fork Nooksack River drainage and beyond. Here we use the physics-based D-Claw software package to model nine lahar scenarios that are initiated as water-saturated landslides between Sherman Crater and the Roman Wall on the Mount Baker edifice and flow down the Middle Fork Nooksack River. The scenarios range in volume from 1 to 260 million cubic meters and have an initial hydraulic permeability from 10 −12 to 10 −10 meters squared. Model output includes data such as flow depth, velocity, runout distance, area inundated, arrival time, and sediment concentration as well as information that allows scientists to calculate other important hydrologic characteristics such as lahar discharge. These data are important to officials who have the responsibility to plan for, or take mitigation measures against, future Mount Baker lahars. To check the validity of the D-Claw results, we compare the scenarios to known geologic information. We also compare D-Claw results with empirical models that have been used in the past to determine potential inundation areas, runout distances, and arrival times. These comparisons highlight similarities and differences between empirical and physics-based models. We also present D-Claw scenario-based animations to help scientists, officials, and lay people alike to visualize how future lahars could affect communities.

Washington↗

Utility of a swath laser rangefinder for characterizing mass movement flow depth and landslide initiation

Mass movements such as debris flows and landslides are some of the deadliest and most destructive natural hazards occurring mostly in alpine and volcanic settings. With ever-growing populations located downslope from known debris flow channels, early warning systems can help prevent loss of life. Geophysical and technological advances have improved monitoring and detection capabilities in recent years; however, they can often be cost prohibitive and resource intensive, making them less accessible to disadvantaged populations. We tested and validated a readily available and cost-effective two-dimensional swath laser rangefinder in a controlled experimental setting against two independent flow-depth lasers. The swath laser successfully recorded cross-sectional changes in flow depth from four debris flows and a water-only flood, in addition to geomorphic changes associated with landslide initiation. The results suggest that a swath laser could be integrated into systems for debris flow detection and characterization of mass movements in natural settings, thus improving the ability to monitor these hazards.

Landslides↗

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗

Forecasting inundation of catastrophic landslides from precursory creep

Forecasting landslide inundation upon catastrophic failure is crucial for reducing casualties, yet it remains a long-standing challenge owing to the complex nature of landslides. Recent global studies indicate that catastrophic hillslope failures are commonly preceded by a period of precursory creep, motivating a novel scheme to foresee their hazard. Here, we showcase an approach to hindcast landslide inundation by linking satellite-captured precursory displacements to modeling of consequent granular-fluid flows. We present its application to the 2021 Chunchi, Ecuador landslide, which failed catastrophically and evolved into a mobile debris flow after four months of precursory creep, destroying 68 homes along its lengthy flow path. Underpinned by uncertainty quantification and in situ validations, we highlight the feasibility and potential of forecasting landslide inundation damage using observable precursors. This forecast approach is broadly applicable for flow hazards initiated from geomaterial failures.

Geophysical Research Letters↗

Numerical modeling of debris flows: A conceptual assessment

Real-world hazard evaluation poses many challenges for the development and application of numerical models of debris flows. In this chapter we provide a conceptual overview of physically based, depth-averaged models designed to simulate debris-flow motion across three-dimensional terrain. When judiciously formulated and applied, these models can provide useful information about anticipated depths, speeds, and extents of debris-flow inundation as well as debris interactions with structures such as levees and dams. Depth-averaged debris-flow models can differ significantly from one another, however. Some of the greatest differences result from simulation of one-phase versus two-phase flow, use of parsimonious versus information-intensive initial and boundary conditions, use of tuning coefficients versus physically measureable parameters, application of dissimilar numerical solution techniques, and variations in computational speed and model accessibility. This overview first addresses these and related attributes of depth-averaged debris-flow models. It then describes model testing and application to hazard evaluation, with a focus on our own model, D-Claw. The overview concludes with a discussion of outstanding challenges for development of improved debris-flow models and suggestions for prospective model users.

Book chapter↗

Forecasting the inundation of postfire debris flows

In the semi-arid regions of the western United States, postfire debris flows are typically runoff generated. The U.S. Geological Survey has been studying the mechanisms of postfire debris-flow initiation for multiple decades to generate operational models for forecasting the timing, location, and magnitude of postfire debris flows. Here we discuss challenges and progress for extending operational capabilities to include modeling postfire debris-flow inundation extent. Analysis of volume and impacted area scaling relationships indicated that postfire debris flows do not conform to assumptions of geometric self-similarity. We documented sensitivity of impacted areas to rainfall intensity using a candidate methodology for generating inundation hazard assessments. Our results emphasize the importance of direct measurements of debris-flow volume, inundated area, and high temporal resolution rainfall intensity.

California↗

Simulating debris flow and levee formation in the 2D shallow flow model D-Claw: Channelized and unconfined flow

Debris flow runout poses a hazard to life and infrastructure. The expansion of human population into mountainous areas and onto alluvial fans increases the need to predict and mitigate debris flow runout hazards. Debris flows on unconfined alluvial fans can exhibit spontaneous self-channelization through levee formation that reduces lateral spreading and extends runout distances compared to unchannelized flows. Here we modify the D-Claw shallow flow model in two ways that are hypothesized to generate levees. We evaluate these modifications with observations from a large-scale flume experiment. We investigate model performance when including the effect of two different friction sub-models, as well as the inclusion of segregation effects on granular permeability. Results show that, for a wide range of plausible model input parameters, simulations including the effects of segregation promoted modeled levee formation, whereas simulations without the effects of segregation did not create levees. Further, using a forward predictive framework, simulations with the effects of segregation were more likely to better model the magnitude of debris flow depth and runout distance, whereas simulation timing of the debris flow was affected by the choice of friction sub-model. Our results indicate that including the effects of segregation on granular permeability can improve the likelihood of better predictions of debris flow depth and runout prior to an event occurring.

Earth and Space Science↗

Insights on multistage rock avalanche behavior from runout modeling constrained by seismic inversions

Inversion of low-frequency regional seismic records to solve for a time series of bulk forces exerted on the earth by a landslide (a force-time function) is increasingly being used to infer volumes and dynamics of large, highly energetic landslides, such as rock avalanches and flowslides, and to provide calibration information on event dynamics and volumes for numerical landslide runout models. Much of the work to date using landslide runout modeling constrained by seismic data has focused on using single-phase models with frictional or velocity-weakening rheologies. Awareness of multistage landslide initiations is increasing, with discrete failures separated in time contributing to the final impact of an event. Our work utilizes a method for incorporating seismic data as a calibration constraint for landslide runout models, considering variable rheologies and different initiation conditions. This study presents a systematic examination of multiple rheologies and initiation conditions, and shows how these factors affect the force-time function derived from the landslide runout model. Our work confirms that, while rheology and fragmenting or initially coherent initiations affect the force-time function, multiple collapses separated by tens of seconds have the greatest impact on the shape and amplitude. We apply this method to the analysis of three real rock avalanches to better constrain plausible initiation conditions and rheology parameters using both seismic and field data. This study provides insights on how assumptions about the initiation dynamics of the source zone and the runout model definition can aid in the interpretation of seismic inversions for multistage rock avalanches.

Journal of Geophysical Research: Solid Earth↗

Modeling the dynamics of lahars that originate as landslides on the west side of Mount Rainier, Washington

Large lahars pose substantial threats to people and property downstream from Mount Rainier volcano in Washington State. Geologic evidence indicates that these threats exist even during the absence of volcanic activity and that the threats are highest in the densely populated Puyallup and Nisqually River valleys on the west side of the volcano. However, the precise character of these threats can be difficult to anticipate. To help predict depths and rates of possible lahar inundation in the area, this report presents the results of simulations of hypothetical future lahars that originate high on the west side of Mount Rainier and travel downstream into the Puyallup and Nisqually River valleys. Many of the results portrayed as still images in the figures of this report are also available as animated files that can be accessed at the web address provided in the figure captions. We simulated eight scenarios, including worst-case scenarios in which the simulated lahars are similar in size and mobility to the approximately 260 million cubic meter (Mm 3 ; 340 million cubic yard) Electron Mudflow lahar that descended from Mount Rainier and inundated the Puyallup River valley about 500 years ago. The other six scenarios place the worst-case scenarios in perspective by simulating lahars that originate from the same source areas but have smaller volumes or lesser mobilities. We perform our simulations using an open-source software package that we developed called D-Claw. The numerical model composing the kernel of D-Claw solves a system of five hyperbolic partial differential equations that describe the depth-averaged dynamics of static or flowing grain-fluid mixtures interacting with three-dimensional topography. In D-Claw, the volume fraction occupied by solid grains is a dependent variable that can freely evolve, enabling simulation of landslide liquefaction and of lahar interaction with static bodies of water. The latter feature facilitates a seamless simulation of a lahar in the Nisqually River valley entering Alder Lake reservoir. In the event of an approximately 260 Mm 3 high-mobility lahar originating on the west side of Mount Rainier, our results point to two areas of pronounced hazard. One area, comprising the densely populated lowlands of Orting, Washington, and environs, could be inundated by lahars originating from either the Sunset Amphitheater or Tahoma Glacier headwall areas. In the worst-case scenario we consider for the Orting lowlands, which involves a 260 Mm 3 high-mobility lahar originating from a landslide in the Sunset Amphitheater, a flow front approximately 4 meters deep and traveling about 4 meters per second reaches the Orting lowlands about 1 hour after the onset of slope failure. After passing through the Orting lowlands, the simulated lahar slows down and comes to rest in the valleys surrounding Sumner and Puyallup. A second area of pronounced hazard is the stretch of the Nisqually River valley beginning in Mount Rainier National Park and extending downstream to Alder Lake reservoir and Alder Dam. This area would be substantially affected in the worst-case scenario that involves a 260 Mm 3 high-mobility lahar originating from the Tahoma Glacier headwall area—the locality identified by a previous study as the sector of Mount Rainier most prone to large-scale gravitational collapse. The simulated lahar passes through the area of Ashford, Washington, within about 20 minutes of the onset of slope failure and reaches the head of Alder Lake within about 50 minutes. The lahar ultimately displaces enough reservoir water to cause overtopping of the 100 meter (330 foot) tall Alder Dam, but consequences of such dam overtopping are not addressed in this report.

Washington↗

Multi-model comparison of computed debris flow runout for the 9 January 2018 Montecito, California post-wildfire event

Hazard assessment for post-wildfire debris flows, which are common in the steep terrain of the western United States, has focused on the susceptibility of upstream basins to generate debris flows. However, reducing public exposure to this hazard also requires an assessment of hazards in downstream areas that might be inundated during debris flow runout. Debris flow runout models are widely available, but their application to hazard assessment for post-wildfire debris flows has not been extensively tested. Necessary inputs to these models include the total volume of the mobilized flow, flow properties (either inherent material properties or calibration coefficients), and site topography. Estimates of volume are possible in post-event (“back calculation”) studies, yet before an event, volume is an uncertain quantity. We simulated debris flow runout for the well-constrained 9 January 2018 Montecito event using three models (RAMMS, FLO2D, and D-Claw) to determine the relative importance of volume and flow properties. We broke the impacted area into three domains, and for each model-domain combination, we performed a numerical sampling study in which volume and flow properties varied within a wide, but plausible range. We assessed model performance based on inundation patterns and peak flow depths. We found all models could simulate the event with comparable results. Simulation performance was most sensitive to flow volume and less sensitive to flow properties. Our results emphasize the importance of reducing uncertainty in pre-event estimates of flow volume for hazard assessment.

California↗

Preliminary assessment of the wave generating potential from landslides at Barry Arm, Prince William Sound, Alaska

We simulated the concurrent rapid motion of landslides on an unstable slope at Barry Arm, Alaska. Movement of landslides into the adjacent fjord displaced fjord water and generated a tsunami, which propagated out of Barry Arm. Rather than assuming an initial sea surface height, velocity, and location for the tsunami, we generated the tsunami directly using a model capable of simulating the dynamics of both water and landslide material. The fjord below most of the landslide source area was occupied by the Barry Glacier until about 2012; therefore, our direct simulation of tsunami generation by landslide motion required new topographic and bathymetric data, which was collected in 2020. The topographic data also constrained landslide geometries and volumes. We considered four scenarios based on two landslide volumes and two landslide mobilities—a more mobile, contractive landslide and a less mobile, noncontractive landslide. The larger of the two volumes is 689 × 10 6 cubic meters (m 3 )—larger than the volume estimate in a previous study—and reflects the largest plausible volume given current observational data. The considered scenario that generated the largest wave heights resulted in forecast wave heights of over 200 meters (m) in the northern part of Barry Arm, adjacent to the landslide source area and runup on the opposite fjord wall in excess of 500 m. Simulated wave heights in excess of 5 m in southern Barry Arm and in Harriman Fjord occurred within 10–15 minutes (min) of landslide motion. The simulated tsunami reached Whittier, Alaska, approximately 20 min after initial rapid landslide motion, with peak heights of just over 2 m in Passage Fjord, 500 m offshore Whittier, occurring 26 min after initial rapid motion. Time of peak wave heights was consistent with previous modeling. Although results are preliminary and can be refined with additional observations and analyses, they provide a refined assessment of the upper bound of the hazard presented by the Barry Arm landslides. The results herein support the National Oceanic and Atmospheric Administration’s National Tsunami Warning Center mission to detect, forecast, and warn for tsunamis in Alaska.

Alaska↗

When hazard avoidance is not an option: Lessons learned from monitoring the postdisaster Oso landslide, USA

On 22 March 2014, a massive, catastrophic landslide occurred near Oso, Washington, USA, sweeping more than 1 km across the adjacent valley flats and killing 43 people. For the following 5 weeks, hundreds of workers engaged in an exhaustive search, rescue, and recovery effort directly in the landslide runout path. These workers could not avoid the risks posed by additional large-scale slope collapses. In an effort to ensure worker safety, multiple agencies cooperated to swiftly deploy a monitoring and alerting system consisting of sensors, automated data processing and web-based display, along with defined communication protocols and clear calls to action for emergency management and search personnel. Guided by the principle that an accelerating landslide poses a greater threat than a steadily moving or stationary mass, the system was designed to detect ground motion and vibration using complementary monitoring techniques. Near real-time information was provided by continuous GPS, seismometers/geophones, and extensometers. This information was augmented by repeat-assessment techniques such as terrestrial and aerial laser scanning and time-lapse photography. Fortunately, no major additional landsliding occurred. However, we did detect small headscarp failures as well as slow movement of the remaining landslide mass with the monitoring system. This was an exceptional response situation and the lessons learned are applicable to other landslide disaster crises. They underscore the need for cogent landslide expertise and ready-to-deploy monitoring equipment, the value of using redundant monitoring techniques with distinct goals, the benefit of clearly defined communication protocols, and the importance of continued research into forecasting landslide behavior to allow timely warning.

Washington↗

Landslide monitoring and runout hazard assessment by integrating multi-source remote sensing and numerical models: An application to the Gold Basin landslide complex, northern Washington

The landslide complex at Gold Basin, Washington, has been drawing considerable attention after a catastrophic runout of the nearby landslide at Oso, Washington, in 2014. To evaluate potential threats of the Gold Basin landslide to the campground down the slope, remote sensing and numerical modeling were integrated to monitor recent landslide activity and simulate hypothetical runout scenarios. Bare-earth LiDAR DEM (digital elevation model) differencing, InSAR (Interferometric Synthetic Aperture Radar), and offset tracking of SAR images reveal that localized collapses at the headscarps have been the primary type of landslide activity at Gold Basin from 2005 to 2019, and currently no signs indicative of movement of a large centralized block or a deep-seated main body were detected. The maximum horizontal deformation rate is 5 m/year occurring primarily from headscarp recession of the middle lobe, and the annual landsliding volume of the whole landslide complex averages 1.03 × 10 5 m 3 . From three-dimensional limit equilibrium analysis of generalized terrace structures, the maximum landslide volume is estimated as 2.0 × 10 6 m 3 . Simulations of hypothetical runout scenarios were carried out using the depth-averaged two-phase model D-claw with above-obtained landslide geometry constraints. The simulation results demonstrate that debris flows with volume less than 10 5 m 3 only pose limited threats to the campground, while volumes over 10 6 m 3 could cause severe damages. Consequently, the estimated maximum landslide volume of 2.0 × 10 6 m 3 suggests a potential risk to the campground nearby. Adaption of our methodology could prove useful for evaluating other similar landslides globally for hazards prevention and mitigation.

Washington↗

The Missoula and Bonneville floods—A review of ice-age megafloods in the Columbia River basin

The Channeled Scabland of eastern Washington State, USA, brought megafloods to the scientific forefront. A 30,000-km2 landscape of coulees and cataracts carved into the region’s loess-covered basalt attests to overwhelming volumes of energetic water. The scarred landscape, garnished by huge boulder bars and far-travelled ice-rafted erratics, spurred J Harlen Bretz’s vigorously disputed flood hypothesis in the 1920s. First known as the Spokane flood, it was rebranded the Missoula flood once understood that the water came from glacial Lake Missoula, formed when the Purcell Trench lobe of the last-glacial Cordilleran ice sheet dammed the Clark Fork valley in northwestern Idaho with ice a kilometer thick. Bretz’s flood evidence in the once-remote Channeled Scabland, widely seen and elaborated by the 1950s, eventually swayed consensus for cataclysmic flooding. Missoula flood questions then turned to some that continue today: how many? when? how big? what routes? what processes? The Missoula floods passed through eastern Washington by a multitude of valleys, coulees and scabland tracts, some contemporaneously, some sequentially. Which routings and their timing depended on the positions of various lobes of the multi-pronged Cordilleran ice sheet and the erosional development of the channels themselves. The first floods mostly followed the big bend of Columbia valley looping through north-central Washington. But the south-advancing Okanogan ice lobe soon blocked that path, forming long-lasting glacial Lake Columbia in the impounded Columbia valley. Missoula floods into this lake were diverted south out of the Columbia valley and into eastern Washington coulees and scabland tracts. At least four floods entered Moses Coulee, but then as the Okanogan lobe advanced over and blocked the head of that coulee, more eastern paths took the water, including Grand Coulee and the Telford-Crab-Creek and Cheney-Palouse scabland tracts. Flood routing also depended on the erosion of the coulees. At some point, headward erosion of upper Grand Coulee lowered the divide saddle between the west-running Columbia valley and the deep and wide Grand Coulee heading southwest. Still uncertain is when this happened and the consequences with respect to the stage and extent of glacial Lake Columbia and to flood access to the other, higher, flood routes. Downstream, all flood routes converged onto Pasco Basin, flowed through Wallula Gap and the Columbia River Gorge into the Pacific Ocean, following submarine canyons and depositing sediment layers on abyssal plains. Stratigraphic studies indicate dozens—likely more than a hundred—separate Missoula floods during the last glacial period. Over the length of the flood route, backwater areas and depositional basins preserve multiple flood beds, many of which are separated by signs of time, including volcanic ash layers and soil development in subaerial environments; and varve-like beds and pelagic mud layers in lacustrine and marine settings. Evidence also comes from the glacial Lake Missoula basin, where stratigraphy indicates dozens of filling and emptying cycles. Varve counts in conjunction of radiocarbon dating and paleomagnetic secular variation show the repeated filling-and-release cycles of glacial Lake Missoula had intervals possibly as long as 100 years early in the lake’s history but diminished to just one or two years for the last few floods. This behavior accords with jökulhlaup-style floods released by subglacial drainage from a self-dumping ice-dammed lake. But not yet clear is whether such a mechanism applies to all the floods or if some emptied more cataclysmically as hypothesized by some. Radiocarbon dating of sparse organic materials remains key to defining flood chronology but has been lately bolstered by analyses of terrestrial cosmogenic nuclides and optically stimulated luminescence. Varve counts and paleomagnetic secular variation studies help to define durations and intervals represented by sequences of flood beds. The ~16 ka Mount St. Helens Set S tephra is commonly interbedded within flood deposits, enabling correlation of deposits among sites. Tephra from the 13.7–13.4 ka eruption of Glacier Peak overlies all glacial Lake Missoula and Missoula flood deposits, defining an end time. Overall conclusions are that glacial Lake Missoula was extant and producing floods for at least 3–4 ky during 20–14 ka. At least ~75 floods preceded Mount St Helens Set S, followed by 30 or more after the tephra fall. Most floods entered glacial Lake Columbia, impounded by the Okanogan lobe, for 2–5 ky between about 18.5 and 15 ka. Glacial Lake Columbia outlived Lake Missoula by >200–400 yr but may have been born later since at least one flood came down the Columbia valley before the Okanogan ice lobe blocked the Columbia valley at 18.5–18 ka. The maximum extent of the Okanogan and Purcell Trench lobes, many Missoula floods, substantial erosion of upper Grand Coulee, and the widespread tephra falls from Mount St. Helens eruptions all happened about 17–15 ka. People, in the area since 16.6–15.3 ka, almost certainly witnessed the last of the Missoula floods and later large floods from other ice-dammed lakes in the Columbia River basin. Quantitative flow analyses give peak discharge estimates and support understanding of erosional and depositional processes. The first flow assessments were simple cross-section calculations but recent assessments employ two-dimensional hydrodynamic models. The general finding is that emplacement of the maximum stage evidence requires about 20 million m3/s near the Lake Missoula outlet and about 5–15 million m3/s through Wallula Gap and downstream in the Columbia River Gorge. These hydraulic analyses raise still-unresolved questions regarding canyon erosion and possible additional water sources. The large Pleistocene Bonneville flood entered the Columbia River system from the southeast from pluvial Lake Bonneville, the Pleistocene predecessor to Great Salt Lake in the eastern Great Basin. During the last glacial, the lake basin filled, covering >50,000 km2 with 10,400 km3 of water before reaching its maximum possible stage governed by Red Rock Pass, the lowest divide separating the basin from the Snake River basin to the north. The overtopping lake rapidly incised 108–125 m into the Red Rock Pass outlet, spilling half of its total lake volume. G.K. Gilbert described the essential sequence in the 1870s, but the flood was mostly forgotten until the late 1950s when Harold Malde linked the spectacular scabland topography and bouldery “melon gravel” on the Snake River Plain to the Lake Bonneville overflow. The Bonneville flood appears to have been a singular event at about 18 ka. No evidence of multiple or pre-last-glacial spillovers has yet been found. Its total volume was about twice that of a maximum Lake Missoula flood yet its peak discharge was ~1 million m3/s, less than a tenth of the largest Missoula floods. Its comparatively simple flow path and much steadier flow make the Bonneville flood ideal for new studies of erosional and depositional processes. At least two floods seem to have passed down the Columbia valley after the last of the Missoula floods, including a large flood about ~14 ka likely from cataclysmic demise of the thinning Okanogan ice lobe dam impounding glacial Lake Columbia. Floods from earlier glacial ages left scant yet clear evidence in the Channeled Scabland and Columbia valley. But their source, timing, and magnitudes are little understood. Some deposits are paleomagnetically reversed, thus older than ~800 ka. Last-glacial floods and perhaps older ones affected the Snake River Plain, some likely sourced in lakes dammed by alpine glaciers in central Idaho.

Idaho, Oregon, Washington↗

Diverse cataclysmic floods from Pleistocene glacial Lake Missoula

In late Wisconsin time, the Purcell Trench lobe of the Cordilleran ice sheet dammed the Clark Fork of the Columbia River in western Montana, creating glacial Lake Missoula. During part of this epoch, the Okanogan lobe also dammed the Columbia River downstream, creating glacial Lake Columbia in northeast Washington. Repeated failure of the Purcell Trench ice dam released glacial Lake Missoula, causing dozens of catastrophic floods in eastern Washington that can be distinguished by the geologic record they left behind. These floods removed tens of meters of pale loess from dark basalt substrate, forming scars along flowpaths visible from space. Different positions of the Okanogan lobe are required for modeled Missoula floods to inundate the diverse channels that show field evidence for flooding, as shown by accurate dam-break flood modeling using a roughly 185 m digital terrain model of existing topography (with control points dynamically varied using automatic mesh refinement). The maximum extent of the Okanogan lobe, which blocked inundation of the upper Grand Coulee and the Columbia River valley, is required to flood all channels in the Telford scablands and to produce highest flood stages in Pasco Basin. Alternatively, the Columbia River valley must have been open and the upper Grand Coulee blocked to nearly match evidence for high water on Pangborn bar near Wenatchee, Washington, and to flood Quincy Basin from the west. Finally, if the Columbia River valley and upper Grand Coulee were both open, Quincy Basin would have flooded from the northeast. In all these scenarios, the discrepancy between modeled flood stages and field evidence for maximum flood stages increases in all channels downstream, from Spokane to Umatilla Basin. The pattern of discrepancies indicates that bulking of floods by loess increased flow volume across the scablands, but this alone does not explain low modeled flow stages along the Columbia River valley near Wenatchee. This latter discrepancy between modeled flood stages and field data requires either additional bulking of flow by sediment along the Columbia reach downstream of glacial Lake Columbia, or coincident dam failures of glacial Lake Columbia and glacial Lake Missoula.

Montana, Washington↗

Basal stress equations for granular debris masses on smooth or discretized slopes

Knowledge of basal stresses is essential for analyzing slope stability and modeling the dynamics and erosive potential of debris flows and avalanches. Here we derive and test new algebraic formulas for calculating the shear stress τ and normal stress σ at the base of variable‐thickness granular debris masses in states of static or dynamic equilibrium on slopes. The formulas include a lateral pressure coefficient κ , but use of a fixed value κ = 0.7 yields predictions of σ that on average err by less than 3% and of τ that on average err by less than 13% in matching basal stresses measured in six large‐scale experiments involving wet debris masses with varying geometries and compositions. Much larger prediction errors result from use of infinite‐slope or shallow‐debris approximations. Specialized versions of the new formulas apply if basal topography is discretized and represented by a “staircase” function in a digital elevation model. Use of these formulas to assess static limiting equilibrium conditions shows that the apparent basal Coulomb friction angle ϕ tread of debris that engages friction acting on the horizontal surfaces (or “treads”) of a staircase sloping at an angle θ is generally described by tan ϕ tread = tan ( ϕ − θ )+ κ tan θ , where ϕ is the true basal friction angle of the same debris in contact with a uniformly sloping bed. Differences between the values of ϕ and ϕ tread can greatly influence the results of numerical simulations that use unsmoothed digital elevation model topography to calculate the stability or dynamics of debris masses on slopes.

Journal of Geophysical Research-Earth Surface↗