Heat as a ground water tracer

Heat carried by ground water serves as a tracer to identify surface water infiltration, flow through fractures, and flow patterns in ground water basins. Temperature measurements can be analyzed for recharge and discharge rates, the effects of surface warming, interchange with surface water, hydraulic conductivity of streambed sediments, and basin-scale permeability. Temperature data are also used in formal solutions of the inverse problem to estimate ground water flow and hydraulic conductivity. The fundamentals of using heat as a ground water tracer were published in the 1960s, but recent work has significantly expanded the application to a variety of hydrogeological settings. In recent work, temperature is used to delineate flows in the hyporheic zone, estimate submarine ground water discharge and depth to the salt-water interface, and in parameter estimation with coupled ground water and heat-flow models. While short reviews of selected work on heat as a ground water tracer can be found in a number of research papers, there is no critical synthesis of the larger body of work found in the hydrogeological literature. The purpose of this review paper is to fill that void and to show that ground water temperature data and associated analytical tools are currently underused and have not yet realized their full potential.     

Variability and comparison of hyporheic water temperatures and seepage fluxes in a small Atlantic salmon stream

Ground water discharge is often a significant factor in the quality of fish spawning and rearing habitat and for highly biologically productive streams. In the present study, water temperatures (stream and hyporheic) and seepage fluxes were used to characterize shallow ground water discharge and recharge within thestreambed of Catamaran Brook, a small Atlantic salmon (Salmo salar) stream in central New Brunswick, Canada. Three study sites were instrumented using a total of 10 temperature sensors and 18 seepage meters. Highly variable mean seepage fluxes, ranging from 1.7 x 10(-4) to 2.5 cm3 m(-2) sec(-1), and mean hyporheic water temperatures, ranging from 10.5 degrees to 18.0 degrees C, at depths of 20 to 30 cm in the streambed were dependent on streambed location (left versus right stream bank and site location) and time during the summer sampling season. Temperature data were usefulfor determining if an area of the streambed was under discharge (positive flux), recharge (negative flux), or parallel flow (no flux) conditions and seepage meters were used to directly measure the quantity of water flux. Hyporheic water temperature measurements and specific conductance measurements of the seepage meter sample water, mean values ranging from 68.8 to 157.9 microS/cm, provided additional data for determining flux sources. Three stream banks were consistently under discharge conditions, while the other three stream banks showed reversal from discharge to recharge conditions over the sampling season. Results indicate that the majority of the water collected in the seepage meters was composed of surface water. The data obtained suggests that even though a positive seepage flux is often interpreted as ground water discharge, this discharging water may be of stream water origin that has recently entered the hyporheic zone.The measurement of seepage flux in conjunction with hyporheic water temperature or other indicators of water origin should be considered when attempting to quantify the magnitude of exchange and the source of hyporheic water.     

Delineating and quantifying ground water discharge zones using streambed temperatures

Streambed temperature mapping, hydraulic testing using minipiezometers, and geochemical analyses of interstitial water of the streambed were used to delineate the pattern of ground water discharge in a sandy streambed and to develop a flux-based conceptual model for ground water/surface water interactions. A new and simple empirical method was used to relate fluxes obtained from minipiezometer data to streambed temperatures. The relationship allowed flux to be calculated at locations where only streambed temperature measurements were made. Slug testing and potentiomanometer measurements at 34 piezometers indicated ground water discharge ranged from 0.03 to 446 L/m2/day (and possibly as high as 7060 L/m2/day) along a 60 m long by 11 to 14 m wide reach of river. Complex but similar plan-view patterns of flux were calculated for both summer and winter using hundreds of streambed temperatures measured on a 1 by 2 m grid. The reach was dominated by ground water discharge and 5% to 7% of the area accounted for approximately 20% to 24% of the total discharge. < 12% of the total area consisted of recharge zones or no-discharge zones. A conceptual model for ground water/surface water interactions consisting of five different behaviors was developed based on the magnitude and direction of flux across the surface of the streambed. The behaviors include short-circuit discharge (e.g., high-flow springs), high discharge (e.g., preferential flowpaths), low to moderate discharge, no discharge (e.g., horizontal hyporheic or ground water flow), and recharge. Geological variations at depth played a key role in determining which type of flow behavior occurred in the streambed.     

Modelling discharge through artesian springs based on a high‐resolution piezometric network

