The Hydraulic Evolution of Groundwater-Fed Pit Lakes After Mine Closure.

Open pit mining frequently requires regional water tables to be lowered to access ore deposits. When mines close, dewatering ceases allowing the water table to recover. In arid and semi-arid mining regions, the developing pit lakes are predominantly fed by groundwater during this recovery phase and pit lakes develop first into "terminal sinks" for the surrounding groundwater system. With time, the re-establishment of regional hydraulic gradients can cause pit lakes to develop into throughflow systems, in which pit lake water outflows into adjacent aquifers. In this study, we use numerical groundwater modeling to aid process understanding of how regional hydraulic gradients, aquifer properties, net evaporation rates, and pit geometry determine the hydraulic evolution of groundwater-fed pit lakes. We find that before the recovery of the regional water table to its new equilibrium, pit lakes frequently transition to throughflow systems. Throughflow from pit lakes to downstream aquifers can develop within two decades following cessation of dewatering even under low hydraulic gradients (e.g., 5 × 10-4) or high net evaporation rates (e.g., 2.5 m/year). Pit lakes remain terminal sinks only under suitable combinations of high evaporation rates, low hydraulic gradients, and low hydraulic conductivities. In addition, we develop an approximate analytical solution for a rapid assessment of the hydraulic status of pit lakes under steady-state conditions. Understanding whether pit lakes remain terminal sinks or transition into throughflow systems largely determines the long-term water quality of pit lakes and downstream aquifers. This knowledge is fundamental for mine closure and planning post-mining land use.

Managed Aquifer Recharge in Mining: A Review.

Managed aquifer recharge (MAR) has been gaining adoption within the mining industry for managing surplus water volumes and reducing the groundwater impacts of dewatering. This paper reviews MAR for mining and includes an inventory of 27 mines using or considering MAR for current or future operations. Most mines using MAR are in arid or semi-arid regions and are implementing it through infiltration basins or bore injection to manage surplus water, preserve aquifers for environmental or human benefit, or adhere to licensing that requires zero surface discharge. Surplus water volumes, hydrogeological conditions, and economics play a pivotal role in the feasibility of MAR for mining. Groundwater mounding, well clogging, and interaction between adjacent mines are common challenges. Mitigation strategies include predictive groundwater modeling, extensive monitoring programs, rotation of infiltration or injection facilities, physical and chemical treatments for clogging, and careful location for MAR facilities in relation to adjacent operations. Should water availability alternate between shortage and excess, injection bores may be used for supply, thus reducing costs and risks associated with drilling new wells. MAR, if applied strategically, also has the potential to accelerate groundwater recovery post-mine closure. The success of MAR for mining is emphasized by mines opting to increase MAR capacity alongside dewatering expansions, as well as prospective mines proposing MAR for future water requirements. Upfront planning is the key to maximizing MAR benefits. Improved information sharing could help increase awareness and uptake of MAR as an effective and sustainable mine water management tool.

Facilitating Open Pit Mine Closure with Managed Aquifer Recharge.

Dewatering of open pit mines can lower the regional water table for distances of several kilometers from the pit. When the mine is closed, dewatering operations usually cease, and the water table near the pit begins to rise. If the pit is backfilled, the water table will eventually recover, but this recovery may take several hundred years. However, if the extracted water is re-injected into the subsurface, then this may accelerate recovery of the water table. We show that there is an optimal distance for re-injection, which is sufficiently far from the mine to minimize the amount of groundwater that flows back to the pit during mine operations (and hence necessitate additional pumping) but is still close enough to speed up the water table recovery post-mine closure. The optimal injection distance increases with the aquifer hydraulic diffusivity and the mine life (duration of dewatering and injection), and typically ranges between about two and nine times the radius of the mine pit. Where the mine pit is not backfilled, the relative reduction in drawdown due to injecting all the pumped water at the optimal distance is between approximately 10% and 50% after a recovery time equal to the mining period, increasing to 30% to 90% after a recovery time five times the mining period. The relative drawdown reduction due to managed aquifer recharge will be even greater for a pit which is backfilled when mining ceases.

The Effect of Undetected Barriers on Groundwater Drawdown and Recovery.

In large-scale pumping projects, such as mine dewatering, predictions are often made about the rate of groundwater level recovery after pumping has ceased. However, these predictions may be impacted by geological uncertainty-including the presence of undetected impermeable barriers. During pumping, an impermeable barrier may be undetected if it is located beyond the maximum extent of the cone of depression; yet it may still control drawdown during the recovery phase. This has implications for regional-scale modeling and monitoring of groundwater level recovery. In this article, non-dimensional solutions are developed to show the conditions under which a barrier may be undetected during pumping but still significantly impact groundwater level recovery. The magnitude of the impact from an undetected barrier will increase as the ratio of pumping rate to aquifer transmissivity increases. The results are exemplified for a hypothetical aquifer with an unknown barrier 3 km from a pumping well. The difference in drawdown between a model with and without a barrier may be <1 m in the 10 years while pumping is occurring, but up to 50 m after pumping has ceased.

