USGS Science⌕ Search

SEARCH · USGS Science

Results for “River Research and Applications”

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

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

1,666 records · Page 56Linked to original sources
USGS2025

Comparative properties of saponitic fault gouge and serpentinite muds cored from mud volcanoes of the Mariana subduction zone

We obtained 12 core samples for physical and chemical characterization from three serpentinite mud volcanoes (Yinazao, Asùt Tesoru, and Fantangisña) located on the forearc of the Mariana subduction system, that were drilled during International Ocean Discovery Program Expedition 366. Two samples from the Fantangisña mud volcano are interpreted to be clay-rich fault gouges derived from the subduction channel. Their bulk compositions are intermediate between the serpentinites and oceanic basalts. The oceanic crustal materials in the gouges have been thoroughly metasomatized and the serpentinites extensively altered to the trioctahedral, Mg-rich smectite clays saponite and corrensite. The only relict phases in clasts of crustal rock are accessory Ti- and P-bearing minerals. The two fault gouge samples have lower frictional strengths ( μ < 0.2) than the serpentinites ( μ = 0.2–0.4), and their measured permeabilities are also somewhat lower. Their physical and compositional properties correspond to saponitic gouges from other faults that juxtapose serpentinite against crustal rocks, in particular gouges from the two creeping traces of the San Andreas Fault recovered in the core from the San Andreas Fault Observatory at Depth. The décollement beneath Fantangisña mud volcano is thus expected to be very weak and likely characterized by stable slip.

Geophysics, Geochemistry, Geosystems↗
USGS2025

Earthquake probabilities and hazards in the U.S. Pacific Northwest

Earthquakes and their cascading consequences pose a significant threat to the people, environment, infrastructure, and economy of the U.S. Pacific Northwest. The Pacific Northwest is susceptible to three types of earthquakes: deep (intraslab) earthquakes, subduction zone (megathrust) earthquakes, and shallow crustal earthquakes. For each of these earthquake types, earth scientists can use a variety of methods to estimate the probability of occurrence for future events, which constrains seismic hazard and informs building codes. The timing of past earthquakes indicates that there is an 85-percent chance of a magnitude 6.5 or greater deep earthquake in the Puget Sound region; a 10-15-percent chance of an approximately magnitude 9 earthquake on the Cascadia Subduction Zone; and a 17-percent chance of a magnitude 6.5 or greater crustal fault earthquake in the Puget Sound region in the next 50 years. Individuals and communities can take simple steps to prepare for and reduce the impact of future earthquakes.

California, Oregon, Washington↗
USGS2026

Tectonic implications of transitional melting regimes from petrological, geochronological, and compositional characterization of the ophiolitic Seventymile terrane, Alaska, USA

