USGS Science⌕ Search

SEARCH · USGS Science

Results for “Landslides”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58Linked to original sources

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report↗

Scoops3D: software to analyze 3D slope stability throughout a digital landscape

The computer program, Scoops3D, evaluates slope stability throughout a digital landscape represented by a digital elevation model (DEM). The program uses a three-dimensional (3D) method of columns approach to assess the stability of many (typically millions) potential landslides within a user-defined size range. For each potential landslide (or failure), Scoops3D assesses the stability of a rotational, spherical slip surface encompassing many DEM cells using a 3D version of either Bishop’s simplified method or the Ordinary (Fellenius) method of limit-equilibrium analysis. Scoops3D has several options for the user to systematically and efficiently search throughout an entire DEM, thereby incorporating the effects of complex surface topography. In a thorough search, each DEM cell is included in multiple potential failures, and Scoops3D records the lowest stability (factor of safety) for each DEM cell, as well as the size (volume or area) associated with each of these potential landslides. It also determines the least-stable potential failure for the entire DEM. The user has a variety of options for building a 3D domain, including layers or full 3D distributions of strength and pore-water pressures, simplistic earthquake loading, and unsaturated suction conditions. Results from Scoops3D can be readily incorporated into a geographic information system (GIS) or other visualization software. This manual includes information on the theoretical basis for the slope-stability analysis, requirements for constructing and searching a 3D domain, a detailed operational guide (including step-by-step instructions for using the graphical user interface [GUI] software, Scoops3D-i) and input/output file specifications, practical considerations for conducting an analysis, results of verification tests, and multiple examples illustrating the capabilities of Scoops3D. Easy-to-use software installation packages are available for the Windows or Macintosh operating systems; these packages install the compiled Scoops3D program, the GUI (Scoops3D-i), and associated documentation. Several Scoops3D examples, including all input and output files, are available as well. The source code is written in the Fortran 90 language and can be compiled to run on any computer operating system with an appropriate compiler.

Techniques and Methods↗

Usoi Dam wave overtopping and flood routing in the Bartang and Panj Rivers, Tajikistan

The Usoi dam was created in the winter of 1911 after an enormous seismogenic rock slide completely blocked the valley of the Bartang River in the Pamir Mountains of southeastern Tajikistan. At present the dam impounds 17 million cubic meters of water in Lake Sarez. Flood volume and discharge estimates were made for several landslide generated floods that could overtop the dam. For landslide volumes of 200, 500, and 1,000 million cubic meters, estimated overtopping flood volumes were 2, 22, and 87 million cubic meters of water, respectively. Estimated peak discharge at the dam for these three flood scenarios were 57,000, 490,000, and 1,580,000 cubic meters per second, based on triangular hydrographs of 70-, 90-, and 110-second durations, respectively. Flood-routing simulations were made for the three landslide-induced overtopping floods over a 530-kilometer reach of the Bartang and Panj Rivers below the Usoi dam. A one-dimensional flow model using a Riemann numerical solution technique was selected for the study. A constant 50-meter wide rectangular channel, which represented the mean channel width, was used for the entire reach. A roughness coefficient of 0.038, appropriate for steep mountainous streams, also was used for the entire reach. For the 87 million cubic meter volume overtopping flood scenario, the peak flows were approximately 1,100, 800, and 550 cubic meters per second at locations 50, 100, and 150 kilometers downstream of the dam, respectively. The model was also used to simulate the less likely scenario of an instantaneous dam breach and draining of the total volume of the lake. Simulated peak flows were approximately 64,000, 52,000, 40,000, and 20,000 cubic meters per second at locations 50, 100, 150, and 530 kilometers downstream of the Usoi dam.

Water-Resources Investigations Report↗

Flood hazard assessment of the Hoh River at Olympic National Park ranger station, Washington

