USGS Science⌕ Search

SEARCH · USGS Science

Results for “Critical Reviews in Environmental 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,679 records · Page 71Linked to original sources
USGS2024

Revision of ModelMuse to support the use of PEST software with MODFLOW and SUTRA models

Executive Summary ModelMuse is a graphical user interface for several groundwater modeling programs. ModelMuse was updated to generate the input files for the parameter estimation software suite PEST. The software is used with MODFLOW or SUTRA models to run PEST-based parameter estimation and display the updated model inputs after parameter estimation. The PEST input files can also be used with the PEST++ version 5 software suite. Parameter estimation typically requires defining the parameters being adjusted during calibration and observations for assessing calibration quality. After a parameter is defined in ModelMuse, it can be applied to all or part of a model dataset. Pilot points—a parameterization device that facilitates higher levels of parameterization—can be used to assign spatially variable distributions of model inputs. Parameters can be applied to temporally varying features, such as boundary conditions, by either applying them to all the values in a series in one step or by applying separate parameters to individual members of a series. ModelMuse allows the definition of many observation types from various model output files. For MODFLOW 6 and SUTRA models, new options were added to ModelMuse to allow it to display the changed input after parameter estimation is complete. For MODFLOW–2005 and MODFLOW–NWT models, ModelMuse can import an entire model for visualization. An example illustrates the use of PEST with a MODFLOW 6 model in ModelMuse.

Techniques and Methods↗
USGS2025

What 25+ years of "Did You Feel It" intensities tell us about shaking in California

“When will the Big One happen?” is a question that people often have for earthquake scientists. But while waiting for the “Big One” to occur, people will usually experience frightening or damaging shaking from multiple relatively smaller‐magnitude earthquakes. Given this context, it raises the question: “Where does most of the damage come from?” Could smaller, yet more frequent, earthquakes account for the majority of reported impactful shaking? To explore this question, we consider reports of earthquake damage and felt shaking experiences from a catalog of community‐collected intensity values from the U.S. Geological Survey’s “Did You Feel It?” system. Comparing these intensities to expectations from a ground‐motion model, we find that earthquakes of magnitudes smaller than expected are responsible for most reported intensities of community decimal intensities (CDI) 4.5 and above (moderate and higher shaking intensity levels). (Here “expected value” is meant in its mathematical sense of the mean or equivalently median shaking intensity.) We also present a regional analysis of observed earthquake shaking for specific areas, detailing the maximum intensity experienced within a predetermined area. We identify several instances of M < 4.5 events that generated maximum intensities of CDI > 5 in regions in California surrounding the cities of Eureka, San Francisco, Los Angeles, and San Diego. Our results motivate the need to include smaller‐magnitude earthquakes in communications about earthquake hazard and risk reduction.

California↗
USGS2024

Projected sea-level rise and high tide flooding at Big Cypress National Preserve, Florida

Introduction National parks and preserves in the South Atlantic-Gulf Region contain valuable coastal habitats such as tidal wetlands and mangrove forests, as well as irreplaceable historic buildings and archeological sites located in low-lying areas. These natural and cultural resources are vulnerable to accelerated sea-level rise and escalating high tide flooding events. Through a Natural Resources Preservation Program-funded project during 2021–23, the U.S. Geological Survey, in collaboration with the National Park Service, estimated the probability of inundation at Big Cypress National Preserve, Florida, and several other parks under various sea-level rise scenarios and contemporary high tide flooding thresholds. The maps produced for this effort can be used to assess potential habitat change and explore how infrastructure and cultural resources within the park may be exposed to future flooding-related hazards.

Florida↗
USGS2024

Developing international standards and guidelines for disseminating and validating simulated ground motions

