SEARCH · Geology Reports

Results for “Journal of Geophysical Research--Earth Surface”

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

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2Linked to original sources

USGS2021

Identifying geomorphic process domains in the synthetic landscapes of West Virginia, USA

Human activities such as mining, agriculture, and urbanization have resulted in severe, large‐scale alteration to landform organization and associated geomorphic processes. The mountaintop mining (MTM) region of West Virginia, USA has experienced dramatic topographic alteration, by removing steep slopes and introducing plateau‐like areas at ridgelines and benches on valley fills. The resulting engineered landforms create synthetic landscapes, disconnected from previous geomorphic processes. Invoking the process domain concept, we compare differences in slope‐area relations, cumulative area distributions (CADs), elevation, slope, upslope accumulated area, and a slope*area product before and after mining to adjacent unmined sub‐catchments in five study basins. Differences in the slope‐area relation include a 42% slope reduction in low drainage areas, corresponding to hillslopes, unchanneled valleys, and debris flow dominated channels, which may fall below thresholds required for debris flow processes. The curved slope‐area relation that represents valley incision by debris flows is replaced by slope‐area relations that resemble basins where gullying and the stream power law dominate. Extremely high chemical weathering of unconsolidated valley fills materials may facilitate process domain shifts from debris flows to gullying and fluvial erosion. The characteristic power law scaling break in CADs that represents the headward limit of the channelized network is subdued in post‐mined sites and may reflect headward channelized network extension in mined basins. Slope‐area relations and CADs present a unique topographic signature of MTM activity, potentially providing an analytical approach to assess impacts on underlying geomorphic processes for other synthetic landscapes such as cities or large‐scale agricultural production.

West Virginia↗

USGS2021

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗

USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗

USGS2021

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

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

California↗

USGS2022

Multi-stage soil-hydraulic recovery and limited ravel accumulations following the 2017 Nuns and Tubbs wildfires in Northern California

Wildfire can impact soil-hydraulic properties by reducing saturated hydraulic conductivity and sorptivity, making recently burned landscapes prone to debris flows and flash floods. The post-fire hazard window can range from years to decades. In Northern California, where wildfire frequency is steadily increasing, the impact and soil-hydraulic recovery from wildfires is unknown. Following the October 2017 Nuns and Tubbs fires in the Northern Bay Area of California, we established 41 monitoring sites for repeat tension-disc infiltrometer measurements of field-saturated hydraulic conductivity ( K fs ) over 3.5 years. Our site arrays, which encompass grasslands, chaparral, and oak and conifer forests across a range in lithology, show a marked decrease in K fs following the wildfires and a swift partial recovery following the initial post-fire rainy season. Our time series reveals a complex path to soil-hydraulic recovery marked by distinct seasonal stages. Analysis of changing K fs , sorptivity, and infiltration model residuals collectively suggests that these stages are related to transitions between soil-hydraulic processes like structural soil sealing from rainsplash, thermal cracking of bare soil, and vegetation regrowth. While soil infiltration rates were strongly impacted by the 2017 fires, dry ravel estimates are an order of magnitude less for similar slopes than the 2009 Station fire in the San Gabriel mountains of Southern California, suggesting that limited ravel flux may insufficiently load channels for debris flows that initiate from within-channel failure. Our analysis suggests that burned landscapes in the Northern Bay Area of California may experience rapid soil-hydraulic recovery and limited pathways toward post-fire debris flow initiation.

California↗

USGS2022

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface↗

USGS2023

Time-lapse seafloor surveys reveal how turbidity currents and internal tides in Monterey Canyon interact with the seabed at centimeter-scale

