Abstract: The origin of groundwater recharge and subsequent flow paths are often difficult to establish in fractured, multi-lithological, and highly compartmentalized aquifers such as the Troodos Fractured Aquifer (TFA). As the conjunctive use of stable isotopes and hydrogeochemical data provides additional information, we established a monitoring network for stable isotopes in precipitation in Cyprus. The local meteoric water line, altitude effect and seasonal variation of stable isotopes in precipitation are derived from monitoring data. Stable isotopes and hydrogeochemical data are combined to model water-rock interactions and groundwater evolution along a complete ophiolite sequence. As a result a generic hydrogeologic description for the observed water types is developed. Isotope hydrology was applied in conjunction with hydrogeochemical modelling in Kargiotis Watershed, a major north-south transect of the TFA. PHREEQC was used for hydrogeochemical modelling to establish generic descriptions for observed water types. Mean precipitation-weighted values from 16 monitoring stations were used to calculate the Local Meteoric Water Line (LMWL), which was found to be equal to δ2H = (6.58 ± 0.13)*δ18O + (12.64 ± 0.91). A general decrease of 1.22‰ for δ2H and 0.20‰ for δ18O in precipitation was calculated per 100 m altitude. A generic groundwater evolution path was established: 1. Na/MgClHCO3, 2. MgHCO3, 3. Ca/MgHCO3, 4. Ca/MgNaHCO3, 4a. MgNa/CaHCO3/Cl, 5. NaMg/CaHCO3/Cl, 6. NaHCO3, 7. Na/MgHCO3SO4, 8. NaSO4Cl/HCO3. Hydrogeologic descriptions, consisting of groundwater origin, flow path and possible active water-rock processes, have been realised for the observed water types. The first two water types occur in serpentine and ultramafic-gabbro springs. Type 3 waters represent early stages of recharge and/or short flow paths, in gabbro whereas types 4 and 5 are typical for further percolating waters in gabbro and diabase. Water types 6 and 7 occur both in diabase and in the basal group and represent the regional flow. Water type 8 is the end member of regional, upwelling groundwater in the basal group. The presented descriptions and methods have practical applications in groundwater exploration, characterization, and protection. The methodology can be applied in other complex aquifer systems.

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.