USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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.

At least 649 records · Page 36Linked to original sources
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS2025

Statewide cumulative human health risk assessment of inorganics-contaminated groundwater wells, Montana, USA

Across the United States, rural residents rely on unregulated and generally unmonitored private wells for drinking water, which may pose serious health risks due to unrecognized contaminants. We assessed the nature, degree, and spatial distribution of cumulative health risks from inorganic contaminants in groundwater. Our analysis included nearly 84,000 data points from 6500+ wells, across 51 of Montana's 98 watersheds, using a public groundwater database. We compared a drinking water screening level cumulative risk assessment (CRA) for inorganics based on the U.S. Environmental Protection Agency (EPA) protective health thresholds (Maximum Contaminant Level Goals, Health Advisories [MCLG-HAs]) to a CRA based on EPA public supply enforceable Maximum Contaminant Levels (MCLs). Based on median concentrations of 19 inorganics (antimony, arsenic, barium, beryllium, boron, cadmium, chromium, copper, fluoride, manganese, molybdenum, nickel, nitrate, lead, selenium, strontium, thallium, uranium, zinc), 75% of watersheds had MCLG-HA-based cumulative risk values > 1.0; arsenic and uranium contributed the most risk, followed by strontium, fluoride, manganese and boron. Hence, this screening level (Tier I) CRA indicated widespread potential for unrecognized human health risk to private well users from inorganic contaminants considering both carcinogenic and non-carcinogenic risks. Sensitivity analysis showed that benchmarks applied (MCLG-HAs versus MCLs) exerted the largest control on results. Our findings identify priority regions for Tier 2 risk assessments to elucidate local sources and distributions of geogenic versus anthropomorphic contaminants. Our study is the first statewide assessment of cumulative health risk from groundwater that we are aware of, and results support increased statewide drinking water education and testing to reduce human health risks from contaminated private well water.

Montana↗
USGS2025

Evaluation of daily stream temperature predictions (1979-2021) across the contiguous United States using a spatiotemporal aware machine learning algorithm

Stream temperature controls a variety of physical and biological processes that affect ecosystems, human health, and economic activities. We used 42 years (1979–2021) of data to predict daily summary statistics of stream temperature across >50,000 stream reaches in the contiguous United States using a recurrent graph convolution network. We comprehensively documented the performance – both across all reaches and by stream type (e.g., reservoir or groundwater influence) – as a baseline for future improvement. The model showed reach-level RMSE of <2 °C with 90 % prediction intervals that contain 90.7 % of observations. We also assessed how the model captured variability in ecologically relevant metrics (e.g., R 2 for annual 7-day maximum = 0.76; R 2 for days exceeding 25 °C = 0.75). This model does not outperform state-of-the-art machine learning efforts (e.g., RMSE ≤1.5 °C) due to a limited input set but does provide the most spatially complete modeling to date to support water availability assessments.

contiguous United States↗
USGS2025

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

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

Earth and Planetary Science Letters↗
USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗
USGS2025

Triple-oxygen isotopic evidence of prolonged direct bioleaching of pyrite with O2

Sulfate is often touted as containing atmospheric oxygen whose isotopic signature can constrain redox, environmental conditions, and biological activity. Yet, the amount and isotopic fractionation associated with air-O 2 incorporation during sulfate formation is still debated, making its verification difficult. In this study, we identify a distinct, microbially dominated environment with the potential to preserve maximum signals of air-O 2 in sulfate. We report triple-oxygen isotope data for sulfate produced from pyrite oxidation in microbial and abiotic experiments, and from natural dissolved sulfate from the Rio Tinto, Spain, an acid mine drainage site. The oxygen isotope systematics of sulfate in these environments define a unique kinetic isotope effect associated with initial stage pyrite oxidation by Acidithiobacillus ferrooxidans that preserves >80 % oxygen from air-O 2 in sulfate. Unlike experiments, which evolve toward water-oxygen dominated sulfate on short time scales, Rio Tinto, Spain hosts a microbe rich environment with distinct geochemistry that maintains high O 2 -oxygen in sulfate. Therefore, in addition to containing isotopic records from water and air, sulfates can also contain a biosignature that is promising for understanding conditions on Mars and early Earth.

Earth and Planetary Science Letters↗
USGS2025

