Abstract: In this paper, a new numerical model has been developed for simulation of Qanat-aquifer system. This model employs quasi-3D mixed-form of Richards’ equation and 1D fully-hydrodynamic form of Saint-Venant equations to simulate subsurface and overland flow, respectively. In order to handle non-orthogonal grids, subsurface flow module benefits from coordinate transformation technique. Using the above-mentioned governing equations, the presented model is able to simulate water flow inside both collection and conveyance sections of the gallery as well as dynamics of groundwater and vadose zone from impermeable bed rock to the soil-air interface. Since measured data corresponding to the hydraulics of Qanats is scarce, the overland and subsurface modules have been validated with analytical, numerical and experimental benchmarks in the literature. Subsequently, the model was employed to simulate ten different hypothetical aquifer-Qanat systems with different properties including the depth of groundwater aquifer, roughness of the gallery and saturated hydraulic conductivity of the gallery-aquifer boundary and the influence of each the parameters was monitored on the outflow rate at the appearance point of each Qanat. Furthermore, the advance of water inside two initially dry galleries were simulated at different time levels up to steady state. Eventually, the streamlines have been shown at the steady state for two Qanat-aquifer systems. Although, the presented study sheds light on some aspects of Qanat-aquifer hydraulics, the validation of the presented model with in-lab or on-field data remains ongoing for the future researches.

Abstract: The development of adaptation and mitigation strategies to tackle anthropic and climate changes impacts is becoming a priority in drought-prone areas. This study examines the capabilities of indigenous rainwater harvesting techniques (RWHT) to be used as a viable solution for flood mitigation. The study analyses the hydraulic performance of the most used micro-catchment RWHT in sub-Saharan regions, in terms of flow peak reduction (FPR) and volume reduction (VR) at the field and basin scale. Parametrized hyetographs were built to replicate the extreme precipitations that strike Sahelian countries during rainy seasons. 2D hydrodynamic simulations showed that half-moons placed with a staggered configuration (S-HM) have the best performances in reducing runoff. At the field scale, S-HM showed a remarkable FPR of 77% and a VR of 70% in case of extreme rainfall. Instead at the basin scale, in which only 5% of the surface was treated, 13% and 8% respectively for FPR and VR were obtained. In addition, the reduction of the runoff coefficient (Rc) between the different configuration was analyzed. The study critically evaluates hydraulic performances of the different techniques and shows how pitting practices cannot guarantee high performance in case of extreme precipitations. These results will enrich the knowledge of the hydraulic behavior of RWHT; aspect marginally investigated in the scientific literature. Moreover, this study presents the first scientific application of HEC-RAS as a rainfall-runoff model. Despite some limitations, this model has the effective feature of using very high-resolution topography as input for hydraulic simulations. The results presented in this study should encourage stakeholders to upscale the use of RWHT in order to lessen the flood hazard and land degradation that oppresses arid and semi-arid areas.

Abstract: The hydrogeological functioning of four different areas in a complex evaporite-karst unit of predominantly aquitard behavior in S Spain was investigated. Environmental dating tracers (3H, 3He, 4He, CFC-12, SF6) and hydrochemical data were determined from spring samples to identify and characterize groundwater flow components of different residence times in the media. Results show a general geochemical evolution pattern, from higher (recharge areas) to lower positions (discharge areas), in which mineralization rises as well as the value of the rCl−/SO42−, evidencing longer water-rock interaction. Ne values show degassing of most of the samples, favored by the high salinity of groundwater and the development of karstification so that the concentration of all the considered gases were corrected according to the difference between the theoretical and the measured Ne. The presence of modern groundwater in every sample was proved by the detection of 3H and CFC-12. At the opposite, the higher amount of radiogenic 4He in most samples also indicates that they have an old component. The 3H/3He dating method does not give reliable ages as a consequence of degassing and the large uncertainty of the 3He/4He ratios of the sources for the radiogenic Helium. The large SF6 concentrations suggest terrigenic production related to halite and dolomite. Binary Mixing and Free Shape Models were created based on 3H and CFC-12 data to interpret the age distribution of the samples. Two parameters (GA50 and >70%) were proposed as an indicator of that distribution, as they provide further information than the mean age. Particularly, GA50 is derived from the median groundwater age and is presented as a new way of interpreting mixed groundwater age data. A greater fraction of old groundwater (3H and CFC-12 free) was identified in discharge areas, while the proportion and estimated infiltration date of the younger fractions in recharge areas were higher and more recent, respectively. The application of different approaches has been useful to corroborate previous theoretical conceptual model proposed for the study area and to test the applicability of the used environmental tracer in dating brine groundwater and karst springs.

