Hydrochemical and bacteriological features of the groundwater: southern border of the Rharb basin (Morocco) / Caractérisation hydrochimique et bactériologique des eaux souterraines: bordure méridionale du bassin du Rharb (Maroc)

Abstract

The western reservoirs represent the principal groundwater system in Morocco. Demographic, industrial and agricultural developments during the last decade have markedly altered groundwater quality. The Mamora coastal aquifer system is among the Atlantic systems which are most heavily threatened by pollution. Agricultural and industrial activities, and rapid urban growth contribute to the pollution of the groundwater. Contamination transport is facilitated by a high permeability of the aquifer formations. In order to assess the actual groundwater quality of the Mamora aquifer and to understand the influence of the factors generating the pollution, an extensive multidisciplinary research programme is in progress, with hydrochemistry and microbiology playing essential roles. The present paper concerns the spatial distribution of physico-chemical parameters in the groundwater, subjected to domestic, industrial and agricultural pollution. Fifty-seven samples were analysed for several parameters (Ca²⁺, Mg²⁺, Na⁺, K⁺, Cl^?, SO₄ ^2?, HCO₃ ^?, NO₃ ^?, pH, electrical conductivity and temperature). The microbiological analysis of 143 samples reveals the presence of four kinds of indicator bacteria in the groundwater resources: faecal Streptococci, faecal coliform, Escherichia coli and Clostridium. The physico-chemical results and bacteriological monitoring show that the nitrate and bacteria concentrations exceed the maximum admissible levels, notably around pumping stations in the sectors of Sidi Taibi, Sidi Ahmed Taleb and Aïn Sbaâ. Contamination is generated by uncontrolled anthropogenic activities and accentuated by the high intrinsic vulnerability of the aquifer system. Several parameters appeared to exceed admissibility standards. Measures are recommended to prevent groundwater pollution in the region.     

Groundwater Quality: Analysis of Its Temporal and Spatial Variability in a Karst Aquifer

This study develops an approach based on hierarchical cluster analysis for investigating the spatial and temporal variation of water quality governing processes. The water quality data used in this study were collected in the karst aquifer of Yucatan, Mexico, the only source of drinking water for a population of nearly two million people. Hierarchical cluster analysis was applied to the quality data of all the sampling periods lumped together. This was motivated by the observation that, if water quality does not vary significantly in time, two samples from the same sampling site will belong to the same cluster. The resulting distribution maps of clusters and box‐plots of the major chemical components reveal the spatial and temporal variability of groundwater quality. Principal component analysis was used to verify the results of cluster analysis and to derive the variables that explained most of the variation of the groundwater quality data. Results of this work increase the knowledge about how precipitation and human contamination impact groundwater quality in Yucatan. Spatial variability of groundwater quality in the study area is caused by: a) seawater intrusion and groundwater rich in sulfates at the west and in the coast, b) water rock interactions and the average annual precipitation at the middle and east zones respectively, and c) human contamination present in two localized zones. Changes in the amount and distribution of precipitation cause temporal variation by diluting groundwater in the aquifer. This approach allows to analyze the variation of groundwater quality controlling processes efficiently and simultaneously.     

Application of Large‐Scale Inversion Algorithms to Hydraulic Tomography in an Alluvial Aquifer

Large‐scale inversion methods have been recently developed and permitted now to considerably reduce the computation time and memory needed for inversions of models with a large amount of parameters and data. In this work, we have applied a deterministic geostatistical inversion algorithm to a hydraulic tomography investigation conducted in an experimental field site situated within an alluvial aquifer in Southern France. This application aims to achieve a 2‐D large‐scale modeling of the spatial transmissivity distribution of the site. The inversion algorithm uses a quasi‐Newton iterative process based on a Bayesian approach. We compared the results obtained by using three different methodologies for sensitivity analysis: an adjoint‐state method, a finite‐difference method, and a principal component geostatistical approach (PCGA). The PCGA is a large‐scale adapted method which was developed for inversions with a large number of parameters by using an approximation of the covariance matrix, and by avoiding the calculation of the full Jacobian sensitivity matrix. We reconstructed high‐resolution transmissivity fields (composed of up to 25,600 cells) which generated good correlations between the measured and computed hydraulic heads. In particular, we show that, by combining the PCGA inversion method and the hydraulic tomography method, we are able to substantially reduce the computation time of the inversions, while still producing high‐quality inversion results as those obtained from the other sensitivity analysis methodologies.     