We are leading an effort to develop international standards and guidelines for curating, disseminating, and validating simulated ground-motion data. This effort is organized as a working group within the Consortium of Organizations for Strong Motion Observation Systems (COSMOS). In 2022 we held online workshops on 7 and 8 June and 20 October. The first workshop focused on curating and disseminating simulated ground-motion data. The second workshop focused on validating simulated ground-motions for engineering applications. About 100 people participated in each of the workshops with strong representation from North America, Europe, and western Asia. In the coming year, we intend to form a technical committee to draft the international guidelines and standards while continuing to engage stakeholders from across the globe through online and in-person meetings. The key points from the workshops include: (1) Numerous groups are generating simulated earthquake ground motions and making them openly available; however there is very little coordination among groups to provide consistent interfaces for searching and retrieving data; (2) Standardizing interfaces for metadata and data access should consider agile approaches that can adapt to changing capabilities and user needs while building upon existing efforts; (3) A distributed architecture is preferred by most participants to allow institutions to host and manage their own data while broadcasting their holdings to a combined catalog; (4) Validation of ground-motion simulations applies to the entire workflow for simulating earthquake ground-motions, including the rupture model, seismic velocity model, and seismic wave propagation software; (5) Metrics to evaluate the validation are application dependent; more research is needed to tie structural response characteristics to ground-motion characteristics; and (6) Validation results should provide a clear, transparent, and quantitative assessment of the simulated ground motions.

Conference Paper↗
USGS2024

Geoelectric evidence for a wide spatial footprint of active extension in central Colorado

Three-dimensional magnetotelluric (MT) imaging in central Colorado reveals a set of north-striking high-conductivity tracks at lower-crustal (50–20 km) depths, with conductive finger-like structures rising off these tracks into the middle crust (20–5 km depth). We interpret these features to represent saline aqueous fluids and partial melt that are products of active extensional tectonomagmatism. These conductors are distributed over a wider region than the narrow corridor along which Rio Grande rift structures are traditionally mapped at the surface, and they consequently demarcate regions of the lower crust where accommodation of bulk extensional strain has concentrated conductive phases. Our observations reveal limitations in existing models of Rio Grande rift activity and may reflect unrecognized spatiotemporal variations in rift system evolution globally.

Colorado↗
USGS2025

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

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

Massachusetts↗
USGS2026

Development of liquefaction-targeted design basis in U.S. seismic provisions

Seismic design criteria for new buildings in the United States have historically focused on life safety performance objectives through minimizing the potential for structural collapse. Development of design criteria to meet this performance objective has evolved over time, leading to the current, risk-targeted maximum considered earthquake (MCE R ) design basis. Corresponding MCE R seismic design loads incorporate the full range of probabilistic ground motion hazard information available from U.S. Geological Survey National Seismic Hazard Models, as well as a representation of the uncertainty in structural collapse capacity. Despite these considerable advances in structural design criteria, guidelines for geotechnical hazards such as liquefaction-induced damage, which continue to be based on uniform ground-shaking hazard (i.e., MCE ground motions), are not well-connected with risk-based, collapse prevention performance objectives. This study, undertaken as part of Building Seismic Safety Council efforts to update recommended seismic provisions through the National Earthquake Hazards Reduction Program, presents probabilistic, liquefaction-targeted design objectives and corresponding ground motion parameters, with a focus on improving consistency in first-level screening criteria for liquefaction hazard assessment. This study identifies a potential maximum acceptable annualized probability of liquefaction triggering ( P L,max ) of 2.5% in 50 years, above which liquefaction consequences and potential mitigation measures would be considered. The P L,max -based design basis is calibrated to effective probabilistic liquefaction hazard levels obtained in practice using current MCE peak ground accelerations, and could be implemented via a new liquefaction-targeted maximum considered earthquake (MCE L ) peak ground acceleration ( PGA L ). The proposed PGA L is derived from probabilistic liquefaction hazard curves and includes uncertainties inherent to liquefaction hazard modeling. The potential design effects of the proposed PGA L are demonstrated via comparison with current MCE PGA levels, effects on first-level liquefaction screening analyses, and improvements in the consistency of P L,max levels across the United States, and for a wide range of site conditions.

conterminous United States↗
USGS2025

Water use across the conterminous United States, water years 2010–20

