A Permanent Multilevel Monitoring and Sampling System in the Coastal Groundwater Mixing Zones

To study the spatial and temporal variability of water dynamics and chemical reactions within the coastal groundwater mixing zones (CGMZs), high‐resolution periodical and spatial groundwater sampling within CGMZs is needed. However, current samplers and sampling systems may require heavy driving machines to install. There is also possible contamination from the metal materials for current samplers and sampling systems. Here, a permanent multilevel sampling system is designed to sample coastal groundwater within CGMZs. This cost‐effective system consists of metal‐free materials and can be installed easily. The system is tested in Po Sam Pai and Tingkok, Tolo Harbor and Hong Kong. Major ions, nutrients, stable isotopes and radium and radon isotopes were analyzed and the data provided scientific information to study the fresh‐saltwater interface fluctuations, and temporal variations and spatial heterogeneity of geochemical processes occurred within CGMZs. The reliable spatial and temporal data from the sampling system demonstrate that the system functions well and can provide scientific data for coastal aquifer studies.     

Spatial and temporal analysis of groundwater recharge with application to sampling design

A numerical experiment of flow in variably saturated porous media was performed in order to evaluate the spatial and temporal distribution of the groundwater recharge at the phreatic surface for a shallow aquifer as a function of the input rainfall process and soil heterogeneity. The study focused on the groundwater recharge which resulted from the percolation of the excess rainfall for a 90-days period of an actual precipitation record. Groundwater recharge was defined as the water flux across the moving phreatic surface. The observed spatial non-uniformity of the groundwater recharge was caused by soil heterogeneity and is particularly pronounced during the stage of recharge peak (substantial percolation stage). During that stage the recharge is associated with preferential flow paths defined as soil zones of locally higher hydraulic conductivity. For the periods of low percolation intensity the groundwater recharge was exhibiting more uniform spatial characteristics. The temporal distribution of the recharge was found to be a function of the frequency and intensity of the rainfall events. Application of sampling design demonstrates the joint influence of the spatial and temporal recharge variability on the cost-effective monitoring of groundwater potentiometric surfaces.     

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.     

Mapping groundwater quality variables using PCA and geostatistics: a case study of Bajo Andarax,southeastern Spain

Abstract

A case study is presented for the application of statistical and geostatistical methods to the problem of estimating groundwater quality variables. This methodology has been applied to the investigation of the detrital aquifer of the Bajo Andarax (Almería, Spain). The use of principal components analysis is proposed, as a first step, for identifying relevant types of groundwater and the processes that bring about a change in their quality. As a result of this application, three factors were obtained, which were used as three new variables (VI: sulphate influence; V2: thermal influence; and V3: marine influence). Analysis of their spatial distribution was performed through the calculation of experimental and theoretical variograms, which served as input for geostatistical modelling using ordinary block kriging. This analysis has allowed a probabilistic representation of the data to be obtained by mapping the three variables throughout the aquifer for each sampling point. In this way, one can evaluate the spatial and temporal variation of the principal physico-chemical processes associated with the three variables VI, V2 and V3 implicated in the groundwater quality of the detrital aquifer.     

Applying a Regional Transport Modeling Framework to Manage Nitrate Contamination of Groundwater

Regional nitrate contamination in groundwater is a management challenge involving multisector benefits. There is always conflict between restricting anthropogenic activities to protect groundwater quality and prioritizing economic development, especially in productive agriculture dominated areas. To mitigate the nitrate contamination in groundwater, it is necessary to develop management alternatives that simultaneously support environmental protection and sustainable economic development. A regional transport modeling framework is applied to evaluate nitrate fate and transport in the Dagu Aquifer, a shallow sandy aquifer that supplies drinking water and irrigation water for a thriving agricultural economy in Shandong Province in east coastal China. The aquifer supports intensive high-value vegetable farms and nitrate contamination is extensive. Detailed land-use information and fertilizer use data were compiled and statistical approaches were employed to analyze nitrogen source loadings and the spatiotemporal distribution of nitrate in groundwater to support model construction and calibration. The evaluations reveal that the spatial distribution and temporal trends of nitrate contamination in the Dagu Aquifer are driven by intensive fertilization and vertical water exchange, the dominant flow pattern derived from intensive agricultural pumping and irrigation. The modeling framework is employed to assess the effectiveness of potentially applicable management alternatives. The predictive results provide quantitative comparisons for the trend and extent of groundwater quality mitigation under each scenario. Recommendations are made for measures that can both improve groundwater quality and sustain productive agricultural development.     

