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

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.

Geological heterogeneity in landslides

Boris Matti

Significant progress has been made these last decades in the development of hydrogeological numerical flow modelling for describing the hydrodynamic behaviour of landslides. However, these new sophisticated methods are still very seldom used in the problems of slope instability in particular because of the hydrogeological complexity which characterizes them; thin aquifers, discontinuous media, succession of saturated and unsaturated zones, low permeabilities, high hydraulic gradients, lithological heterogeneity, strong contrasts of permeabilities and heterogeneous infiltration. Predictive models of flow in the subsurface, which are often based on homogeneous porous media types of representation, are badly adapted to natural systems that are characterized by highly heterogeneous media such as landslides. These models are good and reliable on a landslide scale (regional scale), but their quality may be affected on a local scale by strong geological heterogeneities. Geological heterogeneities of the subsurface take part in determining the hydrodynamical and geomechanical behaviour of landslides. However, their spatial distribution is partially unknown. Thus, the principal objectives of this PhD thesis are: (i) To carry out an integrated multidisciplinary characterization study on the internal structure of landslides in flysch and Quaternary environments, in order to clarify the organisation of the geological heterogeneities and to identify the hydrodynamic implications. (ii) To propose a conceptual model representing the geological architecture and the hydrogeological functioning. (iii) To examine the effects of heterogeneity and anisotropy on flow systems. (iv) To better understand the influence of geological heterogeneities on the mechanical behaviour of large landslides by performing numerical sensitivity analyses, by means of different heterogeneity scenarios on the field parameters. (v) Finally, to test the incidences on slope stabilization techniques; evaluation of the efficiency of a drainage gallery work. The main test site of la Frasse landslide (VD, Switzerland) was chosen, and completed with additional landslide cases. The main results are the following: In most of the case studies, the landslide mass is composed of an old prehistoric stabilized mass, pinched between the active sliding mass and the bedrock, and playing an important hydrologic role. The stabilized mass and the bedrock form the substratum of the landslide. Landslides occurring in these types of media are defined by an organized heterogeneous environment with "fracture" flows and discontinuity porosity. The overall hydraulic conductivity is low, and locally high permeable zones exist. Regional groundwater circulations are limited and form local interconnected aquicludes organised in thin aquifers, and presenting saturated and unsaturated zones. The hydrogeological analyses showed that the system presents a bimodal permeability; (i) Low hydraulic conductivities characterizing the global matrix and defining the capacitive fraction, and (ii) high permeable features, with high hydraulic conductivities defining the conductive fraction, and favouring strong channelling effects. Besides, the observation shows that the aquifer system is generally very reactive with important magnitudes. Often, there is a straight correlation between water level variation and climatic conditions (rainy events). Landslides are characterized by two important inflows namely effective infiltration from the surface and lateral inflows from the neighbouring units. Water transfer between the stabilized mass and the active mass may be important and thus have to be considered. The existence of water transfer between the bedrock and the landslide mass (stabilized and active) is not well established. The bedrock and the landslide mass present a hydrological behavioural independence. Theoretical two- and three-dimensional flow models are used to investigate the effects of the spatial variability of the hydraulic conductivity on the underground flows. The role of the connectivity in generating flow channelling is examined thanks to the observation of close relations between the permeability and the hydraulic pressures. The sensitivity analysis shows clearly that the relation between local permeability and hydraulic pressures is not straight, and that the organization of the flows depends on the heterogeneity of the hydraulic properties and their spatial correlation. Strong channelling effects are observed in highly heterogeneous porous media. The development of flow channelling as a function of the variance of the natural log permeability values and the correlation lengths is demonstrated. The integrated multi-disciplinary geological characterization at the La Frasse test site combined with the hydrogeological and lithological data of several additional case studies led to the proposal of a global conceptual model. The following assumptions are considered to enable a subsequent quantification of flow components: The flow occurs under confined to leaky conditions, with leakage varying in space; The flow framework is controlled by a complex multi-layer system, isolated lenses or perched aquifer; The aquifer system is divided into interconnected hydrological zones presenting various degrees of saturation; Each hydrological zone may function individually from the others; Horizontally and vertically, the flow direction in the porous matrix is affected by prevailing structural patterns generating channeling effects; The flow is multidirectional, free and channelized, and is affected by temporal and spatial changes; The aquifer is under an unsteady flow regime due to seasonal variation of natural gradients; A conceptual model based on a simple reservoir approach is proposed. It allows the representation of most of the field observations and the main characteristics, namely the organized heterogeneity and the duality of the aquifers. The system is represented by various reservoirs more or less connected and saturated. Complex storage capacities and plug-flow effects may record past events and reactive sliding processes several months after the last important rainy event. The analysis shows that function of the capacity and the degree of saturation of the system, an important hydrological event is not necessarily associated to a reactivation. And, according to the degree of complexity of the system (saturation, connectivity...) a localized geological modification (variation of permeability, reservoir burst...) may produce a chain reaction, and generate failures in unexpected places. The conductive fraction favours the drainage of the system, whereas the capacitive fraction controls the distribution of the hydraulic heads. The role of the phreatic nappe, through the conductive fraction, is to drain and control the hydrologic equilibrium of the system. Therefore landslide remediation with the help of a deep drainage gallery is obviously the most valuable method for this type of landslide. It supports and enhances the natural effects of the conductive fraction in draining the system. Finally, in this context the efficiency of civil engineering works was evaluated according to the heterogeneity of the medium. This study describes transient hydrogeological and geomechanical models realized jointly in 2006 by the EPFL and GeoMod SA within the framework of the stabilization work of the La Frasse landslide. These models evaluate the impact of a deep drainage gallery with subvertical pipes towards the surface in terms of reduction of the deformation velocities and increase of the factor of safety of the landslide. Three variants consisting of different inter pipe spacings are tested. Considering the local heterogeneities, the results show that a mean spacing between the pipes of the order of 10 m is able to control the temporal head fluctuations between the wells within a range of some meters. Moreover, this solution induces a strong diminution of the predicted displacements during a specific crisis, from 101cm for the model without drainage to around 14 cm for the drained model, and a significant gain of security (from 1.05 to 1.30).

