USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research: Solid Earth”

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 62Linked to original sources

Projections of multiple climate-related coastal hazards for the US Southeast Atlantic

Faced with accelerating sea level rise and changing ocean storm conditions, coastal communities require comprehensive assessments of climate-driven hazard impacts to inform adaptation measures. Previous studies have focused on flooding but rarely on other climate-related coastal hazards, such as subsidence, beach erosion and groundwater. Here, we project societal exposure to multiple hazards along the Southeast Atlantic coast of the United States. Assuming 1 m of sea level rise, more than 70% of the coastal residents and US \ $ 1 trillion in property are in areas projected to experience shallow and emerging groundwater, 15 times higher than daily flooding. Storms increase flooding exposure by an order of magnitude over daily flooding, which could impact up to ~50% of all coastal residents and US \ $ 770 billion in property value. The loss of up to ~80% of present-day beaches and high subsidence rates that currently affect over 1 million residents will exacerbate flooding and groundwater hazard risks.

Florida, Georgia, North Carolina, South Carolina,

Biodegradation of chlorinated ethenes at a karst site in middle Tennessee

This report presents results of field and laboratory investigations examining the biodegradation of chlorinated ethenes in a karst aquifer contaminated with trichloroethylene (TCE). The study site, located in Middle Tennessee, was selected because of the presence of TCE degradation byproducts in the karst aquifer and available site hydrologic and chlorinated-ethene information. Additional chemical, biological, and hydrologic data were gathered to evaluate whether the occurrence of TCE degradation byproducts in the karst aquifer was the result of biodegradation within the aquifer or simply transport into the aquifer. Geochemical analysis established that sulfate-reducing conditions, essential for reductive dechlorination of chlorinated solvents, existed in parts of the contaminated karst aquifer. Other areas of the aquifer fluctuated between anaerobic and aerobic conditions and contained compounds associated with cometabolism, such as ethane, methane, ammonia, and dissolved oxygen. A large, diverse bacteria population inhabits the contaminated aquifer. Bacteria known to biodegrade TCE and other chlorinated solvents, such as sulfate-reducers, methanotrophs, and ammonia-oxidizers, were identified from karst-aquifer water using the RNA-hybridization technique. Results from microcosms using raw karst-aquifer water found that aerobic cometabolism and anaerobic reductive-dechlorination degradation processes were possible when appropriate conditions were established in the microcosms. These chemical and biological results provide circumstantial evidence that several biodegradation processes are active in the aquifer. Additional site hydrologic information was developed to determine if appropriate conditions persist long enough in the karst aquifer for these biodegradation processes to be significant. Continuous monitoring devices placed in four wells during the spring of 1998 indicated that pH, specific conductance, dissolved oxygen, and oxidation-reduction potentials changed very little in areas isolated from active ground-water flow paths. These stable areas in the karst aquifer had geochemical conditions and bacteria conducive to reductive dechlorination of chlorinated ethenes. Other areas of the karst aquifer were associated with active ground-water flow paths and fluctuated between anaerobic and aerobic conditions in response to rain events. Associated with this dynamic environment were bacteria and geochemical conditions conducive to cometabolism. In summary, multiple lines of evidence developed from chemical, biological, and hydrologic data demonstrate that a variety of biodegradation processes are active in this karst aquifer.

Tennessee

Status of water-quality conditions in the United States, 2010–20

Degradation of water quality can make water harmful or unusable for humans and ecosystems. Although many studies have assessed the effect of individual constituents or narrow suites of constituents on freshwater systems, no consistent, comprehensive assessment exists over the wide range of water-quality effects on water availability. Using published studies, data, and models completed at regional or national scales in the United States during 2010–20, this chapter moves towards a comprehensive assessment by summarizing how selected anthropogenic and geogenic water-quality constituents affect national-scale water availability for human and ecosystem needs. Several types of human health, agricultural, ecological, and beneficial-use standards or thresholds were used to provide context for categorizing surface-water and groundwater quality. Water availability for human and ecological use is limited by elevated concentrations of geogenic and anthropogenic constituents in surface and groundwater. Elevated concentrations of five geogenic constituents (arsenic, manganese, strontium, radium, and adjusted gross alpha) are common in groundwater and collectively affect the drinking water supply to over 30 million people. Surface water sourced drinking water supplies are impaired in about a third of assessed stream miles, most commonly because of non-mercury metals and salinity. Health-based violations at community water systems may disproportionately affect socially vulnerable communities. Ecological water uses are predominantly limited by nutrients, sediment, temperature, pathogens, salinity, and pesticides. Water availability for human and ecological use is adversely affected by human activities including human contaminant sources (for example, wastewater, agriculture), processes (for example, dredging, groundwater pumping), or permanent landscape modifications (for example, dams, urbanization). Primary contaminant sources vary spatially and include fertilizer and manure, atmospheric deposition, wastewater treatment plants, urban land, and a range of natural sources. Contaminants of emerging concern, contaminants without regulatory thresholds, and mixtures of geogenic and anthropogenic water contaminants also contribute to ecological degradation and human exposure.