Multivariate statistical analysis and environmental isotopes of Amman/Wadi Sir (B2/A7) groundwater,Yarmouk River Basin,Jordan

Multivariate statistical techniques, cluster and factor analyses were applied on the Amman/Wadi Sir groundwater chemistry, Yarmouk River basin, north Jordan. The main objective was to investigate the main processes affecting the groundwater chemical quality and its evolution. The k‐means cluster analysis yields three groups with distinct ionic concentrations. Cluster 1 comprises the vast majority of the sampled wells, and the water that belongs to this cluster can be classified as freshwater. Cluster 2 comprises only 2% of the sampled wells; it has the highest ionic concentration. The water of this cluster can be classified as brackish water. Cluster 3 involves 23% of the sampled wells, and it has total ionic concentration intermediate to that of clusters 1 and 2. Factor analysis yields a three‐factor model, which explains 76.77% of the groundwater quality variation. Factor 1 ‘salinity factor’ involves EC, Na⁺, Cl^‐, SO₄^‐2, K⁺ and Mg⁺² and reflects groundwater salinization because of overpumping. Factor 2 ‘hardness factor’ includes Ca⁺², HCO₃^‐ and the pH value and signifies soil–water/rock interaction. Factor 3 ‘nitrate factor’ involves only NO₃^‐ and points to groundwater contamination because of human activities, mainly untreated wastewater, and crops and animal cultivation in the unconfined portion of the aquifer. Factors 1 and 3 can be described as human‐induced factors, whereas factor 2 can be described as geogenic factor. Factors' scores were mapped to deduce the controlling processes on the groundwater chemistry. Stable isotope composition of ¹⁸O and ²H has revealed that the groundwater is a mixture of two water types. The radioactive isotopes tritium and ¹⁴ C were used to evaluate present day recharge to the aquifer and to estimate the groundwater age, respectively. Present day recharge to the groundwater is taking place in the unconfined portion of the aquifer as it is indicated by the measurable tritium content and low groundwater age. Copyright © 2012 John Wiley & Sons, Ltd.     

Groundwater Monitoring Using GRACE and GLDAS Data after Downscaling Within Basaltic Aquifer System

Gravity Recovery and Climate Experiment (GRACE) satellite mission is ground-breaking information hotspot for the evaluation of groundwater storage. The present study aims at validating the sensitivity of GRACE data to groundwater storage variation within a basaltic aquifer system after its statistical downscaling on a regional scale. The basaltic aquifer system which covers 82.06% area of Maharashtra state in India, is selected as the study area. Five types of basaltic aquifer systems with varying groundwater storage capacities, based on hydrologic characteristics, have been identified within the study area. The spatial and seasonal trend analysis of observed in situ groundwater storage anomalies (ΔGWSano) computed from groundwater level data of 983 wells from the year 2002 to 2016, has been performed to analyze the variation in groundwater storages in the different basaltic aquifer system. The groundwater storage anomalies (ΔGWSDano) have been derived from GRACE Release 05 (RL05) after removing the soil moisture anomaly (ΔSMano) and canopy water storage anomaly (ΔCNOano) obtained from Global Land Data Assimilation System (GLDAS) land surface models (NOAH, MOSAIC, CLM and VIC). The artificial neural network technique has been used to downscale the GRACE and GLDAS data at a finer spatial resolution of 0.125°. The study shows that downscaled GRACE and GLDAS data at a finer spatial resolution is sensitive to seasonal groundwater storage variability in different basaltic aquifer systems and the regression coefficient R has been found satisfactory in the range of 0.696 to 0.818.     

Teleconnection and climatic oscillation in aquifer water level in Kumamoto plain,Japan

Relative little is known about the interaction between climate change and groundwater. Analysis of aquifer response to climatic variability could improve the knowledge related to groundwater resource variations and therefore provides guidance on water resource management. In this work, seasonal and annual variations of groundwater levels in Kumamoto plain (Japan) and their possible interactions with climatic indices and El Niño Southern Oscillation (ENSO) were analyzed statistically. Results show the following: (1) The water level in the recharge area mainly fluctuates at 1‐ and 2‐year periods, whereas the significant periodicity for water level oscillation in the coastal aquifer is 0.5 year. (2) The aquifer water levels are possibly influenced by variability in precipitation, air temperature, barometric pressure, humidity variances and ENSO. Relative high correlations and large proportions of similarities in wavelet power patterns were found between these variables and water levels. (3) Aquifer response to climatic variances was evaluated using cross wavelet transform and wavelet coherence. In recharging aquifers, the ENSO‐induced annual variations in precipitation, air temperature, humidity and barometric pressure affect aquifer water levels. The precipitation, air temperature and humidity respond to ENSO with a 4‐, 6‐ and 8‐month time lag, respectively, whereas the ENSO imparts weak influence on the barometric pressure. Significant biennial variation of water levels during 1991–1995 is caused primarily by precipitation and humidity variations. In the coastal aquifer, the 0.5‐year variability in ENSO is transferred by precipitation, barometric pressure and humidity to aquifer water levels, and the precipitation/humidity influence is more significant comparing with the barometric pressure. Copyright © 2014 John Wiley & Sons, Ltd.     

