Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of the Nevada Water Resources Association”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

1,691 records · Page 68Linked to original sources
USGS2024

Methods to evaluate and improve the modeling of rupture directivity in assessment of seismic hazard

In recent years, there have been several advancements related to the modelling of near-source effects of earthquake rupture on strong ground shaking, leading to an improved characterization of ground motions and resulting seismic hazard. Some of these modifications have stemmed from physics-based numerical modelling of the earthquake rupture process, using physics-based dynamic rupture simulations. These contributions have led to a better understanding of how fault rupture characteristics, geometry, and the style of faulting can interact with the hypocenter-dependence on the path from source to site that may ultimately guide the development of seismic directivity models. Moving forward, the application of modern techniques can be used to incorporate these source characteristics and near-fault ground motion behavior that contribute to the azimuthally varying effects that result in rupture directivity. One example is the application of machine learning methods to support more automated integration of new predictor variables in model development and open more evaluation opportunities to access residuals. Here, we utilize several techniques to take advantage of the plethora of synthetic data and its ability to supplement preexisting trends observed in data. We showcase two examples of how models can be either developed, expanded upon, or constrained using artificial neural network model (ANNs). We evaluate the performance of the ANN with existing methods, comparing misfit, potential limitations, and ability to continue to improve upon these methods in the future. One approach uses a set of simulations with corresponding synthetic ground motions from the Southern California Earthquake Center (SCEC) CyberShake study to develop a ground motion model adapted to incorporate seismic directivity information using an ANN. This large database (TBs) enables us to train the model to capture magnitude, period, and distance variations and how these parameters relate to amplification from hypocenters located along finite-faults. In some cases, there is reduced misfit from better representing source features that aren’t included in base ground motion models that neglect hypocenter location (e.g. azimuthal variation, source-to-site terms). Another ANN method uses a shallow-layered neural network model to better fit a hypocenter-independent model. This method adjusts the median and aleatory variability to account for the averaged impact of various hypocenter distributions to fit the underlying directivity adjustment model. This method serves as a template to apply to other directivity models, improving computational efficiency and more readily enabling integration in hazard codes.

California↗
USGS2025

Solution-collapse breccia pipe uranium deposits of the southern Colorado Plateau, northwestern Arizona, USA

Some of the highest-grade uranium deposits in the United States occur in breccia pipes that formed by solution and collapse of sedimentary strata, which occur in the southern portion of the Colorado Plateau in northwestern Arizona. The host breccia pipes are up to 1200 m in vertical extent, average about 90 m in diameter, and can cross-cut strata from their base in the Mississippian Redwall Limestone to as stratigraphically high on some plateaus as the Triassic Chinle Formation. These uranium-base metal deposits are up to 600 m thick and formed within the breccia pipes where they transect the Permian Coconino Sandstone, Hermit Formation, and the Esplanade Sandstone. Of the hundreds of breccia pipes identified across this region, only a small percentage are known to contain mineralization. The main uranium ore mineral is uraninite that is intergrown with at least 20 base-metal sulfide minerals, which contribute Fe, Cu, Co, As, Pb, Zn, Ni, and Ag to the deposits. This study considered regional stratigraphy, sulfur isotope systematics, mineralogy, in situ dating, and compilation and analysis of previous work on the deposits. A comprehensive deposit model has not been published for these deposits. This analysis identified new additions to update the deposit model for these unusual, possibly unique deposits. Proposed modifications to the model include: (1) the source, mechanisms, timing of the base-metal sulfide mineral assemblages, and (2) the source, mechanism, and timing of the uranium mineralization. Sulfide and uranium deposition are shown to be separate mineralization events. The study proposes the possible role of gypsum as a source of sulfur for the sulfide minerals in the deposits. Groundwaters carrying uranium encountered the preexisting sulfides in breccia pipes, reducing the uranyl ions, and precipitating U oxide (as uraninite). Analysis of the regional stratigraphy recognized that numerous beds of gypsum are in the strata that lie only tens of meters above the breccia pipe deposits. In the breccia pipe region, if these stratigraphic units (Toroweap and Kaibab Formations) do not contain gypsum layers then the underlying pipes are not mineralized; where these Permian gypsum layers do occur, breccia pipes can host mineralization. This new understanding should be useful in identifying the prospective region for mineralized pipes.

Arizona↗
USGS2026

Geologic map of pre-Middle Jurassic basement rocks beneath the Atlantic and Gulf Coastal Plains in Florida