Withdrawals of water for human use are fundamental to the evaluation of the Nation’s water availability. This chapter provides an analysis of public supply, crop irrigation, and thermoelectric power water use for the conterminous United States (CONUS) during water years 2010–20. These three categories account for about 90 percent of water withdrawals in the Nation. The values presented here are based on modeling approaches that estimate water use at temporal (monthly) and spatial scales (12-digit hydrologic unit code—small watersheds sized 50–100 square kilometers) compatible for integration into a broader national assessment of water availability. Models also provide an understanding of factors that influence water use. An estimated 244,817 million gallons per day (Mgal/d; 28,677 million cubic meters per month [Mm 3 /mo]) were withdrawn on average within the CONUS during water years 2010–20 from fresh water and saline water for crop irrigation, public supply, and thermoelectric power, with shares of 43, 14.5, and 42.5 percent for each of these categories, respectively. In the same period, estimated withdrawals and consumptive use (1) for public supply were 35,400 and 4,219 Mgal/d (4,081 and 486 Mm 3 /mo), respectively; (2) for crop irrigation were 105,497 and 75,698 Mgal/d (12,147 and 8,716 Mm 3 /mo), respectively; and (3) for thermoelectric power from fresh water were 82,656 and 2,904 Mgal/d (9,952 and 345 Mm 3 /mo), respectively. Withdrawals for these categories of water use are highly spatially variable, with western States dominated by crop irrigation and eastern States dominated by thermoelectric-power water use. Public supply accounts for the largest percentage of water use in several heavily populated northeastern States. Reliance on groundwater compared to surface water depends on the availability of water sources and the type of water use. For public supply, withdrawals from groundwater are greater than withdrawals from surface water in the Western aggregated hydrologic regions, whereas the balance shifts to more surface water for the rest of the CONUS. In all aggregated hydrologic regions, the predominant source of water for crop irrigation is groundwater. Most thermoelectric power facilities in the eastern half of the CONUS use surface water from freshwater and saline sources; most thermoelectric power facilities in the western half of the CONUS use groundwater.

conterminous United States↗
USGS2025

Groundwater budget for the surficial aquifer surrounding Lake Nokomis, Minneapolis, Minnesota

During prolonged periods of above-average precipitation, rising groundwater levels have the potential to cause damage to and interfere with underground infrastructure and building foundations. To understand the relations between precipitation and groundwater in the vicinity of Lake Nokomis, the U.S. Geological Survey, in collaboration with the University of Minnesota, quantified five components of the groundwater budget: groundwater recharge, change in surficial aquifer storage, surficial aquifer groundwater discharge to Lake Nokomis, groundwater evapotranspiration, and groundwater discharge to underlying bedrock aquifers. Field data, geologic records, and empirical calculation methods were used to quantify groundwater budget components for April 2023 through April 2024. Lake water budget data indicate that Lake Nokomis is a flowthrough system during periods with no outflow through the weir, with groundwater inputs equal to outputs. Roughly 40 percent of precipitation that fell in the study area was added to the surficial aquifer as recharge. Uncertainty in the vertical hydraulic conductivity resulted in wide-ranging estimates (spanning three orders of magnitude) of water discharging from the surficial aquifer to the underlying bedrock aquifer. Drought conditions persisted for the duration of this study and were not representative of the conditions that motivated this study. This study is a start towards understanding relations between precipitation, Lake Nokomis levels, and groundwater levels that could affect local underground infrastructure.

Minnesota↗
USGS2026

Investigating the influence of climate and volcanic surface aging on fluvial erosion: A case study of Réunion Island, Indian Ocean

Precipitation is one of the dominant drivers of landscape erosion and evolution; however, the effects of typical rainfall compared with less frequent, high-magnitude precipitation events on erosion remain unclear. Volcanic islands are ideal locations to study such phenomena due to their simple geometries, nontectonic construction, and strong spatiotemporal rainfall gradients. However, spatial variation in surface age, created during their construction, often complicates their degradation histories by introducing temporal changes in erosion rates as drainage networks develop. Réunion Island (western Indian Ocean) presents a clear example of this, with an east–west gradient in both surface age and mean annual precipitation, as well as infrequent cyclones that alter the background rainfall pattern. In this study, we analyze the effects of surface age, average rainfall, and rainfall variability on basin development and fluvial erosion across the island. We calculate basin-averaged values of basin morphology, age, precipitation, river discharge, eroded volumes, and erosion rates, and use these to analyze the dominant drivers of landscape evolution through a series of correlation analyses. Our results indicate a temporal dependence on the influence of precipitation, with young surfaces being dominantly eroded by high-rainfall events and older surfaces eroded by mean annual rainfall patterns. Furthermore, we show that drainage development of shield volcanoes follows similar trends to other volcano types, and suggest that surface permeability and groundwater structure are important controls on runoff-driven erosion on shield volcanoes. These results add new components to the question of how precipitation impacts erosion.

