Academic literature on the topic 'Lez karstic catchment'

Abstract. Rainfall-runoff models are crucial tools for the statistical prediction of flash floods and real-time forecasting. This paper focuses on a karstic basin in the South of France and proposes a distributed parsimonious event-based rainfall-runoff model, coherent with the poor knowledge of both evaporative and underground fluxes. The model combines a SCS runoff model and a Lag and Route routing model for each cell of a regular grid mesh. The efficiency of the model is discussed not only to satisfactorily simulate floods but also to get powerful relationships between the initial condition of the model and various predictors of the initial wetness state of the basin, such as the base flow, the Hu2 index from the Meteo-France SIM model and the piezometric levels of the aquifer. The advantage of using meteorological radar rainfall in flood modelling is also assessed. Model calibration proved to be satisfactory by using an hourly time step with Nash criterion values, ranging between 0.66 and 0.94 for eighteen of the twenty-one selected events. The radar rainfall inputs significantly improved the simulations or the assessment of the initial condition of the model for 5 events at the beginning of autumn, mostly in September–October (mean improvement of Nash is 0.09; correction in the initial condition ranges from −205 to 124 mm), but were less efficient for the events at the end of autumn. In this period, the weak vertical extension of the precipitation system and the low altitude of the 0 °C isotherm could affect the efficiency of radar measurements due to the distance between the basin and the radar (~60 km). The model initial condition S is correlated with the three tested predictors (R2 > 0.6). The interpretation of the model suggests that groundwater does not affect the first peaks of the flood, but can strongly impact subsequent peaks in the case of a multi-storm event. Because this kind of model is based on a limited amount of readily available data, it should be suitable for operational applications.

Abstract. The present study explores the application of a data assimilation (DA) procedure to correct the radar rainfall inputs of an event-based, distributed, parsimonious hydrological model. An extended Kalman filter algorithm was built on top of a rainfall-runoff model in order to assimilate discharge observations at the catchment outlet. This work focuses primarily on the uncertainty in the rainfall data and considers this as the principal source of error in the simulated discharges, neglecting simplifications in the hydrological model structure and poor knowledge of catchment physics. The study site is the 114 km2 Lez catchment near Montpellier, France. This catchment is subject to heavy orographic rainfall and characterised by a karstic geology, leading to flash flooding events. The hydrological model uses a derived version of the SCS method, combined with a Lag and Route transfer function. Because the radar rainfall input to the model depends on geographical features and cloud structures, it is particularly uncertain and results in significant errors in the simulated discharges. This study seeks to demonstrate that a simple DA algorithm is capable of rendering radar rainfall suitable for hydrological forecasting. To test this hypothesis, the DA analysis was applied to estimate a constant hyetograph correction to each of 19 flood events. The analysis was carried in two different modes: by assimilating observations at all available time steps, referred to here as reanalysis mode, and by using only observations up to 3 h before the flood peak to mimic an operational environment, referred to as pseudo-forecast mode. In reanalysis mode, the resulting correction of the radar rainfall data was then compared to the mean field bias (MFB), a corrective coefficient determined using rain gauge measurements. It was shown that the radar rainfall corrected using DA leads to improved discharge simulations and Nash-Sutcliffe efficiency criteria compared to the MFB correction. In pseudo-forecast mode, the reduction of the uncertainty in the rainfall data leads to a reduction of the error in the simulated discharge, but uncertainty from the model parameterisation diminishes data assimilation efficiency. While the DA algorithm used is this study is effective in correcting uncertain radar rainfall, model uncertainty remains an important challenge for flood forecasting within the Lez catchment.

Abstract. The present study explores the application of a data assimilation (DA) procedure to correct the radar rainfall inputs of an event-based, distributed, parsimonious hydrological model. A simplified Kalman filter algorithm was built on top of a rainfall-runoff model in order to assimilate discharge observations at the catchment outlet. The study site is the 114 km2 Lez Catchment near Montpellier, France. This catchment is subject to heavy orographic rainfall and characterized by a karstic geology, leading to flash flooding events. The hydrological model uses a derived version of the SCS method, combined with a Lag and Route transfer function. Because it depends on geographical features and cloud structures, the radar rainfall input to the model is particularily uncertain and results in significant errors in the simulated discharges. The DA analysis was applied to estimate a constant correction to each event hyetogram. The analysis was carried out for 19 events, in two different modes: re-analysis and pseudo-forecast. In both cases, it was shown that the reduction of the uncertainty in the rainfall data leads to a reduction of the error in the simulated discharge. The resulting correction of the radar rainfall data was then compared to the mean field bias (MFB), a corrective coefficient determined using ground rainfall measurements, which are more accurate than radar but have a decreased spatial resolution. It was shown that the radar rainfall corrected using DA leads to improved discharge simulations and Nash criteria compared to the MFB correction.

