USGS Science⌕ Search

SEARCH · USGS Science

Results for “Resources, Conservation, and Recycling”

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

Assessing legacy nitrogen in groundwater using numerical models of the Long Island aquifer system, New York

Nitrogen transported along groundwater flow paths in coastal aquifers can contribute substantially to nitrogen loading into surface water receptors, particularly in hydrologic systems dominated by groundwater discharge. Nitrogen entrained in the aquifer is a function of land use and associated nitrogen sources at the time of groundwater recharge, which may differ considerably from present-day sources. Legacy nitrogen can result in substantial discrepancies between observed present-day nitrogen loading to surface water receptors and loading estimated from present-day sources. Additionally, legacy nitrogen can continue to discharge into surface waters after nitrogen mitigation actions have been undertaken. Here, we use a numerical modeling framework to compare three methods of estimating time-varying historical nitrogen loads to four water bodies (receptors) on eastern Long Island, New York. The methods span a range of data requirements and process complexity, from instantaneous receptor loads calculated from steady-state groundwater contributing areas, to transient loads estimated by explicitly simulating legacy groundwater nitrogen transport over a century with large changes in nitrogen sources and hydrologic conditions. The effects of legacy nitrogen on estimated receptor loads varied temporally and spatially within the study area. Depending on antecedent nitrogen inputs and hydrologic conditions, historical annual nitrogen loads estimated from transient simulations accounting for legacy nitrogen can be quite similar (<10% difference) or substantially different (±100%) from those estimated from simpler instantaneous methods. Continued input of present-day nitrogen sources using methods that account for legacy nitrogen results in asymptotic increases in receptor nitrogen loads over time, indicating that simulated present-day receptor nitrogen loads are not in equilibrium with present-day inputs. For these receptors in disequilibrium, models simulating transient groundwater nitrogen transport could be used to account for legacy nitrogen lag times to help resource managers evaluate the potential effectiveness of proposed nitrogen mitigation actions.

EarthArXiv↗
USGS2026

Birds, breakpoints, and baselines: How citizen science data can reveal ecological boundaries in Kenya’s Upper Tana watershed

Tropical watersheds are increasingly threatened by climate change, land-use conversion, and resource extraction, yet conventional biodiversity monitoring in these systems is often spatially and temporally limited. Citizen science offers a complementary approach, enabling biodiversity data collection over large areas that can supplement professional scientific surveys. We analyzed 10 years (2012–2022) of Upper Tana Watershed bird data from the Kenya Bird Map project, covering 114 pentads (9 × 9 km) within a 17,000 km 2 watershed to assess patterns of bird community composition and distribution across this watershed and to also evaluate the effects of environmental variables, seasonality, and sampling effort to help inform improvements in future citizen science projects. Citizen (or community) scientists recorded 575 species (>50% of Kenya’s total avifauna) in 74 families. Asymptotic species accumulation indicates that most probable species present in the watershed were detected. Threshold indicator taxa analysis revealed distinct ecological boundaries along elevation (∼1,500 m), precipitation (∼1,100 mm), and mean temperature (∼19°C) gradients, corresponding to a turnover from xeric savanna to mesic montane forest assemblages. Notably, bird communities showed little seasonal differentiation between wet and dry periods, consistent with dominance by resident year-round species. Data limitations including uneven survey distribution and frequency, absence of abundance metrics, and coarse representation of local environmental conditions that likely reduced our ability to detect fine-scale species–habitat relationships. Addressing these gaps through spatially balanced sampling at greater resolution, greater survey frequency in underrepresented areas, and improved capture of habitat metrics could strengthen the use of citizen-science bird data for watershed bioassessment. Our findings demonstrate that structured citizen-science initiatives can identify ecological boundaries and inform adaptive management of tropical socio-ecological systems under rapid environmental change.

Upper Tana watershed↗
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↗
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↗
USGS2025

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report↗
USGS2025

Characterization of stream water quality and groundwater levels in the Central Pine Barrens region, Suffolk County, New York, 2017–23