conterminous United States

The projected exposure and response of a natural barrier island system to climate-driven coastal hazards

Accelerating sea level rise (SLR) and changing storm patterns will increasingly expose barrier islands to coastal hazards, including flooding, erosion, and rising groundwater tables. We assess the exposure of Cape Lookout National Seashore, a barrier island system in North Carolina (USA), to projected SLR and storm hazards over the twenty-first century. We estimate that with 0.5 m of SLR, 47% of current subaerial barrier island area would be flooded daily, and the 1-year return period storm would flood 74%. For 20-year return period storms, over 85% is projected to be flooded for any SLR. The modelled groundwater table is already shallow (< 2 m deep), and while projected to shoal to the land surface with SLR, marine flooding is projected to overtake areas with emergent groundwater. Projected shoreline retreat reaches an average of 178 m with 1 m of SLR and no interventions, which is over 60% of the current island width at narrower locations. Compounding these hazards is subsidence, with one-third of the study area currently lowering at > 2 mm/yr. Our results demonstrate the difficulty of managing natural barrier systems such as those managed by federal park systems tasked with maintaining natural ecosystems and protecting cultural resources.

North Carolina

Using geologic mapping to understand temporal and spatial relations of closely clustered to concurrent latest Holocene surface ruptures on two intersecting faults, south-central Mojave Desert, California

The Pinto Mountain Fault Zone (PMFZ) marks a major structural boundary between east-oriented sinistral faults of the eastern Transverse Ranges (to the south) and northwest-oriented dextral faults of the south-central Mojave Desert (to the north). These structural fault systems comprise sinistral and dextral deformational domains of the Eastern California Shear Zone (ECSZ) that intersect one another in the Copper Mountain and Twentynine Palms areas. The U.S. Geological Survey (USGS) is conducting detailed geologic mapping and geochronologic investigations designed to clarify geometric, kinematic, and temporal relations among the two domains, that are focused on the central portion of the left-lateral PMFZ near its intersection with the right-lateral Copper Mountain Fault (CMF) and Mesquite Lake Fault Zone (MLFZ).

California

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas

High‐resolution surface deformation and slip distribution observations for the 2023 Kahramanmaraş, Türkiye, earthquake sequence help constrain the rupture process

Splay, or branch, faults are a common geometric feature of earthquake surface ruptures and may provide constraints on the rupture behavior of an earthquake. The 2023 M w 7.8 Pazarcık and M w 7.5 Elbistan, Türkiye, earthquakes are examples of ruptures with multiple small splays, and the Pazarcık earthquake nucleated on a splay fault, the Narlı fault, before rupturing bilaterally on the East Anatolian fault (EAF). Here, we present 3‐m‐resolution surface displacement from subpixel correlation of Planet Dove optical images for the entirety of both ruptures with corresponding surface slip distributions. For a 30‐km‐long study region spanning the Narlı‐EAF intersection, we compare surface slip derived from five data sets with different resolutions (on‐the‐ground, WorldView, Planet Dove, Sentinel‐2, and Sentinel‐1) to elucidate complementary information. In addition, we integrate information from the surface expression of faulting with published dynamic rupture simulations and rupture process studies to constrain a rupture evolution for the Pazarcık earthquake that is consistent across data sets. This work highlights the complementary nature of disparate surface slip data sets and the role that high‐resolution surface displacement information, including from fault splays, can play in constraining nonunique rupture models and refining understanding of the earthquake rupture process.

Bulletin of the Seismological Society of America

Preferential groundwater discharges along stream corridors are disregarded sources of greenhouse gases

Groundwater delivery of greenhouse gases (GHGs) to stream banks and riparian areas, before mixing with surface waters, has not been well quantified. We measured preferential groundwater delivery of GHGs to stream banks within three stream reaches, and found that stream banks with discharging groundwater emitted more CO 2 and were sources of N 2 O compared to stream banks without actively discharging groundwater, which emitted less CO 2 and were N 2 O sinks. At one of our stream reaches, groundwater CO 2 and N 2 O concentrations were 1.4–19.2 and 1.1–40.6 times higher than those in surface water, respectively, and groundwater delivery rates of CO 2 and N 2 O were 1.5 and 1.6 times higher than surface water emissions per unit area. On average, 21% (range 0%–100%) of CO 2 and N 2 O were emitted at the stream bank before mixing with surface waters. Preferential groundwater GHG emissions may contribute substantially to stream corridor emissions and may be underestimated when using a channel-centric approach to estimate riverine GHG budgets.

