Search USGS⌕ Search

SEARCH · Search USGS

Results for “Natural 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 361 records · Page 20Linked to original sources

Rapid movement of frozen debris-lobes: implications for permafrost degradation and slope instability in the south-central Brooks Range, Alaska

We present the results of a reconnaissance investigation of unusual debris mass-movement features on permafrost slopes that pose a potential infrastructure hazard in the south-central Brooks Range, Alaska. For the purpose of this paper, we describe these features as frozen debris-lobes. We focus on the characterisation of frozen debris-lobes as indicators of various movement processes using ground-based surveys, remote sensing, field and laboratory measurements, and time-lapse observations of frozen debris-lobe systems along the Dalton Highway. Currently, some frozen debris-lobes exceed 100 m in width, 20 m in height and 1000 m in length. Our results indicate that frozen debris-lobes have responded to climate change by becoming increasingly active during the last decades, resulting in rapid downslope movement. Movement indicators observed in the field include toppling trees, slumps and scarps, detachment slides, striation marks on frozen sediment slabs, recently buried trees and other vegetation, mudflows, and large cracks in the lobe surface. The type and diversity of observed indicators suggest that the lobes likely consist of a frozen debris core, are subject to creep, and seasonally unfrozen surface sediment is transported in warm seasons by creep, slumping, viscous flow, blockfall and leaching of fines, and in cold seasons by creep and sliding of frozen sediment slabs. Ground-based measurements on one frozen debris-lobe over three years (2008–2010) revealed average movement rates of approximately 1 cm day −1 , which is substantially larger than rates measured in historic aerial photography from the 1950s to 1980s. We discuss how climate change may further influence frozen debris-lobe dynamics, potentially accelerating their movement. We highlight the potential direct hazard that one of the studied frozen debris-lobes may pose in the coming years and decades to the nearby Trans Alaska Pipeline System and the Dalton Highway, the main artery for transportation between Interior Alaska and the North Slope.

Alaska↗

Variable population exposure and distributed travel speeds in least-cost tsunami evacuation modelling

Evacuation of the population from a tsunami hazard zone is vital to reduce life-loss due to inundation. Geospatial least-cost distance modelling provides one approach to assessing tsunami evacuation potential. Previous models have generally used two static exposure scenarios and fixed travel speeds to represent population movement. Some analyses have assumed immediate departure or a common evacuation departure time for all exposed population. Here, a method is proposed to incorporate time-variable exposure, distributed travel speeds, and uncertain evacuation departure time into an existing anisotropic least-cost path distance framework. The method is demonstrated for hypothetical local-source tsunami evacuation in Napier City, Hawke's Bay, New Zealand. There is significant diurnal variation in pedestrian evacuation potential at the suburb level, although the total number of people unable to evacuate is stable across all scenarios. Whilst some fixed travel speeds approximate a distributed speed approach, others may overestimate evacuation potential. The impact of evacuation departure time is a significant contributor to total evacuation time. This method improves least-cost modelling of evacuation dynamics for evacuation planning, casualty modelling, and development of emergency response training scenarios. However, it requires detailed exposure data, which may preclude its use in many situations.

Natural Hazards and Earth System Sciences↗

On a report that the 2012 M 6.0 earthquake in Italy was predicted after seeing an unusual cloud formation

Several recently published reports have suggested that semi-stationary linear-cloud formations might be causally precursory to earthquakes. We examine the report of Guangmeng and Jie (2013), who claim to have predicted the 2012 M 6.0 earthquake in the Po Valley of northern Italy after seeing a satellite photograph (a digital image) showing a linear-cloud formation over the eastern Apennine Mountains of central Italy. From inspection of 4 years of satellite images we find numerous examples of linear-cloud formations over Italy. A simple test shows no obvious statistical relationship between the occurrence of these cloud formations and earthquakes that occurred in and around Italy. All of the linear-cloud formations we have identified in satellite images, including that which Guangmeng and Jie (2013) claim to have used to predict the 2012 earthquake, appear to be orographic – formed by the interaction of moisture-laden wind flowing over mountains. Guangmeng and Jie (2013) have not clearly stated how linear-cloud formations can be used to predict the size, location, and time of an earthquake, and they have not published an account of all of their predictions (including any unsuccessful predictions). We are skeptical of the validity of the claim by Guangmeng and Jie (2013) that they have managed to predict any earthquakes.

Natural Hazards and Earth System Sciences↗

Discussion and Review of “Seismology of the Oso-Steelhead landslide” by Hibert, Stark, and Ekstrom