Artesian springs are localized aquifer outlets that originate when pressurized ground water is allowed to rise to the surface. Computing artesian discharge directly is often subject to practical difficulties such as restricted accessibility, abundant vegetation or slow flow rates. These circumstances call for indirect approaches to quantify flow. This paper presents a method to estimate ground water discharge through an upwelling spring by means of a three‐layer steady‐state groundwater flow model. Model inputs include on‐site measurements of vertical sediment permeability, sediment temperatures and hydraulic gradients. About 70 spring bed piezometers were used to carry out permeability tests within the spring sediments, as well as to quantify the hydraulic head at different depths below the discharge point. Sediment temperatures were measured at different depths and correlated to permeabilities in order to demonstrate the potential of temperature as a substitute for cumbersome slug tests. Results show that the spatial distribution of discharge through the spring bottom is highly heterogeneous, as sediment permeability varies by several orders of magnitude within centimetres. Sensitivity analyses imply that geostatistical interpolation is irrelevant to the results if field datasets come from a sufficiently high resolution of piezometric records. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling decadal timescale interactions between surface water and ground water in the central Everglades, Florida, USA

Surface-water and ground-water flow are coupled in the central Everglades, although the remoteness of this system has hindered many previous attempts to quantify interactions between surface water and ground water. We modeled flow through a 43,000 ha basin in the central Everglades called Water Conservation Area 2A. The purpose of the model was to quantify recharge and discharge in the basin's vast interior areas. The presence and distribution of tritium in ground water was the principal constraint on the modeling, based on measurements in 25 research wells ranging in depth from 2 to 37 m. In addition to average characteristics of surface-water flow, the model parameters included depth of the layer of ‘interactive’ ground water that is actively exchanged with surface water, average residence time of interactive ground water, and the associated recharge and discharge fluxes across the wetland ground surface. Results indicated that only a relatively thin (8 m) layer of the 60 m deep surfical aquifer actively exchanges surface water and ground water on a decadal timescale. The calculated storage depth of interactive ground water was 3.1 m after adjustment for the porosity of peat and sandy limestone. Modeling of the tritium data yielded an average residence time of 90 years in interactive ground water, with associated recharge and discharge fluxes equal to 0.01 cm d⁻¹. ³H/³He isotopic ratio measurements (which correct for effects of vertical mixing in the aquifer with deeper, tritium-dead water) were available from several wells, and these indicated an average residence time of 25 years, suggesting that residence time was overestimated using tritium measurements alone. Indeed, both residence time and storage depth would be expected to be overestimated due to vertical mixing. The estimate of recharge and discharge (0.01 cm d⁻¹) that resulted from tritium modeling therefore is still considered reliable, because the ratio of residence time and storage depth (used to calculated recharge and discharge) is much less sensitive to vertical mixing compared with residence time alone. We conclude that a small but potentially significant component of flow through the Everglades is recharged to the aquifer and stored there for years to decades before discharged back to surface water. Long-term storage of water and solutes in the ground-water system beneath the wetlands has implications for restoration of Everglades water quality.     

A digital procedure for ground water recharge and discharge pattern recognition and rate estimation

A digital procedure to estimate recharge/discharge rates that requires relatively short preparation time and uses readily available data was applied to a setting in central Wisconsin. The method requires only measurements of the water table, fluxes such as stream baseflows, bottom of the system, and hydraulic conductivity to delineate approximate recharge/discharge zones and to estimate rates. The method uses interpolation of the water table surface, recharge/discharge mapping, pattern recognition, and a parameter estimation model. The surface interpolator used is based on the theory of radial basis functions with thin-plate splines. The recharge/discharge mapping is based on a mass-balance calculation performed using MODFLOW. The results of the recharge/discharge mapping are critically dependent on the accuracy of the water table interpolation and the accuracy and number of water table measurements. The recharge pattern recognition is performed with the help of a graphical user interface (GUI) program based on several algorithms used in image processing. Pattern recognition is needed to identify the recharge/discharge zonations and zone the results of the mapping method. The parameter estimation program UCODE calculates the parameter values that provide a best fit between simulated heads and flows and calibration head-and-flow targets. A model of the Buena Vista Ground Water Basin in the Central Sand Plains of Wisconsin is used to demonstrate the procedure.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.     

A first estimate of ground water ages for the deep aquifer of the Kathmandu Basin, Nepal, using the radioisotope chlorine-36

