Fonctionnement et gestion des aquifères alluviaux de haute altitude

The Haute Sarine alpine valley is located at an average altitude of approximately 1'000 to 1'200 masl. It was formed in its southern part up to the area of Gstaad to the favour of a typical North-South tectonic fault according to the main structures of the Préalpes, then followed normal fault striking East-West to end up on the rock bolt of the Vanel hill. The late and postglacial period was marked by the deposition of glaciolacustrines, glaciofluviatile and flurecent fluviatile sediments which have filled the ice free valley bottom probably starting from a period older than the Bölling-Alleröd interstadial. The depth of this overdeepening varies between 20 and more than 50 m of depth. Drillings which we carried out could document series of typical glaciolacustrine sediments. The implementation of several geophysical survey methods (resistivity, seismic refraction and reflexion, well logs and pumping tests) enabled us to analyse several sources of information which we confronted between them in order to produce a 3D model of the quaternary deposits and their hydrological properties. Most of these sediments are relatively quite permeable and build a free groundwater reservoir that we subdivided into 3 parts. The basin of Gsteig in the south is closed, the valley having been blocked by a massive landslide that probably occurred in several successive stages starting at the first climate improvements. The second part, the basin of Feutersoey has not been retained as a potential important hydrogeologic resource based on our preliminary investigations. The last basin, that of Gstaad-Saanen is the largest and thickest. It is also in this basin that most of the tourist and industrial activities of the area are concentrated. The observation network was established in order to characterize the groundwater flows enabling us at the same time to carry out piezometric measurements and chemical as well as and bacteriological analyses. We also carried out a δ18O isotope analysis campaign. The whole of these approaches highlighted the particular characteristics of the groundwater recharge. The groundwater system is mainly fed by the Sarine and the precipitations infiltration. We also located an underground karstic feeding of evaporitic type. We finally propose a groundwater management and protection concept adapted to the particular situation of the Sarine alluvial mountain aquifer system on both quantitative and qualitative aspects. The high altitude alpine groundwater catchment areas such as the Haute Sarine represent as a whole a quality resource, with high potentials yet to be defined on a larger scale, for example by applying numerical modeling and by organizing and maintaining long term observation networks that will enable us to define a sustainable groundwater resource management within the alpine area.

Typology of recent groundwaters from different aquifer environments based on geogenic tracer elements

