Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geoenergy Science and Engineering”

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

An early Holocene wet period in the southwestern United States

Multiple generations of spring-fed streams traversed ∼800 km 2 of the Las Vegas Valley in southern Nevada between ca. 10.9 ka and 8.5 ka, depositing an extensive tufa network. The scale of this network and diversity of tufa morphologies is novel in North America and offers an opportunity to obtain quantitative paleoclimate data for the region during the early Holocene. We determined isotopic compositions and estimated past temperatures using clumped isotope data from early Holocene tufa on the valley floor (698 m) as well as tufa forming today at higher elevation in the nearby Spring Mountains at Cold Creek Spring (1856 m). Modern and fossil tufa yielded comparably low δ 18 O values, implying that source waters for both were derived from high-elevation winter precipitation. Clumped isotope temperatures of modern tufa average 15.8 ± 2.5 °C, aligning with mean summer temperatures of the emergent spring water, and indicate equilibrium conditions of tufa formation. The early Holocene tufa yielded similar clumped isotope temperatures, averaging 15.2 ± 3.9 °C, meaning it precipitated at temperatures that occur at much higher elevations today. The Las Vegas tufa record, combined with nearby and temporally correlative paleospring and lacustrine records, suggest that cool/wet conditions prevailed throughout the Mojave Desert during the early Holocene. These records also demonstrate that spring ecosystems responded to millennial-scale hydroclimate variations that supersede climate change driven solely by insolation. The previously unrecognized pattern of ecosystem response to hydroclimate documented here may assist in understanding climate drivers for the early Holocene and provide critical information for the fate of groundwater-dependent ecosystems in the southwestern United States.

Nevada↗
USGS2026

Evaluation of submersible pressure transducers for streamflow monitoring in small streams

Compact streamgages requiring minimal infrastructure and equipped with submersible pressure transducers (PTs) are increasingly used to monitor small streams, yet disparate implementations obscure their accuracy under real-world conditions. This study isolated instrumentation-derived uncertainty in stage monitoring by co-locating various combinations of commercial vented and unvented PTs with seven U.S. Geological Survey (USGS) reference streamgages on small streams. Multi-year PT stage records from compact streamgages, collected and corrected following USGS protocols, were compared to concurrent reference observations. Vented PTs demonstrated an average measurement uncertainty of ±0.005 m and mean absolute percent error (MAPE) of ±0.2%. Unvented PTs exhibited higher uncertainty, averaging ±0.009 m and ±0.4% MAPE. Although both sensor types had stage errors up to ±0.3 m, 95% of vented PT and unvented PT observations were within 0.01 and 0.02 m of reference stage, respectively. Analysis of additional unvented PTs revealed stage errors of up to ±1.2% when using barometric sensors within 15 km of the in-water sensor. Propagation of stage error to discharge using reference rating models resulted in cumulative discharge MAPEs of ±4.5% for vented and ±5.5% for unvented PTs. These findings highlight PTs as practical alternatives to reference instrumentation when deployed with standardized procedures, potentially expanding access to reliable streamflow data.

conterminous United States↗
USGS2025

Fossil diatoms in Arctic marine surface sediments

Diatoms are one of the main constituents of marine phytoplankton in the Arctic, and thanks to their siliceous skeletons, diatom fossils are relatively well preserved in sediments. Due to their species-specific sensitivity to different ocean conditions, their abundance and assemblages in sediments are routinely used by paleoceanographers to reconstruct the state of the surface ocean in the past (Koç 2007). By using statistical methods, we can obtain valuable knowledge about their ecological preferences (Oksman et al. 2019) and generate quantitative reconstructions of various parameters, such as sea-surface temperature and sea-ice concentration through time (Krawczyk et al. 2021; Sha et al. 2014). The Marine Arctic Diatoms (MARDI) working group (WG) ( pastglobalchanges.org/mardi ) aims to advance knowledge on marine-diatom ecology and diatom-based reconstructions by compiling and harmonizing data from surface-sediment samples across the Arctic.

Pages Magazine↗
USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗
USGS2024

Integration of rupture directivity models for the US National Seismic Hazard Model

