USGS Science⌕ Search

SEARCH · USGS Science

Results for “Resources, Conservation, and Recycling”

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

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗
USGS2026

Monitoring changes in Landsat thermal features in urban and non-urban interfaces from 1986 to 2023 in two international urban centers: Implications for climate and global issues

Rapid urbanization is reshaping thermal environments worldwide, with the strongest impacts occurring at the interface between urban and non-urban areas. Impervious surfaces, as key indicators of urban expansion, are critical for monitoring urban growth and assessing surface urban heat island (SUHI) effects. Land use and land cover change (LULCC) provides an essential link between urban dynamics and their environmental and societal consequences. Here, we integrated the U.S. Geological Survey (USGS) Climate Global Issues (CGI) Land Cover Product with Landsat thermal time-series to investigate SUHI evolution in two contrasting metropolitan regions: Wuhan, China, and Brasília, Brazil. Using data spanning 1986–2023, we analyzed the relationships between land cover, Landsat-based land surface temperature (LST), and SUHI intensity, and identified persistent thermal hotspots. Results demonstrate that the land cover data utilized increases the accuracy of impervious surface mapping along urban–rural gradients. Average SUHI intensities were 3.4 °C in Wuhan and 3.3 °C in Brasília, with statistically significant warming trends of 0.04 °C/year and 0.01 °C/year, respectively. Maximum temperature proved to be a robust indicator of SUHI intensification, capturing long-term upward trends. Our findings highlight the important role of urban land cover dynamics in shaping temporal SUHI variability and hotspot emergence. This prototype framework demonstrates the scientific and policy value of combining long-term land cover monitoring information with satellite thermal monitoring to quantify and track SUHI at city scale, supporting sustainable urban planning and climate adaptation strategies.

Remote Sensing↗
USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗
USGS2025

Paleomagnetic correlation of surface and subsurface basalt flows in the central and southwestern part of the Idaho National Laboratory, Idaho

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, used paleomagnetic data from 22 coreholes to construct 3 fence diagrams of subsurface basalt flows in the southern part of the Idaho National Laboratory. These diagrams provide comprehensive descriptions of the horizontal and vertical distribution of basalt flows and sediment layers beneath the surface, aiding geological studies and contributing valuable data to numerical models of groundwater flow and contaminant transport. The correlations established though these diagrams include spatial correlations between basalt flows found in multiple coreholes. Correlations were identified by matching average paleomagnetic inclinations and confirming or denying these correlations using petrology, geochemistry and radiometric ages. The fence diagrams aid in identifying potential locations of subsurface vents, volcanic vents that have been buried by more recent volcanic activity, associated to subsurface basalt flows. By tracing the subsurface flows and analyzing where the greatest thickness occurs, the locations of buried vents can be inferred. Some subsurface flows exhibit correlations across several coreholes and may indicate yet unidentified surface or buried vents, thereby enhancing our understanding of the volcanic history and subsurface geology of the region.

Idaho↗
USGS2024

Wildfire, extreme precipitation and debris flows, oh my! Channel response to compounding disturbances in a mountain stream in the Upper Colorado Basin, USA

Compounding changes in climate and human activities stand to increase sediment input to rivers in many landscapes, including via discrete perturbations such as post-fire debris flows. Because sediment supply is a dominant control on river morphology, understanding mountain river responses to sediment regime perturbations is critical to predicting and addressing downstream effects to infrastructure, water security and aquatic habitat. A growing body of literature explores the causes, likelihood, size and composition of post-fire debris flows, but the channel response to these disturbances remains poorly studied. This study used repeat field surveys, time-lapse photographs and pre- and post-disturbance remote sensing datasets to document and analyse space- and time-varying channel response to post-fire debris flows along a steep mountain stream in the Upper Colorado River Basin, USA. Specifically, we evaluated channel morphology and bed composition changes, correlations between channel changes and valley and channel attributes, and the relative importance of spring snowmelt versus summer monsoon events. Several cross-sectional channel change types were observed from lidar a month after post-fire debris-flow events, including channelized and braided incision into deposits, incision into the pre-fire channel bed, bank erosion and no change. Channel changes were most correlated with pre-fire channel width, valley width and unit stream power, and these relationships could be tested in other burned locations to evaluate their transferability. Repeat channel surveys before and after snowmelt indicate rapid recovery and channel narrowing following major sediment disturbances, although sediment deposits remained in the channel margins. Together, these results highlight the importance of field and remote sensing-based channel surveys to improve understanding of, and potential to predict, mountain channel response to compounding climate disturbances.

