Abstract: This study combines historical records of 14C and 3H in the atmosphere and soil with renewal rate and groundwater lumped parameter models to predict the abundance of 14C and 3H in groundwater over time. 624 groundwater samples from numerous studies, over four decades (1978–2019), in South Africa, Namibia, Botswana and Mozambique were collated to compare with predicted groundwater activities of 14C and 3H within the South African Development Community (SADC) region. Spatial datasets of carbonate bearing lithology, C3/C4 vegetation, summer/winter rainfall and coastal proximity were used to apply corrections to 14C and 3H data. Corrected values of 14C and 3H were compared with the theoretical abundance of these tracers, derived from the lumped parameter models, to estimate the general mean residence times and presence of groundwater mixing between modern recharge and older groundwaters. This study found that corrected values produced varying mean residence times derived from 14C ages (∼500–28500 years) and a wide range of potentially mixed waters within each aquifer system (0–100 % of tested wells) across the study area. The largest proportions of mixed groundwater, as well as the youngest mean residence times, were found in alluvial and primary fractured rock aquifers (e.g., western coast of South Africa and southern Mozambique). The smallest proportions of mixed groundwater were predicted in deep confined clay-rich aquifers as well as layered coal bearing carbonate sequences (e.g., Orapa, Malwewe and Serowe, Botswana). Insights into the proportions of mixed groundwater and mean residence times can help assess hydrological resilience on a regional scale. Such information is pertinent in promoting socio-economic development and increased water/food security in the SADC region. By understanding the resilience of groundwater resources, robust and informed strategies for water equality and GDP growth in the SADC region can be envisioned and implemented.

Abstract: Qanat is a passive (unpumped) horizontal well (or a slant well with a very mild inclined angle) that is capable of extracting water from aquifers by gravity. Many qanats are constructed along the radius of the alluvial fan wedge-shaped aquifers. Analytical modeling of such a qanat-aquifer system provides great benefit for quickly screening different designs of qanats and accessing the performance of qanat discharge in the field. The previous analytical modeling of discharge of qanats, however, did not consider the wedge-shaped aquifers. Thus, the goal of this research is to obtain semi-analytical solutions of discharge variations of qanats in alluvial fan aquifers with nearby pumping wells, subjected to areal recharges due to rainfall. The uniform head boundary is considered inside the qanat (because of its enormous permeability in respect to the background aquifer). The influences of the aquifer lateral boundaries on discharge of qanat and its sensitivity to hydraulic and geometric parameters are explored. The influences of the lateral boundaries on the discharge of qanat due to areal recharge and nearby pumping wells discharge are also explored. The results of this study can be utilized for multiple purposes: 1) to predict the discharge of qanat in an alluvial fan aquifer and explore the influences of the areal recharge and nearby pumping well discharge; 2) to estimate the hydraulic parameters of the alluvial fan aquifer depleted by a qanat; 3) to determine the location of the nearby pumping well to minimize its influences on the discharge of a qanat; 4) to calculate the water budgets of aquifers depleted by qanats and pumping wells and replenished by areal recharge among other applications. This paper is an extension to the work presented by Sedghi and Zhan (2020) (which concerns an infinite unconfined aquifer) for an unconfined alluvial fan aquifer setting.

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: 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: Concentrations of uranium (U) \textgreater30 µg/L in groundwater are relatively uncommon in drinking water in the United States but can be of concern in those areas where complex interactions of aquifer materials and anthropogenic alterations of the natural flow regime mobilize U. High concentrations (\textgreater30 µg/L) of U in the southeastern San Joaquin Valley, California, USA, have been detected in 24 percent of 257 domestic, irrigation, and public-supply wells sampled across an approximately 110,000 km2 area. In this study we evaluated mechanisms for mobilization of U in the San Joaquin Valley proposed in previous studies, confirming mobilization by HCO3 and refuting mobilization by NO3 and we refined our understanding of the geologic sources of U to the scale of individual alluvial fans. The location of high concentrations depends on the interactions of geological U sources from fluvial fans that originate in the Sierra Nevada to the east and seepage of irrigation water that contains high concentrations of HCO3 that leaches U from the sediments. In addition, interactions with PO4 from fertilized irrigated fields may sequester U in the aquifer. Principal component analysis of the data demonstrates that HCO3 and ions associated with high total dissolved solids in the aquifer and the percentage of agriculture near the well sampled are associated with high U concentrations. Nitrate concentrations do not appear to control release of U to the aquifer. Age dating of the groundwater and generally increasing U concentrations of the past 25 years in resampled wells where irrigation is prevalent suggests that high U concentrations are associated with younger water, indicating that irrigation of fields over the past 100 years has significantly contributed to increasing concentrations and mobilizing U. In some places, the groundwater is supersaturated with uranyl-containing minerals, as would be expected in roll front deposits. In general, the interaction of natural geological sources high in U, the anthropogenically driven addition of HCO3 and possibly phosphate fertilizer, control the location and concentration of U in each individual fluvial fan, but the addition of nitrate in fertilizer does not appear control the location of high U. These geochemical interactions are complex but can be used to determine controls on anomalously high U in alluvial aquifers.