New geochemical, U-Pb geochronology, and Sr-Nd-Hf isotope data provide evidence for the tectonic evolution of the Seventymile terrane in interior Alaska, USA. Ultramafic and mafic rocks of the Seventymile terrane are thought to represent components of a dismembered ophiolite and provide unique constraints on regional terrane evolution and accretion. The Seventymile ophiolite represents fragments of the Devonian to Permian Slide Mountain Ocean (SMO) that separated allochthonous and parautochthonous continental fragments of western North America. It now occurs as multiple thrust sheets containing Permian mafic and ultramafic rocks overlying and/or possibly imbricated with amphibolite-facies supracrustal rocks of allochthonous Yukon-Tanana terrane and parautochthonous North America. Seventymile klippen contain variably serpentinized peridotite, primarily harzburgite, low-grade meta-mafic rocks, and minor oceanic sedimentary rocks (argillite, chert, limestone, and metasandstone). Mafic rocks include gabbro to diabase, typically as dikes, veinlets, or rare massive stocks intruding peridotite. Mafic rocks also include greenstones of the Seventymile assemblage in klippen structurally underlying, and in shear zone contact with, ultramafic klippen. New trace element and radiogenic isotope data from mafic magmatic rocks associated with the Seventymile ultramafic bodies show evidence for a weakly subduction-modified mantle source, like the mantle source of normal mid-ocean-ridge basalt (N-MORB) or back-arc basin basalt (BABB). Seventymile assemblage greenstones are more heterogeneous. They range from N-MORB to enriched mid-ocean-ridge basalt (E-MORB) and ocean-island basalt (OIB), with a subset of samples indicative of continental arc affinity. These geochemistry results indicate that distinct tectonic environments are represented by at least two, and possibly three, lithological and structural units comprising the Seventymile terrane. Hf-Nd isotope systematics are consistent with a depleted MORB mantle (DMM)−like component that overlaps with Pacific MORB. Primary zircon is rare, but new in situ U-Pb data for gabbro and greenstone indicate ca. 274−272 Ma peak zircon and titanite crystallization. Scattered younger zircons define a ca. 255 Ma zircon peak and correspond to secondary crystallization associated with baddeleyite reaction of high-Si fluids during low-grade metamorphism. If Seventymile suites are contemporaneous, obduction associated with the closure of the SMO resulted in the stacking of ophiolitic packages representing distinct tectonomagmatic settings across the transition from pericontinental, to epicontinental, to distal ocean back-arc. Intrusions hosted in klippe of ultramafic rocks, plus the least subduction-modified greenstones underlying them, geologically and compositionally resemble Slide Mountain rocks of the Campbell Range formation in eastern Yukon and may provide a new piercing point across the Tintina fault.

Alaska, British Columbia↗
USGS2024

Geospatial PDF map of the compilation of GIS data for the mineral industries and related infrastructure of Africa

Introduction In 2021, the U.S. Geological Survey's (USGS) National Minerals Information Center (NMIC) completed the project titled "Compilation of geospatial data for the mineral industries and related infrastructure of Africa." This project aimed to leverage the expertise and capabilities of the NMIC to collect, synthesize, and interpret geospatial data to inform on the extractive resources of the African region and expand the NMIC's understanding on the impact of mineral industry of African nations in the global economy. The African region, which comprises the independent nations that make up the African continent and its associated islands and dependencies, consists of a total of 58 mineral producing countries. The primary objective of this effort was to create a fully attributed Geographic Information System (GIS) portraying existing mining infrastructure, resources, and development capacity across Africa along with the related infrastructure capable of supporting current (for the reference year 2018) and future extractive industry operations in the region. The compiled GIS geodatabase with supporting documentation including comprehensive metadata was published as a USGS data release titled "Compilation of Geospatial Data (GIS) for the Mineral Industries and Related Infrastructure of Africa." This georeferenced portable document format (GeoPDF) map sheet presents a new geographic information product containing a partial representation of the GIS data. This GeoPDF map provides a visual comparison of the distribution of mineral industry and related infrastructure GIS data, which contributes to a deeper understanding of the intersections and complexities of the extractive industries within Africa.

Open-File Report↗
USGS2008

Late Jurassic igneous rocks in south-central Arizona and north-central Sonora: Magmatic accompaniment of crustal extension

In Middle Jurassic time, the region of south-central Arizona and north-central Sonora was part of a continental magmatic arc. In early Late Jurassic time, this magmatic arc gradually gave way to a extensional regime. Late Jurassic rifting was accompanied by magmatism, typically small in volume relative to the preceding arc phase, and with mildly alkaline compositions. These igneous rocks, the 158- to 146-Ma Ko Vaya Suite, make up much or most of the Comobabi, Artesa, Quijotoa, and Brownell Mountains and Sierra del Cobre. The bimodal Ko Vaya Suite constitutes volcanic to shallow plutonic complexes comprising chiefly trachyandesite and volcanic wacke; monzodiorite; rhyolite porphyry; distinctive, compositionally and texturally heterogeneous, quartzpoor granite (the Ko Vaya Granite); A-type perthite granite and quartz syenite; aphyric intrusive rhyolite; and locally abundant hematite veins. Alteration, particularly potassic alteration, is widespread. Minimally altered samples of granite and quartz syenite have the composition of trachyte or trachydacite and alkali rhyolite. The Ko Vaya Suite is nonconformably overlain by the latest Jurassic Sand Wells Formation, approximately correlative with the Glance Conglomerate, the basal unit of the Bisbee Group. The Ko Vaya Suite is in some ways analogous to the coeval Independence dike swarm of eastern California, but less diffuse. We infer that the Ko Vaya Suite represents especially intense or focused Late Jurassic extensional magmatism.

