Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Carbon Capture Science and Technology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,661 records · Page 15Linked to original sources
USGS2026

Seasonal and hydrologic variation influences habitat and functional structure of stream fish assemblages

Introduction: Hydrologic variability is a key driver of ecological structure in lotic systems, shaping habitat conditions, taxonomic diversity, and the functional traits that mediate species’ persistence and performance (e.g., reproductive success). While many studies examine taxonomic responses to variation in flows, few evaluate how spatiotemporal hydrologic variation influences the functional organization within stream fish communities. Methods: We quantified seasonal habitat structure and functional trait diversity of fish assemblages across six Ozark Plateau headwater streams representing two contrasting flow regimes: Groundwater Flashy and Runoff/Intermittent Flashy. Fish and habitat data were collected seasonally during a dry year (2002) and a wet year (2003). Functional space was constructed using PCoA of morphological, ecological, and life-history traits, and functional diversity was measured using community weighted means (CWMs), functional richness (FRic), functional evenness (FEve), and functional divergence (FDiv). Results: We found that habitat structure differed strongly by flow regime and season, with Runoff/Intermittent streams exhibiting pronounced reductions in depth, area, and velocity, while groundwater streams remained structurally stable. Functional identity of assemblages was similar across flow regimes, dominated by benthic, hydrodynamic taxa with opportunistic and periodic life-history strategies. However, functional structure differed significantly: FEve and FDiv were consistently lower in Runoff/Intermittent Flashy streams in both years, indicating assemblage dominance of species with similar trait combinations and reduced trait partitioning under variable flow. FRic and taxonomic richness remained stable across seasons and flow regimes, suggesting high functional redundancy despite species turnover. Discussion: Together, results show that flow regime mediates both habitat structural stability and functional organization. As climatic warming and extreme drought increase hydrologic instability in headwaters, functional trait approaches provide a sensitive tool for detecting losses of functional roles that may not be evident by using taxonomic metrics alone.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2025

From hydrated silica to quartz: Potential hydrothermal precipitates found in Jezero crater, Mars

On Earth, silica-rich phases from opal to quartz are important indicators and tracers of geological processes. Hydrated silica, such as opal, is a particularly good matrix for the preservation of molecular and macroscopic biosignatures. Cherts, a type of silica-dominated rocks, provide a unique archive of ancient terrestrial life while quartz is the emblematic mineral of the Earth's continental crust. On Mars, hydrated silica has been detected in several locations based on remote sensing and rover-based studies. In the present article we report on the detection of cobbles made of hydrated silica (opal or chalcedony), as well as well-crystallized quartz. These detections were made with the SuperCam instrument onboard Perseverance (Mars 2020 mission), using a combination of LIBS, infrared and Raman spectroscopy. Quartz-dominated stones are detected unambiguously for the first time on the Martian surface, and based on grain size and crystallinity are proposed to be of hydrothermal origin. Although these rocks were all found as float, we propose that these detections are part of a common hydrothermal system, and represent different depths / temperatures of precipitation. This attests that hydrothermal processes were active in and around Jezero crater, possibly triggered by the Jezero crater-forming impact. These silica-rich rocks, in particular opaline silica, are very promising targets for sampling and return to Earth given their high biosignature preservation potential.

Earth and Planetary Science Letters↗
USGS2026

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska↗
USGS2025

RioM-1: A new calcite reference material for U-Pb LA-ICP-MS geochronology

Determining absolute ages of carbonate diagenesis, faulting, fossil formation, speleothem growth, carbonate-hosted hydrocarbon deposits, vein mineralisation and hydrothermal alteration has become increasingly accessible through LA-ICP-MS U-Pb dating of calcite, complementing traditional isotope dilution methods still applicable to certain materials (e.g., speleothems via micro-drilling). However, well-calibrated reference materials for LA-ICP-MS calcite geochronology remain scarce. Here, we characterise the Rio Maior calcite, designated ‘RioM-1’, as a potential reference material for U-Pb dating by LA-ICP-MS. Fragments (0.1 to 1 cm 3 ) from a single scalenohedral crystal were analysed by ID-TIMS (U-Pb), LA-ICP-MS (U-Pb and 87 Sr/ 86 Sr), and SIMS (O isotopes). RioM-1 displays high U mass fraction and low, though variable, proportions of common Pb. Combined ID-TIMS analyses from two independent laboratories yielded a Tera-Wasserburg lower intercept date of 63.93 ± 0.11 Ma (2 s , MSWD = 1.3, n = 16). LA-ICP-MS U-Pb data from three independent laboratories are concordant with the TIMS age, producing a pooled date of 64.10 ± 0.12/1.2 Ma (2 s , MSWD = 7.7, n = 708) and an initial 207 Pb/ 206 Pb ratio of 0.85 ± 0.01. Other calcite reference materials, when normalised to RioM-1, yielded dates within uncertainty of their published values. SIMS measurements returned a mean 18 O/ 16 O of 0.002044450 ± 181 (1 s ) and δ 18 O SMOW of 19.57 ± 0.92‰ (1 s ), while LA-ICP-MS strontium isotope measurement yielded a mean 87 Sr/ 86 Sr of 0.708177 ± 9 (2 s ).