Nitrate as a tracer of groundwater flow in a fractured multilayered aquifer

Abstract

Knowledge of the processes that control nitrate migration and its geochemical evolution in the subsurface are fundamental for the regional management of polluted aquifers. In this paper, the spatial distribution and transient variations of nitrate concentrations, associated with manure fertilization, are used to depict hydrogeological dynamics within the sedimentary aquifer system of la Plana de Vic in the Osona region of Spain. Flow systems are identified from geological, hydraulic head, hydrochemical and isotopic data, and by considering nitrate itself as a tracer that indicates how flow paths are modified by human pressures. In this area, nitrates move through fractured aquitards in flows induced by groundwater pumping. Moreover, the lack of casing in the boreholes permits a mixing of groundwater from distinct layers inside the wells, which negates any benefits from the low-nitrate groundwater found in the deepest aquifer layers. Therefore, impacts on groundwater quality are related to intensive manure fertilization as well as to inadequate well construction and exploitation strategies.

Citation Menció, A., Mas-Pla, J., Otero, N. & Soler, A. (2011) Nitrate as a tracer of groundwater flow in a fractured multilayered aquifer. Hydrol. Sci. J. 56(1), 108–122.     

Comparative hydrogeology – reference analysis of groundwater dynamics from neighbouring observation wells

ABSTRACT

Groundwater level fluctuations are caused by spatial and temporal superposition of processes within and outside the aquifer system. Most of the subsurface processes are usually observed on a small scale. Upscaling to the regional scale, as required for future climate change scenarios, is difficult due to data scarcity and increasing complexity. In contrast to the limited availability of system characteristics, high-resolution data records of groundwater hydrographs are more generally available. Exploiting the information contained in these records should thus be a priority for analysis of the chronical lack of data describing groundwater system characteristics. This study analyses the applicability of 63 indices derived from daily hydrographs to quantify different dynamics of groundwater levels in unconfined gravel aquifers from three groundwater regions (Bavaria, Germany). Based on the results of two different skill tests, the study aids index selection for different dynamic components of groundwater hydrographs.     

Application of multivariate statistical procedures to the hydrochemical study of a coastal aquifer: an example from Crete,Greece

A variety of multivariate statistical procedures were applied to three separate sets of quantitative analytical data from a coastal aquifer located in Malia, Crete (Greece), in order to identify the major hydrochemical processes affecting the groundwater quality and to investigate the evolution of groundwater composition in three different sampling periods. Two of them were carried out on October 2001 and September 2002 at the end of the dry season and the third on April 2002 at the end of the wet period. Two factors were found that explained major hydrochemical processes in the aquifer. These factors reveal the existence of an intensive intrusion of seawater and mechanisms of nitrate contamination of groundwater. Bivariate plots of the scores of the two main factors showed that the seawater intrusion and nitrate pollution processes are maintained through three surveys and that the process of nitrate pollution increases from the first to the second dry survey. Q‐mode factor analysis and discriminant analysis of the three sampling periods clearly showed a seasonal variation of the whole chemistry of groundwater samples. This seasonal variation can be attributed to the freshwater recharge and seawater intrusion that affect the groundwater quality of the Malia aquifer. The results of trend surface analysis are in agreement with those of factor analysis. Moreover, the fourth‐order trend surfaces of EC, Cl^? and NO₃^? showed that the salinization process is more intensive during the first dry period and the spatial variation of NO₃^? maxima plumes are strongly affected by the flow regime of the Malia aquifer. Copyright © 2007 John Wiley & Sons, Ltd.     