A key problem in environmental science and in engineering geology is the often incomplete understanding of the origin of dissolved components in groundwater. The dissolved contents of trace elements in groundwater are of special importance for groundwater quality control. The AQUITYP project aims to establish a detailed typology of recent groundwaters based on their geogenic trace element compositions, and to derive a so-called "geo-reference" for groundwaters from five principal aquifer lithologies in the Alpine belt. This geo-reference provides a database for investigations related to groundwater contamination, groundwater resources management and engineering geology. Groundwaters from crystalline, carbonate, and evaporite rocks, as well as molasse and flysch sediments in Switzerland and in neighbouring countries were sampled and documented by previous researchers (Dubois, 1993; Dematteis, 1995; Mandia, 1993; Hesske, 1995; Basabe, 1993). Based on a statistical analysis of the data and examination of the relationship between aquifer lithology and chemical groundwater composition these researchers identified a number of characteristic tracer elements (geogenic tracers). The present study includes (1) a synthesis of the hydrogeology and the hydrochemistry of recent groundwaters in these five aquifer types based on groundwater data acquired within the AQUITYP project, and (2) a hydro-geochemical investigation of the origin and chemical behaviour of the geogenic tracer chromium, based on a comprehensive field and laboratory study. To enable a comparison of the chemical groundwater data gathered and analysed over a time span of 20 years, a rigorous quality control of the entire database was made. An assessment was made of the sampling techniques, the sample treatment and the analytical methods applied in the AQUITYP project since 1981. Different tests were carried out to evaluate the quality and comparability of the data, including geochemical model calculations, comparison of different analytical techniques, and tests to estimate the influence of the filtration procedure. In order to make this large number of quality controlled data accessible also for future investigations, a groundwater data storage system was developed (AQUITYP-DataBase). 1.) Typology of recent groundwaters In the synthesis of groundwater hydrochemistry, the emphasis was laid on the hydro-geochemical evolution leading to the characteristic groundwater composition in each of the five aquifer types. Chemical characteristics and differences between the groundwaters from the different aquifer types were identified and natural concentration ranges for each aquifer type derived. The proposed geogenic tracers were evaluated and the potential sources of these tracers identified. The dominant processes leading to the typical mineralisation of recent groundwaters were investigated using geochemical modelling strategies. Finally, the concentrations of chemical elements in the groundwaters from the different aquifer types were compared to the Swiss drinking water standards in order to assess the quality of the investigated groundwaters. It has been found that each rock type contributes in a characteristic way to the major and trace element composition of the corresponding groundwater: The groundwaters derived from the crystalline Mont-Blanc and Aiguilles-Rouges massifs are characterised by a low total mineralisation (TDS 22 to 158 mg/L) dominated by Ca2+, Na+, Mg2+, alkalinity, SO42-, and F- (Ca-Na-HCO3,-SO4, waters). Elevated amounts of Mo, U, W, and As occur. These groundwaters derive their mineralisation mainly from the interaction with hydrothermal minerals present along fractures. Fractures act as major groundwater flow paths. Minerals relevant for groundwater mineralisation include carbonates (Ca2+, Mg2+, HCO3-), clay minerals (Ca-Na ion exchange), fluorite (F-, Ca2+), Fe-, As-, and Mo-sulphides (SO42-, As, Mo), and U- and W-minerals (U, W). In these crystalline groundwaters the natural concentrations of F- (23% of the investigated springs) and As (7%) exceed the Swiss limits for drinking water. In addition, the WHO guideline values for U are exceeded in 65% of the cases and for Mo in 15% of the cases. The carbonate karst groundwaters obtain their low to intermediate mineralisation (TDS 161 to 547 mg/L) from the dissolution of calcite (Ca-HCO3, waters), as well as in certain regions dolomite (Ca-Mg-HCO3, waters) and gypsum (Ca-Mg-HCO3,-SO4, waters). Together with their very short residence times, the carbonate karst groundwaters generally contain very low trace element contents. Nevertheless, geogenic trace elements occur in specific regions in relation with fossil organic matter (I, V) and accessory minerals such as barite in deep sea limestones (Ba), evaporite minerals (gypsum, celestite: Sr, Li), clay- and Fe-minerals (V), and Mo-sulphides and U-minerals in dolomitic limestones (Mo, U). In 18% of the carbonate karst springs atmospheric derived Pb exceeds the Swiss drinking water quality target value. The groundwaters from Triassic evaporites in the Swiss Rhone basin are characterised by a high total mineralisation (TDS 760 to 2788 mg/L) expressed by elevated amounts of Ca2+, Mg2+, Sr2+, SO42-, and alkalinity (Ca-Mg-SO4-HCO3 waters). Elevated amounts of the trace elements Mn, Ni, Cu, Li, Rb, Y, and Cd occur. The hydrochemical evolution of these groundwaters is governed by incipient dedolomitisation, involving dissolution of gypsum, celestite and dolomite and simultaneous precipitation of calcite. The Na+ and K+ contents probably controlled by ion-exchange reactions on clay minerals. Characteristic trace elements originate mainly from the dissolution of dolomite (Mn, Ni) and small amounts of apatite (Y, Cd), and from the oxidation of sulphide minerals (Cu, Ni, Cd). The elevated concentrations of the highly soluble Li and Rb may be related to brine inclusions in evaporite minerals and eventually to clay minerals. These evaporite groundwaters contain SO42- concentrations exceeding the Swiss quality target for drinking water in all springs, and the concentrations of U and Ni exceed the WHO guideline values in respectively 58% and 2% of the cases. In addition, Mn, Cd, and As concentrations exceed the Swiss quality targets in respectively 11%, 10%, and 7% of the investigated springs. Recent groundwaters circulating in the porous and fissured molasse sandstones and conglomerates acquire their intermediate mineralisation (TDS 48 to 714 mg/L) primarily by dissolution of calcite and minor dolomite (Ca-Mg-HCO3 waters). The particular mineralogy of certain Molasse formations is reflected in specific trace element compositions of corresponding groundwaters: ophiolite detritus in OMM sandstones in western Switzerland (Cr), barite fracture mineralisations in subalpine and folded Molasse units (Ba), granitic detritus containing sulphides (Mo), U-minerals (U) and abundant mica (Li) in the "Glimmersand" (OSM, Ca-Mg-HCO3-SO4 waters), and evaporite minerals (Li, Sr), sulphides (Mo), and U-minerals (U) in the "Gypsum-bearing Molasse" (USM, highly mineralised Ca-Mg-SO4-HCO3 waters, TDS 984 to 1346 mg/L) . In these molasse groundwaters the Swiss quality target values for drinking water of Cr and Pb are exceeded in 36% and 6% of the springs, respectively. The U concentrations of 14% of the molasse groundwaters ("Gypsum-bearing Molasse" and "Glimmersand"), exceed the WHO guideline value. The groundwaters from the "Gypsum-bearing Molasse" display similar quality problems as the evaporite groundwaters. Groundwaters derived from the shallow flysch aquifers in the Niesen and Gurnigel nappes are poorly evolved Ca-(Mg)-HCO3 waters. Their low to intermediate total mineralisation (TDS 160 to 459 mg/L) is acquired primarily by dissolution of calcite and to a lesser degree dolomite. The low trace element content is dominated by Ba originating from barite fracture mineralisations. The poor chemical evolution of recent flysch groundwaters results from (1) their short residence time in the fractured flysch rocks, and (2) the absence of readily dissolving minerals except carbonates and barite. In 32% of the flysch springs Pb derived from atmospheric sources exceeds the Swiss quality target value. 2.) Case study on the chemical weathering of molasse sandstone: Sources and chemical behaviour of Cr For the characterisation of the potentially toxic tracer element Cr a comprehensive field study was carried out on a selected catchment (Lutry spring catchment near Lausanne) situated in a molasse sandstone (OMM, Burdigalian). The investigation of the processes controlling the dissolved Cr content in these groundwaters was based on the groundwater chemical data, as well as on mineralogical, geochemical, and hydrological data. Relic Cr-bearing spinel and pyroxene in the sandstone were identified to be the primary sources of Cr. An electron microscope study showed that in the Burdigalian sandstone and overlying soil these minerals are strongly weathered. The slow weathering of these minerals is the major Cr releasing process. Under the oxidising conditions reigning in the investigated groundwaters, Cr prevails in solution in its highly soluble and toxic hexavalent state (CrO42-). In this case, retention by secondary Cr-hydroxide phases does not occur, as can be shown by geochemical model calculations. Laboratory leaching experiments were carried out with Burdigalian molasse sandstone from the field site, in order to support the field study findings and to quantify the processes responsible for the Cr mineralisation observed in the Lutry groundwater. Two experiments with mountain-wet oxidised and reduced sandstone were carried out over a time span of 2 months each, to obtain information about the influence of the oxidation state of the substratum on the mobilisation of Cr. The experiments clearly showed that the Cr-releasing processes are fast enough to explain the Cr contents found in the groundwater, and that the release of Cr into the groundwater depends on the weathering state of the sandstone. In the oxidised, Fe-hydroxide-coated sandstone, Cr is faster released into solution than in the less altered reduced sandstone. This indicates that the Cr contained in the reduced sandstone is in a more stable state, i.e. mainly incorporated in detrital minerals, while in the oxidised sandstone, the Cr is partly in an unstable state, i.e. adsorbed on the surface of secondary Fe-hydroxides, from where it is more easily leached.

