Delineation of groundwater potential zones using remote sensing and GIS-based data-driven models

The rapid increase in human population has increased the groundwater resources demand for drinking, agricultural and industrial purposes. The main purpose of this study is to produce groundwater potential map (GPM) using weights-of-evidence (WOE) and evidential belief function (EBF) models based on geographic information system in the Azna Plain, Lorestan Province, Iran. A total number of 370 groundwater wells with discharge more than 10 m³s^?1were considered and out of them, 256 (70%) were randomly selected for training purpose, while the remaining114 (30%) were used for validating the model. In next step, the effective factors on the groundwater potential such as altitude, slope aspect, slope angle, curvature, distance from rivers, drainage density, topographic wetness index, fault distance, fault density, lithology and land use were derived from the spatial geodatabases. Subsequently, the GPM was produced using WOE and EBF models. Finally, the validation of the GPMs was carried out using areas under the ROC curve (AUC). Results showed that the GPM prepared using WOE model has the success rate of 73.62%. Similarly, the AUC plot showed 76.21% prediction accuracy for the EBF model which means both the models performed fairly good predication accuracy. The GPMs are useful sources for planners and engineers in water resource management, land use planning and hazard mitigation purpose.     

Hydrogeological studies on the Nubian sandstone aquifer in El-Bahariya Oasis, Western Desert, Egypt

The current study aimed to evaluate hydrogeologically the Nubian sandstone aquifer in El-Bahariya Oasis. It represents the main water-bearing horizon in the study area and consists of continental elastic sediments, mainly sandstone alternating with shale and clays. The general flow lines are directed from SW to NE direction, as detected from the constructed potentiometric head contour map. The piezometric surface reaches 149 m in El-Heiz area at the southern part, while it reaches 90 m at the northern, reflecting higher pressure head of the aquifer in the southern part. The map also illustrates that the southern part is considered as the most promising location for development. The structural elements play an important role in the deposition and distribution of the sedimentary succession of the Nubian sandstone sediments. Consequently, this sedimentary pattern affects the occurrences and movements of the groundwater within the aquifer system. Along the structurally high areas, in the study area, the piezometric head increases, while the reverse is recorded along the structurally low areas. The step-drawdown tests data were carried out by calculating the aquifer loss coefficient (B) and the well loss constant (C). The B values are smaller compared with C values, indicating that the aquifer under pressure has a behavior of leaky aquifer; therefore, it shows hydraulic connection with surrounding formation. The values of well efficiency range from 78.50% to 87.76%. Analysis of 12 pumping test data (constant discharge tests) was carried out in order to calculate the Nubian aquifer hydraulic parameters (transmissivity, hydraulic conductivity, and storage coefficient). The transmissivity values decrease from 3,045 m²/day in the southern part (El-Heiz area) to 236 m²/day in the northeastern part (El-Harra area). Accordingly, the aquifer classified as a high to moderate potentiality. Transmissivity contour map observes gradual increase of transmissivity values from the southern to northeastern direction. This may be due to the increase of shale or clay content in the concerned aquifer in that direction. The storage coefficient values range between 1.04 × 10^?4 and 5.22 × 10^?3, as obtained from the results of pumping test analysis, which ensure that the Nubian sandstone aquifer is classified as semi-confined to confined aquifer type. The S values show a decrease from southwest to northeast direction as detected from S-map. The hydraulic conductivity values vary from to 0.46 m/day in the northern part to 10.88 m/day in the southern part with an average of 5.67 m/day. According to the classification based on K values, the aquifer is mainly composed of coarse sand.     

Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigation of the groundwater modelling component of the Integrated Water Flow Model (IWFM)

ABSTRACT

The Integrated Water Flow Model (IWFM), developed by the California Department of Water Resources, is an integrated hydrological model that simulates key flow processes including groundwater flows, streamflow, stream–aquifer interactions, rainfall–runoff and infiltration. It also simulates the agricultural water demand as a function of soil, crop and climatic characteristics, as well as irrigation practices, and allows the user to meet these demands through pumping and stream diversions. This study investigates the modelling performance of the groundwater module of IWFM using several hypothetical test problems that cover a wide range of settings and boundary conditions, by comparing the simulation results with analytical solutions, field and laboratory observations, or with results from MODFLOW outputs. The comparisons demonstrate that IWFM is capable of simulating various hydrological processes reliably.