A Novel Approach To Quantify Air-Water Gas Exchange in Shallow Surface Waters Using High-Resolution Time Series of Dissolved Atmospheric Gases.

Gas exchange across the air-water interface is a key process determining the release of greenhouse gases from surface waters and a fundamental component of gas dynamics in aquatic systems. To experimentally quantify the gas transfer velocity in a wide range of aquatic settings, a novel method based on recently developed techniques for the in situ, near-continuous measurement of dissolved (noble) gases with a field portable mass spectrometer is presented. Variations in observed dissolved gas concentrations are damped and lagged with respect to equilibrium concentrations, being the result of (a) temperature (and thus solubility) variations, (b) water depth, and (c) the specific gas transfer velocity ( ki). The method fits a model to the measured gas concentrations to derive the gas transfer velocity from the amplitude and the phase lag between observed and equilibrium concentrations. With the current experimental setup, the method is sensitive to gas transfer velocities of 0.05-9 m/day (for N2), at a water depth of 1 m, and a given daily water temperature variation of 10 °C. Experiments were performed (a) in a controlled experiment to prove the concept and to confirm the capability to determine low transfer velocities and (b) in a field study in a shallow coastal lagoon covering a range of transfer velocities, demonstrating the field applicability of the method.

Effects of Intraborehole Flow on Purging and Sampling Long-Screened or Open Wells.

Hydraulic head differences across the screened or open interval of a well significantly influence the sampled water mixture. Sample bias can occur due to an insufficient pumping rate and/or due to native groundwater displacement by intraborehole flow (IBF). Proper understanding of the sampled water mixture is crucial for accurate interpretation of environmental tracers and groundwater chemistry data, and hence groundwater characterization. This paper uses numerical modeling to quantify sample bias caused by IBF in an un-pumped high-yield well, and the influence of pumping rate and heterogeneity on the volume of pumpage required to purge an IBF plume. The results show that (1) the pumping rate must be at least an order of magnitude greater than the IBF rate to achieve permeability-weighted yield, (2) purge volume was 2.2 to 20.6 times larger than the IBF plume volume, with the ratio depending on plume location relative to hydraulic conductivity and head distributions, and (3) after an example 1000-day un-pumped period, purging required removal of at least three orders of magnitude more water than the common practice of three to five well volumes. These results highlight the importance of knowing the borehole flow regime to identify IBF inflow and outflow zones, estimate IBF rates, and to develop a strategic sampling approach.

Towards Quantifying the Likelihood of Water Resource Impacts from Unconventional Gas Development.

Gas production from unconventional reservoirs has led to widespread environmental concerns. Despite several excellent reviews of various potential impacts to water resources from unconventional gas production, no study has systematically and quantitatively assessed the potential for these impacts to occur. We use empirical evidence and numerical and analytical models to quantify the likelihood of surface water and groundwater contamination, and shallow aquifer depletion from unconventional gas developments. These likelihoods are not intended to be exact. They provide a starting point for comparing the probabilities of adverse impacts between types of water resources and pathways. This analysis provides much needed insight into what are "probable" rather than simply "possible" impacts. The results suggest that the most likely water resource impacts are surface water and groundwater contamination from spills at the well pad, which can be as high as 1 in 10 and 1 in 100 for each gas well, respectively. For wells that are hydraulically fractured, the likelihood of contamination due to inter-aquifer leakage is 1 in 106 or lower (dependent on the separation distance between the production formation and the aquifer). For gas-bearing formations that were initially over-pressurized, the potential for contamination from inter-aquifer leakage after production ceases could be as high as 1 in 400 where the separation between gas formation and shallow aquifer is 500 m, but will be much lower for greater separation distances (more characteristic of shale gas).

Groundwater-driven nutrient inputs to coastal lagoons: The relevance of lagoon water recirculation as a conveyor of dissolved nutrients.

