Search USGS⌕ Search

SEARCH · Search USGS

Results for “Arizona Geological Society Digest”

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.

745 records · Page 4Linked to original sources
USGS2025

Update of the 2008 provisional Enhanced Geothermal Systems (EGS) assessment for the Great Basin, USA

In response to the Energy Act of 2020, the U.S. Geological Survey (USGS) is updating the Enhanced Geothermal Systems (EGS) resource assessment for the Great Basin, USA. The previous 2008 provisional assessment estimated how much electricity could be generated from EGS resources of the western United States using models of electric-grade heat, models of heat extraction over time, and estimates of how much rock might be stimulated to produce viable amounts of heat. Herein, a similar conceptual strategy is applied, using updated models of heat extraction as a function of fracture spacing and well distance. Previously used reservoir heat delivery models are updated to have a dependence on fracture and well spacing, potentially improving future estimates of EGS resources as ongoing research provides a better understanding about the success of reservoir stimulation as a function of geology and location. For a range of well distances (250-1000 m) and fracture spacings (1-50 m), heat extraction efficiency ranges from 25-62%, demonstrating the importance of accounting for the most likely results of proven viable fracturing technologies. Although fracturing is important, the biggest uncertainty by far in estimating the EGS resource for the Great Basin is estimating which geologic units at what depths can be stimulated sufficiently to produce geothermal energy economically and efficiently. Uncertainties in these factors yield estimates that range over two orders of magnitude with an upper limit of ~174 terawatts-thermal (TWth) produced for 30 years from the upper 7 km of the crust. This upper limit would require significant technological advances to access most of the electric-grade resource across the Great Basin. Assuming that 1% of this estimate will be accessible in the next few decades gives a resource estimate similar to that made in the 2008 provisional assessment. These estimated EGS heat extraction rates far exceed ( greater than 100x) the natural geothermal heat production rate, thus geothermal electricity production at these rates might not be sustainable unless heat is also recharged from other sources (e.g., excess solar energy when supply exceeds demand). In addition to assessment maps and cumulative estimates, the new models of fractured reservoirs developed herein can be used to estimate steady power production given a set of fractures and well spacing, and estimates can be made for setback distances to ensure no thermal interference with nearby powerplants.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2025

“Leaky weirs” capture alluvial deposition and enhance seasonal mountain-front recharge in dryland streams

“Leaky weirs” are rock structures installed in dryland streams, which are anchored into exposed bedrock, loosely cemented, and designed to allow water to slowly pass through. They are being tested at a ranch in southeastern Arizona, USA, to restore and conserve the historic range and desert wetlands. Data are collected to assess how leaky weirs impact surface water, subsurface water, and groundwater recharge—including stream discharge, timing, and depth of infiltration, and groundwater elevations. Three adjacent watersheds, two with outlets just below leaky weirs and one with leaky weirs farther upstream, were instrumented with water-level loggers, wildlife cameras, and crest stage instruments with temperature sensors in the soil. As most groundwater recharge is assumed to be focused along the mountain fronts in this region, mountain-block recharge is also evaluated to differentiate between the two using isotope analyses. Finally, a single, late-season flood event is scrutinized to consider the leaky weir effect on all monitored components in the water budget. Results indicated groundwater flow is primarily from the mountains to the east via older, regional mountain-block recharge. However, the development of shallow alluvial aquifers is supported by the leaky weirs, that slow flows, capture permeable sediments, and allow infiltration, thus enhancing mountain-front recharge. In turn, these new pockets of water help support the restoration of historic wetlands. Sediment accumulates where leaky weirs are installed, reducing flashy peak flows, and resulting in a series of infiltration ponds along the channel that support vegetation during growing seasons and recharge the shallow aquifer during non-growing seasons.

Arizona↗
USGS2025

A scientific vision and roadmap for earthquake rupture forecast developments, a USGS perspective

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which is particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, un-quantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the United States Geological Survey and the external community. The primary goal of this paper is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early career scientists and engineers. The paper is written so readers can focus on the sections that interest them most (see table of contents), with the Introduction and Discussion providing a stand-alone overview and summary.

Bulletin of the Seismological Society of America↗
USGS2026

Ungulate migrations of the Western United States, volume 6

This report, volume 6 in the “Ungulate Migrations of the Western United States” report series, showcases the migrations of 23 ungulate herds in the Western United States. The report series is produced by the Corridor Mapping Team (CMT). Led by the U.S. Geological Survey, the CMT is a collaboration among 11 State agencies, as well as regional and Federal partners, and an expanding number of Tribal wildlife agencies. The CMT was initiated in response to the U.S. Department of the Interior Secretarial Order 3362, which was signed in 2018 and provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations throughout the Western United States. Including this volume, the report series has detailed the migrations of 237 unique ungulate herds throughout the Western United States and continues to serve as a valuable resource to guide local and regional management, policy, and on-the-ground work necessary to maintain intact and functional ungulate migrations. This report highlights several guiding principles of the CMT that facilitate collaboration among the diverse set of partners and contribute to the program’s continued successes. Notably, raw global positioning system data are not shared among participating agencies and the U.S. Geological Survey, delineating migration corridors and seasonal ranges relies on empirical data, the CMT provides flexible approaches to participating State and Tribal partners, and regular CMT meetings create a framework for open communication among agency partners that supports transboundary mapping of migrations. The 237 ungulate migrations that have been included in the report series are an expanding inventory, which can help maintain ungulate migrations in perpetuity.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS2024

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

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

