Next Article in Journal
Shoreline Solutions: Guiding Efficient Data Selection for Coastal Risk Modeling and the Design of Adaptation Interventions
Previous Article in Journal
Evaluation of the WRF Model to Simulate a High-Intensity Rainfall Event over Kampala, Uganda
Previous Article in Special Issue
Urban Groundwater Contamination by Non-Steroidal Anti-Inflammatory Drugs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural and Anthropogenic Geochemical Tracers to Investigate Residence Times and Groundwater–Surface-Water Interactions in an Urban Alluvial Aquifer

U.S. Geological Survey, Colorado Water Science Center, Denver Federal Center, P.O. Box 25046, MS-415, Denver, CO 80225, USA
*
Author to whom correspondence should be addressed.
Water 2021, 13(6), 871; https://doi.org/10.3390/w13060871
Submission received: 25 February 2021 / Revised: 19 March 2021 / Accepted: 20 March 2021 / Published: 23 March 2021
(This article belongs to the Special Issue Urban Groundwater)

Abstract

:
A multi-component geochemical dataset was collected from groundwater and surface-water bodies associated with the urban Fountain Creek alluvial aquifer, Colorado, USA, to facilitate analysis of recharge sources, geochemical interactions, and groundwater-residence times. Results indicate that groundwater can be separated into three distinct geochemical zones based on location within the flow system and proximity to surface water, and these zones can be used to infer sources of recharge and groundwater movement through the aquifer. Rare-earth-element concentrations and detections of wastewater-indicator compounds indicate the presence of effluent from wastewater-treatment plants in both groundwater and surface water. Effluent presence in groundwater indicates that streams in the area lose to groundwater in some seasons and are a source of focused groundwater recharge. Distributions of pharmaceuticals and wastewater-indicator compounds also inform an understanding of groundwater–surface-water interactions. Noble-gas isotopes corroborate rare-earth-element data in indicating geochemical evolution within the aquifer from recharge area to discharge area and qualitatively indicate variable groundwater-residence times and mixing with pre-modern groundwater. Quantitative groundwater-residence times calculated from 3H/3He, SF6, and lumped-parameter modeling generally are less than 20 years, but the presence of mixing with older groundwater of an unknown age is also indicated at selected locations. Future investigations would benefit by including groundwater-age tracers suited to quantification of mixing for both young (years to decades) and old (centuries and millennia) groundwater. This multi-faceted analysis facilitated development of a conceptual model for the investigated groundwater-flow system and illustrates the application of an encompassing suite of analytes in exploring hydrologic and geochemical interactions in complex systems.

1. Introduction

Groundwater is a major source of water used for domestic and industrial purposes worldwide, particularly in arid and semi-arid regions. Alluvial aquifers are important in many such regions because of their generally high transmissivity, correspondingly large production rates, and accessibility for public-water supplies. Because many alluvial aquifers are unconfined and located along present-day streams, they also can be affected by land use and anthropogenic chemical sources [1], possibly necessitating treatment of groundwater prior to use as a drinking-water supply. In addition to contaminants of anthropogenic origin that have been studied for decades, such as NO3, contaminants of emerging concern (CECs) are becoming increasingly identified in alluvial aquifers and other settings [2,3]. These CECs can include fuel components, solvents, pharmaceutical and wastewater-indicator compounds [2,4], and per- and polyfluoroalkyl substances (PFAS) [5,6].
There are numerous techniques for evaluating the occurrence and distribution of solutes in the environment, including spatiotemporal trend analysis [7], stable isotopes of water and nitrogen [8], and statistical methods [9,10]. Although each of these methods provides insight into geochemical processes occurring in an aquifer, additional knowledge from environmental tracers can provide independent evidence to define recharge conditions, groundwater–surface-water interactions, groundwater mixing, and residence times [11,12,13,14], all of which help to control solute distributions in an aquifer.
This study used a large geochemical dataset (approximately 200 total parameters) to characterize recharge and solute sources, geochemical reactions, flow paths, and groundwater–surface-water interactions in an alluvial aquifer in Colorado, USA. This study was proposed, in part, because of the occurrence of PFAS compounds in the Fountain Creek alluvial aquifer [5]. Geochemical constituents used in the analysis include major ions, trace elements, rare earth elements (REEs), stable isotopes, noble-gas concentrations and isotopes, and pharmaceuticals and wastewater-indicator compounds. Recent studies applying subsets of analytes used in this study have highlighted the benefits of using various tracers in complex hydrologic systems [7,14]. This study builds on previous work, by combining analyses of multiple types of tracers that have not yet been utilized together in the literature (e.g., REEs, noble-gas isotopes, and pharmaceuticals). For example, Reference [15] analyzed noble gases, nitrate, and stable isotopes of water, to investigate sources of water to wells in a mountainous environment, whereas Reference [16] used pharmaceutical compounds and tritium content across the United States, to evaluate occurrence and distribution of anthropogenic compounds in groundwater, as did Reference [7], using REEs in a near-stream aquifer. However, previous studies, such as those noted, used less comprehensive datasets than those described here (e.g., Reference [7] does not include an evaluation of groundwater age). This study is novel in that a set of tracers is applied which allows a more complete understanding of complex aspects of this hydrologic system including the interaction and mixing of groundwater and surface water with a range of ages and compositions. Understanding the ranges in groundwater residence times is critical in this study area because of transient histories of potential contaminant loading from several distinct sources. One end goal of this analysis is to further inform ongoing groundwater modeling in the area, and the results provide an example of the benefit of an integrated approach to investigating sources of recharge for an alluvial aquifer.
The study area is the Fountain Creek alluvial aquifer, located near the city of Colorado Springs, Colorado, USA (Figure 1), on the eastern margin of the Southern Rocky Mountains. The alluvial aquifer is primarily composed of Quaternary-age sand and gravel fluvial sediments derived from weathering of Precambrian granitic rocks to the west and Mesozoic sedimentary rocks to the east [17,18]. The unconfined aquifer has been the subject of previous studies that focused on sources of solutes in the interconnected groundwater–surface-water system, as well as on groundwater-flow directions and water budgets [17,18,19,20,21,22]. Previous studies generally found that NO3 was a constituent of concern and was derived from a variety of sources including fertilizers, effluent from wastewater treatment plants (WWTPs), and septic systems [8]. Work by Reference [22] indicated the influence of WWTPs on solute concentrations and flow volume of Fountain Creek and determined that resource extraction activities in the area contributed dissolved metals (Cu, Mn, Ni, Pb, and Zn) to the stream. In addition to anthropogenic compounds, the aquifer has naturally elevated concentrations of Se and U derived from the Cretaceous-age Pierre Shale, which underlies the alluvial aquifer [17]. Recently, concentrations of PFAS compounds greater than US Environmental Protection Agency (EPA) lifetime health advisory level of 70 nanograms per liter (ng/L) [23] have been found in the alluvial aquifer [5]. Potential sources of PFAS compounds include WWTPs, aqueous-film-forming foam (AFFF) used in firefighting operations, and metal-plating operations [6,10], each of which is documented in the area.
Hydrologically, Fountain Creek interacts with the shallow unconfined groundwater system by both gaining and losing flow depending on the relative elevation difference between the stream and the water table [17,20]. Effluent from WWTPs contributes a substantial portion of flow to the creek during base-flow periods [22]. In the study reach, there are diversions from the creek for agricultural purposes, and these diversions, as well as streamflow, are simulated by using a transit-loss model [29], which is used by water managers for water-supply planning [30]. Streamflow in the area generally peaks in late spring and early summer in response to snowmelt runoff from the mountainous areas to the west (based on streamflow measured at US Geological Survey site number 07105500, Fountain Creek at Colorado Springs, Colorado) [24]. Groundwater flow within the central portion of the alluvial aquifer is generally from north to south consistent with the topographic gradient and streamflow direction [18]. Several small paleochannels are saturated to the northeast of the central portion of the aquifer (and distant from Fountain Creek). These paleochannels have been termed the tributary alluvium and are hydraulically and topographically upgradient from the main valley of the Fountain Creek aquifer (Figure 1). Groundwater flow between the tributary alluvium and the main aquifer appears to be contiguous based on water-level observations (see Supplemental Materials).
The study area is semi-arid with a median annual precipitation rate of approximately 430 millimeters per year, determined based on gridded meteorological datasets based on the period from 1990 to 2020 [31]. Precipitation is seasonal with the greatest quantity generally falling in the summer as monsoon rainstorms. Because of the semi-arid climate, precipitation recharge is minimal, and much of the groundwater recharge is derived from streamflow losses and irrigation return flow [17,18]. Groundwater pumping in the study area has occurred for decades for domestic, industrial, and agricultural uses [20,32] with combined pumping withdrawals up to 19,000 cubic meters per day [20]. Groundwater pumping of this magnitude for prolonged periods of time has likely altered the flow field creating localized groundwater depressions around pumping wells. Patterns in groundwater-level elevations are complicated by pumping in the area, but wells with less pumping-induced disturbances display seasonal variations controlled by snowmelt runoff recharge, streamflow seasonality, and precipitation seasonality [24].

2. Materials and Methods

2.1. Sample Collection, Analysis, and Quality Control

