Geology topics

Sean Richard LaHusen

Publications and source records attributed to Sean Richard LaHusen.

8 recordsLinked to original sources

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↗

USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

Washington↗

USGS2025

Strong shaking from past Cascadia Subduction Zone earthquakes encoded in coastal landforms

Strong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large-magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty-four deep-seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

California, Oregon, Washington↗

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

Complex landslide patterns explained by local intra-unit variability of stratigraphy and structure: Case study in the Tyee Formation, Oregon, USA

Lithology and geologic structure are important controls on landslide susceptibility and are incorporated into many regional landslide hazard models. Typically, metrics for mapped geologic units are used as model input variables and a single set of values for material strength are assumed, regardless of spatial heterogeneities that may exist within a map unit. Here we describe how differences in bedding thickness, grain size, inferred uniaxial compressive strength, and bedding dip control the inherent susceptibility of slopes to deep-seated failure within a single mapped geologic unit - the Tyee Formation of Oregon, USA. The Tyee, which covers over 15,000 km2 and underlies much of the Oregon Coast Range, comprises gently folded alternating beds of sandstone and siltstone deposited as turbidites, forming a 2-km thick Eocene submarine fan which has been uplifted and exhumed through the Cenozoic. Deep-seated landslides are widespread in the Tyee, but form a complex spatial pattern such that landslide density ranges from 0 to 24% of the total landscape area. These slides are often extensive and sufficiently deep to reduce local hillslope gradients, resulting in a strong negative correlation between landslide density and mean local slope. Mean annual precipitation and predicted strong ground motions from Cascadia earthquake scenarios also fail to explain the spatial distribution of deep-seated landslides. Consequently, landslide stability models, which are strongly influenced by landscape slope, pore-water pressure, and seismic acceleration, yield landslide susceptibility maps which are broadly anti-correlated with mapped deep-seated landslide density. Through a multivariable linear regression model, we show that much of the variance in deep-seated landslide density can be explained by variability of intra-unit stratigraphic and structural characteristics, which we measure at 128 sites across two study areas totaling ∼3000 km2. Our results suggest bedding dip is only weakly correlated to landslide density, but strongly influences landslide failure style. Subtle increases in bedding dip, even in the gently folded Tyee Formation, result in a substantially higher likelihood of a landslide being cataclinal, or parallel to bedding. Overall, we find a slight majority of landslides fail within these cataclinal slopes, and that these landslides tend to be larger than non-cataclinal landslides. We also show that the lithological and structural properties that influence landslide susceptibility are distinct for these two populations of landslides. Our results demonstrate how localized, intra-unit, geologic variability can exert strong control on landslide susceptibility and failure style. This suggests that in some locations, landslide hazard models could be significantly improved by incorporating detailed, spatially variable, geologic properties rather than relying solely on generalized geologic map units.

Oregon↗

USGS2023

Evidence of Seattle Fault earthquakes from patterns of deep-seated landslides

Earthquake‐induced landslides can record information about the seismic shaking that generated them. In this study, we present new mapping, Light Detection and Ranging‐derived roughness dating, and analysis of over 1000 deep‐seated landslides from the Puget Lowlands of Washington, U.S.A., to probe the landscape for past Seattle fault earthquake information. With this new landslide inventory, we observe spatial and temporal evidence of landsliding related to the last major earthquake on the Seattle fault ∼1100 yr before present. We find spatial clusters of landslides that correlate with ground motions from recent 3D kinematic models of Seattle fault earthquakes. We also find temporal patterns in the landslide inventory that suggest earthquake‐driven increases in landsliding. We compare the spatial and temporal landslide data with scenario‐based ground motion models and find stronger evidence of the last major Seattle fault earthquake from this combined analysis than from spatial or temporal patterns alone. We also compare the landslide inventory with ground motions from different Seattle fault earthquake scenarios to determine the ground motion distributions that are most consistent with the landslide record. We find that earthquake scenarios that best match the clustering of ∼1100‐year‐old landslides produce the strongest shaking within a band that stretches from west to east across central Seattle as well as along the bluffs bordering the broader Puget Sound. Finally, we identify other landslide clusters (at 4.6–4.2 ka, 4.0–3.8 ka, 2.8–2.6 ka, and 2.2–2.0 ka) in the inventory which let us infer potential ground motions that may correspond to older Seattle fault earthquakes. Our method, which combines hindcasting of the surface response to the last major Seattle fault earthquake, using a roughness‐aged landslide inventory with forecasts of modeled ground shaking from 3D seismic scenarios, showcases a powerful new approach to gleaning paleoseismic information from landscapes.

Bulletin of the Seismological Society of America↗

USGS2022

Limits to coseismic landslides triggered by Cascadia Subduction Zone earthquakes

Landslides are a significant hazard and dominant feature throughout the landscape of the Pacific Northwest. However, the hazard and risk posed by coseismic landslides triggered by great Cascadia Subduction Zone (CSZ) earthquakes is highly uncertain due to a lack of local and global data. Despite a wealth of other geologic evidence for past earthquakes on the Cascadia Subduction Zone, no landslides have been definitively linked to such earthquakes, even in areas otherwise highly susceptible to failure. While shallow landslides may not leave a lasting topographical signature in the landscape, there are thousands of deep-seated landslides in Cascadia, and these deposits often persist for hundreds of years and multiple earthquake cycles. Synthesizing newly developed inventories of dated large deep-seated landslides in the Oregon Coast Range, we use statistical methods to estimate the proportion of these types of landslides that could have been triggered during past great Cascadia Subduction Zone earthquakes. Statistical analysis of high-precision dendrochronology ages of landslide-dammed lakes and surface roughness-dated bedrock landslides reveal Cascadia Subduction Zone earthquakes may have triggered 0–15 % of large deep-seated landslides in the Oregon Coast Range over multiple earthquake cycles. Our results refine estimates from previous studies and further suggest that coseismic triggering accounts for a small fraction of the total deep-seated bedrock landslides mapped in coastal Cascadia. However, if the real rate of coseismic landslide triggering during CSZ earthquakes is near our estimated upper bound for the 1700 CSZ earthquake, we estimate up to 2400 coseismic large deep-seated landslides could occur in the Oregon Coast Range in a single earthquake. These findings suggest Cascadia is consistent with global observations from other subduction zones and that coseismic landslides may still represent a serious geohazard in the region.

Oregon↗

USGS2022

Developing landslide chronologies using landslide-dammed lakes in the Oregon Coast Range

The Oregon Coast Range is a dynamic landscape that is continually shaped by shallow and deep-seated landslides that can have disastrous consequences to infrastructure and human lives. Searching for evidence of potentially coseismic mass wasting is incredibly difficult, particularly when historical observations are limited. Landslide-dammed lakes with submerged “ghost forests” in the Oregon Coast Range present the unique opportunity to establish landslide chronologies with subannual accuracy when dendrochronology is applied. This field guide will visit the unique landslide-dammed Klickitat Lake and explore a drowned ‘ghost forest’ to discuss methods used to establish a prehistoric landslide chronology in western Oregon, USA. After exploring the lake and exposing its geomorphic secrets, the guide will end with a stop on Marys Peak, a mafic volcanic intrusion composed of gabbroic dikes and pillow basalt that forms the highest point in the Oregon Coast Range. With the landscape of western Oregon laid out before us, we will discuss short- and long-term geomorphic evolution of the Oregon Coast Range and Willamette Valley.

Oregon↗