Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geochemical Exploration”

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,228 records · Page 9Linked to original sources
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

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

conterminous United States↗
USGS2026

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗
USGS2026

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2026

Uranium—Deposits, production and resources, market dynamics, and supply chain risks

Introduction Interest in nuclear power for the generation of electricity has risen with the increase in the need for more diverse baseload power, enhanced energy security, and the development of new technologies, such as small modular reactors (SMRs), which could provide power for remote areas, industrial applications, and artificial intelligence (AI) data centers. In 2024, the U.S. Department of Energy received $2.7 billion in congressional funding to bolster the domestic uranium production and nuclear fuel supply chain and address reliance on imports from foreign suppliers. In 2025, the U.S. Government issued several Executive and Secretary’s orders aimed at revitalizing the U.S. nuclear sector. If SMRs are to be as widely deployed in the United States and worldwide as envisioned, demand for uranium (nuclear reactor fuel) will likely increase. After the Fukushima nuclear accident in 2011, the market spot price of uranium began a decline, followed by a decrease in U.S. and global uranium exploration and mine development expenditures that led to a uranium supply deficit until 2020, when prices started to recover, prompting a resurgence in uranium exploration and development. In January of 2024, the uranium spot price rose to a 17-year high $106 (U.S. dollars) per pound of U3O8 (triuranium oxide, commonly known as “yellowcake”), which is expected to increase uranium exploration, mine development, and uranium production domestically and worldwide.

Fact Sheet↗
USGS2026

Validating antler- and movement-based estimates of caribou reproduction using video camera collars

Wildlife agencies invest substantial resources in monitoring ungulate reproductive rates given ongoing concerns about population trends and their implications for management. For barren-ground caribou ( Rangifer tarandus ) in the North American Arctic, data on parturition and neonate survival have traditionally been collected using visual observations of antler retention and calf presence from small, fixed-wing aircraft. However, these surveys are weather-dependent, expensive, and dangerous, motivating the development of movement-based approaches that infer reproductive events from global positioning system (GPS) location data. We had the unique opportunity to use video camera collars to validate both antler- and movement-based methods for estimating reproduction in barren-ground caribou, using data from the Porcupine (2018–2023) and Western Arctic (2021, 2023) herds in Alaska, USA, and Yukon, Canada. For the movement-based approaches, we assessed individual-based and population-based methods, which both inferred parturition following a sharp decline in movement and inferred calf loss following a sharp increase in movement. We investigated these approaches using rarified 1-, 2-, 4-, and 8-hour fix intervals. Based on video collar observations, we found that antler retention at the onset of the calving period was highly accurate for predicting parturition (96% accuracy) but became unreliable as calving progressed, as approximately 80% of parturient females shed their antlers during the calving period. In contrast, both movement methods performed marginally for estimating caribou parturition (~70% accuracy), with estimates from the individual-based method exhibiting greater bias than those from the population-based method. The individual-based method also poorly predicted neonate survival (≤49% accuracy for detecting parturition and calf fate for the Porcupine Herd). Video data revealed that inaccuracies with the movement methods often occurred when caribou bedded during storms, rested after traversing mountainous terrain, lost a calf shortly after birth, or increased movement to evade insects. These results have important implications for management agencies using antler- and movement-based methods to estimate caribou reproduction, particularly when informing assessments of population decline or recovery. Our results demonstrate the utility of video collar data for validating reproductive monitoring approaches and provide insights into how these approaches can be improved.

Alaska, Yukon↗
USGS2026

Overstorey mortality promotes juvenile piñon pine growth during favourable weather at cooler, wetter sites

Hotter droughts have resulted in widespread tree die-off events globally, frequently leading to regeneration failure. Dry forest recovery often depends on the growth and survival of extant juvenile trees. However, it is unclear how microenvironmental changes following overstorey tree die-off affect juvenile trees, particularly in dryland systems where tree recruitment is typically limited by water availability and heat stress. We simulated an overstorey tree die-off event by girdling trees in piñon-juniper woodlands across the south-western United States. We sampled juvenile piñon pine growth from live and dead overstorey treatments across six study sites spanning a regional latitudinal gradient and local elevational gradients. We examined how juvenile branch and needle growth differed between live and dead overstorey treatments, and whether responses varied with weather conditions and juvenile tree size following overstorey mortality. We found greater juvenile branch and needle growth under dead compared with live overstorey trees for 2 years following overstorey mortality at mid- and high-elevation sites which are typically cooler and wetter than the other sites. These observed growth releases were contingent on favourable post-mortality weather conditions. Higher growth under dead overstorey occurred at sites experiencing near-average climatic water deficits compared with sites experiencing above-average climatic water deficits. Growth at all sites increased from the first to second year after overstorey mortality. Across sites, growth was unrelated to juvenile tree size. Synthesis . Our results underscore differentiation in juvenile responses to overstorey tree die-off driven by local site conditions and weather across the range of Pinus edulis . Overstorey mortality resulted in consistently higher juvenile growth only at climatically favourable sites and during favourable weather, while unmeasured microsite differences could help account for variation observed at the hottest and driest site. Results from less climatically favourable sites suggest that overstorey trees neither directly limit nor facilitate juvenile growth, though further study over longer timeframes is needed to resolve the pace and magnitude of potential recovery or decline. Overall, juvenile vigour may be promoted following overstorey mortality only in a narrow spatial (site) and temporal (weather) environmental context, suggesting additional vulnerabilities for piñon populations under more arid conditions.

