USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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,687 records · Page 23Linked to original sources
USGS1999

Water quality of Rob Roy Reservoir and Lake Owen, Albany County, and Granite Springs and Crystal Lake Reservoirs, Laramie County, Wyoming, 1997-98

The water quality of four reservoirs was assessed during 1997 and 1998 as a cooperative project between the Cheyenne Board of Public Utilities and the U. S. Geological Survey. The four reservoirs, Rob Roy, Lake Owen, Granite Springs, and Crystal Lake, provide approximately 75 percent of the public water supply for Cheyenne, Wyoming. Samples of water and bottom sediment were collected and analyzed for selected physical, chemical, and biological characteristics to provide data about the reservoirs. Water flows between the reservoirs through a series of pipelines and stream channels. The reservoirs differ in physical characteristics such as elevation, volume, and depth. Profiles of temperature, dissolved oxygen, specific conductance, and pH were examined. Three of the four reservoirs exhibited stratification during the summer. The profiles indicate that stratification develops in all reservoirs except Lake Owen. Stratification developed in Rob Roy, Granite Springs, and Crystal Lake Reservoirs by mid-July in 1998 and continued until September, with the thickness of the epilimnion increasing during that time. Secchi disk readings indicated Rob Roy Reservoir had the clearest water of the four reservoirs studied. The composition of the phytoplankton community was different in the upper two reservoirs from that in the lower two reservoirs. Many of the species found in Rob Roy Reservoir and Lake Owen are associated with oligotrophic, nutrient-poor conditions. In contrast, many of the species found in Granite Springs and Crystal Lake Reservoirs are associated with mesotrophic or eutrophic conditions. The total number of taxa identified also increased downstream. The chemical water type in the reservoirs was similar, but dissolved-solids concentrations were greater in the downstream reservoirs. Water in all four reservoirs was a calcium-bicarbonate type. In the fall of 1997, Rob Roy Reservoir had the lowest dissolved-solids concentration (19 milligrams per liter), whereas Crystal Lake Reservoir had the highest concentration (63 milligrams per liter). Relatively little differences in the concentrations of major-ion species were noted between samples collected near the surface and near the bottom of the same reservoir. In contrast, iron and manganese concentrations generally were higher in samples collected near the bottom of a reservoir than in near-surface samples collected from the same reservoir. Composite bottom-sediment samples from all four reservoirs contained similar concentrations of bulk constituents such as aluminum, iron, phosphorus and titanium, but varied in concentrations of trace elements. Trace-element concentrations in Rob Roy Reservoir and Lake Owen were similar to the crustal average, whereas in Granite Springs and Crystal Lake Reservoirs the concentrations were similar to granitic rocks.

Wyoming↗
USGS2025

Assessment of coastal and fluvial morphodynamic changes using Structure-for-Motion: A case study of the Sfȃntu Gheorghe Mouth (Danube Delta, Romania)

The ability to accurately map erosion, flooding, and habitat loss in coastal environments is crucial for formulating national strategies aimed at preventing and mitigating the impacts of natural disasters. A fundamental component of this process is the implementation of coastal morphodynamics monitoring through Structure-from-Motion (SfM) techniques, utilizing high-resolution 2D/3D data obtained from aerial photogrammetry. To assess morphodynamic changes over a three-year period (2022 – 2024), several SfM-based photogrammetric studies were conducted, each year, in the Romanian sector of the Danube-Black Sea coastal zone, specifically at the mouth of one of the Danube River distributaries (Sf Gheorghe branch) into the Black Sea, and along the left bank, near Sf Gheorghe locality, located within the Danube Delta Biosphere Reserve (DDBR). The essential equipment for aerial photogrammetry comprises Unmanned Aerial Vehicles (UAVs) and Global Navigation Satellite Systems (GNSS). In this study, the UAV used was a DJI Mavic 3T (Enterprise/Thermal) drone, complemented by two Trimble R12i and R4 GNSS systems, as well as approximately 10 Ground Control Points (GCPs). Data acquisition and processing were carried out using specialized photogrammetric software (Agisoft Metashape) along with various GIS tools (e.g., Blue Marble Geographics Global Mapper and ESRI ArcMap). The photogrammetric products generated for the study, as detailed in this paper, include Digital Elevation Models (DEMs), Digital Terrain Models (DTMs), orthomosaics (orthophotos), and others. At Sfântu Gheorghe beach, a comparison between 2023 and 2024 photogrammetric surveys revealed that the left bank of the Sf. Gheorghe Arm, at the river mouth into the Black Sea, suffered from a twist (erosion) of up to 64 metres. Additionally, on the selected perimetre (total area of 31,910 square meters ) from the beach and dune zone of Sf. Gheorghe, an area of up to 16,202 square meters was eroded between 2023 and 2024. This contrasts with the period between 2022 and 2023, during which deposition predominated. Erosion at the Danube mouths and the adjacent Black Sea coastline is driven by a complex interaction of natural and anthropogenic factors. Natural processes, including subsidence, sea-level rise, and episodic extreme storm events, contribute significantly to coastal dynamics. Meanwhile, human-induced factors, such as upstream hydrotechnical works that limits sediment transport, cutting of navigation canals, as well as the exacerbating effects of climate change, further accelerate erosion. The recent Structure-from-Motion (SfM) surveys provide essential quantitative data, enabling a detailed analysis of both short-term and long-term morphodynamic changes influenced by seasonal variations and extreme hydrometeorological events in this highly dynamic coastal system.

