SEARCH · Geology Reports

Results for “Bulletin Seismological Society America”

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 613 records · Page 34Linked to original sources

USGS2022

The 8 April 1860 Jour de Pâques earthquake sequence in southern Haiti

The grave threat posed by the Enriquillo‐Plantain Garden fault zone (EPGFZ) and other fault systems on the Tiburon Peninsula in southern Haiti was highlighted by the catastrophic M 7.0 Léogâne earthquake on 12 January 2010 and again by the deadly M 7.2 Nippes earthquakes on 14 August 2021. Early Interferometric Synthetic Aperture Radar observations suggest the 2021 earthquake broke structures associated with this fault system farther west of the 2010 event, but the rupture zones of both events are separated by a ∼50 km gap. This sequence provided the impetus to reconsider a nineteenth century earthquake that may have occurred within this gap. Though previous studies identified a single moderately large event on 8 April 1860, original sources describe a complex and distributed seismic sequence to the west of Port‐au‐Prince. These provide evidence for an initial event to the west of Les Cayes, on the southern coast of the Tiburon Peninsula. This was followed on the morning of 8 April 1860 by a damaging earthquake near l’Anse‐à‐Veau along the northern coast of the peninsula, which was succeeded 14 hr later by a larger mainshock to the east. Although locations cannot be determined precisely from extant macroseismic data, our preferred scenario includes an intensity magnitude ( M I ) ∼ 6.2 "> ( M I ) ∼ 6.2 (MI)∼6.2 foreshock, followed by an M I ∼ 6.8 "> M I ∼ 6.8 MI∼6.8 mainshock wholly or partially on the EPGFZ to the east of the foreshock. Our results suggest that, although earthquakes tend to trigger subsequent earthquakes on immediately neighboring fault segments, sequences can include gaps due to strain release by poorly characterized earthquakes in the historic past. We also document 83 aftershocks to the 1860 earthquake sequence, including a previously unknown damaging shock, estimated M I ∼ 6.6 "> M I ∼ 6.6 MI∼6.6 ⁠ , on 10 April 1860 near Jérémie, farther west of l’Anse‐à‐Veau.

Tiburon Peninsula↗

USGS2022

Real-time earthquake detection and alerting behavior of PLUM ground-motion-based early warning in the United States

We examine the real‐time earthquake detection and alerting behavior of the Propagation of Local Undamped Motion (PLUM) earthquake early warning (EEW) algorithm and compare PLUM’s performance with the real‐time performance of the current source‐characterization‐based ShakeAlert System. In the United States (U.S.), PLUM uses a two‐station approach to detect earthquakes. Once a detection is confirmed, observed modified Mercalli intensity (MMI) distributions are forecast onto a regular grid, in which the preferred alert regions are grid cells with MMI 4.0+ forecasts. Although locations of dense station coverage allow PLUM to detect small ( M < 4.5) earthquakes typically not considered for EEW in the U.S., a PLUM detection on a small earthquake does not always generate an alert. This is because PLUM alerts are determined by current shaking distributions. If the MMI 4.0+ shaking subsides prior to detection confirmation by shaking at a second neighboring station, the prior MMI 4.0+ information will not be in the alert forecasts. Of the 432 M 3.0+ U.S. West Coast earthquakes in 2021, 33 produced ground motions large enough to be detected by PLUM. Twenty‐four generated MMI 4.0+ PLUM alerts, whereas ShakeAlert issued public EEW alerts for 13 of these earthquakes. We compare PLUM and ShakeAlert alert regions with ShakeMap and “Did You Feel It?” intensity distributions. Because PLUM alert regions surround stations observed to have strong ground motions (regardless of earthquake magnitude), PLUM alerts reliably include locations that experience significant shaking. This is not necessarily the case for ShakeAlert alert regions when there are large errors in magnitude or epicenter estimates. For two of the largest earthquakes in our real‐time dataset, the M 6.0 Antelope Valley and M 5.1 Petrolia earthquakes, the inclusion of PLUM would have improved real‐time ShakeAlert performance. Our results indicate that incorporation of PLUM into ShakeAlert will improve the robustness of the EEW system.

