Regional estimation of catchment‐scale soil properties by means of streamflow recession analysis for use in distributed hydrological models

The estimation of catchment‐scale soil properties, such as water storage capacity and hydraulic conductivity, is of primary interest for the implementation of distributed hydrological models at the regional scale. This estimation is generally performed on the basis of information provided by soil databases. However, such databases are often established for agronomic uses and generally do not document deep‐weathered rock horizons (i.e. pedologic horizons of type C and deeper), which can play a major role in water transfer and storages. Here, we define the Drainable Storage Capacity Index (DSCI), an indicator that relies on the comparison between cumulated streamflow and precipitation to assess catchment‐scale storage capacities. DSCI is found to be reliable to detect underestimation of soil storage capacities in soil databases. We also use the streamflow recession analysis methodology defined by Brutsaert and Nieber in 1977 to estimate water storage capacities and lateral saturated hydraulic conductivities of the nondocumented deep horizons. The analysis is applied to a sample of 23 catchments (0.2–291 km²) located in the Cévennes‐Vivarais region (south of France). For regionalization purposes, the obtained results are compared with the dominant catchment geology and present a clear hierarchy between the different geologies of the area. Hard crystalline rocks are found to be associated with the thickest and less conductive deep soil horizons. Schist rocks present intermediate values of thickness and of saturated hydraulic conductivity, whereas sedimentary rocks and alluvium are found to be less thick and most conductive. These results are of primary interest in view of the future set‐up of distributed hydrological models over the Cévennes‐Vivarais region. Copyright © 2013 John Wiley & Sons, Ltd.     

Combining isotopic tracers (222Rn and δ13C) for improved modelling of groundwater discharge to small rivers

In regions where aquifers sustain rivers, the location and quantification of groundwater discharge to surface water are important to prevent pollution hazards, to quantify and predict low flows and to manage water supplies. ²²²Rn is commonly used to determine groundwater discharge to rivers. However, using this isotopic tracer is challenging because of the high diffusion capacity of ²²²Rn in open water. This study illustrates how a combination of isotopic tracers can contribute to an enhanced understanding of groundwater discharge patterns in small rivers. The aim of this paper is to combine ²²²Rn and δ¹³C_DIC to better constrain the physical parameters related to the degassing process of these tracers in rivers. The Hallue River (northern France) was targeted for this study because it is sustained almost exclusively by a fractured chalk aquifer. The isotopes ²²²Rn, δ¹³C_DIC, δ²H and δ¹⁸O were analysed along with other natural geochemical tracers. A mass balance model was used to simulate ²²²Rn and δ¹³C_DIC. The results of δ²H and δ¹⁸O analyses prove that evaporation did not occur in the river. The calibration of a numerical model to reproduce ²²²Rn and δ¹³C_DIC provides a best‐fit diffusive layer thickness of 3.21 × 10^?5 m. This approach is particularly useful for small rivers flowing over carbonate aquifers with high groundwater DIC where the evolution of river DIC reflects the competing processes of groundwater inflow and CO₂ degassing. This approach provides a means to evaluate groundwater discharge in small ungauged rivers. Copyright © 2014 John Wiley & Sons, Ltd.     

Exotic pollen as an indicator of variable atmospheric circulation over the Labrador Sea region during the mid to late Holocene

Variability in the abundance of exotic (non‐native) pollen in sediment cores has long been considered as a potential proxy for changing atmospheric circulation, but the difficulty of gaining sufficient total exotic pollen and the incomplete understanding of atmospheric pollen transport patterns has hindered its application. In light of recent advances in the study of pollen transport, we present an exotic pollen record from two fjord sediment cores taken from the west (Placentia Bay, Newfoundland) and east (Narsaq Sund, Greenland) Labrador Sea as a basis for studying variations in regional atmospheric circulation. The two cores cover the last ca. 5500 years and indicate a shift in dominant spring/summer air masses at ca. 2000 (southern Greenland) and 3000 cal a BP (Newfoundland) transporting reduced concentrations of pollen from southerly and south‐westerly vegetation zones. This may suggest a shift away from more dominantly zonal atmospheric circulation (a feature of positive North Atlantic Oscillation years) to more frequent meridional circulation. These results support sea ice/sea‐surface temperature proxy reconstructions from Newfoundland, investigated as part of the same project, which also suggest increased winter atmospheric circulation during the early part of the time period studied. In this region, more positive North Atlantic Oscillation years, and therefore more zonal atmospheric circulation, are associated with increased atmospheric circulation in both the winter and the summer seasons. Copyright © 2011 John Wiley & Sons, Ltd.     

