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.