Réunion Island↗
USGS2025

Road salt collection and redistribution at an urban rain garden on sandy soil, Gary, Indiana

Rain gardens installed as green infrastructure to divert storm runoff from entering combined sewers also collect dissolved constituents and particulates. An urban rain garden in northwestern Indiana, USA, was continuously monitored from November 2019 to May 2021 to evaluate the fate of dissolved constituents entering the rain garden in runoff. Physical and chemical properties of soils in the rain garden were also monitored, along with underlying groundwater. Linear regression models relating specific conductance to chloride concentration indicated that the 0.0371-ha (3998 square feet) rain garden collected approximately 1490 kg (3285 pounds) of road salt from the surrounding 0.2228 ha (24,500 square feet) of impervious surfaces. Soils and groundwater were seasonally affected by road salt application but carryover from year to year was not indicated. Rain garden soil permeability (5.20 × 10 −5 to 9.72 × 10 −5 m/s) remained unchanged during the study period and soil organic carbon generally increased under native vegetation. The results suggest that a rain garden built on sandy soil can divert substantial quantities of runoff and dissolved constituents from combined sewers; however, chloride is transported to sub-infrastructure groundwater that eventually discharges to adjacent waterways with concentrations lower than those observed in runoff.

Indiana↗
USGS2025

False positives in the identification of dynamic earthquake triggering

Dynamic earthquake triggering is commonly identified through the temporal correlation between increased seismicity rates and global earthquakes that are possible triggering events. However, correlation does not imply causation. False positives may occur when unrelated seismicity rate changes coincidently occur at around the time of candidate triggers. We investigate the expected false positive rate in Southern California with global M ≥ 6 earthquakes as candidate triggers. We compute the false positive rate by applying the statistical tests used by DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 to synthetic earthquake catalogs with no real dynamic triggering. We find a false positive rate of ∼3.5%–8.5% when realistic earthquake clustering is present, consistent with the 95% confidence typically used in seismology. However, when this false positive rate is applied to the tens of thousands of spatial-temporal windows in Southern California tested in DeSalvio and Fan (2023), https://doi.org/10.1029/2023jb026487 , thousands of false positives are expected. The expected false positive occurrence is large enough to explain the observed apparent triggering following 70% of large global earthquakes (DeSalvio & Fan, 2023, https://doi.org/10.1029/2023jb026487 ), without requiring any true dynamic triggering. Aside from the known triggering from the nearby El Mayor-Cucapah, Mexico, earthquake, the spatial and temporal characteristics of the reported triggering are indistinguishable from random false positives. This implies that best practice for dynamic triggering studies that depend on temporal correlation is to estimate the false positive rate and investigate whether the observed apparent triggering is distinguishable from the correlations that may occur by chance.

JGR Solid Earth↗
USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗
USGS2024

New methodology for assessing underground natural gas storage resources – Example from Michigan Basin, United States