The area locally known as the “Central Pine Barrens” region, located in Suffolk County, New York, contains most of Long Island’s preserved and undeveloped land. This region overlays an aquifer system that provides potable groundwater for residents of Suffolk County. Between 2017 and 2023, the U.S. Geological Survey, in cooperation with the Central Pine Barrens Joint Planning & Policy Commission and the Town of Brookhaven, monitored groundwater levels and stream water quality in this region. Groundwater levels were measured monthly at five wells and continuously (15-minute intervals) at a sixth well. Water quality was monitored at five locations in the Carmans River and at two locations in the Peconic River, and samples were analyzed for major ions, trace elements, nutrients, pharmaceuticals, and pesticides. The major ion compositions at the sites were mainly sodium-chloride type waters, and compositions varied the most at the furthest upstream sites in both streams. Concentrations above aquatic-life criteria thresholds also occurred most frequently at the furthest upstream sites. The seasonal patterns of nutrient loads and concentrations varied between the Carmans and Peconic Rivers. Several organic compounds including pharmaceuticals, domestic use products, and pesticides were detected at low concentrations in both streams. Metformin was the most frequently detected pharmaceutical compound, and herbicides were the most frequently detected pesticide class. Water-quality conditions influenced by anthropogenic contributions are a result of current and historical land use, and these contributions include onsite wastewater disposal systems, commercial or domestic fertilizers and pesticides, and urban or industrial contaminants in road runoff. This study characterizes and improves understanding of the current hydrologic conditions in the Central Pine Barrens region and the study findings can help inform the development of plans to manage, protect, and restore water resources.

New York↗
USGS2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

California, Oregon, Washington↗
USGS2025

Implications of physics-based M9 ground motions on liquefaction-induced damage in the Cascadia Subduction Zone: Looking forward and backward

Given the likelihood of future M9 Cascadia Subduction Zone (CSZ) earthquakes, various estimates of the resulting, regional ground motions have been made, including a suite of 30 physics-based simulations that reflect key modeling uncertainties. However, because the last CSZ interface rupture occurred in 1700 CE, the shaking expected in such an event is especially uncertain, as are the impacts to the built and living environments. Like other coseismic impacts, soil liquefaction poses a significant threat and must be considered by any scenario study used to inform planning and response, or to focus mitigation resources. Liquefaction is also notable for its potential to “ground truth” ground-motion estimates, given that its presence or absence in the geologic record can provide constraint on the intensities of shaking in past events. It is thus an important phenomenon looking both forward and backward. Accordingly, using recent physics-based simulations, this study (1) predicts liquefaction in M9 CSZ ruptures at 400 locations in Oregon, Washington, and British Columbia (BC) using an array of cone-penetration-test based models and (2) uses paleoliquefaction evidence at ten sites spanning from Southern Oregon to Vancouver, BC to constrain possible ground-motion intensities experienced in the 1700 CE earthquake. The forward predictions indicate that liquefaction in M9 events could be pervasive in the region and affect numerous population hubs, with the potential for damage across hundreds of square kilometers. The backward analyses suggest that 1700 CE ground-motion intensities may have been less than expected from M9 simulations in some northern portions of the CSZ (e.g. Seattle), given the paucity of 1700 CE liquefaction evidence in these areas. Ultimately, further discovery and analysis of CSZ paleoliquefaction, or lack thereof, will confirm or modify this possibility and the conclusions drawn herein.

Oregon, Washington↗
USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗
USGS2025

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama↗
USGS2024

Estimation of reservoir storage capacity and geomorphic change detection analysis from a multibeam bathymetric survey of Randy Poynter Lake, Rockdale County, Georgia

