USGS Science⌕ Search

SEARCH · USGS Science

Results for “Critical Reviews in Environmental Science and Technology”

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,679 records · Page 46Linked to original sources
USGS2025

Preserving and increasing water resources—Natural infrastructure in dryland streams in Baja California Sur, Mexico

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land. To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being. Resume La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas. Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Baja California Sur↗
USGS2026

Water use permits as of July 2024 and reported water use near the North Unit of Theodore Roosevelt National Park, North Dakota, 1980–2023

Starting in the early 2000s, increasing oil and gas development in western North Dakota created a need for additional water resources from surface-water and groundwater sources near the North Unit of Theodore Roosevelt National Park. To summarize the use of water in that area, the U.S. Geological Survey, in cooperation with the National Park Service, developed a map of surface-water and groundwater resources, aquifers, and water-use diversions, and plotted water-use trends from 1980 to 2023. Reported water used from permits in the map area has more than doubled since 2020, increasing from about 750 acre-feet in 2020 to about 2,300 acre-feet in 2022 and 2,000 acre-feet in 2023. Surface water provided the primary source of reported water used for the study period with an average of about 410 acre-feet per year from 1980 through 2017 and about 1,330 acre-feet per year from 2018 through 2023. After 2011, groundwater sourced from the Little Missouri River, Tobacco Garden Creek, Fox Hills, Fort Union, and Dakota aquifers became a larger portion of total annual reported water use from permits in the map area. From 1980 through 2015, water use for irrigation averaged 86 percent of the total annual reported surface-water and groundwater use in the map area. Starting in 2011, however, industrial uses became a proportionally larger total use of water, and in 2015, became the highest reported volume of water use in the map area. From 2011 to 2023, industrial use designated for water depots increased from 50 acre-feet to about 1,370 acre-feet, accounting for about 70 percent of total reported water use in the map area in 2023.

North Dakota↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
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↗
USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗
USGS1996

Magnitude and frequency of floods in the Suwannee River Water Management District, Florida

Flood-frequency statistics for 2-, 5-, 10-, 25-, 50-. 100-, 200-, and 500-year recurrence intervals, based on three methods of analysis, are presented for 25 continuous-record and seven peak flow partial-record gaging stations in the Suwannee River Water Management District. The first method, for gaged stations, utilizes station records; the second method, for ungaged sites, utilizes regional regression analysis; and the third method uses a weighted combination of the station and regional values. Because the weighted values utilize two more or less independent estimates of the peak flow statistic, they are considered more accurate than the station estimates or the regression estimates alone. Also, the use of another weighting scheme to improve estimates of flood frequency statistics at ungaged sites is demonstrated. The karstic nature of much of the Suwannee River Water Management District significantly attenuates flood peaks in some streams by providing substantial subsurface storage when river stages are high. At such times, springs discharging into rivers may reverse flow temporarily and become sinks.

Florida↗
USGS2025

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

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

GCAGS Journal↗
USGS2025

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California↗
USGS2025

Bathymetric contour maps, surface area and capacity tables, and bathymetric change maps for selected water-supply lakes in Missouri, 2022–23

Bathymetric data were collected at 13 water-supply lakes around the periphery of Missouri by the U.S. Geological Survey in cooperation with the Missouri Department of Natural Resources and various local agencies, as part of a multiyear effort to establish or update the surface area and capacity tables for the surveyed lakes. Surveys were carried out during the months of April and May in 2022 and 2023. All but two of the lakes had been surveyed previously by the U.S. Geological Survey, and the recent surveys were compared to the earlier surveys to document the changes in the bathymetric surface and capacity of the lake. Bathymetric data were collected using a high-resolution multibeam mapping system mounted on a boat. Supplemental depth data at three of the lakes were collected in shallow areas with an acoustic Doppler current profiler on a remote-controlled boat. Data points from the various sources were exported at a gridded data resolution appropriate to each lake, either 0.82 foot, 1.64 feet, or 3.28 feet. Data outside the multibeam survey extent and greater than the surveyed water-surface elevation were obtained from data collected using aerial light detection and ranging (lidar) point cloud data. A linear enforcement technique was used to add points to the dataset in areas of sparse data (the upper ends of coves where the water was shallow or aquatic vegetation precluded data acquisition) based on surrounding multibeam and upland data values. The various point datasets were used to produce a three-dimensional triangulated irregular network surface of the lake-bottom elevations for each lake. A surface area and capacity table was produced from the three-dimensional surface for each lake showing surface area and capacity at specified lake water-surface elevations. Various quality-assurance tests were conducted to ensure quality data were collected with the multibeam, including beam angle checks and patch tests. Additional quality-assurance tests were conducted on the gridded bathymetric data from the survey, the bathymetric surface created from the gridded data, and the contours created from the bathymetric survey.

