FIND YOUR NEXT DISCOVERY

Results for “Landslides”

Original records, connected by a shared subject.

SourceResearch connectionFrom yearThrough year

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

62 records · Page 3Linked to original sources

USGS2026

Color-Shaded Relief Bathymetry of Lake Chelan, Washington

In 2023, the U.S. Geological Survey (USGS) completed a 3-week field operation in Lake Chelan, Washington, the third deepest lake in the conterminous United States, to acquire sediment cores, high-resolution bathymetry data, and chirp sub-bottom profiles. The high-resolution bathymetry displayed in this map publication was collected using a Norbit Winghead i67-curved array multibeam system; raw data were cleaned and gridded (3 by 3 meters [m]) in the Teledyne CARIS software. In some locations where the slopes of Lake Chelan are near-vertical, multibeam returns were sparse, resulting in artifacts resembling holes that appear as grey dots within the color-shaded bathymetric relief. Although the 2023 bathymetric mapping of Lake Chelan resulted in similar morphometric characteristics to previous hydrographic campaigns (table 1; Kendra and Singleton, 1987; Pelletier and others, 1989), these new data provide a much more detailed understanding of the two physiographic basins that comprise Lake Chelan: the deeper (as deep as 445.8 m [fig. 1]; previously reported as 453 m) and narrower (1–2 kilometers [km] wide) Lucerne Basin to the north, and the shallower (as deep as 127.8 m [fig. 1]), and wider (2–3 km wide) Wapato Basin to the south. With this detailed imaging of Lake Chelan, we find pervasive morphological evidence of slope failure and mass-transport deposits (MTDs) in at least 89 locations within the lake and 63 subaqueous deltas (fig. 2). The deeper Lucerne Basin is dominated by steep slopes and numerous subaqueous deltas fed by steep fluvial catchments, whereas the Wapato Basin contains numerous landslide scarps, scars, and debris aprons along the more gently sloping lake margins, and abundant relict kettle structures from the Pleistocene deglaciation. Lake Chelan’s geomorphic, geophysical, and sediment data contain important records of glacial retreat and advance, regional volcanic eruptions, and earthquakes that generate ground motions severe enough to induce landslides recorded in the sediment as MTDs. The epicenter of the 1872 Washington State Earthquake (M 6.5–7.2; Bakun and others, 2002; Brocher and others, 2017; Sherrod and others, 2021) is along the Spencer Canyon Fault, near the southern end of Lake Chelan. Shaking from this earthquake was strong enough to affect local aquifers, induce liquefaction of unconsolidated sediments, produce large seiches (Brocher and others, 2018), and likely cause landslides within Lake Chelan.

Washington↗

USGS2025

Geology and geochemistry of Ni-Co laterite occurrences of northern California and southern Oregon

Ni-Co laterite occurrences in the Klamath Mountains of southern Oregon and northern California contain Ni, Co, and Sc values comparable to global deposits. However, they are thin and limited in extent, often occurring as landslides and slumps. Olivine in the parent material contains Ni and Co and enstatite contains Sc. Through pervasive weathering, these elements where concentrated in Fe- and Mn-oxides. Unlike the mined Nickel Mountian deposit near Riddle, Oregon, these occurrences have no hydrous-Mg silicate minerals. Although the region experienced pervasive weathering in the Eocene, preservation of the laterite profiles was limited in the Klamath Mountains due to accelerated uplift and erosion compared to the Coast Range Provence.

Conference Paper↗

USGS2025

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗

USGS2025

Preface to the focus section on earthquake-induced ground deformation

Ground failure due to earthquake-induced liquefaction, surface fault rupture, and landsliding has caused substantial damage in the past and continues to pose risks to critical infrastructure in the future. The consideration of ground failure is central to engineering analysis and design, ranging from site-specific engineering assessments of ground failure potential to broader regional-scale ground failure assessments of near-real-time earthquake impacts or longer-term portfolio risk. Uncertainties in all aspects of the problem complicate the prediction of earthquake-induced ground failure, including the identification of hazard-susceptible geomaterials, characterizing their spatial extent and physical properties, spatio-temporal variability in groundwater conditions, characterization of earthquake ground motions, assessing ground failure severity, and linking ground failure severity to infrastructure damage and loss.