Evaluating the sources of nutrient inputs to coastal lagoons is required to understand the functioning of these ecosystems and their vulnerability to eutrophication. Whereas terrestrial groundwater processes are increasingly recognized as relevant sources of nutrients to coastal lagoons, there are still limited studies evaluating separately nutrient fluxes driven by terrestrial groundwater discharge and lagoon water recirculation through sediments. In this study, we assess the relative significance of these sources in conveying dissolved inorganic nutrients (NO3-, NH4+ and PO43-) to a coastal lagoon (La Palme lagoon; France, Mediterranean Sea) using concurrent water and radon mass balances. The recirculation of lagoon water through sediments represents a source of NH4+ (1900-5500 mol d-1) and PO43- (22-71 mol d-1), but acts as a sink of NO3-. Estimated karstic groundwater-driven inputs of NO3-, NH4+ and PO43- to the lagoon are on the order of 200-1200, 1-12 and 1.5-8.7 mol d-1, respectively. A comparison between the main nutrient sources to the lagoon (karstic groundwater, recirculation, diffusion from sediments, inputs from a sewage treatment plant and atmospheric deposition) reveals that the recirculation of lagoon water through sediments is the main source of both dissolved inorganic nitrogen (DIN) and phosphorous (DIP) to La Palme lagoon. These results are in contrast with several studies conducted in systems influenced by terrestrial groundwater inputs, where groundwater is often assumed to be the main pathway for dissolved inorganic nutrient loads. This work highlights the important role of lagoon water recirculation through permeable sediments as a major conveyor of dissolved nutrients to coastal lagoons and, thus, the need for a sound understanding of the recirculation-driven nutrient fluxes and their ecological implications to sustainably manage lagoonal ecosystems.

Temperature dependency of virus and nanoparticle transport and retention in saturated porous media.

The influence of temperature on virus (PRD1 and ΦX174) and carboxyl-modified latex nanoparticle (50 and 100nm) attachment was examined in sand-packed columns under various physiochemical conditions. When the solution ionic strength (IS) equaled 10 and 30mM, the attachment rate coefficient (katt) increased up to 109% (p<0.0002) and the percentage of the sand surface area that contributed to attachment (Sf) increased up to 160% (p<0.002) when the temperature was increased from 4 to 20°C. Temperature effects at IS=10 and 30mM were also dependent on the system hydrodynamics; i.e., enhanced retention at a lower pore water velocity (0.1m/day). Conversely, this same temperature increase had a negligible influence on katt and Sf values when IS was 1mM or >50mM. An explanation for these observations was obtained from extended interaction energy calculations that considered nanoscale roughness and chemical heterogeneity on the sand surface. Interaction energy calculations demonstrated that the energy barrier to attachment in the primary minimum (∆Φa) decreased with increasing IS, chemical heterogeneity, and temperature, especially in the presence of small amounts of nanoscale roughness (e.g., roughness fraction of 0.05 and height of 20nm in the zone of influence). Temperature had a negligible effect on katt and Sf when the IS=1mM because of the large energy barrier, and at IS=50mM because of the absence of an energy barrier. Conversely, temperature had a large influence on katt and Sf when the IS was 10 and 30mM because of the presence of a small ∆Φa on sand with nanoscale roughness and a chemical (positive zeta potential) heterogeneity. This has large implications for setting parameters for the accurate modeling and transport prediction of virus and nanoparticle contaminants in ground water systems.

Transport and retention of bacteria and viruses in biochar-amended sand.

The transport and retention of Escherichia coli and bacteriophages (PRD1, MS2 and ФX174), as surrogates for human pathogenic bacteria and viruses, respectively, were studied in the sand that was amended with several types of biochar produced from various feedstocks. Batch and column studies were conducted to distinguish between the role of attachment and straining in microbe retention during transport. Batch experiments conducted at various solution chemistries showed negligible attachment of viruses and bacteria to biochar before or after chemical activation. At any given solution ionic strength, the attachment of viruses to sand was significantly higher than that of biochar, whereas bacteria showed no attachment to either sand or biochar. Consistent with batch results, biochar addition (10% w/w) to sand reduced virus retention in the column experiments, suggesting a potential negative impact of biochar application to soil on virus removal. In contrast, the retention of bacteria was enhanced in biochar-amended sand columns. However, elimination of the fine fraction (<60µm) of biochar particles in biochar-amended sand columns significantly reduced bacteria retention. Results from batch and column experiments suggest that land application of biochar may only play a role in microbe retention via straining, by alteration of pore size distribution, and not via attachment. Consequently, the particle size distribution of biochar and sediments is a more important factor than type of biochar in determining whether land application of biochar enhances or diminishes microbial retention.

Using Sealed Wells to Measure Water Levels Beneath Streams and Floodplains.