Arizona, Colorado↗
USGS2025

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

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

South Dakota, Wyoming↗
USGS2026

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

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

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS2026

Linking distribution and return-on-investment models to optimize woody management for prairie grouse in Nebraska

Grasslands in Nebraska, USA, face threats from agricultural conversion, urban development, and woody encroachment, all of which negatively affect prairie grouse ( Tympanuchus spp.) populations. To optimize conservation planning, Nebraska wildlife agencies developed probabilistic area-based surveys for greater prairie-chicken ( T. cupido ) and sharp-tailed grouse ( T. phasianellus ) to sample landscapes across a range of environmental conditions. This design improves historical surveys and enables the development of distribution models that quantitatively define habitat associations and support scenario-based conservation planning. Using survey data collected during 2020–2022, we modeled prairie grouse occurrence and abundance as functions of land cover, topography, and climate using Bayesian logistic and zero-inflated negative binomial models with regularized horseshoe priors. We then conducted a maximum potential return-on-investment analysis of woody cover treatments, assuming sustained treatment success, relative to projected impacts of woody encroachment on prairie grouse populations by 2050. Among modeled associations were a positive association with grasslands having low woody cover and a negative association with grasslands having high woody cover. Across the 3-year period, median estimated annual populations were 142,380 for greater prairie-chicken (range of 95% CIs across years = 68,821–277,615) and 64,154 for sharp-tailed grouse (range of 95% CIs across years = 27,550–146,559). Under projected woody encroachment, mean predicted population declines were 10% for greater prairie-chicken (range of 95% CIs = 7–14%) and 6% for sharp-tailed grouse (range of 95% CIs = 5–7%). Areas with high prairie grouse density and low treatment costs in 2021, and high projected woody encroachment and population loss by 2050, offered the greatest return on investment for woody management. Return on investment was greatest in the northwestern Shortgrass Prairie ecoregion (northwestern Nebraska) for sharp-tailed grouse and the eastern Sandhills ecoregion (central Nebraska) for both species. These models underscore the value of evidence-based, quantitative approaches for prioritizing conservation actions on working lands. Scenario-based modeling could be extended to guide other treatments, such as optimizing restoration (e.g., Conservation Reserve Program) or incentivizing grassland persistence in areas with predicted climate resilience.

Nebraska↗
USGS2026

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li 2 O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li 2 O. At 90% confidence, the region contains at least 130 kt Li 2 O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li 2 O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Alabama, Georgia, Maryland, North Carolina, South ↗
USGS2026

Evolution of arsenic speciation during thermal treatment simulating wildfire heating in arsenic-rich sediments

Understanding arsenic transformations during wildfires can help better constrain the environmental impacts of increased wildfire intensity, frequency, and burned area. Previous studies have monitored the evolution of metal(loid) speciation at specific time/temperature endpoints, hampering the comprehension of the dynamic transformations of metal(loid)s during wildfires. Here, a novel approach, namely in situ time-resolved X-ray absorption near edge structure (TR-XANES) spectroscopy, was used for the first time to monitor the heat-induced redox transformations of arsenic (As) in two As-rich sediments with different characteristics. No arsenic transformations were observed in either sediment at room temperature or when the sediments were heated to a maximum temperature of 400 ºC. When heated to a maximum temperature of 600 or 700 ºC, As underwent complex, dynamic, and partially reversible redox transformations, the extent of which varied with sediment properties, initial As speciation, and heating temperature and duration. In the case of sediments initially containing As(V) and a low sulfide/sulfate ratio, the As(V) was reduced to As(III) over a short period of time, followed by immediate reoxidation of As(III) to As(V). The extent of the transient As reduction increased with increasing maximum temperature. The final As speciation at the end of the heating experiment was not substantially different from that at the beginning of the experiment. In contrast, in sediments that were best described by a more complex combination of reference compounds (As(V) coprecipitated with or adsorbed on iron oxide, arseniosiderite, arsenopyrite) with a high sulfide/sulfate ratio, As(V) phases were reduced to As(III) and As(III) was further reduced to As(-I) (at 700 ºC, only), followed by delayed reoxidation of As(-I) to As(III) (at 700 ºC, only) and reoxidation of As(III) to As(V). Approximately 9% and 26% of all As remained as As(III) ( i.e., was not reoxidized to As(V)) at the end of the experiment when heated to a maximum temperature of 600 and 700 ºC, respectively, which is higher than the relative abundance of As(III) in the native sediment. These fire-induced As transformations can potentially enhance As mobilization to surface water and groundwater, posing risks to environmental and human health.