EDITOR M.C. Acreman; ASSOCIATE EDITOR A. Efstratiadis     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Optimization of a Photo‐assisted Fenton Oxidation Process: A Statistical Model for MDF Effluent Treatment

In the present study, the effects of initial COD (chemical oxygen demand), initial pH, Fe²⁺/H₂O₂ molar ratio and UV contact time on COD removal from medium density fiberboard (MDF) wastewater using photo‐assisted Fenton oxidation treatment were investigated. In order to optimize the removal efficiency, batch operations were carried out. The influence of the aforementioned parameters on COD removal efficiency was studied using response surface methodology (RSM). The optimal conditions for maximum COD removal efficiency from MDF wastewater under experimental conditions were obtained at initial COD of 4000 mg/L, Fe²⁺/H₂O₂ molar ratio of 0.11, initial solution pH of 6.5 and UV contact time of 70 min. The obtained results for maximum COD removal efficiency of 96% revealed that photo‐assisted Fenton oxidation is very effective for treating MDF wastewater.     

Evaluation of the physico-mechanical parameters affecting the deterioration rate of Ahlat ignimbrites (Bitlis,Turkey)

The paper principally focuses on the durability assessment of various stratigraphic levels of Ahlat ignimbrites collected from the eastern region of Turkey. A total of four different ignimbrite types with dissimilar color, texture and particularly welding degree were tested in laboratory. The laboratory tests performed on the ignimbrite specimens indicate that the welding degree as well as the lithic material content mainly controls the strength and capillarity properties of the ignimbrites. In addition, the durability of highly porous ignimbrites strongly depends upon the degree of welding. The effect of several weathering agents on the ignimbrites was evaluated on the basis of decay constant parameter. Accordingly, salt and ice crystallization pressures are a couple of major destructive agents acting within the micropores of the ignimbrites. Conversely, the investigated specimens are relatively durable against cyclic wetting–drying. Statistical evaluations reveal that the pore diameter is the major controlling factor on the deterioration rate of the ignimbrites after specifically recurrent freeze–thaw cycles. Moreover, the dry unit weight of the ignimbrites is more significant than the uniaxial compressive strength considering the deterioration rates during wetting–drying and salt crystallization. A less significant relationship was obtained between pore diameter and salt crystallization decay constant.     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.     

Link between SSZ ophiolite formation, emplacement and arc inception, Northland, New Zealand: U–Pb SHRIMP constraints; Cenozoic SW Pacific tectonic implications

New U–Pb age-data from zircons separated from a Northland ophiolite gabbro yield a mean ²⁰⁶Pb/²³⁸U age of 31.6 ± 0.2 Ma, providing support for a recently determined 28.3 ± 0.2 Ma SHRIMP age of an associated plagiogranite and  29–26 Ma ⁴⁰Ar/³⁹Ar ages (n = 9) of basalts of the ophiolite. Elsewhere, Miocene arc-related calc-alkaline andesite dikes which intrude the ophiolitic rocks contain zircons which yield mean ²⁰⁶Pb/²³⁸U ages of 20.1 ± 0.2 and 19.8 ± 0.2 Ma. The ophiolite gabbro and the andesites both contain rare inherited zircons ranging from 122–104 Ma. The Early Cretaceous zircons in the arc andesites are interpreted as xenocrysts from the Mt. Camel basement terrane through which magmas of the Northland Miocene arc lavas erupted. The inherited zircons in the ophiolite gabbros suggest that a small fraction of this basement was introduced into the suboceanic mantle by subduction and mixed with mantle melts during ophiolite formation.

We postulate that the tholeiitic suite of the ophiolite represents the crustal segment of SSZ lithosphere (SSZL) generated in the southern South Fiji Basin (SFB) at a northeast-dipping subduction zone that was initiated at about 35 Ma. The subduction zone nucleated along a pre-existing transform boundary separating circa 45–20 Ma oceanic lithosphere to the north and west of the Northland Peninsula from nascent back arc basin lithosphere of the SFB. Construction of the SSZL propagated southward along the transform boundary as the SFB continued to unzip to the southeast. After subduction of a large portion of oceanic lithosphere by about 26 Ma and collision of the SSZL with New Zealand, compression between the Australian Plate and the Pacific Plate was taken up along a new southwest-dipping subduction zone behind the SSZL. Renewed volcanism began in the oceanic forearc at 25 Ma producing boninitic-like, SSZ and within-plate alkalic and calc-alkaline rocks. Rocks of these types temporally overlap ophiolite emplacement and subsequent Miocene continental arc construction.