Jean-Luc Probst

The French OZCAR critical zone network offers the opportunity to conduct multi-site studies and to explore the critical zone functioning under contrasted climate, geology, vegetation and land use. In this study, an integrated modeling of the water cycle is performed with the ecohydrological model EcH2O-iso in three long-term observatories: (1) the Naizin watershed characterized by an oceanic climate, a metamorphic bedrock and an intensive agriculture (north-west of France, AgrHyS observatory); (2) the Aurade watershed, a watershed with a warmer semi-continental oceanic climate, a sedimentary geological substratum and a crop cover with a wheat-sunflower rotation (south-west of France, Aurade observatory) and; (3) the Strengbach watershed characterized by a mountain climate, a granitic bedrock, and a beech-spruce forest cover (north-east of France, OHGE observatory).Modeling robustness is evaluated by taking advantage of the large database for critical zone sciences including stream flow, water level in piezometers, and evapotranspiration fluxes measured from climatological stations and flux-towers located in the watersheds. Our comparative study brings these general outcomes: (1) the long term CZ evolution controlling the regolith thickness strongly impacts the total water storage in watersheds; (2) the Quaternary geomorphological evolution influences the current hydrological partitioning and the separation of hydrologically active and inactive water storage; (3) Both internal watershed characteristics and external forcings, such as current atmospheric forcing and recent land use need to be considered to infer stream persistence and to understand hydrological diversity; and (4) the observed hydrological diversity cannot be fully understood without considering a continuum of time scales in CZ evolution. Overall, this work illustrates the strength of critical zone networks, allowing a new level of multi-site and comparative studies that are crossing several observatories and encompassing a wide diversity of geology and climate.

An integrated ecohydrological modeling approach was deployed in three long‐term critical zone (CZ) observatories of the French CZ network (CZ Observatories—Application and Research) to better understand how the CZ heterogeneity modulates the water cycle within territories. Ecohydrological simulations with the physically based model EcH₂O‐iso constrained by a wide range of observations crossing several disciplines (meteorology, hydrology, geomorphology, geophysics, soil sciences, and satellite imagery) are able to capture stream water discharges, evapotranspiration fluxes, and piezometric levels in the Naizin, Auradé, and Strengbach watersheds. In Naizin, an agricultural watershed in northwestern France with a schist bedrock underlying deep weathered materials (5–15 m) along gentle slopes, modeling results reveal a deep aquifer with a large total water storage (1,080–1,150 mm), an important fraction of inactive water storage (94%), and relatively long stream water transit times (0.5–2.5 years). In the Auradé watershed, representative of agricultural landscapes of the southwestern France developed on molasse, a relatively shallow regolith (1–8 m) is observed along hilly slopes. Simulations indicate a shallow aquifer with moderate total water storage (590–630 mm), an important fraction of inactive water storage (91%), and shorter stream water transit times (0.1–1.3 years). In the Strengbach watershed, typical of mid‐mountain forested landscapes developed on granite, CZ evolution implies a shallow regolith (1–5 m) along steep slopes. Modeling results infer a shallow aquifer with the smallest total water storage (475–575 mm), the shortest stream water transit times (0.1–0.7 years), but also the highest fraction of active water storage (18%).

Local waterlogging often occurs on the steep slopes of clayey–calcareous soils in southwestern France, causing nutrients and pollutants transfer to the river bodies and reduced ecosystems services. These soils developed in the Miocene molassic hill formation and are generally impermeable with abundant traces of hydromorphy and heterogenous spatial distribution. This article aims to describe the hydrological functioning of these soils, based on a cross analysis of pedological, hydrological, and geophysical characterizations. Our experimental site is the catchment area located in Auradé (southwestern France). Here, we analyze the flows at the outlet of the studied watershed together with piezometric and climatic monitoring from September 2020 to September 2021. We show that the hydrological year is divided into three phases: first, a soil recharge phase with an effective rainfall of about 100 mm; second, a saturation phase, when 80% of the effective precipitation is drained mostly by runoff and hypodermic flows; third, a drying phase. Soil waterlogging events usually occur during the saturation phase. They are due to several forms of flow: surface runoff associated with return flow, hypodermic flow caused by the presence of soil layers with lower hydraulic conductivity in the subsurface (swelling clays and plowing sole) and groundwater flow with intermittent connection of the soil water table in the hillside to the alluvial groundwater table. We also conducted independent seismic refraction tomography analyses that validate localized waterlogging patterns along the catchment and open the way to spatializing areas with high waterlogging potential at the scale of the study plot.

