USGS Science⌕ Search

SEARCH · USGS Science

Results for “Applied and Environmental Microbiology”

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

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

Arizona↗
USGS2025

Updating regional‐scale geospatial liquefaction models with locally available geotechnical data

We present a method to update the geospatial liquefaction model used by the U.S. Geological Survey’s near‐real‐time ground failure product with subsurface geotechnical data. The geospatial model estimates liquefaction probability from peak ground velocity (via ShakeMap) and geospatial susceptibility proxies. In many regions, additional information relevant to constraining liquefaction likelihood is also available, including surface geology maps and subsurface geotechnical measurements. There is currently no mechanism to use these data in the ground failure product liquefaction model, even though these data could provide more precise constraints on spatial variations in the lithologic character of the soil (surface geology) and direct measurements of the subsurface mechanical properties that affect liquefaction occurrence and severity (geotechnical measurements). In this study, we develop a method to integrate these data with the geospatial model and assess how these data can improve regional‐scale predictions. We develop a Bayesian updating framework and apply it to the 1989 magnitude 6.9 Loma Prieta, California, earthquake, for which mapped observations are available to evaluate performance. We constrain the Bayesian framework with 373 Northern California cone penetration tests and liquefaction susceptibility classes based on the mapped surface geology. This Bayesian model incorporates geotechnical information into the geospatial model and more accurately predicts liquefaction occurrences than the geospatial model, while sacrificing less accuracy in terms of predicting the absence of liquefaction than the geotechnical model. In future applications, this approach could be adapted to update other geospatial models using locally available subsurface data.

California↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2024

Estimating domestic self-supplied water use in Rhode Island, 2014–21

Water withdrawal from private groundwater wells is often unaccounted for in water planning studies, and water from private wells can be a source of exposure to environmental contaminants. The sizes of populations that depend on private wells for domestic water use and the amounts of water that are withdrawn from these wells are generally poorly represented in data collection efforts because of the challenges of locating, metering, or gathering withdrawal information from individual property owners. To address this problem, the U.S. Geological Survey, in cooperation with the Rhode Island Water Resources Board, estimated the volume of water withdrawn from domestic self-supply wells and the populations who use them for the State of Rhode Island at a 30-meter pixel spatial resolution and one-month temporal resolution between July 2014 and June 2021. The number of people reliant on domestic self-supply wells has increased in Rhode Island over the study period; however, the statewide estimate of total water withdrawal has not statistically increased. Withdrawals from private wells are largest in the inland areas of the western part of the State, and the towns of Scituate and Charlestown have the highest estimated withdrawals. Statewide monthly withdrawals ranged from 3.987 million gallons per day in March 2018 to 7.767 million gallons per day in September 2016. The median per capita domestic water use rate was 46.0 gallons per capita per day.

Rhode Island↗
USGS2026

Continuous monitoring and temporal variability of fluorescence of dissolved organic matter in karst groundwater of the Edwards aquifer, south-central Texas, 2019–24

Widespread urbanization on the Edwards aquifer recharge zone in south-central Texas has prompted the development of new approaches and tools for evaluating the current and future status of water quality in the San Antonio and Barton Springs segments of the Edwards aquifer. The U.S. Geological Survey, in cooperation with the San Antonio Water System and the City of Austin, applied continuous monitoring of fluorescence of dissolved organic matter (fDOM) to characterize the sources and transport of dissolved organic matter and associated constituents in the Edwards aquifer. Continuous fDOM time-series data were adjusted using empirically derived temperature, turbidity, and inner-filter-effect correction factors prior to analysis and interpretation. Fully corrected fDOM time-series data from the study sites were evaluated for short- and long-term temporal variability in fDOM in the context of the hydrologic and climatic conditions that occurred during the 5-year study period to understand vulnerability of the aquifer to potential contaminants in recharge. Study results indicate that fDOM is an effective proxy for dissolved organic carbon that, in turn, is indicative of the influx of recent surface water and associated contaminants, and thus, aquifer vulnerability to those contaminants. Although other continuously monitored water-quality parameters provide insights into changes in water quality in response to varying hydrologic conditions, the quantitative relation between fDOM and dissolved organic carbon directly indicates the timing and magnitude of pulses of organic constituents derived from the land surface, whereas other water-quality parameters do not.

