Heat Tracing in a Fractured Aquifer with Injection of Hot and Cold Water.

Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log-log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures.

Site-selection criteria for the Einstein Telescope.

The Einstein Telescope (ET) is a proposed next-generation, underground gravitational-wave detector to be based in Europe. It will provide about an order of magnitude sensitivity increase with respect to the currently operating detectors and, also extend the observation band targeting frequencies as low as 3 Hz. One of the first decisions that needs to be made is about the future ET site following an in-depth site characterization. Site evaluation and selection is a complicated process, which takes into account science, financial, political, and socio-economic criteria. In this paper, we provide an overview of the site-selection criteria for ET, provide a formalism to evaluate the direct impact of environmental noise on ET sensitivity, and outline the necessary elements of a site-characterization campaign.

Parametric assessment of hydrochemical changes associated to underground pumped hydropower storage.

Underground pumped hydropower storage (UPHS) using abandoned mines is an alternative to store and produce electricity in flat regions. Excess of electricity is stored in form of potential energy by pumping mine water to a surface reservoir. When the demand of electricity increases, water is discharged into the mine (i.e., underground reservoir) through turbines producing electricity. During the complete operational process of UPHS plants, hydrochemical characteristics of water evolve continuously to be in equilibrium successively with the atmosphere (in the surface reservoir) and the surrounding porous medium (in the underground reservoir). It may lead to precipitation and/or dissolution of minerals and their associated consequences, such as pH variations. Induced hydrochemical changes may have an impact on the environment and/or the efficiency (e.g., corrosions and incrustations affect facilities) of UPHS plants. The nature of the hydrochemical changes is controlled by the specific chemical characteristics of the surrounding porous medium. However, the magnitude of the changes also depends on other variables, such as hydraulic parameters. The role of these parameters is established to define screening criteria and improve the selection procedure of abandoned mines for constructing UPHS plants. This work evaluates the role of the main hydrogeological factors for three different chemical composition of the porous medium. Results are obtained by means of numerical reactive transport modeling. Potential impacts on the environment (mainly on groundwater and surface water bodies) and on the efficiency of the UPHS plants vary considerably from a hydraulic parameter to another showing the need for a detailed characterization before choosing locations of future UPHS plants.

Model-based classification of CPT data and automated lithostratigraphic mapping for high-resolution characterization of a heterogeneous sedimentary aquifer.

Cone penetration testing (CPT) is one of the most efficient and versatile methods currently available for geotechnical, lithostratigraphic and hydrogeological site characterization. Currently available methods for soil behaviour type classification (SBT) of CPT data however have severe limitations, often restricting their application to a local scale. For parameterization of regional groundwater flow or geotechnical models, and delineation of regional hydro- or lithostratigraphy, regional SBT classification would be very useful. This paper investigates the use of model-based clustering for SBT classification, and the influence of different clustering approaches on the properties and spatial distribution of the obtained soil classes. We additionally propose a methodology for automated lithostratigraphic mapping of regionally occurring sedimentary units using SBT classification. The methodology is applied to a large CPT dataset, covering a groundwater basin of ~60 km2 with predominantly unconsolidated sandy sediments in northern Belgium. Results show that the model-based approach is superior in detecting the true lithological classes when compared to more frequently applied unsupervised classification approaches or literature classification diagrams. We demonstrate that automated mapping of lithostratigraphic units using advanced SBT classification techniques can provide a large gain in efficiency, compared to more time-consuming manual approaches and yields at least equally accurate results.

Factors controlling spatial and temporal patterns of multiple pesticide compounds in groundwater (Hesbaye chalk aquifer, Belgium).

Factors governing spatial and temporal patterns of pesticide compounds (pesticides and metabolites) concentrations in chalk aquifers remain unclear due to complex flow processes and multiple sources. To uncover which factors govern pesticide compound concentrations in a chalk aquifer, we develop a methodology based on time series analyses, uni- and multivariate statistics accounting for concentrations below detection limits. The methodology is applied to long records (1996-2013) of a restricted compound (bentazone), three banned compounds (atrazine, diuron and simazine) and two metabolites (deethylatrazine (DEA) and 2,6-dichlorobenzamide (BAM)) sampled in the Hesbaye chalk aquifer in Belgium. In the confined area, all compounds had non-detects fractions >80%. By contrast, maximum concentrations exceeded EU's drinking-water standard (100 ng L-1) in the unconfined area. This contrast confirms that recent recharge and polluted water did not reach the confined area, yet. Multivariate analyses based on variables representative of the hydrogeological setting revealed higher diuron and simazine concentrations in the southeast of the unconfined area, where urban activities dominate land use and where the aquifer lacks protection from a less permeable layer of hardened chalk. At individual sites, positive correlations (up to τ=0.48 for bentazone) between pesticide compound concentrations and multi-annual groundwater level fluctuations confirm occurrences of remobilization. A downward temporal trend of atrazine concentrations likely reflects decreasing use of this compound over the last 28 years. However, the lack of a break in concentrations time series and maximum concentrations of atrazine, simazine, DEA and BAM exceeding EU's standard post-ban years provide evidence of persistence. Contrasting upward trends in bentazone concentrations show that a time lag is required for restriction measures to be efficient. These results shed light on factors governing pesticide compound concentrations in chalk aquifers. The developed methodology is not restricted to chalk aquifers, it could be transposed to study other pollutants with concentrations below detection limits.