Energy consumption in the United States (U.S.) and across the world is shifting away from traditional fossil fuels like coal and oil, and towards natural gas and renewable sources, including hydrogen. Because gas demand is typically greatest during cold seasons and renewable sources sometimes produce variable supplies, it is important to store energy for use when demand exceeds supply. Whereas batteries and tanks typically store energy above ground, geologic (underground) storage may be able to retain much greater quantities of energy over much longer durations (e.g., Matos et al. 2019; Buursink et al. 2023). Consequently, the U.S. Geological Survey (USGS) is developing a methodology to locate new geologic energy storage sites with an initial focus on depleted hydrocarbon reservoirs followed by estimates of associated pore space or gas storage capacity. Depleted hydrocarbon fields in Michigan Basin in the U.S. already host large quantities of injected natural gas (about 681 billion cubic feet of seasonal gas capacity in 42 gas storage facilities; U.S. Energy Information Administration 2024). Therefore, we used this basin as an example to assess potential new storage resources alongside existing storage capacity and carbon dioxide sequestration resources (e.g., Katz and Coats 1968; Haagsma et al. 2020). In Michigan Basin gas storage is taking place in reservoirs within several key hydrocarbon producing formations with different lithologies, including the informal Stray sandstone of the Michigan Formation and the pinnacle reefs of the Niagara Group (Gautier et al. 1995; Swezey et al. 2015). The assessment methodology introduced here consists of two steps relying on both well information and reservoir data.

Michigan↗
USGS2026

Alternative approaches to dryland reclamation enhance vegetation cover and soil stability at a former uranium mine

Introduction Conventional reclamation methods to recover dryland ecosystems after mineral extraction often have low success. Alternative reclamation techniques may help overcome the many challenges to establishing persistent native vegetation and stable soils. Objectives Compare conventional reclamation approaches that focus primarily on soil stabilization to alternative approaches that may improve native ecosystem restoration at the Kanab North Mine in northern Arizona. Methods In a fully crossed experiment, conventional seed mix (native and non-native species) and seeding methods (drill-seeding) were compared to alternative practices, including a native-only seed mix, connectivity modifiers (ConMods) with hand seeding, and a biocrust inoculation. Restoration metrics like plant establishment, species composition, soil surface cover and stability, biocrust establishment, and near-surface wind-driven sediment were monitored over the first 5 years post-reclamation. Results While all restoration metrics improved with time, compared to traditional methods, ConMods increased total, graminoid, and native plant species cover, and the native-only seed mix increased woody, native, and perennial species cover. Biocrust inoculation increased late-successional biocrust development, but biocrusts only persisted over 5 years in combination with ConMods. Soil stability increased in all treatments, but both ConMods and biocrust inoculation increased stability compared to traditional practices. Conclusion Overall, alternative approaches to reclamation tested here (ConMods, native seed, and biocrust inoculation) were more effective than traditional practices at increasing total plant cover, facilitating biocrust establishment, and promoting soil stability after uranium mine reclamation.

Arizona↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2026

Localization of spatiotemporally heterogeneous subsurface flows using autoencoder-based deep learning framework for time-lapse self-potential tomography

Self-potential (SP) monitoring has emerged as a valuable method for characterizing subsurface hydrogeological features and processes due to its sensitivity to fluid-induced electrokinetic effects. Despite advancements in SP inversion, challenges remain in imaging groundwater dynamics from SP activities due to complex hydrological settings and transient noise. In this study, a deep learning autoencoder (AE)-based framework is proposed for the spatiotemporal localization of subsurface fluid movement from time-lapse SP tomography. Temporal segments of time-lapse numerical inversions were first derived from long-term SP monitoring conducted from a floodplain site in Oak Ridge, Tennessee, known for active hyporheic exchange. Subsequently, AE models based on vision transformer (ViT), convolutional long short-term memory (ConvLSTM), convolutional neural network, and temporal convolutional network were individually trained and compared on the SP tomography segments for reconstruction performance. Finally, the reconstruction error over time serves as an anomaly score to identify moments of active SP variation, whereas spatial distributions of errors within these moments are analyzed to image and localize regions associated with anomalous subsurface fluid movement. The results demonstrate that ConvLSTM- and ViT-AE are most capable for the localization task with contrasting error distributions and consistent delineation of anomalies. Applying the method to both SP arrays parallel and perpendicular to the stream produced consistent anomaly zones near a fault or karst feature, validating the robustness and generalization of the approach. These results demonstrate the potential of the proposed framework as a scalable and interpretable tool for spatiotemporal analysis of subsurface flow dynamics in complex hydrogeological systems.

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