Geostandards and Geoanalytical Research↗
USGS2024

Fish beta diversity associated with hydrologic and anthropogenic disturbance gradients in contrasting stream flow regimes

Understanding the role of hydrologic variation in structuring aquatic communities is crucial for successful conservation and sustainable management of native freshwater biodiversity. Partitioning beta diversity into the additive components of spatial turnover and nestedness can provide insight into the forces driving variability in fish assemblages across stream flow regimes. We examined stream fish beta diversity across hydrologic and anthropogenic disturbance gradients using long-term (1916–2016) site occurrence records ( n = 17,375) encompassing 252 species. We assessed total beta diversity (Sørensen dissimilarity), spatial turnover, and nestedness of fish assemblages in contrasting stream flow regimes across a gradient of decreasing flow stability: groundwater stable ( n = 77), groundwater ( n = 67), groundwater flashy ( n = 175), perennial runoff ( n = 141), runoff flashy ( n = 255), and intermittent ( n = 63) streams. Differences in total beta diversity among the stream flow regimes were driven predominantly (>86 %) by spatial turnover (i.e. species replacement) as opposed to nestedness (i.e. species loss or gain). Total fish beta diversity and spatial turnover were highest in streams with intermediate flow stability (groundwater flashy), while more flow-stable streams (groundwater stable and groundwater) had lower turnover and higher nestedness. Species turnover was also strongly associated with seasonal variation in hydrology across all flow regimes, but these relationships were most evident for assemblages in intermittent streams. Distance-based statistical comparisons showed significant correlations between beta diversity and anthropogenic disturbance variables, including dam density, dam storage volume and water withdrawals in catchments of groundwater stable streams, while hydrologic variables were more strongly correlated with beta diversity in streams with runoff-dominated and flashy flow regimes. The high spatial turnover of species implies that fish conservation actions would benefit from watershed-focused approaches targeting multiple streams with wide spatial distribution, as opposed to simply focusing on preserving sites with the greatest number of species.

Science of the Total Envionrment↗
USGS1976

Water-quality study of a reach of the Merced River in Yosemite National Park and vicinity, California, April 1973 through September 1974

A 25-mile (40-kilometre) reach of the Merced River was studied to provide the National Park Service with water-quality data prior to the operation of a new sewage-treatment plant. Data were collected periodically at four Merced River and two inflow stations. Discharge varied seasonally with highest flow occurring at the time of spring snowmelt runoff. Lowest flow occurred during late summer, coinciding with high park visitor use and elevated water temperature. Water-quality variables investigated were discharge, temperature, major chemical constituents, selected trace metals, specific conductance, suspended sediment, turbidity, dissolved oxygen, pH, alkalinity, carbon dioxide, total organic carbon, nitrogen, phosphorus, algal growth potential, coliform bacteria, periphyton, and benthic invertebrates. Diel measurements were made at some stations to delineate daily fluctuations in selected water-quality variables. Results indicate that water in the reach of the Merced River studied was of good quality. The river had low concentrations of dissolved solids (8 to 41 milligrams per litre) and low alkalinity values (2 to 24 milligrams per litre of alkalinity as CaCO 3 ). Mean concentrations of nitrogen and phosphorus were highest in the Merced River at Rancheria Flat. Algal growth potential tests showed that phosphorus was limiting algal growth at Happy Isles Bridge, El Capitan Bridge, and South Fork Merced River, whereas nitrogen was limiting growth at Rancheria Flat. Input of excessive nitrogen and phosphorus concentrations in the Merced River during the summer low-flow period, and mixing with warm tributary inflow could cause excessive algal growth and production. The algal types found were diatoms, green algae, and blue-green algae. Biomass/chlorophyll a ratios in the Merced River were low (0.3 to 7.7). Benthic invertebrates consisted chiefly of Diptera, Ephemeroptera, and Plecoptera.