Weather-induced temporal variations in nitrate concentrations in shallow groundwater

For the evaluation of policy action programs to improve groundwater quality, research institutes and governments intensively monitor nitrate concentrations in shallow or near surface groundwater. However, trend detection is often hampered by the large seasonal and multi-annual temporal variability in nitrate concentrations, especially in shallow groundwater within 0–5 m below the surface in relatively humid regions. This variability is mainly caused by variations in precipitation excess (precipitation minus evapotranspiration) that results in strong variability in groundwater recharge. The objective of this study was to understand and quantify this weather-induced variability in shallow groundwater nitrate concentrations.We present an example of measured weather related variations in shallow groundwater nitrate concentrations from De Marke, an intensively monitored experimental farm in The Netherlands. For the quantification of the weather-induced variability, concentration-indices were calculated using a 1D model for water and solute transport. The results indicate that nitrate concentrations in the upper meter of groundwater at De Marke vary between 55% and 153% of the average concentration due to meteorological variability. The concentration-index quantification method was successfully used to distinguish weather related variability from human-induced trends in the nitrate concentration monitoring data from De Marke. Our model simulations also shows that sampling from fixed monitoring wells produces less short term variability than measuring from open boreholes. In addition, using larger screen depths and longer screens filters out short term temporal variability at the cost of a more delayed detection of trends in groundwater quality.     

FEUWAnet: a multi-box water level and lateral exchange model for riparian wetlands

Riparian wetlands as typical aquatic-terrestrial interfaces control, in a very specific way, nonpoint water and related chemical fluxes exchanging between catchment areas to their respective water systems (streams, lakes). The existing groundwater and soilwater flow models reveal gaps in dealing with the complex behaviour of processes and the considerable spatial and temporal heterogeneity of riparian wetlands. Based on long-term experience gained through field observations and the interpretation of model produced data, a multi-box aggregation of processes which determines lateral as well as vertical flows and, as a whole, water balance, is used to discretise a generic riparian wetland transect situated between an upland aquifer and a receiving water body.

The resulting mathematical model, FEUWAnet, endowed also with an original methodology to adapt parameters, has been applied to a riparian alder wetland adjacent to Lake Belau (northern Germany). Results of simulations illustrate a good fit between calculated water levels and observed values and an accordance of calculated water balance to previous independent evaluations. This confirms that the sound simplifications of real situations performed by the FEUWAnet mathematical model are a promising way to deal with hydrological complexity of riparian zones. Moreover, FEUWAnet permits, to a certain extent, one to unravel the spatial heterogeneity and temporal variation of lateral (from catchment area to water systems) and vertical (from canopy to groundwater zone) water fluxes typical of riparian ecosystems: this is the necessary step to undertake when developing integrated models capable of assessing the effectiveness of riparian systems in controlling the fluxes of nonpoint pollution discharging in the open water bodies.     

Probabilistic risk and uncertainty analysis for bioremediation of four chlorinated ethenes in groundwater

Groundwater contamination risk assessment for health-threatening compounds should benefit from a stochastic environmental risk assessment which considers the effects of biological, chemical, human behavioral, and physiological processes that involve elements of biotic and abiotic aquifer uncertainty, and human population variability. This paper couples a complex model of chemical degradation and transformation with movement in an aquifer undergoing bioremediation to generate a health risk analysis for different population cohorts in the community. A two-stage Monte Carlo simulation has separate stages for population variability and aquifer uncertainty yielding a computationally efficient and conceptually attractive algorithm. A hypothetical example illustrates how risk variance analysis can be conducted to determine the distribution of risk, and the relative impact of uncertainty and variability in different sets of parameters upon the variation of risk values for adults, adolescents, and children. The groundwater example considers a community water supply contaminated with chlorinated ethenes. Biodegradation pathways are enhanced by addition of butyrate. The results showed that the contribution of uncertainty to the risk variance is comparable to that of variability. Among the uncertain parameters considered, transmissivity accounted for the major part of the output variance. Children were the most susceptible and vulnerable population cohort.     

