Unexpected spatial stability of water chemistry in headwater stream networks.

Understanding how water and solutes enter and propagate through freshwater landscapes in the Anthropocene is critical to protecting and restoring aquatic ecosystems and ensuring human water security. However, high hydrochemical variability in headwater streams, where most carbon and nutrients enter river networks, has hindered effective modelling and management. We developed an analytical framework informed by landscape ecology and catchment hydrology to quantify spatiotemporal variability across scales, which we tested in 56 headwater catchments, sampled periodically over 12 years in western France. Unexpectedly, temporal variability in dissolved carbon, nutrients and major ions was preserved moving downstream and spatial patterns of water chemistry were stable on annual to decadal timescales, partly because of synchronous variation in solute concentrations. These findings suggest that while concentration and flux cannot be extrapolated among subcatchments, periodic sampling of headwaters provides valuable information about solute sources and subcatchment resilience to disturbance.

Flow intermittency, dispersion, and correlated continuous time random walks in porous media.

We study the intermittency of fluid velocities in porous media and its relation to anomalous dispersion. Lagrangian velocities measured at equidistant points along streamlines are shown to form a spatial Markov process. As a consequence of this remarkable property, the dispersion of fluid particles can be described by a continuous time random walk with correlated temporal increments. This new dynamical picture of intermittency provides a direct link between the microscale flow, its intermittent properties, and non-Fickian dispersion.

Anomalous kinetics in diffusion limited reactions linked to non-Gaussian concentration probability distribution function.

We investigate anomalous reaction kinetics related to segregation in the one-dimensional reaction-diffusion system A + B → C. It is well known that spatial fluctuations in the species concentrations cause a breakdown of the mean-field behavior at low concentration values. The scaling of the average concentration with time changes from the mean-field t(-1) to the anomalous t(-1/4) behavior. Using a stochastic modeling approach, the reaction-diffusion system can be fully characterized by the multi-point probability distribution function (PDF) of the species concentrations. Its evolution is governed by a Fokker-Planck equation with moving boundaries, which are determined by the positivity of the species concentrations. The concentration PDF is in general non-Gaussian. As long as the concentration fluctuations are small compared to the mean, the PDF can be approximated by a Gaussian distribution. This behavior breaks down in the fluctuation dominated regime, for which anomalous reaction kinetics are observed. We show that the transition from mean field to anomalous reaction kinetics is intimately linked to the evolution of the concentration PDF from a Gaussian to non-Gaussian shape. This establishes a direct relationship between anomalous reaction kinetics, incomplete mixing and the non-Gaussian nature of the concentration PDF.

Graph-based flow modeling approach adapted to multiscale discrete-fracture-network models.

Fractured rocks are often modeled as multiscale populations of interconnected discrete fractures (discrete fracture network, DFN). Graph representations of DFNs reduce their complexity to their connectivity structure by forming an assembly of nodes connected by links (edges) to which physical properties, like a conductance, can be assigned. In this study, we address the issue of using graphs to predict flow as a fast and relevant substitute to classical DFNs. We consider two types of graphs, whether the nodes represent the fractures (fracture graph) or the intersections between fractures (intersection graph). We introduce an edge conductance expression that accounts for both the portion of the fracture surface that carries flow and fracture transmissivity. We find that including the fracture size in the expression improves the prediction of flow compared to expressions used in previous studies that did not. The two graph types yield very different results. The fracture graph systematically underestimates local flow values. In contrast, the intersection graph overestimates the flow in each fracture because of the connectivity redundancy in fractures with multiple intersections. We address the latter issue by introducing a correction factor into the conductances based on the number of intersections on each fracture. We test the robustness of the proposed conductance model by comparing flow properties of the graph with high-fidelity DFN simulations over a wide range of network types. The good agreement found between the intersection graph and test suite indicates that this representation could be useful for predictive purposes.

Landscape dynamics revealed by luminescence signals of feldspars from fluvial terraces.