California↗
USGS1976

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
USGS2024

Deformation by pressure solution and grain boundary sliding in a retrograde shear zone in southern New England, USA

Alleghanian phyllonites in a shear zone in southern New England were formed by the retrogression and hydration of a high-grade Acadian pelitic schist. The retrogression was locally incomplete, resulting in both heterogeneous mineralogy and mineral compositions, and in many arrested reaction textures. These features, documented by backscattered electron and element map imaging, and by 40 Ar/ 39 Ar dating, provide a record of the processes active as the schist was converted to a phyllonite. On the margins of the shear zone relic biotite and garnet grains survive but are partially replaced by chlorite. Relic high-grade muscovite flakes containing up to 35 % paragonite (pg) also persist but are entrained in swarms of pg-poor muscovite flakes producing folia with an anastomosing network of muscovite grains of variable composition. Single crystals of pg-rich muscovite are truncated by muscovite in successively younger cross-cutting folia that contain decreasing pg contents, falling to <5% pg in the youngest folia. In the core of the shear zone all high-grade minerals are destroyed, but recrystallization continues with earlier greenschist facies folia replaced by lower grade and pg-poorer muscovite in younger folia. Muscovite grains in truncated folia are commonly kinked and selective recrystallization of muscovite in kink bands to pg-poor-compositions demonstrates that strain energy helped drive recrystallization. The boundaries of these new pg-poor muscovite folia are sharp, and truncations indicate that grain boundary sliding (GBS) involved dissolution. GBS also enabled folding by flexural slip along muscovite grain boundaries. Fold amplitudes increased as chlorite dissolved from fold limbs precipitated in fold hinges producing crescent-shaped saddle-reef-like structures. Together these observations of truncation and replacement by dissolution and precipitation demonstrate that pressure solution and GBS facilitated both the formation of the phyllonitic shear zone and slip along it.

Connecticut↗
USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗
USGS2026

Spatio-temporal modeling for assessing geoenergy resources: A workflow applied to gas in place variation in coal beds

The ability to estimate spatio-temporal changes in hydrocarbon reservoir properties and energy resources within pore volumes is essential for optimizing production, reservoir management, geologic energy storage, and safety in underground mining operations. In coal seams, predicting remaining methane gas-in-place (GIP) is critical for quantifying producible gas and improving mine safety and productivity through effective ventilation planning. Although such changes are commonly evaluated using physics-based numerical simulation models, these approaches often require extensive data, calibration effort, and time. This study presents a spatio-temporal geostatistical modeling approach that bridges the gap between purely spatial models and full numerical simulations. The method is applied to a case study of coal seam degasification in the Mary Lee coal group, Black Warrior Basin, Alabama, USA, to estimate GIP evolution over time within a selected mining district. The analysis uses published data from prior natural gas production history-matching of degasification using vertical wells. Empirical spatial and temporal statistics were calculated for reservoir pressure and water saturation, and spatio-temporal variogram models were fitted to experimental variograms. These models provided the structural basis for spatio-temporal kriging, integrated with spatial estimates of time-invariant parameters (porosity, density, and thickness) to estimate GIP. This approach enabled estimation of GIP changes over time, including periods without data. Boxplots of GIP estimates indicated systematic depletion and decreasing spatial variability, reflecting the impacts of degasification. Comparison with cumulative gas production from empirical well records showed approximately 85% agreement based on a relative similarity metric. Spatio-temporal GIP estimates were also used to estimate methane emissions to longwall ventilation systems and compared with reported emissions from the U.S. EPA Greenhouse Gas Reporting Program, showing similar distributions (≈80%) given data limitations. Overall, this integrated modeling approach provides time-dependent GIP estimates with broader implications for resource assessment applications.

Alabama↗
USGS2026

Evaluation of nutrient, alkalinity, and acid-neutralizing capacity stabilities in water samples analyzed by the U.S. Geological Survey National Water Quality Laboratory, 2023–24