Effect of temporal resolution of water level and temperature inputs on numerical simulation of groundwater–surface water flux exchange in a heavily modified urban river

Groundwater interacts with surface water features nearly in all types of landscapes. Understanding these interactions has practical consequences on the quantity and quality of water in either system, because the depletion or contamination of one of the systems will eventually affect the other one. Many studies have shown that the use of heat as natural tracer in conjunction with water level measurements is an effective method for estimating water flow (fluxes) between groundwater and surface water. A number of studies have explored the effects of spatial and temporal variability of groundwater–surface water flux exchanges using temperature and water level measurements; however, the effect of temporal resolution of water level and temperature data on estimating flux remains unexplored. Therefore, this study investigated the effect of temporal resolution of input data on temporal variation of groundwater–surface water flux exchanges. To this end, we calibrated a variably saturated two‐dimensional groundwater flow and heat transport model (VS2DH) at hourly and daily time scales using temperatures measured at multiple depths below the riverbed of the Zenne River, located at a well‐known Belgian brownfield site. Results of the study showed that the computed water flux through the streambed ranged between ?32 mm/day and +25 mm/day using the hourly model and from ?10 mm/day to ?37 mm/day using the daily model. The hourly model resulted in detecting reversal of flow direction inducing short‐term surface water flow into the streambed. However, such events were not captured if daily temperature and water level measurements were used as input. These findings have important implications for understanding contaminant mass flux and their attenuation in the mixing zone of groundwater and surface water. Copyright © 2012 John Wiley & Sons, Ltd.     

Factors Influencing Variability in Groundwater Monitoring Data Sets

“Random” variability in groundwater monitoring data sets reduces the ability to identify long‐term concentration trends. This, in turn, increases the time and cost required to evaluate the effectiveness of natural attenuation and other groundwater remedies. To better understand the factors influencing variability in groundwater monitoring results, we have analyzed three large groundwater monitoring data sets. For the three data sets, the long‐term trend in contaminant concentration in each well accounted for an average of 30% to 40% of the overall variation in contaminant concentration. Understanding the causes of the remaining variability would support the development of improved groundwater monitoring methods. All three data sets show large differences in the temporal monitoring records between individual wells (e.g., coefficient of variation for monitoring results from individual wells ranges from 0.08 to 4.6) indicating that well and aquifer factors are more important contributors to variability than sample collection and analysis factors. However, the depth to groundwater (R² = 0.020) and distance between water level and screened interval (R² = 0.049) accounted for only a portion of the differences in variability between wells and other aquifer characteristics evaluated and were not correlated with the observed variability in monitoring results. Unidentified factors were apparently much more important contributors to variability than these factors. The monitoring data sets exhibited two distinct timescales for variability: Time‐independent variability that was apparent even when wells were re‐sampled within a few days and a long‐term variability likely associated with the long‐term concentration trend. The observation of time‐independent variability suggests that frequent monitoring of contaminated monitoring wells serves primarily to characterize sources of variability unrelated to the long‐term trend of primary interest.     

Patterns and dynamics of river–aquifer exchange with variably-saturated flow using a fully-coupled model

The parallel physically-based surface–subsurface model PARFLOW was used to investigate the spatial patterns and temporal dynamics of river–aquifer exchange in a heterogeneous alluvial river–aquifer system with deep water table. Aquifer heterogeneity at two scales was incorporated into the model. The architecture of the alluvial hydrofacies was represented based on conditioned geostatistical indicator simulations. Subscale variability of hydraulic conductivities (K) within hydrofacies bodies was created with a parallel Gaussian simulation. The effects of subscale heterogeneity were investigated in a Monte Carlo framework. Dynamics and patterns of river–aquifer exchange were simulated for a 30-day flow event. Simulation results show the rapid formation of saturated connections between the river channel and the deep water table at preferential flow zones that are characterized by high conductivity hydrofacies. Where the river intersects low conductivity hydrofacies shallow perched saturated zones immediately below the river form, but seepage to the deep water table remains unsaturated and seepage rates are low. Preferential flow zones, although only taking up around 50% of the river channel, account for more than 98% of total seepage. Groundwater recharge is most efficiently realized through these zones. Subscale variability of K_sat slightly increased seepage volumes, but did not change the general seepage patterns (preferential flow zones versus perched zones). Overall it is concluded that typical alluvial heterogeneity (hydrofacies architecture) is an important control of river–aquifer exchange in rivers overlying deep water tables. Simulated patterns and dynamics are in line with field observations and results from previous modeling studies using simpler models. Alluvial heterogeneity results in distinct patterns and dynamics of river–aquifer exchange with implications for groundwater recharge and the management of riparian zones (e.g. river channel-floodplain connectivity via saturated zones).     