Earth Surface Processes and Landforms↗
USGS2025

Timescales of surface faulting preservation in low-strain intraplate regions from landscape evolution modeling and the geomorphic and historical record

Large surface-rupturing intraplate earthquakes in stable continental regions (SCRs) are uncommon globally and have recurrence intervals of thousands to hundreds of thousands of years based on the paleoseismic and geomorphic record, challenging accurate active fault identification in these regions. To constrain the timescales of preservation for scarps created by surface ruptures from dip-slip earthquakes, we use a two-dimensional scarp diffusion model for typical intraplate settings and explore which parameters influence fault scarp preservation. These parameters include the coseismic vertical surface offset, the recurrence interval of similar magnitude earthquakes, diffusivity (as a proxy for mean annual precipitation rate), and the erodibility of the surficial material. We constrain parameter ranges from a compilation of historical surface ruptures in intraplate settings in a variety of climates, including the Central and Eastern United States, Australia, Europe, Central Asia (Mongolia, China), India, and West Africa. The timescales of scarp preservation from landscape evolution modeling agree well with observations of scarp preservation in low-strain SCR and intraplate tectonic settings, with some notable exceptions for Australian scarps. We find that the erodibility of the surficial material and earthquake recurrence interval have a stronger effect on the timescales of scarp preservation than diffusivity or coseismic vertical surface offset. Our model results may aid in identifying and characterizing subtle, slow-moving active faults in low-strain SCR and intraplate tectonic settings for different tectonic, geomorphic, and climatic characteristics. Accurate fault locations and characterization from the landscape record has implications for both probabilistic seismic and fault displacement hazard analyses.

Journal of Geophysical Research Solid Earth↗
USGS2025

Expression of corticoid-regulatory genes in the gills of Atlantic salmon (Salmo salar) parr and smolt and during salinity acclimation

In teleost fishes, cortisol is the major corticoid and has both glucocorticoid and mineralocorticoid actions. However, how fish tissues discriminate between these distinct corticosteroid actions is unclear. In mammals, the major factors responsible for intracellular corticosteroid regulation are glucocorticoid receptors ( gr s) and the mineralocorticoid receptor ( mr ), but their role in osmoregulation of fish is unclear. 11β-hydroxysteroid dehydrogenases ( hsd11b s) control the levels of intracellular corticosteroids by converting from bioactive forms to inert forms. To investigate how Atlantic salmon ( Salmo salar ) respond to cortisol in different physiological or environmental conditions, we performed comparisons of parr and smolt, and osmotic challenge experiments to examine the physiological responses and gill transcript levels of genes underlying cortisol-signalling, including gr1 , gr2 , mr , hsd11b2 and hsd11b3 . Because cortisol may interact with growth hormone and prolactin during salinity changes, transcript levels encoding growth hormone receptors ( ghr1 , ghr2 ) and the prolactin receptor ( prlr ) were also examined. Hsd11b2 transcript levels in seawater-acclimated fish were consistently lower compared to fish acclimated to fresh water. After transfer to seawater, prlr transcript levels in fish significantly decreased and transcript levels of ghr1 , ghr2 and hsd11b3 showed no change or were slightly higher than those of freshwater control groups. Gr1 , gr2 and mr transcript levels were slightly but consistently higher in fish acclimated to fresh water relative to seawater. Our results indicate that changes in corticosteroid receptor and hsd11b2 transcript levels in the gills may be important mechanisms that regulate corticoid signals to achieve ion homeostasis in Atlantic salmon.