Luminescence signals of quartz and feldspar minerals are widely used to determine the burial age of Quaternary sediments. Although luminescence signals bleach rapidly with sunlight exposure, incomplete bleaching may affect luminescence ages, in particular in fluvial settings where an unbleached remnant signal is commonly encountered in modern alluvium. Here, we use feldspar single-grain post-infrared IR stimulation (pIRIR) dating to show that recent (<11 ka) fluvial terraces of the Rangitikei River (New Zealand) were formed in a context of non-linear incision rate. We relate this pattern to the rapid reinstatement of steady-state incision following the formation of a major, climate-driven, aggradation terrace, causing a phase of accelerated incision. In addition, we show systematic variations in the proportion of unbleached grains in the fluvial sediments over time, mirroring incision rate at the time of deposition. Deposits formed during rapid incision contain fewer bleached grains, which we attribute to large input of unbleached material and limited bleaching opportunities during fluvial transport. This finding demonstrates that the luminescence signals recorded in fluvial terraces not only yield age information, but also inform us on past fluvial transport and ultimately, landscape dynamics.

General database for ground water site information.

In most cases, analysis and modeling of flow and transport dynamics in ground water systems require long-term, high-quality, and multisource data sets. This paper discusses the structure of a multisite database (the H+ database) developed within the scope of the ERO program (French Environmental Research Observatory, http://www.ore.fr). The database provides an interface between field experimentalists and modelers, which can be used on a daily basis. The database structure enables the storage of a large number of data and data types collected from a given site or multiple-site network. The database is well suited to the integration, backup, and retrieval of data for flow and transport modeling in heterogeneous aquifers. It relies on the definition of standards and uses a templated structure, such that any type of geolocalized data obtained from wells, hydrological stations, and meteorological stations can be handled. New types of platforms other than wells, hydrological stations, and meteorological stations, and new types of experiments and/or parameters could easily be added without modifying the database structure. Thus, we propose that the database structure could be used as a template for designing databases for complex sites. An example application is the H+ database, which gathers data collected from a network of hydrogeological sites associated with the French Environmental Research Observatory.

Impact of climate changes during the last 5 million years on groundwater in basement aquifers.

Climate change is thought to have major effects on groundwater resources. There is however a limited knowledge of the impacts of past climate changes such as warm or glacial periods on groundwater although marine or glacial fluids may have circulated in basements during these periods. Geochemical investigations of groundwater at shallow depth (80-400 m) in the Armorican basement (western France) revealed three major phases of evolution: (1) Mio-Pliocene transgressions led to marine water introduction in the whole rock porosity through density and then diffusion processes, (2) intensive and rapid recharge after the glacial maximum down to several hundred meters depths, (3) a present-day regime of groundwater circulation limited to shallow depth. This work identifies important constraints regarding the mechanisms responsible for both marine and glacial fluid migrations and their preservation within a basement. It defines the first clear time scales of these processes and thus provides a unique case for understanding the effects of climate changes on hydrogeology in basements. It reveals that glacial water is supplied in significant amounts to deep aquifers even in permafrosted zones. It also emphasizes the vulnerability of modern groundwater hydrosystems to climate change as groundwater active aquifers is restricted to shallow depths.

Anomalous diffusion exponents in continuous two-dimensional multifractal media.

We study diffusion in heterogeneous multifractal continuous media that are characterized by the second-order dimension of the multifractal spectrum D2, while the fractal dimension of order 0, D0, is equal to the embedding Euclidean dimension 2. We find that the mean anomalous and fracton dimensions, d(w) and d(s), are equal to those of homogeneous media showing that, on average, the key parameter is the fractal dimension of order 0 D0, equal to the Euclidean dimension and not to the correlation dimension D2. Beyond their average, the anomalous diffusion and fracton exponents, d(w) and d(s), are highly variable and consistently range in the interval [1,4]. d(w) can be consistently either larger or lower than 2, indicating possible subdiffusive and superdiffusive regimes. On a realization basis, we show that the exponent variability is related to the local conductivity at the medium inlet through the conductivity scaling.

Persistence of incomplete mixing: a key to anomalous transport.

Anomalous dispersion in heterogeneous environments describes the anomalous growth of the macroscopic characteristic sizes of scalar fields. Here we show that this phenomenon is closely related to the persistence of local scale incomplete mixing. We introduce the mixing scale ε as the length for which the scalar distribution is locally uniform. We quantify its temporal evolution due to the competition of shear action and diffusion and compare it to the evolution of the global dispersion scale σ. In highly heterogeneous flow fields, for which the temporal evolution of σ is superdiffusive, we find that ε evolves subdiffusively. The anomalous evolutions of the dispersion and mixing scales are complementary, εσ ∝ t. This result relates anomalous global dispersion to the dynamics of local mixing.