Comparative study of methods for WHPA delineation

Human activities, whether agricultural, industrial, commercial, or domestic, can contribute to ground water quality deterioration. In order to protect the ground water exploited by a production well, it is essential to develop a good knowledge of the flow system and to adequately delineate the area surrounding the well within which potential contamination sources should be managed. Many methods have been developed to delineate such a wellhead protection area (WHPA). The integration of more information on the geologic and hydrogeologic characteristics of the study area increases the precision of any given WHPA delineation method. From a practical point of view, the WHPA delineation methods allowing the simplest and least expensive integration of the available information should be favored. This paper presents a comparative study in which nine different WHPA delineation methods were applied to a well and a spring in an unconfined granular aquifer and to a well in a confined highly fractured rock aquifer. These methods range from simple approaches to complex computer models. Hydrogeological mapping and numerical modeling with MODFLOW-MODPATH were used as reference methods to respectively compare the delineation of the zone of contribution and the zone of travel obtained from the various WHPA methods. Although applied to simple ground water flow systems, these methods provided a relatively wide range of results. To allow a realistic delineation of the WHPA in aquifers of variable geometry, a WHPA delineation method should ensure a water balance and include observed or calculated regional flow characteristics.     

Hydromechanical reactivation of natural discontinuities: mesoscale experimental observations and DEM modeling

Fracture interaction mechanisms and reactivation of natural discontinuities under fluid pressurization conditions can represent critical issues in risk assessment of caprock integrity. A field injection test, carried out in a damage fault zone at the decameter scale, i.e., mesoscale, has been studied using a distinct element model. Given the complex structural nature of the damage fault zone hydraulically loaded, the contribution of fracture sets on the bulk permeability has been investigated. It has been shown that their orientation for a given in situ stress field plays a major role. Based on these results, a simpler model with a fluid-driven fracture intersecting a second fracture has been set up to perform a sensitivity analysis. It is in presence of a minimum differential stress value with a minimum angle with the maximum principal stress that the second fracture could be both, hydraulically and mechanically reactivated. Results also showed that in the vicinity of the fluid-driven fracture, a natural fracture will offer contrasted hydromechanical responses on each side of the intersection depending on the stress conditions and its orientation with respect to the stress field. In this case, we show that a hydromechanical decoupling can occur along the same plane. These results provide insights into fracture-controlled permeability of fault zones depending on the properties of the fractures and their hydromechanical interactions for a given in situ stress field.

     

A hierarchical generalized linear model with variable selection: studying the response of a representative fish assemblage for large European rivers in a multi-pressure context

Global change calls for an understanding of how temperature and flow regimes influence aquatic ecosystems. Fish assemblages are a major component of river ecosystems and are thought to exhibit more integrative informative responses than single species to environmental variations, whether rare and sudden or gradual and continuous. The use of long-term datasets is thus of primary importance, allied to statistical modeling. For each of three previously identified species clusters, we performed Bayesian variable selection and inference within a hierarchical log Poisson Generalized Linear Model using a spike and slab normal prior to pinpoint which subset of environmental variables is of importance for each fish assemblage. Fish counts from electrofishing experiments are known to provide overdispersed data and, not surprisingly, the contribution of recorded environmental effects is found to be weak compared with those of other intra-assemblage sources of variation. The posterior distribution of the regression parameters is in coherence with what was expected from biological knowledge of the three species clusters. In particular, thermophilic species tend to benefit from warmer waters, whereas the recruitment of cold water species decreases due to global warming effects. Our study provides an example of the advantages of hierarchical modeling for quantifying interspecies ecological effects and selecting common environmental variables of importance.     

Permeability Evolution During Non-linear Viscous Creep of Calcite Rocks