A review of thermal response test analysis using pumping test concepts

The design of ground-coupled heat pump systems requires knowledge of the thermal properties of the subsurface and boreholes. These properties can be measured with in situ thermal response tests (TRT), where a heat transfer fluid flowing in a ground heat exchanger is heated with an electric element and the resulting temperature perturbation is monitored. These tests are analogous to standard pumping tests conducted in hydrogeology, because a system that is initially assumed at equilibrium is perturbed and the response is monitored in time, to assess the system's properties with inverse modeling. Although pumping test analysis is a mature topic in hydrogeology, the current analysis of temperature measurements in the context of TRTs is comparatively a new topic and it could benefit from the application of concepts related to pumping tests. The purpose of this work is to review the methodology of TRTs and improve their analysis using pumping test concepts, such as the well function, the superposition principle, and the radius of influence. The improvements are demonstrated with three TRTs. The first test was conducted in unsaturated waste rock at an active mine and the other two tests aimed at evaluating the performance of thermally enhanced pipe installed in a fully saturated sedimentary rock formation. The concepts borrowed from pumping tests allowed the planning of the duration of the TRTs and the analysis of variable heat injection rate tests accounting for external heat transfer and temperature recovery, which reduces the uncertainty in the estimation of thermal properties.     

Special issue: Marine and river bedforms dynamics,MARID VI

Marine and river bedforms are rhythmic features driven by unidirectional or reversing currents and waves. They are ubiquitous on the floors of rivers, estuaries, coastal and marine settings. Despite a considerable history of dedicated studies, many aspects of their origin, development and dynamics are still the subject of scientific debate in various disciplines. The MARID conferences series hosts experts from geosciences, physics, engineering and other disciplines to provide a platform for the interdisciplinary exchange of fundamental and applied knowledge of marine and river dune dynamics. MARID VI took place in April 2019 in Bremen, Germany. Related scientific contributions are presented in this special issue of Earth Surface Processes and Landforms. An overview of the conference series and the specifics of MARID VI, as well as considerations in the context of equity, diversity and inclusion, are given. This introduction highlights the progress made with the papers published in the special issue.     

First assessment of the Caulerpa racemosa (Caulerpales, Chlorophyta) invasion along the French Mediterranean coast

The introduced green alga Caulerpa racemosa var. cylindracea has been rapidly spreading in the Mediterranean Sea since 1990. It was first observed in France in 1997 (Marseilles). In early 2004, the stretch of the French Mediterranean coastline and the surface area affected by the invasion were estimated at about 83 km and 4014 ha, respectively. The depth range of colonized areas was usually 10-35 m depth. Shallow (0-10 m) and deep (down to 40 m) dense meadows were rarely observed. In contrast to the dead matte of Posidonia oceanica, which constituted the most widely colonized substratum, dense P. oceanica meadows and fine sand with large ripple-marks were not invaded. Few rocky areas were colonized and coarse sand bottoms were usually colonized below 20 m depth. All the colonized areas were exposed to human activities and more than 40% were fishing areas. Mild climate, suitable substrata, presence of vectors of dispersal and absence of efficient biological control make the French Mediterranean coast particularly vulnerable to the further spread of the alga.     

3D geologic framework models for regional hydrogeology and land-use management: a case study from a Quaternary basin of southwestern Quebec, Canada