California↗
USGS2025

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska↗
USGS2025

Characterizing Meteor Crater impact melts through geochemistry and textural analysis

The U.S. Geological Survey Astrogeology Science Center houses the Meteor Crater sample collection, an assemblage of over 2,500 meters of cuttings from 161 drill holes into Meteor Crater’s rim, flanks, and ejecta blanket. We have utilized this unique collection to study the composition and spatial distribution of impact-generated materials from within the ejecta blanket. Meteor Crater has historically been known to have generated only a relatively small amount of impact melt compared to other terrestrial craters of similar size. A detailed compositional and textural dataset of impact-derived melts from this impact can therefore be a useful asset in improving our understanding of crater formation, and in particular impact melt formation. We have characterized 42 impact-melt particles from Meteor Crater using a scanning electron microscope and an electron microprobe for textural and compositional analysis. We analyzed samples from six drill holes in the ejecta blanket, situated to the northwest, southeast, south, and southwest of the crater (ejecta northeast of the crater is devoid of impact melts). Impact melts were collected from drill cuttings at various depths within the ejecta blanket, ranging from a few centimeters below the surface down to ~6.5 meters. Backscattered electron (BSE) images were acquired for each analyzed impact-melt particle. To characterize the various textures and phases present in each impact melt, we also took many detailed BSE images. Our geochemical analyses include full spectral profiles using energy dispersive X-ray spectrometry and well-calibrated wavelength dispersive spectrometry for a number of phases, including minerals (olivine, pyroxene, and so on), pristine glass, and metallic inclusions. The full dataset is available in ScienceBase as a data release (Gullikson and others, 2024), accessible at https://doi.org/10.5066/P9OGAJ8P . Our goal for this Open-File Report is to provide a summary of this immense dataset, details on data collection, descriptions of the different phases observed within impact-melt particles (both geochemically and texturally), and observable trends.

Arizona↗
USGS2025

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

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

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

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

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

Journal of Analytical and Applied Pyrolysis↗
USGS2025

Reduction of red bed sedimentary rocks in connection with energy metal ore formation: A case study from the Sinbad seep, Mesa County, Colorado

The Paradox Basin’s Sinbad seep is a modern analog for ancient bleaching of red bed sediments by introduced alkaline, reducing brines. This bleaching, involving reductive alteration of former red beds, is essential ground preparation that enables the altered rocks to trap Cu, U, and V from later oxidized fluids, forming ore deposits. Study of Sinbad thus offers insights into these metallic mineralization processes in the Paradox and other sedimentary basins. The Sinbad seep occurs where shallow groundwater interacts with organic-rich petroleum source rocks, flows up a fault, and discharges into Salt Creek. Ratios of Na/Cl, Na/Br, and Cl/Br, plus high sulfate concentrations, indicate that the salinity of the seep water originated from dissolution of halite and gypsum during topographically driven flow of meteoric water across a diapir of the Pennsylvanian Paradox Formation, rather than from deep basinal brines. In contact with the organic-rich shales of the Paradox, these SO 4 -rich waters become reduced through bacterial reduction, producing H 2 S. The waters react with red beds of the Permian Cutler Formation, causing pervasive bleaching. The bleaching at Sinbad is characterized by iron-conservative reduction of ferric iron in diagenetic hematite and detrital ilmenite and magnetite to Fe sulfides. The ferrous iron is retained in these sulfides likely as sorbed Fe. Associated alteration includes precipitation of quartz and feldspar overgrowths, partial dissolution of detrital quartz and feldspar caused by pressure solution, formation of clay and authigenic rutile, and precipitation of carbonate and gypsum. The Fe sulfides rapidly degenerate to a mixture of jarosite and iron oxides on weathering at the surface. Among the major-oxide elements and most trace elements, there is no statistical difference between unbleached and bleached samples. The exceptions are uranium and sulfur, which are somewhat greater in bleached samples. Fission track radiography illuminates that uranium is preferentially concentrated in iron oxide cements, iron sorbed into detrital clasts, and finely disseminated iron oxide in illite cement in bleached samples. A leaching experiment suggests that uranium may be more easily available for mobilization from rocks that have undergone bleaching alteration, making them potential U sources. In addition, bleached rocks form effective traps for U, V, and some Cu mineralization.

Colorado↗