Ultramafic float rocks at Jezero crater (Mars): Excavation of lower crustal rocks or mantle peridotites by impact cratering?

Based on observation and data from meteorites and in situ scientific missions, experiments as well as models, the Martian mantle is assumed to share some compositional and mineralogical affinity with the terrestrial mantle. However, there might be subtle differences like the Martian mantle being more ferroan. Yet, we do not have any direct analysis of a Martian mantle rock to confirm this assumption. NASA’s Perseverance rover found olivine-rich boulder-sized float rocks on the upper Jezero fan (Mars). These boulders have an ultramafic composition and their mineralogy is dominantly composed of Fo 73±3 olivine with high-Mg orthopyroxene, Cr-rich Ti-Fe oxides and minor plagioclase and high-Ca pyroxene. Microtextural and petrological analysis reveals that these minerals crystallized at equilibrium. In addition, these boulders are different from all the bedrocks analyzed by Perseverance along its traverse which are crustal igneous rocks and sediments. Comparing our data to Martian meteorites and available Mars bulk silicate models (BSM), we discuss that these boulders could represent primitive melts and/or lower crustal material, and we specifically hypothesize that they could be mantle peridotites. We propose that these putative mantle rocks could have been excavated by the succession of impacts from the shallow mantle or lower crust in the Isidis region where Jezero crater is located. These olivine-rich boulders could thereby constitute the first direct analysis of a Martian mantle rock.

Earth and Planetary Science Letters↗
USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗
USGS2025

Sand provenance boundary in the Mu Us Sandy Land of northern China

Desert dunes are often assumed to have uniform mineral compositions due to extensive mixing during lateral transport, which complicates provenance studies. The Mu Us Sandy Land in north-central China, near the East Asian summer monsoon precipitation boundary, experiences a wetter climate than most deserts. Climate wetting as a result of a warming climate, and the ‘Sand Control Project’ implemented by the Chinese government to decrease the lateral movement of sand dunes in this area provide an opportunity to study surface processes of sand production and transport. Previous studies using zircon U Pb geochronology and heavy mineral composition indicate distinct sand sources for the Mu Us Sandy Land: local basement-derived middle Yellow River sediments and recycled dried-up lacustrine sediments for the eastern part, and northeastern Tibetan Plateau-derived upper Yellow River sediments for the western part. However, zircons and heavy minerals only represent trace amounts of the bulk mineralogy within a sand dune, so broader provenance analysis targeting common minerals is essential. We focus on the area near the proposed provenance difference boundary between the western and eastern Mu Us Sandy Land, using comprehensive sampling and provenance techniques to confirm distinct provenance characteristics and delineate the provenance boundary. Our findings reveal that sand from most of the Mu Us Sandy Land originated from the erosion of local basement by the middle Yellow River and recycled local dried-up lacustrine sediments, whereas the southwestern corner and the neighboring western-central Chinese Loess Plateau received sediments from the distal northeastern Tibetan Plateau.

Mu Us Sandy Land↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2025

Favorability mapping for hydrothermal power resource assessments of the Great Basin, USA

The U.S. Geological Survey (USGS) is updating the 2008 assessment of conventional hydrothermal resources for the Great Basin in the western United States. As part of this work, the workflow for hydrothermal resource favorability maps is being modified to integrate modern data-driven machine learning (ML) methods. Improvements include: [1] using new and refined evidence layers (features); [2] using an order of magnitude more training sites (labeled examples); [3] utilizing simple but non-linear supervised ML algorithms; [4] representing positive training sites (wells with measured heat flow) with their ordinal value proportional to the magnitude of convective upflow (i.e., low, high, or very high convective signals instead of past strategies using positive-negative labels); [5] supplementing training sites with additional sites with low convective signals to represent diverse under-sampled areas where hydrothermal systems are unlikely to exist; [6] comparing with competing approaches; and [7] utilizing Monte Carlo cross-validation to estimate and evaluate prediction uncertainty. For the new favorability map, over half of the power-producing systems (i.e., 15 of 28) are predicted in the 99th percentile of most favorable locations (i.e., the highest 1 % of favorability, corresponding to 1 % of the map area), exceeding the performance of past models that have explicitly used power plants as training sites. Previous favorability maps predicted approximately half of the power-producing hydrothermal systems above the 80th percentile (i.e., 20 % of the map area). For the new favorability map, 93 % of power-producing systems (i.e., 26 of 28) are above the 80th percentile. The power-producing systems for which the new model does not perform well are either comparatively small, low-temperature systems or systems also not predicted well by prior modeling approaches, suggesting that these few systems are unusual when compared with most power-producing systems. Focusing research on these known, seemingly different systems may yield new insights and subsequent discovery of new prospects.

