SEARCH · USGS Science

Results for “Tectonics”

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 1,567 records · Page 87Linked to original sources

USGS2002

Subsurface geometry and evolution of the Seattle fault zone and the Seattle Basin, Washington

The Seattle fault, a large, seismically active, east-west-striking fault zone under Seattle, is the best-studied fault within the tectonically active Puget Lowland in western Washington, yet its subsurface geometry and evolution are not well constrained. We combine several analysis and modeling approaches to study the fault geometry and evolution, including depth-converted, deep-seismic-reflection images, P-wave-velocity field, gravity data, elastic modeling of shoreline uplift from a late Holocene earthquake, and kinematic fault restoration. We propose that the Seattle thrust or reverse fault is accompanied by a shallow, antithetic reverse fault that emerges south of the main fault. The wedge enclosed by the two faults is subject to an enhanced uplift, as indicated by the boxcar shape of the shoreline uplift from the last major earthquake on the fault zone. The Seattle Basin is interpreted as a flexural basin at the footwall of the Seattle fault zone. Basin stratigraphy and the regional tectonic history lead us to suggest that the Seattle fault zone initiated as a reverse fault during the middle Miocene, concurrently with changes in the regional stress field, to absorb some of the north-south shortening of the Cascadia forearc. Kingston Arch, 30 km north of the Seattle fault zone, is interpreted as a more recent disruption arising within the basin, probably due to the development of a blind reverse fault.

Washington↗

USGS2003

Gravity and magnetic expression of the San Leandro gabbro with implications for the geometry and evolution of the Hayward Fault zone, northern California

The Hayward Fault, one of the most hazardous faults in northern California, trends north-northwest and extends for about 90 km along the eastern San Francisco Bay region. At numerous locations along its length, distinct and elongate gravity and magnetic anomalies correlate with mapped mafic and ultramafic rocks. The most prominent of these anomalies reflects the 16-km-long San Leandro gabbroic block. Inversion of magnetic and gravity data constrained with physical property measurements is used to define the subsurface extent of the San Leandro gabbro body and to speculate on its origin and relationship to the Hayward Fault Zone. Modeling indicates that the San Leandro gabbro body is about 3 km wide, dips about 75??-80?? northeast, and extends to a depth of at least 6 km. One of the most striking results of the modeling, which was performed independently of seismicity data, is that accurately relocated seismicity is concentrated along the western edge or stratigraphically lower bounding surface of the San Leandro gabbro. The western boundary of the San Leandro gabbro block is the base of an incomplete ophiolite sequence and represented at one time, a low-angle roof thrust related to the tectonic wedging of the Franciscan Complex. After repeated episodes of extension and attenuation, the roof thrust of this tectonic wedge was rotated to near vertical, and in places, the strike-slip Hayward Fault probably reactivated or preferentially followed this pre-existing feature. Because earthquakes concentrate near the edge of the San Leandro gabbro but tend to avoid its interior, we qualitatively explore mechanical models to explain how this massive igneous block may influence the distribution of stress. The microseismicity cluster along the western flank of the San Leandro gabbro leads us to suggest that this stressed volume may be the site of future moderate to large earthquakes. Improved understanding of the three-dimensional geometry and physical properties along the Hayward Fault will provide additional constraints on seismic hazard probability, earthquake modeling, and fault interactions that are applicable to other major strike-slip faults around the world.

California↗

USGS2004

Frequency-dependent Lg Q within the continental United States