Arizona, Sonora↗
USGS2026

Assessment of groundwater quantity and quality contributions to Lake Huron

Lake Huron, one of the five Great Lakes, borders the United States and Canada, with Michigan as the only U.S. State on its shoreline. Like other freshwater lakes, it faces water-quality challenges from nutrients and chemicals applied across its drainage basin. Although past studies focused on surface-water sources, groundwater contributions remain less understood. To address this gap, the U.S. Geological Survey, as part of the Cooperative Science and Monitoring Initiative, classified drainage basins to Lake Huron into eight hydrogeologic zones based on bedrock rock type and glacial sediment transmissivity. Utilizing existing data and empirical field data, we quantified groundwater discharge and identified areas of concern for loading of chloride and nitrate to Lake Huron. Groundwater contributions, including indirect and shoreline discharge, ranged from 5.8 to 11.5 inches annually, totaling 1.9 cubic miles and 0.09 cubic mile, respectively. Hydrogeologic zones with higher glacial sediment transmissivity yielded greater indirect groundwater discharge. Chloride levels above the U.S. Environmental Protection Agency’s 250-mg/L recommendation were mainly in the Saginaw lowlands, whereas nitrate above the 10-mg/L standard was rare—found in only 11 wells. Together, the analysis of where groundwater discharge is occurring in the Lake Huron Basin and the identification of areas with potential groundwater-quality concerns can help prioritize areas that are critical to protecting the long-term health of Lake Huron.

Michigan↗
USGS2024

Taking heat (downstream): Simulating groundwater and thermal equilibrium controls on annual paired air–water temperature signal transport in headwater streams

Headwater stream temperature often exhibits spatial variation at the kilometer-scale, but the relative importance of the underlying hydrogeological processes and riverine perturbations remains poorly understood. In this study, we investigated the relative importance of groundwater (GW) and other processes on downstream annual stream temperature signal characteristics using deterministic heat budget model (HFLUX) scenarios within an idealized stream reach representative of mountainous forested conditions. We summarized annual stream thermal regimes from the relationship of paired sinusoidal air and water temperature signals (amplitude ratio, phase lag, and mean ratio). Results showed that downstream changes in annual temperature depended on the thermal gradient between water and the hypothetical equilibrium temperature (where all heat fluxes sum to zero). GW inflow, riparian shading, and the boundary input signal were the most significant factors affecting downstream annual water temperature signals, while flow volume and channel dimensions impacted how quickly annual temperature signals changed. Effects of GW were dominated by advective rather than conductive heat exchange processes, but conduction played a larger role when GW input was more spatially diffuse. Our results indicated several mechanisms by which local processes may affect stream thermal resilience to disturbances and can help guide management of wildfire and climate change.

Journal of Hydrology↗
USGS2026

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗
USGS2026

Evaluating groundwater quality influences from oil field operations and other anthropogenic activities in an urban setting, Santa Fe Springs, California

