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-HCO3 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)).