Journal of Fish Biology↗
USGS2025

Physicochemical properties and bioreactivity of sub-10 µm geogenic particles: Comparison of volcanic ash and desert dust

Exposure to ambient particulate matter (PM) with an aerodynamic diameter of <10 μm (PM 10 ) is a well-established health hazard. There is increasing evidence that geogenic (Earth-derived) particles can induce adverse biological effects upon inhalation, though there is high variability in particle bioreactivity that is associated with particle source and physicochemical properties. In this study, we investigated physicochemical properties and biological reactivity of volcanic ash from the April 2021 eruption of La Soufrière volcano, St. Vincent, and two desert dust samples: a standardized test dust from Arizona and an aeolian Gobi Desert dust sampled in China. We determined particle size, morphology, mineralogy, surface texture and chemistry in sub-10 μm material to investigate associations between particle physicochemical properties and observed bioreactivity. We assessed cellular responses (cytotoxic and pro-inflammatory effects) to acute particle exposures (24 hr) in monocultures at the air-liquid interface using two types of cells of the human airways: BEAS-2B bronchial epithelial cells and A549 alveolar type II epithelial cells. In acellular assays, we also assessed particle oxidative potential and the presence of microorganisms. The results showed that volcanic ash and desert dust exhibit intrinsically different particle morphology, surface textures and chemistry, and variable mineralogical content. We found that Gobi Desert dust is more bioreactive than freshly erupted volcanic ash and Arizona test dust, which is possibly linked to the presence of microorganisms (bacteria) and/or nanoscale elongated silicate minerals (potentially clay such as illite or vermiculite) on particle surfaces.

GeoHealth↗
USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗
USGS2025

Characterization of the hydrogeologic framework, groundwater-flow system, geochemistry, and aquifer hydraulic properties of the shallow groundwater system in the Wilcox and Lorraine process areas of the Wilcox Oil Company Superfund site near Bristow, Oklahoma, 2022

The Wilcox Oil Company Superfund site (hereinafter referred to as “the site”) was formerly an oil refinery northeast of Bristow in Creek County, Oklahoma. Historical refinery operations contaminated the soil, surface water, streambed sediments, alluvium, and groundwater with refined and stored products at the site. The Wilcox and Lorraine process areas are where the highest concentrations of volatile organic compounds, semivolatile organic compounds, polycyclic aromatic hydrocarbons, and trace elements (including metals) (collectively hereinafter referred to as “contaminants”) were measured in a local shallow perched groundwater system within the alluvium (hereinafter referred to as the “alluvial aquifer”) at the site during previous site assessments. In order to understand the potential migration of contaminants through the soil and groundwater in these areas, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, investigated aquifer characteristics of the alluvial aquifer in the Wilcox and Lorraine process areas of the site to (1) document hydraulic conductivity and other aquifer characteristics of the alluvial aquifer that govern contaminant fate and transport, (2) describe the geospatial extent and concentration of the contaminants in the alluvial aquifer in the Wilcox and Lorraine process areas, and (3) describe the geochemical controls pertaining to oxidation and reduction governing the fate and transport and the degradation potential of contaminants in the groundwater. Various data were compiled and collected to evaluate the aquifer characteristics at the site including the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer. A total of 20 new (2022) groundwater monitoring wells were installed at the site to collect data used to supplement groundwater-level altitude and groundwater-quality data collected from older, existing groundwater monitoring wells and piezometers. Data compiled and collected for the study were used to evaluate the characteristics of the alluvial aquifer at the site. These aquifer characteristics are defined by the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer.

Oklahoma↗
USGS2025

Hydrologic variability and groundwater age of springs in eastern Oregon and northern Nevada, USA