Groundwater quality is often affected by anthropogenic activities in urban settings. This study examines groundwater quality in and around the Santa Fe Springs Oil Field in Los Angeles County, California, where oil and gas production commonly intersects with high density industrial, commercial and residential land uses. Utilizing a combination of new and historical data, we evaluated potential pathways that would allow for oil field formation fluids to migrate into groundwater and whether mixing may have occurred based on the distribution of groundwater and oil field formation fluid tracers in samples. Samples were analyzed for a wide array of constituents including volatile organic compounds, light hydrocarbons, major ions, and various isotopic compositions. Despite evidence of oil field infrastructure providing potential pathways of migration via uncemented annular spaces, casing breaches and historical disposal of oil field formation water in surface ponds, the distribution and occurrence of stable isotopes of water, chloride, boron, and total dissolved solids do not indicate mixing of oil field formation water and groundwater. However, methane isotopic signatures and the presence of heavier alkanes suggest gas from oil-bearing formations have migrated from depth via oil field well infrastructure. Volatile organic compound detections were mainly from manufactured compounds unrelated to oil and gas production, with a relatively limited number of petroleum hydrocarbons also detected. Volatile organic compounds were generally found in wells tapping shallow, modern aged groundwater, indicating anthropogenic activities occurring at or near land surface as the source. Study results suggest that while oil field infrastructure provides migration pathways for oil field formation fluids to be introduced into groundwater, urban land uses not related to oil and gas production are the primary drivers of groundwater quality degradation.

California↗
USGS2026

End-user needs for remote sensing wetlands of the Prairie Pothole Region of North America

The Prairie Pothole Region (PPR) of North America comprises globally important grassland and wetland ecosystems critical for numerous populations of migratory birds. Due to the importance of this region for migratory birds, and particularly waterfowl, and the threats of habitat loss due to intensifying agriculture, there is a mature and diverse system of conservation organizations, agencies, and partnerships that spends hundreds of millions of dollars annually on habitat conservation to support migratory bird populations. Remote sensing can be a powerful tool for observing and evaluating global change at large scales as well as expanding inferences from field studies to the broader landscape with statistical models. However, development and utilization of these tools has lagged behind their demand for several reasons, including concerns over spatial and temporal resolution and accuracy of products; perception of a misalignment with decision-maker needs; technological barriers such as skill sets of conservation professionals, computing resources, data access, and usability. In this report, we summarize the needs of conservation professionals and scientists who use or want to use remote sensing data products to inform science about wetland change and conservation of wetlands in the PPR. We assembled this information through several methods leading up to, during, and following a January 2026 PPR Wetland Remote Sensing Workshop. The workshop included United States and Canadian scientists, conservation professionals, and policy experts. Our goal was to bring together end-users and remote sensing product developers jointly to explore reducing the lag between product development and utilization of products to inform science and conservation. Specifically, we aimed to identify gaps in wetland remote sensing that limit effective monitoring, management, and conservation in the PPR, and to develop a framework that outlines pathways to address these gaps by fostering collaboration, improving communication networks, encouraging discussion, and building on existing and ongoing efforts. This report summarizes our participants’ descriptions of end-user needs and the outcomes of the workshop.

Prairie Pothole region↗
USGS2025

Estimated ultimate recovery (EUR) Prediction for Eagle Ford Shale using integrated datasets and artificial neural networks

The estimated ultimate recovery (EUR) is an important parameter for forecasting oil and gas production and informing decisions regarding field development strategies. In this study, we combined site-specific geologic, completion, and operational parameters with the predictive capabilities of machine learning (ML) models to predict EURs of the wells for the Eagle Ford Marl Continuous Oil Assessment Unit. We developed an extensive dataset of wells that have produced from the lower and upper Eagle Ford Shale intervals and reduced the model complexity using principal component analysis. We tested the ML models and estimated the sensitivities of ML-predicted EURs to changes in the values of different input variables. The results of applying the optimized ML model to the Eagle Ford suggest that the approach developed in this study could be promising. The ML estimates of the EURs fit the DCA-based values with an R 2 ~ 0.9 and a mean absolute error of ~36 × 10 3 bbl. In the lower Eagle Ford Shale, the EUR estimates were found to be most sensitive to changes in porosity, net thickness of the interval, clay volume, and the API gravity of the oil; and that in the upper Eagle Ford Shale they were most sensitive to changes in the total organic carbon and water saturation, which suggests that it could be important to consider these parameters in assessing these intervals or close analogs.