Missouri↗
USGS2026

Tectonic controls on volcanism and associated hydrothermal activity in a sediment-dominated mid-ocean ridge; Escanaba Trough

Mid-ocean ridges, the Earth's most extensive volcanic system, exhibit unique characteristics in sediment-dominant environments. Thick sediment cover insulates the crust and channels fluid along pathways that can lead to the formation of distinct crustal alteration patterns, exceptionally large mineral deposits, and specialized chemosynthetic ecosystems. This study presents an interdisciplinary investigation into the tectonics of the Escanaba Trough, a heavily sedimented axial valley at the southern Gorda Ridge in the Northeast Pacific Ocean. A primary challenge in such environments is overcoming the masking effect of thick sediments on basement structures that control magmatic and hydrothermal activity. We address this by employing three-dimensional (3D) magnetic modeling of high-resolution near-seafloor magnetic data collected by an autonomous underwater vehicle (AUV). The 2022 surveys with AUV Sentry provided data for 3D magnetic susceptibility models, refining our understanding of the geometry of sub-sediment laccoliths/saucer-shaped sills and hydrothermal alteration. In conjunction with a new 1:100,000 scale lithostratigraphic map, we outline the tectonic controls on the emplacement of Escanaba Trough's three main volcanic centers, characterize the geometry of its spreading segments, and provide volumetric data on the distribution of sub-sediment volcanism in the southern Gorda Ridge.

Escanaba Trough, Pacific Ocean↗
USGS2024

Flood-inundation maps for the Cuyahoga River in and near Independence, Ohio, 2024

Digital flood-inundation maps for a 9.9-mile reach of the Cuyahoga River in and near Independence, Ohio, were created by the U.S. Geological Survey (USGS) in cooperation with the Northeast Ohio Regional Sewer District Board of Trustees. Water-surface profiles were computed for the stream reach by using a one-dimensional steady-state step-backwater model. The model was calibrated to the current (2024) stage-streamflow relation (rating curve 43.0) for the USGS streamgage 04208000, Cuyahoga River at Independence, Ohio. The resulting hydraulic model was then used to compute 13 water-surface profiles for water levels (flood stages) ranging from 14.00 to 26.00 feet. The flood stages range from “action stage” to above “major flood stage” as reported by the National Weather Service. The simulated water-surface profiles were then used in combination with a digital elevation model derived from light detection and ranging data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood mitigation strategies and to define flood hazard areas to protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid postflood recovery efforts.

Ohio↗
USGS2026

Stream sediment sources in Medicine Creek, northern Missouri and southern Iowa

This report presents the results of a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources to quantify sediment transport source contributions in the Medicine Creek drainage basin. Understanding relative source contributions provides valuable information for selecting the conservation practices that may be most effective in reducing sediment and sediment-associated nutrient transport in the Medicine Creek drainage basin and similar areas of the Lower Grand River drainage basin. Sediment samples were collected from potential contributing areas (source samples) and from fluvial-transported samples (target samples). Source sample types included streambanks, row crop fields, and a combined pastures and forests category. Samples were analyzed for particle size and quantity of carbon, nitrogen, stable isotopes of carbon and nitrogen, and 49 mineral elements as potential tracers. Results for the carbon stable isotope ratio of carbon-13/carbon-12 (δ 13 C) and concentrations of total carbon, total nitrogen, calcium, potassium, and copper were selected by discriminant function analysis as the best combination of multiple tracers to differentiate each source type. The discriminant function analysis poorly differentiated pastures and forests, so these source types were combined. The sources defined by the discriminant function analysis were then used in an unmixing model to apportion sources for each target sample. In the study area, transported sediment was predominantly bank sediment, with an overall average of 86.9 percent of suspended-sediment samples and depositional streambed samples attributed to bank material. Suspended-sediment samples from the mainstem of Medicine Creek were dominated by bank sediments (average of 95.8 percent), and depositional streambed samples from throughout the drainage basin had more variable source contributions with an average of 71.1 percent attributed to bank material. The relative importance of upland sources (row crop fields and the combined pastures and forests category) varied seasonally and with streamflow but was not related to land use or drainage basin size. Relative contributions from upland sources were greater in the summer through winter rather than spring and during lower streamflow, though this may be driven by the seasonality of streamflow. These results indicate management practices that reduce bank erosion could be effective strategies for managing the dominant source of sediment and sediment-associated phosphorus.

