USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Investigations Map”

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.

862 records · Page 46Linked to original sources
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↗
USGS2026

Geochemical, mineralogical, and isotopic evidence for multi-stage genesis of the Hicks Dome REE + Y-HFSE-fluorite deposit, Illinois, USA

Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.

Illinois↗
USGS2024

Post-wildfire debris flows

Post-wildfire debris flows pose severe hazards to communities and infrastructure near and within recently burned mountainous terrain. Intense heat of wildfires changes the runoff characteristics of a watershed by combusting the vegetative canopy, litter, and duff, introducing ash into the soil and creating water repellant soils. Following wildfire, rainfall on bare ground is less able to infiltrate into the fire-altered soils and overland flow is less impeded by vegetation. Rainfall runoff in recently burned areas can erode hillslopes owing to the removal of soil binding organic matter near the soil surface by fire. In channels, loose, dry-ravel deposits composed of sand and gravel are readily entrained by concentrated runoff in channels. Entrainment of soil on hillslopes and in channels bulks up the sediment concentration of the rainfall runoff to generate debris flows capable of transporting boulders and large woody debris. Post-wildfire debris flows can be triggered by rainfall conditions that would typically produce little runoff during unburned conditions. The primary rainfall trigger for post-wildfire debris flows is high intensity rainfall during short duration convective rainstorms or periods of high rainfall intensity embedded within a long-duration frontal storm. Numerous observations of debris flows triggered by storms lasting less than an hour following periods of little to no rainfall indicate that antecedent rainfall is not a requirement for initiation of post-wildfire debris flows. Post-wildfire debris-flow hazard assessment entails estimating probability and magnitude of debris flows in the burned area, estimating debris-flow runout and intensity, and defining rainfall intensity-duration thresholds for debris-flow initiation. In the United States, probability and magnitude is estimated using empirically derived models largely based on data collected in southern California. The models provide maps to identify watersheds and drainage paths where post-wildfire hazards are most pronounced. Rainfall intensity-duration thresholds can be incorporated into flood hazard forecasting tools. Currently, work is underway to identify how to best implement debris-flow runout models in burned areas with efficiency and accuracy. Post-wildfire debris flows have been a long-recognized process in the Transverse Ranges of southern California; however, climate change is driving more frequent wildfires to burn more mountainous terrain throughout the western United States and worldwide. As a result, post-wildfire debris flows are becoming a more common threat in areas where they were once infrequent. As the threat of post-wildfire debris flow expands into new areas, evaluating the hazard becomes challenging because the degree to which wildfire increases debris-flow susceptibility varies from region to region. This chapter summarizes the knowledge to date for evaluating post-wildfire debris-flow susceptibility and hazard assessment. We summarize the characteristics of wildfire burn severity, topography, underlying soil and geology, and rainfall conditions that contribute to making a watershed most likely to produce post-wildfire debris flows. Methods for hazard assessment in the United States and other countries are summarized. We highlight knowledge gaps for how post-wildfire debris-flow susceptibility varies throughout the western United States and worldwide and identify research needs to improve hazard assessment methods in different geographies.

Book chapter↗
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗
USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗
USGS2025

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts↗
USGS1976

Water-quality study of a reach of the Merced River in Yosemite National Park and vicinity, California, April 1973 through September 1974