Hydrological variations hold a significant influence over the water chemistry in the karstic critical zone. In this context, the karstic Baget Catchment (BC) has been monitored at a high resolution over two years at the outlet in order to set up a typology of the flood events. The objective was to assess the multiple streamwater physico-chemical patterns in response to hydrological variations, streamflow component (quick-response, subsurface, and baseflow) and lithological contributions, and biogeochemical processes. The karstic catchment exhibited an impulsive response to flood events in relation to the typical structural and morphological characteristics of the karst. In addition, this response was constrained by the magnitude of the rainfall and the preceding hydroclimatic conditions. The variability of the dissolved load in streamflow was closely associated with the characteristics of the weathered rocks and the hydrological conditions throughout the year. Two simple indicators allow to characterize the concentration–discharge relationships with different hysteresis patterns on a set of floods with various intensities and shapes of the hydrograph and under different hydrological conditions before the flood. Almost all elements exhibited either clockwise loops or more complex behaviors, suggesting a higher overall concentration when the major water contribution comes from the quick-response flow (karst and surface runoff fraction). Besides, the epikarst flushing under dry conditions led counterclockwise hysteresis patterns for calcium (Ca2+) and bicarbonate (HCO3−) which revealed an overall chemostatic behavior as a result of carbonate dissolution in the karst. On the contrary, sulfate (SO42−) exhibited the widest relative variation during flooding and showed a significant sensitivity to the dilution process with increasing discharge. For medium flood episodes (Qmax < 4.4 m3·s−1), an overall concentration increase or chemostatic behavior could be observed during the rising limb of the hydrograph. On the opposite, under extreme flood episodes (Qmax > 8.3 m3·s−1) occurring after several rain events, a dilution pattern was noted for all elements originating from rock weathering. Finally, high-frequency sampling during storm events improved the understanding of the factors controlling the hydrochemical dynamic in karstic catchments.

Abstract. Lowland karst aquifers can generate unique wetland ecosystems which are caused by groundwater fluctuations that result in extensive groundwater–surface water interactions (i.e. flooding). However, the complex hydrogeological attributes of these systems, linked to extremely fast aquifer recharge processes and flow through well-connected conduit networks, often present difficulty in predicting how they will respond to changing climatological conditions. This study investigates the predicted impacts of climate change on a lowland karst catchment by using a semi-distributed pipe network model of the karst aquifer populated with output from the high spatial resolution (4 km) Consortium for Small-scale Modelling Climate Lokalmodell (COSMO-CLM) regional climate model simulations for Ireland. An ensemble of projections for the future Irish climate were generated by downscaling from five different global climate models (GCMs), each based on four Representative Concentration Pathways (RCPs; RCP2.6, RCP4.5, RCP6.0 and RCP8.5) to account for the uncertainty in the estimation of future global emissions of greenhouse gases. The one-dimensional hydraulic/hydrologic karst model incorporates urban drainage software to simulate open channel and pressurised flow within the conduits, with flooding on the land surface represented by storage nodes with the same stage volume properties of the physical turlough basins. The lowland karst limestone catchment is located on the west coast of Ireland and is characterised by a well-developed conduit-dominated karst aquifer which discharges to the sea via intertidal and submarine springs. Annual above ground flooding associated with this complex karst system has led to the development of unique wetland ecosystems in the form of ephemeral lakes known as turloughs; however, extreme flooding of these features causes widespread damage and disruption in the catchment. This analysis has shown that mean, 95th and 99th percentile flood levels are expected to increase by significant proportions for all future emission scenarios. The frequency of events currently considered to be extreme is predicted to increase, indicating that more significant groundwater flooding events seem likely to become far more common. The depth and duration of flooding is of extreme importance, both from an ecological perspective in terms of wetland species distribution and for extreme flooding in terms of the disruption to homes, transport links and agricultural land inundated by flood waters. The seasonality of annual flooding is also predicted to shift later in the flooding season, which could have consequences in terms of ecology and land use in the catchment. The investigation of increasing mean sea levels, however, showed that anticipated rises would have very little impact on groundwater flooding due to the marginal impact on ebb tide outflow volumes. Overall, this study highlights the relative vulnerability of lowland karst systems to future changing climate conditions, mainly due to the extremely fast recharge which can occur in such systems. The study presents a novel and highly effective methodology for studying the impact of climate change in lowland karst systems by coupling karst hydrogeological models with the output from high-resolution climate simulations.

