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.

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.

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).

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.

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.

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.

Applied tracer tests in fractured rock: Can we predict natural gradient solute transport more accurately than fracture and matrix parameters?

Applied tracer tests provide a means to estimate aquifer parameters in fractured rock. The traditional approach to analysing these tests has been using a single fracture model to find the parameter values that generate the best fit to the measured breakthrough curve. In many cases, the ultimate aim is to predict solute transport under the natural gradient. Usually, no confidence limits are placed on parameter values and the impact of parameter errors on predictions of solute transport is not discussed. The assumption inherent in this approach is that the parameters determined under forced conditions will enable prediction of solute transport under the natural gradient. This paper considers the parameter and prediction uncertainty that might arise from analysis of breakthrough curves obtained from forced gradient applied tracer tests. By adding noise to an exact solution for transport in a single fracture in a porous matrix we create multiple realisations of an initial breakthrough curve. A least squares fitting routine is used to obtain a fit to each realisation, yielding a range of parameter values rather than a single set of absolute values. The suite of parameters is then used to make predictions of solute transport under lower hydraulic gradients and the uncertainty of estimated parameters and subsequent predictions of solute transport is compared. The results of this study show that predictions of breakthrough curve characteristics (first inflection point time, peak arrival time and peak concentration) for groundwater flow speeds with orders of magnitude smaller than that at which a test is conducted can sometimes be determined even more accurately than the fracture and matrix parameters.

Tracer mass recovery in fractured aquifers estimated from multiple well tests.

Forced-gradient tracer tests in fractured aquifers often report low mass recoveries. In fractured aquifers, fractures intersected by one borehole may not be intersected by another. As a result (1) injected tracer can follow pathways away from the withdrawal well causing low mass recovery and (2) recovered water can follow pathways not connected to the injection well causing significant tracer dilution. These two effects occur along with other forms of apparent mass loss. If the strength of the connection between wells and the amount of dilution can be predicted ahead of time, tracer tests can be designed to optimize mass recovery and dilution. A technique is developed to use hydraulic tests in fractured aquifers to calculate the conductance (strength of connection) between well pairs and to predict mass recovery and amount of dilution during forced gradient tracer tests. Flow is considered to take place through conduits, which connect the wells to each other and to distant sources or sinks. Mass recovery is related to the proportion of flow leaving the injection well and arriving at the withdrawal well, and dilution is related to the proportion of the flow from the withdrawal well that is derived from the injection well. The technique can be used to choose well pairs for tracer tests, what injection and withdrawal rates to use, and which direction to establish the hydraulic gradient to maximize mass recovery and/or minimize dilution. The method is applied to several tracer tests in fractured aquifers in the Clare Valley, South Australia.

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.

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.

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.

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.

Analysis of a vertical dipole tracer test in highly fractured rock.

The results of a vertical dipole tracer experiment performed in highly fractured rocks of the Clare Valley, South Australia, are presented. The injection and withdrawal piezometers were both screened over 3 m and were separated by 6 m (midpoint to midpoint). Due to the long screen length, several fracture sets were intersected, some of which do not connect the two piezometers. Dissolved helium and bromide were injected into the dipole flow field for 75 minutes, followed by an additional 510 minutes of flushing. The breakthrough of helium was retarded relative to bromide, as was expected due to the greater aqueous diffusion coefficient of helium. Also, only -25% of the total mass injected of both tracers was recovered. Modeling of the tracer transport was accomplished using an analytical one-dimensional flow and transport model for flow through a fracture with diffusion into the matrix. The assumptions made include: streamlines connecting the injection and withdrawal point can be modeled as a dipole of equal strength, flow along each streamline is one dimensional, and there is a constant Peclet number for each streamline. In contrast to many other field tracer studies performed in fractured rock, the actual travel length between piezometers was not known. Modeling was accomplished by fitting the characteristics of the tracer breakthrough curves (BTCs), such as arrival times of the peak concentration and the center of mass. The important steps were to determine the fracture aperture (240 microm) based on the parameters that influence the rate of matrix diffusion (this controls the arrival time of the peak concentration); estimating the travel distance (11 m) by fitting the time of arrival of the centers of mass of the tracers; and estimating fracture dispersivity (0.5 m) by fitting the times that the inflection points occurred on the front and back limbs of the BTCs. This method works even though there was dilution in the withdrawal well, the amount of which can be estimated by determining the value that the modeled concentrations need to be reduced to fit the data (approximately 50%). The use of two tracers with different diffusion coefficients was not necessary, but it provides important checks in the modeling process because the apparent retardation between the two tracers is evidence of matrix diffusion and the BTCs of both tracers need to be accurately modeled by the best fit parameters.

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.

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.