Federal regulations require buildings and public facilities on Federal land to be located beyond or protected from inundation by a 100-year flood. Flood elevations, velocities and boundaries were determined for the occurrence of a 100-year flood through a reach, approximately 1-mi-long, of the Hoh River at the ranger station complex in Olympic National Park. Flood elevations, estimated by step-backwater analysis of the 100-year flood discharge through 14 channel and flood-plain cross sections of the Hoh River, indicate that the extent of flooding in the vicinity of buildings or public facilities at the ranger station complex is likely to be limited mostly to two historic meander channels that lie partly within loop A of the public campground and that average flood depths of about 2 feet or less would be anticipated in these channels. Mean flow velocities at the cross sections, corresponding to the passage of a 100-year flood, ranged from about 5 to over 11 ft/sec. Flooding in the vicinity of either the visitors center or the residential and maintenance areas is unlikely unless the small earthen dam at the upstream end of Taft Creek were to fail. Debris flows with volumes on the order of 100 to 1,000 cu yards could be expected to occur in the small creeks that drain the steep valley wall north of the ranger station complex. Historic debris flows in these creeks have generally traveled no more than about 100 yards out onto the valley floor. The potential risk that future debris flows in these creeks might reach developed areas within the ranger station complex is considered to be small because most of the developed areas within the complex are situated more than 100 yards from the base of the valley wall. Landslides or rock avalanches originating from the north valley wall with volumes potentially much larger than those for debris flows could have a significant impact on the ranger station complex. The probability that such landslides or avalanches may occur is unknown. Inspection of aerial photographs of the Hoh River valley revealed the apparent presence, along the ridge crest of the north valley wall, of ridge-top depressions--geologic features that are sometimes associated with the onset of deep-seated slope failures. However, evaluation of the potential landslide hazard associated with these depressions would require an onsite examination of the area by trained personnel. Such an effort was outside the scope of this study.

Washington↗

Investigating geomorphic change using a structure from motion elevation model created from historical aerial imagery: A case study in northern Lake Michigan, USA

South Manitou Island, part of Sleeping Bear Dunes National Lakeshore in northern Lake Michigan, is a post-glacial lacustrine landscape with substantial geomorphic changes including landslides, shoreline and bluff retreat, and sand dune movement. These changes involve interrelated processes, and are influenced to different extents by lake level, climate change, and land use patterns, among other factors. The utility of DEM of Difference (DoD) and other terrain analyses were investigated as a means of understanding interrelated geomorphologic changes and processes across multiple decades and at multiple scales. A 1m DEM was developed from 1955 historical aerial imagery using Structure from Motion Multi-View Stereo (SfM-MVS) and compared to a 2016 lidar-based DEM to quantify change. Landslides, shoreline erosion, bluff retreat, and sand dune movement were investigated throughout South Manitou Island. While the DoD indicates net loss or gain, interpretation of change must take into consideration the SfM-MVS source of the historical DEM. In the case of landslides, where additional understanding may be gleaned through review of the timing of lake high- and lowstands together with DoD values. Landscape-scale findings quantified cumulative feedbacks between interrelated processes. These findings could be upscaled to assess changes across the entire park, informing future change investigations and land management decisions.

Michigan↗

Laboratory simulation of earthquake-induced damage in lava dome rocks

Earthquakes can impart varying degrees of damage and permanent, inelastic strain on materials, potentially resulting in ruptures that may promote hazards such as landslides and other collapse events. However, the accumulation of damage in rocks under the frequency and amplitude of shaking experienced during earthquake events is rarely systematically measured due to technical limitations. Here, we characterize damage evolution during laboratory experiments on a suite of dacitic rocks from Unzen volcano, Japan, to help resolve accumulated damage and landslide susceptibility of lava domes during regional earthquake events. Damage was imparted during slow (time-dependent creep) and fast (stress-oscillation earthquake simulations) uniaxial loading in compression and tension. Damage evolution is approximated from strain during experiments; all samples accumulate strain during earthquake events, but microfracture-dominated samples tend to be more susceptible to damage than vesicle-dominated samples. The orientation of existing fabrics with respect to loading direction dictates the magnitude of strain accumulation under load oscillations. During each “earthquake” experiment of multiple dynamic stress-oscillations, samples accumulate inelastic strain. The strain imparted during each successive event is initially high and then reduces after 5-7 events, except when stressing results in failure. The strain rate during phases of intermittent stressing tends to be higher than prior to them. Understanding the accumulation of damage and the potential for brittle failure of rocks subjected to earthquakes can help define the origin and timing of certain landslides, rockfalls, lava dome collapses, and other failure events.