Geostatistical simulation of spatial variability of convective storms in Mexico City Valley

Precipitation is one of the main components of the hydrological cycle and knowledge of its spatial distribution is fundamental for the prediction of other closely related environmental variables, for example, runoff, flooding and aquifer recharge. Most of the precipitation in Mexico City is due to convective storms characterized by a high spatial variability, implying that modeling its behavior is very complex. In this work stochastic simulation techniques with a geostatistical approach were applied to model the spatial variability of the rainfall of three convective storms. The analysis of the results shows that using the proposed methodology spatial distributions of rain are obtained that reproduce the statistical characteristics present in the available information.     

Saltwater intrusion modelling in Jorf coastal aquifer,South‐eastern Tunisia: geochemical,geoelectrical and geostatistical application

Marine intrusion is the most serious problem facing the coastal Jorf shallow aquifer, located in south‐eastern Tunisia on the Mediterranean Sea. Jorf Aquifer is intensively exploited to supply the growing needs of agriculture and domestic sectors. This work proposes a multidisciplinary investigation, involving hydro‐geochemical, geoelectrical survey and geostatistical techniques for modelling the saltwater intrusion. For this purpose, 36 water samples were conducted and analysed. Electric conductivity, pH, total dissolved solids and major ions were measured and analysed. Pie and Durov Diagrams, Q‐mode hierarchical cluster and geostatistical analysis were considered to identify the main groundwater mineralization processes. Results revealed that the Na‐Cl‐Ca‐SO₄ is the dominant water type suggesting that dissolution of halite and gypsum was the main mineralization source of groundwater in the central and southern part of study area. However, saltwater intrusion was shown to control groundwater quality essentially in coastal areas. Variographic analyses were used to select the variographic model that best fits the spatial development of apparent resistivity. Kriged apparent resistivity profiles showed an abnormal decrease of resistivity values in the coastal zone, implying highly saline water because of seawater intrusion. Apparent resistivity values also decrease considerably in the faulted areas, suggesting a contribution of faults to seawater intrusion. Finally, saltwater mixing ratio was computed for each sample, and a refined seawater intrusion map was developed. Copyright © 2012 John Wiley & Sons, Ltd.     

Acquisition de salinité et qualité des eaux d’une nappe profonde en Tunisie: approche statistique et géochimique

Abstract

The multi-layered Jeffara de Gabes aquifer system is greatly influenced by tectonics. This system is limited at the base and laterally by evaporite layers and has lateral contacts with the sebkhas (salt flats). The groundwater in this aquifer is characterized by high salinity (3–10 g L^-1). Multivariate statistical analysis and a geochemical approach were applied to determine the influence of the evaporite layers and sebkhas on the hydrochemical quality of the Jeffara de Gabes aquifer, and to understand the processes governing its salinity. According to these methods, and based in part on the Sr²⁺/Ca²⁺ ratio, it is demonstrated that the strong salinity of the groundwater is due to interactions between water and the evaporite layers that act as a substratum of this aquifer, as well as saltwater intrusion from the sebkhas. Moreover, the medium- to poor-quality groundwaters are characterized by geochemical interactions: cationic exchange and the precipitation/dissolution process of minerals in the aquifer formations.

Editeur Z.W. Kundzewicz

Citation Ben Alaya, M., Zemni, T., Mamou, A. et Zargouni, F., 2014. Acquisition de salinité et qualité des eaux d’une nappe profonde, Tunisie: approche statistique et géochimique. Hydrological Sciences Journal, 59 (2), 395–419.     

Geostatistical Borehole Image‐Based Mapping of Karst‐Carbonate Aquifer Pores