Several rupture directivity models (DMs) have been developed in recent years to describe the near-source spatial variations in ground motion amplitudes related to propagation of rupture along the fault. We recently organized an effort towards incorporating these directivity effects into the USGS National Seismic Hazard Model (NSHM), by first evaluating the community's work and potential methods to implement directivity adjustments into probabilistic seismic hazard analysis (PSHA). Guided by this evaluation and comparison among the considered DMs, we selected an approach that can be readily implemented into the USGS hazard software, that provides an azimuthally varying adjustment to the median ground motion and its aleatory variability. This method allows assessment of the impact on hazard levels and provides a platform to test the DM amplification predictions using a generalized coordinate system, necessary for consistent calculation of source-to-site distance terms for complex ruptures. We give examples of the directivity-related impact on hazard, progressing from a simple, hypothetical rupture, to more complex fault systems, composed of multiple rupture segments and sources. The directivity adjustments were constrained to strike-slip faulting, where DMs have good agreement. We find that rupture directivity adjustments using a simple median and aleatory adjustment approach can impact hazard both from a site perspective and on a regional scale, increasing shaking off the end of the fault trace up to 30--40\% and potentially reducing it for sites along strike. Statewide hazard maps of California show that the change in shaking along major faults can be a factor to consider for assessing long-period (>ls) near-source effects within the USGS NSHM going forward, reaching up to 10--20\%. Finally, we suggest consideration of minimum parameter ranges and baseline requirements as future DMs are developed to minimize single approach adaptations, to enable more consistent application within both ground motion and hazard studies.

Earthquake Spectra↗
USGS2026

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗
USGS2025

The role of alkali bicarbonate-sulfate brines in the genesis of carbonatite REE resources at the Bear Lodge Alkaline Complex, Wyoming

Rare-earth element (REE) resources in the Bear Lodge Alkaline Complex, Wyoming, are hosted in a variably leached carbonatite dike swarm spatially related to bodies of diatreme breccia. This study examines fluid inclusions in carbonatite dikes, peripheral fluorite breccias, and smoky quartz veins to reconstruct the physiochemical conditions of REE mineralization. Results reveal a multi-stage evolution of carbonatitic fluids: 1) high-temperature (330–432°C), CO 2 -rich magmatic fluids, 2) REE-enriched alkali bicarbonate-sulfate brines, and 3) low-temperature meteoric (117–182°C) or diluted magmatic fluids. Multiphase inclusions (stage-2) contain burbankite, nahcolite, and alkali sulfate daughter crystals, linking the brines to early burbankite [(Na,Ca) 3 (Sr,Ba,REE) 3 (CO 3 ) 5 ] mineralization in carbonatite. REE mobility has been linked to the stability of aqueous complexes with ligands such as Cl-, SO 4 2- , and HCO 3 - . However, light REEs were not initially transported significant distances, rather they crystallized as burbankite from brines or brine-melts associated with the carbonatite magma. Crystallization of burbankite removed Na from the system, leaving residual fluids with high K/Na ratios and promoting potassic fenitization throughout the complex. Fractionation of light and heavy REEs and involvement of low temperature (<300°C), high K/Na bicarbonate brines may have facilitated transport and formation of peripheral occurrences that are more enriched in heavy REEs at Bear Lodge.

Wyoming↗
USGS2026

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin↗
USGS2026

Evaluating harvest liberalization strategies on population dynamics of southern latitude temperate-breeding Canada geese