Frequency-dependent crustal attenuation (1/ Q ) is determined for seven distinct physiographic/tectonic regions of the continental United States using high-quality Lg waveforms recorded on broadband stations in the frequency band 0.5 to 16 Hz. Lg attenuation is determined from time-domain amplitude measurements in one-octave frequency bands centered on the frequencies 0.75, 1.0, 3.0, 6.0, and 12.0 Hz. Modeling errors are determined using a delete-j jackknife resampling technique. The frequency-dependent quality factor is modeled in the form of Q = Q 0 f η. Regions were initially selected based on tectonic provinces but were eventually limited and adjusted to maximize ray path coverage in each area. Earthquake data was recorded on several different networks and constrained to events occurring within the crust (<40 km depth) and at least m b 3.5 in size. A singular value decomposition inversion technique was applied to the data to simultaneously solve for source and receiver terms along with Q for each region at specific frequencies. The lowest crustal Q was observed in northern and southern California where Q is described by the functions Q = 152(±37) f 0.72(±0.16) and Q = 105(±26) f 0.67(±0.16) , respectively. The Basin and Range Province, Pacific Northwest, and Rocky Mountain states also display lower Q and a strong frequency dependence characterized by the functions Q = 200(±40) f 0.68(±0.12) , Q = 152(±49) f 0.76(±0.18) , and Q = 166(±37) f 0.61(±0.14) , respectively. In contrast, in the central and northeast United States Q functions are Q = 640(±225) f 0.344(±0.22) and Q = 650(±143) f 0.36(±0.14) , respectively, show a high crustal Q and a weaker frequency dependence. These results improve upon previous Lg modeling by subdividing the United States into smaller, distinct tectonic regions and using significantly more data that provide improved constraints on frequency-dependent attenuation and errors. A detailed attenuation map of the continental United States can provide significant input into hazard map mitigation. Both scattering and intrinsic attenuation mechanisms are likely to play a comparable role in the frequency range considered in the study.

Bulletin of the Seismological Society of America↗

USGS2005

Stress distribution along the Fairweather-Queen Charlotte transform fault system

Tectonic loading and Coulomb stress transfer are modeled along the right-lateral Fairweather-Queen Charlotte transform fault system using a threedimensional boundary element program. The loading model includes slip below 12 km along the transform as well as motion of the Pacific plate, and it is consistent with most available Global Positioning System (GPS) displacement rate data. Coulomb stress transfer is shown to have been a weak contributing factor in the failure of the southeastern (Sitka) segment of the Fairweather fault in 1972, hastening the occurrence of the earthquake by only about 8 months. Failure of the Sitka segment was enhanced by a combination of cumulative loading from below (95%) by slip of about 5 cm/yr since 1848, by stress transfer (about 1%) from major earthquakes on straddling segments of the Queen Charlotte fault (M 8.1 in 1949) and the Fairweather fault (M 7.8 in 1958), and by viscoelastic relaxation (about 4%) following the great 1964 Alaska earthquake, modeled by Pollitz et al. (1998). Cumulative stress increases in excess of 7 MPa at a depth of 8 km are projected prior to the M 7.6 earthquake. Coulomb stress transferred by the rupture of the great M 9.2 Alaska earthquake in 1964 (Bufe, 2004a) also hastened the occurrence of the 1972 event, but only by a month or two. Continued tectonic loading over the last half century and stress transfer from the M 7.6 Sitka event has resulted in restressing of the adjacent segments by about 3 MPa at 8 km depth. The occurrence of a M 6.8 earthquake on the northwestern part of the Queen Charlotte fault on 28 June 2004, the largest since 1949, also suggests increased stress. The Cape St. James segment of the fault immediately southeast of the 1949 Queen Charlotte rupture has accumulated about 6 MPa at 8 km through loading since 1900 and stress transfer in 1949. A continued rise in earthquake hazard is indicated for the Alaska panhandle and Queen Charlotte Islands region in the decades ahead as the potential for damaging earthquakes increases.

Bulletin of the Seismological Society of America↗

USGS2004

Inverse kinematic and forward dynamic models of the 2002 Denali fault earthquake, Alaska

We perform inverse kinematic and forward dynamic models of the M 7.9 2002 Denali fault, Alaska, earthquake to shed light on the rupture process and dynamics of this event, which took place on a geometrically complex fault system in central Alaska. We use a combination of local seismic and Global Positioning System (GPS) data for our kinematic inversion and find that the slip distribution of this event is characterized by three major asperities on the Denali fault. The rupture nucleated on the Susitna Glacier thrust fault, and after a pause, propagated onto the strike-slip Denali fault. Approximately 216 km to the east, the rupture abandoned the Denali fault in favor of the more southwesterly directed Totschunda fault. Three-dimensional dynamic models of this event indicate that the abandonment of the Denali fault for the Totschunda fault can be explained by the Totschunda fault's more favorable orientation with respect to the local stress field. However, a uniform tectonic stress field cannot explain the complex slip pattern in this event. We also find that our dynamic models predict discontinuous rupture from the Denali to Totschunda fault segments. Such discontinuous rupture helps to qualitatively improve our kinematic inverse models. Two principal implications of our study are (1) a combination of inverse and forward modeling can bring insight into earthquake processes that are not possible with either technique alone, and (2) the stress field on geometrically complex fault systems is most likely not due to a uniform tectonic stress field that is resolved onto fault segments of different orientations; rather, other forms of stress heterogeneity must be invoked to explain the observed slip patterns.