The U.S. Geological Survey evaluated the stability of water-sample chemical analysis of nutrient, alkalinity, and acid-neutralizing capacity constituents with respect to the duration between sample collection and laboratory analysis, also known as the sample holding time. A study began in the spring of 2023 to evaluate the sample stability, between 2 and 180 days after sample collection, of the chemical properties and chemical constituents of alkalinity as calcium carbonate, filtered; acid-neutralizing capacity as calcium carbonate, unfiltered; total ammonia as nitrogen, filtered; total ammonia plus organic nitrogen as nitrogen, filtered and unfiltered; nitrite as nitrogen, filtered; nitrate plus nitrite as nitrogen, filtered; total nitrogen, filtered and unfiltered; orthophosphate as phosphorous, filtered; and total phosphorus as phosphorus (filtered and unfiltered) in water. Both surface water and groundwater matrices were represented. Sample instability varied by observed property and matrix; therefore, providing general guidance for sample holding time is not possible based on matrices alone. No correlations between field measurements of sample characteristics and sample instability were observed. Although observations for some properties indicate sample stability that exceeds the recognized U.S. Geological Survey National Water Quality Laboratory method holding times, this is not necessarily the case for matrices and seasonal characteristics that were not investigated. Based on the limited number of six sample sources used in this study, some patterns emerge for the 12 observed properties studied. Five observed properties generally indicate stability for as many as 180 days after sampling (total nitrogen as nitrogen, both filtered and unfiltered; orthophosphate as phosphorus, filtered; and phosphorus as phosphorus, both filtered and unfiltered). Other observed properties indicate stability for as many as 180 days for some matrices, but not for others. Finally, some observed properties indicate instability well before 180 days.

Scientific Investigations Report↗
USGS2024

Preliminary observations of the April 5th, 2024, Mw4.8 New Jersey earthquake

On 5 April 2024, 10:23 a.m. local time, a moment magnitude 4.8 earthquake struck Tewksbury Township, New Jersey, about 65 km west of New York City. Millions of people from Virginia to Maine and beyond felt the ground shaking, resulting in the largest number (>180,000) of U.S. Geological Survey (USGS) “Did You Feel It?” reports of any earthquake. A team deployed by the Geotechnical Extreme Events Reconnaissance Association and the National Institute of Standards and Technology documented structural and nonstructural damage, including substantial damage to a historic masonry building in Lebanon, New Jersey. The USGS National Earthquake Information Center reported a focal depth of about 5 km, consistent with a lack of signal in Interferometric Synthetic Aperture Radar data. The focal mechanism solution is strike slip with a substantial thrust component. Neither mechanism’s nodal plane is parallel to the primary northeast trend of geologic discontinuities and mapped faults in the region, including the Ramapo fault. However, many of the relocated aftershocks, for which locations were augmented by temporary seismic deployments, form a cluster that parallels the general northeast trend of the faults. The aftershocks lie near the Tewksbury fault, north of the Ramapo fault.

New Jersey↗
USGS2026

Evaluating groundwater quality influences from oil field operations and other anthropogenic activities in an urban setting, Santa Fe Springs, California

Groundwater quality is often affected by anthropogenic activities in urban settings. This study examines groundwater quality in and around the Santa Fe Springs Oil Field in Los Angeles County, California, where oil and gas production commonly intersects with high density industrial, commercial and residential land uses. Utilizing a combination of new and historical data, we evaluated potential pathways that would allow for oil field formation fluids to migrate into groundwater and whether mixing may have occurred based on the distribution of groundwater and oil field formation fluid tracers in samples. Samples were analyzed for a wide array of constituents including volatile organic compounds, light hydrocarbons, major ions, and various isotopic compositions. Despite evidence of oil field infrastructure providing potential pathways of migration via uncemented annular spaces, casing breaches and historical disposal of oil field formation water in surface ponds, the distribution and occurrence of stable isotopes of water, chloride, boron, and total dissolved solids do not indicate mixing of oil field formation water and groundwater. However, methane isotopic signatures and the presence of heavier alkanes suggest gas from oil-bearing formations have migrated from depth via oil field well infrastructure. Volatile organic compound detections were mainly from manufactured compounds unrelated to oil and gas production, with a relatively limited number of petroleum hydrocarbons also detected. Volatile organic compounds were generally found in wells tapping shallow, modern aged groundwater, indicating anthropogenic activities occurring at or near land surface as the source. Study results suggest that while oil field infrastructure provides migration pathways for oil field formation fluids to be introduced into groundwater, urban land uses not related to oil and gas production are the primary drivers of groundwater quality degradation.

California↗
USGS2025

Assessing nonpoint-source uranium pollution in an irrigated stream-aquifer system