Physically based groundwater vulnerability assessment using sensitivity analysis methods.

A general physically based method is presented to assess the vulnerability of groundwater to external pressures by numerical simulation of groundwater flow. The concept of groundwater vulnerability assessment considered here is based on the calculation of sensitivity coefficients for a user-defined groundwater state for which we propose several physically based indicators. Two sensitivity analysis methods are presented: the sensitivity equation method and the adjoint operator method. We show how careful selection of a method can significantly minimize the computational effort. An illustration of the general methodology is presented for the Herten aquifer analog (Germany). This application to a simple, yet insightful, case demonstrates the potential use of this general and physically based vulnerability assessment method to complex aquifers.

Regional transport modelling for nitrate trend assessment and forecasting in a chalk aquifer.

Regional degradation of groundwater resources by nitrate has become one of the main challenges for water managers worldwide. Regulations have been defined to reverse observed nitrate trends in groundwater bodies, such as the Water Framework Directive and the Groundwater Daughter Directive in the European Union. In such a context, one of the main challenges remains to develop efficient approaches for groundwater quality assessment at regional scale, including quantitative numerical modelling, as a decision support for groundwater management. A new approach combining the use of environmental tracers and the innovative 'Hybrid Finite Element Mixing Cell' (HFEMC) modelling technique is developed to study and forecast the groundwater quality at the regional scale, with an application to a regional chalk aquifer in the Geer basin in Belgium. Tritium data and nitrate time series are used to produce a conceptual model for regional groundwater flow and contaminant transport in the combined unsaturated and saturated zones of the chalk aquifer. This shows that the spatial distribution of the contamination in the Geer basin is essentially linked to the hydrodynamic conditions prevailing in the basin, more precisely to groundwater age and mixing and not to the spatial patterns of land use or local hydrodispersive processes. A three-dimensional regional scale groundwater flow and solute transport model is developed. It is able to reproduce the spatial patterns of tritium and nitrate and the observed nitrate trends in the chalk aquifer and it is used to predict the evolution of nitrate concentrations in the basin. The modelling application shows that the global inertia of groundwater quality is strong in the basin and trend reversal is not expected to occur before the 2015 deadline fixed by the European Water Framework Directive. The expected time required for trend reversal ranges between 5 and more than 50 years, depending on the location in the basin and the expected reduction in nitrate application. To reach a good chemical status, nitrate concentrations in the infiltrating water should be reduced as soon as possible below 50mg/l; however, even in that case, more than 50 years is needed to fully reverse upward trends.

A new tracer technique for monitoring groundwater fluxes: the Finite Volume Point Dilution Method.

Quantification of pollutant mass fluxes is essential for assessing the impact of contaminated sites on their surrounding environment, particularly on adjacent surface water bodies. In this context, it is essential to quantify but also to be able to monitor the variations with time of Darcy fluxes in relation with changes in hydrogeological conditions and groundwater - surface water interactions. A new tracer technique is proposed that generalizes the single-well point dilution method to the case of finite volumes of tracer fluid and water flush. It is called the Finite Volume Point Dilution Method (FVPDM). It is based on an analytical solution derived from a mathematical model proposed recently to accurately model tracer injection into a well. Using a non-dimensional formulation of the analytical solution, a sensitivity analysis is performed on the concentration evolution in the injection well, according to tracer injection conditions and well-aquifer interactions. Based on this analysis, optimised field techniques and interpretation methods are proposed. The new tracer technique is easier to implement in the field than the classical point dilution method while it further allows monitoring temporal changes of the magnitude of estimated Darcy fluxes, which is not the case for the former technique. The new technique was applied to two experimental sites with contrasting objectives, geological and hydrogeological conditions, and field equipment facilities. In both cases, field tracer concentrations monitored in the injection wells were used to fit the calculated modelled concentrations by adjusting the apparent Darcy flux crossing the well screens. Modelling results are very satisfactory and indicate that the methodology is efficient and accurate, with a wide range of potential applications in different environments and experimental conditions, including the monitoring with time of changes in Darcy fluxes.

Influence of injection conditions on field tracer experiments.

Calibration of ground water transport models is often performed using results of field tracer experiments. However, little attention is usually paid to the influence, on resulting breakthrough curves, of injection conditions and well-aquifer interactions, more particularly of the influence of the possible trapping of the tracer in the injection wellbore. Recently, a new mathematical and numerical approach has been developed to model injection conditions and well-aquifer interactions in a very accurate way. Using an analytical solution derived from this model, a detailed analysis is made of the evolution of the tracer input function in the aquifer. By varying injection conditions from one simulation to another, synthetic breakthrough curves are generated with the SUFT3D ground water flow and transport finite-element simulator. These tests show clearly that the shape of the breakthrough curves can be dramatically affected by injection conditions. Using generated breakthrough curves as "actual" field results, a calibration of hydrodispersive parameters is performed, neglecting the influence of injection conditions. This shows that neglecting the influence of actual injection conditions can lead to (1) errors on fitted parameters and (2) misleading identification of the active transport processes. Conclusions and guidelines are drawn in terms of proposed methodologies for better controlling the tracer injection in the field, in order to minimize risk of misinterpretation of results.

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation.

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.