Sfȃntu Gheorghe Mouth (Danube Delta)↗
USGS2025

Gravity and magnetic surveys of the Skaergaard intrusion, East Greenland

Aeromagnetic and gravity surveys of the Skaergaard intrusion in East Greenland were carried out in July–August 1971 as part of a grant to the University of Oregon Center for Volcanology to refine the models of crystallization and differentiation of the intrusion, specifically to test whether the intrusion is underlain by dense rocks of a reservoir 20 kilometers (km) thick (referred to as a “hidden zone”). The Skaergaard intrusion is a source of platinum group elements that are critical mineral resources for many technologies, and because no new data have been collected these legacy datasets remain a valuable asset. The total-intensity aeromagnetic survey was flown in early July 1971 with a proton precession magnetometer at a constant barometric altitude of 1.5 km (5,000 feet) with a nominal line spacing of 1 km. Two gravimeters were used to acquire 168 stations of which 86 were at known altitudes (mainly sea level) and 82 had altitudes measured by altimetry in late July–August 1971. Finally, a north-south ground vertical-intensity magnetic traverse was completed across the intrusion together with collection of oriented hand specimens. The hand specimens were measured for remnant magnetization and density, along with density measurements of more specimens collected by expedition geologists for other purposes. The intrusion is composed of layered gabbro with extensive crystal fractionation that is dense and strongly reversely polarized. After terrain correction and standard Bouguer gravity reduction, the gravity anomaly dataset was corrected for all rock above sea level using the density measurements of the various zones of the intrusion and the topographic and geologic maps (variable density Bouguer gravity reduction). A large regional gradient in the gravity anomaly data was removed using orthogonal polynomial fitting to the gridded data. The zonal volumes of rock below sea level were calculated from the dipping polygonal layer gravity model of the intrusion below sea level and combined with elliptic cross–section cylinders for the various zones above sea level to approximate the original zonal volumes of the intrusion. The residual gravity anomaly of 18–20 milligals (mGal) was only about half of the expected anomaly if a large hidden zone proposed from petrologic considerations were present, and both two-dimensional and three-dimensional models imply that the exposed series of intrusion zones explain the gravity anomaly by their down-dip extension below sea level together with a small hidden-zone volume. A three-dimensional model of the exposed rocks and their down-dip extension below sea level also can account for the aeromagnetic anomaly with little or no requirement for hidden-zone rock. The middle and upper zone units of the intrusion contain the most magnetite and account for most of the aeromagnetic anomaly.

Skaergaard intrusion↗
USGS2026

Hyperspectral narrowband imaging spectroscopy: A new paradigm for Earth observation

