Search

International audience; Surface and subsurface flows interact at different spatial and temporal scales through the development of saturated areas occurring when the aquifer reaches the surface. While this interaction exerts a strong control in the partitioning of water between base flow, seepage flow and overland flow, its quantification remains a challenge. Here, we propose a novel modeling approach based on two equivalent hillslopes to capture spatial and temporal variabilities of the main processes. We calibrate their subsurface hydraulic properties based on the temporal dynamics of stream discharge. The model is tested on two pilot catchments located in Brittany (France). For both catchments, the model is successfully calibrated on 40 years of stream discharge data. The results demonstrate that contrasted hydraulic properties are required, with: (1) a relatively low conductive hillslope (conductivities between 3x10–8 m/s and 2x10–6 m/s) enhancing overland flows during recharge periods while also sustaining low flows in the late recession period through slow aquifer discharge, and (2) a highly diffusive hillslope (diffusivity between 5x10–3 m2/s and 3x10–1 m2/s) dominantly shaping the event- to seasonal-scale streamflow recession behavior. The strong contrast of the two hillslopes reveals the fundamental role of heterogeneity in controlling groundwater storage-discharge functions, ruling out any homogeneous equivalent at the catchment scale. Low flows appear to be a non-obvious combination of the contributions of the two hillslopes, the less diffusive hillslope catching up to the more conductive one by the end of the low flow period. Catchment-scale responses integrate complex interactions between recharge, groundwater and overland flows through the volume of subsurface storage and the extent of the saturated area..

International audience; A long duration pumping test conducted over 151 days in a fractured sandstone and shale formation displays a nonstandard drawdown response and anomalous pressure diffusion, which cannot be properly interpreted using existing frameworks (e.g., homogeneous, double porosity, boundary conditions, and fractal models). An alternative framework with simple geometry and more complex hydraulic properties is thus proposed to interpret such kind of drawdown responses. The analytical development allows first to demonstrate all scaling relations in this interpretation framework. Then, and most importantly, the multi-scale hydraulic test provides consistent scalings of transmissivity, T, to storativity, S, over distances ranging from 83 to 383 m in a faulted area. Drawdown analysis in several monitoring wells shows persistent decrease of transmissivity in highly channelized fracture flow structures. In one structure, the cubic dependency of transmissivity to storativity identifies a well-defined fault and also demonstrates the validity of Poiseuille flow at a scale rarely investigated. In the other structure, the linear dependency of transmissivity to storativity indicates that the flow-bearing structure is the surrounding fracture network. Well-designed pumping tests combined with scaling analysis driven by geological evidence thus provide essential information on flow-bearing structures for site characterization and modeling tasks. At least for moderate to low permeable fractured rocks, the scaling of transmissivity to storativity appears to be more informative than any separate interpretation of hydraulic property scaling exponents.

International audience; Hydrological models and the data derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have been widely used to study the variations in terrestrial water storage (TWS) over large regions. However, both GRACE products and model results suffer from inherent uncertainties, calling for the need to make a combined use of GRACE and models to examine the variations in total TWS and their individual components, especially in relation to natural and human-induced changes in the terrestrial water cycle. In this study, we use the results from two state-of-the-art hydrological models and different GRACE spherical harmonic products to examine the variations in TWS and its individual components, and to attribute the changes to natural and human-induced factors over large global river basins. Analysis of the spatial patterns of the long-term trend in TWS from the two models and GRACE suggests that both models capture the GRACE-measured direction of change, but differ from GRACE as well as each other in terms of the magnitude over different regions. A detailed analysis of the seasonal cycle of TWS variations over 30 river basins shows notable differences not only between models and GRACE but also among different GRACE products and between the two models. Further, it is found that while one model performs well in highly-managed river basins, it fails to reproduce the GRACE-observed signal in snow-dominated regions, and vice versa. The isolation of natural and human-induced changes in TWS in some of the managed basins reveals a consistently declining TWS trend during 2002–2010, however; significant differences are again obvious both between GRACE and models and among different GRACE products and models. Results from the decomposition of the TWS signal into the general trend and seasonality indicate that both models do not adequately capture both the trend and seasonality in the managed or snow-dominated basins implying that the TWS variations from a single model cannot be reliably used for all global regions. It is also found that the uncertainties arising from climate forcing datasets can introduce significant additional uncertainties, making direct comparison of model results and GRACE products even more difficult. Our results highlight the need to further improve the representation of human land-water management and snow processes in large-scale models to enable a reliable use of models and GRACE to study the changes in freshwater systems in all global regions.

International audience; Interest in the residence time distribution (RTD) as a comprehensive measure of subsurface hydrologic systems is growing. This focus is resulting from recognition that diverse vadose zone, groundwater flows, and transfer between hydrological compartments, are fundamentally related to the system RTD