Iowa, Missouri↗
USGS1996

Hydrogeologic investigation and simulation of ground-water flow in the Upper Floridan Aquifer of north-central Florida and southwestern Georgia and delineation of contributing areas for selected city of Tallahassee, Florida, water-supply wells

A 4-year investigation of the Upper Floridan aquifer and ground-water flow system in Leon County, Florida, and surrounding counties of north-central Florida and southwestern Georgia began in 1990. The purpose of the investigation was to describe the ground-water flow system and to delineate the contributing areas to selected City of Tallahassee, Florida, water-supply wells. The investigation was prompted by the detection of low levels of tetrachloroethylene in ground-water samples collected from several of the city's water-supply wells. Hydrologic data and previous studies indicate that; ground-water flow within the Upper Floridan aquifer can be considered steady-state; the Upper Floridan aquifer is a single water-bearing unit; recharge is from precipitation; and that discharge occurs as spring flow, leakage to rivers, leakage to the Gulf of Mexico, and pumpage. Measured transmissivities of the aquifer ranged from 1,300 ft 2 /d (feet squared per day) to 1,300,000 ft2/d. Steady-state ground-water flow in the Upper Floridan aquifer was simulated using a three-dimensional ground- water flow model. Transmissivities ranging from less than 5,000 ft2/d to greater than 11,000,000 ft 2 /d were required to calibrate to observed conditions. Recharge rates used in the model ranged from 18.0 inches per year in areas where the aquifer was unconfined to less than 2 inches per year in broad areas where the aquifer was confined. Contributing areas to five Tallahassee water-supply wells were simulated by particle- tracking techniques. Particles were seeded in model cells containing pumping wells then tracked backwards in time toward recharge areas. The contributing area for each well was simulated twice, once assuming a porosity of 25 percent and once assuming a porosity of 5 percent. A porosity of 25 percent is considered a reasonable average value for the Upper Floridan aquifer; the 5 percent porosity simulated the movement of ground-water through only solution-enhanced bedding plains and fractures. The contributing areas were generally elliptical in shape, reflecting the influence of the sloping potentiometric surface. The contributing areas delineated for a 5 percent porosity were always much larger than those determined using a 25 percent porosity. The lowest average ground-water velocity computed within a contributing area, using a 25 percent porosity, was 1.0 ft/d (foot per day) and the highest velocity was 1.6 ft/d. The lowest average ground-water velocity, determined using a 5 percent porosity, was 2.4 ft/d and the highest was 7.4 ft/d. The contributing areas for each of the five wells was also determined analytically and compared to the model-derived areas. The upgradient width of the simulated contributing areas were larger than the upgradient width of the analytically determined contributing areas for four of the five wells. The model could more accurately delineate contributing areas because of the ability to simulate wells as partially penetrating and by incorporating complex, three-dimensional aquifer characteristics, which the analytical method could not.

Florida, Georgia↗
USGS2026

Teach me how to pycap: A high-capacity well decision support tool using analytical solutions in Python

Regulatory agencies in humid temperate environments rely on timely evaluations of streamflow depletion and drawdown to protect aquatic ecosystems and existing water users. Numerical models offer detailed insights, but their complexity and time demands often preclude their practical use in rapid decision-making. We present pycap-dss, an open-source Python package that implements a suite of analytical solutions for estimating streamflow depletion and drawdown. The tool supports superposition of multiple wells and time-varying pumping, enabling cumulative impact assessments in situations with multiple wells and streams. The software is modular and extensible, allowing users to interchange solutions or add new analytical methods. A YAML-based configuration supports batch processing of multiple wells, and an optional AnalysisProject class facilitates integration with regulatory workflows. Rigorous unit and regression testing ensures computational reliability, and continuous integration supports ongoing development. We demonstrate deterministic examples of drawdown where multiple solutions are readily compared and streamflow depletion with multiple wells in the Central Sands region of Wisconsin. We also show the value of Monte Carlo analyses of streamflow depletion in the same Central Sands example, leveraging computational efficiency to evaluate the uncertainty of individual and cumulative streamflow depletion calculations from over 200 high-capacity wells.

Wisconsin↗
USGS1984

Techniques for estimating magnitude and frequency of floods on streams in Indiana

