Stochastic water balance model for rainfall recharge quantification in Ruataniwha Basin,New Zealand

Quantification of groundwater recharge is important for water resources management. Different methods can be used to estimate groundwater recharge. The most suitable approach depends on site characteristics. Water balance model was used in this study to quantify groundwater recharge from rainfall in Ruataniwha Basin, Hawke’s Bay, New Zealand. Because it is a closed basin, this method was determined to be more suitable than any other method. The basin is fed by surface water flow and rainfall, without any lateral groundwater flow. Records of surface water inflow and outflow are available with a certain degree of reliability. To cope with uncertainty in different components of water balance, different sampling methods were used. The Mean Value Latin Hypercube Sampling (MVLHS) was used and compared with Latin Hypercube Sampling (LHS), and Monte Carlo Simulation (MCS). It was found the groundwater recharge from rainfall is about 415 million m³/year. Results of this study have revealed that MVLHS converges faster and with lower variance than LHS and MCS.     

Elasto-plastic Behavior of Raghadan Tunnel Based on RMR and Hoek–Brown Classifications

Lining contact pressure and ground deformation of Raghadan transportation tunnel (Amman, Jordan) were investigated. The tunnel is 1.1 km in length and 13.5 m in diameter. This study was intended to integrate useful relations among the widely used rock classification system (RMR: rock mass rating), Hoek–Brown classification, and lining-ground interaction. The materials encountered along the tunnel alignment were limestone, dolomatic limestone, marly limestone, dolomite, and sillicified limestone. The ground conditions along the tunnel alignment including bedding planes, joint sets and joint conditions, rock quality, water flow, and rock strength were evaluated based on the drilled boreholes and rock exposures. Elasto-plastic finite element analyses were conducted to study the effect of rock mass conditions and tunnel face advance on the behavior of lining-ground interaction. The results of the analyses showed that lining contact pressure decreases linearly with the increase in RMR value. Also the results showed that tunnel lining contact pressure and crown inward displacement decreases with the increase in the unsupported distance (distance between tunnel face and the end of the erected lining). Ground displacement above the tunnel crown was found to be increases in an increasing rate with the decrease in the depth above the crown. This displacement was also found to be affected by the RMR value and the unsupported distance.     

Modelling surface–groundwater interaction in the Ruataniwha basin, Hawke’s Bay, New Zealand

Modelling groundwater and surface water is important for integrated water resources management, especially when interaction between the river and the aquifer is high. A transient groundwater and surface water flow model was built for Ruataniwha basin, New Zealand. The model covers a long-time period; starting in 1990, when water resources development in the area started, to present date. For a better resolution, the simulation period was divided into 59 stress periods, and each stress period was divided to 10 time steps. The model uses data obtained from surface water, and groundwater collected over the last 20 years. Rivers and streams were divided into 28 segments and flow and streambed data at the beginning and end of each segment was used. Parameter estimation and optimisation ‘PEST’ was used for automatic calibration of hydraulic conductivity, groundwater recharge and storativity; whereas riverbed conductance was manually calibrated. Model results show that the rivers gain from the aquifer considerably more than the river losses. The cumulative groundwater abstraction over the last 20 years is approximately 210 million m³. This amount is very low compared to other water budget components; however, the effect of groundwater abstraction on storage is significant. Based on the results of this study, it was found that the loss of storage over the last 20 years is more than 66 million m³. Results also reveal that the effect of groundwater abstraction on rivers and springs flow is significant. The rivers gain from the groundwater system, and the springs flow have been decreasing.     

LES-DEM simulations of sediment transport

In this work, a fully-coupled Computational Fluid Dynamics (CFD) model and Discrete Element Method (DEM) are used to simulate a unidirectional turbulent open-channel flow over the full range of sediment transport regimes. The fluid and particles are computed on separate grids using a dual-grid formulation to maintain consistency and avoid instability issues. The results of coupling the dispersed phase to a multiphase flow solver that uses volume-averaged Navier-Stokes equations are compared to those obtained from coupling through drag to a single flow solver. The current work also examines the applicability and limitations of lumping particles as a representative particle to reduce the cost of simulations. Insight to the impact of different turbulent events to the entrainment of particles is also given. The simulation results of sediment transport from both coupling techniques show good agreement with empirical formulas in the bedload regime, but under-predict sediment transport in the suspended load regime. In the suspended load regime, using partial coupling, the rate of sediment transport was found to be under-predicted as compared to full-coupling. The deviation in results in the suspended load regime was found to increase with increases in the applied shear stress. Both coupling methods revealed the same effect on the friction factor where friction increases in the bedload regime and decreases in the suspended load regime reaching a maximum at the transition between regimes. This result is contrary to past studies which have shown a discrete jump in the friction factor at the transition. Lumping particles as representative particles is shown to reduce the simulation cost by more than a factor of 5 when using a scaling factor of 2. By doing a quadrant analysis on information obtained from particle and flow field results, it was found that most of the particles are entrained by more frequent sweep events.     

Mapping groundwater contamination risk using GIS and groundwater modelling. A case study from the Gaza Strip, Palestine

Increasing pressure on water resources worldwide has resulted in groundwater contamination, and thus the deterioration of the groundwater resources and a threat to the public health. Risk mapping of groundwater contamination is an important tool for groundwater protection, land use management, and public health. This study presents a new approach for groundwater contamination risk mapping, based on hydrogeological setting, land use, contamination load, and groundwater modelling. The risk map is a product of probability of contamination and impact. This approach was applied on the Gaza Strip area in Palestine as a case study. A spatial analyst tool within Geographical Information System (GIS) was used to interpolate and manipulate data to develop GIS maps of vulnerability, land use, and contamination impact. A groundwater flow model for the area of study was also used to track the flow and to delineate the capture zones of public wells. The results show that areas of highest contamination risk occur in the southern cities of Khan Yunis and Rafah. The majority of public wells are located in an intermediate risk zone and four wells are in a high risk zone.     

Vulnerability assessment for the Gaza Strip, Palestine using DRASTIC

The main usefulness of groundwater vulnerability assessment maps is their ability to be an effective preliminary tool for planning, policy, and operational levels of decision-making. DRASTIC is one such assessment method. The DRASTIC index is made up of a calculated sum of products rating and weights for seven hydrogeological parameters that contribute to aquifer vulnerability. With the help of GIS, and based on the available data, maps of DRASTIC parameters were prepared for the Gaza Strip area in a case study. Each map was given a proper rate and a special weight factor developed. The final vulnerability map was obtained as a summation of the seven maps after multiplying each one with the appropriate weight. The vulnerability map was checked against the actual pollution potential in the area and nitrate concentration. The obtained vulnerability map is strongly correlated to known pollution values in the area.