The ecological importance of springs in semiarid regions is far greater than their small size and sparse distribution, yet little is known about the hydrologic functioning of these systems. During 2016–22, 261 springs were visited in the volcanic terrane of eastern Oregon and northern Nevada. When conditions were suitable, measurements of discharge, water temperature, and specific conductance were made, and samples for the analysis of carbon-14, tritium, and water stable isotopes (WSI) were collected. A subset of 60 springs was revisited during different seasons in the same year and during the dry season in multiple years to evaluate variability in discharge, chemistry, and groundwater age. Specific conductance and WSI varied considerably among springs across the study area but were unexpectedly stable across seasons and years at individual springs. Seasonal and interannual variability in spring discharge was related to the residence time of the discharging groundwater. Springs discharging older groundwater (10 3 –10 4 years) had significantly less variability in their discharge compared to springs discharging younger groundwater (10 0 –10 1 years). Variability among springs discharging younger groundwater included cessation of late-summer discharge at 18 % of the repeat-visit springs. A logistic regression model predicted the age of discharging spring water with 89 % accuracy using only the spring latitude, longitude, elevation, and δ 2 H value. This study framework provides a simple, inexpensive, and robust method to provisionally assess the hydrologic behavior of springs having little or no prior information in understudied, semiarid regions across the globe.

Nevada, Oregon↗
USGS2024

Bioconcentration of per- and polyfluoroalkyl substances and precursors in fathead minnow tissues environmentally exposed to aqueous film-forming foam-contaminated waters

Exposure to per- and polyfluoroalkyl substances (PFAS) has been associated with toxicity in wildlife and negative health effects in humans. Decades of fire training activity at Joint Base Cape Cod (MA, USA) incorporated the use of aqueous film-forming foam (AFFF), which resulted in long-term PFAS contamination of sediments, groundwater, and hydrologically connected surface waters. To explore the bioconcentration potential of PFAS in complex environmental mixtures, a mobile laboratory was established to evaluate the bioconcentration of PFAS from AFFF-impacted groundwater by flow-through design. Fathead minnows ( n = 24) were exposed to PFAS in groundwater over a 21-day period and tissue-specific PFAS burdens in liver, kidney, and gonad were derived at three different time points. The ∑PFAS concentrations in groundwater increased from approximately 10,000 ng/L at day 1 to 36,000 ng/L at day 21. The relative abundance of PFAS in liver, kidney, and gonad shifted temporally from majority perfluoroalkyl sulfonamides (FASAs) to perfluoroalkyl sulfonates (PFSAs). By day 21, mean ∑PFAS concentrations in tissues displayed a predominance in the order of liver > kidney > gonad. Generally, bioconcentration factors (BCFs) for FASAs, perfluoroalkyl carboxylates (PFCAs), and fluorotelomer sulfonates (FTS) increased with degree of fluorinated carbon chain length, but this was not evident for PFSAs. Perfluorooctane sulfonamide (FOSA) displayed the highest mean BCF (8700 L/kg) in day 21 kidney. Suspect screening results revealed the presence of several perfluoroalkyl sulfinate and FASA compounds present in groundwater and in liver for which pseudo-bioconcentration factors are also reported. The bioconcentration observed for precursor compounds and PFSA derivatives detected suggests alternative pathways for terminal PFAS exposure in aquatic wildlife and humans.

Environmental Toxicology and Chemistry↗
USGS1976

Guidelines for collection and field analysis of ground-water samples for selected unstable constituents

The unstable nature of many chemical and physical constituents in ground water requires special collection procedures and field analysis immediately after collection. This report describes the techniques and equipment commonly used in the collection and field analysis of samples for pH, temperature, carbonate, bicarbonate, specific conductance, Eh, and dissolved oxygen.

Open-File Report↗
USGS2024

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within northern Medina County, Texas