Quantification of the character and spatial distribution of porosity in carbonate aquifers is important as input into computer models used in the calculation of intrinsic permeability and for next‐generation, high‐resolution groundwater flow simulations. Digital, optical, borehole‐wall image data from three closely spaced boreholes in the karst‐carbonate Biscayne aquifer in southeastern Florida are used in geostatistical experiments to assess the capabilities of various methods to create realistic two‐dimensional models of vuggy megaporosity and matrix‐porosity distribution in the limestone that composes the aquifer. When the borehole image data alone were used as the model training image, multiple‐point geostatistics failed to detect the known spatial autocorrelation of vuggy megaporosity and matrix porosity among the three boreholes, which were only 10 m apart. Variogram analysis and subsequent Gaussian simulation produced results that showed a realistic conceptualization of horizontal continuity of strata dominated by vuggy megaporosity and matrix porosity among the three boreholes.     

Statistical and geostatistical investigations into the effects of the Gabcikovo hydropower plant on the groundwater resources of northwest Hungary / Analyses statistiques et géostatistiques des effets de la centrale hydroélectrique de Gabcikovo sur les ressources en eaux souterraines du Nord-Ouest de la Hongrie

Abstract

Abstract The construction of the Gabcikovo hydropower plant and the diversion of the Danube River over 25 km into an artificial channel in 1992 influenced the groundwater regime of the region considerably. Statistical and geostatistical methods are used to quantify changes of different groundwater characteristics on the Hungarian side of the river based on observations in the time period 1960–2000. External drift kriging was used to interpolate groundwater levels and the other related variables. While mean groundwater levels did not change appreciably, there are significant changes in the variability. Standard deviations of the groundwater levels and the amplitude of the annual cycle decreased near the old river bed of the Danube. The water-level fluctuations of the Danube have a decreased influence on the groundwater dynamics. Interrelationships between water levels in wells have also changed.     

Lithological control of land subsidence induced by groundwater withdrawal in new urban AREAS (Granada Basin,SE Spain). Multiband DInSAR monitoring

Ground subsidence in the southeastern border of the Granada Basin (SE Spain) has been studied using remote sensing techniques. Over the last decades, the region has experienced a huge urban expansion, which has caused a substantial increase in water supply requirements. Water needs are exclusively met by groundwater by means of numerous pumping wells, which exploit a confined detrital aquifer of alluvial fan deposits with a heterogeneous facies distribution. A general piezometric level decline (up to 50 m) has been recorded in the aquifer during the past 30 years that has induced the generation of a subsiding area with oval shape oriented WNW‐ESE just where the new urban areas and pumping wells are located. Subsidence has been monitored by exploiting synthetic aperture radar (SAR) images from ENVISAT (2003–2009) and Cosmo‐SkyMed (2011–2014). A new approach, which combines A‐DInSAR and small‐area persistent scatterer interferometry (PSI) analysis, has been applied obtaining a good accuracy regarding temporal and spatial dimension of the subsidence. ENVISAT data (2003–2009) reveal subsidence rates up to 10–15 mm/year, and Cosmo‐SkyMed (2011–2014) values slightly lower; up to 10 mm/year. Temporal variations in the subsidence velocity are in accordance with the rainfall pattern and piezometric fluctuations in the aquifer. The sector with highest rates of subsidence does not correspond to the area with more intense groundwater exploitation but to the area with greater presence of clays in the confining layer of the aquifer. There is a clear lithological control in the spatial distribution of the ground subsidence. This work integrates detailed geological and hydrogeological data with differential SAR interferometry monitoring with the aim to better understand subsidence processes in detrital aquifers with small‐scale heterogeneity. Copyright © 2016 John Wiley & Sons, Ltd.     

Quality Assessment of Southeast Nile Delta Groundwater for Irrigation