The recovery of the Canada goose ( Branta canadensis ) is one of North America's greatest conservation success stories. Today, continental abundance of temperate-breeding Canada geese (those breeding in southern Canada and lower 48 states in the United States) greatly exceeds historical levels. As a result of increased abundance, human–goose conflicts have also increased, ranging from private and agricultural property damage to human health concerns. Managers have primarily attempted to lower Canada goose populations using hunter harvest via liberalized hunting regulations (increased bag limits and total hunting days). To evaluate the effectiveness of harvest strategies for temperate-breeding Canada geese, managers need a better understanding of how liberalized hunting regulations affect population dynamics. We estimated survival and harvest probabilities, abundance, and recovery distribution of temperate-breeding Canada geese banded in Arkansas, USA, during 2005–2020. We found that adult harvest probabilities declined overall during the study and adult survival probabilities increased. Annual abundances of juvenile geese (≤1 year old) declined during the study, whereas there was no detectable trend in abundance of the adult population over time. Most Canada geese breeding in Arkansas were shot in Arkansas, indicating potential for population dynamics to be influenced largely by state-specific harvest regulations. However, our results suggest that harvest liberalization as a management tool may have limited capacity to further influence population dynamics of Canada geese. Declines in productivity and recruitment likely had a greater influence on the similar observed declines in juvenile and adult abundances, respectively, during our study. More research is needed to better understand the ecological mechanism affecting the population dynamics of Canada geese in both urban and rural environments.

Arkansas↗
USGS2025

40Ar/39Ar geochronology supporting mineral resources research at USGS Denver

The 40 Ar/ 39 Ar geochronology method is used to date potassium-bearing rocks and minerals, based on the decay of 40 K to 40 Ar, which provides important temporal constraints for geological events. The USGS Denver Argon Geochronology Laboratory dates samples from a variety of projects, mainly in the USGS Mineral Resource Program and the National Cooperative Geologic Mapping Program, facilitating in-depth research into the timing of geological events and processes. Recent applications of this method include geochronological studies in the Yellow Pine district, revealing mineralization ages that range from 51 to 70 Ma, and investigations into the Bear Lodge alkaline complex, providing insights into the timing of REE carbonatite formation. Additionally, studies at Alunite Ridge highlight the method's potential for understanding concealed mineral deposits. Ongoing projects continue to leverage 40 Ar/ 39 Ar data to enhance understanding of geologic frameworks and mineral resources, underpinning its importance in modern geochronology and mineral resource assessment.

Conference Paper↗
USGS2025

Waning greenhouse gas emissions from U.S. Federal lease coal production by the mid-21st century

This study presents estimates of future years (2024–2051) United States Federal lease coal production and the resulting greenhouse gas (GHG) emissions from the combustion, transport, and mining of that fuel. Results from the coal production estimate indicate a decline in production from Federal leases; with known production of 240 million short tons (mtn) in 2023 and a projected decline to 34.0 mtn by 2051, which represents a reduction to 14.2% of the 2023 value. In parallel with this projection, total GHG emissions are estimated to decrease from 402.2 million metric tons of carbon dioxide equivalent (MMT CO 2 eq.) in 2024 to 55.0 MMT CO 2 eq. in 2051, a decline to 13.7% of 2024 emissions estimates. The reductions in coal production and emissions are mainly the result of planned coal combustion power plant closures, with major projected closures in 2037 and 2048. However, GHG emissions estimates for future years can be uncertain as they rely heavily on coal production estimates from operators' public business plans and other publicly available resources. Forward looking plans of this type are subject to significant changes if economic and political factors deviate from current information. Results suggest that average GHG emissions over the time series breakout to 95% end point combustion, 3.7% transportation combustion emissions, and 1.3% fugitive emissions, although there is uncertainty associated with these figures. Uncertainty stemming from production projections, sector distributions, and emissions factors on the future emissions estimates increases with time, ranging from −28% to +48% within the 2024–2051 timeframe.

Alabama, Colorado, Montana, North Dakota, Utah, Wy↗
USGS2026

ShakeAlert Earthquake Early Warning System performance during the Mw 7.0 offshore Cape Mendocino earthquake