Alaska↗

USGS2008

Dynamic stresses, Coulomb failure, and remote triggering

Dynamic stresses associated with crustal surface waves with 15-30-sec periods and peak amplitudes < 1 MPa are capable of triggering seismicity at sites remote from the generating mainshock under appropriate conditions. Coulomb failure models based on a frictional strength threshold offer one explanation for instances of rapid-onset triggered seismicity that develop during the surface-wave peak dynamic stressing. Evaluation of the triggering potential of surface-wave dynamic stresses acting on critically stressed faults using a Mohr's circle representation together with the Coulomb failure criteria indicates that Love waves should have a higher triggering potential than Rayleigh waves when incident on vertical, strike-slip faults. That (1) the onset of triggered seismicity often appears to begin during the Rayleigh wave rather than the earher-arriving Love wave, (2) vertical strike-slip faults pervade the crust in most tectonic regimes, and (3) Love-wave amplitudes typically exceed those for Rayleigh waves suggests that the explanation for rapid-onset dynamic triggering may not reside solely with a simple static-threshold friction mode. The results also indicate that thrust faults should be more susceptible to dynamic triggering by Rayleigh-wave stresses than normal faults in the shallow seismogenic crust (< 5 km), while the reverse should be true in the lower seismogenic crust (> 5 km). The latter is consistent with the observation that extensional or transtensional tectonic regimes are more susceptible to remote triggering by Rayleigh-wave dynamic stresses than compressional or transpressional regimes. Locally elevated pore pressures may have a role in the observed prevalence of dynamic triggering in extensional regimes and geothermal/volcanic systems.

Bulletin of the Seismological Society of America↗

USGS2011

LiDAR-Assisted identification of an active fault near Truckee, California

We use high-resolution (1.5-2.4 points/m2) bare-earth airborne Light Detection and Ranging (LiDAR) imagery to identify, map, constrain, and visualize fault-related geomorphology in densely vegetated terrain surrounding Martis Creek Dam near Truckee, California. Bare-earth LiDAR imagery reveals a previously unrecognized and apparently youthful right-lateral strike-slip fault that exhibits laterally continuous tectonic geomorphic features over a 35-km-long zone. If these interpretations are correct, the fault, herein named the Polaris fault, may represent a significant seismic hazard to the greater Truckee-Lake Tahoe and Reno-Carson City regions. Three-dimensional modeling of an offset late Quaternary terrace riser indicates a minimum tectonic slip rate of 0.4 ?? 0.1 mm/yr.Mapped fault patterns are fairly typical of regional patterns elsewhere in the northern Walker Lane and are in strong coherence with moderate magnitude historical seismicity of the immediate area, as well as the current regional stress regime. Based on a range of surface-rupture lengths and depths to the base of the seismogenic zone, we estimate a maximum earthquake magnitude (M) for the Polaris fault to be between 6.4 and 6.9.

Bulletin of the Seismological Society of America↗

USGS2011

The 2007 Nazko, British Columbia, earthquake sequence: Injection of magma deep in the crust beneath the Anahim volcanic belt

On 9 October 2007, an unusual sequence of earthquakes began in central British Columbia about 20 km west of the Nazko cone, the most recent (circa 7200 yr) volcanic center in the Anahim volcanic belt. Within 25 hr, eight earthquakes of magnitude 2.3-2.9 occurred in a region where no earthquakes had previously been recorded. During the next three weeks, more than 800 microearthquakes were located (and many more detected), most at a depth of 25-31 km and within a radius of about 5 km. After about two months, almost all activity ceased. The clear P- and S-wave arrivals indicated that these were high-frequency (volcanic-tectonic) earthquakes and the b value of 1.9 that we calculated is anomalous for crustal earthquakes but consistent with volcanic-related events. Analysis of receiver functions at a station immediately above the seismicity indicated a Moho near 30 km depth. Precise relocation of the seismicity using a double-difference method suggested a horizontal migration at the rate of about 0:5 km=d, with almost all events within the lowermost crust. Neither harmonic tremor nor long-period events were observed; however, some spasmodic bursts were recorded and determined to be colocated with the earthquake hypocenters. These observations are all very similar to a deep earthquake sequence recorded beneath Lake Tahoe, California, in 2003-2004. Based on these remarkable similarities, we interpret the Nazko sequence as an indication of an injection of magma into the lower crust beneath the Anahim volcanic belt. This magma injection fractures rock, producing high-frequency, volcanic-tectonic earthquakes and spasmodic bursts.