The Nile Delta aquifer has deteriorated in the quality of the groundwater due to domestic, agricultural and industrial activities. In order to examine this, a dataset of thirty-one shallow groundwater samples and four surface water samples were collected in May 2014. The objective of our study is to investigate the hydrochemical characteristics of the groundwater at El-Khanka region in El-Qalubia governorate, southern Nile Delta to discuss the possibility of groundwater use for agricultural purpose. Groundwater types were defined, and the suitability for use in irrigation was evaluated. The factor analysis was conducted to investigate the relationship between the thirteen variables for exploring the loading of them in the model. Then, the principal component analysis was performed to identify the linear combination of variables that account for the greatest amount of common variance. Results showed that groundwater samples are mainly alkaline with an average pH value of 8.60. The total dissolved solids (TDS) range from 350 to 1456 mg/L. The highest concentrations of the anions and cations are sulfate (\(\rm{SO}_4^{2-}\)) and sodium (Na⁺) respectively. The residual sodium carbonate (RSC) is less than 1.25 meq/L. Also, all groundwater samples are located in good and permissible salinity with TDS < 1500 mg/L. In addition, all samples are located in the low sodium hazard zone where sodium adsorption ratio (SAR) is less than 10. Therefore, it is concluded that, the groundwater is suitable for irrigation use in El-Qalubia Governorate. Four factors with Eigenvalues above 1.0 which correlated to each other contributed to the model with 81% of the total variance and governed the spatial variability of the aquifer.     

A conceptual–numerical model to simulate hydraulic head in aquifers that are hydraulically connected to surface water bodies

In this paper, we present a conceptual‐numerical model that can be deduced from a calibrated finite difference groundwater‐flow model, which provides a parsimonious approach to simulate and analyze hydraulic heads and surface water body–aquifer interaction for linear aquifers (linear response of head to stresses). The solution of linear groundwater‐flow problems using eigenvalue techniques can be formulated with a simple explicit state equation whose structure shows that the surface water body–aquifer interaction phenomenon can be approached as the drainage of a number of independent linear reservoirs. The hydraulic head field could be also approached by the summation of the head fields, estimated for those reservoirs, defined over the same domain set by the aquifer limits, where the hydraulic head field in each reservoir is proportional to a specific surface (an eigenfunction of an eigenproblem, or an eigenvector in discrete cases). All the parameters and initial conditions of each linear reservoir can be mathematically defined in a univocal way from the calibrated finite difference model, preserving its characteristics (geometry, boundary conditions, hydrodynamic parameters (heterogeneity), and spatial distribution of the stresses). We also demonstrated that, in practical cases, an accurate solution can be obtained with a reduced number of linear reservoirs. The reduced computational cost of these solutions can help to integrate the groundwater component within conjunctive use management models. Conceptual approximation also facilitates understanding of the physical phenomenon and analysis of the factors that influence it. A simple synthetic aquifer has been employed to show how the conceptual model can be built for different spatial discretizations, the parameters required, and their influence on the simulation of hydraulic head fields and stream–aquifer flow exchange variables. A real‐world case was also solved to test the accuracy of the proposed approaches, by comparing its solution with that obtained using finite‐difference MODFLOW code. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical groundwater modelling as an effective tool for management of water resources in arid areas

Abstract

Groundwater, possibly of fossil origin, is used for water supply in some arid regions where the replenishment of groundwater by precipitation is low. Numerical modelling is a helpful tool in the assessment of groundwater resources and analysis of future exploitation scenarios. To quantify the groundwater resources of the East Owienat area in the southwest of the Western Desert, Egypt, the present study assesses the groundwater resources management of the Nubian aquifer. Groundwater withdrawals have increased in this area, resulting in a disturbance of the aquifer’s natural equilibrium, and the large-scale and ongoing depletion of this critical water reserve. Negative impacts, such as a decline in water levels and increase in salinity, have been experienced. The methodology includes application of numerical groundwater modelling in steady and transient states under different measured and abstraction scenarios. The numerical simulation model developed was applied to assess the responses of the Nubian aquifer water level under different pumping scenarios during the next 30 years. Groundwater management scenarios are evaluated to find an optimal management solution to satisfy future needs. Based on analysis of three different development schemes that were formulated to predict the future response of the aquifer under long-term water stress, a gradual increase in groundwater pumping to 150% of present levels should be adopted for protection and better management of the aquifer. Similar techniques could be used to improve groundwater management in other parts of the country, as well as other similar arid regions.

