USGS Science⌕ Search

SEARCH · USGS Science

Results for “Nature Water”

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.

807 records · Page 44Linked 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↗
USGS2024

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report↗
USGS2024

The advantages of electric-grade heat storage in long closed-loop wells

Though long (e.g., tens of kilometers) closed-loop geothermal wells are under consideration for the extraction of subsurface heat, these wells might also serve as an efficient energy storage mechanism for electricity generation. Using a semi-analytic model solution, the potential for electric-grade heat storage as a function of ambient temperature (e.g., corresponds to depth of loop), well length, well diameter, and flow rate is evaluated. The following simplified cases for comparison of standard closed loop operation with energy storage operation were considered: [1] constant flow rate and constant temperature (90 °C) injection to represent the standard closed-loop base-case; and [2] constant flow rate and annual cycle sinusoidal temperature (90-150 °C) to represent seasonal (summer) charging while solar resources are peak. Initial temperatures for all scenarios considered herein are a uniform 175 °C. Electric-grade heat is assumed to be delivered whenever temperatures at the extraction point exceed 90 °C. For calculation purposes only, if temperature falls below 100 °C, it is assumed that heat delivered is sub-economic, so no electricity would be produced. For all scenarios, temperatures at the extraction well asymptotically approach the flow-weighted average injection temperature, but energy storage scenarios exhibited a damped time-varying signal that diminishes in magnitude with length of the loop. The asymptotic approach depends on initial temperatures in the rock and the heat extraction rate (a function of well diameter and flowrate). This analysis demonstrates that shorter closed loops can produce more electricity over time than longer closed-loops previously proposed for electricity production over typical engineering design lifetimes (e.g., 30 years). Although only a high-temperature scenario is considered herein, rock that is initially below boiling temperature would not host a standard closed-loop resource, but injection of hot water seasonally would asymptotically heat this low-temperature system to temperatures capable of electricity production. In other words, regardless of initial temperatures, closed loops could be used to store electricity with no critical minerals in the geothermal battery.

Conference Paper↗
USGS2024

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Earth and Planetary Science Letters↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗
USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗
USGS2025

The δ13C signature of dissolved organic and inorganic carbon reveals complex carbon transformations within a salt marsh

Coastal wetlands have high rates of atmospheric CO 2 uptake, which is subsequently respired back to the atmosphere, stored as organic matter within flooded, anoxic soils, or exported to the coastal ocean. Transformation of fixed carbon occurs through a variety of subsurface aerobic and anaerobic microbial processes, and results in a large inventory of dissolved carbon. Carbon source and the roles of aerobic respiration, sulfate reduction, and methane cycling were evaluated within salt marsh peat and the underlying sandy subterranean estuary. There is a large increase in dissolved inorganic carbon (DIC, 7,350 ± 3,900 μmol L −1 ), dissolved organic carbon (DOC, 1,040 ± 1,480 μmol L −1 ) and CH 4 (14.5 ± 33.3 μmol L −1 ) within the marsh porewaters compared to creek waters. Alkalinity production (5,730 ± 2,170 μeq L −1 ) and sulfate removal (1,810 ± 1,970 μmol L −1 ) indicate anaerobic respiration, however, relative contributions from the various decomposition pathways cannot be identified due to overlapping geochemical signatures. The δ 13 C of the DOC (−29.0 ± 3.7‰) and DIC (−11.2 ± 1.1‰) produced within the marsh differed from the bulk soil organic matter δ 13 C (−14.5 ± 0.2‰). We explore a variety of mechanisms that could result in co-occurring depleted δ 13 C-DOC and enriched δ 13 C-DIC compared to the bulk soil organic carbon pool and salt marsh vegetation, including selective mineralization, production of δ 13 C-depleted bacterial biomass, and methane-derived DOC. While important questions remain about carbon cycling pathways, we found evidence of a cryptic methane cycle. Alteration of the δ 13 C of carbon species complicates source attribution in solid and dissolved phases and careful consideration should be used when carbon is partitioned between in situ salt marsh production and external marine and terrestrial sources.

