Search USGS⌕ Search

SEARCH · Search USGS

Results for “Bulletin, Southern California Academy of Sciences”

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,666 records · Page 2Linked to original sources
USGS2025

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington↗
USGS2025

Relations of groundwater quality to long-term surface disposal of produced water near the Midway-Sunset and Buena Vista Oil Fields, California, USA

Contamination of groundwater by oil-field fluids in proximity to oil and gas development has been an issue of concern to water users and regulators given long histories of development and legacy disposal practices. A robust set of geochemical tracers including petroleum hydrocarbon compounds, thermogenic gases, inorganic ion concentrations, stable isotopes, radioactive isotopes, and noble gases were used to assess if oil-field fluids mixed with groundwater near the Midway-Sunset and Buena Vista Oil Fields in California, USA. Results show evidence of mixing of oil-field fluids with groundwater within the study area from either anthropogenic or natural processes. Produced water plumes associated with modern surface disposal facilities, used since the late 1950s, extend up to 1.5 km and currently remain within the boundaries of the oil fields. Plumes associated with earlier routing of produced water down natural drainages and in large retention structures (Midway Basin and Sunset Basin) near the Buena Vista Lake Bed are present in groundwater east of the oil fields. Based on geochemical tracer evidence, aerial imagery, and aerial electromagnetic surveys, these legacy plumes reach the western portion of the Central Valley aquifer system, an important groundwater resource for agricultural and domestic supply. The legacy plume associated with Sunset Basin may further be detected downgradient in deeper groundwater beneath the southern extent of the Buena Vista Lake Bed based on the presence of thermogenic gases and petroleum hydrocarbon compounds.

California↗
USGS2026

Widespread abyssal turbidites record megathrust earthquake-triggered landslides and coseismic deformation in the Cascadia subduction zone

Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

California, Oregon↗
USGS2025

Near-fault amplification and ground motion variability during the 2019 Ridgecrest, California sequence

We estimate ground-motion variability near the 2019 M 7.1 Ridgecrest earthquake sequence. Accurate seismic hazard estimation requires understanding ground-motion spatial correlations, yet many studies lack the dense station coverage needed to resolve small-scale variability. The 2019 M 7.1 Ridgecrest earthquake sequence presents a unique opportunity to examine ground motions and their spatial correlations at a range of interstation distances. The permanent seismic network was augmented with hundreds of temporary stations including several fault-crossing nodal arrays. We compute the event ( δE i ) and within-event ( δW ij ) residuals from the observed peak ground velocity and peak ground acceleration data to isolate potential sources of ground-motion variability. We then compare δW ij between station pairs that record an event to understand the semivariance of the ground motion versus interstation distance. By fitting an exponential model to the semivariances, we determine a correlation range of 25 km for the Ridgecrest region. Although the exponential model fits the broad-scale increase of semivariance with interstation distance, we also observe smaller-scale trends. We find that ground motions are less correlated for station pairs that are near or across faults that ruptured during the 2019 Ridgecrest sequence. We also find large, positive median δW ij with relative values 2–3 times larger than nearby stations for individual stations’ near-fault traces. Near-fault amplification and greater ground-motion variability can delineate fault zones and may locally increase the seismic hazard.

California↗
USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗
USGS2026

Constraining source and path effects of large magnitude earthquakes using ground motion simulations

The purpose of this study is to use ground‐motion simulations to investigate ways in which source and path effects for large‐magnitude earthquakes can be represented in nonergodic ground‐motion models (GMMs). To achieve this, we designed a ground‐motion study in the San Francisco Bay Area that includes earthquakes with a broad range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of kinematic simulations (magnitude 4–7), we then develop a nonergodic GMM with the simulated data. We find that trends in the within‐site residuals are affected significantly by the earthquake radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within‐site residuals among multiple source realizations. Finally, after removing the source effects from the within‐site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site as the large events, whereas the second approach considers all small events on the fault plane. The results indicate that it is difficult to satisfactorily approximate the path effects of large events with those of small events using either approach, at least in the case of simulations.

California↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2024

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

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

Bulletin of the Seismological Society of America↗
USGS2026

Detecting earthquakes in noisy real-time GNSS data with deep learning for improved PGD magnitude estimation