Bulletin of the Seismological Society of America↗

USGS2012

Contemporary seismicity in and around the Yakima-Fold-and-Thrust Belt in eastern Washington

We examined characteristics of routinely cataloged seismicity from 1970 to the present in and around the Yakima fold‐and‐thrust belt (YFTB) in eastern Washington to determine if the characteristics of contemporary seismicity provide clues about regional‐scale active tectonics or about more localized, near‐surface processes. We employed new structural and hydrologic models of the Columbia River basalts (CRB) and found that one‐third to one‐half of the cataloged earthquakes occur within the CRB and that these CRB earthquakes exhibit significantly more clustered, and swarmlike, behavior than those outside. These results and inferences from published studies led us to hypothesize that clustered seismicity is likely associated with hydrologic changes in the CRB, which hosts the regional aquifer system. While some general features of the regional groundwater system support this hypothesis, seismicity patterns and mapped long‐term changes in groundwater levels and present‐day irrigation neither support nor refute it. Regional tectonic processes and crustal‐scale structures likely influence the distribution of earthquakes both outside and within the CRB as well. We based this inference on qualitatively assessed alignments between the dominant northwest trends in the geologic structure and the seismicity generally and between specific faults and characteristics of the 2009 Wooded Island swarm and aseismic slip, which is the only cluster studied in detail and the most vigorous since regional monitoring began.

Washington↗

USGS2015

A century of induced earthquakes in Oklahoma?

Seismicity rates have increased sharply since 2009 in the central and eastern United States, with especially high rates of activity in the state of Oklahoma. Growing evidence indicates that many of these events are induced, primarily by injection of wastewater in deep disposal wells. The upsurge in activity has raised two questions: What is the background rate of tectonic earthquakes in Oklahoma? How much has the rate varied throughout historical and early instrumental times? In this article, we show that (1) seismicity rates since 2009 surpass previously observed rates throughout the twentieth century; (2) several lines of evidence suggest that most of the significant earthquakes in Oklahoma during the twentieth century were likely induced by oil production activities, as they exhibit statistically significant temporal and spatial correspondence with disposal wells, and intensity measurements for the 1952 El Reno earthquake and possibly the 1956 Tulsa County earthquake follow the pattern observed in other induced earthquakes; and (3) there is evidence for a low level of tectonic seismicity in southeastern Oklahoma associated with the Ouachita structural belt. The 22 October 1882 Choctaw Nation earthquake, for which we estimate Mw 4.8, occurred in this zone.

Oklahoma↗

USGS2017

Influence of lithostatic stress on earthquake stress drops in North America

We estimate stress drops for earthquakes in and near the continental United States using the method of spectral ratios. The ratio of acceleration spectra between collocated earthquakes recorded at a given station removes the effects of path and recording site and yields source parameters including corner frequency for, and the ratio of seismic moment between, the two earthquakes. We determine stress drop from these parameters for 1121 earthquakes greater than M ∼3 in 60 earthquake clusters. We find that the average Brune stress drop for the few eastern United States (EUS) tectonic mainshocks studied (2.6–36 MPa) is about three times greater than that of tectonic mainshocks in the western United States (WUS, 1.0–7.9 MPa) and five times greater than mainshocks potentially induced by wastewater injection in the central United States (CUS, 0.6–5.6 MPa). EUS events tend to be deeper thrusting events, whereas WUS events tend to be shallower but have a wide range of focal mechanisms. CUS events tend to be shallow with strike‐slip to normal‐faulting mechanisms. With the possible exception of CUS aftershocks, we find that differences in stress drop among all events can be taken into account, within one standard deviation of significance, by differences in the shear failure stress as outlined by Mohr–Coulomb theory. The shear failure stress is a function of vertical stress (or depth), the fault style (normal, strike slip, or reverse), and coefficient of friction (estimated here to be, on average, 0.64). After accounting for faulting style and depth dependence, we find that the average Brune stress drop is about 3% of the failure stress. These results suggest that high‐frequency shaking hazard (>∼1 Hz) from shallow induced events and aftershocks is reduced to some extent by lower stress drop. However, the shallow hypocenters will increase hazard within several kilometers of the source.