In this study, Hibert et al. use the seismic waves generated by the disastrous 22 March 2014 landslide near Oso, Washington to investigate the dynamics of the event by combining long-period source inversions and short period analysis. Overall the paper is concise and well written and takes an innovative approach by integrating both high- and low-frequency seismic data together in the interpretation. It also demonstrates the useful prospect of resolving a trajectory of the center of mass of the landslide from a few traces of noisy seismic data. However, the main point of discussion that I would like to address stems from differences between the findings by Hibert et al. and the findings of a recently published multidisciplinary paper on the Oso landslide (Iverson et al., 2015). I am a co-author on the paper and performed the seismic analysis, which was similar to that of Hibert et al., yet the interpretation of how the event unfolded differs significantly between the two papers. In my opinion, this stems largely from a difference in information, some of which Hibert et al. unfortunately did not have access to prior to submission because it was just published in Iverson et al. (2015). However, differences also seem to stem from technical details of the methods and data used. Typically, such details are only of interest to seismologists, but in this case they make a substantial difference in understanding how the event unfolded and how it became so disastrous. For that reason, they are worthy of discussion here. The crux of the difference lies in two points 1) Iverson et al. (2015) found a multi-stage initiation of the main event, while Hibert et al. did not and 2) Iverson et al. (2015) found the second high-frequency event to be from a significantly smaller debris fall whereas Hibert et al. estimate it to be a substantial mass. I detail these two points below, then follow with an additional comment regarding the trajectory estimation made by Hibert et al. and some minor comments.

Natural Hazards and Earth System Sciences↗

Spatial and seasonal responses of precipitation in the Ganges and Brahmaputra river basins to ENSO and Indian Ocean dipole modes: Implications for flooding and drought

We evaluated the spatial and seasonal responses of precipitation in the Ganges and Brahmaputra basins as modulated by the El Niño Southern Oscillation (ENSO) and Indian Ocean Dipole (IOD) modes using Global Precipitation Climatology Centre (GPCC) full data reanalysis of monthly global land-surface precipitation data from 1901 to 2010 with a spatial resolution of 0.5° × 0.5°. The GPCC monthly total precipitation climatology targeting the period 1951–2000 was used to compute gridded monthly anomalies for the entire time period. The gridded monthly anomalies were averaged for the years influenced by combinations of climate modes. Occurrences of El Niño alone significantly reduce (88% of the long-term average (LTA)) precipitation during the monsoon months in the western and southeastern Ganges Basin. In contrast, occurrences of La Niña and co-occurrences of La Niña and negative IOD events significantly enhance (110 and 109% of LTA in the Ganges and Brahmaputra Basin, respectively) precipitation across both basins. When El Niño co-occurs with positive IOD events, the impacts of El Niño on the basins' precipitation diminishes. When there is no active ENSO or IOD events (occurring in 41 out of 110 years), precipitation remains below average (95% of LTA) in the agriculturally intensive areas of Haryana, Uttar Pradesh, Rajasthan, Madhya Pradesh, and Western Nepal in the Ganges Basin, whereas precipitation remains average to above average (104% of LTA) across the Brahmaputra Basin. This pattern implies that a regular water deficit is likely, especially in the Ganges Basin, with implications for the agriculture sector due to its reliance on consistent rainfall for successful production. Historically, major droughts occurred during El Niño and co-occurrences of El Niño and positive IOD events, while major flooding occurred during La Niña and co-occurrences of La Niña and negative IOD events in the basins. This observational analysis will facilitate well-informed decision making in minimizing natural hazard risks and climate impacts on agriculture, and supports development of strategies ensuring optimized use of water resources in best management practice under a changing climate.

Brahmaputra Basin, Ganges Basin↗

Spatial distribution of water level impact to back-barrier bays

Water level in semi-enclosed bays, landward of barrier islands, is mainly driven by offshore sea level fluctuations that are modulated by bay geometry and bathymetry, causing spatial variability in the ensuing response (transfer). Local wind setup can have a secondary role that depends on wind speed, fetch, and relative orientation of the wind direction and the bay. Inlet geometry and bathymetry primarily regulate the magnitude of the transfer between open ocean and bay. Tides and short-period offshore oscillations are more damped in the bays than longer-lasting offshore fluctuations, such as storm surge and sea level rise. We compare observed and modeled water levels at stations in a mid-Atlantic bay (Barnegat Bay) with offshore water level proxies. Observed water levels in Barnegat Bay are compared and combined with model results from the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system to evaluate the spatial structure of the water level transfer. Analytical models based on the dimensional characteristics of the bay are used to combine the observed data and the numerical model results in a physically consistent approach. Model water level transfers match observed values at locations inside the Bay in the storm frequency band (transfers ranging from 70-100%) and tidal frequencies (10-55%). The contribution of frequency-dependent local setup caused by wind acting along the bay is also considered. The approach provides transfer estimates for locations inside the Bay where observations were not available resulting in a complete spatial characterization. The approach allows for the study of the Bay response to alternative forcing scenarios (landscape changes, future storms, and rising sea level). Detailed spatial estimates of water level transfer can inform decisions on inlet management and contribute to the assessment of current and future flooding hazard in back-barrier bays and along mainland shorelines.

Natural Hazards and Earth System Sciences↗

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↗

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↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