USGS Science⌕ Search

SEARCH · USGS Science

Results for “Regional Environmental Change”

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.

955 records · Page 51Linked to original sources
USGS2024

Predicting hydrothermal reservoir depth from chemical geothermometers using a three-dimensional temperature model in the Great Basin, USA

Recent work in the Great Basin region of the western United States has made it possible to predict the depth of hydrothermal reservoirs (i.e., the depth at which heat is accumulated prior to ascent via hydrothermal upflow) identified through geochemistry and to contextualize the spatial patterns of these reservoir depths. Chemical geothermometers use the chemical and mineral constituents of hydrothermal fluids to predict the temperature at which fluids equilibrated with the host rocks at depth. Assuming that most of the Great Basin is dominated by conductive conditions until a vertically connected hydrothermal flow path is created (e.g., by faulting), geothermometers reflect the chemical and thermal conditions at the depth interval that the fluid has conductively equilibrated over a long period before a vertical conduit allows convective upflow. By pairing geothermometer temperature estimates with our recent three-dimensional temperature model of conductive heat flow in the Great Basin, we estimate the corresponding reservoir depths and construct a map of circulation depths. The predicted depths from geothermometers have spatial patterns across the Great Basin that relate to patterns seen in other geologic and geophysical data. Deeper springs generally occur disproportionately in areas with higher strain rates and in basins. We posit that current elevated strain rates reflect patterns of historic deformation where ongoing tectonic activity maintains permeable pathways to deeper reservoirs, some of which are estimated to exceed 6 km depth. Basins, as expected, contain a disproportionate number of these deep systems, because the underlying aquifers are closer to the surface in basins, thus requiring less water pressure to reach the surface than in mountain ranges. Most springs estimated to have their source in a deep reservoir occur at places known to host a hydrothermal system; these refined depth estimates of the source reservoir can help to better constrain the source depth for many known hydrothermal systems across the Great Basin.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
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

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗
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↗
USGS1995

Effects of waste-disposal practices on ground-water quality at five poultry (broiler) farms in north-central Florida, 1992-93

Waste-disposal areas such as chicken-house floors, litter stockpiles, fields that receive applications of litter, and dead-chicken pits are potential sources of nitrates and other chemical constituents in downward-percolating recharge water. Broiler- farms in north-central Florida are concentrated in a region where the Upper Floridan aquifer is unconfined and susceptible to contamination. Eighteen monitoring wells installed at five sites were sampled quarterly from March 1992 through January 1993. Increases in median concentrations of constituents relative to an upgradient well were used to determine the source of the nitrate at two sites. At these sites, increases in the median concentrations of nitrate as nitrogen in ground water in the vicinity of waste-disposal areas at these sites were: 5.4 mg/L for one chicken house; 9.0 mg/L for a second chicken house; 2.0 mg/L for a fallow field that received an application of litter; and, 2.0 mg/L for a dead-chicken pit. At the three remaining sites where the direction of local ground-water flow could not be ascertained, the sources of concentrations of nitrate and other constituents could not be determined. However, median nitrate concentrations in the vicinity of waste-disposal areas at these sites were: 45.5 mg/L for a set of two chicken houses; 3.0 mg/L for a stockpile area; and 2.1 mg/L for a hayfield that received an application of litter. The nitrate concentration in ground water in the vicinity of a field that had previously received heavy applications of litter increased from 3.0 mg/L to 105 mg/L approximately 4 months after receiving an application of commercial fertilizer. Increases in concentrations of organic nitrogen in ground water in the vicinity of waste-disposal areas may be related to the decomposition of litter and subsequent movement with downward percolating recharge water.

Florida↗
USGS2026

The potential impact of three-dimensional distributed slip models derived from real-time GNSS data on the performance of the ShakeAlert earthquake early warning system for slab interface earthquakes

The ShakeAlert® earthquake early warning (EEW) system is designed to warn users of imminent strong ground motion with sufficient time to take protective actions. ShakeAlert currently uses three algorithms to characterize the earthquake source. One estimates the location and magnitude using the first few seconds of the P wave and, while fast, tends to underestimate magnitude for M w 7.0+ earthquakes. A second estimates the location, orientation, length, and corresponding magnitude of a line source using observed peak ground acceleration and contributes primarily to M w 5.5+ earthquakes. The third infers earthquake magnitude from peak ground displacement measured using Global Navigation Satellite System (GNSS) data and offers nonsaturating magnitudes for M w 7.0+ earthquakes. Other EEW algorithms exist that infer temporally evolving spatially variable slip on a 3D fault surface from real‐time GNSS data, information that might enable more accurate and timely alerts in the event of large‐magnitude subduction interface earthquakes. Here, we evaluate the potential contribution of one such algorithm, BEFORES ( Minson et al. , 2014 ), to improve ShakeAlert performance through a simulated real‐time implementation of Bayesian evidence‐based fault orientation and real‐time earthquake slip (BEFORES) and other ShakeAlert algorithms using data for eight M w 7.6+ earthquakes. The test results demonstrate that BEFORES can produce well‐constrained and accurate magnitude estimates as soon as or sooner than other EEW algorithms, in turn enabling it to increase the amount of warning time users receive in many cases. However, with a modified Mercalli intensity (MMI) threshold of 3.5, which is commonly used for issuing alerts, BEFORES would tend to alert large geographic regions that did not feel strong shaking (MMI 6+). This effect can be mitigated using a higher alert threshold of MMI 4.5 without negative impact on the amount of warning time obtainable with BEFORES.

