Abstract. Understanding the variability of the chemical composition of surface waters
is a major issue for the scientific community. To date, the study of
concentration–discharge relations has been intensively used to assess the
spatiotemporal variability of the water chemistry at watershed scales.
However, the lack of independent estimations of the water transit times
within catchments limits the ability to model and predict the water
chemistry with only geochemical approaches. In this study, a dimensionally
reduced hydrological model coupling surface flow with subsurface flow (i.e.,
the Normally Integrated Hydrological Model, NIHM) has been used to constrain
the distribution of the flow lines in a headwater catchment (Strengbach
watershed, France). Then, hydrogeochemical simulations with the code KIRMAT
(i.e., KInectic Reaction and MAss Transport) are performed to calculate the
evolution of the water chemistry along the flow lines. Concentrations of
dissolved silica (H4SiO4) and in basic cations (Na+, K+,
Mg2+, and Ca2+) in the spring and piezometer waters are correctly
reproduced with a simple integration along the flow lines. The seasonal
variability of hydraulic conductivities along the slopes is a key process to
understand the dynamics of flow lines and the changes of water transit times
in the watershed. The covariation between flow velocities and active lengths
of flow lines under changing hydrological conditions reduces the variability
of water transit times and explains why transit times span much narrower
variation ranges than the water discharges in the Strengbach catchment.
These findings demonstrate that the general chemostatic behavior of the
water chemistry is a direct consequence of the strong hydrological control
of the water transit times within the catchment. Our results also show that
a better knowledge of the relations between concentration and mean transit time (C–MTT
relations) is an interesting new step to understand the diversity of C–Q
shapes for chemical elements. The good match between the measured and modeled concentrations while respecting the water–rock interaction times
provided by the hydrological simulations also shows that it is possible to
capture the chemical composition of waters using simply determined reactive
surfaces and experimental kinetic constants. The results of our simulations
also strengthen the idea that the low surfaces calculated from the
geometrical shapes of primary minerals are a good estimate of the reactive
surfaces within the environment.
Concentration-discharge (C-Q) relations can provide insight into the dynamic behavior of the Critical Zone (CZ), as C-Q relations integrate the spatial distribution and timing of watershed hydrogeochemical processes. This study blends geomorphologic analysis, C-Q relations and reactive-transport modeling using a rich dataset from an elevation gradient of eight watersheds in the Southern Sierra Nevada, California. We found that the CZ structure exerts a strong control on the C-Q relations, and on the hydrogeochemical behavior of headwater watersheds. Watersheds with thin regolith, a large stream network, and limited water storage have fast mean transit times along subsurface flow lines, and show limited seasonal variability in ionic concentrations in streamflow (i.e., chemostatic behavior). In contrast, watersheds with thicker regolith, a small stream network and more water storage have longer transit times along subsurface flow lines, and exhibit greater chemical variability (i.e., chemodynamic behavior). Independent estimates of mean transit times and water storage from other isotopic, hydrologic and geophysical studies were consistent with results from modeling C-Q relations. The stream chemistry and its variability were controlled by lateral flow within the regolith, and no mixing with deep groundwater was needed to explain the observed chemical variability. This study opens the possibility to estimate water-storage capacity and mean transit times, and thus drought resistance in watersheds, by using quantitative modeling of C-Q relations.
Abstract. Understanding the spatiotemporal variability of the chemical composition of surface waters is a major issue for the scientific community, especially given the prospect of significant environmental changes for the next decades. To date, the study of concentration-discharge relationships has been intensively used to assess the spatiotemporal variability of the water chemistry at watershed scales; however, the lack of independent estimations of the water transit times within catchments limits our ability to model and predict the water chemistry with only geochemical approaches. This study demonstrates the potential of coupling mathematical hydrology with hydrogeochemical modeling to better understand the spatiotemporal variability of the composition of surface waters. In a first step, a dimensionally reduced hydrological model coupling surface flow with subsurface flow (i.e., the Normally Integrated Hydrological Model, NIHM) has been used to constrain the distribution of the flow lines that are feeding the springs. In a second step, hydrogeochemical simulations with the code KIRMAT (KInectic Reaction and MAss Transport) have been performed to calculate the evolution of the water chemistry along the flow lines. The results indicate that the concentrations of dissolved silica (H4SiO4) and in basic cations (Na+, K+, Mg2+, and Ca2+) in the spring waters are correctly reproduced with a simple integration along the flow lines. The results also show that the modest variabilities of the flow line distribution and of the flow velocity imply that the water transit times only vary from approximately 1.5 to 3 months from floods to drought events. These findings demonstrate that the chemostatic behavior of the spring chemistry is a direct consequence of the strong hydrological control of the water transit times within the catchment. The good matching between the measured and modeled concentrations while respecting the water-rock interaction times provided by the hydrological simulations also shows that it is possible to capture the chemical composition of waters using simply determined reactive surfaces and standard kinetic constants. The results of our simulation strengthen the idea that the low surfaces calculated from the geometrical shapes of minerals are a good estimate of the reactive surfaces within the natural environment and certainly the one to be used for hydrogeochemical modeling such as that performed in this work, in addition to the use of the experimental kinetic constants for mineral dissolution. Overall, this work shows that the hydrogeochemical functioning of an elementary watershed, such as the Strengbach catchment, is relatively simple. The acquisition and variability of the water chemistry can be explained through process-based modeling approaches and by only formulating few hypotheses on the functioning of the watershed.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.