During 2023–24, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised a previous publication of the geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers that was completed during 2018–20 within northern Medina County, Texas. The purpose of this report is to present the updated geologic framework and hydrostratigraphy of the rocks containing the Edwards and Trinity aquifers in northern Medina County from field observations of the surficial expressions of the rocks. The report includes a detailed 1:24,000-scale hydrostratigraphic map with names and descriptions of the geologic and hydrostratigraphic units (HSUs) in the study area. This study includes updates to the geology of the Kainer Formation (or its stratigraphic equivalent, the Fort Terrett Formation) with the addition of a burrowed unit between the basal nodular and dolomitic members. The hydrostratigraphy was also updated with the addition of HSU IIA for the upper part of the Devils River Limestone and the Seco Pass HSU for the burrowed member of the Kainer (or Fort Terrett) Formation. The Cretaceous age rocks (listed in ascending order) in the study area are part of the Trinity Group, Edwards Group and stratigraphically equivalent Devils River Limestone, Washita Group, Eagle Ford Group, Austin Group, and Taylor Group, with isolated areas where Late Cretaceous age igneous rocks have intruded. The groups and formations are composed primarily of relatively thick layers of clays, shales, and limestone. The igneous rocks are coarse-grained ultramafic in composition. Hydrostratigraphically, the rocks exposed in the study area (listed in descending order from land surface) are igneous, the upper confining unit to the Edwards aquifer, the Edwards aquifer, the upper zone of the Trinity aquifer, and the upper part of the middle zone of the Trinity aquifer. The karstic carbonate Edwards and Trinity aquifers developed because of their original depositional history, primary and secondary porosity, diagenesis, fracturing, and faulting. These factors have resulted in development of modified porosity, permeability, and transmissivity within and between the aquifers.

Texas↗
USGS2025

Sources and risk factors for nitrate, pathogens, and fecal contamination of private wells in rural southwestern Wisconsin, USA

Household well water can be degraded by contaminants from the land's surface, but private well owners lack means to protect the source water from neighboring disturbances. Rural residents of southwestern Wisconsin, USA, rely on private well water, and the combination of land use and fractured carbonate bedrock makes groundwater vulnerable to contamination. To identify the extent, sources, and risk factors of private well contamination, randomly selected wells sampled during two-day periods in fall (n = 301) and spring (n = 529) were analyzed for nitrate and indicator bacteria, and a subset (n = 138) was sampled across four seasonal events for analysis of pathogens and microbial source tracking markers by quantitative polymerase chain reaction. Risk factors representing land use, hydrology, geology, and well construction were analyzed for associations with contamination in multivariable models. The importance of risk factors varied by contaminant, illustrating the multifaceted nature of rural groundwater quality. Nitrate contamination was associated with agricultural land use, and wells with casings that extended below a shale aquitard accessed less contaminated water than those drawing water from above it. Human fecal microbes were detected in 64 wells (46%), and rainfall was the key risk factor for contamination, indicating that wastewater from septic systems was available to contaminate wells when transport conditions were favorable. Manure microbes from cattle/ruminants and pigs were detected in 33 and 13 wells, respectively, and concentrations increased with the hectarage of cultivated land near wells. Pathogen genes for viruses, bacteria, and protozoa were detected in 66 wells (48%), including more detections of zoonotic than human-specific pathogens, and human Bacteroides , an indicator of wastewater, was an equivocal predictor of pathogen presence in private wells. Characterizing important elements of the setting, like geology, and identifying sources and risk factors for contaminants can inform landscape-level policies to protect groundwater quality.

Wisconsin↗
USGS2024

SUTRA— A code for simulation of saturated-unsaturated, variable-density groundwater flow with solute or energy transport—Documentation of the version 4.0 enhancements—Freeze-thaw capability, saturation and relative-permeability relations, spatially varying properties, and enhanced budget and velocity outputs