Bed sediments were collected monthly at the outlet of the Ibrahim River between May 2016 and July 2017. Their physical and chemical characteristics were studied to highlight the impact of seasonal variations and discharge on the sediments’ geochemical signatures. Granulometric analysis showed that samples collected after rain events contained abundant sand fractions (> 80%), while the clay fraction was more present at the beginning of the dry season, at low monthly average flows (1 m³/s). Ten major elements, 14 rare earth elements (REE), and 30 trace elements (TE) were analyzed. An excess of CaO reflected the contribution of the carbonate rocks of the Ibrahim karstic springs (30.35 ± 3.91%) but CaO concentration decreased during periods of high water. On another hand, no REE enrichment was detected, both in high and low flows. Most sediments collected at low flow had a pronounced REE depletion that occurred particularly for LREE composition (Nd, Pr, Ce, and La). A negative Ce anomaly (0.992) and a positive Eu anomaly (1.313) were revealed with an average La/Yb ratio of 0.570, reflecting a slight enrichment in HREE. Most studied TE were less concentrated than the averages mentioned in PAAS, UCC, and WSA references. However, a Ze enrichment was mainly due to the regional geochemical background, an As enrichment was associated with anthropogenic contribution, and a Zr enrichment was linked to discharges from pharmaceutical industries located at the river outlet. This study still needs to be complemented binding both spatial and temporal criteria for further fluvial sediments’ monitoring of the entire catchment area.

Cet article met en évidence la répartition saisonnière des précipitations dans le bassin versant de Beni-Bahdel, et son impact sur les écoulements de surface et les transports solides en suspension. Pour ce faire, nous avons exploité les données météorologiques, les débits liquides et les concentrations de matières en suspension fournies par l’Agence Nationale des Ressources Hydriques (ANRH) sur une période de 46 ans (1972-2018). En considérant uniquement les écoulements de surface, nous avons établi les relations directes entre les précipitations et les débits liquides et entre les précipitations et les transports solides pour mettre en évidence l’impact saisonnier des précipitations sur l’érosion et les ruissellements de surface. Les résultats obtenus dans ce travail montrent que les pluies automnales et estivales provoquent directement l’érosion des versants induisant des concentrations en MES importantes à l’exutoire du bassin versant. En hiver et au printemps, elles contribuent aux écoulements de surface tout en apportant des flux d’eau significatifs au barrage. Les calculs basés sur une modélisation linéaire montrent qu’en moyenne 1 mm de pluie génère 0,031 Hm3 d’écoulement à l’exutoire du bassin et 14 x 103 tonnes de transports solides en automne, alors qu’au printemps, il provoque 0,063 Hm3 d’écoulement et 1,64 x 103 tonnes de matières en suspension.

This study simulates carbon dioxide (CO₂) sequestration in 300 major world river basins (about 70% of global surface area) through carbonates dissolution and silicate hydrolysis. For each river basin, the daily timescale impacts under the RCP 2.6 and RCP 8.5 climate scenarios were assessed relative to a historical baseline (1969–1999) using a cascade of models accounting for the hydrological evolution under climate change scenarios. Here we show that the global temporal evolution of the CO₂ uptake presents a general increase in the annual amount of CO₂ consumed from 0.247 ± 0.045 Pg C year⁻¹ to 0.261 and 0.273 ± 0.054 Pg C year⁻¹, respectively for RCP 2.6 and RCP 8.5. Despite showing a general increase in the global daily carbon sequestration, both climate scenarios show a decrease between June and August. Such projected changes have been mapped and evaluated against changes in hydrology, identifying hot spots and moments for the annual and seasonal periods.

River water quality is particularly of concern in semi-arid countries with limited water resources. Increasing anthropogenic activities can lead to the accumulation of trace metals (TM) in bottom sediments, which is a specific storage compartment. The present study aimed to investigate the geochemistry of trace metals (As, Cd, Co, Cr, Cu, Ni, Pb, Zn) and of some physico-chemical parameters in bottom sediments from the Sebou basin, which represents 1/3 of the surface water resources of Morocco. The order of abundance of the metals was Zn > Cr > Cu > Ni > Pb > Co > As > Cd. A major fingerprint of weathering on metal concentration, and point and non-point anthropogenic sources were highlighted. The origin and intensity of the contamination were determined using a combination of geochemical indicators. The contamination was on the whole moderate, with Cr, Zn, Cu, and Pb as the most enriched metals, especially at the A1, S3, and S4 stations located downstream of Fez city, well known for its intensive industrial and tannery activities. A multi-variate analysis evidenced the strong link between natural elements such as Co with clays and Fe oxides, and As with Ca, whereas Cd, Cu, Cr, Ni, Pb, Zn, partly originating from anthropogenic activities (industrial and domestic waste, agricultural inputs), were linked to phosphorus, oxides, carbonates, and/or POC, indicating their anthropic source and/or control by sediment compounds. Cadmium, Pb, and Cu were the most available metals. Finally, in addition to Cd, Pb and Zn were identified as hazardous metals in sediments as evidenced by the positive relationship between the proportion of the labile fraction and the enrichment factor revealing anthropogenic inputs.

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 (Ca²⁺) and bicarbonate (HCO₃⁻) which revealed an overall chemostatic behavior as a result of carbonate dissolution in the karst. On the contrary, sulfate (SO₄²⁻) exhibited the widest relative variation during flooding and showed a significant sensitivity to the dilution process with increasing discharge. For medium flood episodes (Q_max < 4.4 m³·s⁻¹), an overall concentration increase or chemostatic behavior could be observed during the rising limb of the hydrograph. On the opposite, under extreme flood episodes (Q_max > 8.3 m³·s⁻¹) 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.