Propagation of seasonal temperature signals into an aquifer upon bank infiltration

Acesulfame allows the tracing of multiple sources of wastewater and riverbank filtration

Aquifers providing drinking water are increasingly threatened by emerging contaminants due to wastewater inputs from multiple sources. These inputs have to be identified, differentiated, and characterized to allow an accurate risk assessment and thus ensure the safety of drinking water through appropriate management. We hypothesize, that in climates with seasonal temperature variations, the sweetener acesulfame potassium (ACE) provides new pathways to study wastewater inputs to aquifers. Specifically, this study investigates the temperature-driven seasonal oscillation of ACE to assess multiple sources of wastewater inputs at a riverbank filtration site. ACE concentrations in the river water varied from 0.2 to 1 μg L^-1 in the cold season (T < 10 °C) to 0-0.1 μg L^-1 in the warm season (T > 10 °C), due to temperature-dependent biodegradation during wastewater treatment. This oscillating signal could be traced throughout the aquifer over distances up to 3250 m from two different infiltration sources. A transient numerical model of groundwater flow and ACE transport was calibrated over hydraulic heads and ACE concentrations, allowing the accurate calculation of mixing ratios, travel times, and flow-path directions for each of the two infiltration sources. The calculated travel time from the distant infiltration source was of 67 days, while that from the near source was of 20 days. The difference in travel times leads to different potential degradation of contaminants flowing into the aquifer from the river, thus demonstrating the importance of individually assessing the locations of riverbank infiltration. The calibrated ACE transport model allowed calculating transient mixing ratios, which confirmed the impact of river stage and groundwater levels on the mixing ratio of the original groundwater and the bank filtrate. Therefore, continuous monitoring of ACE concentrations can help to optimize the management of the water works with the aim to avoid collection of water with very short travel times, which has important regulative aspects. Our findings demonstrate the suitability of ACE as a transient tracer for identifying multiple sources of wastewater, including riverbank filtration sites affected by wastewater treatment plant effluents. ACE seasonal oscillation tracking thus provides a new tool to be used in climates with pronounced seasonal temperature variations to assess the origins of contamination in aquifers, with time and cost advantages over multi-tracer approaches.

Seasonal biodegradation of the artificial sweetener acesulfame enhances its use as a transient wastewater tracer

The persistence of the artificial sweetener acesulfame potassium (ACE) during wastewater treatment and subsequently in the aquatic environment has made it a widely used tracer of wastewater inputs to both surface water and groundwater. However, the recently observed biodegradation of ACE during wastewater treatment has questioned the validity of this application. In this study, we assessed the use of ACE not only as a marker of wastewater, but also as a transient wastewater tracer that allows both the calculation of mixing ratios and travel times through the aquifer as well as the calibration of transient groundwater flow and mass transport models. Our analysis was based on data obtained in a nearly 8-year river water and groundwater sampling campaign along a confirmed wastewater-receiving riverbank filtration site located close to a drinking water supply system. We provide evidence that temperature controls ACE concentration and thus its seasonal oscillation. River water data showed that ACE loads decreased from 1.5-4 mg·s^-1 in the cold season (December to June; T<10 °C) to 0-0.5 mg·s^-1 in the warm season (July to November; T>10 °C). This seasonal variability of >600% was detectable in the aquifer and preserved >3 km, with ACE concentrations oscillating between <LOQ in the warm season up to 1 μg·L^-1 in the cold season. The large seasonal variation in ACE concentrations during wastewater treatment, compared to the other sweeteners (sucralose, cyclamate, and saccharin) and chloride enables its use as a transient tracer of wastewater inflows and riverbank filtration. In addition, the arrival time of the ACE concentration peak can be used to estimate groundwater flow velocity and mixing ratios, thereby demonstrating its potential in the calibration of groundwater numerical models.

Thermal influences on groundwater in urban environments - A multivariate statistical analysis of the subsurface heat island effect in Munich