Excess nitrogen (N) deposition and gaseous N emissions from industrial, domestic, and agricultural sources have led to increased nitrate leaching, the loss of biological diversity, and has affected carbon (C) sequestration in forest ecosystems. Nitrate leaching affects the purity of karst water resources, which contribute around 50% to Austria’s drinking water supply. Here we present an evaluation of the drivers of dissolved inorganic N (DIN) concentrations and fluxes from a karst catchment in the Austrian Alps (LTER Zöbelboden) from 27 years of records. In addition, a hydrological model was used together with climatic scenario data to predict expected future runoff dynamics. The study area was exposed to increasing N deposition during the 20th century (up to 30 to 35 kg N ha−1 y−1), which are still at levels of 25.5 ± 3.6 and 19.9 ± 4.2 kg N ha−1 y−1 in the spruce and the mixed deciduous forests, respectively. Albeit N deposition was close to or exceeded critical loads for several decades, 70–83% of the inorganic N retained in the catchment from 2000 to 2018, and NO3- concentrations in the runoff stayed <10 mg L−1 unless high-flow events occurred or forest stand-replacing disturbances. We identified tree growth as the main sink for inorganic N, which might together with lower runoff, increase retention of only weakly decreasing N deposition in the future. However, since recurring forest stand-replacement is predicted in the future as a result of a combination of climatically driven disturbance agents, pulses of elevated nitrate concentrations in the catchment runoff will likely add to groundwater pollution.

Because spring water quality from alpine karst aquifers can change very rapidly during event situations, water abstraction management has to be performed in near real-time. Four summer events (2005–2008) at alpine karst springs were investigated in detail in order to evaluate the spectral absorption coefficient at 254 nm (SAC254) as a real-time early warning proxy for faecal pollution. For the investigation Low-Earth-Orbit (LEO) Satellite-based data communication between portable hydrometeorological measuring stations and an automated microbiological sampling device was used. The method for event triggered microbial sampling and analyzing was already established and described in a previous paper. Data analysis including on-line event characterisation (i.e. precipitation, discharge, turbidity, SAC254) and comprehensive E. coli determination (n&gt;800) indicated that SAC254 is a useful early warning proxy. Irrespective of the studied event situations SAC254 always increased 3 to 6 hours earlier than the onset of faecal pollution, featuring different correlation phases. Furthermore, it seems also possible to use SAC254 as a real-time proxy parameter for estimating the extent of faecal pollution after establishing specific spring and event-type calibrations that take into consideration the variability of the occurrence and the transferability of faecal material It should be highlighted that diffuse faecal pollution from wildlife and live stock sources was responsible for spring water contamination at the investigated catchments. In this respect, the SAC254 can also provide useful information to support microbial source tracking efforts where different situations of infiltration have to be investigated.

The chemical composition of throughfall depends on the age of the Norway spruce (Picea abies Karst) stands and season of the year. The pH of throughfall decreased and the amount of hydrogen ion in throughfall deposited to the soil increased with increasing age of spruce stands, especially in the winter season. Concentrations of K+, H+, SO42−, Mn2+, and NH4+in throughfall were higher than bulk precipitation for the whole year and K+, H+, and Mn2+concentrations were higher in throughfall in winter and the growing season. This indicates that these ions were washed out or washed from the surface of needles and/or the bark, and that NO3−, NH4+, Ca2+, Mg2+, Fe2+, and Zn2+were absorbed in the canopy. The effect of high nitrogen deposition, above critical loads, and an increase in the amount of sulfur and in the sum of the strong acids (S-SO42−and N-NO3−) that reached the soil with throughfall may have implications for the vitality of spruce stands, especially in older age classes. The application of Principal Component Analysis (PCA) has led to identification of five factors responsible for the data structure (“mineral dust”, “acidic emissions”, “heavy metals-dust particles”, “ammonium [NH4+]”, and “H+”). They explain more than 60% of the total variance system. The strong positive correlation between stand age class and ionic concentrations in throughfall occurs for all year and the winter period for ions within the following categories: “acidic emissions”, SO42−+ NO3−; “heavy metals-dust particles”, Fe2++ Mn2++ Zn2+; “mineral dust”, Na++ K++ Ca+2+ Mg2+; “NH4+”; and “H+”. The strength of the relationship decreases in the growing period, probably due to processes occurring in the canopy (adsorption, leaching, etc.).