Massachusetts↗
USGS2024

On-orbit calibration and performance of the EMIT imaging spectrometer

The Earth surface Mineral dust source InvesTigation (EMIT) is a remote visible to shortwave infrared (VSWIR) imaging spectrometer that has been operating onboard the International Space Station since July 2022. This article describes EMIT's on-orbit spectroradiometric calibration and validation. Accurate spectroscopy is vital to achieve consistent mapping results with orbital imaging spectrometers. EMIT takes a unique approach to this challenge, with just six optical elements, no shutter, and no onboard calibration systems. Its simple design focuses on uniformity and stability to enable vicarious spectroradiometric calibration. Our experiments demonstrate that this approach is successful, approaching the fidelity of manual field spectroscopy in some cases, and enabling new and more accurate products across diverse Earth science disciplines. EMIT achieves several notable firsts for an instrument of its class. It demonstrates successful on-orbit adjustments of Focal Plane Array (FPA) alignment with sub-micron precision. It offers spectral uniformity better than 98%. Optical artifacts in the measurement channels are at least three orders of magnitude below the primary solar-reflected surface signals. Its noise performance enables percent-level discrimination in the depths of mineral absorption features. In these aspects, EMIT satisfies the stringent performance needs for the next generation of VSWIR imaging spectrometers to observe the Earth's ecosystems, geology , and water resources.

Remote Sensing of Environment↗
USGS2025

Evaluating mass flow meter measurements from chambers for greenhouse gas emissions from orphan wells and other point sources

This study evaluates the performance of a rigid gas flux chamber equipped with a mass flow meter (MFM) for measuring gas emissions from leaking orphan wells and similar pressure-driven gas point sources. We conducted a series of laboratory and field experiments to evaluate the sensitivity, stability, and dynamic range of an MFM chamber system and found an optimal method for sealing the chamber to the ground to isolate the emission source. From these results, we estimate the effects of different soil gas permeabilities on measurements and identify the uncertainty of environmental processes that can impact measurements. Simulations of an MFM chamber are compared to those of a dynamic flux chamber to contrast the data derived with both methodologies and illustrate the potential for measuring high variability leaks with the MFM chamber. Using a low flow resistance MFM and a chamber well-sealed to the ground, it is possible to measure leaks down to 1.08 x 10 -3 cubic meters per hour (m 3 h −1 ) (refenced to 25°/1 atm), corresponding to 0.77 grams per hour (g h −1 ) methane or 2.11 g h −1 carbon dioxide, with a mean uncertainty of 0.89 % relative standard deviation. Environmental processes such as heated gas inside the chamber from solar gain, wind blowing across the chamber vent, and changing humidity in the chamber, can cause variation in MFM measurements. Over 11 d of continuous monitoring under varying weather conditions, the standard deviation of the environmentally sourced signals was found to be 7.40 x 10 -3 m 3 h −1 (equivalent to or 5.27 g h −1 methane or 14.45 g h −1 carbon dioxide). Strategies to obtain the highest quality data from MFM chambers include burying the edges of the chamber below the surface sufficiently deep to seal the chamber edges against gas flow and soaking the dirt with water to lower the chances of escaping gases, while monitoring the gas flow and adjusting the chamber seal to achieve a maximum flow rate.

Atmospheric Measurement Techniques↗
USGS2026

Repeated coral bleaching events have eroded reef accretion potential in the Florida Keys