EPFL2001

The effect of seasonal soil frost on the alpine groundwater recharge including climate change aspects

Daniel Bayard

In alpine areas, the snow cover plays an important role as a water reservoir. Water is stored as snow over the winter and released in spring, recharging mountain aquifers through infiltration. These aquifers are essential, especially for supplying water for human activities during dry seasons. Numerous studies have shown that locally soil frost can drastically reduce the water infiltration. However, we know much less about the hydrological impact of soil frost at a larger scale, in particular with regard to groundwater recharge. This was our motivation for initiating an extensive field experiment in the southern Swiss Alps. Several methods at different spatial scales were adapted to explore (a) the local effect of a partially frozen ground on the snowmelt discharge in an alpine area, (b) the key processes influencing groundwater recharge during snowmelt periods, such as snow cover and soil frost evolution, snowmelt water runoff types, (c) the large-scale effect of soil frost on the aquifer recharge, and (d) winter situations that are critical with respect to flooding. In order to take into account spatial, altitudinal and climatic differences, the field study was run for two winter seasons at Gd St Bernard (2500 m), a site with very high winter precipitation and strong winds, and at Hannigalp above Grächen (2100 m), a ski resort with a rather dry winter climate. The different components of the soil water balance (lateral runoff, deep percolation, liquid soil water content) were measured on delimited plots of 5 m2. Additionally, a dye tracer experiment visualized the snowmelt infiltration patterns into the frozen, respectively non-frozen ground. At Gd St Bernard, investigations of the liquid water content and the soil temperature were additionally conducted at two secondary sites, differing in their orientation, so as to gain a better insight of the spatial soil frost expansion. Finally, an analysis of the water-table variation at the village of Grächen (some 500 m below Hannigalp) during the snowmelt was carried out for a 10-year time period, including our experimental seasons. We complemented our field investigations by using a numerical soil-snow-atmosphere transfer model that simulates the seasonal development of the soil frost and the snow cover, as well as lateral and vertical water flow in the top soil. This model was used as an upper boundary condition to a 2-D groundwater model, to calculate the outlet flow at the base of a simplified alpine aquifer. The two winters investigated had contrasting meteorological conditions, which resulted in intriguing differences and similarities with respect to the soil physical conditions: the first one having a thick snowpack and hardly any soil frost, and the second one with hardly any snow until January and a deep and persistent soil frost. In spite of a very thick snowpack, the water balance measurements from winter 2000/2001 showed for both sites zero to low surface and subsurface flow due to the fact that soil frost was local and the soil permeability high. During the next winter, approximately 25% of the total meltwater run off laterally. The soil infiltration capacity was mainly reduced by the presence of a basal ice sheet, caused by the freezing of meltwater after a mid-winter snowmelt event. The field data were used to calibrate and validate a physically-based one-dimensional model. To allow an accurate parameter setting for the water repellent soil of Hannigalp, laboratory infiltration experiments on soil columns were additionally carried out. Using the resultant information, the model predicted correctly most measured soil parameters, in particular good matching was observed between measured and simulated snowmelt discharge. At Gd St Bernard, the simulated results were less satisfactory, especially during the first winter, when large discrepancies were noted between measured and simulated lateral runoff. Apart from the shortage in the soil parameter characterization, the deviation was a result of the simplified description of the experimental field. The model was, for example, not able to simulate surface runoff generated by a steep slope under unfrozen conditions, as the soil surface is considered by the model to be flat. Nevertheless, results were accurate enough to reproduce the occurrences of snowmelt events and the accumulated values of snow-melt discharge pathways under frozen conditions. Long-term weather data were used to simulate the soil frost depth at both sites. The simulation results showed that the probability of the occurrence of frozen soils was lowest on slopes with south or north exposure. On southerly exposed plots, the heat stored in the soil during the summer inhibited the development of deep and persistent soil frost, while most pore ice thawed due to underneath heating. On north plots, the building of soil frost was rare, as early snowfall insulated the ground from the atmosphere. However, during snow-poor winters, the frost penetrated deep into the ground, and the low underneath heat flux at that location hardly affected the frost depth until snowmelt. The snow depth and average soil temperature were hence the two components affecting the frost depth extent. The altitude soil frost variation was related to the spatial extent of these two components. The higher the altitude, the higher the precipitation and the lower the mean air/surface temperature. Simulating the whole Grächen recharge area, seasonal soil frost was rare between 1800 m and 2000 m, as the snowpack was thick enough to insulate the ground, and the soil was warm enough to inhibit deep soil frost. In lower areas, the shallow snowpack enabled the soil to freeze during most winters, whereas at higher areas, deep and persistent soil frost was simulated due to the low underneath heating. A 10-year water-table depth record was used to examine possible large-scale effects of seasonal soil frost. We noted that the water-table rise at snowmelt was lowest at the end of winters with extensive soil frost. However, such winters were characterized by low precipitation, and reduced snow cover. In contrast, winters with little soil frost were mostly snow rich. Consequently, only a reduced number of winters were directly comparable. Considering these winters, the reduction in the water-table rise at snowmelt due to seasonal soil frost varied only between 10 and 30%, as, during snowmelt, most meltwater was able to re-infiltrate into the very permeable ground in lower unfrozen areas. We may assume that a soil with a smaller infiltration capacity can increase the influence of the soil frost on the aquifer recharge. Finally, we investigated the effect of a changing climate on the hydraulic and thermic regimes at different altitudes in the area of Grächen. We assume an increase in the air temperature of 2°C generates deeper soil frost on lower areas, as the resultant smaller snow cover does not insulate the soil from the atmosphere any more. However, it does reduce the soil frost on higher, snow rich areas, due to the warmer air temperature. As a whole, minor changes were simulated on the mean groundwater recharge, as more/less meltwater was able to infiltrate in higher/lower areas. By increasing additionally the precipitation by 15%, the soil frost diminished slightly over the whole catchment. By combining both climate change scenarios, we may expect an increase in extreme events, as, on the one hand, the rainfall intensity increases, and, on the other hand, rain on snow events over frozen soils occurs more often in lower areas. In practice, this work shows the importance of better investigating the hydrothermal behaviour of alpine regions. Such a knowledge is necessary, as more and more efforts are put into managing efficiently, durably and globally the available water resources. This work stresses also the necessity to consider seasonally soil frost in alpine risk management. As an example, flooding is mostly a consequence of strong rain precipitation at high altitudes combined with heavy snowmelt. Obviously, the risk of extreme runoff events increases during frozen winters, when the soil infiltration capacity is further reduced by the soil ice. Also, it has been shown that unstable slopes are mostly related to specific hydrogeological situations, like strong precipitation, which induces a sharp increase in the underneath water circulation. The melt water in spring is hence a potential activating factor, and a frozen soil may reduce the risk of sliding, as less water infiltrates into the ground.