Much of the southeastern United States, including all of Florida, is covered by flat-lying sedimentary strata of the Atlantic and Gulf Coastal Plains which have accumulated since Middle Jurassic time. The pre-Middle Jurassic rocks that underlie these coastal plains in Florida, here collectively referred to as “basement,” are known only from a relatively small number of boreholes. This scientific investigations map presents an interpretation of the basement geology in a 1:1,000,000-scale subsurface geologic map with supporting text, data, and figures. The subsurface mapping methodology integrates petrographic, geochronological, thermochronological, geochemical, and mineralogical analyses of drill cores and cuttings in the context of regional geophysical data. The pre-Middle Jurassic rocks of Florida consist of the Gondwanan (West African) Suwannee terrane which was accreted to Laurentia during the Alleghanian orogeny and subsequently intruded by Permian granites, superposed by early Mesozoic rift basins, and partially overlain by bimodal Jurassic volcanic rocks. The younger basement components, specifically the Southwest Florida volcanic province, North Florida tholeiites, early Mesozoic rift basins, and Alleghanian granitoids, have correlative and contemporaneous units throughout the Appalachian orogen. In contrast, Florida’s older basement rocks, including Paleozoic siliciclastic strata of the Suwannee basin, North Florida volcanic series, Osceola and Gaskin intrusive complexes, and the St. Lucie Metamorphic Complex, have neither surface exposures nor unequivocal correlates. Major structures include early Mesozoic normal faults and northwest-striking transfer zones such as the Jay fault. Many of these faults define the boundaries of subbasins within the South Georgia rift system. Top-of-basement structure contours show gentle arches and embayments that are also recognized in overlying coastal plain strata.

Florida↗
USGS2024

Triggering the 2022 eruption of Mauna Loa

Distinguishing periods of intermittent unrest from the run-up to eruption is a major challenge at volcanoes around the globe. Comparing multidisciplinary monitoring data with mineral chemistry that records the physical and spatio-temporal evolution of magmas fundamentally advances our ability to forecast eruptions. The recent eruption of Mauna Loa, Earth’s largest active volcano, provides a unique opportunity to differentiate unrest from run-up and improve forecasting of future eruptions. After decades of intermittent seismic and geodetic activity over 38 years of repose, Mauna Loa began erupting on 27 November 2022. Here we present a multidisciplinary synthesis that tracks the spatio-temporal evolution of precursory activity by integrating mineral and melt chemistry, fluid inclusion barometry, numerical modeling of mineral zoning, syn-eruptive gas plume measurements, the distribution and frequency of earthquake hypocenters, seismic velocity changes, and ground deformation. These diverse data indicate that the eruption occurred following a 2-month period of sustained magma intrusion from depths of 3–5 km up to 1–2 km beneath the summit caldera, providing a new model of the plumbing system at this very high threat volcano. Careful correlation of both the geochemistry and instrumental monitoring data improves our ability to distinguish unrest from the run-up to eruption by providing deeper understanding of the both the monitoring data and the magmatic system—an approach that could be applied at other volcanic systems worldwide.

Hawaii↗
USGS2024

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington↗
USGS2025

Fault displacement model for surface principal rupture of strike-slip faults

The probability distribution model for principal displacement accommodated on the surface main trace is a critical input to the fault displacement hazard analysis. This article presents a new model for strike-slip ruptures in the moment magnitude ( M ) range of 6 to 8.3. The new model is the outcome of a multi-year research effort to update the widely used model developed by Petersen and others in 2011. Updates include the adoption of the Fault Displacement Hazard Initiative database and enhancements to rupture and displacement data preparation. Statistical formulation and estimation have also been updated substantially. A three-parameter modified normal distribution that we refer to as the negative Exponentially Modified Gaussian distribution is adopted to model the probability distribution of the natural logarithm of principal displacement. Formulation for the mean parameter of the modified normal includes a random earthquake term, a nonlinear scaling relation with M , and an ellipse function for along-main-trace variation. The aleatory variability of the updated model now depends on M as well as site’s along-main-trace position. These updates not only significantly improve the fit to the distribution of the observed displacements but also yield reasonable 95th percentile predictions for M > 7.5 events. Alternative models representing the estimation uncertainty of the M -scaling relation are also developed. These new models are compared to the previous model in terms of percentile predictions and the calculated hazard curves. The steeper hazard curves from the new models yield a lower exceedance rate than the normal-distribution based model developed previously by Petersen and others.

Earthquake Spectra↗
USGS2026