California, Idaho, Nevada, Oregon, Utah↗
USGS2025

Simulating human behavior under earthquake early warning

Earthquakes are a rapid-onset hazard where advance planning and learning plays a key role in mitigating injuries and death to individuals. Recent advances in earthquake detection have resulted in the development of earthquake early warning (EEW) systems. These systems can provide advance warning to predetermined geographic regions that an earthquake is in progress, which may result in individuals receiving warning seconds before significant shaking is felt at their location. This additional time could allow individuals to take more effective protective actions during the immediate disaster. To demonstrate this, we created an agent-based simulation of a basic apartment that allowed us to randomly and repeatedly simulate an individual receiving and responding to an EEW message. The results of our preliminary simulation show that, in our study environment, earthquake early warning alerts have the potential to allow for sufficient time for individuals to take protective actions.

Washington↗
USGS2025

Pediment formation and subsequent erosion in Gale crater: Clues to the climate history of Mars

Evidence of paleo-rivers, fans, deltas, lakes, and channel networks across Mars has prompted much debate about what climate conditions would permit the formation of these surface water derived features. Pediments, gently sloping erosional surfaces of low relief developed in bedrock, have also been identified on Mars. On Earth, these erosional landforms, often thought to be created by overland flow and shallow channelized runoff, are typically capped by fluvial sediments, and thus in exceptionally arid regions, pediments are interpreted to record past wet periods. Here we document the Greenheugh pediment in Gale crater, exploiting the observational capability of the Curiosity rover. The absence of a fluvial cap suggests that the pediment was likely cut by wind erosion, not fluvial processes. The pediment was then buried by an aeolian deposit (Stimson sandstone) that mantled the lower footslopes of Aeolis Mons (informally known as Mt. Sharp). This burial terminated active wind erosion, preserving the pediment surface (as an angular unconformity). Groundwater was present prior-to, during, and shortly after Stimson deposition, perhaps contributing to lithification and certainly to early diagenesis. Post lithification, wind erosion first cut canyons in the northern most footslopes (north of Vera Rubin ridge). Unlithified gravels were deposited in these canyons, likely due to runoff from Mt. Sharp. Boulder-rich fluvial and debris flow deposits built a > 70 m thick sequence (Gediz Vallis ridge) on the southern Greenheugh pediment. Continued wind erosion left elevated patches of gravel on the northern footslopes, and exposure age dating shows that erosion essentially ceased before 1 Ga (but possibly much earlier). Erosion to the south led to emergence of Vera Rubin ridge, retreat of the Greenheugh pediment, and the formation of Glen Torridon valley. Hence, this footslope environment of Mt. Sharp records climate-driven periods of wind erosion, aeolian deposition (and groundwater activity), surface runoff and sediment deposition, followed by further significant wind erosion that declined to present very slow rates. This likely occurred during the late Hesperian and possibly into the Amazonian.

Icarus↗
USGS2025

Estimated average annualized losses from potential building damage and fatalities due to earthquake-generated tsunamis in the United States

Earthquake-generated tsunamis represent substantial economic threats to states and territories in the United States (U.S.), but we are unaware of any effort to quantify potential impacts at the national level. This gap is partially due to the lack of nationally consistent data on tsunamigenic sources and associated return periods. This study addresses this issue and provides estimates of average annualized losses (AAL) for potential residential fatalities and capital stock losses associated with building damage (i.e., structural, non-structural, contents, and inventory damage) in the U.S. by curating tsunami-hazard information based on deterministic scenarios and probabilistic approaches, calculating potential losses, and estimating return periods where necessary. This assessment was done for the U.S. West Coast, Alaska, Hawaii, U.S. Pacific Territories, and U.S. Atlantic Territories. We estimate that earthquake-generated tsunamis that could affect these states and territories collectively represent $1 billion in potential AAL with 79 % of losses due to residential fatalities and 21 % of losses due to capital stock losses from building damage. We identify AAL variations based on county and county equivalents, states and territories, geographic regions, return periods, and departure-delay assumptions for evacuating residents. Results include high AAL values for potential fatalities in Puerto Rico and the U.S. Pacific Northwest region, high AAL values for potential building-related damage in Hawaii and California, and high building- and population-loss ratios for county equivalents in Alaska and U.S. territories.