Texas↗
USGS2024

Deep structure of Siletzia in the Puget Lowland: Imaging an obducted plateau and accretionary thrust belt with potential fields

Detailed understanding of crustal components and tectonic history of forearcs is important due to their geological complexity and high seismic hazard. The principal component of the Cascadia forearc is Siletzia, a composite basaltic terrane of oceanic origin. Much is known about the lithology and age of the province. However, glacial sediments blanketing the Puget Lowland obscure its lateral extent and internal structure, hindering our ability to fully understand its tectonic history and its influence on modern deformation. In this study, we apply map-view interpretation and two-dimensional modeling of aeromagnetic and gravity data to the magnetically stratified Siletzia terrane revealing its internal structure and characterizing its eastern boundary. These analyses suggest the contact between Siletzia (Crescent Formation) and the Eocene accretionary prism trends northward under Lake Washington. North of Seattle, this boundary dips east where it crosses the Kingston arch, whereas south of Seattle the contact dips west where it crosses the Seattle uplift (SU). This westward dip is opposite the dip of the Eocene subduction interface, implying obduction of Siletzia upper crust at this southern location. Elongate pairs of high and low magnetic anomalies over the SU suggest imbrication of steeply-dipping, deeply rooted slices of Crescent Formation within Siletzia. We hypothesize these features result from duplication of Crescent Formation in an accretionary fold-thrust belt during the Eocene. The active Seattle fault divides this Eocene fold-thrust belt into two zones with different structural trends and opposite frontal ramp dips, suggesting the Seattle fault may have originated as a tear fault during accretion.

Washington↗
USGS2024

A conceptual framework to assess post-wildfire water quality: State of the science and knowledge gaps

Wildfire substantially alters aquatic ecosystems by inducing moderate to catastrophic physical and chemical changes. However, the relations of environmental and watershed variables that drive those effects are complex. We present a Driver-Factor-Stressor-Effect (DFSE) conceptual framework to assess the current state of the science related to post-wildfire water-quality. We reviewed 64 peer-reviewed papers using the DFSE framework to identify drivers, factors, stressors, and effects associated with each study. A total of five drivers were identified and ranked according to their frequency of occurrence in the literature: atmospheric processes > fire characteristics > ecologic processes and characteristics > land surface characteristics > soil characteristics. Commonly reported stressors include increased nutrients, runoff, and sediment transport. Furthermore, although several different factors have been used at least once to explain water-quality effects, relatively few factors outside of precipitation and fire characteristics are frequently studied. We identified several gaps indicating the need for long-term monitoring, multi-factor studies, consideration of organic contaminants, consideration of groundwater, and inclusion of soil characteristics. This assessment expands on other reviews and meta-analyses by exploring causal linkages between influential variables and overall effects in post-wildfire watersheds. Information gathered from our assessment and the framework itself can be used to inform future monitoring plans and as a guide for modeling efforts focused on better understanding specific processes or to mitigate potential risks of post-wildfire water quality.

Water Resources Research↗
USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗
USGS1995

The distribution of bromide in water in the Floridan aquifer system, Duval County, northeastern Florida

Although Duval County, Florida, has ample ground-water resources for public supply, the potential exists for a problem with excessive disinfectant by-products. These disinfectant by-products result from the treatment of raw water containing low concentrations of bromide and naturally occurring organic compounds. Because of this potential problem, the relation of bromide concentrations to aquifer tapped, well location and depth, and chemical characteristics of water in the Floridan aquifer system underlying Duval County were studied to determine if these relations could be applied to delineate water with low-bromide concentrations for future supplies. In 1992, water samples from 106 wells that tap the Floridan aquifer system were analyzed for bromide and major dissolved constituents. A comparison of bromide concentrations from the 1992 sampling with data from earlier studies (1979-80) indicates that higher bromide concentrations were detected during the earlier studies. The difference between the old and new data is probably because of a change in analytical methodology in the analysis of samples. Bromide concentrations exceeded the detection limit (0.10 milligrams per liter) in water from 28 of the 106 wells (26 percent) sampled in 1992. The maximum concentration was 0.56 milligrams per liter. There were no relations between bromide and major dissolved constituents, well depth, or aquifer tapped that would be useful for determining bromide concentrations. Areal patterns of bromide occurrence are not clearly defined, but areas with relatively high bromide concentrations tend to be located in a triangular area near the community of Sunbeam, Florida, and along the St. Johns River throughout Duval County.