Bulletin of the Seismological Society of America↗

USGS2016

Relations between some horizontal‐component ground‐motion intensity measures used in practice

Various measures using the two horizontal components of recorded ground motions have been used in a number of studies that derive ground‐motion prediction equations and construct maps of shaking intensity. We update relations between a number of these measures, including those in Boore et al. (2006) and Boore (2010) , using the large and carefully constructed global database of ground motions from crustal earthquakes in active tectonic regions developed as part of the Pacific Earthquake Engineering Research Center–Next Generation Attenuation‐West2 project. The ratios from the expanded datasets generally agree to within a few percent of the previously published ratios. We also provide some ratios that were not considered before, some of which will be useful in applications such as constructing ShakeMaps. Finally, we compare two important ratios with those from a large central and eastern North American database and from many records from subduction earthquakes in Japan and Taiwan. In general, the ratios from these regions are within several percent of those from crustal earthquakes in active tectonic regions.

Bulletin of the Seismological Society of America↗

USGS2019

On the portability of ML-MC as a depth discriminant for small seismic events recorded at local distances

In this paper we show that M L -M C is a viable and regionally portable depth discriminant and therefore may contribute in nuclear test ban treaty verification. A recent study found that the difference between local magnitude (M L ) and coda duration magnitude (M C ) discriminates shallow seismic events (mining blasts, mining-induced earthquakes, and shallow tectonic earthquakes) from deeper tectonic earthquakes in the Utah region. The shallow seismic events had anomalously high M C values, with increasingly negative M L -M C values as depth decreased. Here we evaluate the performance of M L -M C as a depth discriminant in three new regions, finding that M L -M C increases between 0–9 km depth in all cases. Initially, we investigated M L -M C as a function of depth for naturally occurring earthquakes in the region around Yellowstone National Park, as recorded by the University of Utah Seismograph Stations. For 3,358 Yellowstone earthquakes with well-constrained depths, we found M L -M C increased 0.030 ± 0.007 magnitude units (m.u.) for each 1 km increase in depth up to 10 km depth. Next, we examined M L -M C values for anthropogenic seismicity in northern Oklahoma and southern Kansas, as recorded by the National Earthquake Information Center. For 1,628 events with well-constrained depths, we computed a slope for M L -M C of 0.022 ± 0.010 m.u./km. Finally, we analyzed M L -M C for 28,722 well-located earthquakes in Italy, as recorded by the National Institute of Geophysics and Volcanology, and found an M L -M C slope of 0.018 ± 0.001 m.u./km. In each case, the quoted error bounds represent 95% confidence regions which exclude zero, implying that the depth-dependence of M L -M C is statistically significant. We performed several robustness tests in which we varied the criterion used to define a well-constrained depth and the depth range used in the linear fit. In nearly all cases, we found a positive slope for M L -M C vs. depth at a confidence level above 95%.

Utah, Montana, Idaho, Wyoming, Kansas, Oklahoma↗

USGS2026

How to determine an earthquake rate budget? A comparison of geodetic, geologic, and seismologic moment rates