Version 4.0 of the Saturated-Unsaturated Transport (SUTRA) software code provides the capability to simulate the freezing and thawing of groundwater during energy transport simulations under saturated and unsaturated conditions. In addition to the types of hydrogeologic processes that SUTRA has been able to simulate in the past, this version can be used to study the effects of the freeze-thaw process on the flow and energy dynamics of hydrogeologic systems. The freeze-thaw simulation capability accounts for the latent heat of fusion and allows thermal property values to vary with changing total-water saturation, liquid-water saturation, and ice saturation. It allows the effective permeability of the porous medium to change as a result of freezing and thawing. This version also provides several user-selectable relations for the dependence of total-water saturation on fluid pressure, the dependence of liquid-water saturation on temperature during freezing and thawing, and the dependence of relative permeability on liquid saturation, as well as three user-selectable formulae for defining the bulk thermal conductivity of a mixture of solid grains, liquid water, ice, and air. For unsaturated simulations without freezing, the selectable total-water saturation relations eliminate the need for the user to program these and their associated relative-permeability functions, as had been required in previous SUTRA versions. Optional nonlinear dependence of fluid density on temperature, which covers the range from supercooled (about −50 degrees Celsius) to superheated (about 400 degrees Celsius), is also provided. Additionally, this version makes it possible to spatially vary parameters that, in previous versions of SUTRA, were required to be spatially uniform: solid-matrix properties, adsorption parameters, and parameters for production of solute mass or energy. Spatial variation is also allowed for the newly included freeze-thaw process parameters. Additional enhancements provide (1) output of water-mass and energy budgets that include values of all component terms in the governing balance equations, and (2) output of Darcy velocities (fluid fluxes), in addition to the velocity output provided by previous SUTRA versions. These enhanced outputs allow fuller interpretation of simulation results, especially for freeze-thaw phenomena. The set of processes simulated by this version of SUTRA are useful for studying a wide range of hydrogeologic system types, conditions, and questions. For cryohydrogeologic simulations, however, this version of the code is limited in that (1) it does not simulate thermomechanical effects of freeze-thaw, (2) pressure changes due to water density change during freezing are neglected, (3) ice saturation cannot exceed the initial porosity of the simulated medium, and (4) cryosuction, the migration of liquid water toward freezing fronts, is neglected. Furthermore, this version does not account for air flow or for water vaporization and sublimation under unsaturated conditions.

Techniques and Methods↗
USGS2025

Geochemistry of quartz and fluid inclusions associated with gold, tungsten, and antimony in the Stibnite-Yellow Pine district, Idaho