A 25-mile (40-kilometre) reach of the Merced River was studied to provide the National Park Service with water-quality data prior to the operation of a new sewage-treatment plant. Data were collected periodically at four Merced River and two inflow stations. Discharge varied seasonally with highest flow occurring at the time of spring snowmelt runoff. Lowest flow occurred during late summer, coinciding with high park visitor use and elevated water temperature. Water-quality variables investigated were discharge, temperature, major chemical constituents, selected trace metals, specific conductance, suspended sediment, turbidity, dissolved oxygen, pH, alkalinity, carbon dioxide, total organic carbon, nitrogen, phosphorus, algal growth potential, coliform bacteria, periphyton, and benthic invertebrates. Diel measurements were made at some stations to delineate daily fluctuations in selected water-quality variables. Results indicate that water in the reach of the Merced River studied was of good quality. The river had low concentrations of dissolved solids (8 to 41 milligrams per litre) and low alkalinity values (2 to 24 milligrams per litre of alkalinity as CaCO 3 ). Mean concentrations of nitrogen and phosphorus were highest in the Merced River at Rancheria Flat. Algal growth potential tests showed that phosphorus was limiting algal growth at Happy Isles Bridge, El Capitan Bridge, and South Fork Merced River, whereas nitrogen was limiting growth at Rancheria Flat. Input of excessive nitrogen and phosphorus concentrations in the Merced River during the summer low-flow period, and mixing with warm tributary inflow could cause excessive algal growth and production. The algal types found were diatoms, green algae, and blue-green algae. Biomass/chlorophyll a ratios in the Merced River were low (0.3 to 7.7). Benthic invertebrates consisted chiefly of Diptera, Ephemeroptera, and Plecoptera.

California↗
USGS2026

Assessing environmental drivers of denitrification in restored riverine floodplains

Restoration of impaired floodplains is an increasingly prevalent strategy for alleviating water quality concerns and reducing downstream flooding at watershed scales. Floodplains temporarily store water and slow flow velocity to promote sedimentation during overbank flooding and remove inorganic nitrogen from floodwater and groundwater via denitrification. Evaluating the impacts of different restoration strategies on denitrification can inform more strategic investments into floodplain modifications that improve water quality outcomes. Our research investigates how denitrification rates in floodplains respond to environmental factors that are actionable from an engineering perspective through design and water resources management. We seasonally measured soil denitrification enzyme activity and various environmental characteristics in 4 floodplains with different restoration design and management approaches at the confluence of the Wabash and Tippecanoe Rivers in Indiana, United States. Our results showed that denitrification rates in an agricultural floodplain were significantly lower than in restored floodplains with native vegetation. Certain soil conditions characteristic of floodplain wetlands were associated with higher denitrification, particularly elevated total nitrogen, moisture, silt, and organic matter contents. Vegetation species composition was correlated with denitrification rates. This link may reflect the direct effects of vegetation on soil conditions, such as supplying labile organic carbon, or indirect effects, such as vegetation acting as an indicator of hydrologic regime and land use. Denitrification seasonally varied, peaking in winter when nitrate supply from rivers draining agricultural watersheds in the region is also high. Substrate limitation of soil denitrification enzyme activity was most significant during the summer when overbank flooding, which replenishes soil nitrogen stocks, rarely occurs. Our findings indicate that denitrification capacity will likely be maximized in riverine floodplains that are restored as wetlands with diverse native vegetation and enhanced hydrologic connectivity. Such restoration activities promote higher denitrification rates via elevated moisture, fine sediment deposition, and soil organic matter.

Tippecanoe River, Wabash River↗
USGS2024

A seismic nodal deployment to understand magmatic structure in the vicinity of the Pahala earthquake swarm

In summer-fall 2022, 80 three-component SmartSolo IGU-BD3C-5 nodal seismometers were deployed surrounding the Pāhala seismic swarm on the Island of Hawaiʻi, with the goal of improving seismicity catalogs, and seismic velocity images of the crust and upper mantle in this region. The Pāhala swarm, located south of Mauna Loa and Kīlauea, has been the site of a multiyear sustained swarm of seismicity at depths of ~ 25 – 40 km, with order of magnitude increases in rate in 2015, and then again in 2019. This seismicity is possibly related to the input of magma from the mantle plume below, which may then be subsequently transported to volcanic edifices. However, these processes remain enigmatic, in part due to a lack of precise earthquake locations and seismic velocity models in this region. Here we provide an overview of the deployment, an assessment of the quality of the collected data, and discuss the viability of the dataset for local earthquake relocation, tomography, and teleseismic receiver functions. Through comparisons with proximal permanent broadband and short period instruments, we find that the nodes produce high quality data, particularly at periods shorter than 5 s, although we find, document, and correct discrepancies with the gain and polarities of the instruments. We successfully record signals from teleseismic earthquakes, even at periods longer than 5 s (the corner of the flat response of the nodes). We also record local earthquakes, including details related to source characteristics. This indicates that the data is likely to prove useful for investigations using both local and teleseismic earthquake signals to better understand the connections between the deep and shallow magmatic systems of Hawaiʻi. While this deployment provides a snapshot in time, its success may provide a useful benchmark for future studies as the volcanic systems of Hawaiʻi continue to evolve in the future.