Equations are presented for estimating the magnitude and frequency of floods at ungaged sites on unregulated and nonurban streams in Indiana. The equations were developed by multiple-regression, analysis of basin characteristics and peak-flow statistical data from 242 gaged locations in Indiana, Ohio, and Illinois. The State of Indiana was divided into seven areas on the basis of the regression analysis. A set of equations for estimating peak discharges with recurrence intervals of 2, IO, 25, 50, and 100 years was developed for each area. Significant basin characteristics in the equations are drainage area, channel length, channel slope, mean annual precipitation, storage, precipitation intensity, and a runoff coefficient. Standard errors of estimate for the equations range from 24 to 45 percent. Methods are also presented for estimating flood magnitude and frequency at sites on gaged streams. Flood-frequency data based on observed peaks are given for 270 gaged locations. Twenty of these are on regulated streams, and six are on urban streams. Basin characteristics are also included car 245 of the gaged locations on unregulated and nonurban streams. No techniques are given for estimating flood magnitude and frequency at ungaged sites on regulated or urban streams. A rainfall-runoff model was tlsed to synthesize long-term peak data at 11 gaged locations on small streams. Flood-frequency curves developed from the long-term synthetic data were combined with curves based on short-term observed data to provide weighted estimates of flood magnitude and frequency at the rainfall-runoff stations.

Indiana↗
USGS2024

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah↗
USGS2025

Geochemical and hydrological investigations of historical data collected at the Lee Acres Landfill and Giant Bloomfield Refinery, New Mexico, 1985–2020

The Lee Acres Landfill and Giant Bloomfield Refinery are adjacent properties near the City of Farmington, New Mexico, each having undergone monitoring and remediation related to historical site activities. At the landfill, site cleanup has included the installation of a capillary barrier over former liquid waste lagoons and periodic monitoring of groundwater elevations and groundwater quality. At the refinery, remediation has focused on several petrochemical and crude oil release areas and included soil excavation, groundwater treatment, and regular monitoring of groundwater elevations and quality. Groundwater at both sites has higher concentrations of volatile organic compounds and trace metals than background aquifer concentrations. In 2022, the U.S. Geological Survey compiled the Lee Acres-Giant Bloomfield Refinery Database (LAGBRD), which contains publicly available groundwater-elevation data and organic and inorganic groundwater-quality data from both sites, spanning from 1985 to 2020. Data from the LAGBRD and precipitation data from other sources were used to better understand the cause of relatively high manganese concentrations observed in some groundwater wells at the site through comparison of groundwater chemistry to chemical end members, interpretation of spatial and temporal patterns in the groundwater chemistry, and interpretation of groundwater flow properties. In this study, elevated chloride concentrations in groundwater downgradient from the landfill have been attributed to landfill leachate based on the temporal and spatial variability of chloride concentrations and chloride-to-bromide ratios. Installation of a capillary barrier and surface-water runoff controls at the landfill in 2005 appears to have altered infiltration patterns at that site, resulting in a decrease in chloride at some wells but an increase in chloride and dissolved manganese at others. The timing and relation among groundwater elevation, chloride concentration, and manganese concentration suggest that leachate stored in the vadose zone provides a continued source of contamination to groundwater.

New Mexico↗
USGS2025

Long-term trends in microseismicity during operational shut-ins at the Coso Geothermal Field, California

Pausing injection and production can lead to induced seismicity in a variety of settings, with some of the largest events occurring during these so-called shut-ins. In geothermal fields, shut-ins are periodically conducted for maintenance on wells and surface infrastructure, thereby offering recurring means of estimating stress changes in the subsurface that lead to increased seismicity rates. Here, we examine daily production and microseismicity data from the Coso Geothermal Field (CGF) in California between 1996 and 2010 to investigate the repetitive trends of operational shut-in microseismicity. Using 24 local seismic stations, we first analyze spatial and temporal trends of over 60,000 earthquakes with magnitudes between -0.4 to 3.8. We find that the northern region exhibits no significant seismicity changes during shut-ins, whereas the rest of the field experiences induced seismicity during almost every shut-in with an increasing intensity towards the southern and eastern portions of the field, highlighting local differences in stress within the CGF. Additionally, we cluster the seismicity using waveform cross-correlation, revealing several earthquake clusters primarily occurring during shut-in periods. These observations suggest that certain fracture and fault sections respond quicker to changes in pore pressure and poroelastic stresses within the geothermal system, possibly highlighting main fluid pathways.

California↗