Here we show how ultra-high resolution seabed mapping using new technology can help to understand processes that sculpt submarine canyons. Time-lapse seafloor surveys were conducted in the axis of Monterey Canyon, ∼50 km from the canyon head (∼1,840 m water depth) over an 18-month period. These surveys comprised 5-cm resolution multibeam bathymetry, 1-cm resolution lidar bathymetry, and 2-mm resolution stereophotographic imagery. Bathymetry data reveal centimeter-scale textures that would be undetectable by more traditional survey methods. Upward-looking Acoustic Doppler Current Profilers at the site recorded the flow character of internal tides and the passage of three turbidity currents, while sediment cores collected from the site record flow deposits. Combined with flow and core data, the bathymetry shows how turbidity currents and internal tides modify the seabed. The turbidity currents drape sediment across the site, infilling bedform troughs and smoothing erosional features carved by the internal tides (e.g., rippled scours). Turbidity currents with speeds of 0.9–3.3 m/s failed to cause notable bedform movement, which is surprising given that flows with similar speeds produced rapid bedform migration elsewhere, including the upper Monterey Canyon. The lack of migration may be related to the character of the underlying substrate or indicate that turbidity currents at the site lack dense, near-bed layers. The scale of scours produced by the internal tides (≤0.7 m/s) approaches the scale of features recorded in the ancient rock record. Thus, these results illustrate how the scale gap between seabed mapping technology and the rock record may eventually be bridged.

California↗

USGS2023

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗

USGS2022

The influence of drying on the aeolian transport of river-sourced sand

Transgression and regression of water levels (stages) have impacted the evolution of aeolian landforms and sedimentary deposits throughout geologic history. We studied this phenomenon over a five-day period of reduced flow on the Colorado River in Grand Canyon National Park, AZ, USA, in March 2021. These transient low flows exposed river-channel sand deposits to the air, causing progressive desiccation (drying) and thereby making these deposits susceptible to aeolian transport. We measured aeolian threshold friction velocities ( u *t ) for sand saltation and PM10 dust emissions, as well as other characteristics, on a subaerially exposed sandbar and downwind aeolian dunefield during each day of the low river flow. The sandbar transitioned from supply-limited to transport-limited aeolian sediment transport conditions during the regression in river water stage. A possible tipping point between the two transport conditions occurred approximately 48 hours after the drop in river flow. The empirically measured u *t decreased as the sandbar sediment dried with increased subaerial exposure time. Theoretical estimates and empirical measurements of u *t corresponded closely on the aeolian dunefield and on the sandbar when it was drier during the third and fourth day of the experiment. Eighty-seven percent of the variability in u *t was explained by empirical models that provide practical estimates of aeolian transport potential of subaerial river sediment deposits using monitoring data that are commonly available in this and other river systems. The work provides theoretical insight into the response of aeolian processes to sediment supply changes driven by periods of anthropogenic activity, drought, and climate change.

Arizona↗

USGS2023

Impacts of spontaneous waterfall development on bedrock river longitudinal profile morphology

River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self-form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1-D numerical model to investigate how self-formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km-scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady-state. Consistency between our model predictions and field observations of waterfall-rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

Journal of Geophysical Research - Earth Surface↗

USGS2023

Kinematic evolution of a large paraglacial landslide in the Barry Arm fjord of Alaska

Our warming climate is adversely affecting cryospheric landscapes via glacial retreat, permafrost degradation, and associated slope destabilization. In Prince William Sound, Alaska, the rapid retreat of Barry Glacier has destabilized the slopes flanking the glacier, resulting in numerous landslides. The largest of these landslides (∼500 Mm 3 in volume) is more than 2 km wide and has the potential to generate a tsunami that could affect nearby recreationists, marine traffic, infrastructure, natural and cultural resources, and the community of Whittier, located 60 km from the landslide. Here, we combine landslide structural and kinematic element mapping with data acquired from bi-yearly airborne lidar, multi-week satellite-based synthetic aperture radar (SAR), sub-hourly ground-based SAR, and seismic monitoring from 2020 to 2022 to characterize this landslide and examine its evolution. While some methods serve as a snapshot in time that is a culmination of events, others emphasize the ever-evolving nature of the landslide and associated hazards. Four major kinematic elements define the overall structure of the landslide, which vary in deformation type and rate, from creep (5 mm per day over several months) to episodic movement (2 m in 30 days) and landslide-wide to localized events. In some areas of the landslide, short-term deformation deviates from structures formed by cumulative movement, implying structural and kinematic evolution associated with glacier retreat. These insights are important for assessing landslide hazards and hazard evolution for large, slow-moving bedrock landslides in actively deglaciating environments.

Alaska↗

USGS2023

Tree-ring derived avalanche frequency and climate associations in a high-latitude, maritime climate

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Alaska↗