JGR Biogeosciences

Bathymetric maps, surface areas, and storage capacities of Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma, 2024

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, completed high-resolution multibeam bathymetric surveys to compute new elevation-area-capacity tables for Council Grove Lake and Marion Reservoir, Kansas, and Pine Creek Lake, Oklahoma. Elevation-area-capacity tables identify the relation between the water-surface elevation, surface area, and storage capacity of the lake. The surface areas and storage capacities of each lake were computed from bathymetric surfaces combining multibeam echo sounder data collected in 2024 and light detection and ranging point-cloud data collected in 2016 and 2018.

Kansas, Oklahoma

Latitudinal gradients of snow contamination in the Rocky Mountains associated with anthropogenic sources

Seasonal snow is an important source of drinking water and recreation, and for agriculture in the Rocky Mountain region. Monitoring snow-water quality can inform on the effects to the albedo and energy balance of the snowpack, and the sources of natural and anthropogenic aerosol and gases. This study analyzed metals in the seasonal snowpack from water year (WY) 2018 for 49 sites. Calcium, lanthanum, and cerium concentrations support the importance of mineral dust to the southern Rocky Mountains. Mercury (Hg), zinc (Zn), and cadmium (Cd) concentrations showed a similar spatial pattern to mineral dust, whereas antimony (Sb) concentrations were highest in the northern Rocky Mountains. To assess the relative contributions from dust versus anthropogenic contaminant sources, enrichment factors (EF) were calculated, with values above 10 indicating anthropogenic contamination. For Cd, Hg, Sb, and Zn, EF values exceeded 10 at northern sites. These observations were compared to spatial trends of EF values of Hg from WY2009 to WY2018, regional monitoring networks, and back trajectory analyses. The agreement between these datasets revealed temporally consistent contaminant sources and/or transport processes to the northern Rocky Mountains snowpack. Sources include current and historical mining and smelting in the region. Strategies to limit the emissions of these metals to the Northern Rockies could benefit from focusing on remediation of contaminated sites, and continued monitoring and mitigation of active mining and smelting.

Colorado , Idaho, Montana, New Mexico, Utah, Wyomi

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund Site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

EarthArXiv

Magnitude and frequency of low flows in the Suwannee River Water Management District, Florida

Low-flow frequency statistics for 20 gaging stations having at least 10 years of continuous record and 31 other stations having less than 10 years of continu ous record or a series of at least two low- flow measurements are presented for unregulated streams in the Suwannee River Water Management District in north-central Florida. Statistics for the 20 continuous-record stations included are the annual and monthly minimum consecutive-day average low- flow magnitudes for 1, 3, 7, 14, and 30 consecutive days for recurrence intervals of 2, 5, 10, 20, and, for some long-term stations, 50 years, based on records available through the 1994 climatic year.Only theannual statistics are given for the 31 other stations; these are for the 7- and 30-consecutive day periods only and for recurrence intervals of 2 and 10 years only. Annual low-flow frequency statistics range from zero for many small streams to 5,500 cubic feet per second for the annual 30- consecutive-day average flow with a recurrenceinterval of 2 years for the Suwannee River near Wilcox (station 02323500). Monthly low-flow frequency statistics range from zero for many small streams to 13,800 cubic feet per second for the minimum 30-consecutive-day average flow with a 2-year recurrence interval for the month of March for the same station. Generally, low-flow characteristics of streams in the Suwannee River Water Management District are controlled by climatic, topographic, and geologic fac tors. The carbonate Floridan aquifer system underlies, or is at the surface of, the entire District. The terrane's karstic nature results in many sinkholes and springs. In some places, springs may contribute greatly to low streamflow and the contributing areas of such springs may include areas outside the presumed surface drainage area of the springs. In other places, water may enter sinkholes within a drainage basin, then reappear in springs downstream from a gage. Many of the smaller streams in the District go dry or have no flow for several months in many years. In addition to the low-flow statistics, four synoptic low-flow measurement surveys were conducted on 161 sites during 1990, 1995, and 1996. The measurements were made to provide "snapshots" of flow conditions of streams throughout the Suwannee River Water Management District. Magnitudes of low flows during the 1990 series of measurements were in the range associated with minimum 7-consecutive-day 50-year recurrence interval to the minimum 7-consecutive-day 20-year recurrence interval, except in Taylor and Dixie Counties, where the magnitudes ranged from the minimum 7-consecutive-day 5-year flow level to the7-consecutive-day 2-year flow level. The magnitudes were all greater than the minimum 7- consecutive-day 2-year flow level during 1995 and 1996. Observations of no flow were recorded at many of the sites for all four series of measurements.

Florida

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California

Ages and trace element fertility of porphyry-related mineralization in the Philipsburg polymetallic district, Montana, with a comparison to Butte