Permeability, storage capacity and volumetric strain were measured in situ during deformation of hot-pressed calcite aggregates containing 10, 20, and 30 wt% quartz. Both isostatic and conventional triaxial loading conditions were used. The tests were performed at confining pressure of 300 MPa, pore pressures between 50 to 290 MPa, temperatures from 673 to 873 K and strain rates of 3 × 10⁻⁵ s⁻¹. Argon gas was used as the pore fluid. The initial porosities of the starting samples varied from 5% to 9%, with higher porosity correlated to higher quartz content. Microstructural observations after the experiment indicate two kinds of pores are present: 1) Angular, crack-like pores along boundaries between quartz grains or between quartz and calcite grains and 2) equant and tubular voids within the calcite matrix. Under isostatic loading conditions, the compaction rate covaries with porosity and increases with increasing effective pressure. Most of the permeability reduction induced during compaction is irreversible and probably owes to plastic processes. As has been found in previous studies on hot-pressed calcite aggregates, permeability, k, is nonlinearly related to porosity, ϕ. Over small changes in porosity, the two parameters are approximately related as k ∝ ϕⁿ. The exponent n strongly increases as porosity decreases to a finite value (from about 4 to 6% depending on quartz content), suggesting a porosity percolation threshold. When subjected to triaxial deformation, the calcite-quartz aggregates exhibit shear-enhanced compaction, but permeability does not decrease as rapidly as it does under isostatic conditions. During triaxial compaction the exponent n only varies between 2 and 3. Non-isostatic deformation seems to reduce the percolation threshold, and, in fact, enhances the permeability relative to that at the same porosity during isostatic compaction. Our data provide constraints on the governing parameters of the compaction theory which describes fluid flow through a viscous matrix, and may have important implications for expulsion of sedimentary fluids, for fluid flow during deformation and metamorphism, and melt extraction from partially molten rocks.     

Constantes de formation des complexes hydroxydés de l'aluminium en solution aqueuse de 20 a 70°C

Stability constants of hydroxocomplexes of Al(III):Al(OH)²⁺ and A1(OH)₄^? have been measured in the 20–70°C temperature range by reactions involving only dissolved species. The stability constant ${}^1\text{K}{}_1$ of the first complex ion is studied by measuring pH of solutions of aluminium salts at several concentrations. ${}^1\text{β}{}_4$ of aluminate ion is deduced from equilibrium constants of the reaction between the trioxalato aluminium (III) complex ion and Al³⁺ in acid medium, and between the same complex ion and A1(OH)₄^? in alkaline medium. The K values and the associated ΔH are ${}^1\text{K}{}_1\text{= 10}{}^{\mathrm{?5.00}}$ and ΔH₁ = 11.8 Kcal; ${}^1\text{β}{}_4\text{= 10}{}^{\mathrm{?22.20}}$ and ΔH₄ = 42.45 Kcal. These last results are not in agreement with the values of recent tables for $\text{ΔG}{}^0\text{?}$ and $\text{ΔH}{}^0\text{?}$ of Al³⁺ and Al(OH)₄^?. We suggest a consistent set of data for dissolved and solid Al species and for some aluminosilicates.     

Large-scale GIS-based hydrogeological modeling of Flanders: a tool for groundwater management

For the implementation of the European Union Water Framework Directive (WFD), technological and scientific support are required. This paper presents a methodology to support a first step of the implementation of WFD, which is the delineation of groundwater bodies. The methodology consists of (1) the development of a complete and generally-accepted hydrogeological classification system for Flanders, named the HCOV code, (2) the development of a geographic information systems (GIS)-managed borehole database, and (3) the development of aquifer and aquitard models by means of a solid modeling approach. For each unit of the hydrogeological classification code for Flanders unit, GIS maps are generated for the three basic characteristics of hydrogeological layers: extent, base level and thickness, such that combined, the volume and extent of a hydrogeological layer is unambiguously defined. This GIS-based hydrogeological database has become a useful tool for groundwater management purposes and to provide the input for groundwater modeling.     

Analyse des rapports C/N/P du seston dans la partie orientale de l'atlantique equatorial

A study of the C/N ratios of particulate matter in the eastern part of the Gulf of Guinea demonstrated a statistically significant spatial variation in these ratios. The variation in the ratios ranged between 7.2 and 9.5 and were due in part to differences in the composition of the suspended organic matter and in part to its age and history. Young nutrient rich waters had low ratios (generally less than 7), whereas older waters often had ratios greater than 10. A relation between Chlorophyl a and the C/N ratio was also observed. The distribution of C/P ratios with depth showed a linear correlation with the vertical distribution of oxygen (r = 0.74). C/P values in the neighbourhood of 100 at the surface, increased with depth to around 300 at 1000 m.     

Vertical structure variability in a seasonal simulation of a medium-resolution regional model of the Eastern South Pacific