Rockdale County Department of Water Resources has a directive to update estimates of the reservoir storage capacity of Randy Poynter Lake, located in northern Georgia, and to assess recent sedimentation and associated storage capacity loss. In 2022, the U.S. Geological Survey completed a multibeam bathymetric survey of Randy Poynter Lake to update storage capacity estimates and to quantify storage capacity change since the first multibeam bathymetric survey in 2012 in consideration of estimated errors inherent to bathymetric surveys. Data from the 2022 survey were used to generate contours of the reservoir as well as compute storage capacity at regular increments of water-surface elevation. Storage capacity comparisons between 2012 and 2022 at Randy Poynter Lake show minimal changes that are within the estimated uncertainties, with consistent or slightly increased storage capacities observed at most water-surface elevations and reductions observed at the remaining few elevations. Comparison of the multibeam bathymetric data collected in 2012 with data collected in 2022 further allowed for a formal geomorphic change detection analysis to map, quantify, and infer causation of morphological change over time with respect to a level of detectable change. The volume change in Randy Poynter Lake for the decade between 2012 and 2022 was slightly net-depositional and within the estimated uncertainty. The spatial distribution of sediment deposition was primarily concentrated in the northern portion of the lake, where the principal tributary flows into Randy Poynter Lake. The results of the geomorphic change analysis were used to further understand the future implications to storage capacity change. Despite the challenges of confirming systematic biases because of uncertainties exceeding the observed changes, insights from the study help predict long-term reservoir sediment accumulation, indicating a reservoir half-life extending about 650 years from 2022 on the basis of the current sediment yield estimates.

Georgia↗
USGS2026

Using a temporary emigration model to estimate abundance of stream fishes from hybrid removal surveys with and without block nets

Monitoring programs are often faced with a decision to allocate resources into either robust spatiotemporal coverage to estimate a population index (e.g., not true abundance) or confirming closed sampling conditions (e.g., with block nets) for an unbiased population estimate at the cost of spatiotemporal coverage. However, making accurate and precise abundance estimates at robust spatiotemporal scales is possible when combining open and closed sampling designs with integrated modeling techniques. We used simulations and a case study of backpack electrofishing surveys in the Santa Ana River, California to test the efficacy of an integrated abundance model (temporary emigration model, TE) to estimate abundance of fishes using removal sampling methods with a hybrid sampling design (sampling with and without block nets during removal sampling). We found that the TE model performed well under most modeling scenarios (sample size, amount of closure violation, number of samples collected during closure), although at least a few samples with block nets were necessary for all parameters to be estimable. When applied to fish surveys in the Santa Ana River, we found that catch of the fishes fit to the TE model (Santa Ana Sucker, Arroyo Chub, Channel Catfish, Largemouth Bass, Yellow Bullhead) showed little evidence that the closure assumption was violated when block nets were not used. Additionally, we found that the abundance of non-native fishes negatively affected the abundance of the native Santa Ana Sucker, which was also found to adversely affect the native fish’s access to critical habitat consisting of gravel and cobble substrate. Our results indicate that the TE model presents a viable solution to common sampling problems that impact many monitoring programs, where precise and accurate population estimates can be made at large spatiotemporal scales even when most samples violate the closure assumption.

California↗
USGS2026

Groundwater and surface-water interactions in the Waihe‘e-Kahalu‘u watershed, O‘ahu, Hawai‘i—Analysis of historical data and numerical groundwater-model simulations

The Waihe‘e-Kahaluʻu watershed and surrounding areas lie in a dike-impounded groundwater setting of the Ko‘olau Range, O‘ahu, Hawai‘i, where groundwater withdrawal from wells and tunnels have resulted in reductions in streamflow. Resource managers and surface-water users seek information that can be used to balance the needs of competing uses of groundwater and streamflow in the watershed. This study assesses the effects of four groundwater developments in the Waihe‘e-Kahalu‘u watershed—Kahaluu Tunnel, Waihee Tunnel, Kahaluu Well, and Waihee Incline Wells—on Kahalu‘u and Waihe‘e Streams. Analyses of historical streamflow, withdrawal, and rainfall data indicate that base flow (the component of streamflow that is supplied from groundwater discharge) in Kahalu‘u and Waihe‘e Streams decreased in response to withdrawals from wells and tunnels within and outside the watershed. Groundwater-model simulations indicate that shutdown of selected tunnels and wells in the Waihe‘e-Kahaluʻu watershed will result in increases of groundwater discharge, but not all of the increase will occur in streams within the watershed; increases will also occur in streams outside the basin, the ocean, and tunnels other than those that were shut down. The analysis of historical data and model simulations is consistent in showing that the effects of withdrawal changes on streams are not constrained by the topographic boundaries of watersheds. Stream base flows in the Waihe‘e-Kahalu‘u watershed are subject to a complex interaction of multiple withdrawals from wells and tunnels installed at various times within and outside the watershed. The correspondence between withdrawals and base-flow changes is consistent with the conceptual model and understanding of the connection between groundwater and streamflow in the dike-impounded groundwater setting of the Waihe‘e-Kahalu‘u watershed and the Ko‘olau Range.

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