International audience; Nitrogen pollution of freshwater and estuarine environments is one of the most urgent environmental crises. Shallow aquifers with predominantly local flow circulation are particularly vulnerable to agricultural contaminants. Water transit time and flow path are key controls on catchment nitrogen retention and removal capacity, but the relative importance of hydrogeological and topographical factors in determining these parameters is still uncertain. We used groundwater dating and numerical modeling techniques to assess transit time and flow path in an unconfined aquifer in Brittany, France. The 35.5 km2 study catchment has a crystalline basement underneath a ∼60 m thick weathered and fractured layer, and is separated into a distinct upland and lowland area by an 80 m-high butte. We used groundwater discharge and groundwater ages derived from chlorofluorocarbon (CFC) concentration to calibrate a free-surface flow model simulating groundwater flow circulation. We found that groundwater flow was highly local (mean travel distance = 350 m), substantially smaller than the typical distance between neighboring streams (∼1 km), while CFC-based ages were quite old (mean = 40 years). Sensitivity analysis revealed that groundwater travel distances were not sensitive to geological parameters (i.e. arrangement of geological layers and permeability profile) within the constraints of the CFC age data. However, circulation was sensitive to topography in the lowland area where the water table was near the land surface, and to recharge rate in the upland area where water input modulated the free surface of the aquifer. We quantified these differences with a local groundwater ratio (rGW-LOCAL), defined as the mean groundwater travel distance divided by the mean of the reference surface distances (the distance water would have to travel across the surface of the digital elevation model). Lowland, rGW-LOCAL was near 1, indicating primarily topographical controls. Upland, rGW-LOCAL was 1.6, meaning the groundwater recharge area is almost twice as large as the topographically-defined catchment for any given point. The ratio rGW-LOCAL is sensitive to recharge conditions as well as topography and it could be used to compare controls on groundwater circulation within or between catchments.

International audience; This review presents the physical mechanisms generating residence time distributions (RTDs) in hydrologic systems with a focus on steady-state analytical solutions. Steady-state approximations of the RTD in hydrologic systems have seen widespread use over the last half-century because they provide a convenient, simplified modeling framework for a wide range of problems. The concept of an RTD is useful anytime that characterization of the timescales of flow and transport in hydrologic systems is important, which includes topics like water quality, water resource management, contaminant transport, and ecosystem preservation. Analytical solutions are often adopted as a model of the RTD and a broad spectrum of models from many disciplines has been applied. Although these solutions are typically reduced in dimensionality and limited in complexity, their ease of use makes them preferred tools, specifically for the interpretation of tracer data. Our review begins with the mechanistic basis for the governing equations, highlighting the physics for generating a RTD, and a catalog of analytical solutions follows. This catalog explains the geometry, boundary conditions and physical aspects of the hydrologic systems, as well as the sampling conditions, that altogether give rise to specific RTDs. The similarities between models are noted, as are the appropriate conditions for their applicability. The presentation of simple solutions is followed by a presentation of more complicated analytical models for RTDs, including serial and parallel combinations, lagged systems, and non-Fickian models. The conditions for the appropriate use of analytical solutions are discussed, and we close with some thoughts on potential applications, alternative approaches, and future directions for modeling hydrologic residence time.

International audience; The detection of preferential flow paths and the characterization of their hydraulic properties are major challenges in fractured rock hydrology. In this study, we propose to use temperature as a passive tracer to characterize fracture connectivity and hydraulic properties. In particular, we propose a new temperature tomography field method in which borehole temperature profiles are measured under different pumping conditions by changing successively the pumping and observation boreholes. To interpret these temperature- depth profiles, we propose a three step inversion-based framework. We consider first an inverse model that allows for automatic permeable fracture detection from borehole temperature profiles under pumping conditions. Then we apply a borehole-scale flow and temperature model to produce flowmeter profiles by inversion of temperature profiles. This second step uses inversion to characterize the relationship between temperature variations with depth and borehole flow velocities (Klepikova et al., 2011). The third inverse step, which exploits cross-borehole flowmeter tests, is aimed at inferring inter-borehole fracture connectivity and transmissivities. This multi-step inverse framework provides a means of including temperature profiles to image fracture hydraulic properties and connectivity. We test the proposed approach with field data obtained from the Ploemeur (N.W. France) fractured rock aquifer, where the full temperature tomography experiment was carried out between three 100 m depth boreholes 10 m apart. We identified several transmissive fractures and their connectivity which correspond to known fractures and corroborate well with independent information, including available borehole flowmeter tests and geophysical data. Hence, although indirect, temperature tomography appears to be a promising approach for characterizing connectivity patterns and transmissivities of the main flow paths in fractured rock.