The design of wells beneath streams and floodplains has often employed with tall standpipes to prevent incursion of surface water into the well during flood events. Here, an approach has been presented to minimise the infrastructure demands in these environments by sealing the well top (e.g., prevent water entering the well) and monitor the total pressure in the water column using an absolute (non-vented) pressure transducer. The sealed well design was tested using a laboratory experiment where the total pressure responses were monitored in both an unsealed and sealed well, while the water level was varied. It is observed that, whether the well is sealed or not, the total pressure at a given depth in the aquifer will be equal to that within the well. This indicates that the sealed well design is a viable alternative to tall standpipes and also facilitates installation of wells beneath streams and floodplains.

Limitations of the use of environmental tracers to infer groundwater age.

Apparent ages obtained from the measured concentrations of environmental tracers have the potential to inform recharge rates, flow rates, and assist in the calibration of groundwater models. A number of studies have investigated sources of error in the relationships between the apparent ages, and the age assumed by models to relate this quantity to an aquifer property (e.g., recharge). These studies have also provided a number of techniques for correcting the known biases of apparent ages. In this paper, we review some of the concepts of age bias. We then demonstrate this bias through the use on four numerical examples, and assess the accuracy of correction methods in overcoming this bias. We examine this for CFCs, SF6, 3H/3He, 39Ar, and 14C. We demonstrate that in our four simple steady-state aquifer examples, bias occurs for a wide range of environmental tracers and flow configurations. When applying correction methods, we found that the values obtained are limited by the model assumptions. Models accounting for exchange with aquitards represent whole mobile zones and not discrete well screens. Mean transit times (comparable to mean ages) obtained from lumped parameter models deviate from actual values as the assumed distribution varies from the actual distribution. Methods that use multiple tracer ages are limited to ranges where both tracers report apparent ages. Our findings suggest that the incorporation of environmental tracer data into the understanding of groundwater systems requires approaches such as the direct use of concentrations, or the simulation of full age distributions.

Influence of Groundwater Hydraulic Gradient on Bank Storage Metrics.

The hydraulic gradient between aquifers and rivers is one of the most variable properties in a river/aquifer system. Detailed process understanding of bank storage under hydraulic gradients is obtained from a two-dimensional numerical model of a variably saturated aquifer slice perpendicular to a river. Exchange between the river and the aquifer occurs first at the interface with the unsaturated zone. The proportion of total water exchanged through the river bank compared to the river bed is a function of aquifer hydraulic conductivity, partial penetration, and hydraulic gradient. Total exchange may be estimated to within 50% using existing analytical solutions provided that unsaturated zone processes do not strongly influence exchange. Model-calculated bank storage is at a maximum when no hydraulic gradient is present and increases as the hydraulic conductivity increases. However, in the presence of a hydraulic gradient, the largest exchange flux or distance of penetration does not necessarily correspond to the highest hydraulic conductivity, as high hydraulic conductivity increases the components of exchange both into and out of an aquifer. Flood wave characteristics do not influence ambient groundwater discharge, and so in large floods, hydraulic gradients must be high to reduce the volume of bank storage. Practical measurement of bank storage metrics is problematic due to the limitations of available measurement technologies and the nested processes of exchange that occur at the river-aquifer interface. Proxies, such as time series concentration data in rivers and groundwater, require further development to be representative and quantitative.

Induced temperature gradients to examine groundwater flowpaths in open boreholes.

Techniques for characterizing the hydraulic properties and groundwater flow processes of aquifers are essential to design hydrogeologic conceptual models. In this study, rapid time series temperature profiles within open-groundwater wells in fractured rock were measured using fiber optic distributed temperature sensing (FO-DTS). To identify zones of active groundwater flow, two continuous electrical heating cables were installed alongside a FO-DTS cable to heat the column of water within the well and to create a temperature difference between the ambient temperature of the groundwater in the aquifer and that within the well. Additional tests were performed to examine the effects of pumping on hydraulic fracture interconnectivity around the well and to identify zones of increased groundwater flow. High- and low-resolution FO-DTS cable configurations were examined to test the sensitivities of the technique and compared with downhole video footage and geophysical logging to confirm the zones of active groundwater flow. Two examples are presented to demonstrate the usefulness of this new technique for rapid characterization of fracture zones in open boreholes. The combination of the FO-DTS and heating cable has excellent scope as a rapid appraisal tool for borehole construction design and improving hydrogeologic conceptual models.

Bias of apparent tracer ages in heterogeneous environments.