This editorial introduces the Special Issue entitled “Hyperspectral Narrowband Imaging Spectroscopy: A New Paradigm for Earth Observation” in the August 2026 issue of Photogrammetric Engineering & Remote Sensing (PE&RS), the flagship journal of the American Society for Photogrammetry and Remote Sensing (ASPRS). This volume represents the fourth dedicated hyperspectral special issue published in PE&RS, following earlier contributions by Thenkabail et al. (2025, 2024a, 2024b), and continues ASPRS’s commitment to advancing cutting‑edge imaging spectroscopy research and its applications across Earth system science. Remote sensing is undergoing one of the most profound transformations in its history. The emergence of hyperspectral narrowband (HNB) imaging spectroscopy data, capable of acquiring hundreds of contiguous, narrow spectral bands, has shifted the discipline from observing Earth in a handful of broad spectral windows to capturing continuous spectral signatures of the Planet. This transition marks a decisive break from the multispectral paradigm that has dominated satellite remote sensing for nearly five decades, driven by the advent of new orbital imaging spectrometers such as EnMAP, PRISMA, and NASA’s EMIT, and by the forthcoming Surface Biology and Geology (SBG) mission ((Pires Silva et al., 2026; Bourriz et al., 2025; Thenkabail et al., 2025; Chabrillat et al., 2024; Aneece et al., 2024; Dave et al., 2024; Thenkabail et al., 2024a; Thenkabail et al., 2024b; Thenkabail, 2024a; Thenkabail, 2024b; Thompson et al., 2022; Kokaly et al., 2022; Aneece & Thenkabail, 2022; Cawse Nicholson et al., 2021; Guanter et al., 2021; Vangi et al., 2021; Thenkabail et al., 2021). These missions (e.g., Table 1) deliver unprecedented spectral fidelity, improved signal to noise ratios, and global coverage capabilities, enabling a new era of quantitative, spectroscopy based Earth observation. Where multispectral broadbands (MBBs) provide only a few discrete measurements along the electromagnetic spectrum, HNB systems deliver rich, diagnostic information that enables scientists to characterize Earth’s surface with unprecedented biochemical, biophysical, and structural detail (Figure 1a, 1b). The implications for environmental monitoring, agriculture, water resources, and mineral exploration are profound. Several overarching themes emerge: • Spectral fidelity matters. The ability to preserve subtle absorption features is essential for mineral mapping, vegetation trait retrieval, and biochemical modeling. • AI and deep learning are indispensable. From destriping to classification, modern analytics must be scalable, label‑efficient, and capable of exploiting the full spectral–spatial richness of HNB data. • Physics‑based and data‑driven approaches must converge. Radiative transfer models such as PROSAIL, enhanced with localized soil parameterizations, remain foundational for biophysical retrievals and model‑based inference. • Dimensionality reduction and feature extraction are critical. Techniques such as L1‑ISOMAP demonstrate that intelligent manifold learning can unlock the structure of fused, high‑dimensional datasets. • Next‑generation architectures must be interactive and multimodal. ICTNet exemplifies the future of hyperspectral classification: hybrid, synergistic, and capable of modeling both local textures and global spectral dependencies.

Photogrammetric Engineering and Remote Sensing (PE↗
USGS2025

Unlocking ecological insights from sub-seasonal visible-to-shortwave infrared imaging spectroscopy: The SHIFT campaign

We stand at the threshold of a transformative era in Earth observation, marked by space-borne visible-to-shortwave infrared (VSWIR) imaging spectrometers that promise consistent global observations of ecosystem function, phenology, and inter- and intra-annual change. However, the full value of repeat spectroscopy, the information embedded within different temporal scales, and the reliability of existing algorithms across diverse ecosystem types and vegetation phenophases have remained elusive due to the absence of suitable sub-seasonal spectroscopy data. In response, the Surface Biology and Geology (SBG) High-Frequency Time Series (SHIFT) campaign was initiated during late February 2022 in Santa Barbara County, California. SHIFT, designed to support NASA's SBG mission, addressed mission scoping, scientific advancement, applications development, and community building. This ambitious endeavor included weekly Airborne Visible InfraRed Imaging Spectrometer-Next Generation (AVIRIS-NG) imagery acquisitions for 13 weeks (spanning February 24 to May 29, 2022), accompanied by coordinated terrestrial vegetation and coastal aquatic data collection. We describe the rich datasets collected and illustrate how the complex sub-seasonal patterns of change can be linked to biological science and applications, surpassing insights from multispectral observations. Leveraging open-source processing methods and cloud-based analysis tools, the SHIFT campaign showcases the readiness of the scientific community to harness ecological insights from remotely sensed hyperspectral time series. We provide an overview of SHIFT's goals, data collections, preliminary results, and the collaborative efforts of early career scientists committed to unlocking the transformative potential of high-frequency time series data from space-borne VSWIR imaging spectrometers.

Ecosphere↗
USGS2025

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia↗
USGS1993

Ground-water quality near an inactive landfill and sludge-spreading area, Tallahassee, Florida

Ground-water quality of a 120-acre area encompassing an inactive landfill southwest of Tallahassee, Florida, was assessed by installing and sampling 21 monitoring wells. Sludge from a municipal sewage-treatment plant was also applied to the land surface near the landfill site. Water samples were collected from the 21 monitoring wells from June 1987 through September 1990 and analyzed for concentrations of major ions, nitrogen species, trace metals, field measurements, and more than 140 organic constituents. Water quality in the Upper Floridan aquifer at the study site was compared to background water quality in other areas of Leon County. Ground-water quality of the study site was evaluated by the location of wells sampled in relation to the landfill and the sludge-spreading areas. Results of nonparametric statistical analysis indicated that potassium and nitrite plus nitrate as nitrogen concentrations in water samples from wells completed in the Upper Floridan aquifer at the study site were significantly different at the 5 percent probability level from those in samples from background wells completed in this aquifer in Leon County. Median potassium concentrations were 0.4 mg/L (milligrams per liter) in water from wells at the study site and 0.7 mg/L in water from background wells in Leon County. Median nitrite plus nitrate concentrations were 6.48 mg/L in water samples from wells at the study site and 0.51 mg/L in water samples from background wells. Graphical comparison of the water quality characteristics in water from six categories of wells; upgradient, landfill, adjacent to the landfill, downgradient onsite, downgradient offsite, and background Leon County, indicated that sodium, bicarbonate, sulfate, iron, manganese, and dissolved solids concentrations and specific conductance values were highest within the landfill. Upgradient and downgradient well categories include wells in areas affected by sludge spreading. Localized reducing conditions were indicated in a downgradient area by median iron concentrations of 930 and 30 ug/L (micrograms per liter) in samples from wells only 5 feet apart. Nitrite plus nitrate concentrations were lowest in water from wells in the landfill compared to water from the other categories. Concentrations of most trace elements and organic compounds were mostly below detection limits, although State maximum contaminant levels of 1 ug/L for benzene and vinyl chloride and 3 ug/L for tetrachloroethene were exceeded in 14 samples. Concentrations of sodium, potassium, chloride, sulfate, nitrite plus nitrate, total organic and ammonia nitrogen, and ammonia nitrogen in ground water were statistically analyzed in areas affected by sludge-spreading and other areas. Only nitrite plus nitrate and chloride concentrations were significantly different between the two areas. Median nitrite plus nitrate and chloride concentrations in sludge-spreading areas were 6.9 mg/L and 2.9 mg/L, respectively, compared to 1.1 mg/L and 1.8 mg/L in other areas.

Florida↗
USGS1999

Simulation of the effects of development of the ground-water flow system of Long Island, New York

Extensive development on Long Island since the late 19th century and projections of increased urbanization and ground-water use makes effective water-resource management essential for preservation of the island's hydrologic environment and maintenance of a reliable source of water supply. This report presents results of a ground-water flow simulation analysis of the effects of development on the Long Island ground-water system. It describes ground-water levels, stream-flow, and the ground-water budget for the predevelopment period (pre-1900), the 1960's drought, and a more recent (1968-83) period with significant hydrologic stress. The report also presents estimated effects of a proposed water-supply strategy for the year 2020. Long Island has three major aquifers-the upper glacial (water-table), the Magothy, and the Lloyd aquifers-that are separated to varying degrees by confining units. Before development, recharge from precipitation entered the ground-water system at a rate of more than 1.1 billion gallons per day. An equal amount discharged to streams (41 percent), the shore (52 percent), and subsea boundaries (7 percent) . Urbanization and withdrawal of more than 400 Mgal/d (million gallons per day) from wells have resulted in local effects that include declines in ground-water levels, drying up and burial of streams and wetlands, reduction of ground-water recharge by increased overland flow to the ocean, a general decrease in ground-water discharge, and salt water intrusion. In some areas, the reduction in recharge is mitigated by leakage from water-supply and wastewater disposal lines, and infiltration of storm water through recharge basins. During 1968-83, a net loss of 240 Mgal/d from the ground-water system caused a decrease in ground-water discharge to streams (135 Mgal/d), to the shore (82 Mgal/d), and to subsea boundaries (23Mgal/d).The greatest adverse effects have been in western Long Island, where the most severe development has occurred. This analysis shows stream base flow to be highly sensitive to water-table fluctuations, and long streams to be more sensitive than short ones. A water-supply scenario for the year 2020 was simulated that employs redistribution of pumping centers to mitigate extreme local effects . Although the net stress on the ground-water system was projected to increase 57 Mgal/d (24 percent) above that of 1968-83, redistribution of ground-water withdrawals across the island would allow recovery of cones of depression in western Long Island, thereby reducing the threat of salt water intrusion and increasing base flow of some streams . The increased stress would cause a net decrease in base flow island wide of 44 Mgal/d; total base flow would be 281 Mgal/d - 39 percent below predevelopment levels or 14 percent below 1968-83 levels. The most severe effects would be in Nassau and western Suffolk Counties.

New York↗
USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗
USGS2026

Conventional hydrothermal power-producing systems of the Great Basin, USA

As part of the update to the electric-grade conventional hydrothermal assessment of the Great Basin, USA, Monte Carlo analyses of identified resources within explored regions will be performed to make estimates of discovered resources and associated uncertainty. Analyses use conditional statistics where estimates are conditioned upon a hydrothermal favorability map, allowing for the likelihood that more resources exist in regions of higher hydrothermal favorability. For these analyses, a dataset of identified hydrothermal systems is compiled, and the new compilation is described herein. Recognizing that a single hydrothermal system may be developed with multiple power plants, and that the hydrothermal upflow zone may be several kilometers across with many measurements characterizing a single hydrothermal system, a procedure was developed and employed to create clusters of points (power plants, measurements, etc.) that are associated with a single system, and a new central point was defined as the best estimator of the center of the hydrothermal system. Hydrothermal systems were uniquely identified by grouping electric-grade hydrothermal measurements and operating power plants within a distance of 10 km. Groups that are >10 km apart are assumed to be different electric-grade hydrothermal systems. While 10 km was used as the threshold, most systems were significantly further apart, and most points within groups were typically within 5 km of each other. A well measurement was considered an electric-grade measurement of a hydrothermal system if it had two properties: a measured temperature of >85 °C and evidence of hydrothermal convection. Other points that were added to the dataset are locations of operating powerplants or locations that have been classified as an electric-grade hydrothermal resource by either the U.S Geological Survey (USGS) or the Great Basin Center for Geothermal Energy. After all points are assigned to systems, new points were computed with the goal of identifying the center of the throat of the hydrothermal upflow zone. If operating powerplants exist for a system, then the arithmetic average of all power plant locations is used. Otherwise, if USGS made an estimate, that location is used. In the absence of both powerplants or USGS estimates, the arithmetic average of all electric-grade measurement locations is used. An example is shown of how these newly compiled locations might be ranked for uncertainty analyses, where higher confidence is assumed if measured temperature is higher and there are many supporting measurements indicating an electric-grade resource. In summary, 28 systems have operating power plants, an additional 78 systems are known identified electric-grade hydrothermal resources, and 99 new systems were identified as probable electric-grade systems with varying levels of confidence. These 205 locations are shown as a function of a recent hydrothermal favorability map, conceptually illustrating the conditional statistics that can be used to make estimates of the undiscovered resources of the Great Basin. An accompanying data release provides summaries of developed capacity by system and USGS estimates of likely total capacity and associated uncertainty.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS1995

Streamflow, ground-water recharge and discharge, and characteristics of surficial deposits in Buzzards Bay Basin, southeastern Massachusetts

Streamflow measurements at 14 low-flow partial-record stations and two discontinued streamflow-gaging stations were related to concurrent streamflows at six long-term streamflow-gaging stations to estimate streamflows at selected flow durations from 50 to 99 percent for water years 1967-91. At low flows, stream discharge per square mile generally increased with increasing percentage of subbasin underlain by stratified-drift deposits. At the 70-, 90-, and 99-percent flow durations, subbasins underlain primarily by stratified-drift deposits had stream discharges per square mile at least 2, 4, and 8 times greater, respectively, than subbasins underlain primarily by till and bedrock deposits. Streamflow measured three times at the Paskamanset River showed that streamflows downstream of municipal pumped wells increased at a lower rate per unit drainage area than streamflows upstream of the pumped wells. When pumpage from wells near the river was added to measured streamflow at the most downstream station, the sum was similar to the potential streamflow at the station if its streamflow per unit area was the same as that upstream from the pumped wells. Mean ground-water recharge and discharge rates were computed from continuous records of daily mean discharge during water years 1967-91 for six streamflow-gaging stations in southeastern Massachusetts and Rhode Island. Estimates of mean ground-water recharge were 19.7 to 22.6 inches per year for stations with drainage areas primarily underlain by till and bedrock deposits, and 23.8 to 25.2 inches per year for stations with drainage areas primarily underlain by stratified-drift deposits. During drought years, annual ground-water recharge can be less than 50 percent of the mean ground-water recharge rate for water years 1967-91. Estimates of ground-water discharge rates generally were 1 to 4 inches per year less than estimates of ground-water recharge rates. Ground-water discharge-duration curves were calculated from hydrograph separation of streamflow records at the six continuous streamflow-gaging stations for water years 1967-91. When streamflows were less than the median streamflow (50-percent duration), ground-water discharges averaged 85.7 percent of the total streamflow at the six streamflow-gaging stations. Ground-water discharges were estimated at the 14 low-flow partial-record stations and two discontinued streamflow-gaging stations by multiplying streamflows at the selected flow durations from 50 to 99 percent for water years 1967-91 by 85.7 percent. Data on surficial deposits in the basin were obtained by drilling 10 wells, by conducting 12 seismic-refraction surveys, and by collecting well log data for 313 wells from consulting firms and the Massachusetts Department of Environmental Management. These data indicated that a few minor corrections could be made to existing saturated thickness maps for stratified-drift deposits in Buzzards Bay Basin. Few well logs were detailed enough to estimate transmissivities; thus existing transmissivity maps were not updated for stratified-drift deposits in the basin.

Massachusetts↗
USGS2026

Multiple-well monitoring site adjacent to the Midway- Sunset and Buena Vista Oil Fields, Kern County, California

Groundwater quality in and around oil fields in the Southern San Joaquin Valley is of interest to many California residents that rely heavily on groundwater for domestic, commercial, and agricultural use. To help assess the effects of historical oil-field activities and natural geologic sources on groundwater near the southwest margins of the Kern County Groundwater Subbasin, a multiple-well monitoring site was installed near the administrative boundary between the Midway-Sunset and Buena Vista Oil Fields in Kern County, California. The installation of the Midway-Sunset Buena Vista multiple-well monitoring site (MSBV) supports regional analysis of the relations of oil and gas sources to groundwater quality by providing information about the geology, hydrology, geophysical properties, and water quality of the alluvial and upper Tulare aquifers in areas where groundwater data were limited. Data collected from the site included drill cuttings, whole core samples, sidewall core samples, mud-gas analysis, borehole geophysical logs, depth to water measurements, and water quality samples. Whole cores were scanned using dual energy computed tomography. Subsamples of selected cores were analyzed for density, porosity, specific retention, and bulk minerology. Thin sections of the subsamples were prepared, photographed, and examined. Two samples were analyzed using scanning electron microscope technology to examine the microporosity of diatomite laden sediment. Instrumentation installed in the wells collect hourly depth to water measurements. Analysis of the data show there is 355 feet of alluvium overlying the Tulare Formation at the well site. The contact between the two formations is an aquitard resulting in a perched aquifer in the alluvium and unconfined aquifer in the Tulare Formation. The alluvium is more heterogenous and finer grained than the Tulare Formation resulting in markedly higher porosity in the alluvium compared to the Tulare Formation. Higher specific retention observed in the alluvium is attributed to the finer grained sediment and greater abundance of reworked diatomite (as represented by opal-CT [cristobalite-tridymite]) compared to the Tulare Formation. Total dissolved solids (TDS) approached or exceeded 10,000 milligrams per liter (mg/L) in the alluvium from approximately 176 to 242 feet below land surface and at the top of the Amnicola clay at approximately 670 feet below land surface within the Tulare Formation. Elevated TDS, chloride, and boron concentrations in the alluvium and on top of the Amnicola clay likely reflect groundwater that is mixed with oil-field water. Water chemistry and modern-aged groundwater in the alluvial monitoring well (MSBV #3) are consistent with the oil-field water in the alluvium being derived from documented historical surface disposal of oil-field water upslope (northwest) of the site. Water chemistry and pre-modern groundwater age in the deeper Tulare monitoring well (MSBV #1) on top of the Amnicola clay are consistent with oil-field fluids derived from upslope natural geologic sources or old oil wells that leak in the subsurface. Shallow groundwater in the Tulare (MSBV #2) is not affected by mixing with oil-field sources.

EarthArXiv↗
USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗
USGS2024

Three-dimensional temperature maps of the Williston Basin, USA: Implications for deep hot sedimentary and enhanced geothermal resources

As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.

Montana, North Dakota, South Dakota↗
USGS1999

Environmental setting and water-quality issues in the lower Tennessee River basin

The goals of the National Water-Quality Assessment Program are to describe current water-quality conditions for a large part of the Nation's water resources, identify water-quality changes over time, and identify the primary natural and human factors that affect water quality. The lower Tennessee River Basin is one of 59 river basins selected for study. The water-quality assessment of the lower Tennessee River Basin study unit began in 1997. The lower Tennessee River Basin study unit encompasses an area of about 19,500 square miles and extends from Chattanooga, Tennessee, to Paducah, Kentucky. The study unit had a population of about 1.5 million people in 1995. The study unit was subdivided into subunits with relatively homogeneous geology and physiography. Subdivision of the study unit creates a framework to assess the effects of natural and cultural settings on water quality. Nine subunits were delineated in the study unit; their boundaries generally coincide with level III and level IV ecoregion boundaries. The nine subunits are the Coastal Plain, Transition, Western Highland Rim, Outer Nashville Basin, Inner Nashville Basin, Eastern Highland Rim, Plateau Escarpment and Valleys, Cumberland Plateau, and Valley and Ridge.The lower Tennessee River Basin consists of predominantly forest (51 percent) and agricultural land (40 percent). Activities related to agricultural land use, therefore, are the primary cultural factors likely to have a widespread effect on surface- and ground-water quality in the study unit. Inputs of total nitrogen and phosphorus from agricultural activities in 1992 were about 161,000 and 37,900 tons, respectively. About 3.7 million pounds (active ingredient) of pesticides was applied to crops in the lower Tennessee River Basin in 1992. State water-quality agencies identified nutrient enrichment and pathogens as water-quality issues affecting both surface and ground water in the lower Tennessee River Basin. Water-quality data collected by State and Federal agencies between 1980 and 1996 were summarized to characterize surface- and ground-water quality of the subunits with respect to these issues. Median concentrations of nitrogen species generally were less than 1 milligram per liter in surface and ground water in all subunits, and were highest throughout the subunits that had the largest percentages of agricultural land use. Median phosphorus concentrations also were less than 1 milligram per liter in all subunits. Phosphatic limestones present in two subunits had a larger effect on phosphorus concentrations in surface and ground water than did the amount of agricultural land use in these subunits. Median counts of fecal coliform were higher in surface water than in ground water in all subunits. The highest median counts in surface water were in the Valley and Ridge (7,500 colonies per 100 milliliters) and the Outer Nashville Basin subunits (5,000 colonies per 100 milliliters). Highest median counts in ground water were in the Inner and Outer Nashville Basin subunit. Natural setting likely has an important effect with respect to fecal contamination of surface and ground water in the lower Tennessee River Basin.

Alabama, Georgia, Kentucky, Mississippi, Tennessee↗
USGS2025

Cancer risk and estimated lithium exposure in drinking groundwater in the US

Importance Lithium is a naturally occurring element in drinking water and is commonly used as a mood-stabilizing medication. Although clinical studies have reported associations between receiving lithium treatment and reduced cancer risk among patients with bipolar disorder, to our knowledge, the association between environmental lithium exposure and cancer risk has never been studied in the general population. Objectives To evaluate the association between exposure to lithium in drinking groundwater and cancer risk in the general population. Design, Setting, and Participants This cohort study included participants with electronic health record and residential address information but without cancer history at baseline from the All of Us Research Program between May 31, 2017, and June 30, 2022. Participants were followed up until February 15, 2023. Statistical analysis was performed from September 2023 through October 2024. Exposure Lithium concentration in groundwater, based on kriging interpolation of publicly available US Geological Survey data on lithium concentration for 4700 wells across the contiguous US between May 12, 1999, and November 6, 2018. Main Outcome and Measures The main outcome was cancer diagnosis or condition, obtained from electronic health records. Stratified Cox proportional hazards regression models were used to estimate the hazard ratios (HRs) and 95% CIs for risk of cancer overall and individual cancer types for increasing quintiles of the estimated lithium exposure in drinking groundwater, adjusting for socioeconomic, behavioral, and neighborhood-level variables. The analysis was further conducted in the western and eastern halves of the US and restricted to long-term residents living at their current address for at least 3 years. Results A total of 252 178 participants were included (median age, 52 years [IQR, 36-64 years]; 60.1% female). The median follow-up time was 3.6 years (IQR, 3.0-4.3 years), and 7573 incident cancer cases were identified. Higher estimated lithium exposure was consistently associated with reduced cancer risk. Compared with the first (lowest) quintile of lithium exposure, the HR for all cancers was 0.49 (95% CI, 0.31-0.78) for the fourth quintile and 0.29 (95% CI, 0.15-0.55) for the fifth quintile. These associations were found for all cancer types investigated in both females and males, among long-term residents, and in both western and eastern states. For example, for the fifth vs first quintile of lithium exposure for all cancers, the HR was 0.17 (95% CI, 0.07-0.42) in females and 0.13 (95% CI, 0.04-0.38) in males; for long-term residents, the HR was 0.32 (95% CI, 0.15-0.66) in females and 0.24 (95% CI, 0.11-0.52) in males; and the HR was 0.01 (95% CI, 0.00-0.09) in western states and 0.34 (95% CI, 0.21-0.57) in eastern states. Conclusions and Relevance In this cohort study of 252 178 participants, estimated lithium exposure in drinking groundwater was associated with reduced cancer risk. Given the sparse evidence and unknown mechanisms of this association, follow-up investigation is warranted.

contiguous United States↗
USGS2026

Coastal barrier resilience and resistance: Analysis and metrics for characterizing coastal state

Barrier islands are shaped by a variety of short- and long-term environmental processes such as storms and relative sea-level rise. These islands, found along the estuarine-marine interface, provide ecosystem services including storm surge and wave attenuation, erosion protection to inland marshes, habitat for fish and wildlife, and recreation. Natural resource managers require actionable information on how barrier island resilience and resistance changes over time to understand how an island’s current state relates to past conditions and to inform restoration prioritization and implementation. The U.S. Geological Survey and The Water Institute collaborated on a study to develop indicators of resilience and resistance for barrier islands in Louisiana. Here, resilience captures island persistence on yearly to decadal time scales, and resistance captures persistence on event time scales of days to weeks. The indicators fall in two categories: Tier 1 Screening Metrics, that can be readily calculated from available data, are easily interpretable as an evaluation of barrier condition, and provide a high-level snapshot of overall barrier resilience and resistance; and Tier 2 Analysis Metrics, which are detailed metrics that required specialized analysis and interpretation and are more applicable to answering specific questions managers may have about barrier state. The research team derived Tier 1 resilience indicators from subaerial land and vegetation cover calculated from publicly available maps and products based on satellite imagery. By benchmarking the total land and vegetation extent against their respective historical maxima, this metric provides a snapshot of an island’s current state in the context of its long-term trajectory. The research team developed Tier 1 resistance indicators based on subaerial island configuration and water level recurrence as a proxy for evaluating island resistance to storms, which are the primary driver of short-term change. These Tier 1 metrics can be analyzed over time to provide a high-level assessment of how an island’s resistance decreases because of elevation loss or sea-level rise or increases due to restoration or natural recovery. The research team developed Tier 2 resilience and resistance indicators and associated analyses to provide detailed information for specific time periods or applications (e.g., wildlife management). These metrics include habitat coverage from high-resolution maps, which show composition changes over time to capture the evolving resilience of specific habitat types; high tide flooding analysis, which evaluate island area relative to specified flooding thresholds to characterize resistance in the short-term or, if analyzed over time, indicate changes in resilience; and hypsometric curve analysis, which allows managers to evaluate island area changes above their own elevation benchmarks of interest and similarly characterize resistance in the short-term or indicate changes in resilience if assessed over time. The research team calculated Tier 1 metrics of the barrier islands and headlands along the coast of Louisiana for the period of 1984 through 2021 and Tier 2 metrics for select times during that period depending on data available and quality. The results were captured in a report card for each barrier, which also includes an overview of the metrics and their interpretation; a restoration and storm history; and Tier 1 and Tier 2 metric analysis, including benchmarking against coastwide and regional values as well as to an island’s pre-restoration trajectory. These report cards provide a readily digestible synthesis of barrier condition and trajectory that coastal managers can use to support restoration prioritization and other decisions.

Louisiana↗
USGS2024

Hydrogeologic conceptual model of groundwater occurrence and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River is a major concern in the southwestern United States where the river provides water to about 40 million people for municipal and industrial use and is used to irrigate about 5.5 million acres of land. Much of the salinity in the Colorado River Basin is derived from natural interactions of surface water and groundwater with various geologic materials (rocks, soils, and alluvial deposits). The Dolores River in southwest Colorado is a major tributary of the Colorado River that historically accounts for about 6 percent of the salinity load to the Upper Colorado River Basin with the Paradox Valley being the primary source of salinity to the Dolores River. The Paradox Valley, one of several salt-anticline valleys in the region, is a fault-bounded topographic basin aligned with and exposing an underlying salt-anticline core. Salt deposits in the Pennsylvanian Paradox Formation of the Hermosa Group form an elongated salt diapir oriented northwest to southeast that is up to 12,000 feet (ft) thick beneath the present valley floor. Surface erosion, groundwater circulation, and weathering during Tertiary and Quaternary valley formation contributed to development of a cap rock, collapse features, breccia, and brine at the top of the exposed salt diapir. Today (2023), brine occurring in the brecciated cap rock and underlying salt deposits is in hydraulic connection with an overlying freshwater alluvial aquifer, and depending on seasonal river stage and hydrologic conditions, the brine discharges to the Dolores River causing the observed increase in salinity as the river crosses the Paradox Valley. To reduce salinity concentrations in the Dolores River, the Bureau of Reclamation (Reclamation) operates the Paradox Valley Unit (PVU). The PVU project consists of nine shallow brine pumping wells near the Dolores River and one deep disposal well where the brine is injected for disposal. When operational, the PVU pumping wells extract brine from the base of the alluvial aquifer that is piped and injected into a deep disposal well about 3 miles southwest of the PVU. The PVU became fully operational July 1, 1996, and by 2015, operation of the PVU had reduced salinity concentrations in the Dolores River by as much as 70 percent compared to pre-PVU conditions. In response to a 4.5 magnitude earthquake, injection operations, and thus PVU pumping, were ceased from March 2019 to June 2022. A trial period of PVU operation began in June 2022 with a reduced injection rate, and thus PVU pumping rate, of about two-thirds capacity to gather additional information and guide future operational decisions. In cooperation with Reclamation, the U.S. Geological Survey (USGS) developed this report to present the current (2023) understanding of groundwater and brine occurrence and discharge to the Dolores River in the Paradox Valley. Results from the compilation of spatial datasets, groundwater sampling and age dating, and aquifer tests are presented to provide improved understanding of the Paradox Valley hydrogeology, to supply datasets for a numerical groundwater-flow and brine-transport model, and to support future operations of the PVU. The hydrogeologic data provided herein, along with the most recent loading analysis for the Dolores River in the Paradox Valley, and a previous conceptual model for brine discharge to the river are used to present a conceptual understanding of groundwater occurrence in the Paradox Valley.

Colorado↗