Florida↗
USGS1976

Non-radiometric borehole geophysical detection of geochemical halos surrounding sedimentary uranium deposits

Roll-type uranium deposits are formed by the concentration of uranium by ground water in geochemical cells. Non-uranium minerals having different solubilities may be deposited ahead of or behind the uranium minerals, forming halos that surround the ore. In addition, oxidizing and reducing environmental conditions may cause zones of mineral alteration to develop beyond the limits of the uranium deposit. Certain physical property anomalies that are commonly associated with halos can be detected by relatively fast and inexpensive borehole geophysical measurements made either in individual holes, or between two adjacent holes. Borehole measurements that have been found to be useful include electrical resistivity, induced polarization, and magnetic susceptibility. Electrical resistivity is increased by the presence of calcite and other cementing minerals that sometimes create permeability barriers in the neighborhood of uranium deposits. Induced polarization (IP) response is increased by sulfide and clay minerals that are commonly found in anomalous concentrations near roll-type deposits. Magnetic susceptibility is usually decreased by the oxidation of magnetite to hematite or limonite in the zone of chemical alteration that is left as a trail behind roll fronts. Borehole measurements of electrical resistivity, induced polarization and magnetic susceptibility were made in the vicinity of a uranium roll-type deposit in south Texas. Results indicate that mineral halos can be detected by borehole measurements made. in wide-spaced drill holes, and that the total amount of drilling needed to find a deposit can be reduced substantially by this exploration approach.

Open-File Report↗
USGS2025

Advancing current understanding of Martian impact-generated hydrothermal systems through novel coupled modeling: Insights from Gale, Jezero, and other craters

Impact-generated hydrothermal systems are considered potentially habitable environments on Mars, Earth, and other planetary bodies for microbial life. However, there is an ongoing debate regarding what geological features on Mars provide definitive evidence for such systems. Although earlier studies have modeled hydrothermal processes in Martian craters, they often lacked integration with shock physics hydrocodes to constrain initial impact conditions. The importance of this two-code coupling was demonstrated by successfully replicating alteration signatures in the Earth's Haughton impact structure. In this study, we use a similar two-code approach, combining the iSALE hydrocode with the HYDROTHERM hydrothermal model to simulate the full evolution of impact-generated hydrothermal systems. We apply this method to craters the size of Jezero (∼50 km) and Gale (∼154 km) in diameter. Although Jezero's interior is largely buried, our results align with hypothesized hydrothermal vents and alteration minerals near central uplifts in similarly sized exposed craters, such as Toro and Auki. Furthermore, our models correspond to alteration patterns observed by the Curiosity in the lower layers of Mount Sharp, which may represent remnants of impact-driven hydrothermal activity. A key finding is that these systems may persist much longer than previously estimated. Our simulations suggest that a Jezero-sized system could remain habitable for thermophiles for approximately 720,000 years, whereas a Gale-sized system could persist for nearly 2 million years. Additionally, simulations under unsaturated crustal conditions reveal that air-dominated near-surface layers can suppress vertical fluid flow, enabling deep subsurface alteration without producing detectable mineral signatures at the surface.

Journal of Geophysical Research - Planets↗
USGS2024

Triggering the 2022 eruption of Mauna Loa