Mt. Unzen↗

Source processes for the probabilistic assessment of tsunami hazards

The importance of tsunami hazard assessment has increased in recent years as a result of catastrophic consequences from events such as the 2004 Indian Ocean and 2011 Japan tsunamis. In particular, probabilistic tsunami hazard assessment (PTHA) methods have been emphasized to include all possible ways a tsunami could be generated. Owing to the scarcity of tsunami observations, a computational approach is used to define the hazard. This approach includes all relevant sources that may cause a tsunami to impact a site and all quantifiable uncertainty. Although only earthquakes were initially considered for PTHA, recent efforts have also attempted to include landslide tsunami sources. Including these sources into PTHA is considerably more difficult because of a general lack of information on relating landslide area and volume to mean return period. The large variety of failure types and rheologies associated with submarine landslides translates to considerable uncertainty in determining the efficiency of tsunami generation. Resolution of these and several other outstanding problems are described that will further advance PTHA methodologies leading to a more accurate understanding of tsunami hazard.

Oceanography↗

Geologic hazards of the Wasatch Front, Utah

The results of recent and ongoing research into six significant geologic hazards of the Wasatch Front region will be summarized on this field trip, including: (1) surface fault rupture on the Salt Lake City segment of the Wasatch fault zone; (2) seismic site response in the Salt Lake Valley, including ground shaking and liquefaction; (3) liquefaction-induced landsliding at the Farmington Siding landslide complex; (4) lake flooding along the shores of Great Salt Lake; (5) debris-flow deposition on alluvial fans at the base of the Wasatch Range; and (6) landsliding in the Ogden area. The trip will provide an opportunity to discuss the scientific, engineering, and administrative aspects involved in geologic-hazard evaluation in this rapidly growing region.

Brigham Young University Geology Studies↗

Remote sensing of rainfall for debris-flow hazard assessment

Recent advances in remote sensing of rainfall provide more detailed temporal and spatial data on rainfall distribution. Four case studies of abundant debris flows over relatively small areas triggered during intense rainstorms are examined noting the potential for using remotely sensed rainfall data for landslide hazard analysis. Three examples with rainfall estimates from National Weather Service Doppler radar and one example with rainfall estimates from infrared imagery from a National Oceanic and Atmospheric Administration satellite are compared with ground-based measurements of rainfall and with landslide distribution. The advantages and limitations of using remote sensing of rainfall for landslide hazard analysis are discussed. ?? 2003 Millpress,.

Conference Paper↗

Morphology of chasma walls, Mars

The landforms developed on the walls of the Valles Marineris system of chasmas are of three major types, which are locally transitional. The most common type is composed of steep spurs and gullies. The dominant process in the formation or modification of this type appears to be the downslope movement of. material under the influence of gravity, resulting in the accumulation of extensive talus deposits. The type is morphologically similar to high, steep terrestrial scarps in desert or alpine environments. The second morphologic type consists of walls dissected by tributary canyons with characteristic V-shaped cross profiles and blunt canyon heads that locally contain lobate deposits. The tributary canyons may be relict features of the time, when the existence of running water was possible on the surface of Mars. The third morphologic type consists of landslide scars forming broad curved or straight recessed sections of chasma wall. This type is accompanied by landslide deposits that form hummocky floors at the base of the recessed sections. The landslides developed at the expense of other wall morphologies. Chains of rimless depressions and craters that parallel the main structural trends of the chasmas are best interpreted as collapse holes. The origin of the chasmas on Mars is conjectural and may have been structural (grabens), but, on the basis of morphologic studies of their walls, it is suggested that most of the present wall configuration is the result of erosional scarp retreat, where erosion follows preestablished structural planes of weakness.

Journal of Research of the U.S. Geological Survey↗

Tsunami flooding