Groundwater Quality in the Yucatan Peninsula: Insights from Stable Isotope and Metals Analysis

High surface water-groundwater connectivity characterizes watersheds underlain by karsts, increasing contaminant transport risks. However, karsts are highly complex, making research necessary to understand the transport of contaminants from the surface, through the aquifer, to discharge areas. In Yucatan, the lack of waste water treatment raises the risk of groundwater contamination. We monitored stable isotopes (δ¹⁸O-NO₃ and δ¹⁵N-NO₃), cadmium, and lead to document waste water contamination and transport during the rainy and dry seasons, using water samples collected along the Ring of Cenotes during each season. Specific conductance and pH showed no consistent seasonality, with conductance ranging from 0.5 to 55 mS/cm and pH ranging from 6.6 to 8.6 for most samples. Nitrate concentrations in the cenotes averaged 205 ± 260 μM and no seasonal pattern was observed. Cd and Pb concentrations were 0.1 to 37.9 μg/L and 0.2 to 243.2 μg/L, respectively. Nitrate stable isotope values were 2.6 to 27.2‰ for δ¹⁸O and 1.2 to 20.7‰ for δ¹⁵N. The statistical relationship between δ¹⁵N and δ¹⁸O, in dry season samples, indicated that denitrification was occurring. A scale measure for waste water recognition showed: (1) high variability among sites probably related with dry/rainy seasons and/or diverse anthropogenic activities; and (2) specific water quality variables that contribute to contamination at each site during each season. Importantly, our analyses indicate that in the area surrounding the Ring of Cenotes, waste water exhibits spatial and temporal patterns related to complex transport and dilution processes, as is the case in karsts in general.     

Characterization of Flow Dynamics and Vulnerability in a Coastal Aquifer System

Traditional aquifer vulnerability techniques primarily rely on spatial property data for a region and are limited by their ability to directly or indirectly assess flow and transport processes occurring from the surface to depth within an aquifer system. The main objective of this study was to investigate groundwater vulnerability in terms of aquifer interconnectivity and flow dynamics. A combination of stable isotopes, groundwater age‐dating (radiocarbon), and geomorphic/geogenic spatial analyses was applied to a regional, highly developed coastal aquifer to explain the presence of nitrate at depth. The average δ¹³C value (?17.3 ± 2‰ VPDB, n = 27) is characteristic of groundwater originating from locally infiltrated precipitation through extensively cultivated soils. The average δ¹⁸O and δD values (?4.0 ± 0.1‰ VSMOW, n = 27; δD: ?19.3 ± 1‰ VSMOW, n = 27, respectively) are similar to precipitation water derived from maritime sources feeding the region's surface water and groundwater. Stable and radioactive isotopes reveal significant mixing between shallow and deep aquifers due to high velocities, hydraulic connection, and input of local recharge water to depths. Groundwater overdevelopment enhances deeper and faster modern water downward flux, amplifying aquifer vulnerability. Therefore, aquifer vulnerability is a variable, dependent on the type and degree of stress conditions experienced by a groundwater system as well as the geospatial properties at the near surface.     

Net groundwater inflow in an enclosed lake: from synoptic variations to climatic projections

Using lake Stechlin in northeastern Germany as an example of a small groundwater‐feed lake without surface inflows and outflows, we estimated the temporal scales and the variability ranges of the net groundwater contribution to the lake water budget. High‐resolution water level measurements by a bottom‐mounted pressure logger provided the background for the estimation of the total lake water budget. This method has demonstrated reliability for estimation of lake level variations during periods ranging from subdiurnal to perennial. The typical amplitudes of the synoptic‐to‐perennial variability characterizing the groundwater climate of lake Stechlin are estimated by comparing the two subsequent years 2006 and 2007; one of these years shows an extremely high, and the other an extremely low, annual precipitation–evaporation balance. The net groundwater flow, estimated as the difference between the total water budget and the precipitation–evaporation balance at the surface, revealed synoptic effects of lake water exfiltration into the groundwater aquifer following strong precipitation events. Perennial variations between wet and dry years superimposed seasonal oscillations. The probable origin of the latter is seasonality in the groundwater level on the watershed, although the exact amplitudes are subject to further quantification on account of seasonality in the evaporation estimation error. The results emphasize the non‐stationary behaviour of groundwater flow on timescales shorter than climatic ones. The analysis yielded a net quantitative relationship between groundwater flow and water balance at the lake surface: The water level changes in the lake due to evaporation and precipitation are damped to 60% because of the lake–groundwater exchange by means of intermittent infiltration and exfiltration events. Assuming the remaining 40% of the surface water budget may potentially result in perennial water level variability, we estimated an effect of the precipitation decrease on the lake water budget as predicted by the regional climate scenarios for the next century. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterization of the groundwater flow regime and hydrochemistry of groundwater from the Buem formation,Eastern Ghana