Distinguishing periods of intermittent unrest from the run-up to eruption is a major challenge at volcanoes around the globe. Comparing multidisciplinary monitoring data with mineral chemistry that records the physical and spatio-temporal evolution of magmas fundamentally advances our ability to forecast eruptions. The recent eruption of Mauna Loa, Earth’s largest active volcano, provides a unique opportunity to differentiate unrest from run-up and improve forecasting of future eruptions. After decades of intermittent seismic and geodetic activity over 38 years of repose, Mauna Loa began erupting on 27 November 2022. Here we present a multidisciplinary synthesis that tracks the spatio-temporal evolution of precursory activity by integrating mineral and melt chemistry, fluid inclusion barometry, numerical modeling of mineral zoning, syn-eruptive gas plume measurements, the distribution and frequency of earthquake hypocenters, seismic velocity changes, and ground deformation. These diverse data indicate that the eruption occurred following a 2-month period of sustained magma intrusion from depths of 3–5 km up to 1–2 km beneath the summit caldera, providing a new model of the plumbing system at this very high threat volcano. Careful correlation of both the geochemistry and instrumental monitoring data improves our ability to distinguish unrest from the run-up to eruption by providing deeper understanding of the both the monitoring data and the magmatic system—an approach that could be applied at other volcanic systems worldwide.

Hawaii↗
USGS2025

Factors regulating the potential for freshwater mineral soil wetlands to function as natural climate solutions

There are increasing global efforts and initiatives aiming to tackle climate change and mitigate its impacts via natural climate solutions (NCS). Wetlands have been considered effective NCS given their capacity to sequester and retain atmospheric carbon dioxide (CO 2 ) while also providing a myriad of other ecosystem functions that can assist in mitigating the impacts of climate change. However, wetlands have a dual impact on climate, influencing the atmospheric concentrations of both CO 2 and methane (CH 4 ). The cooling effect associated with wetland CO 2 sequestration can be counterbalanced by the warming effect caused by CH 4 emissions from wetlands. The relative ability of wetlands to sequester CO 2 versus emit CH 4 is dependent on a suite of interacting physical, chemical, and biological factors, making it difficult to determine if/which wetlands are considered important NCS. The fact that wetlands are embedded in landscapes with surface and subsurface hydrological connections to other wetlands (i.e., wetlandscapes) that flow over and through geochemically active soils and sediments adds a new layer of complexity and poses further challenges to understanding wetland carbon sequestration and greenhouse gas fluxes at large spatial scales. Our review demonstrates how additional scientific advances are required to understand the driving mechanisms associated with wetland carbon cycling under different environmental conditions. It is vital to understand wetland functionality at both wetland and wetlandscape scales to effectively implement wetlands as NCS to maximize ecological, social, and economic benefits.

Wetlands↗
USGS2024

Mapping Arundo donax (Arundo cane) with multispectral imagery before, during, and after herbicide treatment along the Rio Grande in Webb County, Texas, 2020–21

Arundo donax , commonly called Arundo cane, giant reed, or Carrizo cane, is an invasive bamboo-like perennial grass common in riparian areas throughout the southwestern United States. In Texas, not only does it negatively affect riparian ecosystems, but it has also become a problem for border security because it reduces visibility along the Rio Grande. To address these problems, in 2015 the Texas State Soil and Water Conservation Board was authorized by the Texas State Legislature to develop a program to eradicate Arundo cane along the Rio Grande. In 2020, the Texas State Soil and Water Conservation Board applied imazapyr and glyphosate herbicides along a 19.3-kilometer reach of the Rio Grande, northwest of Laredo, Texas. The U.S. Geological Survey, in cooperation with the Texas State Soil and Water Conservation Board and the Webb Soil and Water Conservation District, used WorldView-3 Standard high-resolution satellite imagery to map Arundo cane extent along the reach before , during , and after the herbicide-treatment period on June 30, 2020, September 26, 2020, and May 7, 2021, respectively. A maximum likelihood supervised classification analysis was computed on the images to map the spatial extent and estimate the area covered by Arundo cane. The estimated area covered by Arundo cane in the before classification was 1,282,000 square meters, in the during classification was 1,064,000 square meters, and in the after classification was 1,108,000 square meters. The qualitative comparison of the three images shows that there was an overall decrease in vegetation classified as Arundo cane throughout the study area.

Texas↗
USGS2025

Factors influencing landslide occurrence in low-relief formerly glaciated landscapes: Landslide inventory and susceptibility analysis in Minnesota, USA