Les crues « éclair » parfois dévastatrices qui touchent les bassins versants méditerranéens du Sud de la France sont difficiles à anticiper. Leur prévision passe par l'utilisation de modèles pluie-débit, dont l'efficacité est encore limitée par les incertitudes liées notamment à la variabilité spatiale des pluies méditerranéennes et à la caractérisation de l'état hydrique initial des hydrosystèmes. Dans le cas de bassins karstiques, à ces incertitudes s'ajoutent celles liées à la dynamique des aquifères et à leur rôle sur la formation des crues. La première partie de ce travail de thèse propose un modèle pluie-débit horaire, distribué, événementiel et parcimonieux pour reproduire les crues « éclair » à l'exutoire du bassin karstique du Lez (Montpellier) de 114 km2. Le modèle est évalué non seulement sur la qualité des simulations de débits mais aussi sur la qualité de son initialisation obtenu grâce à une relation entre sa condition initiale et divers indicateurs de l'état hydrique de l'hydrosystème. Calibré sur 21 épisodes de crues, le modèle fournit des simulations satisfaisantes, et sa condition initiale est significativement corrélée à l'indice d'humidité Hu2 du modèle SIM de Météo-France ou à la piézométrie dans l'aquifère du Lez. Les pluies mesurées par radar en début d'automne sont de bonne qualité et conduisent à une amélioration des simulations de débit et de l'estimation de la condition initiale du modèle. En revanche, les pluies mesurées par radar en fin d'automne sont de moindre qualité et n'améliorent pas les simulations. Face aux incertitudes liées à la paramétrisation du modèle ou à l'estimation des pluies radar, la deuxième partie du travail de thèse analyse l'apport de l'assimilation des débits observés pour corriger en temps réel les paramètres les plus sensibles du modèle et notamment sa condition initiale ou les pluies radar en entrée du modèle. La procédure d'assimilation de données a été mise en place à l'aide du coupleur PALM, qui permet de relier modèle hydrologique à l'algorithme d'assimilation. La correction de la condition initiale du modèle permet généralement d'améliorer les prévisions (sous hypothèse de pluie future connue); la correction de la pluie a des effets similaires. Néanmoins les limites de cette correction sont atteintes dans le cas où le modèle ne reproduit pas de façon satisfaisante la partie initiale de montée des eaux, ce qui pourra être amélioré par la suite. Finalement, ce travail de thèse montre que la complexité d'un bassin karstique peut être représentée efficacement à l'aide d'un nombre réduit de paramètres, pour simuler les débits, et contribue à l'amélioration des outils opérationnels pour la prévision des crues
The sometimes devastating flash floods which affect the Mediterranean watersheds of the South of France are difficult to anticipate. Flood forecasting requires the use of rainfall-runoff models which are limited in their efficiency by uncertainty related to the spatial variability of Mediterranean rainfall and the characterization of the initial hydric state of the system. In karstic catchments, these uncertainties are added to those due to aquifer dynamics and their role in flood genesis. The first part of this work will present a distributed event-based parsimonious hourly rainfall-runoff model in order to reconstruct flash flood events at the outlet of the 114 km2 Lez Catchment (Montpellier). The model is evaluated not only for the quality of the simulations produced, but for the quality of its parameter initialization obtained using a relationship between the initial condition and various hydric state indicators of the system. Calibrated using 21 flood episodes, the model produces satisfactory simulations and its initial condition is significantly correlated with the Hu2 soil humidity index of the Météo-France model or piezometers measuring the Lez aquifer. Radar rainfall data measured in early fall are of good quality and lead to improved discharge simulations and an improved estimation of the model initial condition. However, rainfall measured by radar in late fall are of poor quality and do not improve the simulations. Confronted with the uncertainty related to model parametrization or the estimation of radar rainfall, the second part of this dissertation analyzes improvements achieved by assimilating observed discharge measurements in order to perform real-time corrections to the most sensitive model parameters and notably the initial condition and the radar rainfall input to the model. The data assimilation procedure was implemented with the help of the PALM coupling software which allows for the linking of the hydrological model with the assimilation algorithm. Correcting the initial condition allowed for, on average, the improvement of forecasting (under a known future rainfall hypothesis); correcting the rainfall had similar effects. Nevertheless, the limits of this approach are reached when the model is unable to satisfactorily reproduce the rising limb of the hydrograph, a problem which may be addressed by future research. Finally, this body of work demonstrates that the complexity of a karstic catchment can be efficiently represented with a reduced number of parameters in order to simulate discharges and contribute to the improvement of operational tools for flood forecasting