The Philipsburg mining district is a Mo-Cu porphyry system with associated Cordilleran polymetallic veins. Geochronology was employed to date the porphyry (~66 Ma, U/Pb in zircon) and molybdenite mineralization from the veins (~76 Ma, Re-Os). Age results suggest that the two-mineralization events model proposed by Lund et al. (2018) for the Butte district can be applied to the Philipsburg deposit. Furthermore, common fertility indicators from zircon trace elements were analyzed with variable success to characterize the existing Mo-Cu mineralization of the Philipsburg porphyry. This raises concerns about the benefits of this method and/or the importance of sampling, which could impact exploration for similar porphyry deposits.

Montana

High-resolution landscape-scale maps of nutrients, phytoplankton, and related water-quality constituents in Suisun Bay and the Sacramento–San Joaquin Delta, California, during 2018—Results from high-resolution underway surveys

We examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun, Grizzly, and Honker Bays in California by completing three spatial surveys that used continuous underway high frequency (1 hertz) measurements and sampling onboard a high-speed boat. Surveys were conducted in May, July, and October 2018. Water was sampled continuously during surveys to simultaneously collect information about the concentration and spatial distribution of all major nutrient forms together with information about the major classes of phytoplankton and associated standard field measurements of water quality, such as dissolved oxygen, dissolved organic matter, pH, salinity, specific conductance, turbidity, and water temperature. The results show a greater than 50-fold variation in nutrient and phytoplankton concentrations across space and time, providing evidence of the dynamic environmental processes that shape the ways nutrients interact with and affect Delta aquatic habitats.

Data Report

Polystyrene microplastics alter the accumulation and elimination dynamics of silver nanoparticles in Daphnia magna

Rationale Microplastics (MPs) can interact with engineered nanomaterials and alter their environmental fate and bioavailability. However, their influence on the bioaccumulation dynamics of silver nanoparticles (AgNPs) in aquatic filter feeders under environmentally relevant conditions remains poorly understood. We hypothesized that polystyrene microplastics (PS-MPs) alter the uptake, elimination and overall bioaccumulation dynamics of AgNPs in Daphnia magna . Methodology Adult D. magna were exposed to isotopically labeled citrate-coated 109 AgNPs in the presence and absence of PS-MPs. Waterborne uptake, dietary uptake, assimilation efficiency, food ingestion and elimination were quantified experimentally and incorporated into a biodynamic model to predict steady-state silver (Ag) concentrations under environmentally relevant exposure scenarios. Results PS-MPs increased waterborne Ag accumulation by approximately eightfold compared with AgNPs exposure alone. However, Ag elimination was substantially faster, with 98% of accumulated Ag eliminated after five days of depuration compared with 66% without PS-MPs. PS-MPs also reduced food ingestion rates (IRs), while assimilation efficiency remained largely unchanged. Biodynamic modeling predicted that steady-state Ag accumulation was approximately fivefold greater in the presence of PS-MPs, with waterborne exposure becoming the dominant accumulation pathway. Discussion These findings suggest that PS-MPs alter bioaccumulation dynamics of AgNPs in filter-feeding organisms by enhancing organism-associated Ag during waterborne exposure while accelerating Ag elimination, likely through particle-associated transport and gut egestion. Collectively, these results underscore the complexity of Ag bioaccumulation processes when MPs and AgNPs interact and highlight the importance of studies conducted under environmentally relevant conditions.

Environmental Chemistry

Rewilding the upper Klamath River Basin: Rapid recolonization of Chinook salmon following the world’s largest dam removal

Habitat fragmentation is a global driver of decline in biodiversity, yet evaluations of ecological response following re-connection remain rare. In 2024, four hydroelectric dams on the Klamath River were removed, reconnecting over 640 km of anadromous habitat and marking the largest dam removal project completed to date. We used a stationary high-resolution multi-beam imaging SONAR to quantify adult Chinook salmon ( Oncorhynchus tshawytscha ) passage into newly reopened reaches during the first two migratory seasons following dam removal. An estimated 7,742 (95% CI: 7,702–7,778) Chinook salmon migrated upstream in 2024 and 13,310 (95% CI: 12,876–13,733) in 2025, representing 18–19% of fish returning to the Klamath River Basin. Within this period, Chinook salmon recolonized 88% of their documented historical distribution, expanding into 345 km of reconnected habitats and reaching elevations of up to 1,250 m. Recolonization required fish to swim through dam removal sites, complex hydraulic corridors and engineered fishways, demonstrating the propensity for rapid recolonization at basin scale despite these challenges. These findings provide the first quantitative assessment of salmon response to the world’s largest dam removal and demonstrate how reconnection can catalyze large-scale recolonization in historically fragmented river networks.

California, Oregon

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado