USGS ScienceSearch

SEARCH · USGS Science

Results for “The Seismic Record”

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,669 records · Page 25Linked to original sources

Preliminary ground and airborne-based geophysical mapping and modelling of an active hydrothermal system at Mammoth Lakes, California

Mammoth Lakes, California hosts a productive hydrothermal system within the seismically active south moat of Long Valley Caldera. Surficial evidence of the shallow hydrothermal system includes discrete zones of tree-kill dispersed between Shady Rest Park and the Casa Diablo Geothermal Power Plant (40 MW), as well as east of the power plant. The tree-kill areas are associated with elevated diffuse CO₂ emissions, heated ground, hydrothermal alteration, diffuse soil H₂S emissions, and gas vents. Previous mapping delineates prominent north and northwest trending structures within the south moat along the southwestern edge of the resurgent dome that may accommodate gas and fluid flow at the Shady Rest Park and Basalt Canyon Tree Kill Areas (SRTKA and BCTKA, respectively). Both tree-kill areas are also located along contacts between resurgent rhyolite, mafic lavas, and surficial deposits which may provide additional pathways for gas and fluid migration in the shallow subsurface. Characterizing structure and lithology using geophysical anomalies is critical to determining primary structural controls on the hydrothermal system and the extent of subsurface alteration at these sites. We conducted ground and airborne-based potential field geophysical surveys to map gravity and magnetic anomalies. These anomalies are then used to model subsurface geology, structure, and hydrothermal alteration. Here we present our preliminary geophysical mapping and modelling results at both tree-kill locations. Gravity and magnetic data suggest complex structural intersections are coincident with heated ground and gas emissions at the SRTKA and BCTKA. Hydrothermal systems are often observed or interpreted to exploit fault intersections which can serve as highly permeable pathways for hydrothermal fluid and gas discharge, enabling economic geothermal energy production. Geophysical mapping and modelling are an effective means of investigating such structural complexity at Mammoth Lakes due to the presence of unidentified and concealed structures.

California

Refinements to the attenuated propagation of local earthquake shaking (APPLES) ground-motion-based earthquake early warning algorithm

We refined the Attenuated ProPagation of Local Earthquake Shaking (APPLES) ground-motion-based earthquake early warning (EEW) approach, and directly compare APPLES performance with that of the source-characterization-based U.S. ShakeAlert EEW system for a suite of historical earthquakes in the U.S. West Coast and Japan. APPLES is an extension of the Propagation of Local Undamped Motion (PLUM) algorithm in which observed shaking intensity at seismic stations is used to forward-predict intensity distributions to surrounding areas using an attenuation model derived from an intensity prediction equation. We test new configuration options within APPLES, such as using the second highest estimated ground motion rather than the maximum, to better match median ground-motion observations and reduce alerts for small magnitude earthquakes, both of which are key alerting priorities within ShakeAlert. We evaluate these configurations alongside ShakeAlert by comparing the ground-motion estimation accuracy and available warning times relative to station observations and ShakeMap distributions. Our preferred APPLES configuration produces accurate ground-motion estimates and corresponds better with median observations compared to ShakeAlert’s estimates. This preferred configuration substantially reduces alert issuance for M < 5.0 earthquakes compared to the previous APPLES configuration, and alert-release criteria can further restrict alerts to primarily M ≥ 5.5 earthquakes without requiring magnitude estimation. Prioritizing matching median-observed ground motions may reduce APPLES warning times compared to configurations that were tuned to avoid missed alerts (such as those that use the maximum estimated ground motions), which can lead to shorter warning times compared to ShakeAlert for the same alert threshold. However, station-based warning time assessments demonstrate that APPLES can outperform ShakeAlert for high target thresholds. APPLES is a simple, independent EEW approach that may improve the robustness of EEW for the West Coast of the U.S.

California, Oregon, Washington

Site response and wave propagation effects in the eastern United States

Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.

eastern United States

Corwin Springs Known Geothermal Resources Area, Park County, Montana