The interpretation of apparent ages often assumes that a water sample is composed of a single age. In heterogeneous aquifers, apparent ages estimated with environmental tracer methods do not reflect mean water ages because of the mixing of waters from many flow paths with different ages. This is due to nonlinear variations in atmospheric concentrations of the tracer with time resulting in biases of mixed concentrations used to determine apparent ages. The bias of these methods is rarely reported and has not been systematically evaluated in heterogeneous settings. We simulate residence time distributions (RTDs) and environmental tracers CFCs, SF6 , (85) Kr, and (39) Ar in synthetic heterogeneous confined aquifers and compare apparent ages to mean ages. Heterogeneity was simulated as both K-field variance (σ(2) ) and structure. We demonstrate that an increase in heterogeneity (increase in σ(2) or structure) results in an increase in the width of the RTD. In low heterogeneity cases, widths were generally on the order of 10 years and biases generally less than 10%. In high heterogeneity cases, widths can reach 100 s of years and biases can reach up to 100%. In cases where the temporal variations of atmospheric concentration of individual tracers vary, different patterns of bias are observed for the same mean age. We show that CFC-12 and CFC-113 ages may be used to correct for the mean age if analytical errors are small.

When can inverted water tables occur beneath streams?

Decline in regional water tables (RWT) can cause losing streams to disconnect from underlying aquifers. When this occurs, an inverted water table (IWT) will develop beneath the stream, and an unsaturated zone will be present between the IWT and the RWT. The IWT marks the base of the saturated zone beneath the stream. Although a few prior studies have suggested the likelihood of an IWT without a clogging layer, most of them have assumed that a low-permeability streambed is required to reduce infiltration from surface water to groundwater, and that the IWT only occurs at the bottom of the low-permeability layer. In this study, we use numerical simulations to show that the development of an IWT beneath an unclogged stream is theoretically possible under steady-state conditions. For a stream width of 1 m above a homogeneous and isotropic sand aquifer with a 47 m deep RWT (measured in an observation point 20 m away from the center of the stream), an IWT will occur provided that the stream depth is less than a critical value of 4.1 m. This critical stream depth is the maximum water depth in the stream to maintain the occurrence of an IWT. The critical stream depth decreases with stream width. For a stream width of 6 m, the critical stream depth is only 1 mm. Thus while theoretically possible, an IWT is unlikely to occur at steady state without a clogging layer, unless a stream is very narrow or shallow and the RWT is very deep.

An analysis of river bank slope and unsaturated flow effects on bank storage.

Recognizing the underlying mechanisms of bank storage and return flow is important for understanding streamflow hydrographs. Analytical models have been widely used to estimate the impacts of bank storage, but are often based on assumptions of conditions that are rarely found in the field, such as vertical river banks and saturated flow. Numerical simulations of bank storage and return flow in river-aquifer cross sections with vertical and sloping banks were undertaken using a fully-coupled, surface-subsurface flow model. Sloping river banks were found to increase the bank infiltration rates by 98% and storage volume by 40% for a bank slope of 3.4° from horizontal, and for a slope of 8.5°, delay bank return flow by more than four times compared with vertical river banks and saturated flow. The results suggested that conventional analytical approximations cannot adequately be used to quantify bank storage when bank slope is less than 60° from horizontal. Additionally, in the unconfined aquifers modeled, the analytical solutions did not accurately model bank storage and return flow even in rivers with vertical banks due to a violation of the dupuit assumption. Bank storage and return flow were also modeled for more realistic cross sections and river hydrograph from the Fitzroy River, Western Australia, to indicate the importance of accurately modeling sloping river banks at a field scale. Following a single wet season flood event of 12 m, results showed that it may take over 3.5 years for 50% of the bank storage volume to return to the river.

Disconnected surface water and groundwater: from theory to practice.

When describing the hydraulic relationship between rivers and aquifers, the term disconnected is frequently misunderstood or used in an incorrect way. The problem is compounded by the fact that there is no definitive literature on the topic of disconnected surface water and groundwater. We aim at closing this gap and begin the discussion with a short introduction to the historical background of the terminology. Even though a conceptual illustration of a disconnected system was published by Meinzer (1923), it is only within the last few years that the underlying physics of the disconnection process has been described. The importance of disconnected systems, however, is not widely appreciated. Although rarely explicitly stated, many approaches for predicting the impacts of groundwater development on surface water resources assume full connection. Furthermore, management policies often suggest that surface water and groundwater should only be managed jointly if they are connected. However, although lowering the water table beneath a disconnected section of a river will not change the infiltration rate at that point, it can increase the length of stream that is disconnected. Because knowing the state of connection is of fundamental importance for sustainable water management, robust field methods that allow the identification of the state of connection are required. Currently, disconnection is identified by showing that the infiltration rate from a stream to an underlying aquifer is independent of the water table position or by identifying an unsaturated zone under the stream. More field studies are required to develop better methods for the identification of disconnection and to quantify the implications of heterogeneity and clogging processes in the streambed on disconnection.