Predicting saltwater intrusion into aquifers in vicinity of deserts using spatio-temporal kriging

The primary objective of this study was to provide a detailed framework to use the spatio-temporal kriging to model the spatio-temporal variations of salinity data and predict saltwater intrusion into freshwater aquifers in the vicinity of deserts. EC data, measured in extraction wells in the Mahvelat plain located in the Northeastern part of Iran, available from 2007 to 2013, were used to demonstrate the developed framework. The source of data was not a well-designed measurement network. Therefore, to homogenize the data, spatial analysis was used to find EC distribution in the area in each year of study. To conduct the spatial analysis, a guideline and a systematic process were developed to select an appropriate kriging method and optimize its parameters. This process can be applied to different variables. After spatial analysis of EC data for all the years of the analysis period using empirical Bayesian kriging (EBK) method with manually optimized parameters, spatio-temporal and corresponding variogram analysis was conducted using R software. This process was based on a separable product-sum model applied to the data from 2007 to 2012. The data of 2013 and the data available for the years 1999 and 2006 were used for evaluating the performance of the spatio-temporal model. The EC distribution maps, developed for different years until 2021, show a high level of EC in the north, south, and west of the study area and growing saltwater intrusion into the central freshwater aquifer. This result can be attributed to the over-exploitation of the aquifer and hydraulic head and gradient distribution in the area. The framework provided in this study for spatio-temporal analysis of unstructured EC data is useful for groundwater managers in making proper decisions.