Bulletin of the Seismological Society of America↗
USGS2022

Waters divided: A history of alluvial fan research and a view of its future

Flows exiting confined valleys tend to deposit sediment in fan-shaped landforms. Where deposition is wholly or largely by the tractive forces of flowing water, these landforms are called alluvial fans. They are the product of the progressive division of water and sediment downfan, from slopes that may exceed 0.10 to distal slopes that may be below 0.01. Channel depths also tend to decline, from values that approach one to several meters at steep fanheads, to a few decimeters at distal fan margins. The result is a radiating, depositional ramp where confined or unconfined flows transport sediment from source basins to bounding streams, subsiding basins, or stable platforms. Where streams or subsiding basins consume the sediment supply from the source basin, the fan may approach a steady form whose extent and distal slope are set by stream location or subsidence rate. Where boundary conditions do not remove sediment, the fan may prograde out to long distances and low slopes (<0.01). Theoretical and experimental work over the past several decades support the notion that alluvial fan long-profiles become steeper as sediment supply increases or transport capacity decreases, and increasingly concave upward as the rate of bed material deposition decreases downfan. Grainsize distributions of alluvial fans seem to span the range observed in alluvial rivers, with no processes that uniquely identify them, apart from the distributary pattern of deposition. Bed sand cover tends to increase downfan in arid-region fans, with an absence of systematic downfan fining of coarser grain sizes. Surficial mapping and geochronology have demonstrated that fan deposition varies greatly through time, arguably from climate variations that alter hillslope sediment supply. The combination of surficial mapping and hydraulic modeling with high-resolution topography can now produce detailed flood susceptibility maps. The effective use of these maps to protect lives and property, however, depends on answering many of the enduring questions about the mechanics of how water and sediment divide down alluvial fans.

Book chapter↗
USGS2026

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-​central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico↗
USGS2025

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

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

California↗
USGS2025

Sources of water and salts for the Zuni Salt Lake in west-central New Mexico

The Zuni Salt Lake is located in a maar in west-central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. There have been several investigations suggesting different sources for the water and salt to the lake. Springs, seeps, and ephemeral streamflow have all been observed to contribute freshwater to the lake, and brackish to hypersaline seeps have been documented along the banks of the lake. This report summarizes the findings of a study that characterizes the lake’s hydrology, its water and salinity sources, and the hydrogeologic conceptual model. Regional groundwater levels indicate that each of the aquifers in the area have the potential to discharge groundwater to the lake. There is also evidence of vertical groundwater flow pathways at the maar that were likely created by the igneous intrusion that fractured the intersecting aquifers. A detailed water budget was constructed from continuous lake stage, precipitation, and evaporation data to estimate the groundwater inflow to the Zuni Salt Lake. It was determined that groundwater inflow to the lake is 441 ±94 acre-feet per year, which composes as much as 77 percent of the total inflows. The high sodium and chloride concentrations measured in two hypersaline samples collected near the lake indicate that the majority of the dissolved solids entering the lake are from a hypersaline groundwater source. The geochemical and isotopic compositions measured in the lake and surrounding features support the interpretation that hypersaline groundwater is the primary source of salts to the lake, which is likely sourced from the older (and deeper) Permian units. The hypersaline groundwater samples collected during this investigation have a unique aqueous chemistry relative to each of the mapped aquifers, and variability in groundwater compositions is interpreted to result from differences in minerology and residence time.

New Mexico↗
USGS2023

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia↗
USGS2026

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

South Dakota↗
USGS2024

Three-dimensional hydrogeologic framework model of the Russian River watershed, California

The Russian River watershed is in northern Sonoma County and southern Mendocino County, California, in the northern part of the California Coast Ranges. The Russian River serves as a supply for agricultural irrigation and for municipal, domestic, and commercial uses. Through a cooperative agreement with the California State Water Resources Control Board and Sonoma County Water Agency, the U.S. Geological Survey has completed studies to better understand the hydrogeologic system and develop numerical hydrologic modeling tools to evaluate and aid in managing groundwater resources. This report focuses on the development of a digital three-dimensional hydrogeologic framework model of the Russian River watershed for use in groundwater resource assessment and numerical models. The digital three-dimensional hydrogeologic framework model of the Russian River watershed portrays the altitude, thickness, and extent of five hydrogeologic units. These five hydrogeologic units include (1) a basement unit, (2) the Sonoma Volcanics, (3) a consolidated sedimentary rock unit, (4) an unconsolidated sediment unit, and (5) channel alluvium. Model input data were compiled from published geologic maps, interpreted well data, and a model of the top of basement derived from gravity data. These data were used to construct surfaces that represent the upper and lower subsurface boundaries of each hydrogeologic unit. Top surfaces were created for the five hydrogeologic units and then stacked in three dimensions to create a solid-volume digital model. The digital three-dimensional hydrogeologic framework model described in this report and the corresponding data represent the generalized geometry of the subsurface geologic units; the model reproduces the input geologic data with reasonable accuracy and is consistent with previously published subsurface conceptualizations of the region. The model indicates the overall geometry of the basement within the watershed and the spatial extent, altitude, and thickness of the basin-filling units. The hydrogeologic framework model is at a scale and resolution appropriate for use as the foundation for a numerical hydrologic model of the study area.