Bulletin of the Seismological Society of America↗
USGS2025

Ungulate migrations of the Western United States, volume 5

Many ungulates migrate between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. In 2018, the U.S. Department of the Interior established Secretarial Order 3362, which provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations by enhancing habitat quality for ungulates across the Western United States. In response to the order, the U.S. Geological Survey (USGS) created the Corridor Mapping Team, which is a collaboration among 11 State agencies, regional and Federal partners, and an expanding number of Tribal wildlife agencies. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States and publishes them in the USGS “Ungulate Migrations of the Western United States” report series. This report details migrations and seasonal ranges from 36 additional herds and includes 2 herd updates detailed in previous reports. The Corridor Mapping Team has mapped the migrations and seasonal ranges of 218 unique herds for the report series, including this report. The report series serves as a map-based inventory of the ungulate migrations across the Western United States for biologists, managers, policymakers, and conservation practitioners. Building on the previous report volumes in the series, volume 5 additionally describes some of the local and national initiatives that are incorporating the products, tools, and information from this growing USGS report series.

Navajo Nation↗
USGS2024

Slip rate for the Rose Canyon fault through San Diego, California, based on analysis of GPS data: Evidence for a potential Rose Canyon–San Miguel-Vallecitos fault connection?

The Rose Canyon fault is the southern extension of the larger Newport–Inglewood–Rose Canyon fault system, which represents a major structural boundary in the Inner Continental Borderland (ICB) offshore of southern California. Ten to fifteen percent of total plate boundary motion in southern California is thought to be accommodated by the faults of the ICB, but the exact distribution of slip is uncertain. With an onshore segment, the Rose Canyon fault offers an opportunity to measure the slip rate using traditional geodetic methods. In this study, we use Global Positioning System (GPS) surface velocities from a combined campaign and continuous GPS network to constrain elastic models of the Rose Canyon fault. We then compare the observed surface velocities with proposed conceptual models of regional fault connections that facilitate the transfer of slip into the Rose Canyon fault to assess how well the observations are explained by the models. The results of elastic half‐space models suggest that the Rose Canyon fault may be slipping toward the higher end of geologic estimates, with the preferred model indicating a slip rate of 2.4 ± 0.5 mm/yr. Although limited in terms of near‐fault benchmarks, we find an improved model fit using an asymmetrical elastic half‐space model and a higher slip rate, suggesting a potential rheological contrast across the Rose Canyon fault, similar to observations from the northern Newport–Inglewood fault segments. Observed GPS surface velocities, background seismicity, and gravity anomalies south of San Diego Bay point toward a more easterly trace for the Rose Canyon fault, suggesting a possible connection with the San Miguel–Vallecitos fault system. Such a connection could increase the potential rupture lengths of future earthquakes and have important consequences for regional seismic hazards.

California↗
USGS2025

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

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

British Columbia, California, Oregon, Washington↗
USGS2024

Precariously balanced rocks in northern New York and Vermont, U.S.A.: Ground-motion constraints and implications for fault sources

Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake‐source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground‐based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time‐independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site‐specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground‐motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M ∼6.5 earthquake‐source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake‐source faults in regions of cryptic faults.

Bulletin of the Seismological Society of America↗
USGS2026

Resurvey of the Marble Canyon and Bridge Canyon dam sites in Grand Canyon National Park—Changes in sediment storage and evidence supporting the occurrence of bedrock incision through the mid-20th century

The Bureau of Reclamation developed an extensive plan for a network of dams, water tunnels, and hydropower plants in and around Grand Canyon, Arizona, in the 1940s through 1960s. The two largest of these planned dams were the Marble Canyon and Bridge Canyon Dams on the Colorado River. Though these dams were ultimately never built, Reclamation conducted extensive topographic, bathymetric, and subsurface exploration work at the sites proposed for these dams in the 1940s and 1950s. Resurveys of these dam sites were conducted between 1998 and 2021 to determine the changes in sediment storage at these dam sites caused by the upstream construction and operation of Glen Canyon Dam and by the recession of Lake Mead, the reservoir impounded by Hoover Dam. The resurveys of the Marble Canyon dam sites indicate that the post-1950s changes in sediment storage at these dam sites are broadly consistent with flux-based estimates of voluminous sand erosion from Marble Canyon since the 1963 closure of Glen Canyon Dam. These resurveys also suggest that the pre-dam longitudinal variation in sediment thickness over bedrock played a key role in determining the locations of the sand erosion induced by Glen Canyon Dam; more sand eroded from locations where more sand was present in the 1950s. The resurvey of the Bridge Canyon dam sites indicates that the Colorado River’s incision of the Lake Mead delta is regulated both by bed-sediment grain size and downstream hydraulic controls. Finally, analyses of bed-sediment thickness and sedimentological data at the dam sites, and observations of bed scour and gravel transport, suggest that sufficient bedrock was exposed to allow bedrock incision during commonly recurring pre-dam snowmelt floods that entrained small boulders into transport.