Water-quality samples were collected from groundwater wells and surface-water features, illustrated in Figure 1, to characterize aquifer conditions, as well as the composition of end-member sources of recharge (i.e., recharge sources with distinctive geochemistry). Groundwater wells sampled for this study included monitoring, public-supply, agricultural-supply, and domestic-supply wells. Groundwater sites were selected based on location within the alluvial aquifer and accessibility. Where possible, dedicated monitoring wells from previous investigations [21] were selected for sampling. Monitoring wells were sampled, using either a portable Grundfos stainless steel pump or a Proactive Mini-Monsoon 12V low-flow pump, and pumps were thoroughly cleaned between wells following standard US Geological Survey (USGS) techniques [33]. Where dedicated monitoring wells were not available, samples were collected from public-supply, agricultural-supply, or domestic-supply wells. For wells with dedicated pumps (i.e., pumps that remain in place), samples were collected from available sample ports prior to any water treatment, thus providing complementary data to samples collected by using removable pumps which were also not subjected to any treatment prior to sample collection. Surface-water sampling sites were selected to characterize the surface-water chemistry of Fountain Creek, as well as Canal 4, which diverts water from Fountain Creek to eastern parts of the study area (Figure 1). Both Fountain Creek and Canal 4 are considered potential sources of groundwater recharge to the alluvial aquifer.
All samples were collected according to procedures described in the USGS National Field Manual [33]. Briefly, groundwater wells were pumped at a low rate until temperature, specific conductance, pH, dissolved oxygen (DO), and turbidity had stabilized to within respectively +/−0.2 ℃, +/−3%, +/−2 standard units (SU), +/−0.3 milligrams per liter (mg/L), and +/−0.5 turbidity units, or until three casing volumes were removed, whichever took longer. This ensured that water being sampled was representative of water within the aquifer. Surface-water samples were collected, using methods for isokinetic sampling [33], although necessary velocity to achieve true isokinetic conditions was not obtained in all instances. Alkalinity and pH were measured in the field for groundwater samples and at the USGS National Water Quality Laboratory (NWQL) in Lakewood, Colorado, for surface-water samples [34].
Samples for analysis of major and trace elements, pharmaceuticals, and wastewater-indicator compounds, and stable isotopes of NO3 were collected in polyethylene bottles and filtered to 0.45 micrometers (µm). Nitrogen-isotope samples were then filtered a second time, using a 0.20 µm filter and frozen. Samples for metals analysis were preserved in the field with nitric acid (HNO3) and refrigerated, and samples for anions were unpreserved but refrigerated. Samples for stable isotopes of water were collected, unfiltered, in glass bottles. Quality-control samples (field blanks and replicates) were collected for trace elements, pharmaceuticals and wastewater-indicator compounds, and stable isotopes of nitrate (replicate only, no blank) to evaluate sampling and/or analytical bias. Tritium (3H) samples were collected, unfiltered, in polyethylene bottles with no headspace by bottom filling after flushing with three volumes of sample water. Samples of the major dissolved gases O2, N2, Ar, CO2, and CH4 were collected, unfiltered, in replicate, in glass bottles sealed with no headspace by bottom filling. Samples of noble gases (He, Ne, Ar, Kr, and Xe) were collected, unfiltered, in replicate, in copper tubes [35] and flushed of all air bubbles before being sealed under positive pressure. Samples of sulfur hexafluoride (SF6) were collected, unfiltered, in replicate, in glass bottles with no headspace by bottom filling after flushing with three volumes of sample water.
Anions were analyzed by ion-exchange chromatography, and cations were analyzed by inductively coupled plasma-atomic emission spectrometry at the NWQL [36,37]. Trace elements were analyzed by inductively coupled plasma-optical emission spectrometry and inductively coupled plasma-mass spectrometry at the USGS Analytical Trace Elements Chemistry Laboratory located in Boulder, Colorado [38]. Pharmaceutical compounds and wastewater-indicator compounds were analyzed by high-performance liquid chromatography tandem mass spectrometry [39] at the NWQL. Stable isotopes of nitrate (δ15N, δ18ONO3) were analyzed by continuous-flow isotope-ratio mass spectrometry (IRMS) [40] at the USGS Reston Stable Isotope Laboratory (RSIL) located in Reston, Virginia, and are reported in standard stable-isotopic units of per mil (‰) [41] with uncertainties of +/−0.5‰. Stable isotopes of water (δ2H, δ18OH2O) were analyzed by dual-inlet IRMS [42,43] at the RSIL with uncertainties of +/−2‰ and +/−0.2‰, respectively. Tritium analysis was conducted by using distillation and electrolytic enrichment, followed by liquid scintillation, at the USGS Menlo Park Tritium Laboratory, located in Menlo Park, California, and are reported in tritium units (TU) with analytical uncertainty ranging from 0.09 to 0.17 TU. Concentrations of major gases were analyzed, using gas chromatography at the USGS Reston Groundwater Dating Laboratory located in Reston, Virginia, with uncertainties of 0.002 mg/L (O2), 0.001 mg/L (N2), 0.003 mg/L (Ar), 0.04 (CO2), and 0.0005 mg/L (CH4). Concentrations and isotopes of noble gases were analyzed at the USGS Noble Gas Laboratory located in Denver, Colorado, using a magnetic-sector mass spectrometer and ultralow vacuum extraction line [44], and are reported in units of cubic centimeters at standard temperature and pressure per gram of water (ccSTP/gH2O) and as ratios (e.g., 20Ne/22Ne), respectively, with analytical uncertainties of 1% (He), 2% (Ne), 2% (Ar), 3% (Kr), and 3% (Xe). Concentrations of SF6 were analyzed at the USGS Groundwater Dating Laboratory, using purge and trap gas chromatography followed by an electron capture detector according to the method of Reference [45] and are reported in units of femtomoles per kilogram (fm/kg) with uncertainty ranging from 20% at the minimum reporting limit to 3% at the maximum reporting limit (https://water.usgs.gov/lab/sf6/lab/analytical_procedures/, accessed on 1 January 2021). All water-quality data are available from the USGS National Water Information System (NWIS) database [24], using the USGS site IDs provided in Supplementary Materials Table S1.
Data quality was evaluated through quality-control samples, calculation of charge balance for major-ion data, and comparison to historical data. This dataset included more than 10 percent field blanks and replicates. Source solution, equipment, and field blanks were collected to test for possible contamination [46]. For the blanks, there were few detections for trace elements where only 7% of these detections were greater than 1 microgram per liter (µg/L). These detections greater than 1 µg/L occurred for Al, Cu, and Zn. The potential bias in samples for these elements was assessed by comparing the maximum concentration reported in blanks to the average environmental concentration. Potential bias of Cu and Zn are generally low given that the maximum concentration detected in the blanks are 6% and 16%, respectively, of the average environmental concentration. For Al, the maximum concentration in blanks was 103% of the average environmental concentration, indicating greater potential for bias and resulting in Al being removed from subsequent statistical analyses. Variability between environmental and replicate samples, quantified by the median of relative percent difference (RPD), for elements with potential bias (Al, Cu, Zn) were respectively 17%, 18%, and 9%, indicating that there is relatively little variability in the dataset. Of the 11 pharmaceutical and wastewater-indicator compounds detected in blanks, none had values above the reporting level. Rare earth elements did not have detections in the blanks with exception of 3 samples, which were greater than 0.002 µg/L for Ce. The maximum concentration of detections in blanks for Ce are 9% of the average environmental concentration, meaning samples are likely not substantially biased for Ce. RPD calculated from field replicates was less than 10 percent for nutrient and major-ion data with exception of 1 replicate pair that had an alkalinity RPD of 14%. Charge balance for the dataset ranged from −10.1 to 12%, where all but three samples had charge balances within the +/−10% range suggested for many water-quality studies [47]. Given the infrequent detections in blanks and generally low RPDs, the dataset is considered representative of environmental conditions in the study area.

2.2. Multivariate Evaluation

Laboratory analytical results for many constituents (particularly pharmaceuticals, wastewater-indicator compounds, and trace elements) were commonly below laboratory reporting limits (LRLs), which inherently limits data analysis [16]. In subsequent interpretive analyses, all censored values below the LRL were imputed, using regression-on-order (ROS) statistical methods [48], implemented in the NADA package of R [49]. The ROS method allows for robust estimation of censored values by using the structure of the dataset above the LRL and projecting this data structure to represent the censored concentrations. Results of imputation by ROS were further filtered by removing all parameters for which greater than 80% of the results were below the LRL (i.e., when more than 80% of the dataset was censored for a given parameter). This was done to limit the multivariate evaluation to those parameters with more distributed spatial information and to limit the effects of the ROS imputation on the overall results of the study.
Principal component analysis (PCA) was applied on the modified dataset to reduce the dimensionality of the dataset and aid in interpretation of geochemical signatures of different water types [50]. Geochemical data types were grouped together as follows: major ions, trace elements, REEs, noble gases, noble-gas isotopes, and pharmaceutical and wastewater-indicator compounds. Reference [51] suggests that multivariate compositional data analysis (CoDA) may not be useful for isotopes, so isotope ratios of noble gases were included in the PCA to test that assumption. The group-specific dataset was scaled, using the centered-log ratio transform (clr), which allows closed compositions (i.e., compositions summing to a fixed value [52]) to be reliably evaluated, using statistical methods [50]. There are a variety of transformations that could be applied (isometric-log ratio, ilr; allometric-log ratio, alr; centered-log ratio, clr), and Reference [52] suggests that the ilr may be superior for multivariate evaluation. However, application of the results of ilr with PCA may be difficult because back-transformation is required prior to interpretation [9,50]. As an alternative the clr was used, as it has been established as a reasonable method for scaling prior to PCA [9,51]. PCA was then applied, using the prcomp function of R [49].
Iterative PCA using different combinations of parameters may be helpful to further define groupings and understand processes [9]. Therefore, in addition to PCA of the previously defined geochemical groups, a combined PCA was also completed, using a subset of parameters. These parameters were included based on previous group-specific PCA, prior geochemical knowledge of the system, and other research in the Fountain Creek alluvial aquifer [17,22]. All parameters were converted to the same units of µg/L prior to regrouping because this is required by PCA [50]. Following unit conversions, the dataset was again processed using the clr prior to completion of PCA.

2.3. Noble Gases and Groundwater Residence-Time Tracers

2.3.1. Dissolved Gas Modeling

Estimation of conditions during groundwater recharge (specifically water temperature, atmospheric pressure, excess air, and salinity) are required for attributing proportions of dissolved gases 3H, He, and SF6, which are used in lumped parameter models [53] and apparent-age calculations [54,55]. Inversion of dissolved-gas concentrations are commonly used to estimate recharge conditions, and the inclusion of noble gases in addition to major gases allows for greater confidence because of the specific solubility characteristics of noble gases [35,56].
In this analysis, Ne, Ar, Kr, Xe, O2, and CO2 were used because they are primarily of atmospheric origin. Helium and N2 may have non-atmospheric sources [8,12] and were only used when no viable models were obtained from the set of gases with primarily atmospheric origin. Helium and N2 were used in six of the 19 models for sampled groundwater wells where noble-gas data were available. Although Ar can have terrigenic sources, use in gas modeling is reasonable if samples do not show terrigenic 40Ar production [35]. Evaluation of the ratio of 36Ar/40Ar in samples versus that in the air (36Arair/40Arair = 0.003378) [35] indicates that subsurface radiogenic production of 40Ar can be ignored (i.e., mean 36Arsample/40Arsample / 36Arair/40Arair = 0.998).
For the nine groundwater wells that were not sampled for noble gases, recharge temperature and excess air were computed by the USGS Groundwater Dating Laboratory from concentrations of O2, N2 (corrected for excess gas), Ar, and CO2 according to methods described in References [57,58] and References [59,60].
Modeling with noble-gas data utilized an inverse approach similar to those described by References [61,62]. Input parameters included the recharge salinity from groundwater sampling and the elevation of the land surface in feet above the North American Vertical Datum of 1988 (NAVD88) at the well. Using these inputs, an inverse-error fitting procedure was applied wherein sets of simulated gas concentrations are produced and compared to the observed values, for each of the three most commonly used models of excess air: closed-system equilibration (CE), partial re-equilibration (PR), and unfractionated excess air (UA), as described by Reference [35]. The error in these sets of simulations is quantified by using the χ2 parameter, which accounts for model errors in each gas. The model returns a recharge temperature, excess air concentration, χ2 value, and probability for each of the CE, PR, and UA models. Recharge temperature and excess air are crucial for groundwater-age estimates discussed in the next section, and in of themselves, are useful for hydrologic characterization [63].
A sensitivity analysis was conducted for several wells to evaluate the influence of changing recharge elevation on calculated recharge characteristics. Three different situations were simulated: (1) recharge elevation set equal to wellhead elevation (minimum elevation possible; this is the base-case model described previously), (2) recharge elevation set equal to the median elevation within the watershed, and (3) recharge elevation set equal to the maximum possible elevation within the upgradient watershed. Results of the sensitivity analysis indicate that simulated recharge conditions are relatively insensitive to the input recharge elevation (see results in Supplementary Materials).

2.3.2. Apparent Ages and Lumped-Parameter Modeling

Concentrations of 3H, 3He, and SF6 in groundwater were used as tracers to evaluate the presence of modern water [55,64], as well as to estimate residence-time distributions [65]. Tracer concentrations were normalized to the groundwater-recharge conditions (temperature and excess air) calculated from dissolved gases as described previously.
A full discussion of the complexities and nomenclature relating to estimation of the time since groundwater recharge has occurred (generally referred to as “groundwater age”) is not discussed herein but is provided by References [66,67]. Generally, groundwater age is conceptualized as having a distribution rather than a discrete value because of processes such as mixing and dispersion. Methods were applied in this work to estimate age distributions of modern groundwater (i.e., groundwater recharged since the 1950s). For consistency, the terminology of Reference [67] is used with respect to apparent tritium-helium age (3H/3He), idealized age (for SF6), and mean-residence time (for lumped-parameter models). Apparent 3H/3He ages were calculated according to methods described in References [54,56]. The apparent 3H/3He age relies on the known decay rate of 3H (derived from atmospheric nuclear weapons testing) to 3He, combined with an understanding of He sources in the aquifer [54]. Piston-flow ages (i.e., idealized ages) were calculated from SF6 concentrations in water samples compared to SF6 concentrations in the atmosphere through time based on methods described in Reference [45]. Finally, residence-time distributions were evaluated, using the program TracerLPM [53] for the piston-flow model (PFM), exponential-mixing model (EMM), and dispersion model (DM). These models are consistent with simplified understanding of the groundwater-flow system based on aquifer thickness, groundwater-flow directions, and residence-time distribution functions [68].
In addition to quantitative estimates of groundwater age, the fractions of old and young groundwater can be calculated, using 3H data alone (herein referred to as the 3H thresholds method), which may be less ambiguous than application of lumped parameter models (LPMs) because the conceptual mixing relations need not be known [69]. The 3H thresholds method [64,69] was applied to calculate the fraction of water recharged since the advent of atmospheric nuclear testing in 1953 (Fpost-1953), using the measured 3H concentrations and estimated 3H in precipitation [70,71,72]. The 3H threshold method allows three categories of groundwater to be differentiated; modern (recharged after the 1950s), pre-modern (recharged prior to the 1950s), or mixed. The method consists of comparing decayed concentrations of 3H in precipitation for the study area to the 3H concentrations in water samples. Based on different thresholds that are established, the sample is then categorized into one of the three groups. Differences between the precipitation input datasets of References [70,71] are assessed in Supplementary Material Section S4.

3. Results and Discussion

All data described in the following sections are available through the USGS NWIS database [24]. Results are compiled according to USGS site numbers and site names, provided in Supplementary Materials Table S1. Groundwater wells were categorized into three groups based on their spatial location (Figure 1). Tributary alluvium wells are in the paleochannels to the northeast and upgradient from the main aquifer, wells proximal to the creek are within approximately 200 meters of Fountain Creek, and the remainder of the wells were classified as central groundwater.

3.1. Major Ions and Trace Elements

The major-ion composition of groundwater and surface water covers a broad range but also shows some spatial associations. Based on the Piper diagram [73] for the study area surface water is typically Na+K-Cl type, whereas groundwater is typically Ca-SO4 type or Ca-HCO3 type (Supplementary Materials Figure S1). Fountain Creek and Canal 4 have appreciable Cl possibly indicative of road salting effects. The composition of central groundwater (Supplementary Materials Figure S1) is between tributary alluvium and surface water, indicating a mixed origin. Ratios of Cl/Br and concentrations of Cl in groundwater may be used to indicate mixing and inputs from anthropogenic sources [74,75,76]. An analysis of Cl/Br ratios and Cl concentrations (Supplementary Materials Figure S2) indicates that groundwater well 04-009, located in the tributary alluvium (Figure 1), is representative of native groundwater recharge and that the majority of groundwater at other wells can be explained by variable mixing between native groundwater recharge with WWTP effluent, septic tank discharge, and road salting. Similarly, surface waters have compositions that include WWTP effluent, road salting, and meteoric precipitation. This analysis indicates that mixing relations between surface water, groundwater, and anthropogenic sources are important for understanding solute distributions in the aquifer. These results are further corroborated and strengthened by REE data and PCA, as described below.
Oxidation-reduction (redox) classification of the waters, using the method of Reference [77] indicates that 26 of 31 groundwater sites are characterized by oxic conditions whereas the remaining five display mixed oxic-anoxic conditions (Supplementary Materials Table S2). Mixed conditions of this type are possibly consistent with NO3 reduction or Mn reduction [77].
Trace-element (Al, Cu, Fe, Mn, Pb, and Zn) data show spatial patterns in groundwater and surface water (Supplementary Materials Figure S5) potentially linked to variable processes or sources controlling transport, though Al may be marginally affected by sample bias and is thus not interpreted further. For instance, Cu, Fe, Mn, and Zn tended to have greater concentrations in surface water than in groundwater (Supplementary Materials Figure S5), possibly attributed to runoff from of mine tailings in the area as documented by Reference [22]. In contrast, Pb tended to have greater concentrations in groundwater potentially linked to a subsurface source or redox conditions in the aquifer.

3.2. Stable Isotopes

The δ2H and δ18OH2O composition of groundwater (Figure 2) is generally consistent with mixing of water sources from a snowpack-driven Rocky Mountain meteoric water line (RMWL-Snow) as described by Reference [78] and annual variation in precipitation for the site estimated by using the online isotopes in precipitation calculator (OIPC) [79]. The composition is shifted slightly to the right of the snow-dominated sites likely because the area receives both snowpack accumulation and recharge from summer rainstorms (as reflected in the relative positions of the data and the OIPC meteoric water line (MWL). Tributary alluvium groundwater wells have a large range in isotopic composition, covering nearly the entire spread in the dataset (Figure 2a). This range is most likely related to seasonality of recharge, though the datasets required to thoroughly analyze such trends are not available for the study site (e.g., see Reference [80]). There are no clear spatial associations of wells with more depleted versus enriched compositions. Surface-water samples have a slightly shallower slope (slope = 7.16) versus groundwater (slope = 8.19), consistent with some evaporation of surface water [81].
Stable isotopes of nitrate (δ15N, δ18ONO3) provide evidence that denitrification is occurring in the aquifer, as several wells plot along a line sloping at 1:2 (δ15N:δ18ONO3) in Figure 2c, a trend typical of denitrification [41,82]. These wells also show ammonia and nitrate + nitrite (NO3 + NO2) concentrations consistent with isotopic fractionation by microbial reduction (Figure 3). Several groundwater wells and streams plot in a region with generally elevated NO3 + NO2 and ammonia (circled sites on Figure 3b). These wells are in the general proximity to current or former WWTPs (Figure 1), and thus the coincidence of multiple elevated nitrogen species is likely linked to WWTP effluent [22]. It is interesting to note that the indication of some nitrate-reducing conditions by N isotopes and concentrations of NO3 + NO2 and ammonia is slightly inconsistent with primarily oxic conditions as indicated by the method of Reference [77]. Specifically, of the wells in Figure 3 noted to have stable-isotopic evidence of denitrification, only well U-7 was noted to have mixed oxic/anoxic redox conditions by the method of Reference [77]. Isotopes of N and O are ideally suited to identifying nitrate reduction [41,83] whereas the idealized method of redox-condition assignment is subject to a variety of limitations [77]. The prevalence of oxidizing conditions, as indicated in Supplementary Materials Table S2, is believed to be generally representative of the system, whereas results in Figure 3 indicate that nitrate-reducing conditions may exist in some locations of the aquifer. An alternative is that the observed isotopic fractionation occurred prior to the water entering the hydrologic system (such as within the WWTP). This seems somewhat unlikely given that denitrification is indicated in all three zones of the aquifer (Figure 2c) meaning that water from similarly denitrified sources would have had to migrate throughout the aquifer and enter via recharge zones that are not near one another.

3.3. Rare Earth Elements

Rare earth elements are useful indicators of redox processes [84], groundwater mixing [85], water-rock reaction [86], and for tracing anthropogenic water sources [87]. Groups of elements amongst the REEs are also useful for tracing geochemical processes, specifically gadolinium (Gd) may be indicative of anthropogenic sources because it is commonly used in magnetic resonance imaging [7,88,89], and cerium (Ce) and europium (Eu) are redox sensitive and more likely to be derived from rock weathering [86,90]. In this study, three REE anomalies are used: Gd anomaly calculated according to the method of Reference [89]; Ce anomaly and Eu anomaly, both calculated according to the method of Reference [90]. The calculation of each anomaly has a specific purpose, with the Gd anomaly being indicative of the influence of WWTP effluent [87,89] and anomalies of Ce and Eu being indicative of redox potential and geochemical reaction between water and the aquifer matrix, respectively [86,90].
The Gd anomaly and visual inspection of REE concentrations normalized to the North American Shale Composite (NASC) [91] (Figure 4) indicate that surface water is affected by WWTP effluent. The NASC is a reasonable proxy for use in this instance because the Pierre Shale forms the basal surface of the alluvial aquifer. Gadolinium anomalies (unitless) in surface water range from 0.99 to 124 and are evident in the sharp pattern of increased NASC-normalized Gd concentrations relative to nearby elements in Figure 4, specifically in groundwater wells proximal to Fountain Creek (Figure 4a) and in Fountain Creek (Figure 4c). Note however that the calculated Gd anomaly values are not displayed on Figure 4, only the visual pattern arising from Gd enrichment. For comparison, Gd anomalies of greater than 1.5 are typically indicative of WWTP inputs [87], and Reference [89] noted Gd anomalies up to nearly 140 in WWTP effluent. The range of Gd anomalies in groundwater and surface water in this study indicates that there are a range of WWTP effects, from non-affected (Gd anomaly less than 1.5 in a variety of groundwater and surface-water locations) to highly-affected (Gd anomaly up to 124 in Canals 4 and 46 in Fountain Creek). Groundwater wells in proximity to creeks also display a subdued Gd anomaly (Figure 4a), whereas groundwater in the tributary alluvium and central groundwater generally do not (Figure 4b,d). The Gd anomaly in wells proximal to the creeks indicates that streamflow losses likely contribute to recharge in these areas, whereas distal groundwater from surface-water bodies is not influenced by streamflow losses from the creek or WWTP effluent.
The Ce anomalies range from 0.02 to 9.9 and are substantially greater in surface water than groundwater (Supplementary Materials Figure S3). Lesser Ce anomalies in groundwater may be caused by precipitation of dissolved Ce upon transformation of Ce3+ to Ce4+ [86,90], which is in agreement with the generally oxidizing character of groundwater in the study area. Greater Ce anomalies in surface water are somewhat unexpected given the assumed oxic nature of creeks in the area (which would presumably scavenge Ce from the aqueous phase via adsorption and decrease Ce anomaly). Detailed research on REEs in surface waters has indicated that a complex set of in-stream processes may control Ce anomalies [84], and the dataset in this study does not allow for mechanistic understanding of these in-stream processes such as adsorption, mineral precipitation, or speciation. The Eu anomalies have a smaller range from 0.007 to 2.11 and display the opposite pattern from Ce anomalies in that groundwater wells have greater Eu anomalies than surface-water sites (Supplementary Materials Figure S3). Reference [90] illustrated that Eu anomalies were poorly correlated with redox potential, and Reference [86] found that Eu anomalies were useful as an indicator of geochemical reactions within the aquifer matrix. Given these previous research findings, the pattern of Eu anomalies in the Fountain Creek aquifer is interpreted as indicative of weathering reactions within the aquifer, whereas the Eu anomaly of surface water is likely imprinted from anthropogenic sources (as corroborated by Gd anomaly). The imprint of human activity on Eu anomaly is also evident given that NASC normalized Gd concentrations are used in calculating Eu anomaly, meaning these Gd and Eu anomalies are not entirely independent.
Although general REE patterns and anomalies of Gd, Ce, and Eu are helpful for evaluating the interconnectivity of groundwater and surface water and the presence of WWTP effluent [87], the visual relations in NASC normalized diagrams do not allow for quantitative mixing or other processes to be evaluated. The REE normalization plot [85,90,92], however, may be useful for such analyses (Figure 4e). This plot helps to further contextualize the results graphically displayed in NASC normalized diagrams. For this analysis, the following classifications of REEs were used; light rare earth elements (LREEs) = La, Ce, Pr, and Nd; middle rare earth elements (MREEs) = Sm, Eu, Gd, Tb, and Dy; and heavy rare earth elements (HREEs) = Ho, Er, Tm, Yb, and Lu. See References [85,90,92] for further information on the construction of such plots.
Results shown in the normalized REE plot indicate that surface water of Fountain Creek and Canal 4 with pronounced Gd anomalies tend to plot near the bottom center, whereas other surface water does not plot in a tight grouping. Groundwater plots in a wide range of compositions with tributary alluvium groundwater tending to plot closer to the MREE = LREE line and other groundwater plotting farther to the right. To illustrate the utility of normalized REE plots in understanding mixing of water sources and geochemical reactions, a synthetic mixing line was calculated between two locations with near end-member compositions of REEs, T02-MW006 (tributary alluvium groundwater) and Canal 4 at Headgate (surface water with apparent anthropogenic effects). These water-quality sampling locations were chosen based on their position in space on the normalized REE plot in Figure 4e, which is a distinct mixing analysis from that presented earlier for Cl/Br ratios (Supplementary Materials Figure S2). Points along the synthetic mixing line are labelled with the fraction of groundwater in the mixture (x). The mixing line is not linear but takes an arcuate path, and mixing fractions of equal interval are not plotted at equal intervals along the mixing line.
The purpose of the synthetic mixing line is to show that central groundwater falls outside of the zone that may be explained by simple conservative mixing between water types with distinctive REE chemistry (REE end members). This is somewhat unexpected considering central groundwater wells are located downgradient from the tributary alluvium sites and Canal 4 (Figure 1). There are two possible explanations for this: (1) A third unsampled end-member water composition exists that contributes to the chemistry of the central groundwater, or (2) water–rock interaction produces a shift in the REE composition. The second option appears likely given that the compositions of central groundwater are becoming more enriched in terms of REE concentrations compared to the NASC, which was noted by Reference [85] to indicate geochemical evolution. This analysis provides evidence that groundwater in the Fountain Creek alluvial aquifer is evolving geochemically along the flow path from recharge to discharge zones and that zones of focused groundwater recharge exist. Specifically, based on Gd anomaly in U-7 (Figure 4a), it appears that Fountain Creek loses flow to groundwater recharge in the northern part of the study area in the reach downstream from the gage at Fountain Creek below Janitell Road below Colorado Springs, Colorado (USGS site number 07105530), and that the tributary alluvium was recharged by flow lost from discharge of the WWTP on the East Fork Sand Creek (based on Gd anomaly of 1.75 in MW1-1). This REE mixing analysis is complementary to the mixing analysis of Cl/Br presented previously in that both analyses indicate that groundwater in the tributary alluvium and proximal to Fountain Creek have end-member compositions, whereas central groundwater has a composition that is a mixture between the two end members with the addition of external solutes (both from anthropogenic and natural sources).

3.4. Pharmaceutical Compounds and Wastewater-Indicator Compounds

Previous studies on CECs have commonly used detection frequency as opposed to absolute concentrations [4], and that approach is used here, along with spatial analysis of concentrations. The pharmaceutical and wastewater-indicator compounds most commonly detected from this study are listed in Table 1. Typical uses and sources of many of these compounds are summarized in References [2,93,94].
Several compounds show differing abundance in groundwater versus surface water whereas others are found in similar distributions of samples. Compounds that were more prominent in surface water than groundwater are meprobamate, gabapentin, tramadol, metformin, bupropion, sulfamethoxazole, hexamethylenetetramine, and lidocaine. Other compounds in Table 1 generally had similar incidences of detection between groundwater and surface water. The overall high incidence of detection of these compounds in both groundwater and surface water show that the Fountain Creek alluvial aquifer is affected by WWTP effluent, consistent with REE analysis.
Many of the most commonly detected compounds in the Fountain Creek alluvial aquifer have been noted to occur in other settings [2,4]. An exception is hexamethylenetetramine, which was not noted to occur in other studies. One blank sample had a value greater than the LRL for hexamethylenetetramine of 36 ng/L. Compared to the average environmental concentration of 720 ng/L, this small potential bias of 5% is considered unlikely to affect results and interpretation. Hexamethylenetetramine may be derived from agricultural operations and can be associated with some pesticides [95] or from a variety of pharmaceutical uses [96]. Hexamethylenetetramine is believed to be rapidly leached from soils to groundwater and to be both abiotically and biotically degraded in the environment [95]. The greater prevalence of hexamethylenetetramine in surface water compared to groundwater (Table 1) supports the assumption that rapid degradation of the compound is occurring in groundwater.
Maps of CEC compounds in groundwater and surface water (Supplementary Materials Figure S6) support the interconnected nature of the system. Discharge points of CECs into surface waters are the WWTPs shown as potential solute sources on Figure 1, both adjacent to Fountain Creek in the northern portion of the study reach and on the east fork of Sand Creek near groundwater monitoring well T02-MW006. Canal 4 is a conveyance for water partially sourced from WWTPs and is believed to lose flow to the groundwater system based on field observations. The spatial distribution of pharmaceutical compounds indicates that solutes may be added to the aquifer from seepage from Canal 4 (Supplementary Materials Figure S6). Based on the spatial distribution of CEC concentrations (Supplementary Materials Figure S6), different CECs may have different sources and controlling processes in the hydrologic system. For example, methylbenzotriazole and carbamazepine are relatively widely distributed in groundwater and are found nearly ubiquitously in surface waters, whereas other compounds, such as bupropion and gabapentin, are sparingly observed in groundwater but are abundant in surface water. The variable occurrences are likely linked to different sources, transport properties, and degradation of some compounds in the environment [2].

3.5. Principal Component Analysis

Several iterations of PCA were completed [9], first on groups of parameters (Supplementary Materials Figures S7–S9), then on a subset determined to be most useful for discriminating between water sources and geochemical reactions (Figure 5). The final subset of parameters included the following: Hexamethylenetetramine, methylbenzotriazole, carbamazepine, and fluconazole were included because of their high frequencies of detection and contrasting behavior between groundwater and surface water (see Section 3.4); U and Se were included because they are known geologically derived constituents in Fountain Creek [17]; Gd and Eu were included because Gd is associated with WWTP effluent and streams, and groundwater displays spatially variable Eu-anomaly (see Section 3.3 and Supplementary Materials Figure S3); Cd and Mn were included because of previous investigations showing point sources in the study area [22] and trace-element-specific PCA (Supplementary Materials Figure S8); Rb was included because of trace-element-specific PCA (Supplementary Materials Figure S8) and because it may be indicative of weathering reactions in the aquifer matrix [86]; and finally NO3 + NO2 was included because it is primarily of anthropogenic (WWTP) origin and because of the potential influence of denitrification in the aquifer.
Results of PCA are shown in Figure 5. Rays (black lines) correspond to each parameter, and the lengths of the rays are proportional to the parameter variability within the two plotted components PC1 and PC2; the angles between any two rays is a measure of the correlation between the two variables [50]. The position of the symbols relating to each water-quality sample is a function of the score of that sample against the two PCs. Rays for constituents with a clear anthropogenic source tend to plot in the lower left quadrant with rays of nearly equal length. Constituents of known geologic origin U and Se [17] plot nearly on top of one another and are located near the NO3 + NO2 ray on the right side of the plot. Manganese, Cd, and Rb plot in the upper quadrant whereas Eu plots close to the CEC rays in the lower left, but with a shorter-length ray.
The location of samples in component space supports previous analyses in the identification of distinct end members and mixing. Most surface waters (especially those affected by WWTP effluent) plot between the Cd ray and the CEC rays. Three of the five groundwater wells in proximity to the creek also plot in this region. Reference [22] noted the presence of some trace metals including Cd and Mn in stormwater stream samples, and this analysis indicates that these metals along with pharmaceutical and wastewater- indicator compounds are a tracer of surface waters and are indicative of surface-water infiltration where they are found in groundwater samples. Tributary alluvium groundwater all plot on the right half of the plot in the arc between the Rb and fluconazole rays, with the exception of T01-MW002. This well is located downgradient from the former WWTP on Sand Creek (Figure 1) and appears affected by denitrification and nitrate reduction (Figure 3) such that it plots opposite from the NO3 + NO2 ray (because the well is devoid of NO3 + NO2 and thus this parameter has little control on the composition of the well) and closer to the rays representative of WWTP effluent. Aside from T01-MW002, PCA for all other tributary alluvium groundwater indicates correlation with geologically derived constituents Rb, Se, and U. Central groundwater tend to plot in an intermediate zone between tributary alluvium groundwater and groundwater proximal to the creek and surface water. This is consistent both with mixing of water sources (as indicated previously by Cl/Br and REE analyses) and with the occurrence of anthropogenically and geologically derived constituents in the central zone of the aquifer.

3.6. Noble Gases and Groundwater Residence-Time Tracers

3.6.1. Noble-Gas Isotopes

Isotopes of noble gases can be informative as to residence times and mixing of distinct waters [12,13,35,97]. Helium isotopes (Figure 6a) indicate that groundwater proximal to the creek and in the tributary alluvium are generally close to equilibrium with air-saturated water (ASW) for the study area (calculated using the method of Reference [98] and assuming 14 °C and 1700 meters elevation NAVD88) based on ratios near 1 of 3He/4He in the sample to 3He/4He in air (R/Ra). Several samples in the tributary alluvium have R/Ra values greater than 1, indicating 3He production. Given the abundance of tritium in all groundwater samples (discussed below), the R/Ra values greater than 1 are likely attributable to tritiogenic 3He (3Hetrit; [54]). Samples in the central groundwater zone display lesser R/Ra values (0.75 to 0.90; Figure 6a) and greater ratios of 4He/20Ne. These patterns are consistent with production of 4He in the subsurface relative to atmospherically derived 20Ne likely caused by increasing residence time and interaction with aquifer materials [13]. A notable exception to the pattern of lesser R/Ra values in central groundwater wells is NONUM-2, which although centrally located in the aquifer, is also downgradient from Johnson Reservoir (Figure 1). Noble-gas isotopes indicate that this water is not geochemically evolved and may be sourced from the reservoir.
Possible recharge from the mountain block to the west is indicated by the calculated noble-gas recharge temperature (NGRT) of 4.3 °C for well TH-5, substantially colder than other wells in the aquifer (Supplementary Materials Table S4). Well TH-5 is located near Fountain Creek and is aligned with a paleochannel sourced from higher elevations along the mountain front on the western part of the study area (Figure 1). The presence of abnormally cold NGRTs has been used in the past to recognize mountain-block recharge [62,63,97]. In this instance, the alluvial aquifer may have a recharge source to the west, which would have gone undetected without the application of noble-gas modeling.
Calculation of apparent 3H/3He ages and application of lumped-parameter models are subject to uncertainty regarding the ratio of 3He/4He in water that is no longer in equilibrium with the atmosphere [12]. To evaluate subsurface He production, the method of Reference [99] to derive site-specific ratios for the addition of terrigenic He components was applied (Figure 6b). These results corroborate those of 4He/20Ne in indicating geochemical evolution in the aquifer. A full sensitivity analysis focused on subsurface He inputs is included in Supplementary Materials Section S5 of the (Supplementary Materials Table S3).
Interaction with aquifer materials is also illustrated by Xe isotopes (132Xe/20Ne; Figure 6c). Several wells near Fountain Creek (specifically U-7, A1, and TH-5) have 132Xe concentrations above that explainable by equilibrium with ASW, consistent with either interaction of water with the underlying Pierre Shale, or inflow of groundwater sourced from the shale [13]. These wells are near a reach of the creek known to have variable groundwater–surface-water interactions [20], and based on the Xe enrichments in groundwater wells in this area, it is considered plausible that there is some groundwater circulation within the Pierre Shale at the base of the alluvial aquifer. Upward leakage from the shale (which is typically considered a confining layer) could be attributable to sustained groundwater pumping in the area, which is known to cause leakage from confining layers [100].
Figure 6. Noble-gas isotope compositions, namely (a) R/Ra versus 4He/20Ne (Ra = 1.4 × 10−6), (b) helium evolution diagram (e.g., see Reference [99]), and (c) 132Xe/20Ne versus 4He/20Ne. The helium evolution diagram shown in (b) illustrates the ratio of 3He/4He in samples corrected for air-bubble entrainment (Rc) versus 4He derived from solubility (4Heeq) over total corrected 4He (4Hec). Wells of interest are labeled with site alias, for discussion. The composition of air-saturated water (ASW) was calculated, using the method of Reference [98] and the atmospheric composition from Reference [101], solubility at 14 degrees Celsius (°C), and 1700 m above North American Vertical Datum of 1988 (NAVD88). Extents of analytical uncertainties are within the points on all plots.
Figure 6. Noble-gas isotope compositions, namely (a) R/Ra versus 4He/20Ne (Ra = 1.4 × 10−6), (b) helium evolution diagram (e.g., see Reference [99]), and (c) 132Xe/20Ne versus 4He/20Ne. The helium evolution diagram shown in (b) illustrates the ratio of 3He/4He in samples corrected for air-bubble entrainment (Rc) versus 4He derived from solubility (4Heeq) over total corrected 4He (4Hec). Wells of interest are labeled with site alias, for discussion. The composition of air-saturated water (ASW) was calculated, using the method of Reference [98] and the atmospheric composition from Reference [101], solubility at 14 degrees Celsius (°C), and 1700 m above North American Vertical Datum of 1988 (NAVD88). Extents of analytical uncertainties are within the points on all plots.
Water 13 00871 g006

3.6.2. Apparent Ages and Presence of Atmospheric Tracer Gases

Apparent 3H/3He ages (e.g., see Reference [67]) were calculated for all wells in the dataset with applicable data and assuming He isotope-abundances derived from terrigenic interaction (3He/4He = 2 × 10−8). In groundwater with appreciable terrigenic He (as indicated by noble-gas isotopes), the 3H/3He apparent-age calculations can be biased [102]. Moreover, some of the wells are likely subject to denitrification (Figure 3), which could cause degassing and stripping of noble gases, 3H, 3He, and SF6 in unequal quantities proportional to their solubility [103], which would also bias ages in a variety of ways. Therefore, these results should be viewed with discretion and compared with complementary datasets as discussed below. A summary of the results of all applicable methods of estimating groundwater-residence times is presented in Table 2. Additional details related to environmental-tracer modeling are summarized in Reference [104].
All wells contain measurable 3H, indicating that some recharge has occurred since the 1950s and the advent of atmospheric nuclear testing [54]. All but one well have 3H concentrations in the range of 4–10 tritium units (TU). The CAD1 well, however, has a substantially greater value of nearly 20 TU, indicating that this well may have recharged during the bomb pulse of atmospheric nuclear testing in the 1960s [72]. The apparent 3H/3He age for this well is only 3.2 years, thus indicating the complexity introduced into these calculations by uncertainty in terrigenic He-isotope systematics [99,102]. These isotope systematics may be influenced by upward leakage from the Pierre Shale, which could impart a unique He signature to waters in the vicinity of high-volume pumping wells. Apparent 3H/3He ages for wells with applicable data range from about 0 to 21 years, although a majority of wells have 3H/3He ages of less than 10 years. These results have less precision from wells with noble-gas isotopes indicative of longer residence times (Figure 6) because of uncertainties arising from input of terrigenic He with an unknown 3He/4He ratio. The JCC well has the oldest apparent 3H/3He age of 21.5 years and a noble-gas composition indicative of interaction with the aquifer matrix (Figure 6a). Although it can be stated with some certainty that this well is in contact with older water in the aquifer, the exact age is uncertain, and true advective age for a groundwater sample does not exist because groundwater samples are mixtures of water molecules arising from converging flow paths [67].
Although all wells clearly contain a component of modern water based on the presence of 3H, this component has the ability to mask older water fractions by producing a young bias on the apparent ages [65,67]. One way of evaluating this age bias is through the calculation of the fraction of water recharged since 1953 [69]. Mixing-fraction calculations of this type indicate a range of recharge after 1953 (Fpost-1953) of 0.36 to greater than 1 (Table 2). Theoretically, the Fpost-1953 values should have a maximum at 1, but values greater than 1 may be obtained depending on the times series of 3H in precipitation. This study used the function from Reference [70], as detailed in Supplementary Materials Section S4. The well with the lowest Fpost-1953 (and therefore the highest fraction of water recharged prior to 1953) is 04-009, which is located at the head of a paleochannel in the central portion of the study area (Figure 1) and has a Cl/Br composition representative of meteoric precipitation (Supplementary Materials Figure S2). The high proportion of older water in this well combined with the juvenile composition may indicate that well 04-009 is representative of pre-anthropogenic hydrologic conditions for the study area.
Concentrations of SF6 (both in terms of molality and in equilibrium with recharge conditions in the study area in parts per trillion by volume, pptv) vary widely. Idealized SF6 piston-flow ages range from about 2 to 15 years (Table 2). Comparison of these idealized ages with 3H/3He apparent ages shows that the idealized SF6 piston-flow age tends to be greater, except in the JCC well. Several wells contain concentrations that are several orders of magnitude greater than that explainable by simple atmospheric equilibrium with recharge parameters (temperature, excess air) derived from noble-gas modeling. This phenomenon has been noted to occur both in specific aquifer lithologies [105,106] and in spatial association with urban areas [107]. Given that sediments in the alluvial aquifer are partially derived from weathering of granitic rock, it is considered likely that local contamination occurred from terrigenic SF6 due to U-series decay [108]. However, local variations are extreme, and additional analysis is warranted to separate terrigenic from potential other sources.

3.6.3. Lumped Parameter Models

Results of lumped parameter models (LPMs) can indicate the mean residence time of groundwater [65,67]. The piston-flow model (PFM), exponential-mixing model (EMM), and dispersion model (DM) were applied in this study. Selected results according to different LPMs are illustrated in Figure 7. Additional details related to environmental-tracer modeling are summarized in Reference [104].
Comparison of tracer-tracer plots indicate variable agreement with simplified LPMs. A number of wells in the central groundwater zone and several in the tributary alluvium are consistent with both the EMM or the DM and have a range of mean residence times from 5 to 10 years (CO259-25, CO259-26, NONUM-2, T13-MW004, TH-22, U-12, U-14B). These mean residence times are generally in agreement with 3H/3He apparent ages and idealized SF6 piston-flow ages (Table 2). Data plot closer to the LPM lines when considering 3Hetrit versus 3H (Figure 7a) as opposed to 3Hetrit versus SF6 (Figure 7b) or SF6 versus 3H (Figure 7c). This is because of likely terrigenic and external SF6 sources, which have caused excess SF6 to be present in the system.
Several wells are consistent with a PFM, namely A1, MW2-4, and TH-5, with mean residence times of generally less than five years. Two of these wells are located proximal to Fountain Creek (A1, TH-5) and have Gd anomaly consistent with infiltration of surface water (Figure 4). The short mean residence times, agreement with the PFM, and presence of surface-water REE signatures indicate that groundwater in the immediate vicinity of Fountain Creek receive a majority of their recharge from creek losses and that the flow paths from the creek to the well are simple enough such that they approximate piston flow.
Several wells cannot be explained by any of the applied LPMs, specifically CAD1 and U-7. Further analysis illustrates several possible reasons for these wells not fitting LPMs. For instance, U-7 has 3H contents consistent with recent recharge, no 3Hetrit, and excess SF6. The water from this well also has one of the largest Gd anomalies of any groundwater well (Figure 4). A possible explanation for the lack of agreement with any LPM is that this well was very recently recharged by water lost from Fountain Creek (which is in agreement with the REE chemistry and the 3H/3He apparent age of zero years, Table 2). The excess SF6 is more enigmatic, but one possible explanation is that a more complex model for the formation of excess air beneath the losing stream is required (i.e., aside from the CE, PR, and UA models considered) [35]. Reference [109] found that reasonable results of dissolved-gas modeling near areas of focused stream losses were not always properly linked to physical processes in the aquifer, and that more detailed analyses may be required under these conditions. Considering CAD1 (which has non-zero 3Hetrit and bomb pulse 3H), the inconsistency may be caused by differential diffusion within the aquifer. Diffusion of 3Hetrit is more rapid than 3H [103], and thus it is likely that some 3Hetrit has been lost via diffusion during groundwater flow, offsetting the point in Figure 7a for CAD1 down and to the right, off the DM line. An alternative explanation for lack of fit for U-7 and CAD1 to LPMs is that binary mixing between old and young waters of substantially different ages is occurring [11]. Although this may be possible, it cannot be explored with the dataset given than no tracers of “old” groundwater (e.g., 14C) were collected as part of this study. Taken as a whole, the LPM results affirm the findings of previous research that non-uniqueness and uncertainty may be limited by applying multiple tracers [56,67,110].
A comparison of the three methods of estimating groundwater age (idealized SF6 piston-flow age, 3H/3He apparent age, and mean residence time from LPMs) is shown in Figure 8, which allows age bias within the dataset to be explored and groundwater-residence times to be visualized. Based on Figure 8, 3H/3He apparent age may be biased low, as a majority of wells have apparent ages less than five years old. This may be caused by somewhat complicated He mass balance in the aquifer (see Figure 6), introducing uncertainty in the amount of 3Hetrit used in the apparent-age calculations. Idealized SF6 piston-flow ages tend to be greater than 3H/3He apparent age but span a smaller range. Mean residence time estimates from TracerLPM have a greater number of possible analyses (because TracerLPM may produce multiple age estimates for a single sample based on combinations of different tracers), and these residence-time estimates span a much greater range than either other method (yielding a maximum mean residence time of 62 years). It is important to consider that the range of mean residence-time estimates from TracerLPM is not a measure of the full age distribution of each sample (which could include waters much older than illustrated); rather, the estimates are the distribution of possible mean ages for the entire dataset based on the tracers used in this study.

3.7. Conceptual Hydrologic Model

Combining the insights gained by the array of tracers used in this study allows for a general conceptual model of groundwater flow and solute sourcing and transport (Figure 9) that can serve to guide future studies and modeling. The influence of anthropogenic solutes is clear throughout surface waters in the study area based on a Gd anomaly and the occurrence of pharmaceuticals and wastewater-indicator compounds. Active groundwater recharge occurs both in the tributary alluvium and near Fountain Creek, although recharge in these two locations is derived from different sources. Focused recharge from streamflow loss occurs near Fountain Creek and is recognized by a Gd anomaly in both surface water and groundwater, the distribution of pharmaceuticals and wastewater-indicator compounds, and the agreement of near-creek groundwater with the PFM. Contrastingly, distributed aerial recharge in the tributary alluvium is recognized by Cl/Br ratios, groundwater ages, and noble gases nearly in equilibrium with the atmosphere. Distributed aerial recharge also fits the PFM because the sampled wells are close to the recharge zone.
Between the areas of active recharge, there is a zone of groundwater mixing and flow path convergence. This central zone of the aquifer is recognized based on agreement with the DM, older groundwater ages, and noble gases that are out of equilibrium with the atmosphere. Mixing is also supported by PCA, showing intermediate compositions of central groundwater and Cl/Br ratios. Longer residence times in this zone of the aquifer promote enhanced water-rock reaction and biogeochemical reactions as evidenced by REE relations and stable isotopes of nitrogen.
Dissolved solutes display spatial patterns that are controlled by the presence of anthropogenic point sources, as well as groundwater-residence times. Concentrations of a subset of major ions and trace elements (Cl, Cr, and Pb), a subset of REEs (Ce and Gd), and wastewater-indicator compounds are related primarily to proximity to surface-water bodies that lose flow to groundwater. Other constituents are derived from water-rock reactions and weathering within the aquifer matrix, such as HREE, U, Se, and 4He. This latter group of constituents tends to be enriched in the central portion of the aquifer where groundwater-residence times are longest.
There are several implications of these relations for future studies focused on hydrologic characterization in complex settings. First, a comprehensive dataset may allow for the full range of hydrologic conditions and geochemical processes to be recognized (e.g., see Reference [14]), but investigators may wish to tailor analytes if specific processes are key to the investigation. For example, without the use of REEs in this study, it may have been possible to recognize recharge from streamflow loss (by pharmaceutical and wastewater-compound indicators), but recognition of geochemical evolution in the central part of the aquifer would have been less evident. Moreover, the use of REEs is much more straightforward for recognizing the influence of WWTP effluent throughout the hydrologic system as the Gd anomaly alone is an indicator of effluent migration whereas use of pharmaceuticals and wastewater-indicator compounds is costlier analytically and may be complicated by complex quality-assurance aspects [16]. A second implication is that tracers of groundwater-residence time (be they quantitative or qualitative) are highly useful. Quantitative groundwater-residence time estimates allow for direct calculation of recharge rate [111] and are indicators of transience in the hydrologic system [110,112]. More qualitative estimates of residence times can be used to discriminate different flow-path configurations (i.e., piston flow, dispersion, and convergence) and may indicate the presence of old groundwater via noble-gas isotopes. Noble gases are also unique in their applicability for investigating recharge conditions such as groundwater recharge temperature [63], which was useful in this study for indicating a possible groundwater recharge source to the west.
Ideally, this conceptual model would be used to evaluate the distribution of specific compounds to attribute point sources of these compounds. Potential compounds of interest in this aquifer include pharmaceutical and wastewater-indicator compounds and PFAS [5]. This study includes data on pharmaceutical and wastewater-indicator compounds but does not directly include data for PFAS. Instead potential sources of PFAS are evaluated based on groundwater ages and records on industrial activity.
The most easily discriminated potential solute source for pharmaceutical and wastewater-indicator compounds is the WWTPs throughout the area, given that these compounds occur ubiquitously in surface waters and are found in groundwater of varying ages (i.e., Table 2 and Supplementary Materials Figure S6). Sources of PFAS are more difficult to attribute based on groundwater-age data alone because of the range of calculated groundwater-residence times. All groundwater observed in this study has at least a proportion of modern recharge, but in order to attribute certain compounds to a single source area the active period of solute loading from the given source area must be known [113]. PFAS compounds in the area have at least three possible sources: WWTP effluent [22], former metal-plating operations [114], and AFFF use [5]. These solute sources generally overlap in time (from the 1970s to 1990s for metal-plating operations [114], from the 1970s through the present for AFFF use [27], and for WWTP discharge [28]). Because most groundwater is less than 10 years old, it is plausible that PFAS compounds could be derived from either WWTP or AFFF, though metal-plating operations are less likely. To more fully investigate PFAS sourcing, it would be necessary to evaluate PFAS distributions and compositions, using statistical techniques (e.g., see Reference [10]).

4. Conclusions

An extensive geochemical dataset was collected, to investigate the hydrology and geochemistry of the Fountain Creek alluvial aquifer. This dataset included major and trace elements, stable isotopes, REEs, pharmaceutical and wastewater-indicator compounds, noble gases, and groundwater-age tracers. Each of these constituents has shown use in tracing geochemical reactions or resolving interactions between groundwater and surface water [7,11,14,63].
This study supports previous studies in indicating that groundwater recharge in the Fountain Creek alluvial aquifer is influenced by losses from surface water [18]. Groundwater can be separated into three distinct geochemical zones, which are spatially associated with preponderance of aerial recharge, proximity to surface waters, or groundwater mixing. The tributary alluvium, located in upgradient recharge areas of the aquifer, has a juvenile geochemical composition reflective of recent equilibrium with the atmosphere, lack of interaction with the aquifer materials, and generally young groundwater. The central groundwater zone displays geochemical evolution from the recharge zone and contains groundwater with the longest residence time. Groundwater bodies near Fountain Creek, the primary surface-water feature in the area, were also recently recharged, but they have the geochemical character of creek water that is influenced by WWTP effluent, and thus are partially derived from streamflow losses. The use of REEs is essential to understanding groundwater–surface-water interactions affected by WWTP effluent.
Most groundwater in the hydrological system appears to have been recharged within the past 20 years. However, noble-gas isotopes, lumped-parameter modeling, and the fraction of water recharged since 1953 (Fpost-1953) are consistent with a proportion of older groundwater in the system. The overall groundwater ages are likely biased young by the presence of recent recharge [65,67], and the presence of older groundwater may have gone undetected without the application of noble gases and their isotopes. More detailed evaluation of mixing of old and young waters would require an intermediate-age tracer, such as carbon-14, with the ability to define an old end member based on solute concentrations and detailed modeling (e.g., see Reference [11]).

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/13/6/871/s1, additional details on geochemical analyses within the Fountain Creek alluvial aquifer, Colorado. Specifically, details are included pertaining to major-ion and trace-element geochemistry, principal component analysis (PCA), tritium in precipitation, 3H/3He apparent ages, groundwater recharge temperature, and groundwater-flow direction.

Author Contributions

Conceptualization, S.S.P.; methodology, S.S.P., C.P.N., and G.K.; software, C.P.N.; validation, G.K. and C.P.N.; formal analysis, C.P.N.; investigation, S.S.P. and C.P.N.; resources, S.S.P. and G.K.; data curation, G.K. and C.P.N.; writing—original draft preparation, C.P.N.; writing—review and editing, C.P.N., G.K., and S.S.P.; visualization, C.P.N.; supervision, S.S.P.; project administration, S.S.P. and C.P.N.; funding acquisition, S.S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the US Air Force Civil Engineering Center.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data pertaining to this article are available from the USGS National Water Information Service (NWIS) database, at https://waterdata.usgs.gov/nwis, accessed on 1 January 2021.

Acknowledgments

The authors appreciate the work of Zach Kisfalusi, Mike Holmberg, Evan Gohring, and Nancy Bauch in collecting samples. This manuscript was substantially improved by comments provided by Kim Beisner and Zeno Levy. Comments from three anonymous reviewers were also helpful. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the US Government.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Böhlke, J.K.; Verstraeten, I.M.; Kraemer, T.F. Effects of surface-water irrigation on sources, fluxes, and residence times of water, nitrate, and uranium in an alluvial aquifer. Appl. Geochem. 2007, 22, 152–174. [Google Scholar] [CrossRef]
  2. Bai, X.; Lutz, A.; Carroll, R.; Keteles, K.; Dahlin, K.; Murphy, M.; Nguyen, D. Occurrence, distribution, and seasonality of emerging contaminants in urban watersheds. Chemosphere 2018, 200, 133–142. [Google Scholar] [CrossRef]
  3. Focazio, M.J.; Kolpin, D.W.; Barnes, K.K.; Furlong, E.T.; Meyer, M.T.; Zaugg, S.D.; Barber, L.B.; Thurman, M.E. A national reconnaissance for pharmaceuticals and other organic wastewater contaminants in the United States—II) Untreated drinking water sources. Sci. Total. Environ. 2008, 402, 201–216. [Google Scholar] [CrossRef]
  4. Elliott, S.M.; Erickson, M.L.; Krall, A.L.; Adams, B.A. Concentrations of pharmaceuticals and other micropollutants in groundwater downgradient from large on-site wastewater discharges. PLoS ONE 2018, 13, e0206004. [Google Scholar] [CrossRef]
  5. Hu, X.C.; Andrews, D.Q.; Lindstrom, A.B.; Bruton, T.A.; Schaider, L.A.; Grandjean, P.; Lohmann, R.; Carignan, C.C.; Blum, A.; Balan, S.A.; et al. Detection of Poly- and Perfluoroalkyl Substances (PFASs) in U.S. Drinking Water Linked to Industrial Sites, Military Fire Training Areas, and Wastewater Treatment Plants. Environ. Sci. Technol. Lett. 2016, 3, 344–350. [Google Scholar] [CrossRef]
  6. Weber, A.K.; Barber, L.B.; Leblanc, D.R.; Sunderland, E.M.; Vecitis, C.D. Geochemical and Hydrologic Factors Controlling Subsurface Transport of Poly- and Perfluoroalkyl Substances, Cape Cod, Massachusetts. Environ. Sci. Technol. 2017, 51, 4269–4279. [Google Scholar] [CrossRef] [PubMed]
  7. Keefe, S.; Barber, L.; Hubbard, L.; Bradley, P.; Roth, D.; Kolpin, D. Behavior of major and trace elements in a transient surface water/groundwater system following removal of a long-term wastewater treatment facility source. Sci. Total. Environ. 2019, 668, 867–880. [Google Scholar] [CrossRef] [PubMed]
  8. Böhlke, J.K.; Denver, J.M. Combined Use of Groundwater Dating, Chemical, and Isotopic Analyses to Resolve the History and Fate of Nitrate Contamination in Two Agricultural Watersheds, Atlantic Coastal Plain, Maryland. Water Resour. Res. 1995, 31, 2319–2339. [Google Scholar] [CrossRef]
  9. Bern, C.R.; Walton-Day, K.; Naftz, D.L. Improved enrichment factor calculations through principal component analysis: Examples from soils near breccia pipe uranium mines, Arizona, USA. Environ. Pollut. 2019, 248, 90–100. [Google Scholar] [CrossRef]
  10. Zhang, X.; Lohmann, R.; Dassuncao, C.; Hu, X.C.; Weber, A.K.; Vecitis, C.D.; Sunderland, E.M. Source Attribution of Poly- and Perfluoroalkyl Substances (PFASs) in Surface Waters from Rhode Island and the New York Metropolitan Area. Environ. Sci. Technol. Lett. 2016, 3, 316–321. [Google Scholar] [CrossRef] [Green Version]
  11. Jurgens, B.C.; Bexfield, L.M.; Eberts, S.M. A Ternary Age-Mixing Model to Explain Contaminant Occurrence in a Deep Supply Well. Ground Water 2014, 52, 25–39. [Google Scholar] [CrossRef]
  12. Kulongoski, J.T.; Hilton, D.R. Applications of groundwater helium. In Handbook of Environmental Isotope Geochemistry; Baskaran, M., Ed.; Springer: Berlin, Germany, 2012; pp. 285–304. [Google Scholar] [CrossRef]
  13. McMahon, P.B.; Lindsey, B.D.; Conlon, M.D.; Hunt, A.G.; Belitz, K.; Jurgens, B.C.; Varela, B.A. Hydrocarbons in Upland Groundwater, Marcellus Shale Region, Northeastern Pennsylvania and Southern New York, U.S.A. Environ. Sci. Technol. 2019, 53, 8027–8035. [Google Scholar] [CrossRef]
  14. Moeck, C.; Radny, D.; Popp, A.; Brennwald, M.; Stoll, S.; Auckenthaler, A.; Berg, M.; Schirmer, M. Characterization of a managed aquifer recharge system using multiple tracers. Sci. Total. Environ. 2017, 609, 701–714. [Google Scholar] [CrossRef]
  15. Levy, Z.F.; Fram, M.S.; Faulkner, K.E.; Alpers, C.N.; Soltero, E.M.; Taylor, K.A. Effects of montane watershed development on vulnerability of domestic groundwater supply during drought. J. Hydrol. 2020, 583, 124567. [Google Scholar] [CrossRef]
  16. Bexfield, L.M.; Toccalino, P.L.; Belitz, K.; Foreman, W.T.; Furlong, E.T. Hormones and Pharmaceuticals in Groundwater Used As a Source of Drinking Water Across the United States. Environ. Sci. Technol. 2019, 53, 2950–2960. [Google Scholar] [CrossRef] [Green Version]
  17. Arnold, L.R.; Ortiz, R.F.; Brown, C.R.; Watts, K.R. Groundwater and Surface-Water Interaction, Water Quality, and Processes Affecting Loads of Dissolved Solids, Selenium, and Uranium in Fountain Creek, Pueblo County, Colorado, 2012–2014. In U.S. Geological Survey Scientific Investigation Report 2016–5134; U.S. Geological Survey: Reston, VA, USA, 2016; 78p. [Google Scholar] [CrossRef] [Green Version]
  18. Cain, D.; Edelmann, P. A reconnaissance water-quality appraisal of the Fountain Creek alluvial aquifer between Colorado Springs and Pueblo, Colorado, including trace elements and organic constituents. In U.S. Geological Survey Water-Resources Investigations Report 86-4085; U.S. Geological Survey: Reston, VA, USA, 1986; 45p. [Google Scholar]
  19. Chafin, D.T. Effects of land use on water quality of the Fountain Creek alluvial aquifer, east-central Colorado. In U.S. Geological Survey Water-Supply Paper 2381-D; U.S. Geological Survey: Reston, VA, USA, 1996; 99p. [Google Scholar]
  20. Edelmann, P.; Cain, D. Sources of water and nitrogen to the Widefield aquifer, southwestern El Paso County, Colorado. In U.S. Geological Survey Water-Resources Investigations Report 85-4162; U.S. Geological Survey: Reston, VA, USA, 1985; 81p. [Google Scholar]
  21. Lewis, M.E. Quality of water in the alluvial aquifer and tributary alluvium of the Fountain Creek valley, southwestern El Paso County, Colorado, 1991–1992. In U.S. Geological Survey Water-Resources Investigations Report 94-4118; U.S. Geological Survey: Reston, VA, USA, 1995; 39p. [Google Scholar]
  22. Mau, D.P.; Stogner, R.W.; Edelmann, P. Characterization of stormflows and wastewater treatment-plant effluent discharges on water quality, suspended sediment, and stream morphology for Fountain and Monument Creek watersheds, Colorado, 1981–2006. In U.S. Geological Survey Scientific Investigations Report 2007–5104; U.S. Geological Survey: Reston, VA, USA, 2007; 76p. [Google Scholar]
  23. U.S. Environmental Protection Agency. Lifetime Health Advisories and Health Effects Support Documents for Per-fluorooctanoic Acid and Perfluorooctane Sulfonate. Fed. Regist. 2016, 81, 33250–33251. [Google Scholar]
  24. U.S. Geological Survey. USGS water data for the Nation. In U.S. Geological Survey National Water Information System Database; U.S. Geological Survey: Reston, VA, USA, 2020. [Google Scholar] [CrossRef]
  25. AECOM. Peterson Air Force Base Conceptual Site Model: Report Prepared for US; Air Force Civil Engineering Center by AECOM: Denver, CO, USA, 2017; 24p. [Google Scholar]
  26. Topper, R.; Spray, K.L.; Bellis, W.H.; Hamilton, J.L.; Barkmann, P.E. Ground Water Atlas of Colorado: Colorado Geological Survey Special Publications-53. 2003. Available online: https://coloradogeologicalsurvey.org/water/colorado-groundwater-atlas/ (accessed on 1 January 2021).
  27. Aerostar. Final Preliminary Assessment Report for Perfluorinated Compounds at Peterson Air Force Base El Paso County, Colorado; Air Force Civil Engineering Center, Aerostar SES LLC.: Jacksonville, FL, USA, 2016; 145p. [Google Scholar]
  28. U.S. Environmental Protection Agency. Facility Registry Service—Geospatial Downloads. 2019. Available online: https://www.epa.gov/frs/geospatial-data-download-service (accessed on 1 April 2019).
  29. Kuhn, G.; Arnold, L.R. Application of a stream-aquifer model to Monument Creek for development of a method to estimate transit losses for reusable water, El Paso County, Colorado. In U.S. Geological Survey Scientific Investigations Report 2006–5184; U.S. Geological Survey: Reston, VA, USA, 2006; 111p. [Google Scholar]
  30. U.S. Geological Survey. Transit-Loss Accounting Program for Monument and Fountain Creeks, El Paso and Pueblo Counties, Colorado. 2019. Available online: https://www.usgs.gov/centers/co-water/science/transit-loss-accounting-program-monument-and-fountain-creeks-el-paso-and?qt-science_center_objects=0#qt-science_center_objects (accessed on 1 October 2019).
  31. Huntington, J.L.; Hegewisch, K.C.; Daudert, B.; Morton, C.G.; Abatzoglou, J.T.; McEvoy, D.J.; Erickson, T. Climate Engine: Cloud Computing and Visualization of Climate and Remote Sensing Data for Advanced Natural Resource Monitoring and Process Understanding. Bull. Am. Meteorol. Soc. 2017, 98, 2397–2410. [Google Scholar] [CrossRef]
  32. Jenkins, E.D. Ground water in Fountain and Jimmy Camp Valleys, El Paso County, Colorado. In U.S. Geological Survey Water-Supply Paper 1583; U.S. Geological Survey: Reston, VA, USA, 1964; 66p. [Google Scholar]
  33. U.S. Geological Survey; Variously Dated; National Field Manual for the Collection of Water-Quality Data. U.S. Geological Survey Techniques of Water-Resources Investigations, Book 9, Chaps A1–A10. Available online: http://pubs.water.usgs.gov/twri9A (accessed on 10 November 2019).
  34. Rounds, S.A. Alkalinity and acid neutralizing capacity (ver. 4.0, September 2012). In U.S. Geological Survey Techniques of Water-Resources Investigations; U.S. Geological Survey: Reston, VA, USA, 2012; Book 9, Chapter 6.6. [Google Scholar] [CrossRef]
  35. Aeschbach-Hertig, W.; Solomon, D.K. Noble gas thermometry in groundwater hydrology. In The Noble Gases as Geochemical Tracers; Burnard, P., Ed.; Springer: Berlin, Germany, 2013; pp. 81–122. [Google Scholar]
  36. Fishman, M.J.; Friedman, L.C. (Eds.) Methods for determination of inorganic substances in water and fluvial sediments. In U.S. Geological Survey Techniques of Water-Resources Investigations; U.S. Geological Survey: Reston, VA, USA, 1989; Book 5, Chapter A1; 545p. Available online: http://pubs.usgs.gov/twri/twri5-a1/html/pdf.html (accessed on 27 February 2014).
  37. Fishman, M.J. (Ed.) Methods of analysis by the U.S. Geological Survey National Water Quality Laborato-ry—Determination of inorganic and organic constituents in water and fluvial sediments. In U.S. Geological Survey OpenFile Report 93–125; U.S. Geological Survey: Reston, VA, USA, 1993; 217p. Available online: http://pubs.usgs.gov/of/1993/0125/report.pdf (accessed on 27 February 2014).
  38. Garbarino, J.R.; Struzeski, T.M. Methods of analysis by the U.S. Geological Survey National Water Quality Labor-atory—Determination of elements in whole-water digests using inductively coupled plasma-optical emission spectrometry and inductively coupled plasma-mass spectrometry. In U.S. Geological Survey Open-File Report 98–165; U.S. Geological Survey: Reston, VA, USA, 1998; 101p. Available online: http://nwql.usgs.gov/OFR-98-165.shtml (accessed on 27 February 2014).
  39. Furlong, E.T.; Noriega, M.C.; Kanagy, C.J.; Kanagy, L.K.; Coffey, L.J.; Burkhardt, M.R. Determination of human-use pharmaceuticals in filtered water by direct aqueous injection–high-performance liquid chromatography/tandem mass spec-trometry. In U.S. Geological Survey Techniques and Methods; U.S. Geological Survey: Reston, VA, USA, 2014; Book 5, Chapter B10; 49p. [Google Scholar]
  40. Coplen, T.B.; Qi, H.; Révész, K.; Casciotti, K.; Hannon, J.E. Determination of the δ15N and δ18O of Nitrate in Water; RSIL lab code 2900, chap. 17 of Stable Isotope-Ratio Methods, Methods of the Reston Stable Isotope Laboratory (Slightly revised from Version 1.0 Released in 2007). In U.S. Geological Survey Techniques and Methods; Révész, K., Coplen, T.B., Eds.; U.S. Geological Survey: Reston, VA, USA, 2012; 35p. Available online: https://pubs.usgs.gov/tm/2006/tm10c17/ (accessed on 1 January 2021).
  41. Kendall, C.; Doctor, D.H.; Young, M.B. Environmental isotope applications in hydrologic studies. In Treatise on Geochemistry, 2nd ed.; Holland, H.D., Turekian, K.K., Eds.; Elsevier Science: Amsterdam, The Netherlands, 2015; pp. 273–327. [Google Scholar]
  42. Révész, K.; Coplen, T.B. Determination of the δ(2H/1H) of water: RSIL lab code 1574. Chapter C1, Methods of the Reston Stable Isotope Laboratory. In U.S. Geological Survey Techniques and Methods 10–C1; Révész, K., Coplen, T.B., Eds.; U.S. Geological Survey: Reston, VA, USA, 2008; 27p. [Google Scholar]
  43. Révész, K.; Coplen, T.B. Determination of the δ(18O/16O) of water: RSIL lab code 489. Chapter C2, Methods of the Reston Stable Isotope Laboratory. In U.S. Geological Survey Techniques and Methods, 10–C2; Révész, K., Coplen, T.B., Eds.; U.S. Geological Survey: Reston, VA, USA, 2008; 28p. [Google Scholar]
  44. Hunt, A.G. Noble Gas Laboratory’s standard operating procedures for the measurement of dissolved gas in water samples. In U.S. Geological Survey Techniques and Methods; U.S. Geological Survey: Reston, VA, USA, 2015; Book 5, Chapter A11; 22p. [Google Scholar] [CrossRef] [Green Version]
  45. Busenberg, E.; Plummer, L.N. Dating young groundwater with sulfur hexafluoride: Natural and anthropogenic sources of sulfur hexafluoride. Water Resour. Res. 2000, 36, 3011–3030. [Google Scholar] [CrossRef]
  46. Mueller, D.K.; Schertz, T.L.; Martin, J.D.; Sandstrom, M.W. Design, analysis, and interpretation of field quali-ty-control data for water-sampling projects. In U.S. Geological Survey Techniques and Methods; U.S. Geological Survey: Reston, VA, USA, 2015; Book 4, Chapter C4; 54p. [Google Scholar] [CrossRef]
  47. Appelo, C.A.J.; Postma, D. Geochemistry, Groundwater, and Pollution, 2nd ed.; CRC Press: Amsterdam, The Netherlands, 2005; 649p. [Google Scholar]
  48. Lee, L.; Helsel, D. Baseline models of trace elements in major aquifers of the United States. Appl. Geochem. 2005, 20, 1560–1570. [Google Scholar] [CrossRef]
  49. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: http://www.R-project.org/ (accessed on 5 June 2019).
  50. Reimann, R.C.; Filzmoser, P.; Garrett, R.G.; Dutter, R. Statistical Data Analysis Explained: Applied Environmental Statistics; John Wiley & Sons, Ltd.: Chichester, UK, 2008. [Google Scholar]
  51. Blondes, M.S.; Engle, M.A.; Geboy, N.J. A practical guide to the use of major elements, trace elements, and isotopes in compositional data analysis: Applications for deep formation brine geochemistry. In Compositional Data Analysis; Martín-Fernández, J.A., Thió-Henestrosa, S., Eds.; Springer: Cham, Switzerland, 2016; pp. 13–29. [Google Scholar]
  52. Filzmoser, P.; Hron, K.; Reimann, C. Principal component analysis for compositional data with outliers. Environmetrics 2009, 20, 621–632. [Google Scholar] [CrossRef]
  53. Jurgens, B.C.; Böhlke, J.K.; Eberts, S.M. TracerLPM (Version 1): An Excel® workbook for interpreting groundwater age distributions from environmental tracer data. In U.S. Geological Survey Techniques and Methods Report 4-F3; U.S. Geological Survey: Reston, VA, USA, 2012; 60p. [Google Scholar]
  54. Solomon, D.K.; Cook, P.G. 3H and 3He. In Environmental Tracers in Subsurface Hydrology; Cook, P.G., Herczeg, A.L., Eds.; Springer Science: New York, NY, USA, 2000; pp. 397–424. [Google Scholar]
  55. Chambers, L.; Gooddy, D.; Binley, A. Use and application of CFC-11, CFC-12, CFC-113 and SF6 as environmental tracers of groundwater residence time: A review. Geosci. Front. 2019, 10, 1643–1652. [Google Scholar] [CrossRef]
  56. Cartwright, I.; Cendón, D.; Currell, M.; Meredith, K. A review of radioactive isotopes and other residence time tracers in understanding groundwater recharge: Possibilities, challenges, and limitations. J. Hydrol. 2017, 555, 797–811. [Google Scholar] [CrossRef]
  57. Weiss, R.F. The solubility of nitrogen, oxygen and argon in water and seawater. Deep Sea Res. 1970, 17, 721–735. [Google Scholar] [CrossRef]
  58. Weiss, R. Carbon dioxide in water and seawater: The solubility of a non-ideal gas. Mar. Chem. 1974, 2, 203–215. [Google Scholar] [CrossRef]
  59. Benson, B.B.; Krause, D. The concentration and isotopic fractionation of gases dissolved in freshwater in equilibrium with the atmosphere. 1. Oxygen. Limnol. Oceanogr. 1980, 25, 662–671. [Google Scholar] [CrossRef] [Green Version]
  60. Benson, B.B.; Krause, D. The concentration and isotopic fractionation of oxygen dissolved in freshwater and seawater in equilibrium with the atmosphere1. Limnol. Oceanogr. 1984, 29, 620–632. [Google Scholar] [CrossRef]
  61. Peeters, F.; Beyerle, U.; Aeschbach-Hertig, W.; Holocher, J.; Brennwald, M.S.; Kipfer, R. Improving noble gas based paleoclimate reconstruction and groundwater dating using 20Ne/22Ne ratios. Geochim. Cosmochim. Acta 2003, 67, 587–600. [Google Scholar] [CrossRef] [Green Version]
  62. Manning, A.H.; Solomon, D. Using noble gases to investigate mountain-front recharge. J. Hydrol. 2003, 275, 194–207. [Google Scholar] [CrossRef]
  63. Manning, A.H. Mountain-block recharge, present and past, in the eastern Española Basin, New Mexico, USA. Hydrogeol. J. 2011, 19, 379–397. [Google Scholar] [CrossRef]
  64. Lindsey, B.D.; Jurgens, B.C.; Belitz, K. Tritium as an indicator of modern, mixed, and premodern groundwater age. Sci. Investig. Rep. 2019, 18. [Google Scholar] [CrossRef]
  65. McCallum, J.L.; Cook, P.G.; Simmons, C.T. Limitations of the Use of Environmental Tracers to Infer Groundwater Age. Ground Water 2015, 53, 56–70. [Google Scholar] [CrossRef]
  66. Cook, P.G.; Böhlke, J.K. Determining timescales for groundwater flow and transport. In Environmental Tracers in Subsurface Hydrology; Cook, P.G., Herczeg, A.L., Eds.; Springer Science: New York, NY, USA, 2000; pp. 1–30. [Google Scholar]
  67. Suckow, A. The age of groundwater—Definitions, models and why we do not need this term. Appl. Geochem. 2014, 50, 222–230. [Google Scholar] [CrossRef]
  68. Leray, S.; Engdahl, N.B.; Massoudieh, A.; Bresciani, E.; McCallum, J. Residence time distributions for hydrologic systems: Mechanistic foundations and steady-state analytical solutions. J. Hydrol. 2016, 543, 67–87. [Google Scholar] [CrossRef] [Green Version]
  69. Jasechko, S. Partitioning young and old groundwater with geochemical tracers. Chem. Geol. 2016, 427, 35–42. [Google Scholar] [CrossRef]
  70. Jurgens, B.C. Data for Tritium deposition in precipitation in the United States, 1953–2012. In U.S. Geological Survey Data Release; U.S. Geological Survey: Reston, VA, USA, 2018. [Google Scholar] [CrossRef]
  71. Jasechko, S.; Taylor, R.G. Intensive rainfall recharges tropical groundwaters. Environ. Res. Lett. 2015, 10, 124015. [Google Scholar] [CrossRef]
  72. Michel, R.L.; Jurgens, B.C.; Young, M.B. Tritium deposition in precipitation in the United States, 1953–2012. Sci. Investig. Rep. 2018, 11. [Google Scholar] [CrossRef] [Green Version]
  73. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  74. Alcalá, F.J.; Custodio, E. Using the Cl/Br ratio as a tracer to identify the origin of salinity in aquifers in Spain and Portugal. J. Hydrol. 2008, 359, 189–207. [Google Scholar] [CrossRef]
  75. Davis, S.N.; Whittemore, D.O.; Fabryka-Martin, J. Uses of Chloride/Bromide Ratios in Studies of Potable Water. Ground Water 1998, 36, 338–350. [Google Scholar] [CrossRef]
  76. McArthur, J.; Sikdar, P.; Hoque, M.; Ghosal, U. Waste-water impacts on groundwater: Cl/Br ratios and implications for arsenic pollution of groundwater in the Bengal Basin and Red River Basin, Vietnam. Sci. Total. Environ. 2012, 437, 390–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Jurgens, B.C.; McMahon, P.B.; Chapelle, F.H.; Eberts, S.M. An Excel® workbook for identifying redox processes in ground water. In U.S. Geological Survey Open-File Report 2009–1004; U.S. Geological Survey: Reston, VA, USA, 2009; 8p. [Google Scholar]
  78. Anderson, L.; Berkelhammer, M.; Mast, M.A. Isotopes in North American Rocky Mountain Snowpack 1993–2014. Quat. Sci. Rev. 2016, 131, 262–273. [Google Scholar] [CrossRef]
  79. Bowen, G.J. The Online Isotopes in Precipitation Calculator, Version 3.1. 2020. Available online: http://www.waterisotopes.org (accessed on 2 July 2020).
  80. Solder, J.E.; Beisner, K.R. Critical evaluation of stable isotope mixing end-members for estimating groundwater recharge sources: Case study from the South Rim of the Grand Canyon, Arizona, USA. Hydrogeol. J. 2020, 28, 1–17. [Google Scholar] [CrossRef]
  81. Gammons, C.H.; Poulson, S.R.; Pellicori, D.A.; Reed, P.J.; Roesler, A.J.; Petrescu, E.M. The hydrogen and oxygen isotopic composition of precipitation, evaporated mine water, and river water in Montana, USA. J. Hydrol. 2006, 328, 319–330. [Google Scholar] [CrossRef]
  82. McMahon, P.B.; Böhlke, J.K. Regional Patterns in the Isotopic Composition of Natural and Anthropogenic Nitrate in Groundwater, High Plains, U.S.A. Environ. Sci. Technol. 2006, 40, 2965–2970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Kendall, C.; Aravena. Nitrate isotopes in groundwater systems. In Environmental Tracers in Subsurface Hydrology; Cook, P.G., Herczeg, A.L., Eds.; Springer Science: New York, NY, USA, 2000; pp. 261–297. [Google Scholar]
  84. Leybourne, M.I.; Johannesson, K.H. Rare earth elements (REE) and yttrium in stream waters, stream sediments, and Fe–Mn oxyhydroxides: Fractionation, speciation, and controls over REE+Y patterns in the surface environment. Geochim. Cosmochim. Acta 2008, 72, 5962–5983. [Google Scholar] [CrossRef]
  85. Duvert, C.; Cendón, D.I.; Raiber, M.; Seidel, J.-L.; Cox, M.E. Seasonal and spatial variations in rare earth elements to identify inter-aquifer linkages and recharge processes in an Australian catchment. Chem. Geol. 2015, 396, 83–97. [Google Scholar] [CrossRef]
  86. Göb, S.; Loges, A.; Nolde, N.; Bau, M.; Jacob, D.E.; Markl, G. Major and trace element compositions (including REE) of mineral, thermal, mine and surface waters in SW Germany and implications for water–rock interaction. Appl. Geochem. 2013, 33, 127–152. [Google Scholar] [CrossRef]
  87. Rabiet, M.; Togola, A.; Brissaud, F.; Seidel, J.-L.; Budzinski, H.; Elbaz-Poulichet, F. Consequences of Treated Water Recycling as Regards Pharmaceuticals and Drugs in Surface and Ground Waters of a Medium-sized Mediterranean Catchment. Environ. Sci. Technol. 2006, 40, 5282–5288. [Google Scholar] [CrossRef]
  88. Barber, L.B.; Paschke, S.S.; Battaglin, W.A.; Douville, C.; Fitzgerald, K.C.; Keefe, S.H.; Roth, D.A.; Vajda, A.M. Effects of an Extreme Flood on Trace Elements in River Water—From Urban Stream to Major River Basin. Environ. Sci. Technol. 2017, 51, 10344–10356. [Google Scholar] [CrossRef]
  89. Verplanck, P.L.; Furlong, E.T.; Gray, J.L.; Phillips, P.J.; Wolf, R.E.; Esposito, K. Evaluating the Behavior of Gadolinium and Other Rare Earth Elements through Large Metropolitan Sewage Treatment Plants. Environ. Sci. Technol. 2010, 44, 3876–3882. [Google Scholar] [CrossRef]
  90. Noack, C.W.; Dzombak, D.A.; Karamalidis, A.K. Rare Earth Element Distributions and Trends in Natural Waters with a Focus on Groundwater. Environ. Sci. Technol. 2014, 48, 4317–4326. [Google Scholar] [CrossRef]
  91. Gromet, L.; Haskin, L.A.; Korotev, R.L.; Dymek, R.F. The “North American shale composite”: Its compilation, major and trace element characteristics. Geochim. Cosmochim. Acta 1984, 48, 2469–2482. [Google Scholar] [CrossRef]
  92. Stolpe, B.; Guo, L.; Shiller, A.M. Binding and transport of rare earth elements by organic and iron-rich nanocolloids in Alaskan rivers, as revealed by field-flow fractionation and ICP-MS. Geochim. Cosmochim. Acta 2013, 106, 446–462. [Google Scholar] [CrossRef]
  93. Barnes, K.K.; Kolpin, D.W.; Furlong, E.T.; Zaugg, S.D.; Meyer, M.T.; Barber, L.B. A national reconnaissance of pharmaceuticals and other organic wastewater contaminants in the United States—I) Groundwater. Sci. Total. Environ. 2008, 402, 192–200. [Google Scholar] [CrossRef]
  94. Schaider, L.A.; Rodgers, K.M.; Rudel, R.A. Review of Organic Wastewater Compound Concentrations and Removal in Onsite Wastewater Treatment Systems. Environ. Sci. Technol. 2017, 51, 7304–7317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. U.S. Environmental Protection Agency. Reassessments of Two Exemptions from the Requirement of a Tolerance for Hexamethylenetetramine (HMTA; CAS Reg. No. 100-97-0). 2006. Available online: https://www.epa.gov/sites/production/files/2015-04/documents/hexamethyl.pdf (accessed on 19 March 2020).
  96. National Center for Biotechnology Information. PubChem Database, Methenamine, CID=4101. 2020. Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Methenamine (accessed on 2 July 2020).
  97. Gardner, P.M.; Heilweil, V.M. A multiple-tracer approach to understanding regional groundwater flow in the Snake Valley area of the eastern Great Basin, USA. Appl. Geochem. 2014, 45, 33–49. [Google Scholar] [CrossRef]
  98. Kipfer, R.; Aeschbach-Hertig, W.; Peeters, F.; Stute, M. Noble Gases in Lakes and Ground Waters. Rev. Miner. Geochem. 2002, 47, 615–700. [Google Scholar] [CrossRef]
  99. Kulongoski, J.T.; Hilton, D.R.; Cresswell, R.G.; Hostetler, S.; Jacobson, G. Helium-4 characteristics of groundwaters from Central Australia: Comparative chronology with chlorine-36 and carbon-14 dating techniques. J. Hydrol. 2008, 348, 176–194. [Google Scholar] [CrossRef]
  100. Zinn, B.A.; Konikow, L.F. Potential effects of regional pumpage on groundwater age distribution. Water Resour. Res. 2007, 43, 06418. [Google Scholar] [CrossRef]
  101. Porcelli, D.; Ballentine, C.J.; Wieler, R. An Overview of Noble Gas Geochemistry and Cosmochemistry. Rev. Miner. Geochem. 2002, 47, 1–19. [Google Scholar] [CrossRef]
  102. Mahara, Y.; Ohta, T.; Morikawa, N.; Nakano, T.; Tokumasu, M.; Hukutani, S.; Tokunaga, T.; Igarashi, T. Effects of terrigenic He components on tritium–helium dating: A case study of shallow groundwater in the Saijo Basin. Appl. Geochem. 2014, 50, 142–149. [Google Scholar] [CrossRef] [Green Version]
  103. Visser, A.; Schaap, J.D.; Broers, H.P.; Bierkens, M.F. Degassing of 3H/3He, CFCs and SF6 by denitrification: Measurements and two-phase transport simulations. J. Contam. Hydrol. 2009, 103, 206–218. [Google Scholar] [CrossRef]
  104. Newman, C.P. Environmental-tracer modeling to support hydrogeochemical evaluation of the Fountain Creek alluvial aquifer, El Paso County, Colorado. In U.S. Geological Survey Data Release; U.S. Geological Survey: Reston, VA, USA, 2020. [Google Scholar] [CrossRef]
  105. Koh, D.-C.; Plummer, L.N.; Busenberg, E.; Kim, Y. Evidence for terrigenic SF6 in groundwater from basaltic aquifers, Jeju Island, Korea: Implications for groundwater dating. J. Hydrol. 2007, 339, 93–104. [Google Scholar] [CrossRef]
  106. Von Rohden, C.; Kreuzer, A.; Chen, Z.; Aeschbach-Hertig, W. Accumulation of natural SF6in the sedimentary aquifers of the North China Plain as a restriction on groundwater dating. Isot. Environ. Health Stud. 2010, 46, 279–290. [Google Scholar] [CrossRef]
  107. Santella, N.; Ho, D.T.; Schlosser, P.; Stute, M. Widespread elevated atmospheric SF6 mixing ratios in the Northeastern United States: Implications for groundwater dating. J. Hydrol. 2008, 349, 139–146. [Google Scholar] [CrossRef]
  108. Friedrich, R.; Vero, G.; von Rohden, C.; Lessmann, B.; Kipfer, R.; Aeschbach-Hertig, W. Factors controlling terrigenic SF6 in young groundwater of the Odenwald region (Germany). Appl. Geochem. 2013, 33, 318–329. [Google Scholar] [CrossRef]
  109. Cey, B.D.; Hudson, G.B.; Moran, J.E.; Scanlon, B.R. Evaluation of Noble Gas Recharge Temperatures in a Shallow Unconfined Aquifer. Ground Water 2009, 47, 646–659. [Google Scholar] [CrossRef] [PubMed]
  110. Manning, A.H.; Clark, J.F.; Diaz, S.H.; Rademacher, L.K.; Earman, S.; Plummer, L.N. Evolution of groundwater age in a mountain watershed over a period of thirteen years. J. Hydrol. 2012, 460-461, 13–28. [Google Scholar] [CrossRef]
  111. McMahon, P.B.; Plummer, L.N.; Böhlke, J.K.; Shapiro, S.D.; Hinkle, S.R. A comparison of recharge rates in aquifers of the United States based on groundwater-age data. Hydrogeol. J. 2011, 19, 779–800. [Google Scholar] [CrossRef]
  112. McCallum, J.L.; Dogramaci, S.; Bai, A.; Cook, P.G.; Engdahl, N.B.; Simmons, C.T.; Skrzypek, G.; Grierson, P.F. Assessing Temporal Changes in Groundwater Recharge Using Spatial Variations in Groundwater Ages. Water Resour. Res. 2020, 56. [Google Scholar] [CrossRef]
  113. Kaown, D.; Koh, D.-C.; Solomon, D.K.; Yoon, Y.-Y.; Yang, J.; Lee, K.-K. Delineation of recharge patterns and contaminant transport using 3H–3He in a shallow aquifer contaminated by chlorinated solvents in South Korea. Hydrogeol. J. 2014, 22, 1041–1054. [Google Scholar] [CrossRef]
  114. Colorado Department of Public Health and the Environment. Evaluation of Tetrachloroethylene Vapor Intrusion into Buildings Located above a Contaminated Aquifer, Schlage Lock Company, Security, El Paso County, Colo-Rado, Epa Facility ID: COD082657420. 2006; 33p. Available online: https://www.atsdr.cdc.gov/HAC/pha/SchlageLockCompany/SchlageLockCompanyHC113006.pdf (accessed on 25 March 2020).