The Kathmandu Basin in Nepal contains up to 550 m of Pliocene-Quaternary fluvio-lacustrine sediments which have formed a dual aquifer system. The unconfined sand and gravel aquifer is separated by a clay aquitard, up to 200 m thick, from the deeper, confined aquifer, comprised of Pliocene sand and gravel beds, intercalated with clay, peat, and lignite. The confined aquifer currently provides an important water supply to the central urban area but there are increasing concerns about its sustainability due to overexploitation. A limited number of determinations of the radioisotope 36Cl have been made on bore waters in the basin, allowing us to postulate on the age of ground water in the deeper, confined aquifer. Ground water evolution scenarios based on radioisotope decay, gradual dissolution of formational salts as the ground waters move downgradient, and flow velocity estimations produce comparable ground water ages for the deep waters, ranging from 200,000 to 400,000 years. From these ages, we deduce a mean ground water flow velocity of only 45 mm/year from recharge in the northeast to the main extraction region 15 km to the southwest. We thus estimate current recharge at about 5 to 15 mm/year, contributing 40,000 to 1.2 million m3/year to the ground water system. Current ground water extraction is estimated to be 20 times this amount. The low specific discharge confirms that the resource is being mined, and, based on current projections, reserves will be used up within 100 years.     

Simulation of fluid, heat transport to estimate desert stream infiltration

In semiarid regions, the contribution of infiltration from intermittent streamflow to ground water recharge may be quantified by comparing simulations of fluid and heat transport beneath stream channels to observed ground temperatures. In addition to quantifying natural recharge, streamflow infiltration estimates provide a means to characterize the physical properties of stream channel sediments and to identify suitable locations for artificial recharge sites. Rates of winter streamflow infiltration along stream channels are estimated based on the cooling effect of infiltrated water on streambed sediments, combined with the simulation of two-dimensional fluid and heat transport using the computer program VS2DH. The cooling effect of ground water is determined by measuring ground temperatures at regular intervals beneath stream channels and nearby channel banks in order to calculate temperature-depth profiles. Additional data inputs included the physical, hydraulic, and thermal properties of unsaturated alluvium, and monthly ground temperatures measurements over an annual cycle. Observed temperatures and simulation results can provide estimates of the minimum threshold for deep infiltration, the variability of infiltration along stream channels, and also the frequency of infiltration events.     

On the interpretation of recharge estimates from steady-state model calibrations

Ground water recharge is often estimated through the calibration of ground water flow models. We examine the nature of calibration errors by considering some simple mathematical and numerical calculations. From these calculations, we conclude that calibrating a steady-state ground water flow model to water level extremes yields estimates of recharge that have the same value as the time-varying recharge at the time the water levels are measured. These recharge values, however, are a subdued version of the actual transient recharge signal. In addition, calibrating a steady-state ground water flow model to data collected during periods of rising water levels will produce recharge values that underestimate the actual transient recharge. Similarly, calibrating during periods of falling water levels will overestimate the actual transient recharge. We also demonstrate that average water levels can be used to estimate the actual average recharge rate provided that water level data have been collected for a sufficient amount of time.     

Estimating deep recharge rates beneath an interlobate moraine using temperature logs