EPFL2008

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

Typologie géochimique des eaux des aquifères carbonatés des chaînes alpines d'Europe centrale et méridionale

This study describes the spatial heterogeneousness of the chemical composition of underground waters, particularly the trace elements, and the carbon-13 isotope in several carbonate aquifers of the alpine belt in Europe. The research forms part of the AQUITYP project which have been worked on the last 14 years by the Geology laboratory of the Swiss Federal Institute of Technology in Lausanne (GEOLEP). The main aim of the project is to define the typology of the aquifers, based on the chemical characterisation of underground waters. The carbonate aquifers has been divided in several groups chosen after geological, geographical and hydrological criterias. The geological criterias are petrography and depositional environments of carbonates. The geographical criteria is the altitude of the catchment area. The criteria hydrology is the type of flow in karstic systems. The first stage of this research consisted of a selection of 87 carbonate aquifers divided in 13 regions of study in the alpine belt between France and Greece. Aquifers were selected in order to guarantee the absence of external contamination in groundwaters. Sampling were carried out in a period of base-flow, when groundwaters are equilibrated with rocks. The analyse of 112 samples of water provided the basis of the observations of the various chemical and isotopic composition of groundwaters of carbonate aquifers. Lixiviation of rocks forming aquifers provides information about geogenic origin of some trace elements. The geological and hydrogeological context of aquifers The aquifers studied are composed of platform and basin limestones, that form part of the paleogeographic domain of the Tethys. One example of Devonian limestone belongs to the complex of the Palaeozoic of Graz is also included. Each aquifer has homogenous petrographic composition, but from one aquifer to another the petrographic composition may differ. In this research are included pure limestones, marly limestones and dolomites. From a hydrogeological point of view the 87 aquifers show different degrees of karstification, they are placed either in regions of mountains (caracterised by an important hydraulic gradient), either in region of plateau. The surface of the catchment areas is variable between ten km2 (for example the karstic plateau in Slovenia) and a few km2 (for example the Parmelan massif from subalpine belt in France). Discharges of springs are also variable, between 1 and 1000 l/s. The chemical characterisation of water The study of the chemical typology of the waters of carbonate aquifers are carried out on non thermal water rich in oxygen (median values: temperature 8.8°C, dissolved oxygen 11.8 mg/l, Eh 0.4 V, pH 7.2). The total dissolved solids differs between 100 and 600 mg/l. Underground waters are calcium bicarbonate type. Within these facies in some cases, the Mg and S04 content can be superior, thus defining two sub-facies. The iodine content in the groundwaters of carbonate aquifers usually shows higher values than the other types of aquifers of the project AQUITYP. Concentration in the groundwaters of carbonate aquifers is variable between 10 and 60 µg/l. The origin of iodine is caused by the oxidation of fossil organic matter present in the carbonate rocks. Trace elements allows to distinguish groups of aquifers composed of pure limestones, dolomite or deep sea limestones : pure limestones. The aquifers composed of pure limestones are marked by a very small concentration of trace elements. This can be explained by the poorness of the detrital materials (clays, Fe-Mn oxides, heavy minerals) and the compact lithology that reduces the surface of water-rock interaction. dolomite. The typical elements of dolomite aquifers are Mg, U and Mo. The high primary porosity and the brittle deformation of this type of rock favorise a better water rock interaction compared to poor reef limestone aquifers. Examples of this group of aquifers have been selected in the Italian Dolomites and in the Zone of Karavanke in Slovenia deep sea limestones. The barium is a geogenic tracer in some aquifers composed of basin limestones. The aquifers of the Maiolica and of the Scaglia in the Apennine and the Malm aquifer (spring of Grandchamp) in the Prealpes in Switzerland have a high level of Ba in groundwater (between 60 and 220 µg/l). The limiting factor of Ba mobility in this carbonate environment is the concentration of SO42- in the water. The SO42- makes the Ba precipitate in shape of barytes. This process explains the presence of barytes in fractures found in the aquifer of the Calcare Massiccio in the Apennines. The content of carbonate-13 in groundwaters allows to make a distinction between waters infiltrated in catchment areas with or without vegetation and soils. This isotope is also interpreted as a tracer of the type of underground flow in Karstic system (rapid flow in Karstic channels or slow flow in little fractures of the rock). For the project an example has been taken in Tanneben karst, in Austria. This karst shows that water with a slow flow from blocks has a higher 13C content in comparison with fast flow channel water. Finally some analysis of waters from the bottom of soils on the Jura belt in Switzerland showed that during little rain periods the content of isotope is comparable to the one of slow flow waters of the unsaturated zone in the karst. This demonstrate that the isotopic exchange with rock is completed at low horizon in the soil. The vanadium is a geogenic tracer in the Malm aquifer of the Jura belt in Switzerland and in the Cretaceous aquifer in the region of Trieste, in Italy. The origin of vanadium differs in these two examples. In the Jura belt the vanadium comes probably from marly layers of limestones of the low Portlandien. In the region of Trieste the element vanadium is present in bituminous carbonate sediments from inner platform (Calcari bituminosi from Paleocene). In conclusion the importance and the limits of the chemical criterion for the typology of aquifers are shown. This method of distinction of groups of limestone aquifers is one of the tools which permit the determination of the origin of water sources, particularly in underground and borehole construction projects.

EPFL1995