How much tectonic geodetic strain is accommodated as earthquakes? Can estimates of tectonic geodetic strain be used in probabilistic seismic hazard analyses? Understanding the balance of moment accumulation (i.e., geodetic moment) and accommodation (i.e., geologic and seismologic moment) is a necessary step toward leveraging corrected geodetic strain rates to quantify earthquake budgets. We investigate these questions using updated geodetic, geologic, and seismologic datasets compiled for the western United States (WUS) portion of the 2023 U.S. National Seismic Hazard Model (NSHM23‐WUS) to provide new views of moment rates across the WUS. We compare moment rates using a floating footprint approach across the extent of the NSHM23‐WUS fault sections database. We make these comparisons using “moment rate deviation (MRD),” which is defined as a surplus or deficit of geodetic moment rate relative to the sum of geologic and seismologic moment rates. We find that, in general, higher strain rate areas of the WUS (i.e., plate boundary) have smaller values of MRD, meaning that geodetic moment rates approximate the earthquake record preserved by geologic and seismologic moment rates well. In contrast, the lower strain rate, eastern portions of our analysis domain (i.e., intraplate) have very high values of MRD, with geodetic moment rates exceeding measured geologic and seismologic moment rates. We find that geodetic moment rates appear to be representative of the earthquake budget in relatively well‐characterized, interplate regions, and geodetic estimates of off‐fault deformation could be considered in seismic hazard analyses where MRD is minimized.

western United States↗

USGS2018

Ground motions from induced earthquakes in Oklahoma and Kansas

Improved predictions of earthquake ground motions are critical to advancing seismic hazard analyses and earthquake response. The high seismicity rate from 2009 to 2016 in Oklahoma and Kansas provides an extensive data set for examining the ground motions from these events. We evaluate the ability of three suites of ground‐motion prediction equations (GMPEs)—appropriate for modeling tectonic earthquakes in active crustal and stable continental regions—to reproduce the observed ground motions. Mixed‐effects regressions are used to separate the ground‐motion residuals into bias, between‐event, and within‐event terms. Although the residuals depict differing accuracies in the ability of the three GMPE suites to predict the ground motions, some consistent trends emerge in the period, magnitude, and distance dependence. The trends suggest that aspects of the ground motions from these induced earthquakes are not well modeled by current tectonic GMPEs. Most important, we find evidence for relatively poor overall fit to the ground motions, by all of the GMPE suites, at periods less than about 0.2 s and above 3 s, greater‐than‐predicted magnitude scaling for small to moderate‐magnitude events ( ⁠ M ≲ 5 "> M ≲ 5 ⁠ ), higher‐than‐predicted within‐event variability above 3 s, and an apparent geometric attenuation that is stronger than average predictions at close distances ( ⁠ R ≲ 20    km "> R ≲ 20 km ) and short periods ( ⁠ T ≲ 1 s "> T ≲ 1 s ). Our results indicate that regionally appropriate GMPEs for induced earthquakes should be used for predicting ground motions in Oklahoma and Kansas.

Kansas, Oklahoma↗

USGS2024

(Re)discovering the seismicity of Antarctica: A new seismic catalog for the southernmost continent

We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on‐continent temporary and long‐term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (⁠Mw ⁠) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice‐covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Seismological Research Letters↗

USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Cascadia Subduction Zone science: Call for the next generation community seismic velocity model

The Cascadia subduction zone (CSZ) hosts major seismic and tsunami hazards, yet key questions persist about the relationship between margin structure, fluid distribution, episodic tremor and slip, shallow megathrust behavior, shaking and tsunamigenesis, and the resulting hazard estimates. Addressing these problems requires an empirically grounded, three‐dimensional seismic velocity model to illuminate subsurface structure and properties and to provide a basis for geophysical studies such as earthquake simulations and ground‐motion estimation. In May 2024, the National Science Foundation‐funded Cascadia Region Earthquake Science Center (CRESCENT) community velocity model (CVM) working group, with U.S. Geological Survey and regional partners, convened a workshop to identify priorities for such a model. Participants emphasized the features necessary for addressing key science questions, including implementing findability, accessibility, interoperability, and reusability (FAIR) access, capturing along‐strike and along‐dip structural heterogeneity, resolving shallow offshore–onshore structure, constraining elastic properties and quantifying their uncertainties for numerical wave propagation simulations, their validation benchmarks, and supporting associated accurate earthquake ground‐motion simulations and hazard assessments. This article describes the priorities defined in the workshop, and a description of how, guided by these needs, CRESCENT plans to develop multiple generations of a CVM to advance CSZ science and improve seismic and tsunami hazard modeling across the Pacific Northwest. The CVM will span the CSZ from the surface to ∼100 km depth, offshore and east of the Cascades into Idaho (∼132°–110° W) and the southern and northern tectonic regime transitions (∼36°–52° N) to capture the entire tectonic system as well as its surroundings.

Cascadia Subduction Zone↗