Arizona↗
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

New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States

Four new 10 Be- 26 Al isochron burial ages ranging from 4.4 to 2 Ma on ancestral Colorado River deposits in the lower Colorado River corridor (LCRC) help constrain the river’s evolution during the Pliocene and early Pleistocene. They help fill a gap between previous work that focused on older and younger deposits: Older dated deposits include the 5 Ma Bouse Formation, which records the integration of the Colorado River through a series of preexisting basins to the Gulf of California and the ca. 4.5 - 3.5 Ma Bullhead Alluvium, a 200 to 300 m thick aggregational package that immediately followed integration. The much younger, 100 - 70 ka, Chemehuevi Formation is another major aggradation package mapped throughout the LCRC. The new burial ages on the facies of Santa Fe Railway (4.37 ± 0.71 Ma), boulder conglomerate of Bat Cave Wash (2.12 ± 0.26 and 2.05 ± 0.31 Ma), and the Palo Verde alluvium (3.03 ± 0.26 Ma) partially fill in a 3.5 M.y. gap between the deposition of the Bullhead Alluvium and the Chemeheuvi Formation and document the timescales over which the Colorado River was able to remove the Bullhead aggradational package and initiate newer and smaller aggradational pulses.

Arizona, California, Nevada↗
USGS2024

The impact of source time function complexity on stress drop estimates

Earthquake stress drop—a key parameter for describing the energetics of earthquake rupture—can be estimated in several different, but theoretically equivalent, ways. However, independent estimates for the same earthquakes sometimes differ significantly. We find that earthquake source complexity plays a significant role in why theoretically (for simple rupture models) equivalent methods produce different estimates. We apply time‐ and frequency‐domain methods to estimate stress drops for real earthquakes in the SCARDEC (Seismic source ChAracteristics Retrieved from DEConvolving teleseismic body waves, Vallée and Douet, 2016 ) source time function (STF) database and analyze how rupture complexity drives stress‐drop estimate discrepancies. Specifically, we identify two complexity metrics—Brune relative energy (BRE) and spectral decay—that parameterize an earthquake’s complexity relative to the standard Brune model and strongly correlate with the estimate discrepancies. We find that the observed systematic magnitude–stress‐drop trends may reflect underlying changes in STF complexity, not necessarily trends in actual stress drop. Both the decay and BRE parameters vary systematically with magnitude, but whether this magnitude–complexity relationship is real remains unresolved.

Bulletin of the Seismological Society of America↗
USGS2025

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

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

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

A comprehensive fault system inversion approach: Methods and application to NSHM23

We present updated inversion‐based fault‐system solutions for the 2023 update to the National Seismic Hazard Model (NSHM23), standardizing earthquake rate model calculations on crustal faults across the western United States. We build upon the inversion methodology used in the Third Uniform California Earthquake Rupture Forecast (UCERF3) to solve for time‐independent rates of earthquakes in an interconnected fault system. The updated model explicitly maps out a wide range of fault recurrence and segmentation behavior (epistemic uncertainty), more completely exploring the solution space of viable models beyond those of UCERF3. We also improve the simulated annealing implementation, greatly increasing computational efficiency (and thus inversion convergence), and introduce an adaptive constraint weight calculation algorithm that helps to mediate between competing constraints. Hazard calculations show that ingredient changes (especially fault and deformation models) are the primary driver of hazard changes between NSHM23 and UCERF3. Updates to the inversion methodology are also consequential near faults in which the slip rate in UCERF3 was poorly fit or was satisfied primarily using large multifault ruptures that are now restricted by explicit b ‐value and segmentation constraints.

Bulletin of the Seismological Society of America↗
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

Introduction to the special section on improving measurements of earthquake source parameters

Earthquake source parameters such as magnitude, seismic moment, source dimension, stress drop, and radiated energy are fundamental to understanding earthquake physics, and are also key ingredients in earthquake ground‐motion modeling, rupture simulation, and statistical seismology. However, the uncertainties in these parameters estimated from the radiated seismic wavefield are large due to variability in approaches, including site and attenuation characterizations, and so estimates for an individual earthquake made by different studies can vary greatly. Estimating spectral source parameters remains a popular topic, due to a combination of their intrinsic importance and their apparent ease of measurement, but also a controversial one, due to many sources of variability and large uncertainty. The available methods coupled with necessary parameter choices and assumptions in the analysis make it challenging for researchers to apply methods or understand the reliability in results or reported source parameters. This Special Section on Improving Measurements of Earthquake Source Parameters showcasing comparisons between methods and studies seeks to alleviate some of these difficulties to help the community identify the important components and trade‐offs of decomposing recorded seismograms into their source, path, and site components.

Bulletin of the Seismological Society of America↗