A seasonal simulation from a medium-resolution ocean general circulation mode (OGCM) is used to investigate the vertical structure variability of the Southeast Pacific (SEP). The focus is on the extra-tropical Rossby wave (ETRW) variability and associated forcing mechanism. Some aspects of the model mean state are validated from available observations, which justifies a vertical mode decomposition of the model variability. The analysis of the baroclinic mode contributions to sea level indicates that the gravest mode is dominant over most of the domain at all frequencies. Annual variability is on average twice as large as the semi-annual variability which is confined near the coast for all the modes. The first baroclinic mode contribution to the annual cycle exhibits a clear westward propagation north of the critical latitude. The higher-order modes only contribute near the coast where they are associated with vertically propagating energy. The residual variability, which is the energy at all timescales other than annual and semi-annual periods peaks offshore between 20°S and 30°S for all baroclinic modes. The third baroclinic mode also exhibits a relative maximum variability off the coast of Peru south of the critical latitude of the annual cycle (13°S), where the Peru–Chile Undercurrent is the most intense. Sensitivity experiments to the atmospheric and boundary forcing suggest that the residual variability results from the non-linear interaction between annual Rossby waves and the mean flow, while the annual ETRWs in the model result from the summed-contribution from both the local wind stress and remote equatorial forcing. Overall the study extends the classical analysis of sea level variability in the SEP based on linear theory, and suggests that the peculiarities of the baroclinic modes need to be taken into account for interpreting the sea level variability and understanding its connection with the equatorial variability.     

Comparison of two snowmelt modelling approaches in the Dudh Koshi basin (eastern Himalayas,Nepal)

Abstract

The glaciers in the Nepalese Himalayas are retreating due to rising temperatures. Lack of data and information on Nepal’s cryosphere has impeded scientific studies and field investigations in the Nepalese Himalayas. Therefore, IRD France and Ev-K2 CNR Italy have conducted the PAPRIKA (CryosPheric responses to Anthropogenic PRessures in the HIndu Kush-Himalaya regions: impacts on water resources and society adaptation in Nepal) project in Nepal with the financial support of the French and Italian scientific agencies. This project aims to address the current and future evolution of the cryosphere in response to overall environmental changes in South Asia, and its consequences for water resources in Nepal. Thus, two hydrological models, the GR4J lumped precipitation–runoff model and the snowmelt runoff model (SRM), were used in the Dudh Koshi basin. The GR4J model has been successfully applied in different parts of Europe. To obtain better results in such a harsh and rugged topography, modifications needed to be made, particularly in the snow module. The runoff pattern is analysed herein both for past years and, in a sensitivity analysis, for possible future climatic conditions (i.e. precipitation and temperature) using the SRM and GR4J modelling approaches. The results reveal a significant contribution of snow- and glacier-melt to runoff, and the SRM model shows better performance in Nepalese catchments than the GR4J model.

Editor D. Koutsoyiannis; Associate editor D. Gerten     

Chemical and K–Ar systematics of a mylonitic mica after natural and experimental interactions with varied fluids

Flow-through and batch-leaching experiments combining mineralogical, chemical and K–Ar isotopic analyses of mica separates from a mylonitic sample of the Tyndrum faulting system in Scotland (U.K.) were conducted to evaluate and simulate the natural interaction with H₂SO₄-loaded river water at pHs of 3–4. The flow-through experiments with H₂SO₄-loaded and pure deionized water completed by a batch-leaching experiment with 1 M HCl at room temperature had varied effects on soluble mineral phases, such as Fe-sulfates and Ca/Mg-carbonates, that were present in the varied mica size fractions, but none caused the mineralogical and K–Ar characteristics of the mica to differ from separates of the natural environment.Despite the limited number of K–Ar ages, size fractionation of the mylonitic mica identified two generations of mica with different REE patterns. The coarser (2–8 μm) fraction yielded a high La/Yb ratio of 4.3 and crystallized at 359 ± 6 Ma at a temperature of 250–300 °C, probably during a major tectonic-thermal activity. The smaller (< 0.5 μm) fraction yielded a La/Yb ratio of 2.1, a younger K–Ar age at 315 ± 5 Ma and a lower crystallization temperature of about 200 °C.