Coral bleaching is now the primary driver of coral-reef degradation globally, threatening the persistence of coral-reef functions and the invaluable ecosystem services they support. The Florida Keys have been impacted by four severe, regional-scale coral-bleaching events, which occurred in 1997–1998, 2005, 2014–2015, and 2023. Our aim in this study was to quantify the effects of those thermal-stress events on the persistence of coral-reef structures. We leveraged long-term coral-reef monitoring data to quantify the impacts of coral bleaching on coral cover, carbonate production, and reef-accretion potential throughout the Florida Keys from 1996–2024. Our study indicates that coral bleaching drove significant declines of reef-building corals, which successively diminished regional-scale reef-accretion potential (from 1.65 ± 0.33 mm y -1 in 1996) and culminated with a regional-scale transition to net reef erosion (-0.12 ± 0.08 mm y -1 ) following the 2023 event. Coral bleaching was the primary driver of declining reef-building capacity, with the cumulative impacts of the four coral-bleaching events accounting for >75% of total declines since 1996. The 1997–1998 event alone drove ~50% of overall declines, largely due to losses of the shallow-water ecosystem engineer, Acropora palmata . More moderate impacts of thermal stress on regional reef-accretion potential in 2005 (6%) and 2014–2015 (9%) could lend some support to the hypothesis that the remaining corals were more thermally tolerant; however, the unprecedented thermal stress event in 2023 caused widespread coral-bleaching and bleaching-related mortality, and caused an additional 12% of the total decline in reef-accretion potential. Stony-coral tissue loss disease accounted for an additional 6% of the decline, with the remaining 19% unexplained. Coral bleaching has already dramatically reshaped the structure and function of coral reefs in the Florida Keys over the last three decades. With the high likelihood that thermal stress will accelerate moving forward, the persistence of Florida’s essential coral-reef functions and services is uncertain.

Florida↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2024

Peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. This work is part of that reexamination. The stationarity assumption is a concern because flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done to assess potential nonstationarity in peak flows in the north-central United States. This report summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Scientific Investigations Report↗
USGS2025

Characterization of the hydrogeologic framework, groundwater-flow system, geochemistry, and aquifer hydraulic properties of the shallow groundwater system in the Wilcox and Lorraine process areas of the Wilcox Oil Company Superfund site near Bristow, Oklahoma, 2022

The Wilcox Oil Company Superfund site (hereinafter referred to as “the site”) was formerly an oil refinery northeast of Bristow in Creek County, Oklahoma. Historical refinery operations contaminated the soil, surface water, streambed sediments, alluvium, and groundwater with refined and stored products at the site. The Wilcox and Lorraine process areas are where the highest concentrations of volatile organic compounds, semivolatile organic compounds, polycyclic aromatic hydrocarbons, and trace elements (including metals) (collectively hereinafter referred to as “contaminants”) were measured in a local shallow perched groundwater system within the alluvium (hereinafter referred to as the “alluvial aquifer”) at the site during previous site assessments. In order to understand the potential migration of contaminants through the soil and groundwater in these areas, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, investigated aquifer characteristics of the alluvial aquifer in the Wilcox and Lorraine process areas of the site to (1) document hydraulic conductivity and other aquifer characteristics of the alluvial aquifer that govern contaminant fate and transport, (2) describe the geospatial extent and concentration of the contaminants in the alluvial aquifer in the Wilcox and Lorraine process areas, and (3) describe the geochemical controls pertaining to oxidation and reduction governing the fate and transport and the degradation potential of contaminants in the groundwater. Various data were compiled and collected to evaluate the aquifer characteristics at the site including the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer. A total of 20 new (2022) groundwater monitoring wells were installed at the site to collect data used to supplement groundwater-level altitude and groundwater-quality data collected from older, existing groundwater monitoring wells and piezometers. Data compiled and collected for the study were used to evaluate the characteristics of the alluvial aquifer at the site. These aquifer characteristics are defined by the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer.

Oklahoma↗
USGS2024

Assessment of the interconnection between Tampa Bay and the Floridan aquifer system: Historical groundwater data compilation and analysis, 1976–2022