During a regional hydrogeologic survey in the St. Lawrence Lowlands, Canada, a computer-based 3D Geologic Framework Model (GFM) was constructed to obtain a consistent representation of this typical Quaternary glaciated basin over a 1,400 km² area. Such a detailed stratigraphic reconstruction was needed because the Quaternary sediments control the recharge to the underlying regional fractured rock aquifer and also because buried granular aquifers are partly connected to the regional system. The objectives of this geomodeling effort are 1) to improve understanding of subsurface conditions above the regional aquifer and; 2) to provide a common stratigraphic framework for hydrogeologic applications. The method draws on knowledge-driven discrete modeling using gOcad, as well as standardization and quality control procedures to maximize the use of a multisource database. The resulting model represents the bedrock topography and the complex stratigraphic architecture of overlying sediments. The regional till aquitard, the marine clay aquiclude and the buried granular aquifers have been modeled with unprecedented details thus providing a well-constrained 3D hydrostratigraphic framework. The recharge zones of the rock aquifer represent about 35% of the study area. Buried granular aquifers are directly connected to the regional aquifer system over about 10% of the area. The model allows several applications such as assessing aquifer vulnerability and areal groundwater recharge rates; improving the GFM inter-operability with groundwater modeling systems would be the next logical step.

---

Resumen Se construyó un Modelo del Marco Geológico (GFM) 3D basado en computadora durante un levantamiento hidrogeológico regional en las Tierras Bajas de St. Lawrence, Canada para obtener una representación consistente de esta cuenca glaciar Cuaternaria típica en un área de 1,400 km². Se necesitó tal grado de reconstrucción estratigráfica debido a que los sedimentos Cuaternarios controlan la recarga del acuífero rocoso fracturado regional subyacente y también porque los acuíferos granulares enterrados están parcialmente conectados con el sistema regional. Los objetivos de este esfuerzo de modelizado geológico son: 1) mejorar el entendimiento de las condiciones subsuperficiales por encima del acuífero regional y; 2) aportar un marco estratigráfico común para aplicaciones hidrogeológicas. El método se basa en el conocimiento de modelizado discreto utilizando gOcad, así como también en estandarización y procedimientos de control de calidad para maximizar el uso de bancos de datos de fuentes múltiples. El modelo obtenido representa la topografía del macizo rocoso y la arquitectura de la compleja estratigrafía de los sedimentos sobreyacentes. El acuitardo regional de tillita, el acuicludo arcilloso marino y los acuíferos granulares enterrados se han modelizado con detalles sin precedentes aportando de este modo un marco hidroestratigráfico 3D bien definido. Las zonas de recarga del acuífero rocoso representan casi el 35% del área de estudio. Los acuíferos granulares enterrados están conectados directamente al sistema acuífero regional sobre aproximadamente el 10% del área. El modelo permite varias aplicaciones tal como evaluar la vulnerabilidad de acuíferos y ritmos de recarga areales de agua subterránea; el paso lógico siguiente sería mejorar la inter-operabilidad del GFM con los sistemas de modelizado de agua subterránea.

---

Résumé Dans le cadre dune étude dhydrogéologie régionale dans les basses terres du Saint-Laurent (Canada), un Modèle surfacique géologique en 3D (MSG) a été développé pour obtenir une représentation cohérente par ordinateur de ce bassin quaternaire sur une superficie de plus de 1,400 km². Une telle reconstitution stratigraphique était nécessaire en raison du contrôle quexercent les sédiments quaternaires sur la recharge de laquifère régional fracturé sous-jacent et aussi parce que les aquifères granulaires enfouis sont partiellement connectés au système régional. Les objectifs de cet effort de géomodélisation sont: 1) daméliorer la compréhension des conditions de sous-surface au-dessus de laquifère régional et; 2) de fournir un cadre stratigraphique commun à des fins de caractérisation hydrogéologique. La méthode repose sur lintégration de linterprétation géologique dans la modélisation discrète effectuée à laide du logiciel gOcad, ainsi que sur des procédures duniformisation et de contrôle de la qualité des données afin doptimiser lutilisation dune base de données multisources. Le MSG représente la topographie du roc et larchitecture stratigraphique des sédiments sus-jacents. Laquitard régional (till), laquiclude dargile marine et les aquifères granulaires ont été modélisés à un niveau de détail sans précédent, fournissant ainsi un cadre hydrostratigraphique solidement établi. Les zones de recharge de laquifère rocheux représentent environ 35% de la zone détude. Les aquifères granulaires enfouis sont connectés directement au système aquifère régional sur au moins 10% du territoire. Le MSG permet de multiples usages tels lévaluation de la vulnérabilité des aquifères et lestimation de la distribution de la recharge; la prochaine étape logique serait lamélioration de linteropérabilité avec les systèmes de modélisation numérique de lécoulement.

     