Panel 5 focused on tsunami flooding with an emphasis on Probabilistic Tsunami Hazard Analysis (PTHA) as derived from its counterpart, Probabilistic Seismic Hazard Analysis (PSHA) that determines seismic ground-motion hazards. The Panel reviewed current practices in PTHA and determined the viability of extending the analysis to extreme design probabilities (i.e., 10 -4 to 10 -6 ). In addition to earthquake sources for tsunamis, PTHA for extreme events necessitates the inclusion of tsunamis generated by submarine landslides, and treatment of the large attendant uncertainty in source characterization and recurrence rates. Tsunamis can be caused by local and distant earthquakes, landslides, volcanism, and asteroid/meteorite impacts. Coastal flooding caused by storm surges and seiches is covered in Panel 7. Tsunamis directly tied to earthquakes, the similarities with (and path forward offered by) the PSHA approach for PTHA, and especially submarine landslide tsunamis were a particular focus of Panel 5.

Conference Paper↗

Recent damaging earthquakes in Japan, 2003-2008

During the last six years, from 2003-2008, Japan has been struck by three significant and damaging earthquakes: The most recent M6.6 Niigata Chuetsu Oki earthquake of July 16, 2007 off the coast of Kashiwazaki City, Japan; The M6.6 Niigata Chuetsu earthquake of October 23, 2004, located in Niigata Prefecture in the central Uonuma Hills; and the M8.0 Tokachi Oki Earthquake of September 26, 2003 effecting southeastern Hokkaido Prefecture. These earthquakes stand out among many in a very active period of seismicity in Japan. Within the upper 100 km of the crust during this period, Japan experienced 472 earthquakes of magnitude 6, or greater. Both Niigata events affected the south-central region of Tohoku Japan, and the Tokachi-Oki earthquake affected a broad region of the continental shelf and slope southeast of the Island of Hokkaido. This report is synthesized from the work of scores of Japanese and US researchers who led and participated in post-earthquake reconnaissance of these earthquakes: their noteworthy and valuable contributions are listed in an extended acknowledgements section at the end of the paper. During the Niigata Chuetsu Oki event of 2007, damage to the Kashiwazaki-Kariwa nuclear power plant, structures, infrastructure, and ground were primarily the product of two factors: (1) high intensity motions from this moderate-sized shallow event, and (2) soft, poor performing, or liquefiable soils in the coastal region of southwestern Niigata Prefecture. Structural and geotechnical damage along the slopes of dunes was ubiquitous in the Kashiwazaki-Kariwa region. The 2004 Niigata Chuetsu Earthquake was the most significant to affect Japan since the 1995 Kobe earthquake. Forty people were killed, almost 3,000 were injured, and many hundreds of landslides destroyed entire upland villages. Landslides were of all types; some dammed streams, temporarily creating lakes threatening to overtop their new embankments and cause flash floods and mudslides. The numerous landslides resulted, in part, from heavy rain associated with Typhoon Tokage. The earthquake forced more than 100,000 people into temporary shelters, and as many as 10,000 displaced from their upland homes for several years. Total damages was estimated by Japanese authorities at US$40 billion, making this the second most costly disaster in history, after the 1995 Kobe earth-quake. The 2003 Tokachi-Oki earthquake was the third event of magnitude 8.0+ to strike the southeastern portion of Hokkaido in the last 50 years. The event produced tsunami run-ups along the shoreline of southern Hokkaido that reached maximum heights of 4 meters. Accelerations recorded by seismic networks of Hokkaido indicated a high intensity motion region from Hiroo area to Kushiro City, with a PGA values in the range of 0.35 to 0.6g. Despite high acceleration levels, the observed ground failure, liquefaction, structural, port, and lifeline damages were remarkably light.

Conference Paper↗

Modelling root reinforcement in shallow forest soils