Seismological Research Letters↗

USGS2026

Only very strong shaking can break a tree

A handful of credible accounts, both recent and historical, describe trees being broken during earthquakes. Apart from landslides and tsunamis, such accounts are uncommon, suggesting that trees break only when ground acceleration is extreme or when trees are especially vulnerable to breakage. Here, I consider the question: what ground acceleration is required to break a living tree? I present simplified first‐principles calculations, considering the properties of common trees, focusing on tall, thin trees that can be approximated as uniform vertical cantilevers. The results suggest that, possibly excepting tall trees with heavy canopies, standing trees can only be broken by shaking approaching or exceeding 1 g , as anticipated by the rarity of documented accounts of snapped trees. I briefly consider several well‐documented instances of tree damage during strong earthquakes, drawing inferences about local ground motions.

The Seismic Record↗

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↗

USGS2025

Colored shaded-relief bathymetry and acoustic backscatter of Lake Sammamish, Washington

Evidence of strong earthquakes (such as underwater landslides and associated deposits) may be recorded within the lacustrine sediments of Pacific Northwest lakes. The floor of Lake Sammamish, Wash., an approximately 11 kilometer (6.8 mile) long, 2 kilometer (1.2 mile) wide, and 35 meter (114.8 feet) deep lake located in a populated region just east of Seattle, was mapped by the U.S. Geological Survey in November of 2021 to search for evidence of past earthquakes. Mapping was conducted using a SWATHplus-M 234-kHz interferometric side-scan sonar system was pole-mounted on the U.S. Geological Survey research vessel Parke Snavely, and the system collected full-coverage bathymetric and acoustic backscatter data, which were processed to 2-meter spatial resolution. Two maps were created, a colored shaded-relief bathymetric map showing lake floor morphology (sheet 1), and an acoustic-backscatter map showing backscatter intensities (sheet 2). The results may then be utilized together to investigate past earthquake activity.

Washington↗

USGS2024

Neotectonic mapping of Puerto Rico

Puerto Rico is part of the Puerto Rico-Virgin Islands microplate, along the Caribbean–North American plate boundary between the Puerto Rico trench subduction zone and the Muertos Trough incipient subduction zone. Despite recent seismicity and geodetically constrained deformation of ~3 mm/yr of left-lateral shear across the island, Quaternary fault locations remain largely uncertain. Preservation of recent faulting in the landscape is masked by distributed faulting, high weathering rates of the tropical climate, steep topography, frequent landsliding, and extensive agriculture and urbanization along coastlines and river valleys. We present remote neotectonic mapping of historical aerial imagery from the 1930s-60s and <1-m lidar topography, integrated with field observations, to create an updated active fault map of Puerto Rico. We focus on faults that offset younger geomorphic surfaces, with ages estimated by geologic mapping, OSL and radiocarbon dating, and morphologic interpretations. We present new evidence for Quaternary activity on seven faults, including the South Lajas, Salinas, Punta Montalva, Great Southern Puerto Rico, Cerro Goden, Parguera, and San Marcos faults. We find that active faulting occurs preferentially along the southern and western coasts and does not always spatially coincide with preexisting bedrock faults, possibly reflecting partial reactivation of older faults in the modern strain field.

Puerto Rico↗

USGS2025

Late Pleistocene kinematics of the Great Southern Puerto Rico Fault Zone, Puerto Rico

Several onshore faults in southern Puerto Rico have recently been recognized as Quaternary active. However, the kinematics of these faults, particularly any lateral component, remain largely unconstrained. It is difficult to characterize low strain‐rate faults, partially due to extensive erosional and anthropogenic landscape modification, steep relief, and frequent landsliding, limiting the preservation of geomorphic features that could serve as recorders of fault motion. Here, we constrain the kinematics along sections of the Great Southern Puerto Rico Fault Zone (GSPRFZ) on the southern coastal plain of Puerto Rico. We integrate ∼1‐m‐resolution light detection and ranging (lidar)‐derived topography, historical air photos, and field mapping to identify a series of ∼50–1200‐m‐long fault scarps and lineaments that trend northwest–southeast and extend for ≥25 km across the southern coastal plain. Fault scarps are primarily south facing, cut across topography, and displace Quaternary deposits and landforms. We document multiple offset geomorphic markers, including channel thalwegs and interfluves formed in deposits previously mapped as Quaternary piedmont alluvial plain. We observe both vertical (south‐side‐down) and right‐lateral meter‐scale displacements, which indicate that the GSPRFZ accommodates right‐lateral oblique motion in the late Pleistocene, consistent with northeast motion of the Puerto Rico and the Virgin Islands microplate away from the Hispaniola block.

