USGS Science⌕ Search

SEARCH · USGS Science

Results for “Earth Science”

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,189 records · Page 66Linked to original sources
USGS2021

A regional spatio-temporal analysis of large magnitude snow avalanches using tree rings

Snow avalanches affect transportation corridors and settlements worldwide. In many mountainous regions, robust records of avalanche frequency and magnitude are sparse or non-existent. However, dendrochronological methods can be used to fill this gap and infer historical avalanche patterns. In this study, we developed a tree-ring-based avalanche chronology for large magnitude avalanche events (size ≥ ∼ D 3 "> ≥ ∼ D 3 ) using dendrochronological techniques for a portion of the US northern Rocky Mountains. We used a strategic sampling design to examine avalanche activity through time and across nested spatial scales (i.e., from individual paths, four distinct subregions, and the region). We analyzed 673 samples in total from 647 suitable trees collected from 12 avalanche paths from which 2134 growth disturbances were identified over the years 1636 to 2017 CE. Using existing indexing approaches, we developed a regional avalanche activity index to discriminate avalanche events from noise in the tree-ring record. Large magnitude avalanches, common across the region, occurred in 30 individual years and exhibited a median return interval of approximately 3 years (mean = 5.21 years). The median large magnitude avalanche return interval (3–8 years) and the total number of avalanche years (12–18) varies throughout the four subregions, suggesting the important influence of local terrain and weather factors. We tested subsampling routines for regional representation, finding that sampling 8 random paths out of a total of 12 avalanche paths in the region captures up to 83 % of the regional chronology, whereas four paths capture only 43 % to 73 %. The greatest value probability of detection for any given path in our dataset is 40 %, suggesting that sampling a single path would capture no more than 40 % of the regional avalanche activity. Results emphasize the importance of sample size, scale, and spatial extent when attempting to derive a regional large magnitude avalanche event chronology from tree-ring records.

Montana↗
USGS2023

The influence of large woody debris on post-wildfire debris flow sediment storage

Debris flows transport large quantities of water and granular material, such as sediment and wood, and this mixture can have devastating impacts on life and infrastructure. The proportion of large woody debris (LWD) incorporated into debris flows can be enhanced in forested areas recently burned by wildfire, because wood recruitment into channels accelerates in burned forests. In this study, we explored how LWD transported in post-wildfire debris flows acts to retain sediment within small headwater channels. We found that the largest amount of debris flow sediment is stored in channel reaches with a ratio of LWD length to channel width between 0.25 and 1. Moreover, the largest debris flow deposits forced by LWD were found in low-angle channel reaches. We examined two models for debris flow volume estimation: (1) the current volume prediction model used in U.S. Geological Survey debris flow hazard assessments, and (2) a regional model developed to predict the sediment yield associated with debris-laden flows. We found that the regional model better matched the magnitude of the observed sediment at the terminal fan, suggesting the utility of regionally calibrated parameters for debris flow volume prediction. However, large wood created sediment storage upstream of the terminal fan, and this volume was of the same magnitude as the total volume at the terminal fans. Finally, we demonstrate a method for estimating debris flow velocity based on estimates of the critical velocity required to break wood, which can be used in future field studies to estimate minimum debris flow velocity values.

New Mexico↗
USGS2023

Landslide initiation thresholds in data-sparse regions: Application to landslide early warning criteria in Sitka, Alaska, USA

Probabilistic models to inform landslide early warning systems often rely on rainfall totals observed during past events with landslides. However, these models are generally developed for broad regions using large catalogs, with dozens, hundreds, or even thousands of landslide occurrences. This study evaluates strategies for training landslide forecasting models with a scanty record of landslide-triggering events, which is a typical limitation in remote, sparsely populated regions. We evaluate 136 statistical models trained on a precipitation dataset with five landslide-triggering precipitation events recorded near Sitka, Alaska, USA, as well as > 6000 d of non-triggering rainfall (2002–2020). We also conduct extensive statistical evaluation for three primary purposes: (1) to select the best-fitting models, (2) to evaluate performance of the preferred models, and (3) to select and evaluate warning thresholds. We use Akaike, Bayesian, and leave-one-out information criteria to compare the 136 models, which are trained on different cumulative precipitation variables at time intervals ranging from 1 h to 2 weeks, using both frequentist and Bayesian methods to estimate the daily probability and intensity of potential landslide occurrence (logistic regression and Poisson regression). We evaluate the best-fit models using leave-one-out validation as well as by testing a subset of the data. Despite this sparse landslide inventory, we find that probabilistic models can effectively distinguish days with landslides from days without slide activity. Our statistical analyses show that 3 h precipitation totals are the best predictor of elevated landslide hazard, and adding antecedent precipitation (days to weeks) did not improve model performance. This relatively short timescale of precipitation combined with the limited role of antecedent conditions likely reflects the rapid draining of porous colluvial soils on the very steep hillslopes around Sitka. Although frequentist and Bayesian inferences produce similar estimates of landslide hazard, they do have different implications for use and interpretation: frequentist models are familiar and easy to implement, but Bayesian models capture the rare-events problem more explicitly and allow for better understanding of parameter uncertainty given the available data. We use the resulting estimates of daily landslide probability to establish two decision boundaries that define three levels of warning. With these decision boundaries, the frequentist logistic regression model incorporates National Weather Service quantitative precipitation forecasts into a real-time landslide early warning “dashboard” system ( https://sitkalandslide.org/ , last access: 9 October 2023). This dashboard provides accessible and data-driven situational awareness for community members and emergency managers.

Alaska↗
USGS2024

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)↗
USGS2024

Characteristics of debris-flow-prone watersheds and debris-flow-triggering rainstorms following the Tadpole Fire, New Mexico, USA

Moderate- or high-severity fires promote increases in runoff and erosion, leading to a greater likelihood of extreme geomorphic responses, including debris flows. In the first several years following fire, the majority of debris flows initiate when runoff rapidly entrains sediment on steep slopes. From a hazard perspective, it is important to be able to anticipate when and where watershed responses will be dominated by debris flows rather than flood flows. Rainfall intensity averaged over a 15 min duration, I 15 , in particular, has been identified as a key predictor of debris flow likelihood. Developing effective warning systems and predictive models for post-fire debris flow hazards therefore relies on high-temporal resolution rainfall data at the time debris flows initiate. In this study, we documented the geomorphic response of a series of watersheds following a wildfire in western New Mexico, USA, with an emphasis on constraining debris flow timing within rainstorms to better characterize debris-flow-triggering rainfall intensities. We estimated temporal changes in soil hydraulic properties and ground cover in areas burned at different severities over >2 years to offer explanations for observed differences in spatial and temporal patterns in debris flow activity. We observed 16 debris flows, all of which initiated during the first several months following the fire. The average recurrence interval of the debris-flow-triggering I 15 is 1.3 years, which highlights the susceptibility of recently burned watersheds to runoff-generated debris flows in this region. All but one of the debris flows initiated in watersheds burned primarily at moderate or high soil burn severity. Since soil hydraulic properties appeared to be relatively resilient to burning, we attribute reduced debris flow activity at later times to decreases in the fraction of bare ground. Results provide additional constraints on the rainfall characteristics that promote post-fire debris flow initiation in a region where fire size and severity have been increasing.

New Mexico↗
USGS2024

Evaluation of debris-flow building damage forecasts

Reliable forecasts of building damage due to debris flows may provide situational awareness and guide land and emergency management decisions. Application of debris-flow runout models to generate such forecasts requires combining hazard intensity predictions with fragility functions that link hazard intensity with building damage. In this study, we evaluated the performance of building damage forecasts for the 9 January 2018 Montecito postfire debris-flow runout event, in which over 500 buildings were damaged. We constructed forecasts using either peak debris-flow depth or momentum flux as the hazard intensity measure and applied each approach using three debris-flow runout models (RAMMS, FLO-2D, and D-Claw). Generated forecasts were based on averaging multiple simulations that sampled a range of debris-flow volume and mobility, reflecting typical sources and magnitude of pre-event uncertainty. We found that only forecasts made with momentum flux and the D-Claw model could correctly predict the observed number of damaged buildings and the spatial patterns of building damage. However, the best forecast only predicted 50 % of the observed damaged buildings correctly and had coherent spatial patterns of incorrectly predicted building damage (i.e., false positives and false negatives). These results indicate that forecasts made at the building level reliably reflect the spatial pattern of damage but do not support interpretation at the individual building level. We found the event size strongly influences the number of damaged buildings and the spatial pattern of debris-flow depth and velocity. Consequently, future research on the link between precipitation and the volume of sediment mobilized may have the greatest effect on reducing uncertainty in building damage forecasts. Finally, because we found that both depth and velocity are needed to predict building damage, comparing debris-flow models against spatially distributed observations of building damage is a more stringent test for model fidelity than comparison against the extent of debris-flow runout.

Natural Hazards and Earth System Sciences↗
USGS2024

Assessing locations susceptible to shallow landslide initiation during prolonged intense rainfall in the Lares, Utuado, and Naranjito municipalities of Puerto Rico

Hurricane Maria induced about 70 000 landslides throughout Puerto Rico, USA, including thousands each in three municipalities situated in Puerto Rico's rugged Cordillera Central range. By combining a nonlinear soil-depth model, presumed wettest-case pore pressures, and quasi-three-dimensional (3D) slope-stability analysis, we developed a landslide susceptibility map that has very good performance and continuous susceptibility zones having smooth, buffered boundaries. Our landslide susceptibility map enables assessment of potential ground-failure locations and their use as landslide sources in a companion assessment of inundation and debris-flow runout. The quasi-3D factor of safety, F 3 , showed strong inverse correlation to landslide density (high density at low F 3 ). Area under the curve (AUC) of true positive rate (TPR) versus false positive rate (FPR) indicated success of F 3 in identifying head-scarp points (AUC = 0.84) and source-area polygons (0.85 ≤ AUC ≤ 0.88). The susceptibility zones enclose specific percentages of observed landslides. Thus, zone boundaries use successive F 3 levels for increasing TPR of landslide head-scarp points, with zones bounded by F 3 at TPR = 0.75, very high; F 3 at TPR = 0.90, high; and the remainder moderate to low. The very high susceptibility zone, with 118 landslides km −2 , covered 23 % of the three municipalities. The high zone (51 landslides km −2 ) covered another 10 %.

Puerto Rice↗
USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗
USGS2024

Rapid simulation of wave runup on morphologically diverse, reef-lined coasts with the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process model

Low-lying, tropical, coral-reef-lined coastlines are becoming increasingly vulnerable to wave-driven flooding due to population growth, coral reef degradation, and sea-level rise. Early-warning systems (EWSs) are needed to enable coastal authorities to issue timely alerts and coordinate preparedness and evacuation measures for their coastal communities. At longer timescales, risk management and adaptation planning require robust assessments of future flooding hazard considering uncertainties. However, due to diversity in reef morphologies and complex reef hydrodynamics compared to sandy shorelines, there have been no robust analytical solutions for wave runup to allow for the development of large-scale coastal wave-driven flooding EWSs and risk assessment frameworks for reef-lined coasts. To address the need for fast, robust predictions of runup that account for the natural variability in coral reef morphologies, we constructed the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process modeling system. We developed this meta-process model using a training dataset of hydrodynamics and wave runup computed by the XBeach Non-Hydrostatic process-based hydrodynamic model for 440 combinations of water level, wave height, and wave period with 195 representative reef profiles that encompass the natural diversity in real-world fringing coral reef systems. Through this innovation, BEWARE-2 can be applied in a larger range of coastal settings than meta-models that rely on a parametric description of the coral reef geometry. In the validation stage, the BEWARE-2 modeling system produced runup results that had a relative root mean square error of 13% and relative bias of 5% relative to runup simulated by XBeach Non-Hydrostatic for a large range of oceanographic forcing conditions and for diverse reef morphologies (root mean square error and bias 0.63 and 0.26 m, respectively, relative to mean simulated wave runup of 4.85 m). Incorporating parametric modifications in the modeling system to account for variations in reef roughness and beach slope allows for systematic errors (relative bias) in BEWARE-2 predictions to be reduced by a factor of 1.5–6.5 for relatively coarse or smooth reefs and mild or steep beach slopes. This prediction provided by the BEWARE-2 modeling system is faster by 4–5 orders of magnitude than the full, process-based hydrodynamic model and could therefore be integrated into large-scale EWSs for tropical, reef-lined coasts and used for large-scale flood risk assessments.

Natural Hazards and Earth System Sciences↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗
USGS2025

Topographic controls on landslide mobility: Modeling hurricane-induced landslide runout and debris-flow inundation in Puerto Rico

In 2017, Hurricane Maria triggered more than 70 000 landslides in Puerto Rico. After initiation, these predominantly shallow landslides were mobilized to varying extents – some landslides only traveled partway downslope, whereas others reached drainage channels and were mobilized into long-traveled debris flows that could severely impact roads and infrastructure. Thus, forecasting potential landslide runout and inundation zones is critical for estimating landslide and debris-flow hazards. Here we conduct an in-depth topographic analysis of landslide-affected areas from nine study areas and apply a linked modeling technique to estimate locations susceptible to varying degrees of landslide runout in the Lares, Utuado, and Naranjito municipalities. We find that the longest runout lengths are observed on high-relief escarpments, although highly mobile long-runout debris flows also occurred in lower-relief dissected uplands. These topographic differences indicate that landslides that are initiated under similar conditions and possess equal potential to be mobilized as debris flows may not travel the same distances or affect the same areal extent. Our modeling approach allows the local topography to automatically control the implementation of two runout methods: (1) H/L runout zones are assigned directly downslope of landslide source zones, and (2) debris-flow inundation zones are estimated in the presence of a channel network. Debris-flow volumes are calculated as a function of area-integrated growth factors, estimated as a function of the upstream areas susceptible to shallow landslides. Applying our empirical modeling scheme over an area of 560 km 2 , our results highlight the efficacy of our methods for the assessment of the potential for landslide runout and debris-flow inundation over diverse terrains with varied susceptibility.

Puerto Rico↗
USGS2025

Social sensing a volcanic eruption: Application to Kīlauea, 2018

Protecting lives and livelihoods during volcanic eruptions is the key challenge in volcanology, conducted primarily by volcano monitoring and emergency management organisations, but it is complicated by scarce knowledge of how communities respond in times of crisis. Social sensing is a rapidly developing practice that can be adapted for volcanology. Here we use social sensing of Twitter (currently known as X) posts to track changes in social action and reaction throughout the 2018 eruption of Kīlauea on the island of Hawai`i. The volume of relevant posts very rapidly increases in early May, coincident with the beginning of the eruption; automated sentiment analysis shows a simultaneous shift towards more negative emotions being expressed in post text. Substantial negative trends in sentiment are evident in reaction to high-impact events, including the destruction of a popular residential area and injuries sustained by tourists viewing the eruption. Topics of local Twitter conversation reveal societal actions, including the sharing of hazard warnings, mitigation actions, and aid announcements. Temporal trends in societal actions reflect patterns in volcanic activity (e.g. the peak and waning of eruptive activity), civil protection actions (e.g. risk mitigation actions and the communication of official warnings), and socioeconomic pressures (e.g. the destruction of homes). Local tweets detailing eruption damage and disruption display a similar temporal trend to independent estimates of the number of buildings in contact with lava. We show how hazard and risk information is discussed and reacted to on Twitter, which helps inform our understanding of community response actions and aids situational awareness, and outline how our approach could be adapted for use in real time.

Hawaii↗
USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗
USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗
USGS2026

Changing drivers of regional large magnitude avalanche frequency throughout Colorado, USA

Large magnitude snow avalanches (destructive size ≥ D3) impact settlements, transportation corridors, and public safety worldwide. In Colorado, United States, avalanches have killed more people than any other natural hazard since 1950. In March 2019, a large magnitude avalanche cycle occurred throughout the entire mountainous portion of Colorado resulting in more than 1000 reported avalanches during a two-week period. Nearly 200 of these avalanches were size D4 or larger with at least three D5 avalanches. However, placing this 2019 large magnitude avalanche cycle in historic context requires data prior to the instrumental record. Here, we paired tree disturbance data from dendrochronology (1698 to 2020) with meteorological data from the modeled and instrumental record (1901 to 2020) to understand the frequency and climate drivers of large magnitude snow avalanche cycles. The extensive number of downed trees from the 2019 avalanche cycle allowed us to collect 1,188 cross-sections and cores from 1023 individual trees within 24 avalanche paths across the state. From these samples we identified 4135 avalanche-related growth disturbances. We employed a strategic nested sampling design to spatially aggregate avalanche frequency from individual avalanche paths, to counties, to three major sub-regions (i.e., north, central, and south), and across the entire region (i.e., state of Colorado). Over a period spanning more than three centuries (1698 to 2020), we identified 76 avalanche years within 24 individual avalanche paths. Large magnitude avalanche event frequency varied across paths and sub-regions with several notable region-wide avalanche cycles. Both tree-ring and historical written records highlighted 1899 as a year with widespread and large magnitude avalanche activity similar to the March 2019 avalanche cycle. Since the early-20th century (1900 to 2020) regional avalanche probability declined significantly in parallel with decreasing snowpack throughout Colorado. Similarly, dominant avalanche regimes shifted from large magnitude regional cycles driven by above average snowfall years over most of the record, to regional avalanche cycles occurring more commonly in average to low snow years since 1988. In recent decades, a lack of December precipitation and above average March precipitation characterized years with regional large magnitude avalanche activity. Even with declining snow water equivalent, truly extreme regional large magnitude avalanche cycles remain possible – as demonstrated by the 2019 cycle. This underscores that rare but high-impact events are not eliminated by long-term trends. Understanding the changing snow and weather drivers and subsequent behavior of large magnitude avalanche cycles across multiple spatial scales may improve avalanche forecasting and the products and mitigations strategies developed by structural engineers to mitigate avalanche danger. This can decrease the avalanche risk to the public and improve infrastructure design in avalanche terrain.

Colorado↗
USGS2026

The damability function: A probabilistic approach to regional landslide dam susceptibility analysis applied to the Oregon Coast Range, USA

Landslides can dam rivers and require rapid response to mitigate catastrophic outburst floods. Here, we present a workflow to map landslide dam formation susceptibility at a regional scale. We define a probabilistic function that combines river valley width and landslide volume to efficiently determine the likelihood of a landslide dam or “damability”. We combine damability values with landslide susceptibility to estimate landslide dam susceptibility. The valley width measurements are automated using a new elevation threshold-based algorithm. Landslide volume is represented as a statistical distribution from mapped landslides. We validate and apply our approach to the Oregon Coast Range, USA and find that 36 % of river stretches exceed a dam potential threshold; these are in river headwaters and steeper terrain, which in this case correlate with more resistant lithologies. We also estimate volumes of the potential dammed lakes and find that most rivers with high dam susceptibility are less likely to impound large lakes because they have low drainage areas. However, widespread susceptibility, and the potential impacts from exceptionally large landslides, suggest that this hazard should be considered in the Pacific Northwest. The damability function workflow can ingest new data and be applied more broadly to assess future landslide dam hazards.

Oregon↗
USGS2026

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W↗
USGS2026

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York↗