A hypothesis used to explain the relationship between timber harvesting and landslides is that tree roots add mechanical support to soil, thus increasing soil strength. Upon harvest, the tree roots decay which reduces soil strength and increases the risk of management -induced landslides. The technical literature does not adequately support this hypothesis. Soil strength values attributed to root reinforcement that are in the technical literature are such that forested sites can't fail and all high risk, harvested sites must fail. Both unstable forested sites and stable harvested sites exist, in abundance, in the real world thus, the literature does not adequately describe the real world. An analytical model was developed to calculate soil strength increase due to root reinforcement. Conceptually, the model is composed of a reinforcing element with high tensile strength, i.e. a conifer root, embedded in a material with little tensile strength, i.e. a soil. As the soil fails and deforms, the reinforcing element also deforms and stretches. The lateral deformation of the reinforcing element is treated analytically as a laterally loaded pile in a flexible foundation and the axial deformation is treated as an axially loaded pile. The governing differential equations are solved using finite-difference approximation techniques. The root reinforcement model was tested by comparing the final shape of steel and aluminum rods, parachute cord, wooden dowels, and pine roots in direct shear with predicted shapes from the output of the root reinforcement model. The comparisons were generally satisfactory, were best for parachute cord and wooden dowels, and were poorest for steel and aluminum rods. A parameter study was performed on the root reinforcement model which showed reinforced soil strength increased with increasing root diameter and soil depth. Output from the root reinforcement model showed a strain incompatibility between large and small diameter roots. The peak increase in soil strength attributed to roots was controlled by the small (<4mm) diameter root fraction. These results were used to calculate the effect of timber harvesting on a small, approximately 7.6 m 3 (10 yd 3 ), hypothetical landslide in a shallow, cohesionless, forest soil. The root reinforcement model predicted a post-harvest reduction in soil strength of 14 and 19 percent for a soil with and without 5 kPa (105 lbs/ft 2 ) of cohesion, respectively.

Thesis↗

Constraining parameter uncertainty in modeling debris-flow initiation during the September 2013 Colorado Front Range storm

The occurrence of debris flows during the September 2013 northern Colorado floods took the emergency management community by surprise. The September 2013 debris flows in the Colorado Front Range initiated from shallow landslides in colluvium. Most occurred on south- and east-facing slopes on the walls of steep canyons in crystalline rocks and on sedimentary hogbacks. Previous studies showed that most debris flows occurred in areas of high storm-total rainfall and that strength added by tree roots accounts for the low number of landslides in densely forested areas. Given the lack of rainfall thresholds for debris flow occurrence in northern Colorado, we want to parameterize a numerical model to assess potential for debris flows in advance of heavy rainfall. Natural Resources Conservation Service (NRCS) soil mapping of the area, supplemented by laboratory testing and field measurements, indicates that soil textures and hydraulic properties of landslide source materials vary considerably over the study area. As a step toward modeling storm response, available soil and geologic mapping have been interpreted to define zones of relatively homogeneous properties. A new, simplified modeling approach for evaluating model input parameters in the context of slope and depth of observed debris flow source areas and recorded debris-flow inducing rainfall helps narrow the range of possible parameters to those most likely to produce model results consistent with observed debris flow initiation. Initial results have narrowed the strength parameters to about one third of possible combinations of cohesion and internal friction angle and narrowed hydraulic conductivity to a range spanning slightly more than one order of magnitude.

Colorado↗

Debris-flow growth in Puerto Rico during Hurricane Maria: Preliminary results from analyses of pre- and post-event lidar data

On September 20, 2017, Hurricane Maria triggered widespread debris flows in Puerto Rico. We used field observations and pre- and post-Maria lidar to study the volumetric growth of long-travelled (>400 m) debris flows in four basins. We found overall growth rates that ranged from 0.7 to 30.4 m 3 per meter of channel length. We partitioned the rates into two growth mechanisms, aggregation of multiple landslides, or erosion and entrainment of channel sediment. In three basins, landslides accounted for more than 80% of the total debris-flow volumes. In one basin, entrainment accounted for 96% of the volume. These results indicate that forecasting volumes for regional debris-flow inundation modeling is more complicated than estimating the number and volume of contributing landslide source areas, although this task is difficult by itself. In this preliminary analysis, we did not find geologic, topographic, or morphometric variables that correlated with the growth observations. We suspect that the observed growth rates were heavily influenced by local variations in environmental conditions, including antecedent soil moisture conditions, duration and intensity of rainfall, and availability of channel material. Given these considerations, regional debris-flow inundation modeling may be best achieved by using a suite of scenarios that capture possible mechanisms of debris-flow growth.

Puerto Rico↗

A coccidioidomycosis outbreak following the Northridge, Calif, earthquake

Objective. —To describe a coccidioidomycosis outbreak in Ventura County following the January 1994 earthquake, centered in Northridge, Calif, and to identify factors that increased the risk for acquiring acute coccidioidomycosis infection. Design. —Epidemic investigation, population-based skin test survey, and casecontrol study. Setting. —Ventura County, California. Results. —In Ventura County, between January 24 and March 15, 1994, 203 outbreak-associated coccidioidomycosis cases, including 3 fatalities, were identified (attack rate [AR], 30 cases per 100 000 population). The majority of cases (56%) and the highest AR (114 per 100 000 population) occurred in the town of Simi Valley, a community located at the base of a mountain range that experienced numarous landslides associated with the earthquake. Disease onset for cases peaked 2 weeks after the earthquake. The AR was 2.8 times greater for persons 40 years of age and older than for younger persons (relative risk, 2.8; 95% confidence interval [CI], 2.1-3.7; P <.001). Environmental data indicated that large dust clouds, generated by landslides following the earthquake and strong aftershocks in the Santa Susana Mountains north of Simi Valley, were dispersed into nearby valleys by northeast winds. Simi Valley case-control study data indicated that physically being in a dust cloud (odds ratio, 3.0; 95% CI, 1.6-5.4; P<.001) and time spent in a dust cloud ( P <.001) significantly increased the risk for being diagnosed with acute coccidioidomycosis. Conclusions. —Both the location and timing of cases strongly suggest that the coccidioidomycosis outbreak in Ventura County was caused when arthrospores were spread in dust clouds generated by the earthquake. This is the first report of a coccidioidomycosis outbreak following an earthquake. Public and physician awareness, especially in endemic areas following similar dust cloud-generating events, may result in prevention and early recognition of acute coccidioidomycosis.

Journal of the American Medical Association↗

Rapid reservoir erosion, hyperconcentrated flow, and downstream deposition triggered by breaching of 38 m tall Condit Dam, White Salmon River, Washington

Condit Dam on the White Salmon River, Washington, a 38 m high dam impounding a large volume (1.8 million m 3 ) of fine-grained sediment (60% sand, 35% silt and clay, and 5% gravel), was rapidly breached in October 2011. This unique dam decommissioning produced dramatic upstream and downstream geomorphic responses in the hours and weeks following breaching. Blasting a 5 m wide hole into the base of the dam resulted in rapid reservoir drawdown, abruptly releasing ~1.6 million m 3 of reservoir water, exposing reservoir sediment to erosion, and triggering mass failures of the thickly accumulated reservoir sediment. Within 90&thinsp;min of breaching, the reservoir's water and ~10% of its sediment had evacuated. At a gauging station 2.3&thinsp;km downstream, flow increased briefly by 400&thinsp;m 3 &thinsp;s &minus;1 during passage of the initial pulse of released reservoir water, followed by a highly concentrated flow phase&mdash;up to 32% sediment by volume&mdash;as landslide-generated slurries from the reservoir moved downstream. This hyperconcentrated flow, analogous to those following volcanic eruptions or large landslides, draped the downstream river with predominantly fine sand. During the ensuing weeks, suspended-sediment concentration declined and sand and gravel bed load derived from continued reservoir erosion aggraded the channel by&thinsp;>1&thinsp;m at the gauging station, after which the river incised back to near its initial elevation at this site. Within 15&thinsp;weeks after breaching, over 1 million m 3 of suspended load is estimated to have passed the gauging station, consistent with estimates that&thinsp;>60% of the reservoir's sediment had eroded. This dam removal highlights the influence of interactions among reservoir erosion processes, sediment composition, and style of decommissioning on rate of reservoir erosion and consequent downstream behavior of released sediment.

Washington↗

Model simulations of flood and debris flow timing in steep catchments after wildfire

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most post-wildfire debris flows are generated from water runoff. The majority of existing debris-flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's $n$) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Water Resources Research↗