Abstract: The quantification of groundwater resources is essential especially in water scarce countries like Namibia. The chloride mass balance (CMB) method and isotopic composition were used in determining groundwater recharge along a precipitation gradient at three sites, namely: Tsumeb (600 mm/a precipitation); Waterberg (450 mm/a precipitation) and Kuzikus/Ebenhaezer (240 mm/a precipitation). Groundwater and rainwater were collected from year 2016–2017. Rainwater was collected monthly while groundwater was collected before, during and after rainy seasons. Rainwater isotopic values for δ18O and δ2H range from −10.70 to 6.10‰ and from −72.7 to 42.1‰ respectively. Groundwater isotopic values for δ18O range from −9.84 to −5.35‰ for Tsumeb; from −10.85 to −8.60‰ for Waterberg and from −8.24 to −1.56‰ for Kuzikus/Ebenhaezer, while that for δ2H range from −65.6 to −46.7‰ for Tsumeb; −69.4 to −61.2‰ for Waterberg and −54.2 to −22.7‰ for Kuzikus/Ebenhaezer. Rainwater scatters along the GMWL. Rainwater collected in January, February and March are more depleted in heavy isotopes than those in November, December, April and May. Waterberg groundwater plots on the GMWL which indicates absence of evaporation. Tsumeb groundwater plots on/close to the GMWL with an exception of groundwater from the karst Lake Otjikoto which is showing evaporation. Groundwater from Kuzikus/Ebenhaezer shows an evaporation effect, probably evaporation occurs during infiltration since it is observed in all sampling seasons. All groundwater from three sites plot in the same area with rainwater depleted in stable isotopic values, which could indicates that recharge only take place during January, February and March. CMB method revealed that Waterberg has the highest recharge rate ranging between 39.1 mm/a and 51.1 mm/a (8.7% – 11.4% of annual precipitation), Tsumeb with rates ranging from 21.1 mm/a to 48.5 mm/a (3.5% – 8.1% of annual precipitation), and lastly Kuzikus/Ebenhaezer from 3.2 mm/a to 17.5 mm/a (1.4% – 7.3% of annual precipitation). High recharge rates in Waterberg could be related to fast infiltration and absence of evaporation as indicated by the isotopic ratios. Differences in recharge rates cannot only be attributed to the precipitation gradient but also to the evaporation rates and the presence of preferential flow paths. Recharge rates estimated for these three sites can be used in managing the savannah aquifers especially at Kuzikus/Ebenhaezer where evaporation effect is observed that one can consider rain harvesting.

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: Qanat is a type of drain that extract water from aquifers by gravity. Significant amount of fresh water used in Middle East and other parts of the world are supplied by qanats. Despite their importance, discharge variation of these type of wells received almost no attention. The aim of this research is to obtain a Laplace domain solution of discharge variation of a qanat installed in an anisotropic unconfined aquifer subjected to arbitrary areal recharge and nearby pumping well(s) discharge. A new semi-analytical solution of drawdown is obtained first to implement the effects of arbitrary areal recharge and nearby pumping well(s) using the principle of superposition. Then, the discharge variation solution of the qanat is obtained from the drawdown solution. To establish a constant-head boundary condition at the qanat periphery, the qanat is discretized into several segments. The results of this study are presented in dimensionless discharge-dimensionless time curves. The effects of hydraulic as well as geometric parameters on the discharge variation of the qanat due to arbitrary areal recharge, falling of water table from its initial position and discharge of nearby wells are explored. We also investigate the influences of distance and screen depth and location of the nearby well on the discharge variation of the qanat. The results of this study can be utilized for multiple purposes: 1) to predict discharge of qanat in response to rainfall and nearby pumping well(s); 2) to estimate the aquifer parameters using hydrograph of the qanat; 3) to determine optimal location and pumping pattern of the nearby wells to minimize their influences on the discharge of the qanat; 4) to calculate water budget of aquifers drained by a qanat. The equation presented in this work can also be used to estimate discharge of a horizontal drain installed in cropland subjected to arbitrary irrigation pattern.

Abstract: In Situ Recovery (ISR) is the most important process of uranium production in the world (50%). It consists of an injection of a leaching solution into a permeable mineralized aquifer (sandstone), pumping of the solution after dissolution of the ore minerals and recovery of the uranium from the pregnant solution in a surface plant. In this context, the fate of swelling clay minerals such as smectites is of main importance due to their role in the mobility of cationic elements by diverse geochemical processes such as ion-exchange reactions or dissolution. The present study details analysis of dioctahedral smectites before and after in-situ leaching by acidic (H2SO4) ISR solutions. Samples were collected from two sedimentary basins hosting some of the main uranium roll front deposits exploited by acidic ISR (Tortkuduk deposit, Shu-Saryssu basin, Kazakhstan, and Dulaan Uul and Zoovch Ovoo deposits, Sainshand basin, Mongolia). Scanning Electron Microscope and X-Ray Diffraction analysis revealed that dioctahedral smectite is a ubiquitous mineral in all analyzed samples, before and after acidification, and revealed a difference of crystal chemistry of the smectites between deposits of Kazakhstan (beidellite type) and Mongolia (montmorillonite type). Chemical analysis and semi-quantification of the smectites before and after acidification also revealed a difference in chemical reactivity, with a higher dissolution of montmorillonite layers compared to beidellite ones, and the importance of ion-exchange reactions. These findings are consistent with literature data obtained on model systems. The persistence of dioctahedral smectites after several years of acidification is crucial for the understanding of geochemical processes during uranium production or remediation of the aquifers. Finally, based on the analysis of samples from U-deposits hosted in both sedimentary basins, a schematic model of the impact of acid solutions on dioctahedral smectite was proposed.

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.