Deposits in the historic Stibnite-Yellow Pine district in east-central Idaho were mined for Ag, Au, Hg, Sb, and W. Ore is hosted in breccia and fault zones along NE-striking faults that transect the Atlanta lobe of the Late Cretaceous Idaho batholith and metamorphosed Neoproterozoic to lower Paleozoic rocks. Geochronological studies have delineated five superimposed and zoned hydrothermal events in the mining district which range in age from Late Cretaceous to the Eocene. They formed the following consecutive assemblages: (1) vein and disseminated pre-ore milky quartz (± calcite), muscovite, and brassy pyrite with minor molybdenite and base metal sulfides; (2) vein and disseminated quartz, calcite, sericite, auriferous-pyrite, and arsenopyrite; (3) potassium feldspar, quartz, dolomite, and calcite with vein and disseminated scheelite ± pyrite and arsenopyrite; (4) Au-Ag veins that contain adularia; and (5) breccia, vein, and disseminated quartz, calcite, and stibnite with minor Ag sulfosalts and distal zones of cinnabar and realgar at Fern and Hermes mines. In this study, the pressure, temperature, and composition of ore-forming fluids were measured in Au-, W-, and Sb-bearing samples collected from drill core and surface exposures. Using optical petrography, scanning electron microscopy cathodoluminescence (SEM-CL) and electron microprobe cathodoluminescence (EMPA-CL) imaging, as well as electron microprobe (EMP) and laser ablation-inductively coupled plasma-mass spectrometry trace element analyses, four major overprinting generations of hydrothermal quartz are defined. Fluid inclusion assemblages (FIAs) in each generation of quartz were characterized by microthermometry and laser Raman spectroscopy. Quartz and fluid inclusion paragenesis are as follows: (Qi) altered and overprinted Qi with no coeval FIAs that formed before or during event 1; (Q1) early veins as well as disseminated quartz with dull gray SEM-CL and H 2 O-NaCl-CO 2 FIAs that formed at minimum temperatures and pressures of 281° to 358°C and 1346 to 1789 bar during event 2 (Au); (Q2) vein and disseminated quartz, often euhedral, with H 2 O-NaCl ± CO 2 ± CH 4 FIAs that formed at minimum temperatures and pressures of 183° to 213°C and 72 to 602 bar during event 3 (W); and (Q3) euhedral and disseminated quartz with H 2 O-NaCl FIAs that formed at minimum temperature and pressures of 156° to 161°C and 5 to 6 bar during event 5 (Sb). Stibnite hosts H 2 O-NaCl FIAs trapped at minimum temperatures and pressures of 150° to 189°C and 5 to 12 bar that are similar to those of Q3-hosted FIAs. The quartz paragenetic history in the district follows a progressive temporal path from Tienriched quartz deposited by a relatively deep and higher temperature aqueous-carbonic fluid to a Li, Al, and Sb enriched quartz precipitated from a lower temperature aqueous fluid. Gases in fluid inclusion extracts from ore and gangue minerals were analyzed by mass spectrometry. Noble gas isotope analysis of fluid inclusion extracts revealed three signatures. Pre-ore Qi has R/R A (the isotopic ratio of the sample 3 He/ 4 He to atmospheric 3 He/ 4 He (1.384 × 10 −6 )) values that reflect the composition of the crustal melt-derived Atlanta lobe of the Idaho batholith. Quartz (Q1) and auriferous arsenian pyrite of hydrothermal event 2 have 20 Ne/ 4 He and R/R A values that plot along the air-crust mixing line, and these gases were likely introduced by metamorphic fluids derived from metasedimentary rocks. Stibnite with associated quartz (Q3 ± Q2) have 20 Ne/ 4 He and R/R A values that extend from the air-crust mixing line to the air-mantle mixing line, showing a shift from a crustal signature to an overprinting mantle signature. Late calcite has 20 Ne/ 4 He and R/R A values that extend from air toward the mantle end member, which suggests that magmatic volatiles condensed into meteoric groundwater late in the paragenesis. A mantle volatile source is also supported by high 20 Ne/ 22 Ne ratios. The Stibnite-Yellow Pine district is similar to other Au-W-Sb districts around the world in ore fluid chemistries, with late-stage stibnite being a product of cooling. The fluid pressure-temperature-composition (P-T-X) estimates and the crustal source of He during hydrothermal event 2 gold mineralization are characteristic of orogenic Au deposits, and the observed inclusion reequilibration textures indicate that trapped fluid inclusions were deformed during transitions between the brittle-ductile environment. Scheelite was deposited during event 3 with Q2 quartz and is overgrown by Q3 quartz and stibnite. Our interpretation of the quartz texture, chemistry, and fluid inclusion P-T-X data suggest that Au, W, and Sb mineralization occurred as hydrothermal fluids evolved from H 2 O-NaCl-CO 2 at >300°C and depths of ~4 to 6 km near the ductile to brittle transition to H 2 O-NaCl at ~150°C and epithermal depths as the Idaho batholith was exhumed. Geochronology indicates that from the Late Cretaceous to early Eocene, overprinting hydrothermal events formed deep orogenic Au, intermediate W, and shallow epithermal Au-Ag and Sb. Without comprehensive geochronologic constraints, our fluid inclusion data and textural observations could be interpreted as a single evolving hydrothermal system, as has been done in many other mining districts with similar mineralization.

Idaho↗