Shallow aquifers beneath cities are highly influenced by anthropogenic heat sources, resulting in the formation of extensive subsurface urban heat islands. In addition to anthropogenic factors, natural factors also influence the subsurface temperature. However, the effect of individual factors is difficult to capture due to high temporal dynamics in urban environments. Particularly in the case of shallow aquifers, seasonal temperature fluctuations often override the influence of existing heat sources or sinks. For the city of Munich, we identify the dominant anthropogenic and natural influences on groundwater temperature and analyse how the influences change with increasing depth in the subsurface. For this purpose, we use depth temperature profiles from 752 selected groundwater monitoring wells. Since the measurements were taken at different times, we developed a statistical approach to compensate for the different seasonal temperature influences using passive heat tracing. Further, we propose an indicator to spatially assess the thermal stress on the aquifer. A multiple regression analysis of four natural and nine anthropogenic factors identified surface sealing as the strongest and the district heating grid as a weak but significant warming influence. The natural factors, aquifer thickness, depth-to-water and Darcy velocity show a significant cooling influence on the groundwater temperature. In addition, we show that local drivers, like thermal groundwater uses, surface waters and underground structures do not significantly contribute to the city-wide temperature distribution. The subsequent depth-dependent analysis revealed that the influence of aquifer thickness and depth-to-water increases with depth, whereas the influence of Darcy velocity decreases, and surface sealing and the heating grid remain relatively constant. In conclusion, this study shows that the most critical factor for the mitigation of future groundwater warming in cities is to minimize further sealing of the ground, to restore the permeability of heavily sealed areas and to retain open landscapes.

Estimating Travel Time in Bank Filtration Systems from a Numerical Model Based on DTS Measurements

An approach is presented to determine the seasonal variations in travel time in a bank filtration system using a passive heat tracer test. The temperature in the aquifer varies seasonally because of temperature variations of the infiltrating surface water and at the soil surface. Temperature was measured with distributed temperature sensing along fiber optic cables that were inserted vertically into the aquifer with direct push equipment. The approach was applied to a bank filtration system consisting of a sequence of alternating, elongated recharge basins and rows of recovery wells. A SEAWAT model was developed to simulate coupled flow and heat transport. The model of a two-dimensional vertical cross section is able to simulate the temperature of the water at the well and the measured vertical temperature profiles reasonably well. MODPATH was used to compute flowpaths and the travel time distribution. At the study site, temporal variation of the pumping discharge was the dominant factor influencing the travel time distribution. For an equivalent system with a constant pumping rate, variations in the travel time distribution are caused by variations in the temperature-dependent viscosity. As a result, travel times increase in the winter, when a larger fraction of the water travels through the warmer, lower part of the aquifer, and decrease in the summer, when the upper part of the aquifer is warmer.

Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution, spreadsheet tool, and field applications

Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction-dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson (2004), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater-surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux-LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi-decadal surface signals have disturbed subsurface thermal regimes.

Tracer tomography: design concepts and field experiments using heat as a tracer

Numerical and laboratory studies have provided evidence that combining hydraulic tomography with tomographic tracer tests could improve the estimation of hydraulic conductivity compared with using hydraulic data alone. Field demonstrations, however, have been lacking so far, which we attribute to experimental difficulties. In this study, we present a conceptual design and experimental applications of tracer tomography at the field scale using heat as a tracer. In our experimental design, we improve active heat tracer testing by minimizing possible effects of heat losses, buoyancy, viscosity, and changing boundary conditions. We also utilize a cost-effective approach of measuring temperature changes in situ at high resolution. We apply the presented method to the 8 m thick heterogeneous, sandy gravel, alluvial aquifer at the Lauswiesen Hydrogeological Research Site in Tübingen, Germany. Results of our tomographic heat-tracer experiments are in line with earlier work on characterizing the aquifer at the test site. We demonstrate from the experimental perspective that tracer tomography is applicable and suitable at the field scale using heat as a tracer. The experimental results also demonstrate the potential of heat-tracer tomography as a cost-effective means for characterizing aquifer heterogeneity.

Argon concentration time-series as a tool to study gas dynamics in the hyporheic zone

The oxygen dynamics in the hyporheic zone of a peri-alpine river (Thur, Switzerland), were studied through recording and analyzing the concentration time-series of dissolved argon, oxygen, carbon dioxide, and temperature during low flow conditions, for a period of one week. The argon concentration time-series was used to investigate the physical gas dynamics in the hyporheic zone. Differences in the transport behavior of heat and gas were determined by comparing the diel temperature evolution of groundwater to the measured concentration of dissolved argon. These differences were most likely caused by vertical heat transport which influenced the local groundwater temperature. The argon concentration time-series were also used to estimate travel times by cross correlating argon concentrations in the groundwater with argon concentrations in the river. The information gained from quantifying the physical gas transport was used to estimate the oxygen turnover in groundwater after water recharge. The resulting oxygen turnover showed strong diel variations, which correlated with the water temperature during groundwater recharge. Hence, the variation in the consumption rate was most likely caused by the temperature dependence of microbial activity.