Unraveling protracted modification of Archean and Paleoproterozoic crust in central Laurentia, Penokean orogen, with garnet and accessory mineral geochronology and microstructural analysis

Proterozoic metamorphism and deformation of the southern margin of the Superior craton in the Lake Superior region is attributed to the Penokean orogeny (1890−1830 Ma). This model includes a period of crustal inversion in which Archean basement blocks were exhumed through overlying Paleoproterozoic strata, producing the corridor of gneiss domes that parallels the trend of the Penokean orogen across the northern Midcontinent, USA. However, recent geologic mapping and 40 Ar/ 39 Ar geochronology challenge this interpretation, suggesting instead that the gneiss dome structures reflect younger episodes of tectonic activity along the southern margin of Laurentia. In absence of integrated pressure-temperature-time-deformation constraints for these rocks, interpretations are largely limited to their final cooling history, making it difficult to both identify the tectonic forces that shaped the architecture of the Penokean orogenic belt and assess the extent to which later Proterozoic tectonism modified the southern Superior craton. We address this problem with an approach joining thermodynamic modeling, garnet and accessory mineral geochronology, and microstructural analysis for several metamorphic rocks across the gneiss dome corridor. The U-Pb ages of titanite reveal that the Proterozoic geometries of exhumed basement gneiss domes are governed by preexisting Archean structures. Garnet Lu-Hf geochronology constrains the timing of prograde-to-peak metamorphism in the Penokean orogenic belt. Granulite facies metamorphism is related to the final stages of the Penokean orogeny at 1837 Ma and localized in a belt of high-grade rocks near a major Penokean suture. Garnet Lu-Hf ages of samples adjacent to gneiss domes reflect regional metamorphism following the accretionary phase of the Penokean orogeny, between 1825 Ma and 1782 Ma, which we suggest reflects continued crustal thickening related to convergence farther south during this time interval. Combination of garnet microstructures and Sm-Nd ages reflects later exhumation of gneiss domes and buried metasedimentary rocks by ca. 1750 Ma, consistent with previously published 40 Ar/ 39 Ar cooling ages across the region. Reset Lu-Hf and Sm-Nd garnet ages and U-Pb ages of syn-kinematic titanite reflect reactivation of primary Penokean structures during this period of basement uplift. These data document significant modification of the Penokean orogen and the Archean crust of the southern Superior province between 1800 Ma and 1700 Ma. Tectonic activity during this interval coincides with collisional events recognized in western Laurentia, suggesting that the period immediately following the Penokean orogeny may be a broadly important time for crustal growth and modification in proto-North America.

Michigan, Wisconsin↗
USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗
USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗
USGS2025

Summary of results from monitoring the Geysers with continuous passive seismic and repeat magnetotelluric measurements (2021-2023)

Understanding temporal variations in a geothermal field can support operators in decision making that pertains to optimizing production and mitigating hazards. Between 2021 and 2023, The Geysers geothermal field in northern California was monitored with an array of continuous passive seismic sensors and annual repeat magnetotelluric (MT) measurements. Each of these data sets were analyzed and modelled separately to understand the data, sensitivity, and any observable changes. Then, the data were inverted jointly using a crossgradient method to further constrain temporal changes in geophysical properties within the geothermal field. Multiple permutations of annual datasets were used as inputs to the joint inversion. Results demonstrate seismic data constrain smooth inversion of the MT data, and the MT data provide supplementary information about the location of temporal fluid changes. Estimating relative changes in steam saturation for various time intervals of the joint models shows compartmentalized changes in the field, and good spatial correlation with the location of injection wells. These results demonstrate that collecting both passive seismic and MT measurements then modeling them jointly provide complementary information and a relatively inexpensive method for monitoring temporal changes in an active geothermal field that provides results to support operators.

California↗
USGS2026

Cambrian and Ordovician stratigraphy, conodont biostratigraphy, and microfacies analysis to support 1:24,000-scale geologic mapping of the southern Lake Champlain valley, New York and Vermont