The Corwin Springs Known Geothermal Resources Area (KGRA) is contiguous to Yellowstone National Park along a part of the northern boundary of the park near Gardiner, Park County, Mont. The area contains two known sites of hot-spring activity--LaDuke Spring, 2.8 km southeast of the small resort community of Corwin Springs, and Bear Creek Spring, 2.6 km east of Gardiner. LaDuke Spring issues from brecciated quartzite and has a flow rate of 380 1/min, a surface water temperature of 65°C, a silica-geothermometer temperature of 68.7°C, and a Na-K-Ca geothermometer temperature of 76.8°C. Bear Creek Spring issues from limestone and has a flow rate of 4 1/min, a surface water temperature of 32°C, a silica temperature of 46.7°C, and a Na-K-Ca temperature of 87.2°C. The springs, which are actively depositing travertine, are on or near the trace of the Gardiner fault-a high-angle reverse fault which forms the southwestern boundary of the Beartooth uplift. Features that make the Corwin Springs KGRA potentially significant include: (1) Proximity of the Corwin Springs area to significant Pliocene and Pleistocene volcanism in Yellowstone National Park and related geothermal activity, (2) existence of a large potential geothermal reservoir in the northern extension of the Sepulcher Mountain graben, (3) localization of thermal activity along the Gardiner fault in the Corwin Springs area, (4) occurrence of negative gravity anomaly centered over the northern end of the Sepulcher Mountain graben, (5) existence of negative magnetic anomaly in the Corwin Springs area, (6) location of the area within the Intermountain seismic belt, (7) evidence of recent tectonism in the area, as demonstrated by Pleistocene and Holocene faulting, and (8) observed surface temperatures and estimated geochemical temperatures of the two known hot springs (features which do not indicate high subsurface temperatures).

Montana

Risk implications of Poisson assumptions and declustering inferred from a fully time-dependent earthquake forecast

We use the Third Uniform California Earthquake Rupture Forecast Epidemic Type Aftershock Sequence model, which is fully time-dependent in terms of including spatiotemporal clustering, to evaluate the effects of the Poisson assumption and declustering algorithms on statewide loss exceedance curves. The model is simulation based, meaning it produces synthetic catalogs that exhibit realistic behavior with respect to aftershocks and multi-fault earthquakes. A Poisson version of the model was constructed by randomizing event times, and the influence of two declustering algorithms was examined as well. We demonstrate that the probability of one-or-more loss exceedances (occurrence exceedance probability) is greater for the Poisson model because it has fewer seismically quiet time windows. The discrepancy between dollar loss estimates with a given exceedance probability is up to a factor of 32% but varies depending on the loss threshold (the x-axis value) and the forecast duration (we examined a range between 24 h and 50 years, with the discrepancy for the latter being negligible). We discuss how the one-or-more loss exceedance metric is questionable because it ignores all but the maximum loss experienced in each timeframe. An alternative metric based on total aggregate loss in each time window (aggregate exceedance probability) was therefore also examined, for which the Poisson model again implies higher risk at intermediate losses but lower risk at higher losses (because large, triggered events now contribute to total aggregate losses for the fully time-dependent model). We also argue that declustering is not a scientifically justifiable way to deal with full time dependence, in agreement with a chorus from other recent studies. It is difficult to draw generally applicable conclusions from our study, in part because application specific details will likely be important, but our results highlight how full time dependence can be reckoned with once authoritative forecast models are made available.

California

Streams, springs, and volcanic lakes for volcano monitoring