USGS2024

Predicting characteristic length scales of barrier island segmentation in microtidal environments

Segmented barrier islands can be found in regions with small tidal ranges. In contrast to tidally dominated barriers, where inlet dynamics are thought to control island length scales, the controls on barrier island length scales in wave-dominated environments have not been quantified. These microtidal barriers typically have a curved shoreline, suggesting the influence of wave-driven alongshore sediment transport. Microtidal barriers are also typically hydrodynamically isolated from one another, as weak tidal flows limit interactions between adjoining barriers. To better understand the controls on and scales of barrier segmentation in the relative absence of tides, here we develop a theoretical framework to estimate the alongshore length scales at which a barrier will either breach or heal following a disturbance in the barrier morphology. The non-dimensional framework compares the timescales of overwash (advective) and alongshore sediment transport (diffusive) processes along barrier island chains. We then apply this framework to modern barrier islands in the microtidal Gulf of Mexico using wave hindcast data and the lengths, widths, heights, and lagoon depths measured from remotely sensed geospatial data and topobathymetric data. We find that most of these barriers are currently longer than their critical length scale, often as a result of coastal restoration efforts. Our critical length scale analysis suggests that most of the Gulf of Mexico barriers are vulnerable to segmentation despite coastal restoration efforts intended to protect fisheries and the mainland coasts.

Florida, Louisiana, Mississippi↗

USGS2024

Vortex trapping of suspended sand grains over ripples

Coastal hydrodynamics and morphodynamics integrate the effects of small-scale fluid-sediment interactions; yet, these small-scale processes are not well understood. To investigate sediment trapping by turbulent coherent structures or vortices, the transport of coarse sand over ripples was analyzed in a small-oscillatory flow tunnel with phase-separated Particle Image and Tracking Velocimetry. Results from one of the first direct measurements of vortex-trapped sand grains under oscillatory flows are presented. The vortices mobilized sand grains along the ripple slopes just prior to flow reversal and transported the suspended sediment grains. During several flow cycles, some sand grains were temporarily trapped in the vortex, prescribing semi-circular trajectories off-center from the vortex core in quadrants of the vortex that were closest to the ripple slope, as illustrated by Nielsen (1992, https://doi.org/10.1142/1269 ). Comparisons of the horizontal sediment grain velocity with the horizontal fluid velocity yielded a linear relationship with a slope of 0.87. The vertical grain velocities also varied linearly with the vertical fluid velocity with a slope of approximately 1 and an offset of −0.08 m s −1 . The offset is close to the still water settling velocity for coarse sand grains, as hypothesized during vortex trapping. Additionally, estimates of the off-center distance, between the centers of the semi-circular sediment paths and vortex cores, compared well with the ratio of the settling velocity to the radian frequency of the vortex yielding a linear regression slope of 0.99. Improved understanding of vortex trapping effects on sediment dynamics may decrease uncertainty in model predictions of large-scale coastal hydrodynamics and sediment transport.

Journal of Geophysical Research: Earth Surface↗

USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗

USGS2024

Postfire sediment mobilization and its downstream implications across California, 1984 – 2021

Fire facilitates erosion through changes in vegetation and soil, with major postfire erosion commonly occurring even with moderate rainfall. As climate warms, the western United States (U.S.) is experiencing an intensifying fire regime and increasing frequency of extreme rain. We evaluated whether these hydroclimatic changes are evident in patterns of postfire erosion by modeling hillslope erosion following all wildfires larger than 100 km 2 in California from 1984 to 2021. Our results show that annual statewide postfire hillslope erosion has increased significantly over time. To supplement the hillslope erosion modeling, we compiled modeled and measured postfire debris-flow volumes. We find that, in northern California, more than 50% of fires triggering the top 20 values of sediment mass and sediment yield occurred in the most recent decade (between 2011 and 2021). In southern California, the postfire sediment budget was dominated by debris flows, which showed no temporal trend. Our analysis reveals that 57% of postfire sediment erosion statewide occurred upstream of reservoirs, indicating potential impacts to reservoir storage capacity and thus increased risk to water-resource security with ongoing climate change.

California↗

USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗

USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

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

Alaska↗