International audience; Groundwater in shallow weathered and fractured crystalline rock aquifers is often the only perennial water resource, especially in semi-arid region such as Southern India. Understanding groundwater flows in such a context is of prime importance for sustainable aquifer management. Here, we describe a detailed study of fracture properties and relate the hydraulic connectivity of fractures to groundwater flows at local and watershed scales. Investigations were carried out at a dedicated Experimental Hydrogeological Park in Andhra Pradesh (Southern India) where a large network of observation boreholes has been set up. Twenty-height boreholes have been drilled in a small area of about 18,000 m2 in which borehole loggings and hydraulic tests were carried out to locate the main flowing fractured zones and investigate fractures connectivity. Several hydraulic tests (nineteen slug tests and three pumping tests) performed under two water level conditions revealed contrasting behavior. Under high water level conditions, the interface including the bottom of the saprolite and the first flowing fractured zone in the upper part of the granite controls groundwater flows at the watershed-scale. Under low water level conditions, the aquifer is characterized by lateral compartmentalization due to a decrease in the number of flowing fractures with depth. Depending on the water level conditions, the aquifer shifts from a watershed flow system to independent local flow systems. A conceptual groundwater flow model, which includes depth-dependent fracture connectivity, is proposed to illustrate this contrasting hydrological behavior. Implications for watershed hydrology, groundwater chemistry and aquifer vulnerability are also discussed.

International audience; Peatlands are complex ecosystems driven by many physical, chemical, and biological processes. Peat soils have a significant impact on water quality, ecosystem productivity and greenhouse gas emissions. However, the extent of peatlands is decreasing across the world, mainly because of anthropogenic activities such as drainage for agriculture or groundwater abstractions in underlying aquifers. Potential changes in precipitation and temperature in the future are likely to apply additional pressure to wetland. In this context, a methodology for assessing and comparing the respective impacts of groundwater abstraction and climate change on a groundwater-fed wetland (135 km2) located in Northwest France, is presented. A groundwater model was developed, using flexible boundary conditions to represent surface-subsurface interactions which allowed examination of the extent of the wetland areas. This variable parameter is highly important for land management and is usually not considered in impact studies. The model was coupled with recharge estimation, groundwater abstraction scenarios, and climate change scenarios downscaled from 14 GCMs corresponding to the A1B greenhouse gas (GHG) scenario over the periods 1961-2000 and 2081-2100. Results show that climate change is expected to have an important impact and reduce the surface of wetlands by 5.3-13.6%. In comparison, the impact of groundwater abstraction (100% increase in the expected scenarios) would lead to a maximum decrease of 3.7%. Results also show that the impacts of climate change and groundwater abstraction could be partially mitigated by decreasing or stopping land drainage in specific parts of the area. Water management will require an appropriate compromise which encompasses ecosystem preservation, economic and public domain activities.

International audience; Several recharge processes may be defined for karst systems. However, their respective influence is difficult to determine from spring-water chemistry without a good knowledge of the input signal, i.e. rainfall chemistry. Given the little data published concerning rainfall chemistry in the French Mediterranean area, two rainfall gauges networks were installed – one along a transect from the Mediterranean coast to the most elevated part of the Hérault region, which was monitored over two hydrological cycles (1996–1998), and a second on the catchment area of the main karst aquifer of the Gard region (Fontaine de Nîmes karst system), which was monitored over four hydrological cycles (1998–2002). The main objective of this study was to characterize the input signal of Mediterranean karst systems by analysing the chemical and isotopic compositions of rainfall, in particular with reference to spatial, intra-seasonal and inter-annual variations. Rainwater samples were collected and analysed for major ions, trace elements (B, Br and Sr) and stable water isotopes (δ18O, δ2H). The major results of this study are the following: the degree of fluctuation in concentrations of major elements and stable water isotopes in the rainwater depend strongly on elevation and distance from sea. Higher rainfall concentrations are systematically associated with areas nearest to the Mediterranean coast. Rainfall of Atlantic and Mediterranean origins may be differentiated – rainfall of Mediterranean origin is characterized by a higher deuterium excess and generally enriched in heavy isotopes. This rainfall is predominant during the autumn.

International audience; A monitoring of spring and rain waters in the South of France during two hydrological cycles is presented. Rain waters were sampled after each precipitation event at 3 rain-gauge stations. Four karstic springs located in the same area have also been daily (high discharge events) to monthly (low flow periods) sampled. This paper focuses on compositional changes in the Cl− and Br− ions and the oxygen-18 (δ18O) and hydrogen (δD) isotopes during the high discharge events. The responses of the different karstic systems are quite homogeneous and reflect the hydrological state of the system. The waters discharged during the major autumn and winter high discharge events originate from the epikarstic reservoir and show characteristics chemical variations related to residence in the unsaturated zone close to the surface. Their residence time in the order of 1-3 months. Correlations between the composition of the spring-water and the rainwater during three successive high discharge events during the summer of 1998 indicate that the water for high discharge event "n" is derived from water from precipitation event "n − 1" via a piston-type mechanism with residence time of 2 weeks. These results are interpreted as an indication of the major role of the epikarst reservoir in the karst recharge functioning. The similar behaviour of the four springs, although located in different geological contexts allows to think that the epikarst role could be more important than previously thought.