L'hétérogénéité spatiale est un trait caractéristique des aquifères calcaires. Conditionnée par le karst et ses morphologies, cette hétérogénéité explique la fréquente sous-exploitation des réserves de ces hydrosystèmes mais également le degré lacunaire de connaissance de leurs fonctionnements. Nous avons voulu, dans cette recherche, approfondir cette connaissance en scrutant le fonctionnement des aquifères calcaires qui arment les bas plateaux du Nord-Est de la France. Trois de ces bas plateaux ont été retenus, ceux du Dogger et de l'Oxfordien (Jurassique) et celui du Muschelkalk (Trias). Afin de comprendre comment s'organisent les écoulements sur et dans ces trois systèmes, notre démarche s'articule autour de deux principaux axes ; le premier concerne l'identification des modalités de l'écoulement et le second est consacré à leur quantification. En pratique, la connaissance des modalités de circulation a été recherchée à travers l'étude des différents paramètres structurels (géologique, tectonique, karstique) et conjoncturels (précipitation, évapotranspiration) qui interagissent dans le cycle de l'eau. Nous montrons que la variété de fonctionnements hydrologiques observés, à l'échelle du domaine d'étude, est conditionnée principalement par les paramètres structurels ; l'hétérogénéité lithologique, notamment, apparaît comme une clef fonctionnelle déterminante, d'une part, par l'existence d'une multiplicité de nappes d'eau et d'autre part, par son rôle dans la karstogenèse. Les morphologies du karst s'observent principalement sur les secteurs de contact lithostratigraphique et nous avons pu établir, à travers une répartition spatiale des phénomènes du karst (en surface et en souterrain), que ces trois plateaux sont majoritairement déconnectés du karst et dominés par des processus de circulation lents, de l'infiltration à la vidange aquifère. Afin d'apporter des éléments de quantification sur les circulations, une analyse hydroclimatique et géochimique d'exutoires de petites dimensions a été ensuite menée. C'est ainsi qu'une petite dizaine de sources et de bassins versants, localisés dans des contextes physiographiques variés, ont fait l'objet de suivis sur plusieurs cycles hydrologiques. L'analyse nous a permis de préciser les caractères hydrodynamiques de deux grands types de comportement observés au sein des bas plateaux étudiés : le comportement fissural et le comportement karstique. Le comportement fissural se manifeste, à l'échelle annuelle, par un filtrage important des précipitations par les systèmes (forte inertie des valeurs de débits et de conductivités électriques) et également par de lents processus de mise en charge à l'échelle des évènements de crue. Le comportement karstique, quant à lui, est attesté, à l'échelle annuelle, par des pluies moyennement filtrées par les systèmes mais aussi par des variations hydrogéochimiques accusées lors des sollicitations pluviométriques majeures, en lien avec la transmissivité importante de ces milieux
The spatial heterogeneity is a feature of limestone aquifers. Conditioned by the karst and its morphologies, this heterogeneity explains the frequent under-exploitation of these hydrosystem reserves but it also explains the incomplete degree of knowledge of their functionings. In this research we wanted to deepen this knowledge by scrutinizing the functioning of the limestone aquifers which arm the low plateaus of the North-East of France. Three of these low plateaus held our attention, those of Dogger and Oxfordian (Jurassic) and the one of Muschelkalk (Triassic). In order to understand how the flows on and into these three systems are organized, our approach focuses on two main axes ; the first one concerns the identification of the modes of flow and the second one is devoted to quantification. In practice, the knowledge of the flow modes was investigated through the study of various structural parameters (geology, tectonics, karst) and conjunctural parameters (precipitation, evapotranspiration) interacting the water cycle. We show that the variety of hydrological functionings which are observed, on the scale of the study area, is mainly dependent on the structural parameters ; the lithological heterogeneity, in particular, appears as a determining functional key, on the one hand, by the existence of a multiplicity of groundwater and on the other hand, by its role in the karstogenesis. The morphologies of the karst are primarily obvious on the areas of lithostratigraphic contact and we were able to establish, through a spatial distribution of karst phenomena (on and underground), that these three plateaus are mainly disconnected from the karst and run by slow-flow processes, from infiltration to ground water discharge. In order to provide some quantifications of the circulations, a hydroclimatic and geochemical analysis of small outlets was then conducted. Thus a dozen springs and catchments, located in various physiographic settings have been surveyed for several water years. The analysis allowed us to specify the hydrodynamic characteristics of two types of behavior observed in the low plateaus we studied: the fissural behavior and the karstic behavior. The fissural behavior appears, on an annual scale, through a significant rainfall filtering of the systems (high inertia of the values of discharges and electrical conductivities) and also by slow charging processes when flood events occur. The karst behavior, meanwhile, is evidenced on an annual scale, by a moderately rainy filtering of the systems but also by hydrogeochemical variations accused during major rainfall solicitations in connection with the important transmissivity of these environments