The U.S. Geological Survey used existing data collected after the last major navigational channel modification in the mid-1980s to investigate groundwater levels and chloride concentrations in wells in the Floridan aquifer system and other aquifers beneath and near Tampa Bay. Tampa Bay is located on the west-central coast of Florida and provides access for commercial shipping. In 2021, the U.S. Army Corps of Engineers began to investigate alternatives to improve the efficiency of the deep-draft navigation channels within Tampa Bay, reduce costs, improve safety, and analyze the implications of modifying navigational channels. The Floridan aquifer system underlies Tampa Bay and is the primary source of public water supply in the region. Modifications to the channels have the potential to expose the Floridan aquifer system to the saltwater in Tampa Bay, with the potential to increase salinity in wells in the region. Other factors affecting the groundwater levels and the location of the freshwater/saltwater interface were also examined, including changes in sea level, groundwater extraction, and variations in climate. Groundwater levels and well-construction reports were used to identify whether different aquifer units are well-connected. Twelve wells had available data before and after the last major channel modifications, which took place in the 1980s, with six datasets of chloride concentration available in areas along the northern and eastern coastline of Tampa Bay, which is nearest to historical dredging activities. Of these six, Kendall’s τ and p -values indicated increasing trends in chloride concentration for three datasets (TR 11-2, TR 10-2, 51), no trend in chloride concentration for two datasets (TR 9-3, 50), and a decreasing trend in chloride concentration for one dataset (TR 9-1). The upward trends in chloride concentration observed for TR 10-2 and 51 are likely the result of changes in local groundwater withdrawals. Well TR 11-2 had a gradual increasing trend in chloride concentration, fresh groundwater throughout the period of record, and a 3- to 4-foot increase in hydraulic head during the period of record, possibly caused by the construction and control of the Tampa Bypass Canal, resulting in changes to the regional potentiometric surface.

Florida↗
USGS2025

Rare earth element-mineralized carbonatite in the Bear Lodge Alkaline Complex, USA—Ore genesis implications from fluid inclusion characterization

Rare earth element (REE) resources of the Bear Lodge Alkaline Complex, Wyoming, are hosted in variably leached carbonatite spatially related to diatreme breccia pipes. We investigated the genesis of REE and lesser-known gold resources through fluid inclusion analysis of carbonatite, fluorite breccia, and smoky quartz vein samples. Physicochemical characteristics of inclusion-trapped fluids were evaluated using petrography, microthermometry, Raman spectroscopy, decrepitate mound analysis, energy-dispersive spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and noble gas isotope analysis. Microthermometry results reveal three fluid types that affected carbonatite dikes within deeper zones that escaped near-surface, ore-grade REE enrichment: (1) high-temperature (330–432°C) magmatic fluid captured in fine-grained calcite; (2) REE-enriched alkali bicarbonate-sulfate brine; and (3) low-temperature (117–182°C) diluted magmatic or meteoric water. Multiphase brine-like inclusions contain burbankite, nahcolite, strontianite, celestine and alkali sulfate daughter crystals, linking them to early burbankite mineralization. Peripheral smoky quartz and fluorite occurrences at Smith Ridge, 1.5 km from the central carbonatite dike swarm, contain primary inclusions that are Cl-poor and rich in Na-HCO 3 -SO 4 , similar to secondary and pseudosecondary inclusions in carbonatites. Helium isotopes reveal a MORB-like source for carbonatite samples and an older crust signature at Smith Ridge, consistent with the proximal ridge-top exposures of Archean granite. Results from this fluid inclusion study coupled with previous studies of carbonatite mineral paragenesis, show that light REEs (LREEs) were not mobilized great distances. Instead, burbankite crystallized within carbonatite from alkali bicarbonate fluids. With sodium retained in early burbankite, outward-emanating fluids enriched in potassium relative to sodium (higher K:Na) resulted in potassium–ferric iron metasomatism of silicate host rocks. This alkali fractionation was accompanied by fractionation of LREEs and heavy REEs (HREEs), with LREEs dominating the central carbonatite resources. In contrast, areas of peripheral REE mineralization at Bear Lodge are commonly characterized by higher HREE:LREE ratios. The K:Na ratio of associated fenites or alteration assemblages could be indicative of early crystallized burbankite in carbonatites and REE fractionation processes potentially leading to areas of concentrated HREEs with greater supply vulnerabilities.

South Dakota, Wyoming↗
USGS2026

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas↗