Editor D. Koutsoyiannis; Associate editor X. Chen     

Influence of precipitation and deep saline groundwater on the hydrological systems of Mediterranean coastal plains: a general overview

Abstract

The increasing water demand is a concern affecting many regions in the Mediterranean Basin. To overcome this situation rim countries resorted during the last decades to a massive mobilization of their water resources, often resulting in excessive water exploitation. In such a context, understanding the effects of present recharge and aquifer salinization is crucial for correct water management. Understanding the present hydrogeological situation of coastal plains requires the knowledge of both their past morphologic conditions and their recent geological evolution. Within this framework, this paper presents a review of water related problems in the Mediterranean Basin. It suggests a conceptual model for groundwater resources in Mediterranean coastal plains, deriving from the present and past recharge processes. Special attention is paid to providing a better understanding of climate change impacts on water quantity and quality, and conservation of ecological diversity.

Citation Re, V. & Zuppi, G. M. (2011) Influence of precipitation and deep saline groundwater on the hydrological systems of Mediterranean coastal plains: a general overview. Hydrol. Sci. J. 56(6), 966–980.     

Recharge sources and hydrogeological effects of irrigation and an influent river identified by stable isotopes in the Motril‐Salobreña aquifer (Southern Spain)

Measurement of the stable isotopes oxygen‐18 and deuterium in water is an important tool to characterize aquifer recharge sources. In the driest areas of the Mediterranean, this application is of special interest due to the scarcity of water and the resulting common incidence of human influence on natural hydrological systems. The Motril‐Salobreña detrital aquifer (southern Spain) is a clear example of such an impact as inhabitants have designed irrigation systems and a dam was recently built across the course of the Guadalfeo River, which feeds the aquifer. The sampling of (river or ground) water has allowed the determination of stable isotope contents (oxygen‐18 and deuterium), both temporally and spatially, and the relative importance of the main recharge sources in certain sectors. In addition, we were able to infer seasonal trends and to improve existing knowledge of the main flow paths and the position of a seasonal groundwater divide. Data analysis shows evaporation plays a minor role (despite the high temperatures in the zone), scarce rainwater influence, and the overwhelming contribution of recharge from the Guadalfeo River and from the carbonate aquifer (Escalate aquifer) in contact with the Motril‐Salobreña aquifer. Irrigation return flow during the summer months comprises the main recharge due to the significant volumes of water that infiltrate. The construction of the dam will almost certainly entail great changes in the current dynamics of the hydrogeology of the Motril‐Salobreña aquifer; therefore, knowledge of its behaviour is crucial in order to carry out sustainable use of its groundwater resources. Copyright © 2011 John Wiley & Sons, Ltd.     

Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases,southern Tunisia: implications of return flow

Abstract

Major ions and stable isotopes in groundwaters of the Plio-Quaternary shallow aquifer of the Djerid oases, southern Tunisia, were investigated to elucidate the origin of groundwater recharge and the mineralization processes. It has been demonstrated that the groundwater composition is mainly controlled by the water–rock interaction, the encroachment of brines from the Chotts and the return flow of irrigation waters. The isotopically depleted groundwater samples suggest that the recharge waters derive from an old palaeoclimatic origin. However, the enriched groundwater samples reflect the presence of evaporated recharge water. Furthermore, the large negative deuterium-excess values indicate the effect of secondary evaporation processes, probably related to the return flow of irrigation waters pumped from the underlying aquifer.

Editor D. Koutsoyiannis; Associate editor E. Custodio

Citation Tarki, M., Dassi, L. and Jedoui, Y., 2012. Groundwater composition and recharge origin in the shallow aquifer of the Djerid oases, southern Tunisia: implications of return flow. Hydrological Sciences Journal, 57 (4), 790–804.