Comparison of methods for the estimation of pyrite oxidation rate in a waste rock pile at Mine Doyon site, Quebec, Canada

Several methods were evaluated and compared for the estimation of pyrite oxidation rates (POR) in waste rock at Mine Doyon, Quebec, Canada. Methods based on data collected in situ, such as the interpretation of temperature and oxygen concentration profiles (TOP) measured in the waste rock pile and pyrite mass balance (PMB) on solid phase samples were compared with the oxygen consumption measurements (OCM) in closed chamber in the laboratory. A 1-D analytical solution to a gas and heat transport equation used temperature and oxygen profiles (TOP) measured in the pile for the preliminary POR estimates at a site close to the slope of the pile (Site 6) and in the core of the pile (Site 7). Resulting POR values were 1.1 × 10^− 9 mol(O₂) kg^− 1 s^− 1 and 1.0 × 10^− 10 mol(O₂) kg^− 1 s^− 1 for the slope site and the core site, respectively. Oxidation rates based on pyrite mass balance (PMB) calculations for solid samples were 2.21 × 10^− 9 mol(O₂) kg^− 1 s^− 1 and 2.03 × 10^− 9 mol(O₂) kg^− 1 s^− 1, respectively, for the same slope and core sites, but the difference between sites was within the error margin. The OCM measurements in the laboratory on fresh waste rock samples yielded higher POR values than field methods, with average oxidation rate of 6.7 × 10^− 8 mol(O₂) kg^− 1 s^− 1. However, the OCM results on weathered and decomposed material from the rock stockpile (average oxidation rate 3.4 × 10^− 9 mol(O₂) kg^− 1 s^− 1) were consistent with results from the field-based estimates. When POR values based on fresh material are excluded, the remaining POR values for all methods range from 1.0 × 10^− 10 to 3.4 × 10^− 9 mol(O₂) kg^− 1 s^− 1. The lowest estimated value (1.0 × 10^− 10 mol(O₂) kg^− 1 s^− 1) was based on TOP estimates in the interior of the pile where oxygen transport was limited by diffusion from the surface. These results suggest that small-scale OCM laboratory experiments may provide relatively representative values of POR in the zones of waste rock piles in which oxygen transport is not dominated by diffusion.     

Laboratory Study of Polymer Solutions Used for Mobility Control During In Situ NAPL Recovery

The use of surfactant solutions for the in situ recovery of residual NAPL in aquifers is increasingly considered as a viable remediation technique. The injection of a few pore volumes of high concentration surfactant solutions can mobilize most of the residual NAPL contacted by the solutions. However, the washing solutions'physico-chemical properties (low density and high viscosity), combined with the natural porous media heterogeneity, can prevent a good sweep of the entire contaminated volume. From the petroleum industry, it is well-known that polymer solutions can be injected following a surfactant solution slug to act as a mobility buffer and increase the overall sweep efficiency. The objective of our laboratory study is first to select and characterize polymers that would be suitable for aquifer restoration. Our experiments showed that among several polymers, xanthan gum solution rheology was made in order to predict shear rates, xanthan gum concentrations, salinity, and temperature effects on solution viscosity. The second set of experiments were made with a sand box which was designed to reproduce a simple heterogeneous media consisting of layers of sand with different permeability. These tests illustrate the xanthan gum solution's ability to increase surfactant solution's sweep efficiency and limit viscous fingering. The tests established that: (1) the injection of xanthan solution behind a surfactant solution slug decreases fluid velocity in high permeability layers and increases it in low-permeability ones, thus increasing the sweep efficiency (2) xanthan solutions eliminate viscous fingering at the polymer/surfactant solution front; (3) a xanthan solution preflush is desirable to limit surfactant solution mobility and prevent surfactant adsorption on solids; and (4) depending on site heterogeneity injection strategies should be applied to limit overriding by low-density surfactant solution.