Unsaturated Zone Flow Processes and Aquifer Response Time in the Chalk Aquifer,Brighton, South East England

The Chalk aquifer is one of the main sources of water in South East England. The unsaturated zone in the aquifer plays an important role controlling the time and magnitude of recharge and is major pathway for contaminant transport to the water table. A range of previous work has addressed flow processes in the Chalk unsaturated zone, but physical understanding is still incomplete. Here we present the results of a study on flow mechanism in the Chalk unsaturated zone using a combination of statistical analysis and novel laboratory methods. The study was undertaken at three sites (North Heath Barn [NHB], Pyecombe East [PE], and Preston Park [PP]) on the Chalk of the Brighton block, South East England. Daily and hourly time series data of groundwater level and rainfall were correlated. The results show that a slower groundwater level response to rainfall occurs during dry seasons (summer and autumn) when the amount of effective rainfall is less than 4 mm/day, with a thicker and drier unsaturated zone. A faster response occurs during wet seasons (winter and spring) when the daily effective rainfall exceeds 4 mm/day with a thinner and wetter unsaturated zone. Periods of very rapid response (within 15 h) were observed during wet seasons at NHB and PE sites, with unsaturated hydraulic conductivity (K_u) inferred to reach 839 mm/day. A slower response was observed at an urbanized site (PP) as a result of reduction in direct recharge due to reduced infiltration, due to presences of impermeable infrastructure covering the area around PP borehole. Laboratory measurements of K_u of the Chalk matrix using a geotechnical centrifuge show variation from 4.27 to 0.07 mm/day, according to the level of saturation. Thus, the rapid responses cannot be linked to matrix flow only but indicate the contribution of fracture and karstic flow processes in conducting water.     

Modeling Gaseous CO2 Flow Behavior in Layered Basalts: Dimensional Analysis and Aquifer Response

Particular attention is paid to the risk of carbon dioxide (CO₂) leakage in geologic carbon sequestration (GCS) operations, as it might lead to the failure of sequestration efforts and to the contamination of underground sources of drinking water. As carbon dioxide would eventually reach shallower formations under its gaseous state, understanding its multiphase flow behavior is essential. To this aim, a hypothetical gaseous leak of carbon dioxide resulting from a well integrity failure of the GCS system in operation at Hellisheiði (CarbFix2) is here modeled. Simulations show that migration of gaseous carbon dioxide is largely affected by formation stratigraphy, intrinsic permeability, and retention properties, whereas the initial groundwater hydraulic gradient (0.0284) has practically no effect. In two different scenarios, about 18.3 and 30.6% of the CO₂ that would have been injected by the GCS system for 3 days could be potentially released again into the atmosphere due to a sustained leakage of the same duration. As the gaseous leak occurs, the aquifer experiences high pressure buildups, and the presence of a less conductive layer further magnifies these. Strikingly, the dimensional analysis showed that buoyant and viscous forces can be comparable over time within the predicted gaseous plumes, even far from the leakage source. Local pressure gradients, buoyant, viscous, and capillary forces all play an important role during leakage. Therefore, neglecting one or more of these contributions might lead to a partial prediction of gaseous CO₂ flow behavior in the subsurface, giving space to incorrect interpretations and wrong operational choices.     

Integrated Hydrological Modeling of Climate Change Impacts in a Snow-Influenced Catchment

The potential impact of climate change on water resources has been intensively studied for different regions and climates across the world. In regions where winter processes such as snowfall and melting play a significant role, anticipated changes in temperature might significantly affect hydrological systems. To address this impact, modifications have been made to the fully integrated surface-subsurface flow model HydroGeoSphere (HGS) to allow the simulation of snow accumulation and melting. The modified HGS model was used to assess the potential impact of climate change on surface and subsurface flow in the Saint-Charles River catchment, Quebec (Canada) for the period 2070 to 2100. The model was first developed and calibrated to reproduce observed streamflow and hydraulic heads for current climate conditions. The calibrated model was then used with three different climate scenarios to simulate surface flow and groundwater dynamics for the 2070 to 2100 period. Winter stream discharges are predicted to increase by about 80, 120, and 150% for the three scenarios due to warmer winters, leading to more liquid precipitation and more snowmelt. Conversely, the summer stream discharges are predicted to fall by about 10, 15, and 20% due to an increase in evapotranspiration. However, the annual mean stream discharge should remain stable (±0.1 m³/s). The predicted increase in hydraulic heads in winter may reach 15 m and the maximum decrease in summer may reach 3 m. Simulations show that winter processes play a key role in the seasonal modifications anticipated for surface and subsurface flow dynamics.     