Introduction Volcanic unrest can trigger appreciable change to surface waters such as streams, springs, and volcanic lakes. Magma degassing produces gases and soluble salts that are absorbed into groundwater that feeds streams and lakes. As magma ascends, the amount of heat and degassing will increase, and so will any related geochemical and thermal signal. Subsurface magma movement can cause pressurization that alters hydrostatic head and may induce groundwater discharge. Fluid-pressure changes have been linked to distal volcano-tectonic earthquakes (White and McCausland, 2016; Coulon and others, 2017) and phreatic eruptions (for example, Yamaoka and others, 2016). Clearly, changes in groundwater and surface waters are both indicators of unrest and clues to how and where magma is rising toward the surface. Where possible, it is prudent to incorporate real-time hydrologic data into multiparameter monitoring of restless volcanoes. Hydrologic dynamics can also be tracked by changes in groundwater levels that are commonly measured in shallow boreholes (see of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume, on boreholes; Hurwitz and Lowenstern, 2024). Although inferred to be common, relatively few volcano-hydrology anomalies are well documented, and many are essentially anecdotal (Newhall and others, 2001), reflecting the fact that high-resolution time series remain rare. Extreme examples include the 2008 eruption of Nevado del Huila, Colombia, where relatively minor phreatomagmatic eruptions were accompanied by expulsion of as much as 300 million cubic meters of groundwater from fissures high on the volcano (Worni and others, 2011), generating large lahars. Substantial decreases in flow rate from springs about 8 kilometers from the summit of Mayon Volcano, Philippines, have been noted before most eruptions in the 20th century (Newhall and others, 2001). Stream monitoring at Redoubt Volcano in 2009 allowed Werner and others (2012) to recognize that groundwater was unable to absorb (or scrub) the high flux of volcanic gas and that a high CO 2 /SO 2 precursor signal had been evident for 5 months prior to the eruption. A key to better interpreting hydrologic anomalies—or even identifying them—is therefore obtaining adequate baseline data. Most hydrologic monitoring at U.S. volcanoes has been accomplished by intermittent sampling surveys with annual or less frequent sampling (for example, https://hotspringchem.wr.usgs.gov/index.php ). More frequent sampling, however, generally is needed to establish reliable baselines. A recent hydrologic and hydrothermal monitoring experiment at 25 sites and 10 of the 12 level 4 (very high threat) volcanoes in the U.S. portion of the Cascade Range demonstrated that there is sufficient temporal variability in hydrothermal fluxes, even during quiescent periods, that one-time measurements will commonly have limited interpretive value (Crankshaw and others, 2018). Thus, surveys are best augmented with data from streamgages (for example, Evans and others, 2004; Bergfeld and others, 2008). Streamflow (water discharge) data allow measured temperature and specific conductance to be converted to heat and solute mass fluxes, which could be insightful parameters for detecting anomalous activity (McCleskey and others, 2012). At the Yellowstone Caldera, long-term monitoring of river solutes has allowed calculation of the chloride flux, a proxy for heat discharge (Hurwitz and others, 2007; McCleskey and others, 2016) from the subsurface magma. This is readily accomplished because data from streamgages are continuously recorded and archived by the U.S. Geological Survey (USGS) National Water Information System (NWIS) (USGS, 2024). Similar studies on stratovolcanoes or shield volcanoes would be scientifically useful, and yet are logistically challenging, requiring streamgages on numerous radial drainages complemented by either frequent manual sampling or numerous deployments of equipment to measure water temperature and specific conductance as a proxy for water chemistry. Another challenge is that some volcanic areas, especially shield volcanoes, are characterized by near-surface porous rocks and soils, such that surface streams are rare and replaced by distant, dilute large-volume springs with only a trace of any original volcanically sourced water (Manga, 2001; Hurwitz and others, 2021). Volcanic lakes are worthy of special attention for monitoring efforts, as their temperature and composition can provide evidence of increased flux of volatile-rich fluids from below. Quantifying changes in volatile and heat release from magma can be simpler in lakes than for volcanoes with radial drainages and no major lakes. Moreover, volcanic lakes pose a range of hazards themselves, including phreatomagmatic eruptions, debris flows, flank collapse, tsunamis, and toxic gas release (Mastin and Witter, 2000; Delmelle and others, 2015; Manville, 2015; Rouwet and others, 2015)—hazards that have historically been responsible for substantial loss of life at many volcanoes worldwide (Manville, 2015). Catastrophic CO 2 release at Lake Nyos, Cameroon, in 1986 suffocated about 1,750 people and about 3,500 livestock and was probably triggered by a large landslide into the gas-saturated lake (Kling and others, 1987; Evans and others, 1993). Gas-charged springs in Soda Bay within Clear Lake (California) have caused almost a dozen deaths to bathers in the past hundred years (ABC News, 2000). A 2005 example of lake overturn and abundant gas release was documented at Mount Chiginagak in Alaska (Schaefer and others, 2008) but did not result in any human casualties. Although thermally stratified lakes, which promote trapping of exsolved magmatic gas, tend to develop in tropical regions, the phenomenon can also arise where salinity creates meromixis (a condition in which a lake does not mix completely), as occurs in Mono Lake, California (Jellison and Melack, 1993; Jellison and others, 1998). If magma erupts or flows into a lake, the interaction between hot magma and cold water can be explosive (Mastin and others, 2004; Zimanowski and others, 2015) and substantially expand the area affected by the eruption. Another hazard is the breaching of crater rims by landslides triggered by volcanic and (or) seismic activity. Under some circumstances, substantial volumes of water can be displaced, leading to large floods and lahars. Late Holocene lake flooding from Aniakchak Crater in the Alaska Peninsula (Waythomas, 2022) and from Paulina Lake in Newberry Crater, Oregon (Chitwood and Jensen, 2000), caused by the failure of outlet sills, testify to the substantial hazards at lake-filled calderas. Several volcanic systems in the United States host lakes known to receive heat and gas from underlying magma. These lakes vary widely in area, depth, and chemical composition. Lakes are present at level 4 volcanoes, including Crater Lake and Newberry Volcano in Oregon; Yellowstone Caldera in Wyoming; Long Valley Caldera, Clear Lake volcanic field, Medicine Lake, and Salton Buttes in California; and Aniakchak Crater, Mount Katmai, Fisher Caldera, Mount Okmok, and Kaguyak Crater, among others, in Alaska. A water lake was present in Halemaʻumaʻu, the crater of Kīlauea, Hawai‘i (fig. F1), from October 2019 to December 2020. Level 3 volcanoes with lakes include Mono Lake volcanic field (Calif.), Mount Bachelor (Ore.), Ukinrek Maars and Mount Chiginagak (Alaska), and Soda Lake (Nevada). In addition, there are lakes at many levels 1 and 2 volcanoes. In the United States, there are no strongly acidic lakes that receive abundant input of magmatic gas, such as those found at Mount Ruapehu (New Zealand), Ijen and Kelud (Indonesia), and Poás (Costa Rica). Nevertheless, many contain fluids that provide clues to magmatic processes below. Since publication of a previous report on recommended instrumentation for volcano monitoring (Moran and others, 2008), continuous hydrologic monitoring has become increasingly feasible. However, changes in water pressure, temperature, and chemistry remain, in general, poorly studied phenomena at volcanoes (Sparks, 2003; National Academies of Sciences, Engineering, and Medicine, 2017). Recent efforts by the USGS have included the temporary study of Cascade Range volcanoes, which included frequent (15 minute to hourly) temporal sampling of temperature, depth, and conductivity (Crankshaw and others, 2018; Ingebritsen and Evans, 2019). At Yellowstone Caldera, many streamgages have now added thermistors and specific conductance sensors, allowing estimation of time-dependent chloride flux as a proxy for variations in subsurface heat flux (McCleskey and others, 2012, 2016). Efforts to better understand lakes have also accelerated, with bathymetric mapping and sampling carried out at several locations in the United States. Especially thorough work was done at Yellowstone Lake thanks to the Hydrothermal Dynamics of Yellowstone Lake (HD-YLAKE, https://hdylake.org ) project, funded primarily by the National Science Foundation. In addition to geophysical surveys and recovery of cores and other samples, HD-YLAKE investigations included remotely operated vehicle (ROV) investigations of hydrothermal vents on the lake floor (fig. F2). Data collected by the ROV provided a better understanding of the thermal and chemical influx from lake-bottom hydrothermal systems (Sohn and others, 2017). In this chapter, we focus on detecting changes in the chemistry, temperature, discharge, or water levels of streams, springs, and lakes that can be caused by seismicity, volumetric strains, or increases in gas flux associated with ascending magma. There is unavoidable overlap with other chapters of this report. Samples of water and gas can also be obtained in boreholes ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–K, 5 p., https://doi.org/10.3133/sir20245062k. ">chapter K , this volume; Hurwitz and Lowenstern, 2024), both shallow and deep. Gas monitoring ( of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–E, 11 p., https://doi.org/10.3133/sir20245062e.">chapter E , this volume; Lewicki and others, 2024) relies in part on samples from springs and wells, particularly where measurable gas plumes are absent. Water acts as a trigger and lubricant for landslides and sediment-rich floods, and so hydrology has obvious relevance for lahar monitoring, as discussed in of Flinders, A.F., Lowenstern, J.B., Coombs, M.L., and Poland, M.P., eds., Recommended capabilities and instrumentation for volcano monitoring in the United States: U.S. Geological Survey Scientific Investigations Report 2024–5062–H, 6 p., https://doi.org/10.3133/sir20245062h. ">chapter H (this volume; Thelen and others, 2024). Shared situational awareness among scientists engaged in geophysical, gas, and hydrologic monitoring will improve overall understanding of the volcanic hazard.

Scientific Investigations Report

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay

Ten new insights in climate science 2025

Interdisciplinary understanding is vital for delivering sound climate policy advice. However, navigating the ever-growing and increasingly diverse scholarly literature on climate change is challenging for any individual researcher. This annual synthesis highlights and explains recent advances across a variety of fields of climate change research. This year, the 10 insights focus on: (1) the record-warmth of 2023/2024 and the elevated Earth energy imbalance; (2) acceleration of ocean warming and intensifying marine heatwaves; (3) northern land carbon sinks under strain; (4) reinforcing feedback between biodiversity loss and climate change; (5) accelerated depletion of groundwater; (6) global dengue incidence; (7) global income losses and labour productivity declines; (8) strategic scaling of CDR; (9) integrity challenges in carbon credit markets and emerging responses; and (10) effective policy mixes for emissions reductions. The insights have been written to be accessible to researchers from different fields, serving as entry-points to specific topics, as well as providing an overview of the evolving landscape of climate change research. In the final section, the insights are used to develop overarching policy-relevant messages. This paper provides the basis for a science-policy report that was shared with all Party delegations ahead of COP30 in Belém, Brazil.

Global Sustainability

Geologic models underpinning the 2024 U.S. Geological Survey assessment of undiscovered oil and gas resources in the Hosston and Travis Peak Formations of the onshore Gulf Coast region, U.S.A.

The Early Cretaceous (Berriasian–Hauterivian) Hosston Formation in Louisiana and eastward is time correlative to the Travis Peak Formation of Texas and southern Arkansas. The formation is a first-order clastic sequence with a regional carbonate transgressive surface as an upper contact. The Hosston and Travis Peak formations contain conventional natural gas and oil accumulations that have been produced for nearly a century. These mature reservoirs contain terrigenous fluvial-deltaic, shore-zone, and paralic deposits across the productive trend; organic-lean mudstone and siltstone lithologies are found outboard of the Lower Cretaceous shelf margin. Producing reservoirs exhibit normal pressure gradients from 0.43 to 0.55 psi/ft (9.7 to 12.4 kpa/m), depths from 4000 to over 20,000 ft (1220 to 6100 m), and temperatures from 150 to 385°F (65 to 196°C). Wells are primarily vertical completions. The number of new field wildcats has been declining since the late 1990s. This paper presents comprehensive geologic models, which include lithofacies maps, structure and isopach maps, burial history models, regional seismic interpretations, and events charts that underpin the recently completed U.S. Geological Survey assessment of undiscovered, technically recoverable hydrocarbons within the Hosston and Travis Peak formations. This study also provides geographic and stratigraphic distributions of Hosston–Travis Peak reservoir properties, including geopressure, reservoir temperature, porosity, permeability, API gravity, and gas-oil ratios. Results indicate estimated undiscovered, technically recoverable mean resources of 28 million barrels of oil and 35.8 trillion cubic ft of gas in conventional and continuous accumulations within the Lower Cretaceous Hosston and Travis Peak formations of the onshore U.S. Gulf Coast region. Quantitative assessment results are detailed in U.S. Geological Survey Fact Sheet 2025–3021 and associated Data Release.

GCAGS Journal

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California

Insights on gas hydrate formation and growth within an interbedded sand reservoir from well logging at the Qiongdongnan Basin, South China Sea

Although variable well log resolution and its control on saturation estimation has been studied, it has not been directly applied to a specific location to explore the nature of gas hydrate within a sand reservoir. We applied in-situ measurements of resistivities, neutron porosity, and gamma ray at two sites in the Qiongdongnan Basin, South China Sea (QDN-W05–2021 and QDN-W08–2021) to investigate the reservoir parameters of a hydrate-bearing sand reservoir. Our results show that gas hydrate is distributed in 5 zones with a total thickness of 10.7 m and an average saturation of 69% at the QDN-W05–2021 site, while they are distributed in 2 zones with a total thickness of 4.3 m and an average saturation of 49% at the QDN-W08–2021 site. We found that variances in saturations estimated from lateral-extra deep button (RX), phase shift (P40H-P40L), and attenuation (A40H-A40L) resistivities within the laterally mapped continuous sand body were affected by the nature of gas hydrate occurrences. Results indicate gas hydrate forms and accumulates at the center of the sand layer and tends to be less or not present toward the top and base. Integrated with seismic data, the in-situ measurements provide insights in the evolution of a mushroom-shaped, hydrate-gas reservoir system. In the system, free gas is likely horizontally transported from the top-center of the gas chimney to the surrounding areas in the early stage dominated by a warm-gas environment, whereas hydrate forms in the opposite pathway starting from the surrounding areas in the following stage with temperature reducing. Our study suggests that high-resolution in-situ measurements not only are a tool to identify the physical properties, but also can be used to help explain the physical process of hydrate growth and accumulation.

Marine Geology

Analysis of factors affecting plume remediation in a sole-source aquifer system, southeastern Nassau County, New York

Several plumes of dissolved, chlorinated solvents, including trichloroethylene, have been identified in a sole-source aquifer near the former Northrop Grumman Bethpage Facility and Naval Weapons Industrial Reserve Plant sites in southeastern Nassau County, New York. Past investigations have documented that the groundwater contamination originated from this industrial area and now extends to the south, in the direction of groundwater flow. The intermixed plumes are commonly referred to as the “Navy Grumman groundwater plume.” Detailed groundwater-flow modeling was needed for the New York State Department of Environmental Conservation (NYSDEC) to evaluate design options necessary for the construction, operation, optimization, maintenance, and monitoring of a groundwater extraction and treatment cleanup plan selected in a December 2019 Amended Record of Decision by the NYSDEC to comprehensively address these plumes. Consequently, the NYSDEC began a cooperative study with the U.S. Geological Survey in 2020 to better understand the local hydrogeologic framework using two independent approaches to characterize aquifer heterogeneity and update an existing regional groundwater-flow model to provide transient boundary conditions for new inset groundwater-flow models of the plume area. We developed these detailed inset models for the two independent aquifer characterizations using history-matching techniques coupled with a novel approach to risk-based management optimization of the remedial design. We also used the updated regional model to assess this optimized groundwater extraction and treatment design for potential saltwater intrusion. The ensembles of parameters resulting from history matching provided a platform with which to evaluate capture by water-supply and remedial wells using particle-tracking techniques. Using the ensemble to select a risk stance, we performed multiobjective optimization to identify various configurations of remedial pumping that are consistent with external constraints and that favor potentially competing objectives. Multiple solutions provide tradeoffs that NYSDEC can consider. In general, pumping redistribution may help to prevent further contamination migration downgradient. These and other study results are intended to support decisions for the remedial design focused on the local area encompassing the full extent of the Navy Grumman groundwater plume.

New York

Unprecedented burning in tropical peatlands during the 20th century compared to the previous two millennia

Tropical peatland wildfire incidence has risen in recent decades, driven by drainage for land use and intensified by severe droughts with global climate change. These disturbances have altered vegetation structure, disrupted ecosystem functioning, and increased carbon emissions, particularly in Southeast Asia. However, the long-term history and characteristics of wildfires in tropical peatlands remain largely unknown. Here, we compiled fifty-eight macro-charcoal records from peatlands across the tropics, ranging from lowland forested to montane peatlands, to assess millennia-scale changes and controlling factors of tropical peatland burning. We divided the datasets into four main sub-regions: Neotropical, Afrotropical, Indomalayan and Australasian ecoregions to explore regional variability. Tropical peatlands had high burning levels between 0 and 850 ce , followed by a relatively low and stable period until a marked increase during the 20th century. The general trend in tropical peatland burning follows changes in global temperature, and climate variables that control the length and severity of drought events have a notable influence on peat burning before 1900 ce . During the 20th century, regional differences were observed, with declining fire trends in the Neotropical and Afrotropical regions and increasing fire trends in the Indomalayan and Australasian regions. This difference is likely attributable to human activities, and such intervention is also evident in palm swamps and hardwood swamps under similar wet, weakly seasonal climates. With the increase in anthropogenic pressures on peatlands and greater climate variability, future wildfires in peatlands are likely to become more frequent and widespread across all tropical ecoregions. Conservation and sustainable land-use practices could be used to mitigate and control peatland burning and protect these carbon-rich sinks.

Global Change Biology

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