EPFL2003

Méthodes d'analyse du fonctionnement hydrogéologique des versants instables

Pierre Tullen

Significant progress was made these last decades in the development of hydrogeological methods. However, these new methods are still very little used in the problematics of unstable slopes, mainly because of the hydrogeological complexity, which characterizes them (thin aquifers, discontinuous media, succession of saturated and unsaturated zones, lower permeability, higher hydraulic gradients, evolution of the permeability in time in relation to slope movements, etc.). Consequently, the approach usually developed to study slope instabilities is based on the setting-up of relations between hydroclimatologic parameters and slope movements. This approach exclusively concentrates on causes and effects, without developing the processes induced by these causes and responsible for these effects. This "black-box" type of relation is in some cases satisfactory, but in others, the analysis of the induced hydrogeological processes becomes essential to understand the functioning of the instability phenomenon and to predict the movements. Our project aims at introducing this missing link to allow a better understanding of the hydrogeological processes caused by particular hydroclimatologic conditions. It is a matter of determining the relation between hydroclimatologic phenomena and a critical hydrogeological state, which can lead a phenomenon of instability to reactivation or to a state of rupture. Unfortunately, this relation cannot be approached in the same way for all instabilities throughout the world. Owing to the specificity of each slope, the relation must be established individually, which may lead to a hydrogeological typology of unstable slopes. Accordingly, we propose a methodology comprising a succession of stages that are useful for the understanding of the hydrogeological processes underlying slopes. This approach aims at establishing a hydrogeological assessment on the scale of the slope or, if necessary, on a regional scale, in order to determine groundwater infiltration zones, trajectory and velocity of groundwater run-off as well as outlet areas. The usual ways of doing this are generally the measurement of the river flow rate, a pluviometric follow-up or direct hydrogeological measurements via piezometers. These traditional approaches generally remain limited by problems of space representativeness or life span; on top of this come problems of cost due to the often chaotic morphology of the slope instabilities and the fact that they are not easily accessible. To make up for these various problems, we propose a complete methodological approach to improve the understanding of the groundwater run-off systems for the unstable slopes: A characterization of the geological domain based on geophysical measurements in addition to the borehole data. This approach, original in the specific context of unstable slopes, makes it possible to obtain a three-dimensional representation of the heterogeneity which composes the subsurface grounds in a fast, efficient and economic way; A spatialization of the climatologic observations in order to take into account the distribution of the precipitations while considering the various states (rain and/or snow) that constitute the water contribution to a slope in relation to altitude; A hydrogeochemical characterization of the outlets, which makes it possible to assess the general groundwater circulation system within the slope, and to come to a decision about one or several possible trajectories of groundwater run-off; An isotopic characterization of the outlets, which allows estimating the altitude of infiltration of the groundwater. This characterization can also give some information on the infiltration zones. These methods also make it possible to characterize the general groundwater run-off system. In some cases, it is also possible to come to a conclusion, quantitatively speaking, concerning the mixing between several reservoirs and to estimate the time of residence of the water in the aquiferous system; A qualitative validation of the assumptions on the basis of artificial tracing tests which will attest the existence of hydrogeological connections between a catchment area and an outlet. This approach will also inform about the time of residence of water; A validation of the conceptual model on the basis of a bi- or three-dimensional numerical simulation of the groundwater run-off in order to validate the hydrological and hydrogeological relations. The calibration of these numerical models is based on a set of hydrogeological data, in particular the outlet discharge, the isotopic composition of water and/or their hydrochemical composition. These strictly hydrogeological numerical simulations can then be translated into a mechanical state to evaluate the stability of the slope. The whole of the points suggested by this methodological approach is not always essential to determine the hydrogeological processes that govern each slope. According to the characteristics of each study site and the means at disposal, some could be treated only partly, or not at all, depending on the established objectives. In the end, a hydrogeological typology of the unstable slopes should make it possible to propose a standard methodological approach and specific tools for analysis in direct relationship to the hydrogeological context which characterizes each slope. The methodological approach suggested in this PhD thesis was tested and validated on three representative sites through Switzerland and the alpine arc that present quite distinct characteristics of hydrogeological functioning. They are the slopes of Hohberg in the canton of Fribourg, Triesenberg in the Principality of Liechtenstein and Peney in the canton of Geneva. This validation phase allowed to assess the hydrogeological functioning of the three studied slopes and to evaluate the relation between the hydrogeological processes and the phenomena of instability. An approach such as that presented in this study opens up prospects to carry out means of sustainable stabilization with an aim of solving problems of slope instabilities.

EPFL2002