International Journal of Disaster Risk Reduction↗
USGS2025

Evaluation of solid bitumen created from marine oil shale bituminite under hydrous and anhydrous pyrolysis conditions

To test the influence of environmental conditions on aromaticity of solid bitumen generated during petroleum generation, four organic-rich (26–36 wt% total organic carbon) oil shale samples collected from the Neoproterozoic–Lower Cambrian restricted marine Salt Range Formation in the upper Indus Basin, Pakistan, were pyrolyzed under hydrous and anhydrous conditions. Experiments used closed system batch reactors at subcritical water temperatures between 300 and 370°C for 72 h. Thermal conversion of bituminite in the Salt Range oil shales created a newly formed solid bitumen, similar to previous observations from experiments on the Eocene lacustrine Green River Formation Mahogany zone oil shale. Newly formed solid bitumen in the Salt Range Formation oil shales generally has higher reflectance (R o ) in hydrous residues compared to anhydrous experiments subjected to the same time-temperature conditions, also similar to prior observations. This finding supports the hypothesis that radical disproportionation is favored in the presence of hydrogen contributed by water, promoting aromatization in the solid bitumen residue with concomitant increased expulsion of saturated hydrocarbons. Indigenous solid bitumen (and vitrinite in a comparison sample) also shows higher reflectance values in hydrous versus anhydrous residues, indicating that additional aromatization in the presence of exogenous hydrogen occurs both in newly formed organic matter and in organic matter that is present throughout the experiment. Despite similarities in their bulk rock geochemical screening parameters, R o evolution shows different trajectories amongst the four Salt Range oil shales, suggesting as-yet undetermined differences in kinetic properties which are probably related to differences in a priori chemical composition. These results have implications for the use of solid bitumen reflectance (BR o ) as a thermal proxy, suggesting BR o values and appearance could vary as a function of the concentration of water. Variation in water concentration may be present at the reservoir or formation scale, but may also be present at a much finer scale in tight oil shales where permeability is several orders of magnitude lower than conventional reservoirs. Therefore, local variations in the presence of water potentially could explain substantial variation in BR o values and appearance in closely spaced source rock samples and even within an individual microscope field.

Journal of Analytical and Applied Pyrolysis↗
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

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2025

Using hydrogeologic context and water budgets to evaluate the potential for groundwater contributions to contaminants in Lake Superior

This study presents a synthesis of the hydrogeology in the U.S. Lake Superior watershed and the contribution of groundwater to the water budget of the U.S. Lake Superior basin. Much of the shoreline of Lake Superior in Minnesota and Michigan is composed of hydrogeologic units contributing very little direct groundwater discharge to the lake. Groundwater in watersheds adjacent to the lake typically flows in short, local flow systems characterized by thin glacial sediments with active groundwater flow in fractured bedrock within the top 60–90 m below land surface. The complex groundwater system in Wisconsin’s Bayfield Peninsula has the largest groundwater reservoir near the lake, characterized by thick sand and gravel glacial deposits and underlying sandstone aquifer. Although these thick sandy glacial deposits are not in direct contact with the lake at the shoreline, groundwater discharge may still be significant via subsurface exposures beyond the shoreline or flow through the underlying sandstone aquifer. Overall, most groundwater in the watershed is contributed as indirect base flow in streams around the lakeshore and comprises about 60 percent of the upland water budget. Direct groundwater flow to the shoreline contributes 2 to 9 percent of the inflow. Identifying possible contamination sources through direct sampling of groundwater would be an inefficient way to detect problems if sources are unknown, particularly for some chemicals of concern such as PFAS, pesticides, PCBs, chloride, and nutrients. Evaluating the chemical characteristics of contamination is also important to consider in evaluating how groundwater may contribute to pollution in Lake Superior.

Lake Superior↗