California, Oregon, Washington↗

USGS2022

Estimation of site terms in ground-motion models for California using horizontal-to-vertical spectral ratios from microtremor

The horizontal‐to‐vertical spectral ratios from microtremor (mHVSR) data obtained at 196 seismic stations in California are used to evaluate three alternative microtremor‐based proxies for site amplification for use in ground‐motion models (GMMs): the site fundamental period ( ⁠ f 0 "> f 0 ⁠ ), the period‐dependent amplitude of the mHVSR( T ), and the normalized amplitude of the mHVSR( T ). The alternative parameters are evaluated for the sites with and without measurements of V S 30 "> V S 30 ⁠ . If a V S 30 "> V S 30 measurement is not available for a site, then f 0 "> f 0 has the highest correlation with the site amplification for short periods ( T <1 s) and the normalized amplitude of the mHVSR( T ) has the highest correlation for long periods ( T ≥1 s). If a measurement of the V S 30 "> V S 30 is available for a site, then the normalized amplitude of the mHVSR( T ) has the highest correlation for the site amplification not explained by V S 30 "> V S 30 for all periods. For both cases, the correlations are strongest at the longer periods as mHVSR( T ) measurements excel at providing valuable information for sites with long‐period amplification due to the deeper velocity structure. In particular, for sites with a V S 30 "> V S 30 measurement, the normalized mHVSR( T ) amplitude provides more information about the long‐period site terms than the basin depth currently used in GMMs. Empirical models of the median and standard deviation of the site terms based on the normalized mHVSR( T ) curves are developed for the two cases. These models can be used directly in the ASK14 GMM to modify the median and aleatory standard deviation or they can be used to estimate the site‐specific site term in the context of a partially nonergodic GMM. Including the mHVSR( T ) measurement can have a significant effect on estimates of the ground motion at a site: the range 5%–95% on the observed HVSR( T ) values corresponds to factors of 0.6–1.6 for the median spectral acceleration for periods between 0.5 and 4 s.

California↗

USGS2022

Estimates of k0 and effects on ground motions in the San Francisco Bay area

Ground‐motion studies are a key component of seismic hazard analyses and often rely on information of the source, path, and site. Extensive research has been done on each of these parameters; however, site‐specific studies are of particular interest to seismic hazard studies, especially in the field of earthquake engineering, as near‐site conditions can have a significant impact on the resulting ground motion at a site. There has been much focus on the constraint of site parameters and their application to seismic hazard studies, especially in the development of ground‐motion models (GMMs). Kappa is an observational parameter describing the high‐frequency attenuation of spectra, and its site contribution ( ⁠ κ 0 "> κ 0 ) has shown to be a good predictor of high‐frequency ground motions; however, measurements are often limited. In this study, we develop a κ 0 "> κ 0 dataset for the San Francisco Bay area (SFBA) by estimating κ 0 "> κ 0 for 228 stations, and we produce a continuous regional map of κ 0 "> κ 0 ⁠ . We find κ 0 "> κ 0 to range between 0.003 and 0.072 s, with larger values concentrating on the east, north, and south sides of the bay, and lower values concentrating on the west side. We also evaluate the robustness of κ 0 "> κ 0 as a site parameter and find it to correlate with peak ground acceleration. These estimates of κ 0 "> κ 0 can add predictive power to GMMs, thus increasing the accuracy of predicted ground motion and improving the robustness of ground‐motion studies in the SFBA.

California↗

USGS2023

Comparison of nonergodic ground-motion components from CyberShake and NGA-West2 datasets in California

In this study, we compare the Southern California Earthquake Center CyberShake platform against the Next Generation Attenuation‐West2 empirical datasets. Because the CyberShake and empirical datasets cover very different magnitude ranges and site conditions, we develop ground‐motion models (GMMs) for CyberShake datasets to compare trends with empirical GMMs and decompose the residuals for further analysis. We apply mixed effects regression to four CyberShake datasets in southern, central, and northern California at 2, 3, 5, and 10 s periods, and compare the results with the empirical datasets using the same approach. CyberShake captures the total variability of ground motions in the empirical datasets but tends to predict larger median ground motions relative to the empirical GMMs. We then calculate and compare the repeatable source‐specific location, site, and path effects between CyberShake and empirical datasets. We find that the correlations of site effects between the CyberShake and empirical datasets are generally satisfactory, but the variability of site effects is slightly smaller for CyberShake datasets. There is no apparent correlation of source‐specific location effects between the CyberShake and empirical datasets. Comparison of path effects shows a wide range of correlation coefficients. Finally, we investigate the source of observed differences between the CyberShake and empirical datasets. We attribute the larger median ground‐motion levels in CyberShake to a combination of the homogeneous slip patterns of the earthquake ruptures, the low resolution of near‐surface materials in the velocity models, and strong reflections at high‐contrast boundaries in the velocity models. These factors also impact the correlations of site and path effects between the CyberShake and empirical datasets. Moreover, the leakage from location effects into site and path terms further weakens the correlations. In summary, we find that CyberShake could be improved, but it is still very useful to supplement empirical datasets for ground‐motion studies, especially to inform their nonergodic components.

California↗

USGS2022

The Geodetic Centroid (gCent) Catalog: Global earthquake monitoring with satellite imaging geodesy

Remote sensing geodetic observations (Interferometric Synthetic Aperture Radar [InSAR] and optical correlation [“pixel tracking”]) serve an increasingly diverse and important role in earthquake monitoring and response. This study introduces the Geodetic Centroid (gCent) catalog—an earthquake catalog derived solely from space‐based geodetic observations—and analysis of 74 earthquakes ( ⁠ M W 4.3–7.4) imaged from 1 August 2019 to 01 February 2022. For gCent, we use InSAR and optical correlation observations derived from the Sentinel‐1 satellites and various publicly available optical satellites to systematically image all global earthquakes M W 5.5 or larger and shallower than 25 km, M W 7.0 or larger at any depth, and other high‐impact earthquakes or seismic events of special interest. We invert surface displacements from successfully imaged earthquakes for the location, orientation, and dimensions of a single slipping fault patch that describes the centroid characteristics of the earthquake. These centroid models, in turn, are compiled into a catalog and used in U.S. Geological Survey/Advanced National Seismic System (ANSS) operational earthquake response products such as ShakeMaps and finite‐fault models. We provide a comparison of the gCent catalog to the ANSS Comprehensive Catalog and Global Centroid Moment Tensor (Global CMT) catalog to compare reported locations, depths, and magnitudes. We find that global earthquake catalogs not only generally provide reasonably comparable locations (within 10 km on average), but also they systematically overestimate depth that may have implications for earthquake shaking predictions based solely on earthquake origin information. Geodetic magnitudes are comparable to seismically inferred magnitudes, indicating that gCent models are unlikely to be systematically biased by the presence of postseismic deformation. We additionally highlight limitations of the gCent catalog induced by both the limitations of remote sensing imaging of earthquakes and our imposition of a simplified earthquake source description that does not include spatially distributed slip.

Bulletin of the Seismological Society of America↗

USGS2022

Rupture scenarios for the 3 June 1770 Haiti earthquake

The 2010 M 7.0 Haiti earthquake provided the impetus to reconsider historical earthquakes in Hispaniola ( Bakun et al. , 2012 ). That earthquake also shed new light on complex fault systems along Haiti’s southern peninsula ( Douilly et al. , 2013 ; Saint Fleur et al. , 2015 ). Recently, the 2021 M 7.2 Nippes earthquake ( Calais et al. , 2022 ; Douilly et al. , 2022 ), and a recent study reconsidering the 1860 sequence ( Martin et al. , 2022 ) further underscored the complexity of fault systems and large earthquake ruptures along the peninsula. Motivated by these studies and recent geological investigations ( Prentice et al. , 2010 ; Saint Fleur et al. , 2020 ), we reconsider the 3 June 1770 Haiti earthquake to explore the conventional assumption that it was the last major ( M ≥7.5) earthquake along the Enriquillo–Plantain Garden fault (EPGF). Accounts provide compelling evidence for substantial liquefaction in the Cul‐de‐Sac plain, one or more likely landslide‐driven tsunami in Gonaïves Bay, and extensive landsliding that created at least three documented landslide dams. We consider three end‐member rupture scenarios that are consistent with available constraints: two scenarios with M 7.7 and rupture lengths of 150–170 km, and one scenario with a ∼90 km rupture and M 7.5. Absent future work to identify and date paleoevents along the southern peninsula, none of these scenarios can be ruled out. Our preferred rupture model extends from the Miragoâne pull‐apart to near la Selle mountain, with a rupture length of 127 km, M 7.6, and a high stress drop. Rupture could have been on the EPGF or on an oblique thrust fault associated with overthrusting of the Massif de la Selle. The results do support the conclusion that the 1770 earthquake was the last major earthquake in southern Haiti, with a magnitude upward of M 7.5 and significantly more severe shaking in southern Haiti than during the 2010 earthquake.

Bulletin of the Seismological Society of America↗

USGS2023

A methodology to combine shaking and ground failure models for forecasting seismic damage to buried pipeline networks

How does an earthquake affect buried pipeline networks? It is well known that the seismic performance of buried pipelines depends on ground failures (GFs) as well as strong ground shaking (SGS), but it is unclear how the various types of earthquake hazards should be collectively combined, as existing methodologies tend to examine each of the earthquake hazards separately. In this article, we develop a probability‐based methodology to consistently combine SGS with three types of GF (surface faulting, liquefaction, and landslide) for forecasting seismic damage in buried pipeline networks from a given earthquake rupture scenario. Using a gas transmission pipeline example, we illustrate how the proposed methodology enables others (e.g., researchers, pipeline operators who manage distribution lines, and consultants) to modularly combine various models such as those for estimating probability of GF, permanent ground displacements, and pipeline fragility. Finally, we compare the proposed methodology against the Hazus methodology to explore implications from considering each hazard one at a time.

Bulletin of the Seismological Society of America↗

USGS2023

Applications of nonergodic site response models to ShakeAlert case studies in the Los Angeles area

In this study, we explore whether the Parker and Baltay (2022) site response models for the Los Angeles (LA) basin region can improve ground‐motion forecasts in the U.S. Geological Survey ShakeAlert earthquake early warning system (hereafter ShakeAlert). We implement the peak ground acceleration and peak ground velocity site response models of Parker and Baltay (2022) in ShakeAlert via the earthquake information to ground‐motion (hereafter eqinfo2GM) module, which predicts ground motions from the estimated earthquake parameters of magnitude, rupture length, and location. The nonergodic site response models for the greater LA area were developed using ground motions from 414 M 3–7.3 earthquakes in southern California. We test nonergodic ground‐motion forecasts for five earthquakes in the LA area: the 1994 M 6.7 Northridge earthquake, the 2008 M 5.4 Chino Hills earthquake, the 2019 M 7.1 Ridgecrest earthquake, the 2020 M 4.5 South El Monte earthquake, and a synthetic M 7.8 earthquake on the southern San Andreas fault from the ShakeOut scenario, which was the basis of a statewide emergency response exercise. From the test results, we find that with the nonergodic site response applied, ShakeAlert not only alerts larger areas but can also result in longer warning times in LA region. In addition, the modified Mercalli intensity (MMI) ground‐motion predictions generated by the ShakeAlert eqinfo2GM module are improved in accuracy when compared with the corresponding ShakeMap ground‐truth MMI when the nonergodic site response model is applied.

California↗

USGS2023

Seafloor observations eliminate a landslide as the source of the 1918 Puerto Rico Tsunami

The 11 October 1918 devastating tsunami in northwest Puerto Rico had been used as an example for earthquake‐induced landslide tsunami hazard. Three pieces of evidence pointed to a landslide as the origin of the tsunami: the discovery of a large submarine landslide scar from bathymetry data collected by shipboard high‐resolution multibeam sonar, reported breaks of submarine cable within the scar, and the fit of tsunami models to flooding observations. Newly processed seafloor imagery collected by remotely operated vehicle (ROV) show, however, pervasive Fe–Mn crust (patina) on the landslide walls and floor, indicating that the landslide scar is at least several hundred years old. C 14 "> 14 C dates of sediment covering the landslide floor verify this interpretation. Although we have not searched the region systematically for an alternative tsunami source, we propose a possible source—a two‐segment normal‐fault rupture along the eastern wall of Mona rift. The proposed fault location matches the published normal faults with steep bathymetry and is close to the International Seismological Center–Global Earthquake Model catalog locations of the 1918 mainshock and aftershocks. The ROV observations further show fresh vertical slickensides and rock exposure along the proposed fault trace. Hydrodynamic models from an M w "> M w 7.2 earthquake rupture along the eastern wall of the rift faithfully reproduce the reported tsunami amplitudes, polarities, and arrival times. Our analysis emphasizes the value of close‐up observations and physical samples to augment remote sensing data in natural hazard studies.

Puerto Rico↗

USGS2023

Combinatorial optimization of earthquake spatial distributions under minimum cumulative stress constraints

We determine optimal on‐fault earthquake spatial distributions using a combinatorial method that minimizes the long‐term cumulative stress resolved on the fault. An integer‐programming framework was previously developed to determine the optimal arrangement of a millennia‐scale earthquake sample that minimizes the misfit to a target slip rate determined from geodetic data. The resulting cumulative stress from just slip‐rate optimization, however, can greatly exceed fault strength estimates. Therefore, we add an objective function that minimizes cumulative stress and broad stress constraints to limit the solution space. We find that there is a trade‐off in the two objectives: minimizing the cumulative stress on a fault within fault strength limits concentrates earthquakes in specific areas of the fault and results in excursions from the target slip rate. Both slip‐rate and stress objectives can be combined in either a weighted or lexicographic (hierarchical) method. Using a combination of objectives, we demonstrate that a Gutenberg–Richter sample of earthquakes can be arranged on a constant slip‐rate finite fault with minimal stress and slip‐rate residuals. We apply this method to determine the optimal arrangement of earthquakes on the variable slip‐rate Nankai megathrust over 5000 yr. The sharp decrease in slip rate at the Tokai section of the fault results in surplus cumulative stress under all scenarios. Using stress optimization alone restricts this stress surplus to the northeast end of the fault at the expense of decreasing the slip rate away from the target slip rate at the southwest end of the fault. A combination of both slip‐rate and stress objectives provides an adequate fit to the data, although alternate model formulations for the fault are needed at the Tokai section to explain persistent excess cumulative stress. In general, incorporating stress objectives and constraints into the integer‐programming framework adds an important aspect of fault physics to the resulting earthquake rupture forecasts.

Bulletin of the Seismological Society of America↗

USGS2023

Incorporation of real-time earthquake magnitudes estimated via peak ground displacement scaling in the ShakeAlert Earthquake Early Warning system

The United States earthquake early warning (EEW) system, ShakeAlert®, currently employs two algorithms based on seismic data alone to characterize the earthquake source, reporting the weighted average of their magnitude estimates. Nonsaturating magnitude estimates derived in real time from Global Navigation Satellite System (GNSS) data using peak ground displacement (PGD) scaling relationships offer complementary information with the potential to improve EEW reliability for large earthquakes. We have adapted a method that estimates magnitude from PGD ( Crowell et al. , 2016 ) for possible production use by ShakeAlert. To evaluate the potential contribution of the modified algorithm, we installed it on the ShakeAlert development system for real‐time operation and for retrospective analyses using a suite of GNSS data that we compiled. Because of the colored noise structure of typical real‐time GNSS positions, observed PGD values drift over time periods relevant to EEW. To mitigate this effect, we implemented logic within the modified algorithm to control when it issues initial and updated PGD‐derived magnitude estimates ( ⁠ M PGD "> M PGD "> M PGD PGD ⁠ ), and to quantify M PGD "> M PGD PGD uncertainty for use in combining it with estimates from other ShakeAlert algorithms running in parallel. Our analysis suggests that, with these strategies, spuriously large M PGD "> M PGD PGD will seldom be incorporated in ShakeAlert’s magnitude estimate. Retrospective analysis of data from moderate‐to‐great earthquakes demonstrates that the modified algorithm can contribute to better magnitude estimates for M w > 7.0 "> M w > 7.0 w>7.0 events. GNSS station distribution throughout the ShakeAlert region limits how soon the modified algorithm can begin estimating magnitude in some locations. Furthermore, both the station density and the GNSS noise levels limit the minimum magnitude for which the modified algorithm is likely to contribute to the weighted average. This might be addressed by alternative GNSS processing strategies that reduce noise.

Bulletin of the Seismological Society of America↗

USGS2023

On the ratio of full‐resonance to square‐root‐impedance amplifications for shear‐wave velocity profiles that are a continuous function of depth

Amplifications of seismic waves traveling upward through a continuous, interface‐free velocity profile are consistently smaller when computed using the square‐root‐impedance (SRI) method than when computed using full‐resonance (FR) calculations. This was found for a wide range of velocity profiles. For realistic profiles, for which the gradient of velocity decreases with depth, the differences are not large, with the ratio of FR/SRI amplifications ranging from about 1.05 to 1.3. Comparisons of the amplifications from a continuous velocity profile with those from approximations to that profile using a stack of constant‐velocity layers give some support to the hypothesis that the difference between FR and SRI amplifications for gradient profiles is because the former is controlled by the ratio of seismic impedances, whereas the latter is based on the square root of the seismic impedance ratios. This implies that gradient profiles will always have FR amplifications greater than SRI amplifications. A model‐independent, easy‐to‐implement modification of the SRI amplifications is proposed that shows promise in bringing the SRI amplifications closer to the FR amplifications.

Bulletin of the Seismological Society of America↗

USGS2023

Remote characterization of the 12 January 2020 eruption of Taal Volcano, Philippines, using seismo-acoustic, volcanic lightning, and satellite observations

On 12 January 2020, an eruption began on the shores of the Main Crater Lake (MCL) of Taal Volcano—a caldera system on the southern end of Luzon Island in the Philippines. Taal, one of the most active volcanoes in the Philippines, is located 30 km south of Manila—a major metropolitan area with a population of 13.5 million people. Eruptive activity intensified throughout the day on 12 January, producing prolific volcanic lightning, ashfall, and a sustained plume that reached 16–17 km altitude. The chronology of events was well documented by the Philippine Institute of Volcanology and Seismology and the Tokyo Volcanic Ash Advisory Center. The wealth of data collected during the eruption provides a unique opportunity to investigate how the combination of different remote sensing methods may complement local observations and monitoring. Remote systems tend to provide lower resolution data but are also less likely to be compromised by the eruptive activity, thus providing continuous records of eruptive processes. Here, we present a postevent analysis of the 12 January activity, including data from long‐range lightning, infrasound, and seismic arrays located at distances up to several thousands of kilometers from the volcano. By combining these datasets, we distinguish five phases of activity and infer a major shift in eruption behavior around 12:00 on 12 January (UTC). The remote observations suggest that the most of the water within the MCL ( ⁠∼42 million m 3⁠ ) was vaporized and incorporated into the volcanic plume within the first 12 hr of the eruption.

Taal Volcano↗

USGS2023

A detailed view of the 2020-2023 southwestern Puerto Rico seismic sequence with deep learning

The 2020–2023 southwestern Puerto Rico seismic sequence, still ongoing in 2023, is remarkable for its multiple‐fault rupture complexity and elevated aftershock productivity. We applied an automatic workflow to continuous data from 43 seismic stations in Puerto Rico to build an enhanced earthquake catalog with ∼180,000 events for the 3+ yr sequence from 28 December 2019 to 1 January 2023. This workflow contained the EQTransformer (EQT) deep learning model for event detection and phase picking, the EikoNet‐Hypocenter Inversion with Stein Variational Inference probabilistic earthquake location approach with a neural network trained to solve the eikonal wave equation, and relocation with event‐pair waveform cross correlation. EQT increased the number of catalog events in the sequence by about seven times, though its performance was not quite as good as thorough analyst review. The enhanced catalog revealed new structural details of the sequence space–time evolution, including sudden changes in activity, on a complex system of many small normal and strike‐slip faults. This sequence started on 28 December 2019 with an M 4.7 strike‐slip earthquake followed by 10 days of shallow strike‐slip foreshocks, including several M 5+ earthquakes, in a compact region. The oblique normal fault M w "> M w 6.4 mainshock then happened on 7 January 2020. Early aftershocks in January 2020, with several M 5+ earthquakes, quickly expanded into two intersecting fault zones with diffuse seismicity: one extending ∼35 km on a northward‐dipping normal fault and the other ∼60‐km‐long and oriented west‐northwest–east‐southeast on strike‐slip faults. Months to years later, aftershocks moved westward, deeper, and to outer reaches of the active fault zones, with abrupt rapid seismicity migration following larger M 4.7+ aftershocks in May, July, and December 2020. The observed seismicity evolution indicates cascading failure from stress transfer on multiple critically stressed faults. High aftershock productivity results from the complex multiple‐fault network hosting the sequence, which is characteristic of an immature fault system in the diffuse deformation zone around Puerto Rico, at the complicated North American–Caribbean plate boundary region.

Puerto Rico↗

USGS2023

Slip deficit rates on southern Cascadia faults resolved with viscoelastic earthquake cycle modeling of geodetic deformation

The fore‐arc of the southern Cascadia subduction zone (CSZ), north of the Mendocino triple junction (MTJ), is home to a network of Quaternary‐active crustal faults that accumulate strain due to the interaction of the North American, Juan de Fuca (Gorda), and Pacific plates. These faults, including the Little Salmon and Mad River fault (LSF and MRF) zones, are located near the most populated parts of California’s north coast and show paleoseismic evidence for three slip events of several‐meter scale in the past 1700 yr. However, the geodetic slip rates of these faults are poorly constrained. In this work, we analyze a new compilation of interseismic geodetic velocities from Global Navigation Satellite Systems, leveling, and tide gauge data near the MTJ to constrain present‐day slip deficit rates on upper‐plate faults and coupling on the megathrust. We construct Green’s functions for interseismic slip deficit for discrete faults embedded in an elastic plate overlying a viscoelastic mantle. We then use a constrained least‐squares inversion to determine best‐fitting slip rates on the major faults and investigate slip rate trade‐offs between faults. Results indicate that the LSF and MRF systems together accumulate 4–5 mm/yr of reverse‐slip deficit, although their separate slip rates cannot be determined independently. Modeling of the horizontal and vertical velocities suggests that the southernmost CSZ is coupled interseismically to deeper than 25 km depth. We also find that 6–17 mm/yr of right‐lateral slip deficit extends north of the MTJ and into the southern Cascadia fore‐arc. These results reinforce the notion that both the southernmost Cascadia megathrust and the smaller fore‐arc faults above it contribute to regional seismic hazard.

California, Oregon, Washington↗

USGS2023

Introduction to the special section on seismoacoustics and seismoacoustic data fusion

A variety of geophysical hazards (e.g., volcanic activity, earthquakes, mass movements, marine storms, bolides) and anthropogenic sources (e.g., chemical and nuclear explosions, mining blasts, rocket launches) can release energy as mechanical waves in the ground, ocean, and atmosphere (Arrowsmith et al., 2010; Campus and Christie, 2009). Due to the mechanical coupling between a planetary body, its ocean, and its atmosphere, waves can propagate across these interfaces (Ben-Menahem and Singh, 1981) and carry information about the source and the media they propagated through. The field of seismoacoustics, driven by geophysical observations of both seismic and low-frequency acoustic (infrasound) waves, has several interdisciplinary applications. Observations of both seismic and infrasonic waves can be used to discriminate between atmospheric and subsurface events, such as sonic booms and earthquakes. Moreover, seismoacoustic analyses can provide useful information for the source characterization of shallow anthropogenic events, such as underground or surface explosions, volcanic, and tectonic events (e.g., Arrowsmith et al., 2020; Assink et al., 2016; de Groot-Hedlin and Hedlin, 2019; Matoza et al., 2009). Similarly, remote observations (e.g., on regional seismic and infrasonic arrays) can help monitor natural events such as volcanic eruptions and provide additional details about eruption dynamics. Recent works additionally suggest that infrasound can be used to discriminate between different volcanic processes (Watson et al., 2022). Finally, looking outwards, the study of seismo-acoustics has been a particularly valuable tool for planetary science (e.g., Krishnamoorthy et al., 2018; Martire et al., 2018; Krishnamoorthy et al., 2019; Martire et al., 2020; Garcia et al., 2020; Brissaud et al., 2021).

Bulletin of the Seismological Society of America↗

USGS2023

A population-based performance evaluation of the ShakeAlert earthquake early warning system for M 9 megathrust earthquakes in the Pacific Northwest, U.S.A.

We evaluate the potential performance of the ShakeAlert earthquake early warning system for M 9 megathrust earthquakes in the Pacific Northwest (PNW) using synthetic seismograms from 30 simulated M 9 earthquake scenarios on the Cascadia subduction zone. The timeliness and accuracy of source estimates and effectiveness of ShakeAlert alert contours are evaluated with a station‐based alert classification scheme using an alert threshold equal to the target threshold. We develop a population‐based alert classification method by aligning a population grid with Voronoi diagrams computed from the station locations for each scenario. Using raster statistics, we estimate the PNW population that would receive timely accurate alerts during an offshore M 9 earthquake. We also examine the range of expected warning times with respect to the spatial distribution of the population. Results show that most of the population in our evaluation region could receive alerts with positive warning times for an alert threshold of modified Mercalli intensity (MMI) III, but that late and missed alerts increase because the alert threshold is increased. An average of just under 60% of the population would be alerted for MMI V prior to the arrival of threshold level shaking. Large regions of late and missed alerts for thresholds MMI IV and V are caused by delays in alert updates, inaccurate FinDer source estimates, and undersized alert contours due to magnitude underestimation. We also investigate an alerting strategy where ShakeAlert sends out an alert to the entire evaluation region when the system detects at least an M 8 earthquake along the coast. Because large magnitude offshore earthquakes are rare in Cascadia, overalerting is most likely to occur from an overestimated M 7+ on the Gorda plate. With appropriate criteria to minimize overalerting, this strategy may eliminate all missed and late alerts except at sites close to the epicenter.

California, Oregon, Washington↗