Puerto Rico↗

USGS2024

Participatory engagement to reduce communication gaps

Underserved communities, especially those in coastal areas in Puerto Rico, face significant threats from natural hazards such as hurricanes and rising sea levels. Limited funding hinders the investment in costly mitigation measures, increasing exposure to natural disasters. Providing coastal resources and data products through effective communication mechanisms is fundamental to improving the well-being of these underserved coastal communities. The overall objectives of the pilot effort to engage and connect with underserved coastal communities in Puerto Rico were the following: (1) compile a comprehensive database of the projects and resources relevant to natural hazards in Puerto Rico; (2) foster connections with Puerto Rican interested parties to better understand their priorities regarding coastal hazards and provide them with pertinent U.S. Geological Survey (USGS) resources; and (3) identify knowledge gaps to guide future USGS projects in Puerto Rico. Here we outline our participatory engagement framework and process, along with two specific resources developed with the information collected from this effort. These resources are available in English and Spanish and consist of user-friendly, non-technical information products. Among them are: (1) a website where users can learn about USGS research on landslides, hurricanes, earthquakes, water resources, coastal hazards, tsunamis, and ecosystem hazards and environmental contaminants, and (2) a geonarrative highlighting shoreline changes in Puerto Rico with sections on historical shoreline trends, hurricane impacts, and potential solutions that could help protect communities and mitigate coastal hazards. Continuing participatory engagement in future projects could enhance the accessibility and usability of natural hazards resources within the community.

San Juan↗

USGS2025

Using the D-Claw software package to model lahars in the Middle Fork Nooksack River drainage and beyond, Mount Baker, Washington

Lahars, or volcanic mudflows, are the most hazardous eruption-related phenomena that will affect communities living along rivers that originate on Mount Baker. In the past 15,000 years, the largest lahars from Mount Baker have affected the Middle Fork Nooksack River drainage and beyond. Here we use the physics-based D-Claw software package to model nine lahar scenarios that are initiated as water-saturated landslides between Sherman Crater and the Roman Wall on the Mount Baker edifice and flow down the Middle Fork Nooksack River. The scenarios range in volume from 1 to 260 million cubic meters and have an initial hydraulic permeability from 10 −12 to 10 −10 meters squared. Model output includes data such as flow depth, velocity, runout distance, area inundated, arrival time, and sediment concentration as well as information that allows scientists to calculate other important hydrologic characteristics such as lahar discharge. These data are important to officials who have the responsibility to plan for, or take mitigation measures against, future Mount Baker lahars. To check the validity of the D-Claw results, we compare the scenarios to known geologic information. We also compare D-Claw results with empirical models that have been used in the past to determine potential inundation areas, runout distances, and arrival times. These comparisons highlight similarities and differences between empirical and physics-based models. We also present D-Claw scenario-based animations to help scientists, officials, and lay people alike to visualize how future lahars could affect communities.

Washington↗

USGS2025

Using subducting plate motion to constrain Cascadia slab geometry and interface strength

Subduction zones are home to multiple geohazards driven by the evolution of the regional tectonics, including earthquakes, volcanic eruptions and landslides. Past evolution builds the present-day structure of the margin, while the present-day configuration of the system determines the state-of-stress in which individual hazardous events manifest. Regional simulations of subduction zones provide a tool to synthesize the tectonic history of a region and investigate how geologic features lead to variations in the state of stress across the subduction system. However, it is challenging to design regional models that provide a force-balance that is consistent with the large-scale motion of surrounding tectonic plates while also not over-constraining the solution. Here, we present new models for the Cascadia subduction zone that meet these criteria and demonstrate how the motion of the subducting Juan de Fuca plate can be used to determine the along-strike variations in the viscous (long-term) coupling across the plate boundary. All successful models require lower viscous coupling in the northern section of the trench compared to the central and southern sections. However, due to uncertainties in the geometry of the Cascadia slab, we find that there is a trade-off between along-strike variation in viscous coupling and slab shape. Better constraints on the slab shape, and/or use of other observations are needed to resolve this trade-off. The approach presented here provides a framework for further exploring how geologic features in the overriding plate and the properties of the plate boundary region affect the state-of-stress across this and other subduction zones.

Cascadia subduction zone↗

USGS2026

Interplay between tectonics and submarine mass transport deposits in Cortes Basin: New high-resolution geophysics in the Outer California Borderland

The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.

California↗

USGS2026

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California↗

USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗

USGS2026

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗

USGS2025

Cascading land surface hazards as a nexus in the Earth system

Earth’s surface is sculpted by numerous processes that move sediment, ranging from gradual and benign to abrupt and catastrophic. Although infrequent, high-magnitude sediment mobilization events can be hazardous to people and infrastructure, leaving topographic imprints on the landscape and remarkable narratives in the historical record. Hazardous events such as fires, storms, and earthquakes accelerate erosion and sediment transport, increasing landscape sensitivity to subsequent perturbations, thus forming a cascading hazard. Although the redistribution of sediment across Earth’s landscape can result in higher risks to vulnerable populations, cascading processes are commonly unaccounted for in hazard assessments. Cascading hazards can occur almost immediately after triggering events, such as coseismic landslides, or over months, years, or even decades after an initial perturbation, such as debris flows after wildfires or flooding in channels alluviated by volcanic debris. Sediment cascades span Earth’s surface, from mountaintops to river valleys, where erosion, deposition, and aggradation can lead to a myriad of hazardous processes, including decreased river conveyance capacity, which increases the likelihood of downstream flooding. An improved understanding of the magnitude, frequency, and persistence of cascading hazards is critical given the rapid changes in the frequency and severity of storms, fires, sea-level change, and cryospheric melting, as well as the expansion of high-population-density urban footprints in regions susceptible to solid Earth hazards. Understanding the full consequences and underlying physics of Earth’s cascading land surface hazards can help minimize future human and economic losses.

California↗

USGS2025

Site-specific, extended ShakeMaps for earthquake engineering applications

The U.S. Geological Survey (USGS) routinely produces ShakeMaps of shaking intensity across the globe. Due to practical constraints, the number of response spectral periods was limited to three standard periods (0.3, 1.0, and 3.0 sec). We have recently developed the tools that are necessary to expand this functionality to include 22 periods (matching the current U.S. National Seismic Hazard Model periods) as well as the orientation-independent components (e.g., “RotD50”). We refer to ShakeMap products that include these extensions as “extended ShakeMaps.” The added level of complexity motivated us to also develop a user-friendly tool called the “ShakeMap Sampling Tool” (SST) that gives all the estimated shaking metrics for a specific location (or list of locations). Additionally, we develop a web application where users can input locations of interest and view/download the SST results. We further familiarize users with the concept of “Composite ShakeMaps.” For earthquakes sequences such as a mainshock and larger foreshocks and aftershocks, this provides a map of the maximum value of each shaking metric, which is useful for overall loss estimates, the full extent of ground failure triggering potential, and a better portrayal of the repeated shaking levels at a given point for a series of earthquakes. Such a site-specific shaking history facilitates earthquake forensics at building or infrastructure sites for which damage may be of concern, as described in the Disproportionate Damage Earthquake trigger specified in the IEBC (2018, Section 405.2.2) and in developing ATC-145 guidelines (Guidelines for Post-Earthquake Assessment, Repair, and Retrofit of Buildings). The composite ShakeMap can be combined with the SST for a variety of earthquake-hazard applications, such as systematically inferring triggering shaking estimates at specific sites of geotechnical interest for landsliding, liquefaction, and lateral-spreading hazards.

Kahramanmaraş↗