Hawaii↗
USGS2024

Debris-flow entrainment modelling under climate change: Considering antecedent moisture conditions along the flow path

Debris-flow volumes can increase along their flow path by entraining sediment stored in the channel bed and banks, thus also increasing hazard potential. Theoretical considerations, laboratory experiments and field investigations all indicate that the saturation conditions of the sediment along the flow path can greatly influence the amount of sediment entrained. However, this process is usually not considered for practical applications. This study aims to close this gap by combining runout and hydrological models into a predictive framework that is calibrated and tested using unique observations of sediment erosion and debris-flow properties available at a Swiss debris-flow observation station (Illgraben). To this end, hourly water input to the erodible channel is predicted using a simple, process-based hydrological model, and the resulting water saturation level in the upper sediment layer of the channel is modelled based on a Hortonian infiltration concept. Debris-flow entrainment is then predicted using the RAMMS debris-flow runout model. We find a strong correlation between the modelled saturation level of the sediment on the flow path and the channel-bed erodibility for single-surge debris-flow events with distinct fronts, indicating that the modelled water content is a good predictor for erosion simulated in RAMMS. Debris-flow properties with more complex flow behaviour (e.g., multiple surges or roll waves) are not as well predicted using this procedure, indicating that more physically complete models are necessary. Finally, we demonstrate how this modelling framework can be used for climate change impact assessment and show that earlier snowmelt may shift the peak of the debris-flow season to earlier in the year. Our novel modelling framework provides a plausible approach to reproduce saturation-dependent entrainment and thus better constrain event volumes for current and future hazard assessment.

Illgraben basin↗
USGS2024

Extending the Boore and Abrahamson (2023) modified square-root-impedance method for the development of site amplifications consistent with the full-resonance approach to a range of VS30 values

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions

Bulletin of the Seismological Society of America↗
USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗
USGS2026

Effects of human activity on Khumbu Glacier: Towards a sustainable Everest Base Camp

The Everest Base Camp (EBC), a critical staging area for mountaineers attempting to summit Mount Everest, has seen a significant increase in human activity over the past few decades. The increasing number of climbers, trekkers, and support staff, particularly at the campsite area, has intensified environmental impacts on the Khumbu Glacier. This study investigates the effects of human activity at EBC, focusing on three contributors to glacier melt: global climate change, local fossil fuel consumption, and human urinary discharge. Through field data collected during the spring 2023 climbing season, this study examines the energy released through the use of liquefied petroleum gas (LPG), kerosene, and petrol for cooking, heating, and electricity, as well as the heat generated by human urinary discharge. The total heat energy released from these activities at EBC in the spring of 2023 was 849,174 ± 179,774 MJ, which would be sufficient to melt 2492 ± 528 tons of glacier ice and snow. Trends in land surface temperature at EBC and adjacent environments during the 1991–2023 period were calculated from Landsat 5, 7, and 9 satellite data. The 32-year Landsat record reveals that EBC surface temperature increased by 0.28 °C year −1 , which was roughly twice the rate of warming at the surface of Khumbu Glacier adjacent to EBC and about 15% higher than the temperature increase of the debris-covered terrain immediately to the north of EBC. The findings suggest that addressing these anthropogenic influences could help to preserve the Khumbu Glacier and support the sustainability of mountaineering in the region. This paper also takes a transformative approach and explores stable and safer locations in case there is a necessity to relocate the current EBC. Two sites have been identified southwest of the existing EBC, which, unlike the current supraglacial site, are situated on stable ground.

Everest Base Camp, Khumbu Glacier↗
USGS2026

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

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
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↗