To disseminate accurate and useful warnings, earthquake early warning (EEW) systems must quickly determine the size and location of an earthquake to estimate expected shaking. Traditional seismic‐based algorithms tend to underestimate the true magnitudes of large earthquakes, a phenomenon known as magnitude saturation. This limitation motivated the recent inclusion of Global Navigation Satellite Systems (GNSS) data into the U.S. Geological Survey’s ShakeAlert EEW system with the Geodetic First Approximation of Size and Time (GFAST) algorithm because GNSS data do not saturate with large ground motions. However, the noise levels of GNSS data are very high compared with traditional seismic data, which obscures P ‐wave arrivals and can result in less accurate magnitude estimations if displacement amplitudes are low, such as for lower magnitude earthquakes or large source–station distances. In this study, we develop a deep‐learning model that detects earthquakes in GNSS data and use the Ridgecrest, California, earthquake sequence as a case study to demonstrate how the model could act as a filter to reduce the amount of low‐quality data that enters an algorithm like GFAST. To preserve our limited real earthquake data for model inference, we generated a training dataset composed of >700,000 synthetic displacement waveforms. We combined the synthetic waveforms with real‐time GNSS noise to produce realistically noisy training waveforms and then tested our model on additional synthetic data and performed inference using the real data that were held back. We discuss the performance of our trained model on both the unseen synthetic data and real inference data. Our model can be used to selectively filter only high‐quality data where an earthquake signal is observed for input into an algorithm like GFAST (outperforming a simple signal‐to‐noise ratio–based filter) to reduce the error in GFAST’s real‐time earthquake magnitude estimations.

California↗
USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗
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↗
USGS2025

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California↗
USGS2023

Preliminary geologic map of the Southern Santa Rosa Mountains and Borrego Badlands, San Diego County, Southern California

This investigation delineates the geologic framework of an area of 75 square kilometers (km 2 ) located west of the Salton Sea in southern California (fig. 1, on sheet 1). The study area encompasses the south flank of the Santa Rosa Mountains and the eastern part of the Borrego Badlands (sheet 1). In this study area, regionally important stratigraphic and structural elements collectively inform the late Cenozoic geologic evolution of the Anza-Borrego sector of the Salton Trough province. Critical stratigraphic and structural elements in the map area include the following: The well exposed sequence of late Cenozoic, nonmarine sedimentary rocks that filled the Anza-Borrego subbasin (fig. 1) of the Salton Trough; A tectonic boundary that—in the southern Santa Rosa Mountains—separates the sedimentary strata from underlying crystalline rocks of Peninsular Ranges type. This tectonic boundary, named the West Salton Detachment Fault System by Axen and Fletcher (1998), is projected to underlie all late Cenozoic sedimentary strata in the Anza-Borrego subbasin of the Salton Trough; A variety of transpressional, transtensional, and strike-slip structures that have deformed the late Cenozoic sedimentary strata and collectively guided syntectonic and posttectonic depositional events within the Anza-Borrego subbasin of the Salton Trough; and The southeasternmost surface expression of the Clark Fault, a major strand of the dextral San Jacinto Fault Zone. Geologic mapping and analysis for this investigation focused on clarifying geologic relations among these four stratigraphic and structural aspects in the map area.

California↗
USGS2025

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California↗
USGS2025

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California↗
USGS2026

ShakeAlert Earthquake Early Warning System performance during the Mw 7.0 offshore Cape Mendocino earthquake

The 5 December 2024 M w 7.0 Offshore Cape Mendocino earthquake was a challenging test of the U.S. West Coast ShakeAlert earthquake early warning system due to its offshore epicenter and limited near‐source station coverage. We analyzed real‐time performance of all components of the ShakeAlert system, including the seismic algorithms (earthquake point‐source integrated code [EPIC] and Finite‐fault rupture Detector [FinDer]), the geodetic algorithm (Geodetic First Approximation of Size and Time–peak ground displacement [GFAST‐PGD]), and network telemetry during the event. EPIC created the first solution for this earthquake 15 s after origin time with an initial magnitude estimate of M 5.6 and location error of 10 km from the Advanced National Seismic System epicenter. An early spurious trigger from station CE.89101 fortuitously maintained location accuracy and, correspondingly, magnitude accuracy. FinDer contributed its first solution at 18 s with a location estimate closer to the seismic network and produced two distinct rupture geometries, leading to minor fluctuations in estimated intensity contours. GFAST‐PGD did not meet alerting thresholds but otherwise performed as expected. Network latencies were <2 s for most stations, supporting the rapid detection of this earthquake by the system. Roughly five million alerts were delivered to cell phone devices in California and Oregon during this event. This was also the first instance of a school district‐wide ShakeAlert‐powered system being activated. Comparisons to recorded seismograms demonstrate that the maximum warning times before potentially damaging shaking (intensity 6+) were in the range of 5–55 s. Although the ShakeAlert system provided accurate solutions and useful alert delivery, this earthquake raised awareness of potential issues within the system, including the need for improved offshore location estimates, a combination of solutions from ShakeAlert servers, and handling of spurious triggers.

California↗
USGS2024

SCEC/USGS Community Stress Drop Validation Study: How spectral fitting approaches influence measured source parameters

Spectral source parameters used to estimate an earthquake’s stress drop (Δσ) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress-Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low-frequency displacement (Ω0) and corner frequency ( f c ) by fitting a smooth function to benchmark displacement spectra. Participants used linear- or log-sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate Ω0, f c , and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in Ω0 and f c can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade-off between Ω0 and f c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω0 and somewhat higher f c compared to those assuming a Brune model, although resulting Δσ estimates are similar. When compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ.

California↗
USGS2025

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

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

Conference Paper↗