Uranium (U) in rocks and soils of arid and semi-arid environments can be mobilized by irrigation and fertilization, posing environmental and health risks. Elevated U, along with selenium (Se) and nitrate (NO 3 ) co-constituents, necessitates careful monitoring and management. We developed a distributed-parameter numerical model to assess U pollution in an irrigated stream-aquifer system, applying it to a 552 km 2 region in Colorado's Lower Arkansas River Valley (LARV) over 14 years. A MODFLOW model, describing groundwater and stream flow, was coupled with an RT3D-OTIS model to portray reactive U transport. Calibration using the PESTPP-iES iterative ensemble smoother (iES) software indicated good agreement with observed U concentrations. The model revealed substantial and variable U levels across the LARV, highlighting potential hotspots and possible contributing factors, such as geological composition of the bedrock and near-surface shale and aquifer sediments derived from them, irrigation practices, and riparian landscape. U levels exceed the chronic standard (85th percentile = 30 μg/L, set by the US Environmental Protection Agency), which is the permissible regulatory threshold, in groundwater across 44 % of the region and along the river by an average factor of 2.9. Simulated average U concentrations in the non-riparian aquifer and river are 124 μg/L and 60 μg/L, respectively, compared with 112 μg/L and 62 μg/L for measured values. The average 85th percentile U concentration is 222 μg/L in the aquifer and 82 μg/L in the river. Average simulated U mass loading to the river is 0.17 kg/day per km, compared to an estimated 0.23 kg/day per km. Findings provide a baseline for comparing future simulated outcomes of alternative best management practices (BMPs) for U pollution mitigation and offer a methodology applicable to other irrigated regions.

Colorado↗
USGS2025

Understanding market sensitivity: Estimation of supply and demand elasticities for non-fuel minerals

In today’s rapidly changing economic landscape, understanding market responsiveness to price changes and the factors influencing commodity prices has become increasingly relevant. Price elasticities serve as indicators of how variations in market conditions affect supply and demand, providing insights into the sensitivity of commodity markets to price fluctuations. This paper presents a comprehensive analysis of price elasticities of supply and demand for 74 non-fuel mineral commodities including precious metals, base metals, minor metals, and industrial minerals that are utilized across various industries. We employ various econometric techniques, including fixed effects models for panel data and two-stage dynamic ordinary least squares (2S-DOLS) alongside autoregressive distributed lag (ARDL) models for time series analysis, to derive robust estimates of price elasticities. Our findings reveal variability in elasticities among different commodities and indicate that all studied mineral commodities exhibit price inelastic supply and demand in the short run, which we define as one year for the purposes of our analysis, given that the data is all annual. This research provides original estimates of price responsiveness for a wide range of commodities that have not been previously addressed in the literature, thereby enhancing the understanding of market dynamics in the mineral sector. Given that price elasticities can be influenced by factors such as market structure, technological advancements, mining costs, and industry-specific demand drivers, we use variables that serve as proxies for these factors.

Mineral Economics↗
USGS2024

Produced water geochemistry from hydraulically stimulated Niobrara Formation petroleum wells: Origin of salinity and temporal perspectives on treatment and reuse

Produced water (i.e., a mixture of returned injection fluids and geologic formation brines) represents the largest volumetric waste stream associated with petroleum production in the United States. As such, produced water has been the focus of intense study with emphasis on understanding the geologic origin of the fluids, environmental impacts of unintended or intentional release, disposal concerns, and their commodity (e.g., lithium) potential. However, produced water geochemistry from many active petroleum plays remain poorly constrained leading to knowledge gaps associated with the origin of brine salinity and parameters (e.g., radium levels) that can impact treatment, disposal, and possible reuse. Here we evaluate the major ion geochemistry, radium concentrations, and stable water isotope composition of ~120 produced water samples collected from 17 producing unconventional petroleum wells in Weld County, Colorado from the Late Cretaceous Niobrara Formation. This sample set encompasses eight produced water time series from four new wells across production days 0 to ~365 and from four established wells across production days ~1000 to ~1700. Additionally, produced water from nine other established Niobrara Formation wells were sampled at discrete time points ranging from day 458 to day 2256, as well as hydraulic fracturing input fluids. These results expand the available Niobrara Formation produced water geochemical data, previously limited to few wells sampled within the first year of production, allowing for the heterogeneity of major ions and radium to be evaluated. Furthermore, we explore the geochemical relationships between major ion ratios and stable water isotope composition to understand the origin of salinity in Niobrara Formation brines from the Denver-Julesburg Basin. These findings are discussed with perspective toward potential treatment and reuse of Niobrara produced water prior to disposal.

Colorado, Nebraska, Wyoming↗
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↗