Aeromagnetic and magnetotelluric imaging of west-central Idaho and the Stibnite-Yellow Pine mining district: A regional to district perspective

Aeromagnetic and magnetotelluric (MT) data are used to better understand the geology and mineral resources near the Stibnite-Yellow Pine mining district in central Idaho. The reduced-to-pole (RTP) transformation of regional-scale aeromagnetic data shows that allochthonous island-arc rocks west of the Salmon River suture are significantly more magnetic than the Laurentian continental rocks east of the suture and that the granitoids of the Idaho batholith have moderate to low magnetization in both early, metaluminous, and late, peraluminous phases. Application of tilt derivative to aeromagnetic data highlights major crustal-scale structures. The 5-km upward continued magnetic data indicate island-arc rocks have deep magnetic sources. The 110-km-long MT profile images resistivity structure to depths around 30 km. At shallow depths, resistivity corresponds to mapped geologic units, with moderate resistivities underlying volcanic and roof-pendant metasedimentary rocks and moderate to high resistivities occurring beneath the Idaho batholith. Crustal-scale moderate resistivities beneath the suture image the results of tectonomagmatic processes that accompanied suturing and translating allochthonous terranes. Low resistivity values beneath and fringing the batholith are derived from metasedimentary rocks that may have served as a melt source and reductant during melt generation and provided metals during later ore formation. In the Stibnite-Yellow Pine mining district, a high-resolution aeromagnetic compilation is shown to correlate with mapped lithologies and mineral deposit-related structures. The RTP transform distinguishes magnetic and nonmagnetic granitoid phases of the Idaho batholith. The tilt derivative highlights metasedimentary rocks, some of which are favorable ore hosts. The Meadow Creek fault hosts the Stibnite and Hangar Flats deposits and is imaged as a magnetic low due to hydrothermal alteration. Reconstructions of magnetic anomaly offsets and orebodies indicate around 3 km of post-95 Ma dextral separation, with some or all of the offset inferred to postdate the main Au mineralization episode (61–66 Ma).

Idaho↗
USGS2026

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period. Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Colorado, Kansas, Oklahoma, Texas↗
USGS2026

Groundwater quality near an oil field in a stream-dominated recharge setting, California, USA

Alluvial valley aquifers are important sources of water supply in many areas but effects of co-located oil and gas development on these resources have not been widely reported, especially in settings where recharge is dominated by stream infiltration. Interpreting the presence of geochemical indicators in the context of hydrology, geology, and other factors provides a more complete understanding of the relations between groundwater and sources of oil-field fluids and aids in identifying risks associated with oil and gas development. Groundwater and Salinas River water samples were collected in an alluvial valley near the San Ardo Oil Field in Monterey County, California and analyzed for a wide range of dissolved chemical, gas, and isotopic constituents to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with fresh groundwater used for supply. Hydraulic gradients, age-dating tracers, and other geochemical indicators show that recharge from the Salinas River has the potential to dilute oil-field fluids that might migrate or seep into the aquifer. Groundwater and Salinas River water collected downgradient of the San Ardo Oil Field showed little or no evidence of mixing with oil-field fluids. Some samples within the oil field contained trace amounts of hydrocarbons or elevated temperatures, indicating that any potential effects from oil-field activities are minor or have been diluted by recharge from the Salinas River. The two samples with the most geochemical evidence of potential mixing with oil-field fluids (SP-18 and GW-17) were collected west of or along the Los Lobos fault, where naturally occurring hydrocarbons are near the land surface. Those samples are also near active or inactive oil-field wells, and so anthropogenic activities and pathways cannot be ruled out as a cause of trace detections of hydrocarbons and elevated temperatures in the aquifer.

California↗