The Sandilands area of southeastern Manitoba contains an interlobate moraine that is a major ground water recharge area. Underlying the highly permeable sediments of the moraine are up to 100 m of till and the subcrop of the Winnipeg Formation, which contains a major sandstone aquifer. Ground water flow within the till is examined using high-resolution temperature profiles and solutions to the differential equation for heat flow in porous media. These analyses indicate that recharge to the sandstone aquifer is occurring at a rate of approximately 2 x 10(-8) m/sec beneath the moraine, which is in agreement with recharge rates determined by conventional ground water hydraulics (10(-7) to 10(-10)(m/sec) and another study using multiple environmental tracers (1 x 10(-9) to 6 X 10(-9) m/sec). The use of temperature to determine ground water flux is not limited by half-lives as many environmental tracers are, and this allows for cost-effective estimation of recharge and discharge rates over longer periods.     

The water budget myth revisited: why hydrogeologists model

Within the ground water community, the idea persists that if one can estimate the recharge to a ground water system, one then can determine the size of a sustainable development. Theis addressed this idea in 1940 and showed it to be wrong-yet the myth continues. The size of a sustainable ground water development usually depends on how much of the discharge from the system can be "captured" by the development. Capture is independent of the recharge; it depends on the dynamic response of the aquifer system to the development. Ground water models were created to study the response dynamics of ground water systems; it is one of the principal reasons hydrogeologists model.     

Effects of afforestation on groundwater recharge and water budgets in the western region of Uruguay

Uruguay has stimulated the development of its forest sector since the promulgation of Forest Law N° 15 939 in December of 1987. Nevertheless, the substitution of natural grasslands with forest plantations for industrial use has raised concerns regarding hydrological processes of groundwater recharge and water consumption involving evapotranspiration. The purpose of this study is to assess the effects of this substitution approach on water resources. Input data were collected from two small experimental watersheds of roughly 100–200 hectares located in western Uruguay. The watersheds are characterized by Eucalyptus Globulus ssp. Maidenni and natural grasslands for cattle use. Total rainfall, stream discharge, rainfall redistribution, soil water content and groundwater level data were collected. Groundwater recharge was estimated from water table fluctuations and from groundwater contributions to base flows. Seasonal and annual water budgets were computed from October of 2006 to September of 2014 to evaluate changes in the hydrological processes. The data show a decrease in annual specific discharge of roughly 17% for mean hydrological years and no conclusive effects on annual groundwater recharge in the forested watershed relative to the reference pasture watershed. Reduced annual specific discharge is equivalent to the mean annual interception. The computed actual annual evapotranspiration is consistent with international catchment measurements. Reduction rates vary seasonally and according to accumulated rainfall and its temporary distribution. The degree of specific discharge decline is particularly high for drier autumns and winters (32 to 28%) when the corresponding rainfall varies from 275 to 400 mm. These results are of relevance for water resources management efforts, as water uses downstream can be affected. These findings, based on a study period dominated by anomalous wet springs and summers and by dry autumns and winters, oppose earlier results based on 34 years of rainfall and discharge data drawn from Uruguayan large basins. Copyright © 2016 John Wiley & Sons, Ltd.     

Recharge estimation for transient ground water modeling

Reliable ground water models require both an accurate physical representation of the system and appropriate boundary conditions. While physical attributes are generally considered static, boundary conditions, such as ground water recharge rates, can be highly variable in both space and time. A practical methodology incorporating the hydrologic model HELP3 in conjunction with a geographic information system was developed to generate a physically based and highly detailed recharge boundary condition for ground water modeling. The approach uses daily precipitation and temperature records in addition to land use/land cover and soils data. The importance of the method in transient ground water modeling is demonstrated by applying it to a MODFLOW modeling study in New Jersey. In addition to improved model calibration, the results from the study clearly indicate the importance of using a physically based and highly detailed recharge boundary condition in ground water quality modeling, where the detailed knowledge of the evolution of the ground water flowpaths is imperative. The simulated water table is within 0.5 m of the observed values using the method, while the water levels can differ by as much as 2 m using uniform recharge conditions. The results also show that the combination of temperature and precipitation plays an important role in the amount and timing of recharge in cooler climates. A sensitivity analysis further reveals that increasing the leaf area index, the evaporative zone depth, or the curve number in the model will result in decreased recharge rates over time, with the curve number having the greatest impact.     

Measuring Artificial Recharge with Fiber Optic Distributed Temperature Sensing

Heat was used as a tracer to measure infiltration rates from a recharge basin. The propagation of diurnal oscillation of surface water temperature into the basin bed was monitored along a transect using Fiber Optic Distributed Temperature Sensing (FODTS). The propagation rate was related to downward specific discharge using standard theory of heat advection and dispersion in saturated porous media. An estimate of the temporal variation of heat propagation was achieved using a wavelet transform to find the phase lag between the surface temperature diurnal oscillation and the correlated oscillation at 0.33 and 0.98 m below the bed surface. The wavelet results compared well to a constant velocity model of thermal advection and dispersion during periods of relatively constant discharge rates. The apparent dispersion of heat was found to be due primarily to hydrodynamic mechanisms rather than thermal diffusion. Specific discharge estimates using the FODTS technique also compared well to water balance estimates over a four month period, although there were occasional deviations that have yet to be adequately explained. The FODTS technique is superior to water balance in that it produces estimates of infiltration rate every meter along the cable transect, every half hour. These high resolution measurements highlighted areas of low infiltration and demonstrated the degradation of basin efficiency due to source waters of high suspended solids. FODTS monitoring promises to be a useful tool for diagnosing basin performance in an era of increasing groundwater demand.     

Impact of Recharge Through Residual Oil Upon Sampling of Underlying Ground Water

At an aviation gasoline spill site in Traverse City, Michigan, historical records indicate a positive correlation between significant rainfall events and increased concentrations of slightly soluble organic compounds in the monitoring wells of the site. To investigate the recharge effect on ground water quality due to infiltrating, water percolating past residual oil and into the saturated zone, an in situ infiltration experiment was performed at the site. Sampling cones were set at various depths below a circular test area, 13 feet (4 meters) in diameter. Rainfall was simulated by sprinkling the test area at a rate sufficiently low to prevent runoff. The sampling cones for soil-gas and ground water quality were installed in the unsaturated and saturated zones to observe the effects of the recharge process. At the time of the test, the water table was below the residual oil layer. The responses of the soil-gas and ground water quality were monitored during the recharge and drainage periods, which resulted from the sprinkling.
Infiltrated water was determined to have transported organic constituents of the residual oil, specifically benzene, toluene, ethylbenzene, and ortho-xylene (BTEX), into the ground water beneath the water table, elevating the aqueous concentrations of these constituents in the saturated zone. Soil-gas concentrations of the organic compounds in the unsaturated zone increased with depth and time after the commencement of infiltration. Reaeration of the unconfined aquifer via the infiltrated water was observed. It is concluded that water quality measurements are directly coupled to recharge events for the sandy type of aquifer with an overlying oil phase, which was studied in this work. Ground water sampling strategies and data analysis need to reflect the effect of recharge from precipitation on shallow, unconfined aquifers where an oil phase may be present.     

Spatial and temporal variability of ground water recharge in central Australia: a tracer approach

Two environmental tracer methods are applied to the Ti-Tree Basin in central Australia to shed light on the importance of recharge from floodouts of ephemeral rivers in this arid environment. Ground water carbon-14 concentrations from boreholes are used to estimate the average recharge rate over the interval between where the ground water sample first entered the saturated zone and the bore. Environmental chloride concentrations in ground water samples provide estimates of the recharge rate at the exact point in the landscape where the sample entered the saturated zone. The results of the two tracer approaches indicate that recharge rates around one of the rivers and an extensive floodplain are generally higher than rates of diffuse recharge that occurs in areas of lower topographic relief. Ground water 2H/1H and 18O/16O compositions are all depleted in the heavier isotopes (delta2H = -67 per thousand to -50 per thousand; delta18O = -9.2 per thousand to -5.7%o) compared with the long-term, amount-weighted mean isotopic composition of rainfall in the area (delta2H = -33.8 per thousand; delta18O = -6.3 per thousand). This indicates that recharge throughout the basin occurs only after intense rainfall events of at least 150 to 200 mm/month. Finally, a recharge map is developed to highlight the spatial extent of the two recharge mechanisms. Floodout recharge to the freshest ground water (TDS <1,000 mg/L) is approximately 1.9 mm/year compared with a mean recharge rate of approximately 0.2 mm/year to the remainder of the basin. These findings have important implications for management of the ground water resource.     

Use of temperature profiles and stable isotopes to trace flow lines: Nagaoka area, Japan

In this study, we use borehole temperature data and stable isotopes to delineate the flow system and estimate the effect of urbanization in the Nagaoka area of Japan. Temperature profiles were measured four times in observation wells during the period 2000-2001 and compared with those measured in the same wells during the period 1977-1983 (Taniguchi 1986). Water was sampled in both observation and pumping wells during the same period. The temporal and spatial variability in temperature indicate clearly the effect of urban warming and heavy pumping on the ground water system. Urban warming caused higher temperatures recently as compared to the older values, and pumping caused induced recharge from the river to the ground water. The stable isotope data show the ground water flow system is divided into shallow, intermediate, and deep systems, and that land use and infiltration rate are affecting the shallow flow system.     

Estimating ground water recharge using flow models of perched karstic aquifers

The fraction of rain that is annually recharged to ground water is a function of the transient quantities of precipitation (wet vs. dry years) as well as other meteorological and geologic factors, and thus it is very difficult to estimate. In this study, we have used long records (20 to 30 years) of precipitation and spring discharge to reconstruct the transient character of yearly recharge. These data sets were used to calibrate numerical ground water flow models on the less than 3 km(2) scale for four separate perched karstic aquifers in the Judean and Samarian Mountains of Israel. The stratification and karstic character of the local carbonate rock aquifers cause ground water to flow through discrete dissolution channels and to discharge at isolated springs. An innovative, dual-porosity approach was used where a finite-difference solution simulates flow in the rock matrix, while the karstic channels are simulated using computationally simple drains. Perched conditions are also simulated innovatively using MODFLOW by treating the bottom unsaturated layer as if it is saturated, but by assuming zero pressure head throughout the "unsaturated" layer. Best fitting between measured and computed spring hydrograph data has allowed us to develop a set of empirical functions relating measured precipitation to recharge to the aquifer. The generic methodology presented gives insight into the suspected changes in aquifer recharge rates between particularly wet or dry years.