Figure 1. Map of study site with sampling locations, referenced according to the National Water Information System (NWIS) database [24], surface-water features, and extent of the alluvial aquifer [25,26]. Locations of potential solute sources are from References [27,28] and include wastewater treatment plants (WWTPs), areas of aqueous-film-forming foam (AFFF) use, and former metal-plating facilities. SW and GW refer to surface water and groundwater sample locations, respectively.
Figure 1. Map of study site with sampling locations, referenced according to the National Water Information System (NWIS) database [24], surface-water features, and extent of the alluvial aquifer [25,26]. Locations of potential solute sources are from References [27,28] and include wastewater treatment plants (WWTPs), areas of aqueous-film-forming foam (AFFF) use, and former metal-plating facilities. SW and GW refer to surface water and groundwater sample locations, respectively.
Water 13 00871 g001
Figure 2. (a) Stable-isotopic composition and regression lines for groundwater and surface-water samples (snowpack-driven Rocky Mountain Meteoric Water Line (RMWL-Snow) from Reference [78] and isotopes in precipitation calculator (OIPC) meteoric water line (MWL) from Reference [79]), (b) δ18OH2O vs. δ18ONO3 with specific wells labeled with site alias and lines representing theoretical sources of oxygen during denitrification reactions [41,82], and (c) δ15N vs. δ18ONO3 with specific wells labeled with site alias and lines for ranges of fractionation trends applicable to nitrate reduction [82].
Figure 2. (a) Stable-isotopic composition and regression lines for groundwater and surface-water samples (snowpack-driven Rocky Mountain Meteoric Water Line (RMWL-Snow) from Reference [78] and isotopes in precipitation calculator (OIPC) meteoric water line (MWL) from Reference [79]), (b) δ18OH2O vs. δ18ONO3 with specific wells labeled with site alias and lines representing theoretical sources of oxygen during denitrification reactions [41,82], and (c) δ15N vs. δ18ONO3 with specific wells labeled with site alias and lines for ranges of fractionation trends applicable to nitrate reduction [82].
Water 13 00871 g002
Figure 3. (a) NO3 + NO2 concentrations versus δ15N and (b) NO3 + NO2 concentrations versus ammonia concentrations in groundwater and surface water with specific wells labeled with site alias. Isotopes of N and concentrations of aqueous species are consistent with denitrification reactions in water from wells. The circled area on panel (b) corresponds to generally elevated ammonia and NO3 + NO2 concentrations likely influenced by WWTPs. The two plots do not show all the same wells, because not all wells that were sampled for ammonia and NO3 + NO2 were also sampled for nitrogen isotopes.
Figure 3. (a) NO3 + NO2 concentrations versus δ15N and (b) NO3 + NO2 concentrations versus ammonia concentrations in groundwater and surface water with specific wells labeled with site alias. Isotopes of N and concentrations of aqueous species are consistent with denitrification reactions in water from wells. The circled area on panel (b) corresponds to generally elevated ammonia and NO3 + NO2 concentrations likely influenced by WWTPs. The two plots do not show all the same wells, because not all wells that were sampled for ammonia and NO3 + NO2 were also sampled for nitrogen isotopes.
Water 13 00871 g003
Figure 4. (a) Evaluation of rare earth element (REE) concentrations in groundwater (GW) and surface water (SW). Panels (ad) show REE concentrations normalized to North American Shale Composite (NASC). Prominent Gd anomalies are present in Fountain Creek and groundwater proximal to the creek. Panel (e) shows a REE normalization plot (e.g., see Reference [92]) and conceptual mixing relations. A synthetic mixing line between groundwater composition (T02-MW006) and WWTP-affected surface water (Canal 4 at Headgate) is labeled with the fraction of groundwater in the mixture (x). Arrows at top right of REE normalization diagram indicate the general patterns expected by mixing relations and water–rock reaction (rxn).
Figure 4. (a) Evaluation of rare earth element (REE) concentrations in groundwater (GW) and surface water (SW). Panels (ad) show REE concentrations normalized to North American Shale Composite (NASC). Prominent Gd anomalies are present in Fountain Creek and groundwater proximal to the creek. Panel (e) shows a REE normalization plot (e.g., see Reference [92]) and conceptual mixing relations. A synthetic mixing line between groundwater composition (T02-MW006) and WWTP-affected surface water (Canal 4 at Headgate) is labeled with the fraction of groundwater in the mixture (x). Arrows at top right of REE normalization diagram indicate the general patterns expected by mixing relations and water–rock reaction (rxn).
Water 13 00871 g004
Figure 5. Results of principal component analysis (PCA) for subset of geochemical parameters. Rays are plotted against rotations (top and right axes), and sample results are plotted against principal component loadings (bottom and left axes).
Figure 5. Results of principal component analysis (PCA) for subset of geochemical parameters. Rays are plotted against rotations (top and right axes), and sample results are plotted against principal component loadings (bottom and left axes).
Water 13 00871 g005
Figure 7. Selected results of lumped parameter modeling: (a) 3Hetrit versus 3H, (b) 3Hetrit versus SF6, and (c) SF6 versus 3H. Ages are marked by symbols along lines for piston-flow model (PFM), exponential-mixing model (EMM), and dispersion model (DM).
Figure 7. Selected results of lumped parameter modeling: (a) 3Hetrit versus 3H, (b) 3Hetrit versus SF6, and (c) SF6 versus 3H. Ages are marked by symbols along lines for piston-flow model (PFM), exponential-mixing model (EMM), and dispersion model (DM).
Water 13 00871 g007
Figure 8. Histogram of estimates of mean residence time from lumped-parameter modeling, 3H/3He apparent ages, and idealized SF6 piston-flow ages. Different numbers of analyses are available because multiple age estimates are produced from the program TracerLPM, using the lookup mean ages function (see Reference [53], for method details), and because some model calculations fail for idealized SF6 piston flow (see Table 2).
Figure 8. Histogram of estimates of mean residence time from lumped-parameter modeling, 3H/3He apparent ages, and idealized SF6 piston-flow ages. Different numbers of analyses are available because multiple age estimates are produced from the program TracerLPM, using the lookup mean ages function (see Reference [53], for method details), and because some model calculations fail for idealized SF6 piston flow (see Table 2).
Water 13 00871 g008
Figure 9. Conceptual hydrologic model highlighting processes occurring in the system (indicated by black type and solid black lines) and the means by which those processes were recognized (indicated by underlined type and dashed black lines).
Figure 9. Conceptual hydrologic model highlighting processes occurring in the system (indicated by black type and solid black lines) and the means by which those processes were recognized (indicated by underlined type and dashed black lines).
Water 13 00871 g009
Table 1. Summary of percent detections of contaminants of emerging concern in water-quality dataset from groundwater and surface water associated with the Fountain Creek alluvial aquifer and occurrence of compounds from similar studies within the United States.
Table 1. Summary of percent detections of contaminants of emerging concern in water-quality dataset from groundwater and surface water associated with the Fountain Creek alluvial aquifer and occurrence of compounds from similar studies within the United States.
Compound NamePercent Detections in Groundwater
Samples
Percent Detections in Surface-Water
Samples
Occurrence in Groundwater
Downgradient from Wastewater
Effluent in Minnesota, USA 1
Occurrence in Surface Water near Denver, Colorado, USA 2
Methylbenzotriazole70%92%X
Carbamazepine57%83%XX
Sulfamethoxazole13%83%XX
Fluconazole40%83%X
Lidocaine20%83%XX
Hexamethylenetetramine 37%92%
Bupropion13%75%XX 4
Gabapentin7%83% X
Tramadol10%75%XX
Meprobamate7%83%XX
Metformin3%33%XX
1 Reference [4]. 2 Reference [2]. 3 One blank sample had a reported concentration above the laboratory reporting limit (LRL). The concentration in the blank sample was 5% of the average environmental concentration, indicating little possible bias from blank samples for this compound. 4 Specifically noted presence of metabolite of this compound, not the compound itself.
Table 2. Results of apparent-age calculations, concentrations of atmospheric tracer gases, and idealized piston-flow ages of groundwater wells in the Fountain Creek alluvial aquifer.
Table 2. Results of apparent-age calculations, concentrations of atmospheric tracer gases, and idealized piston-flow ages of groundwater wells in the Fountain Creek alluvial aquifer.
USGS Site IDSite AliasR/Ra3H (TU)3H-3He Age (Years) 1Fpost-1953SF6 Concentration (fmol/kg)SF6
Concentration (pptv) 1
SF6 Piston-Flow Age
(Years) 2,3
38482410440510103-0021.0116.793.30.82.568.035.4
38471010443120104-0090.9423.1100.362.77.267.9
384233104425801A10.9486.0700.7235.4110.4MF
3847181044637013DAA--8.95--1.0628.0890.79MF
384509104435901CO259-250.7464.9110.80.5810.4122.41MF
384719104444701CO259-260.9985.81.20.682.749.131.9
384437104422601DBA--5.8--0.68------
384112104421301FTN10.9955.911.10.73.8211.9MF
384639104461401BAC1--8.49--1.0183.07196.38MF
384949104424501MW1-11.0016.184.60.732.718.134.9
384848104413901MW2-41.0214.214.90.492.286.849.9
384408104424701NONUM-21.0154.952.70.584.688.652.9
384617104455901CAD10.84419.93.22.38369.71714.58MF
384929104431101T01-MW002--8.33--0.993.547.885.9
384956104422801T02-MW006--4.58--0.54------
384917104422701T04-MW004--4.94--0.582.315.5814.4
384818104415701T07-MW004--4.22--0.493.487.656.9
384758104422301T07-MW006--4.31--0.52.976.819.9
384732104430901T13-MW0040.9824.623.20.542.576.3311.4
384648104454501TH-220.8955.321.10.633.6912.25MF
384540104453601TH-46--4.79--0.56219.2616.75MF
384503104451601TH-50.9985.270.60.6225.2862.35MF
384636104465401TH-520.976.480.70.7713.5751.94MF
384534104450302U-11--4.52--0.533.588.613.4
384513104445302U-120.8084.965.60.5833.8472.45MF
384433104440701U-14B0.8914.576.90.5498.94248.12MF
384420104432601U-150.9684.950.60.588.9331.76MF
384652104465101U-70.9688.9801.077.7820.17MF
384604104451502U-9--6.66--0.793.548.463.9
384217104402901JCC0.8025.1621.50.612.666.0512.4
-- Indicates no applicable data. USGS, US Geological Survey. TU, tritium units. 1 Calculation used components that were corrected for excess air and recharge temperature. 2 “MF” indicates a failure of the piston-flow model to explain the data because SF6 concentration is greater than that explainable by the piston-flow model. 3 When available noble gases were used to compute excess air and recharge temperature, which was used in calculating piston-flow age.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Newman, C.P.; Paschke, S.S.; Keith, G. Natural and Anthropogenic Geochemical Tracers to Investigate Residence Times and Groundwater–Surface-Water Interactions in an Urban Alluvial Aquifer. Water 2021, 13, 871. https://doi.org/10.3390/w13060871

AMA Style

Newman CP, Paschke SS, Keith G. Natural and Anthropogenic Geochemical Tracers to Investigate Residence Times and Groundwater–Surface-Water Interactions in an Urban Alluvial Aquifer. Water. 2021; 13(6):871. https://doi.org/10.3390/w13060871

Chicago/Turabian Style

Newman, Connor P., Suzanne S. Paschke, and Gabrielle Keith. 2021. "Natural and Anthropogenic Geochemical Tracers to Investigate Residence Times and Groundwater–Surface-Water Interactions in an Urban Alluvial Aquifer" Water 13, no. 6: 871. https://doi.org/10.3390/w13060871

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop