Geology topics

Faults & earthquakes

Explore USGS earthquake and fault publications. Search fault names, locations and seismic hazards; follow the original report for maps, methods and publication context.

Follow this topic with RSS

453 records · Page 15Linked to original sources

USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗

USGS2024

Status and performance of the ShakeAlert® earthquake early warning system: 2019-2023

The U.S. Geological Survey (USGS)‐operated ShakeAlert® system is the United States West Coast earthquake early warning system ( Given et al. , 2018 ). In this study we detail ShakeAlert’s performance during some of the largest events seen by the system thus far. Statewide public alerting using ShakeAlert messages was authorized in California in October 2019. Over the next few years, public alerts were expanded into Oregon and then into Washington ( U.S. Geological Survey, 2024 ). ShakeAlert source results are routinely compared to the USGS Comprehensive Catalog (ComCat; Guy et al. , 2015 ; U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which contains the earthquake location and magnitude determined using complete waveform data. M 4.5 and larger is the threshold used for public alerting and was deliberately set below the level where damage is likely to compensate for cases where the system underestimates the magnitude. Between 17 October 2019 and 1 September 2023, the ShakeAlert system created 95 events with maximum magnitude estimates of M ≥4.5, the public alerting threshold. 94 of the 95 events were due to real earthquakes. Seven were categorized “false” per ShakeAlert’s internal definition that there was no matching catalog event within 100 km and 30 s of origin time; however, all but one of these were real earthquakes that were poorly located, primarily because they were at the edges of the seismic network. Three detected events were labeled “missed” because they were very poorly located (>100 km location error). In addition, the system did not produce solutions for four ComCat events M ≥4.5 ( U.S. Geological Survey, Earthquake Hazards Program, 2017 ), which were all at the edge of the alerting and network boundaries. The ShakeAlert system has accurately detected the majority of earthquakes that have occurred within the operational region since completing the public rollout, and alerts from the system have been delivered to millions of cell phone users throughout the West Coast.

California, Oregon, Washington↗

USGS2024

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii↗

USGS2025

The GorDAS Distributed Acoustic Sensing experiment above the Cascadia locked zone and subducted Gorda Slab

The southernmost portion of the Cascadia Subduction zone in Northern California produces high rates of moderate and large earthquakes owing to subduction of the Gorda slab and deformation associated with the Mendocino Triple Junction. Distributed Acoustic Sensing (DAS) is rapidly advancing as a method for detecting earthquakes and imaging crustal structure. We have begun a long-term DAS monitoring experiment on buried telecom fiber in Arcata, California, with the goal of increasing the available recordings of moderate to large earthquakes as well as imaging seismogenic structures. We have recorded over a year's worth of data, including most aftershocks of the 2022 M w 6.4 Ferndale earthquake, though not the mainshock itself. The dataset includes numerous magnitude 3.5 and larger earthquakes including the 2023/01/01 M w 5.4 Rio Dell earthquake. Here we present initial results comparing an earthquake detection algorithm, run in real-time on the processing unit of the interrogator system, with both the ShakeAlert earthquake early warning system as well as a post-processed earthquake catalog developed with deep-learning phase-picker algorithms. The rapid onboard processing of the detector demonstrates the potential utility of DAS-based edge computing for earthquake early warning. We also verify the quality of the strain waveforms both in terms of peak amplitudes and waveform similarity using about five months of nodal seismometer data. These instruments were deployed roughly every 300 m along the ~15km long cable and validate large variations in peak strain over short distances that are seen in the DAS data. All data from time windows surrounding both the local and teleseismic earthquakes are publicly available, which will improve our understanding of both the performance of DAS systems in moderate earthquakes and earthquake hazards associated with the Gorda subduction zone.

California↗

USGS2025

National population exposure and evacuation potential in the United States to earthquake-generated tsunami threats

Previous efforts to characterize tsunami threats to people have focused primarily on individual scenarios in specific areas but have not recognized multiple scenarios across an entire country. This study addresses this gap by quantifying population exposure and evacuation potential in the United States to 102 earthquake-related, tsunami-hazard zones, including 92 local scenarios, 8 distant scenarios, and 2 probabilistic products. Geospatial path-distance modeling quantified evacuation potential and the influence of departure delays. We focused on residents to support other national, multi-hazard risk analyses. Millions of residents are in distant-tsunami zones, and hundreds of thousands of residents are in local-tsunami zones. In 41 scenarios, there is at least one resident that may have insufficient time to evacuate before wave arrival. Tens of thousands of residents may have insufficient time to evacuate from local tsunamis that impact the U.S. Pacific Northwest or Puerto Rican coastlines. The largest improvements in evacuation potential may come from reducing departure delays in some areas but may involve vertical-evacuation structures or changing land use in other areas.

International Journal of Disaster Risk Reduction↗

USGS1976

Digital telemetry system network installation for the National Earthquake Information Service

A Digital Telemetry System Network, developed by the Albuquerque Seismological Laboratory, was installed for the USGS National Earthquake Information Service (NEIS) during January 1975. NEIS is located in Golden, Colorado. The initial Digital Telemetry System Network consists of four independent digital telemetry systems that provide wide area coverage over most of the United States. The four systems were installed between: Albuquerque, New Mexico to NEIS; St. Louis, Missouri to NEIS; Salt Lake City, Utah to NEIS; and Menlo Park, California to NEIS. Each digital telemetry system will transmit fourteen channels of high accuracy seismic data in a digital format over a standard C-1 grade telephone line. Digitized data words will have a resolution of eleven digital bits plus one sign bit. Each digital telemetry receiver, located at NEIS, will have the capacity to provide fourteen channels of seismic data in a digital format suitable for either on-line digital computer analysis or magnetic tape recording. In addition the data can be presented in an analog format suitable for recording on helicorders or develocorders. This network is the first large scale digital telemetry system used by the U.S. Geological Survey for the transmission of real time seismic data.

Open-File Report↗

USGS2026

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places. We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

California↗

USGS2024

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W–trending Central China Orogenic Belt has been overprinted by the N-S–trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau's weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W–trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene–Oligocene) and reworking of the deep lithosphere (Stage II, middle–late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

GSA Bulletin↗

USGS2024

Debris avalanches in the northern California Coast Range triggered by plate boundary earthquakes

Determining the timing and cause for ancient hillslope failures proves difficult in the western United States, yet critical as it ties directly into groundmotion estimates for hazardous events. This knowledge gap is important to confront as hillslope failures are candidates to be triggered by earthquakes along active plate boundaries. We identify two prehistoric, i.e., preinstrumental history, debris avalanches (3 10^6 and 6 10^6 m^3) in the Coast Ranges of northern California. These debris avalanches are well suited for studying the timing (to the exact year) and trigger as the densely forested landscape enables effective use of dendrochronology and high-resolution radiocarbon to pinpoint time of failure, and the legacy of landslide studies in the region provides context for evaluating a seismic trigger as the most likely triggering mechanism. Neither debris avalanche was triggered by the largest precipitation event of the instrumental record in 1964 CE, and the two debris avalanche sites are physiographically suited to accommodate topographic amplification of seismic shaking. Through a suite of geologic and dendrochronologic evidence, we establish the likely time of failure of the two debris avalanches as 1906 CE and 933 CE. The failure of the younger landslide coincides with the 1906 CE San Francisco earthquake on the Pacific/North American plate boundary, while the older landslide dates to 933 CE falling within the broad age range (850-966 CE) of a likely megathrust earthquake on the Cascadia subduction zone as recorded in coastal marshes in Humboldt Bay, California. The precise age on the 933 CE debris avalanche could trim the uncertainty on a Cascadia subduction zone event from 116 years to one year. Utilizing debris avalanche records from sites suitable to record seismic shaking improves understanding of plate-boundary earthquake timing and extent of shaking.

Bulletin of the Seismological Society of America↗

USGS2025

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters↗

USGS2025

Assessing earthquake risks to lifeline infrastructure systems in the United States

The security and economic stability of the United States rely heavily on robust lifeline infrastructure systems and yet the risks to such systems are seldom quantified at the national scale. For example, while earthquake risks to buildings in the United States have been investigated at the national scale regularly, such risks to gas pipelines have rarely been investigated nationally. In this paper, we use examples from two critical infrastructure sectors to demonstrate (1) the nature of earthquake risks to lifeline infrastructure systems, (2) complexities involved in regional seismic risk assessments, and (3) how such risks change with time. We found that bridge risks can be underestimated by at least 64 % when viewed from repair costs instead of traffic demands and that regional risks can be underestimated by 19 % when spatial correlations of ground motion are ignored. Further, exceedance of traffic demand can be 50 times more likely to occur when viewed at the regional scale than when viewed at an individual bridge. Similarly, exceedance of repairs can be 180 times more likely to occur when viewed at the pipeline network level than at a segment-specific level. Finally, sensitivity analyses with the 2018 and 2023 USGS National Seismic Hazard Models indicate an increase in bridge risk of at least 24 % and an increase in exposed gas pipeline mileage of 43 %. The evolution of risks, complexities involved in assessments, and limited resources jointly underscore the need for more routine updates to nationwide seismic risk assessments of lifeline systems in the United States.

International Journal of Critical Infrastructure P↗

USGS2025

Remote single-station seismic monitoring of the July–October 2022 earthquake swarm at Tau volcano, American Samoa

From July to October 2022, a non-eruptive volcanic earthquake swarm occurred within ~15 km of Taʻū Island, located in eastern American Samoa. Felt reports from local residents were the only available information about the swarm when it started, as American Samoa lacked a seismic monitoring network. We developed a consistent single-station catalog for the entire swarm, using seismic data from the nearest station IU.AFI, ~250 km away. We applied the EQTransformer deep-learning model (Mousavi et al., 2020), automatically picking Pn and Sn arrivals on IU.AFI continuous data. We retained only events with Sn-Pn times of 22.5–25 seconds, consistent with the expected locations based on felt reports, then detected smaller swarm events with subsequent template-matching. This single-station catalog characterized the swarm’s onset and escalation to peak activity before a multi-agency field response team installed a local seismic network in mid-August 2022. This permanent seismic network captured the swarm’s decline. EQTransformer identified short S-P times on the first two locally deployed seismometers, both Raspberry Shake sensors, to constrain the swarm’s distance from Taʻū Island. Modern seismological processing methods, combined with basic observations such as felt reports, can quickly contribute useful information during an earthquake response in a poorly monitored region.

American Somoa↗

USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗

USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗

USGS2025

Estimated average annualized losses from potential building damage and fatalities due to earthquake-generated tsunamis in the United States

Earthquake-generated tsunamis represent substantial economic threats to states and territories in the United States (U.S.), but we are unaware of any effort to quantify potential impacts at the national level. This gap is partially due to the lack of nationally consistent data on tsunamigenic sources and associated return periods. This study addresses this issue and provides estimates of average annualized losses (AAL) for potential residential fatalities and capital stock losses associated with building damage (i.e., structural, non-structural, contents, and inventory damage) in the U.S. by curating tsunami-hazard information based on deterministic scenarios and probabilistic approaches, calculating potential losses, and estimating return periods where necessary. This assessment was done for the U.S. West Coast, Alaska, Hawaii, U.S. Pacific Territories, and U.S. Atlantic Territories. We estimate that earthquake-generated tsunamis that could affect these states and territories collectively represent $1 billion in potential AAL with 79 % of losses due to residential fatalities and 21 % of losses due to capital stock losses from building damage. We identify AAL variations based on county and county equivalents, states and territories, geographic regions, return periods, and departure-delay assumptions for evacuating residents. Results include high AAL values for potential fatalities in Puerto Rico and the U.S. Pacific Northwest region, high AAL values for potential building-related damage in Hawaii and California, and high building- and population-loss ratios for county equivalents in Alaska and U.S. territories.

International Journal of Disaster Risk Reduction↗

USGS2025

Preliminary depth to basement modeling at Salton Sea, California

The San Andreas Fault – Imperial Fault (SAF-IF) transtensional step-over zone along the southern margin of the Salton Sea hosts substantial geothermal production and lithium brine resources. Recent volcanism at the Salton Buttes and active seismicity along the SAFIF fault system highlight active tectonic and magmatic processes that pose natural hazards and may impact energy and mineral production. Characterizing the subsurface architecture and extent of concealed alteration associated with this tectono-magmatic system enhances understanding of these active processes, associated hazards, and resources. We have compiled a gravity database, consisting of new and re-processed existing data, from which we have constructed a new isostatic residual gravity anomaly map of the Salton trough. We have used this new gravity dataset together with a compilation of publicly available borehole data to develop new depth to basement inversion models for the region. These depth to basement models help to constrain basin geometries, inform alteration mapping, and reveal variations in basement rocks. Due to the concealed nature of the complex tectonic framework at the Salton trough, it is necessary to utilize geophysical methods for subsurface characterization. These new depth to basement models are a first step toward constructing 2D and 3D geophysical and geologic models of the Imperial Valley and Salton Sea geothermal area. This analysis complements other geophysical initiatives, including magnetotelluric (MT) modeling (Tokmakoff et al., 2024), magnetic mapping (Glen and Earney, 2023, 2024) and potential field modeling, and seismic studies focused on hazard and resource investigations in the Imperial Valley.

California↗

USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗

USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