The 5 December 2024 M w 7.0 Offshore Cape Mendocino earthquake was a challenging test of the U.S. West Coast ShakeAlert earthquake early warning system due to its offshore epicenter and limited near‐source station coverage. We analyzed real‐time performance of all components of the ShakeAlert system, including the seismic algorithms (earthquake point‐source integrated code [EPIC] and Finite‐fault rupture Detector [FinDer]), the geodetic algorithm (Geodetic First Approximation of Size and Time–peak ground displacement [GFAST‐PGD]), and network telemetry during the event. EPIC created the first solution for this earthquake 15 s after origin time with an initial magnitude estimate of M 5.6 and location error of 10 km from the Advanced National Seismic System epicenter. An early spurious trigger from station CE.89101 fortuitously maintained location accuracy and, correspondingly, magnitude accuracy. FinDer contributed its first solution at 18 s with a location estimate closer to the seismic network and produced two distinct rupture geometries, leading to minor fluctuations in estimated intensity contours. GFAST‐PGD did not meet alerting thresholds but otherwise performed as expected. Network latencies were <2 s for most stations, supporting the rapid detection of this earthquake by the system. Roughly five million alerts were delivered to cell phone devices in California and Oregon during this event. This was also the first instance of a school district‐wide ShakeAlert‐powered system being activated. Comparisons to recorded seismograms demonstrate that the maximum warning times before potentially damaging shaking (intensity 6+) were in the range of 5–55 s. Although the ShakeAlert system provided accurate solutions and useful alert delivery, this earthquake raised awareness of potential issues within the system, including the need for improved offshore location estimates, a combination of solutions from ShakeAlert servers, and handling of spurious triggers.

California↗
USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗
USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS2025

Framework for assessing intranational risk of critical mineral sectors: Case study of India

Recent geopolitical actions regarding the supply of critical minerals have highlighted both the threats posed by countries dominating a critical mineral supply chain, and the potential mitigations served by alternative supplier countries. However, no systematic quantification exists for the intranational factors that determine such threats or mitigations. This study outlines a methodological framework for quantifying the intranational risk profile of critical mineral sectors and compares them with the sector's international potential for leverage. The four intranational factors of concentration, synchronization, private sector involvement, and regulatory obstacles were analyzed across multiple scales in order to determine an intranational risk score that quantifies the resilience or fragility of the critical mineral sector in question. Network analysis was also separately conducted to identify the most important entities operating across all critical mineral sectors. Finally, this study assessed the critical mineral sector's global share of reserves and production to determine an international risk score that quantifies the country's critical mineral sector's leverage potential. A case study of India showed that two critical mineral sectors with high leverage potential, chromite and barite, were respectively intranationally resilient and fragile. Of the remaining sectors with minimal leverage, five sectors such as aluminum were intranationally resilient and may serve as viable mitigation options for future critical mineral supply. Three sectors such as titanium were intranationally fragile, but recent regulatory developments may both decrease intranational risk and increase international risk, potentially making them major critical mineral sectors in the future.

Resources Policy↗
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↗
USGS2026

Evidence for a biological origin of uranium-rich carbon masses within the Ediacaran Salt Range Formation of Pakistan

Thucholites are unique organic structures found in igneous and sedimentary rocks composed of a U-C-rich interior enclosed by an organic outer shell. Their formation and occurrence have perplexed scientists for over 100 years. Typically, thucholites are sparse in sedimentary rocks but where found in abundance, they may be the result of rapid paleoecological disruptions, e.g., volcanic ashfall. Here, we evaluated thucholites from the Ediacaran Salt Range Formation of the Indus Basin in Pakistan using field emission scanning electron microscopy (FESEM) and nanoscale secondary ion mass spectroscopy (NanoSIMS) to propose a two-stage mechanism for thucholite development. NanoSIMS results suggest organic matter in thucholite cores formed by biological fractionation (represented by the presence of lighter 13 C/ 12 C and 34 S/ 32 S) while the outer organic mantle formed via radiolysis-induced polymerization. FESEM elemental analysis confirms compositional differences between the two thucholite components (core and mantle), further implying their contrasting origins. Dimensional comparison of thucholite cores and their U-bearing mineral morphologies to ancient and modern U-biomineralized microbes suggests that the thucholite cores formed from metal-biological interactions. The presence of volcanogenic biotite, alkali feldspar, and clay spherules (interpreted as devitrified volcanic glass) with thucholite suggests that the thucholite cores in the Salt Range Formation are biological responses to rapid paleoenvironmental change from volcanic eruption which preserved the fossilized morphologies of Proterozoic microorganisms.

Salt Range↗