California↗
USGS2025

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

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

California, Idaho, Nevada, Oregon, Washington↗
USGS2025

Submarine groundwater discharge creates cold‐water refugia that can mitigate exposure of heat stress in nearshore corals

Coral reef mortality around the world is accelerating due to human activities and rising sea temperatures that cause bleaching, which is expected to become more frequent. Our ability to predict which corals will be most resilient, however, remains limited due to insufficient information characterizing nearshore temperature and habitat conditions. In this study, we examine how submarine groundwater discharge (SGD) reduces nearshore water temperatures and exposure of corals to heat stress, complementing the understanding that SGD can adversely affect coral when it contains elevated nutrient concentrations. Data from fixed nearshore sensors and vertical depth profiles along ~100 km of the western shoreline of the Island of Hawai’i from 2003 to 2014 demonstrate that submarine groundwater discharge (SGD) can reduce nearshore water temperatures by 1 °C–5°C and create estuarine-like conditions with salinities as low as 20 PSU, where the prevalent coral species, Pocillopora meandrina, Porites lobata , and Montipora capitata , thrive. Time-series temperature records reveal that exposure to high ambient ocean temperatures, which are known to initiate bleaching events, are reduced up to 5%–46% of the time. Coral health surveys indicated coral bleaching in response to moderately high annual temperatures in 2010 and 2011, with more colonies affected farther from cold, SGD-fed waters. Synthesis of these results, along with coral response data following the more extreme marine heat wave of 2014–2015, demonstrates lower coral loss and greater coral recovery near groundwater seeps, particularly those with higher flux and influence on reducing nearshore water temperatures. Our results demonstrate that SGD may therefore provide a beneficial ecosystem service and enhance coral reef resilience, particularly where human-related nutrient additions to groundwater can be mitigated. The implications of our findings are relevant across tropical coasts where groundwater inputs can be substantial, such as the Caribbean and Indo-Pacific, and contribute to improving our understanding of coral sensitivity to gradients in temperature and nutrient stress. Improved management of groundwater resources could thus be vital to local–regional strategies for mitigating future heat stress.

Hawaii↗
USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗
USGS2025

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2026

Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals

The Mountain Pass carbonatite stock hosts a world-class rare earth element deposit and may be classified as a carbonate-sulfate igneous rock, as it contains on average > 50 volume percent carbonate minerals and 20 to 30 volume percent sulfate minerals. The sulfates range in composition from barite to celestine and locally occur with sparse sulfide minerals. We investigate the origin of sulfur enrichment and the occurrence of sulfur-bearing minerals in the Mountain Pass carbonatite with in-situ sulfur isotope and mineral chemistry. Barite cores with δ 34 S of 1 to 3‰ do not coexist with sulfides, whereas celestine rims with δ 34 S of > 3‰ are associated with sulfides with δ 34 S < -10‰. We propose a model in which sulfur-bearing sediments were subducted during episodes of plate convergence in the Mojave Province that preceded Mountain Pass magmatism. Metasomatism of the overlying mantle by melts derived from the subducted sediments generated an unusually carbon- and sulfur-rich source to yield carbonatite magmas. Sulfur from primary carbonatite magmas and ~ 1 to 7% sulfur from subducted sediment melts yielded a slightly enriched δ 34 S composition (relative to depleted mantle δ 34 S of -1‰) for early crystallizing barite. Celestine rims on magmatic barite cores formed at low, hydrothermal temperatures (< 350 °C) based on S isotope thermometry for equilibrium celestine-galena and celestine-pyrite pairs. The sparse sulfides in the carbonatite stock are not in equilibrium with the primary barite cores and therefore do not permit S isotope thermometry estimates of magmatic temperatures. The S/Se ratios of sulfide minerals (> 3,400) typically exceed primitive mantle values (S/Se of 3,340), also consistent with their derivation from hydrothermal fluids. Trace occurrences of sulfide and sulfate minerals in alkaline silicate stocks related to the carbonatite stock have similar δ 34 S compositions and yield similarly low formation temperatures, suggesting regionally extensive and chemically similar sulfur-bearing hydrothermal fluids that imparted lithologically diverse rocks with a consistent sulfur isotope fingerprint.

California, Nevada↗