In landscapes recently impacted by continental glaciation, landslides may occur where topographic relief has been generated by the drainage of glacial lakes and ensuing post-glacial fluvial network development into unconsolidated glacially derived sediments and exhumed bedrock. To investigate linkages among environmental variables, post-glacial landscape development, and landslides, we created a landslide inventory of nearly 10,000 landslides in five regions of the formerly glaciated low-relief state of Minnesota, USA. Multivariate logistic regression indicates the importance of slope angle, lithology, and the development of stream valleys to landslide distribution. Areas underlain by fine-grained glaciolacustrine and nearshore deposits that are incised by streams are particularly prone to shallow (<1-2 m depth) landslides. Landslides also occur in a wide range of glacial and fluvial deposits, and as rockfall in layered Paleozoic sedimentary rocks in central and southern Minnesota and Precambrian igneous and sedimentary rocks in northeastern Minnesota. Although no more than 1-2% of the studied regions are susceptible to landslides, they can pose risk to life and safety, damage infrastructure, and impact water quality. The combination of recently generated low-relief steep slopes, extensive unconsolidated sediments, and layered sedimentary bedrock make this formerly glaciated landscape more susceptible to landslides than current national-scale models indicate.

Minnesota↗
USGS2025

Speleothem evidence for Late Miocene extreme Arctic amplification – An analogue for near-future anthropogenic climate change?

The Miocene provides an excellent climatic analogue for near-future runaway anthropogenic warming, with atmospheric CO 2 concentrations and global average temperatures similar to those projected for the coming century under extreme-emissions scenarios. However, the magnitude of Miocene Arctic warming remains unclear due to the scarcity of reliable proxy data. Here we use stable oxygen isotope and trace element analyses, alongside clumped isotope and fluid inclusion palaeothermometry of speleothems to reconstruct palaeo-environmental conditions near the Siberian Arctic coast during the Tortonian (8.68 ± 0.09 Ma ). Stable oxygen isotope records suggest warmer-than-present temperatures. This is supported by temperature estimates based on clumped isotopes and fluid inclusions giving mean annual air temperatures between + 6.6 and + 11.1 °C , compared with − 12.3 °C today. Trace elements records reveal a highly seasonal hydrological environment. Our estimate of > 18 °C of Arctic warming supports the wider consensus of a warmer-than-present Miocene and provides a rare palaeo-analogue for future Arctic amplification under high-emissions scenarios. The reconstructed increase in mean surface temperature far exceeds temperatures projected in fully coupled global climate models, even under extreme-emissions scenarios. Given that climate models have consistently underestimated the extent of recent Arctic amplification, our proxy data suggest Arctic warming may exceed current projections.

Lena River delta↗
USGS2026

Escherichia coli monitoring and assessment in 2022 and 2023 after beach restoration at Lake St. Clair Metropark Beach, Macomb County, Michigan

Lake St. Clair Metropark Beach in Michigan has a history of closures because of elevated Escherichia coli ( E. coli ) concentrations in its recreational waters. To reduce closures, restoration projects were implemented in 2021 to deter waterfowl from congregating on the beach. In this study, the U.S. Geological Survey, in cooperation with the Michigan Department of the Environment, Great Lakes, and Energy and in collaboration with Huron-Clinton Metroparks and the Macomb County Health Department, monitored E. coli from 2022–23 in surface water, shallow groundwater, and sediment at Lake St. Clair Metropark Beach. Results were compared to data from a prerestoration (2018–19) study. A significant decrease in daily geometric mean E. coli concentrations in surface water was observed postrestoration, but the number of high concentration events increased. This resulted in more frequent beach closures postrestoration. Surface-sediment E. coli concentrations significantly decreased after restoration, and waterfowl populations generally decreased from 2021 to 2023, suggesting that the deterrence measures could be influencing E. coli concentrations in surface sediments and surface water. Groundwater E. coli concentrations were orders of magnitude higher than those in surface water and revealed no change correlated with restoration. Seepage measurements indicated that groundwater occasionally discharges into surface water, potentially providing a transport mechanism for E. coli to reach the lake. Continued monitoring and consideration of environmental factors could help to better understand the beach system.

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