Louisiana, Mississippi, Texas↗
USGS2024

Cross section N–N' through the Valley and Ridge province of the southern Appalachian basin, from Greene County, west-central Alabama, to Bibb County, central Alabama

Introduction Geologic cross section N–N′ is the sixth in a series of geologic cross sections constructed by the U.S. Geological Survey to document and improve understanding of the geologic framework and petroleum systems of the Appalachian basin. Cross section N–N′ provides a regional view of the structural and stratigraphic framework of the Appalachian basin in the Valley and Ridge province in western and central Alabama; it spans approximately 69 miles (mi) (111 kilometers [km]). This geologic cross section is a companion to geologic cross sections E–E′ , D–D′ , C–C′ , I–I′ , and A–A′ that are located approximately 350 to 550 mi (563 to 885 km) to the northeast. Cross section N–N' complements earlier geologic cross sections through the Alabama part of the Appalachian basin. Although some of the other cross sections show more structural and stratigraphic detail, they are of more limited extent geographically and stratigraphically. Cross section N–N′ contains information that is useful for evaluating energy resources in the Appalachian basin. Although the Appalachian basin petroleum systems are not shown on the cross section, many of their key elements (such as source rocks, reservoir rocks, seals, and traps) can be inferred from lithologic units, unconformities, and geologic structures shown on the cross section. Other aspects of petroleum systems (such as the timing of petroleum generation and petroleum migration pathways) may be evaluated by burial history, thermal history, and fluid flow models based on what is shown on the cross section. In addition, cross section N–N′ may be used as a reconnaissance tool to identify plausible geologic structures and strata for the subsurface storage of liquid waste or for the sequestration of carbon dioxide.

Alabama↗
USGS2024

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

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗
USGS2023

3D Dynamic rupture modeling of the 6 February 2023, Kahramanmaraş, Turkey Mw 7.8 and 7.7 earthquake doublet using early observations

The 2023 Turkey earthquake sequence involved unexpected ruptures across numerous fault segments. We present 3D dynamic rupture simulations to illuminate the complex dynamics of the earthquake doublet. Our models are constrained by observations available within days of the sequence and deliver timely, mechanically consistent explanations of the unforeseen rupture paths, diverse rupture speeds, multiple slip episodes, heterogeneous fault offsets, locally strong shaking, and fault system interactions. Our simulations link both earthquakes, matching geodetic and seismic observations and reconciling regional seismotectonics, rupture dynamics, and ground motions of a fault system represented by 10 curved dipping segments and embedded in a heterogeneous stress field. The M w 7.8 earthquake features delayed backward branching from a steeply branching splay fault, not requiring supershear speeds. The asymmetrical dynamics of the distinct, bilateral M w 7.7 earthquake are explained by heterogeneous fault strength, prestress orientation, fracture energy, and static stress changes from the previous earthquake. Our models explain the northward deviation of its eastern rupture and the minimal slip observed on the Sürgü fault. 3D dynamic rupture scenarios can elucidate unexpected observations shortly after major earthquakes, providing timely insights for data‐driven analysis and hazard assessment toward a comprehensive, physically consistent understanding of the mechanics of multifault systems.

The Seismic Record↗
USGS2026

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Telecommunications fiber for sensing earthquake aftershocks: Progress and hurdles

Aftershocks offer valuable clues to earthquake behavior. The challenge: quickly deploying sensors to capture the early details of earthquake ruptures within the zone of aftershocks. Telecommunication fibers might be an answer, providing denser networks in otherwise difficult areas, potentially faster than traditional methods.

Conference Paper↗