Introduction Geologic mapping in the southern Lake Champlain valley of New York and Vermont (fig. 1) has required evaluation of stratigraphic nomenclature used since the early 1900s. The paleogeography of the Cambrian and Ordovician adjacent to the Adirondack Highlands, the high-relief mountains of Proterozoic igneous and metamorphic rocks of northern New York, has impacted the lithostratigraphic succession, facies changes, and several hiatuses that occur within the package of clastic and carbonate rocks that demonstrates overall deepening consistent from the southern to northern Appalachians. The stratigraphic nomenclature used for the geologic mapping of the Cambrian and Ordovician strata of the study area is derived from evaluation of names put forth by previous workers from areas of northern New York and western Vermont. Descriptions of type localities or type sections and other criteria set forth by the North American Stratigraphic Code (North American Commission on Stratigraphic Nomenclature [NACSN], 2021) were consulted. The North American Stratigraphic Code states that lithostratigraphic units, such as the fundamental unit, the formation, be defined by lithic character and be mappable (NACSN, 2021, Article 24). However, some early workers (Walcott, 1912; Kay, 1937; Fisher and Hanson, 1951; Oxley and Kay, 1959) have defined stratigraphic units based on fossil content, which may or may not be facies dependent, and hence may not be mappable as lithostratigraphic units. Although preservation of older names takes priority in stratigraphic nomenclature (NACSN, 2021, Article 7), some do not have adequately defined stratotypes or are located a distance away with different lithologies than that in the study area. Biostratigraphy can be a helpful tool for understanding correlations and facies changes. Herein, conodont biostratigraphy is utilized to date and correlate separate lithologic units. Further, microfacies analyses of units through inspection of petrographic thin sections provide paleogeographic information that helps in understanding the genesis of geologic units.

New York, Vermont↗
USGS2025

Relationship of basin structure and bedrock lithology to faulting in the 2019 Ridgecrest earthquake region, California, from gravity and aeromagnetic data

We investigate patterns of cumulative offsets on the faults that ruptured in 2019 and along the Garlock Fault in the Ridgecrest region, California using recently published gravity and aeromagnetic data. We also examine the relationship of basin structure and bedrock structure to the 2019 M7.1 Ridgecrest earthquake ruptures (Fig. 1A), which were primarily along a dextral northwest-striking fault system, and along a sinistral northeast-striking fault, which ruptured hours earlier with a M6.4 event.

California↗
USGS2024

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp / Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ~2-km-wide structure with high Vp / Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ~90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments ( Vp ~1600–3000 m/s, Vs ~500 m/s to ~1000 m/s), underlain by higher velocity basement rocks ( Vp ~3000–5800 m/s; Vs ~1000–1500 m/s). (3) Between ~1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ~3 km west and ~1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp / Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp / Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003) . (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (~3000 m/s to ~4200 m/s) and Vs (~1200–2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity ( Vp ~2500–3000 m/s) Jurassic–Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

California↗
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2025

Rupture process of the Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake

The Mw7.0 December 5, 2024 Offshore Cape Mendocino earthquake ruptured a km long portion of the east-west trending Mendocino fault zone (MFZ). In order to clarify the rupture process, we assemble three-component seismograms from regional seismic stations, horizontal coseismic displacement vectors derived from Global Navigation Satellite System (GNSS) time series, and a Sentinel-1 ascending interferogram. These data are interpreted with a model of slip distributed on two vertical fault planes representative of the eastern MFZ and spanning the ~70 km length of the aftershock zone. Assuming right-lateral strike slip, we find that the rupture initiates in the oceanic mantle at 20-30 km depth and proceeds unilaterally updip and toward the east. Early aftershocks locate adjacent to the peak slip areas, tracking the coseismic rupture propagation from oceanic mantle to shallower depth and implying a significant role of static stress transfer in driving aftershocks in an ocean plate environment.

California↗
USGS2024

The importance of nodal plane orientation diversity for earthquake focal mechanism stress inversions

Inversions of earthquake focal mechanisms are among the most accessible and reliable methods for determining crustal stress. However, the use of this method varies widely, and assumptions that underpin it are often violated, potentially compromising stress estimates. We investigate the consequences of violating the little-studied assumption that the focal mechanisms have diverse orientations. Our approach is to employ data-informed synthetic mechanisms, with nodal plane orientations defined by recent earthquake lineaments in the Midland Basin, western Texas, and rakes consistent with slip in the mapped stress field. Using both the traditional stress inversion method that assumes constant shear stress magnitudes on the causative faults as well as a recently published variable shear stress method, we show that low fault plane diversity can cause maximum horizontal stress ( S Hmax ) orientation and relative principal stress magnitude (faulting regime) estimates to differ markedly from the true values. This problem is compounded for catalogs with even modest amounts of noise (≤15°) or few (e.g., 20) mechanisms. Significantly, traditional approaches for quantifying uncertainty such as the bootstrap can severely underestimate the true uncertainty under these circumstances. To remedy this, we provide simple tools to quantify nodal plane orientation diversity and stress inversion reliability.

Geological Society of London Special Publications↗