Abstract: It is well known that rainwater harvesting (RWH) can augment water supply and reduce stormwater pollutant discharges. Due to the lack of continuous 2D modelling of urban flood coverage and its associated damage, the ability of RWH to reduce urban flood risks has not been fully evaluated. Literature suggests that small distributed storage spaces using RWH tanks will reduce flood damage only during small to medium flooding events and therefore cumulative assessment of their benefits is needed. In this study we developed a new integrated modelling framework that implements a semi-continuous simulation approach to investigate flood prevention and water supply benefits of RWH tanks. The framework includes a continuous mass balance simulation model that considers antecedent rainfall conditions and water demand/usage of tanks and predicts the available storage prior to each storm event. To do so, this model couples a rainfall-runoff tank storage model with a detailed stochastic end-use water demand model. The available storage capacity of tanks is then used as a boundary condition for the novel rapid flood simulation model. This flood model was developed by coupling the U.S. EPA Storm Water Management Model (SWMM) to the Cellular-Automata Fast Flood Evaluation (CA-ffé) model to predict the inundation depth caused by surcharges over the capacity of the drainage network. The stage-depth damage curves method was used to calculate time series of flood damage, which are then directly used for flood risk and cost-benefit analysis. The model was tested through a case study in Melbourne, using a recorded rainfall time series of 85 years (after validating the flood model against 1D-2D MIKE-FLOOD). Results showed that extensive implementation of RWH tanks in the study area is economically feasible and can reduce expected annual damage in the catchment by up to approximately 30 percent. Availability of storage space and temporal distribution of rainfall within an event were important factors affecting tank performance for flood reduction.

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: 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.