Abstract: Complex pyrite textures associated with large changes in isotopic and trace element compositions are routinely assumed to be indicative of multi-faceted processes involving multiple fluid and sulfur sources. We propose that the features of ore-stage pyrite from roll front deposits across the world, revealed in exquisite detail via high-resolution trace element mapping by LA-ICP-MS, reflect the dynamic internal evolution of the biogeochemical processes responsible for sulfate reduction, rather than externally driven changes in fluid or sulfur sources through time. Upon percolation of oxidizing fluids into the reduced host-sandstones, roll front systems become self-organized, with a systematic reset of their activity cycle after each translation stage of the redox interface down dip of the aquifer. Dominantly reducing conditions at the redox interface favor the formation of biogenic framboidal pyrite (δ34S from −30.5 to −12.5‰) by bacterial sulfate reduction and the genesis of the U mineralization. As the oxidation front advances, oxidation of reduced sulfur minerals induces an increased supply of sulfate and metals in solution to the bacterial sulfate reduction zone that has similarly advanced down the flow gradient. Hence, this stage is marked by increased rates of the bacterial sulfate reduction associated with the crystallization of variably As-Co-Ni-Mo-enriched concentric pyrite (up to 10,000′s of ppm total trace contents) with moderately negative δ34S values (from −13.7 to −7.5‰). A final stage of pyrite cement with low trace element contents and heavier δ34S signature (from −6.9 to +18.8‰) marks the end of the roll front activity cycle and the transition from an open to a predominantly closed system behavior (negligible advection of fresh sulfate). Blocky pyrite cement is formed using the remaining sulfate, which now becomes quickly heavy according to a Rayleigh isotope fractionation process. This ends the cycle by depleting the nutrient supplies for the sulfate-reducing bacteria and cementing pore spaces within the host sandstone, effectively restricting fluid infiltration. This internally-driven roll front activity cycle results in systematic, large S isotope and trace element fractionation. Ultimately, the long-time evolution of the basin and fluid sources control the metal endowment and evolution of the system; these events, however, are unlikely to be preserved by the roll front, as a direct result of its hydrodynamic nature.

Abstract: A uranium-mineralized sandy aquifer, planned for mining by means of uranium in situ recovery (U ISR), harbors a reservoir of bacterial life that may influence the biogeochemical cycles surrounding uranium roll-front deposits. Since microorganisms play an important role at all stages of U ISR, a better knowledge of the resident bacteria before any ISR actuations is essential to face environmental quality assessment. The focus here was on the characterization of bacteria residing in an aquifer surrounding a uranium roll-front deposit that forms part of an ISR facility project at Zoovch Ovoo (Mongolia). Water samples were collected following the natural redox zonation inherited in the native aquifer, including the mineralized orebody, as well as compartments located both upstream (oxidized waters) and downstream (reduced waters) of this area. An imposed chemical zonation for all sensitive redox elements through the roll-front system was observed. In addition, high-throughput sequencing data showed that the bacterial community structure was shaped by the redox gradient and oxygen availability. Several interesting bacteria were identified, including sulphate-reducing (e.g. Desulfovibrio, Nitrospira), iron-reducing (e.g. Gallionella, Sideroxydans), iron-oxidizing (e.g. Rhodobacter, Albidiferax, Ferribacterium), and nitrate-reducing bacteria (e.g. Pseudomonas, Aquabacterium), which may also be involved in metal reduction (e.g. Desulfovibrio, Ferribacterium, Pseudomonas, Albidiferax, Caulobacter, Zooglea). Canonical correspondence analysis (CCA) and co-occurrence patterns confirmed strong correlations among the bacterial genera, suggesting either shared/preferred environmental conditions or the performance of similar/complementary functions. As a whole, the bacterial community residing in each aquifer compartment would appear to define an ecologically functional ecosystem, containing suitable microorganisms (e.g. acidophilic bacteria) prone to promote the remediation of the acidified aquifer by natural attenuation. Assessing the composition and structure of the aquifer’s native bacteria is a prerequisite for understanding natural attenuation and predicting the role of bacterial input in improving ISR efficiency.