Intrinsic groundwater vulnerability determination at the aquifer scale: a methodology coupling travel time estimation and rating methods

Groundwater vulnerability has been subject of much research due to the valuable information it provides concerning groundwater protection and exploitation potential. Up to now, most groundwater vulnerability studies adopt subjective systems of rating the various factors and subsequently, their results are often ambiguous and contradicting. Within the present study a methodology for the estimation of intrinsic groundwater vulnerability at the aquifer scale is presented. The methodology is based on travel time estimation from specified sources of pollution to the aquifer. Besides the deterministic calculation of travel times, the methodology provides a rating system for each pollution source, based on its relative severity and the estimated threat that it poses to the aquifer. Therefore, it can be regarded as a hybrid method that couples the advantages provided by the physically based methods with those of the subjective rating systems. The methodology is applied to the Neon Sidirochorion aquifer, Northeastern Greece, an overexploited aquifer where river waters, sea waters and lake waters interact, causing groundwater quality deterioration to the aquifer. The results indicated that the proposed groundwater vulnerability assessment methodology is well capturing pollution related to saltwater intrusion and agricultural activities, while it is concluded that the conceptual model is significantly affecting the vulnerability assessment results and therefore has to be previously developed.     

Estimation of air-overpressure produced by blasting operation through a neuro-genetic technique

Blasting operations usually produce significant environmental problems which may cause severe damage to the nearby areas. Air-overpressure (AOp) is one of the most important environmental impacts of blasting operations which needs to be predicted and subsequently controlled to minimize the potential risk of damage. In order to solve AOp problem in Hulu Langat granite quarry site, Malaysia, three non-linear methods namely empirical, artificial neural network (ANN) and a hybrid model of genetic algorithm (GA)–ANN were developed in this study. To do this, 76 blasting operations were investigated and relevant blasting parameters were measured in the site. The most influential parameters on AOp namely maximum charge per delay and the distance from the blast-face were considered as model inputs or predictors. Using the five randomly selected datasets and considering the modeling procedure of each method, 15 models were constructed for all predictive techniques. Several performance indices including coefficient of determination (R ²), root mean square error and variance account for were utilized to check the performance capacity of the predictive methods. Considering these performance indices and using simple ranking method, the best models for AOp prediction were selected. It was found that the GA–ANN technique can provide higher performance capacity in predicting AOp compared to other predictive methods. This is due to the fact that the GA–ANN model can optimize the weights and biases of the network connection for training by ANN. In this study, GA–ANN is introduced as superior model for solving AOp problem in Hulu Langat site.     

Textural evidences for dissolution of silica-rich rocks of the Ypresian phosphatic series,Gafsa-Metlaoui basin,southwestern Tunisia: implication of biogenic silica supply on genesis of fibrous clays

The present work aimed to determine the mineralogical composition of Ypresian series and to clarify the influence of the dissolution of siliceous frustules on the genesis of fibrous clay minerals. The specimens sampled from CPG trench are mainly constituted of silica-rich rocks at Mides area located at the western part of Gafsa-Metlaoui basin. The samples were investigated using X-ray diffraction (XRD) and scanning electron microscopy (SEM) in order to determine texture of constituents. The data obtained indicate that the bulk rock samples are mainly made up of opal CT and clay minerals. The latter consists of palygorskite-sepiolite minerals associated with smectite and few amount of illite. The trend of mineralogical composition indicates that fibrous clay minerals are more concentrated at the upper part. SEM observations indicated that palygorskite mineral appears as thread-like facies, which surround foliated texture of smectite in the lower part of the Mides section, although with the low Mg activity confirmed by the absence of dolomite. But, at the upper part of the Mides section, SEM observations revealed the occurrence of siliceous frustules, which have numerous dissolved areas and replacement of carbonate tests by silica. The dissolution saturated the depositional environment with silica which is required for the formation of palygorskite and sepiolite minerals, in addition to high Mg activity confirmed by the presence of dolomite in the bulk rock, which is required basically for the formation of sepiolite. Although the genesis mode of palygorskite and sepiolite is similar with very little difference, the genesis of sepiolite needs a high alkalinity than the formation of palygorskite.     

Interior boundary-aligned unstructured grid generation and cell-centered versus vertex-centered CVD-MPFA performance

Grid generation for reservoir simulation must honor classical key constraints and be boundary aligned such that control-volume boundaries are aligned with geological features such as layers, shale barriers, fractures, faults, pinch-outs, and multilateral wells. An unstructured grid generation procedure is proposed that automates control-volume and/or control point boundary alignment and yields a PEBI-mesh both with respect to primal and dual (essentially PEBI) cells. In order to honor geological features in the primal configuration, we introduce the idea of protection circles, and to generate a dual-cell feature based grid, we construct halos around key geological features. The grids generated are employed to study comparative performance of cell-centred versus cell-vertex control-volume distributed multi-point flux approximation (CVD-MPFA) finite-volume formulations using equivalent degrees of freedom. The formulation of CVD-MPFA schemes in cell-centred and cell-vertex modes is analogous and requires switching control volume from primal to dual or vice versa together with appropriate data structures and boundary conditions. The relative benefits of both types of approximation, i.e., cell-centred versus vertex-centred, are made clear in terms of flow resolution and degrees of freedom required.     

A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

Seismic stratigraphy of the Mianwali and Bannu depressions,north-western Indus foreland basin

International Journal of Earth Sciences - Regional seismic reflection profiles, deep exploratory wells, and outcrop data have been used to study the structure and stratigraphic architecture of the...