This study demonstrates the application of multivariate statistical methods in definition of groundwater recharge and discharge areas in a sedimentary basin in Ghana. Q‐mode hierarchical cluster analysis (HCA) was applied to 57 hydrochemical data from the Buem formation in the northern part of the Volta Region in Ghana. R‐mode HCA and R‐mode factor analysis were then applied to the same dataset to reveal the processes controlling the hydrochemistry of groundwater from this hydrogeological formation. Results of both the Q‐ and R‐mode analyses were backed by graphical methods. The analyses revealed two major water types, differentiated by salinity levels into four spatial groundwater associations. The characteristics of the four groundwater types are discussed. The recharge areas are characterized by Ca? HCO₃ low salinity waters which evolve through rock–water interactions to Na? HCO₃ high salinity waters in the discharge areas. This study finds that the hydrochemistry of groundwater from this formation is mainly controlled by the weathering of minerals, principally silicates in the aquifer matrix. The effects of the chemistry of recharging precipitation are higher in the recharge areas, while mineral weathering tends to be severe close to the discharge areas in the groundwater flow regime. All the four spatial groundwater associations have low sodium content, but salinity levels increase towards the discharge areas, such that some of wells in the discharge areas may not be acceptable for irrigation on grounds of high salinities which might affect the osmotic potentials of plants. Copyright © 2011 John Wiley & Sons, Ltd.     

Assessment of seasonal groundwater quality and potential saltwater intrusion: a study case in Urmia coastal aquifer (NW Iran) using the groundwater quality index (GQI) and hydrochemical facies evolution diagram (HFE-D)

Coastal aquifers are at threat of salinization in most parts of the world. This work investigated the seasonal hydrochemical evolution of coastal groundwater resources in Urmia plain, NW Iran. Two recently proposed methods have been used to comparison, recognize and understand the temporal and spatial evolution of saltwater intrusion in a coastal alluvial aquifer. The study takes into account that saltwater intrusion is a dynamic process, and that seasonal variations in the balance of the aquifer cause changes in groundwater chemistry. Pattern diagrams, which constitute the outcome of several hydrochemical processes, have traditionally been used to characterize vulnerability to sea/saltwater intrusion. However, the formats of such diagrams do not facilitate the geospatial analysis of groundwater quality, thus limiting the ability of spatio-temporal mapping and monitoring. This deficiency calls for methodologies which can translate information from some diagrams such Piper diagram into a format that can be mapped spatially. Distribution of groundwater chemistry types in Urmia plain based on modified Piper diagram using GQI_Piper(mix) and GQI_Piper(dom) indices that Mixed Ca–Mg–Cl and Ca-HCO₃ are the dominant water types in the wet and dry seasons, respectively. In this study, a groundwater quality index specific to seawater intrusion (GQI_SWI) was used to check its efficiency for the groundwater samples affected by Urmia hypersaline Lake, Iran. Analysis of the main processes, by means of the Hydrochemical Facies Evolution Diagram (HFE-Diagram), provides essential knowledge about the main hydrochemical processes. Subsequently, analysis of the spatial distribution of hydrochemical facies using heatmaps helps to identify the general state of the aquifer with respect to saltwater intrusion during different sampling periods. The HFE-D results appear to be very successful for differentiating variations through time in the salinization processes caused by saltwater intrusion into the aquifer, distinguishing the phase of saltwater intrusion from the phase of recovery, and their respective evolutions. Both GQI and HFE-D methods show that hydrochemical variations can be read in terms of the pattern of saltwater intrusion and groundwater quality status. But generally, in this case (i.e. saltwater and not seawater intrusion) the HFE-D method was presented better efficiency than GQI method (including GQI_Piper and GQI_SWI).