Analytical Solution for Modeling Discharge into a Tunnel Drilled in a Heterogeneous Unconfined Aquifer

Predicting transient inflow rates into a tunnel is an important issue faced by hydrogeologists. Most existing analytical solutions overestimate the initial discharge due to the assumption that drilling was instantaneous over the entire tunnel length. In addition, they assume a homogeneous system. An alternative model was recently developed for tunnels intersecting heterogeneous formations, but its application was reduced to the case of confined flow to deep tunnels in weakly diffusive aquifers. In this paper, we adapt existing analytical solutions for drainage systems to the specific case of a tunnel progressively drilled in a highly diffusive heterogeneous unconfined aquifer. The case of a tunnel overlying an impervious layer is analytically solved by applying the superposition principle, while the case of a tunnel constructed some distance above an impervious layer is solved by discretizing the tunnel length into subsectors. Both models can simulate transient discharge into a tunnel drilled at various speeds through a heterogeneous unconfined aquifer, and allow the prediction of discharge rates in shallow tunnels located in highly diffusive aquifers. We successfully applied this approach to a tunnel in heterogeneous volcanic rock.     

Predicting Pesticide Attenuation in a Fractured Aquifer Using Lumped‐Parameter Models

Neighboring springs draining fractured‐rock aquifers can display large differences in water quality and flow regime, depending on local variations of the connectivity and the aperture size distribution of the fracture network. Consequently, because homogeneous equivalent parameters cannot be assumed a priori for the entire regional aquifer, the vulnerability to pollution of such springs has to be studied on a case by case basis. In this paper, a simple lumped‐parameter model usually applied to estimate the mean transit time of water (or tracer) is presented. The original exponential piston‐flow model was modified to take land‐use distribution into account and applied to predict the evolution of atrazine concentration in a series of springs draining a fractured sandstone aquifer in Luxembourg, where despite a nationwide ban in 2005, atrazine concentrations still had not begun to decrease in 2009. This persistence could be explained by exponentially distributed residence times in the aquifer, demonstrating that in some real world cases, models based on the groundwater residence time distribution can be a powerful tool for trend reversal assessments as recommended for instance by current European Union guidelines.     

Aquifer Testing, Mathematical Modeling, and Regulatory Risk

Mathematical modeling founded on a strong field data base can be a valuable tool for the analysis of ground-water contamination problems. The purposes of this paper are threefold: (1) we demonstrate the dilemma of a knowledgeable ground-water quality regulator whose regulatory decision-making process is confronted with the output of a mathematical model that is based on very limited field test data; (2) we demonstrate a method available to a knowledgeable regulator for assessing approximately a range of possible performances of a contaminated ground-water recovery well field using a range of input data derived from a very limited data base; and (3) we present a strong case for presenting mathematical model outputs as ranges of values rather than as unique solutions. A range is determined by an examination of the level of sophistication of the field data base. Our experience with 12 field sites wherein ground-water contamination has occurred has led us to conclude that field data are seldom, if ever, adequate to defend a unique solution from a mathematical model. Regulatory decisions generally can be reduced to minimization of risks based on the smallest range of model outputs that can be defended on the basis of the field data base. The more limited the field data base, the greater must be the range of defensible model outputs, and consequently, the greater the risk inherent in subsequent regulatory decisions. The knowledgeable regulator can assess the risks in the regulatory decision-making process only if he is able to assess the extent to which the field data base for the mathematical model output reflects state-of-the-art data collection and analysis technologies and methodologies. If an applicant for a permit or license submits a less than adequate data base and concomitantly a broad range of defensible model outputs, he inherently requests that the knowledgeable regulator accept a risk greater than that required if adequate aquifer testing techniques had been employed.