# GIS-SWIAS: Tool to Summarize Seawater Intrusion Status and Vulnerability at Aquifer Scale

Authors: Leticia Baena-Ruiz; David Pulido-Velazquez
Journal: Scientific Programming
 (2021)

---

## Abstract

In this paper, we introduce GIS-SWIAS, a novel generalized ArcGIS ArcToolbox that helps to analyze seawater intrusion (SWI) status and vulnerability at aquifer scale (SWIAS). It is a user-friendly tool that can be applied to any aquifer and is fully integrated in the ArcGIS environment, which is a widely available software tool. It is the first ArcGIS tool with these characteristics focusing on SWI analyses that we can find in the literature. GIS-SWIAS is able to deal with georeferenced information; it is easy to introduce the required data (inputs) and to efficiently perform the demanding computational operations required. Its outputs are in the form of shapes, reports, and images (maps, conceptual cross sections, and time series of lumped indices) to summarize the magnitude, intensity, and temporal evolution of SWI within an aquifer for specific dates or by showing statistics for a chosen time period. It can be applied to assess historical SWI dynamic in cases where there is no groundwater flow model. In those cases, the spatial distribution is assessed by applying simple interpolation techniques. Nevertheless, if we want a rational quantitative analysis of the sustainability of alternative management scenarios to the SWI problem, the GIS-SWIAS tool requires that information on hydraulic head and chloride concentration distribution is generated from simulations of their impacts by a calibrated density-dependent flow model. In such cases, adaptation strategies to potential future scenarios—whose distributed impacts have to be propagated within the previously calibrated models—could usefully be analyzed and compared using this tool. Given all these ways that the GIS-SWIAS tool can be applied, it provides a valuable tool for both the researcher and technician to assess SWI dynamics and aquifer resilience under different scenarios. It can support the decision-making process by helping to make a rational selection of sustainable management strategies. Its performance for the analyses of historical and potential future scenarios has been tested and confirmed in two case studies described in previous research works.


---

## Body

## 1. Introduction

Seawater intrusion (SWI) in coastal aquifers is a worldwide problem affecting groundwater-dependent ecosystems and human health. In recent decades, society’s awareness of this issue has grown and this has been reflected in the legal framework of many countries. For example, the European Union’s Water Framework Directive (WFD) requires that river basin plans address the achievement of a good qualitative and quantitative status of groundwater bodies [1]. In coastal groundwater bodies, intrusion is one of the main issues that need to be considered to achieve or maintain good groundwater status.The impacts of SWI on groundwater have a heterogeneous distribution. Analyses of spatiotemporal distribution of SWI require salinity or chloride concentration to be mapped for different dates. Depending on the issue to be addressed and the available information, SWI can be approximated using various models. SWI can be mapped either by applying simple interpolation models [2, 3] to existing point data or by simulating the physical processes using transient density-dependent groundwater flow models [4, 5]. The results obtained with these physical process models can be applied to assess sustainable management strategies, i.e., strategies that prevent deterioration of the aquifer resource due to SWI [6]. They can even be employed to propagate impacts of potential local climate change (CC) [7] or global change (GC) scenarios and to identify adaptation strategies [8].Based on the salinity or chloride concentration maps at different dates, some authors have defined indices to summarize SWI [9–12]. These indices provide an overview of the intensity and spatial distribution or percentage SWI at aquifer scale. Such indices need to offer descriptive and synthetic results so that the status of SWI in different aquifers and over different periods can be compared. These index-methods [9, 13] establish threshold values for chloride to define the basis of SWI that have been defined in various ways: by referencing natural background levels and/or by taking into account the concentrations required to protect dependent ecosystems or human health [14].When simple interpolation is used to draw the maps used to define the indices, analyses must be limited to the historical period for which there are data. In contrast, physical process models can be applied to propagate various potential conditions and so maps can be obtained for different scenarios (e.g., alternative management scenarios or future potential CC scenarios); in this latter case, the output can be used to determine the optimum strategy and therefore support the decision-making process [15].In addition to the SWI status and dynamics, another important issue to take into account in identifying sustainable management strategies for coastal aquifers is the aquifer’s vulnerability to SWI. In the literature we found various methods for mapping SWI vulnerability, such as the GALDIT method [16]. Then, by applying a method to express vulnerability as an index, we can also get an overview of the intensity and spatial distribution or percentage of SWI at aquifer scale [13].


## 2. Related Works

In the literature we find many tools to assess and analyze water resource problems [17–19]. The success of these software tools lies in their usability. A user-friendly environment and the implementation in commonly used software are key factors for their success and popular use. For example, groundwater studies usually employ Geographic Information Systems (GIS) since they are powerful, widely available tools that can deal with large amounts of spatial georeferenced information [20] and make calculations in an efficient way to provide quick results [21]. The analysis and mapping of hydrogeological data provide useful spatiotemporal information to decision makers [22].GIS tools have been widely used for different purposes related to groundwater issues [23, 24]. Several authors have developed different source codes within GIS environment (Alcaraz et al. [25]; Bhatt et al. [26]) for hydrogeological modelling. A free and open source module included in FREEWAT was developed by Criollo et al. [19] to analyze hydrochemical and hydrogeological data in order to simplify the characterization of the groundwater bodies at chemical risk. Almeida et al. [27] coupled a groundwater flow model into a GIS environment to simulate transient flow in a confined aquifer. Akbar et al. [28] and Ríos et al. [29] presented GIS-based models to simulate contaminants leaching into groundwater.In coastal areas, a three-dimensional GIS-based groundwater flow model was developed [30] to simulate the aquifer’s response to past climate changes. A new ArcGIS tool for groundwater flow simulation and visualization of results was also implemented by [19]. Other authors (De Filippis et al. [31]) applied a previously developed GIS-based tool (AkvaGIS), in addition to a groundwater flow model, to study the impacts of pumping on seawater intrusion in coastal aquifers in Malta and Italy. This tool was used in other studies (Perdikaki et al. [32]) to analyze hydrochemical parameters in a coastal aquifer that presented seawater intrusion problems.Nevertheless, as far as we know, there is no ArcGIS tool focusing on analyzing seawater intrusion (SWI) status and vulnerability at aquifer scale.In this paper, we describe the development of a new ArcGIS tool called GIS-SWIAS, which is the implementation of the index-based method for assessing aquifer status and vulnerability to SWI proposed by [13, 15]. It helps to analyze SWI status and/or vulnerability at aquifer scale using a mixed lumped-distributed analysis. It is a user-friendly ArcGIS® toolbox that performs all the required calculations for specific dates or temporal periods inside a GIS environment. The data inputs to the model are hydraulic head and chloride concentration maps. The tool provides two options to map these variables: the first is to use point data by applying interpolation techniques integrated within GIS-SWIAS, while the second is to take these data from existing external distributed models. The second approach allows different climatic and/or management scenarios to be assessed and compared. From those maps, extensive calculations have been fully automatized in GIS-SWIAS to display the results as distributed maps of affected and nonaffected volumes (at a specific moment or over a period of time), mean conceptual cross sections, and a lumped index (Ma and L_Vul) to analyze the global intensity and the dynamics of SWI.Although there are many GIS-based tools in the literature that allow simulating groundwater flow and analyzing groundwater quality, none of them perform spatial and temporal analysis on groundwater quality and vulnerability issues. Moreover, this new tool provides simple pictures that summarize the proportional affected area within the aquifer according to a chloride threshold. For this purpose, GIS-SWIAS has been applied to analyze the seawater intrusion problem, but this tool could be applied to represent the global status of an aquifer to any contaminant. The main objective of GIS-SWIAS is to provide an easy-to-use tool through a user-friendly interface that can be used by users at different levels of expertise to summarize the SWI problem at aquifer scale. It allows analyzing long time periods with a low computing cost.


## 3. Description of GIS-SWIAS Tool: Models, Inputs, and Outputs

GIS-SWIAS is an ArcGIS ArcToolbox that contains the models required to analyze SWI status and vulnerability at aquifer scale, according to the methodology described in previous works [13, 15]. Figure 1 shows the structure of the tool, which includes inputs, steps, and models, as well as the outputs generated.Figure 1
Workflow of the GIS-SWIAS ArcToolbox.To determine the overall status of the aquifer, theinputs to the tool include variables (to characterize the historical evolution of hydraulic head and chloride concentration) and parameters (to define aquifer geometry and hydrodynamic behaviour). Data describing the historical evolution can come from direct observations (monitoring network) or other techniques (geophysical applications, etc.). For the vulnerability assessment, an additional input is required: a distributed vulnerability index map, which comes from other intrinsic information (aquifer type, conductivity, distance from the shoreline, and bicarbonate concentration).The results/outputs produced to summarize status and vulnerability to SWI through visual displays and time series are as follows: (1) maps of aquifer volumes affected by SWI; (2) 2D conceptual cross sections (with mean penetration inland and mean thickness on specific dates, or mean values over a period of time); (3) lumped index (mass of chloride that causes the concentration in some areas to exceed the SWI threshold and lumped vulnerability index) to summarize the global dynamic of SWI within the aquifer.


### 3.1. Description of the Outputs: Theoretical Background

In order to assess the maps of SWI-affected aquifer volumes for different dates, we need to compile (A) maps of chloride concentration or vulnerability to SWI; (B) maps of groundwater volumes for specific dates; (C) threshold of chloride concentration or vulnerability that will be used to tag which parts of the aquifer are impacted (areas with chloride concentration or vulnerability index exceeding the threshold). The tool provides two options with respect to the input maps (A) (chloride or vulnerability maps) and (B), either calculating the maps internally from point data by applying interpolation techniques integrated within GIS-SWIAS or taking the maps from existing external distributed models; the second option means that various potential climatic and/or management scenarios can be compared and assessed. Maps of groundwater volume are calculated by combining hydraulic head, geometry, and the storage coefficient. The maps of groundwater volume and chloride concentration are combined to assess the aquifer volume affected by using a chloride threshold (Vr). This threshold is assumed to be equal to the natural background level of the aquifer, or the reference quality standard determined by authorities in order to maintain a good groundwater status. The affected volume is defined as the groundwater volume of resource whose chloride concentration is above the established threshold.2D conceptual cross section depicts the magnitude of the intrusion process in the aquifer at a specific moment, or the mean values in a period of time. The cross section is defined orthogonal to the shoreline. It summarizes the mean geometry of the affected volume, i.e., the mean thickness and inland penetration of the aquifer volume with chloride concentration above the threshold. The average affected thicknessTha(m) and inland penetration P(m) of the intrusion can be calculated by summing the values in each cell i of the aquifer mesh where the chloride concentration exceeds the threshold:(1)Tham=∑Vi>Vr∑Si>Vr,Pm=∑Vi>VrTha∗Lcoast,Vi>Vrm3=Sim2∗bim∗αi,where Vi(>Vr) is the groundwater volume (m3) in the cell i with a chloride concentration (or vulnerability) exceeding Vr; Si is the surface area (m2) of the cell i with chloride concentration (or vulnerability) exceeding Vr; bi is the saturated thickness (m) within the cell i with Cl concentration (or vulnerability) above Vr; αi is the storage coefficient in the cell i; Lcoast is the length of coastline (m).The mean chloride concentration (C) of the affected volume is(2)Cmgl=∑Ci>Vr∗Vi>VrV>Vr,where Ci>Vr is the chloride concentration (mg/l) in cell i; Vi(>Vr) is the groundwater volumes (m3) in cell i with a concentration exceeding Vr; V>Vr is the total groundwater volume (m3) of the affected area.The increment of chloride concentration (IC) above the threshold (Vr) in the affected volume is(3)ICmgl=C−Vr.Both variables, the conceptual cross section and IC index, give an overview of the magnitude and intensity of the intrusion process per linear metre of coast at a specific moment in time.The lumped index Ma (mass of chloride that causes the concentration in some areas to exceed the threshold) to summarize the global dynamic of SWI within the aquifer is obtained multiplying the increment of concentration (IC) by penetration (P) and affected thickness (Tha) from (1) and (3).(4)Makgm=Pm∗ICmgl∗10−3∗Tham.The vulnerability to SWI (or vulnerability to contamination in general) is assessed and summarized following the same steps to assess the SWI status. In this case, instead of chloride concentration values, a distributed map of groundwater vulnerability has to be generated by applying any index-based method (e.g., GALDIT) and the threshold applied to identify the affected area is defined by a specific vulnerability class (e.g., high or very high vulnerability).The affected volume corresponds to the groundwater that presents vulnerability values above the threshold (e.g., very high vulnerability). The average affected thickness Tha(m) and inland penetration P(m) are calculated by applying (1).The lumped index to assess vulnerability is(5)L_Vul−∑Vuli>Vr∗Vi>VrV>Vr,where Vuli>Vr is the value of the vulnerability index (-) in cell i.The concept of Ma and L_Vul involves some simplifications, in accordance with the definition of the conceptual cross sections. While 2D maps and cross sections summarize the extent and magnitude of SWI and vulnerability in an aquifer at a specific time, Ma and L_Vul indices show the lumped intensity and temporal dynamic of the SWI and vulnerability to SWI at aquifer scale.


### 3.2. Tool Programming in ArcGIS

GIS-SWIAS is an ArcGIS ArcToolbox composed of a chain of models programmed in ModelBuilder. ModelBuilder is a visual programming language that allows chaining and sequencing several geoprocessing ArcGIS tools through a user-friendly interface. ModelBuilder is available within the tool bar in ArcGIS. It allows the addition of any geoprocessing tool of ArcGIS by linking and providing its output and transferring it to another tool as input.Programming in ModelBuilder allows us to automate a workflow to create a model, which can be documented and shared as a ArcGIS tool. ModelBuilder contains a script tool to link with Python scripts and other models. It also allows iteration of a workflow, so it can be very useful to analyze the evolution of the historical hydrogeological processes.The three models that compose GIS-SWIAS have been compiled by adding different tools from ArcToolbox to create shapes from point data, to interpolate, etc. Figure2 shows the workflow of one of the three models.Figure 2
Workflow of “Chloride concentration map” model.Although ModelBuilder is an intuitive and easy-to-use tool, the integration of lots of geoprocesses in the same model may be difficult. Because several geoprocesses have dynamic parameters and the user interaction is necessary, GIS-SWIAS was divided into three steps (models) that have to be executed following the workflow shown in Figure1.


### 3.3. Description of the Models within GIS-SWIAS

GIS-SWIAS contains three ModelBuilder models (Figure2) that can be applied individually or used sequentially to produce a complete lumped assessment of the SWI at aquifer scale. GIS-SWIAS can be shared with other users and it can be added as a toolbox as shown in Figure 3. The run sequence follows the order shown in the workflow (Figure 1): “Chloride concentration map”, “Hydraulic head map,” and “Summarizing SWI”. These models can be run within the ArcToolbox window providing a user-friendly graphical user interface.Figure 3
GIS-SWIAS in the ArcToolbox.


#### 3.3.1. “Chloride Concentration map” Model

The “Chloride concentration map” model generates a classified chloride concentration shapefile from a point feature table in text format by using Inverse Distance Weighting (IDW) interpolation technique (other interpolation techniques could be implemented in this tool). It can also import chloride concentration fields from Visual MODFLOW files. The dialog box shown in Figure4 is opened by double-clicking the model tool on the ArcToolbox window.Figure 4
Dialog box of the “Chloride concentration map” model.The model requires a polygon shapefile of the aquifer and the workspace containing a text file for each date to be analyzed. The text files have to includeX and Y UTM coordinates that define the locations of the point features and chloride concentration values (mg/l) in the observation wells. Text files also have to include column header, as shown in Figure 4. The text filename cannot exceed eight characters nor include blank spaces or special characters (underscore can be used as a substitute). It is not necessary that all points contain data every date of the period to be analyzed.The user has to indicate the fields (columns) in the input table that contain theX and Y coordinates and the chloride concentration value for each point (Figure 5). Optional settings concerning IDW interpolation techniques can be changed by the user. A reclassification of the values after interpolation is also required. Finally, the user has to choose a folder where the output shapefiles will be saved.Figure 5
Setting options defined by the user.When all the required parameters are filled out, the model can be executed by clicking “OK” at the bottom of the dialog box. The execution screen (Figure6) shows the running processes and it can be closed when the execution has successfully completed.Figure 6
Execution screen of the model.The “Chloride concentration map” model provides the following shapefiles for each date analyzed: (1) point shapefile of chloride concentration data; (2) raster from data interpolation; (3) polygon shapefile from interpolation covering the default extension; and (4) polygon shapefile from interpolation clipped to the shape of the aquifer.


#### 3.3.2. “Hydraulic Head Map” Model

The “Hydraulic head map” model generates a classified hydraulic head shapefile from a point feature table in text format. It also generates a shapefile containing aquifer variables (chloride concentration and hydraulic head values) and aquifer parameters (storage coefficient and bottom of the aquifer). The dialog box is shown in Figure7.Figure 7
Dialog box of the “Hydraulic head map” model.This model has the same input requirements as for the “Chloride concentration map” model but focuses on hydraulic head data (m.a.s.l.). It also allows hydraulic head fields to be imported from Visual MODFLOW. The name of hydraulic head text files must be the same as the chloride concentration text files for each time period analyzed.The user has to indicate the location of the chloride shapefiles generated in the previous model (“Chloride concentration map” model). It also requires polygon shapefiles of storage coefficient and bottom (m) of the aquifer as inputs.The model generates shapefiles of hydraulic head data in an analogous way to the “Chloride concentration map” model. Moreover, it provides a polygon shapefile containing variables for each date analyzed (chloride concentration and hydraulic head values) and parameters (bottom and storage coefficient) of the aquifer. This shapefile is named “union_%name of the hydraulic head text file%_hh.shp”, where “% name of the hydraulic head text file%” is variable if different dates are analyzed.


#### 3.3.3. “Summarizing SWI” Model

For the “Summarizing SWI” model, the methodology proposed in Baena-Ruiz et al. [13,15] and described in Section 3.1 has been implemented in the ArcGis environment. This tool generates Excel® tables containing statistics that summarize SWI at aquifer scale. It also generates conceptual cross sections (.shp), where the mean affected and nonaffected volumes are drawn for the aquifer (average values over a time period or instantaneous values on a specific date). If different dates are analyzed, it shows graphs representing the temporal evolution of Pa and Ta variables, percentage of affected volume, chloride concentration within the aquifer, and Ma index (or lumped vulnerability index). The dialog box for global status assessment is shown in Figure 8.Figure 8
Dialog box of the “Summarizing SWI” model.The “Summarizing SWI” model requires the folder path where the results of the “Hydraulic head model” have been previously saved to be specified. In this folder, the shapefile named “union_%name of the hydraulic head text file%_hh.shp” contains chloride concentration, hydraulic head, bottom of the aquifer, and storage coefficient fields. The user has to select from the pull-down list the corresponding column in the input shapefile for each field, as shown in Figure9.Figure 9
Selection of chloride field from the input table.The next required parameter in this tool is the “Chloride threshold.” It is defined as the chloride concentration value above which the aquifer is considered to be affected by SWI. This threshold may be set as the natural background level of the aquifer or as the relevant environmental quality standards. Hinsby et al. [14] proposed a method to calculate groundwater thresholds values based on these criteria.The shoreline length (m) is also required for subsequent calculations.X and Y axes establish the coordinate system of the conceptual cross sections. The GIS-SWIAS tool provides polyline shapefiles for X and Y axes located at (0,0), but the user can translate them to another coordinate origin or create new ones.The vertical scale factor is used to rescale the vertical magnitude (Ta) of the conceptual cross section if the factor Ta/Pa is too small. If vertical scale factor = 1, the conceptual cross section will maintain the real size ratio.Finally, two paths where the output results will be saved are required. “Output workspace statistic” will contain lumped variables reports in Excel table format for each date analyzed (Figure10) and mean statistics for the entire period. Four graphs will be also saved in this path: (1) temporal evolution of Pa and Ta variables; (2) percentage volume affected; (3) Ma index; and (4) chloride concentration within the aquifer (mean chloride concentration in the aquifer, mean chloride concentration in the affected volume, and the increment of concentration within the affected volume above the threshold).Figure 10
Summary statistics for a date.“Output workspace results” will contain the polygon shapefiles that allow the (1) mean affected and (2) nonaffected conceptual cross section within the aquifer for each date analyzed to be drawn, (3) the mean affected and (4) nonaffected conceptual cross section within a time period, and the (5) maximum affected cross section for a time period. These two paths can be the same for all results, but they have to be different from the output paths of the previous models.Figure11 and Table 1 show the graphical and statistical summaries, respectively, from the GIS-SWIAS tool.Figure 11
Graphical output from the GIS-SWIAS tool.Table 1
Lumped variables output (Excel® table format) from the GIS-SWIAS tool.
Lumped variables (Excel table)At a specific moment in timeStatistics for a time periodTotal aquifer volumeAverage aquifer volumeTotal aquifer affected volumeAverage aquifer affected volumeTotal chloride concentration∗aquifer volumeAverage chloride concentration∗aquifer volumeTotal chloride concentration∗aquifer volume in the affected volumeAverage chloride concentration∗aquifer volume in the affected volumeTotal aquifer areaAverage aquifer areaTotal aquifer affected areaAverage aquifer affected areaShoreline lengthShoreline lengthAverage aquifer thicknessAverage aquifer thickness within the periodAverage affected aquifer thicknessAverage affected aquifer thickness within the periodAverage orthogonal distance of the aquifer to the shorelineAverage orthogonal distance of the aquifer to the shoreline within the periodAverage orthogonal distance of the affected aquifer to the shorelineAverage orthogonal distance of the affected aquifer to the shoreline within the periodAverage chloride concentration within the aquiferAverage chloride concentration within the aquifer within the periodAverage chloride concentration within the affected volumeAverage chloride concentration within the affected volume within the periodIncrement of chloride concentration above the threshold valueIncrement of chloride concentration above the threshold value within the periodMa indexAverage Ma index within the periodPercentage of affected volumeAverage percentage of affected volume within the periodThe GIS-SWIAS tool can provide results for each date where information is available; these are obtained by iterative application of the described method. GIS-SWIAS allows historical [13] and future periods [15] to be analyzed if the hydraulic head and chloride maps come from a density-dependent flow model. By this means, GIS-SWIAS can be used to analyze adaptation strategies [15] in terms of reducing SWI, taking into account future potential scenarios that might include CC and/or GC, also considering projected land use change scenarios (new urbanized areas, crop transformations) [15].Moreover, this tool may be also used to summarize SWI vulnerability, for any index method applied to assess it. In this case, instead of the chloride concentration maps, generated by executing the “Chloride concentration map” model, polygon shapefiles of the vulnerability index (previously prepared by the user) would be used as inputs of the model “Summarizing SWI vulnerability” (Figure12), which also will require the “Hydraulic head map” shapefile generated by the tool, as previously described.Figure 12
Dialog box of the “Summarizing SWI vulnerability” model.The vulnerability index maps must be included as numerical fields (values obtained before defining the vulnerability classes). In order to generate the conceptual cross sections that summarize the “affected” aquifer volume, i.e., where the vulnerability to SWI is identified, the tool requires a vulnerability threshold to be input that represents the reference value chosen to distinguish between affected and nonaffected volumes. This threshold will be also used to assess the lumped vulnerability index.Just as in the definition of the Ma index, the lumped global value of vulnerability in the aquifer on a specific date is obtained by weighting the vulnerability score in each cell with its water storage. This lumped index also allows an analysis of the dynamic of SWI vulnerability of the system at aquifer scale to be performed. The lumped index can be also obtained using different threshold values [13, 15].


## 3.1. Description of the Outputs: Theoretical Background

In order to assess the maps of SWI-affected aquifer volumes for different dates, we need to compile (A) maps of chloride concentration or vulnerability to SWI; (B) maps of groundwater volumes for specific dates; (C) threshold of chloride concentration or vulnerability that will be used to tag which parts of the aquifer are impacted (areas with chloride concentration or vulnerability index exceeding the threshold). The tool provides two options with respect to the input maps (A) (chloride or vulnerability maps) and (B), either calculating the maps internally from point data by applying interpolation techniques integrated within GIS-SWIAS or taking the maps from existing external distributed models; the second option means that various potential climatic and/or management scenarios can be compared and assessed. Maps of groundwater volume are calculated by combining hydraulic head, geometry, and the storage coefficient. The maps of groundwater volume and chloride concentration are combined to assess the aquifer volume affected by using a chloride threshold (Vr). This threshold is assumed to be equal to the natural background level of the aquifer, or the reference quality standard determined by authorities in order to maintain a good groundwater status. The affected volume is defined as the groundwater volume of resource whose chloride concentration is above the established threshold.2D conceptual cross section depicts the magnitude of the intrusion process in the aquifer at a specific moment, or the mean values in a period of time. The cross section is defined orthogonal to the shoreline. It summarizes the mean geometry of the affected volume, i.e., the mean thickness and inland penetration of the aquifer volume with chloride concentration above the threshold. The average affected thicknessTha(m) and inland penetration P(m) of the intrusion can be calculated by summing the values in each cell i of the aquifer mesh where the chloride concentration exceeds the threshold:(1)Tham=∑Vi>Vr∑Si>Vr,Pm=∑Vi>VrTha∗Lcoast,Vi>Vrm3=Sim2∗bim∗αi,where Vi(>Vr) is the groundwater volume (m3) in the cell i with a chloride concentration (or vulnerability) exceeding Vr; Si is the surface area (m2) of the cell i with chloride concentration (or vulnerability) exceeding Vr; bi is the saturated thickness (m) within the cell i with Cl concentration (or vulnerability) above Vr; αi is the storage coefficient in the cell i; Lcoast is the length of coastline (m).The mean chloride concentration (C) of the affected volume is(2)Cmgl=∑Ci>Vr∗Vi>VrV>Vr,where Ci>Vr is the chloride concentration (mg/l) in cell i; Vi(>Vr) is the groundwater volumes (m3) in cell i with a concentration exceeding Vr; V>Vr is the total groundwater volume (m3) of the affected area.The increment of chloride concentration (IC) above the threshold (Vr) in the affected volume is(3)ICmgl=C−Vr.Both variables, the conceptual cross section and IC index, give an overview of the magnitude and intensity of the intrusion process per linear metre of coast at a specific moment in time.The lumped index Ma (mass of chloride that causes the concentration in some areas to exceed the threshold) to summarize the global dynamic of SWI within the aquifer is obtained multiplying the increment of concentration (IC) by penetration (P) and affected thickness (Tha) from (1) and (3).(4)Makgm=Pm∗ICmgl∗10−3∗Tham.The vulnerability to SWI (or vulnerability to contamination in general) is assessed and summarized following the same steps to assess the SWI status. In this case, instead of chloride concentration values, a distributed map of groundwater vulnerability has to be generated by applying any index-based method (e.g., GALDIT) and the threshold applied to identify the affected area is defined by a specific vulnerability class (e.g., high or very high vulnerability).The affected volume corresponds to the groundwater that presents vulnerability values above the threshold (e.g., very high vulnerability). The average affected thickness Tha(m) and inland penetration P(m) are calculated by applying (1).The lumped index to assess vulnerability is(5)L_Vul−∑Vuli>Vr∗Vi>VrV>Vr,where Vuli>Vr is the value of the vulnerability index (-) in cell i.The concept of Ma and L_Vul involves some simplifications, in accordance with the definition of the conceptual cross sections. While 2D maps and cross sections summarize the extent and magnitude of SWI and vulnerability in an aquifer at a specific time, Ma and L_Vul indices show the lumped intensity and temporal dynamic of the SWI and vulnerability to SWI at aquifer scale.


## 3.2. Tool Programming in ArcGIS

GIS-SWIAS is an ArcGIS ArcToolbox composed of a chain of models programmed in ModelBuilder. ModelBuilder is a visual programming language that allows chaining and sequencing several geoprocessing ArcGIS tools through a user-friendly interface. ModelBuilder is available within the tool bar in ArcGIS. It allows the addition of any geoprocessing tool of ArcGIS by linking and providing its output and transferring it to another tool as input.Programming in ModelBuilder allows us to automate a workflow to create a model, which can be documented and shared as a ArcGIS tool. ModelBuilder contains a script tool to link with Python scripts and other models. It also allows iteration of a workflow, so it can be very useful to analyze the evolution of the historical hydrogeological processes.The three models that compose GIS-SWIAS have been compiled by adding different tools from ArcToolbox to create shapes from point data, to interpolate, etc. Figure2 shows the workflow of one of the three models.Figure 2
Workflow of “Chloride concentration map” model.Although ModelBuilder is an intuitive and easy-to-use tool, the integration of lots of geoprocesses in the same model may be difficult. Because several geoprocesses have dynamic parameters and the user interaction is necessary, GIS-SWIAS was divided into three steps (models) that have to be executed following the workflow shown in Figure1.


## 3.3. Description of the Models within GIS-SWIAS

GIS-SWIAS contains three ModelBuilder models (Figure2) that can be applied individually or used sequentially to produce a complete lumped assessment of the SWI at aquifer scale. GIS-SWIAS can be shared with other users and it can be added as a toolbox as shown in Figure 3. The run sequence follows the order shown in the workflow (Figure 1): “Chloride concentration map”, “Hydraulic head map,” and “Summarizing SWI”. These models can be run within the ArcToolbox window providing a user-friendly graphical user interface.Figure 3
GIS-SWIAS in the ArcToolbox.


### 3.3.1. “Chloride Concentration map” Model

The “Chloride concentration map” model generates a classified chloride concentration shapefile from a point feature table in text format by using Inverse Distance Weighting (IDW) interpolation technique (other interpolation techniques could be implemented in this tool). It can also import chloride concentration fields from Visual MODFLOW files. The dialog box shown in Figure4 is opened by double-clicking the model tool on the ArcToolbox window.Figure 4
Dialog box of the “Chloride concentration map” model.The model requires a polygon shapefile of the aquifer and the workspace containing a text file for each date to be analyzed. The text files have to includeX and Y UTM coordinates that define the locations of the point features and chloride concentration values (mg/l) in the observation wells. Text files also have to include column header, as shown in Figure 4. The text filename cannot exceed eight characters nor include blank spaces or special characters (underscore can be used as a substitute). It is not necessary that all points contain data every date of the period to be analyzed.The user has to indicate the fields (columns) in the input table that contain theX and Y coordinates and the chloride concentration value for each point (Figure 5). Optional settings concerning IDW interpolation techniques can be changed by the user. A reclassification of the values after interpolation is also required. Finally, the user has to choose a folder where the output shapefiles will be saved.Figure 5
Setting options defined by the user.When all the required parameters are filled out, the model can be executed by clicking “OK” at the bottom of the dialog box. The execution screen (Figure6) shows the running processes and it can be closed when the execution has successfully completed.Figure 6
Execution screen of the model.The “Chloride concentration map” model provides the following shapefiles for each date analyzed: (1) point shapefile of chloride concentration data; (2) raster from data interpolation; (3) polygon shapefile from interpolation covering the default extension; and (4) polygon shapefile from interpolation clipped to the shape of the aquifer.


### 3.3.2. “Hydraulic Head Map” Model

The “Hydraulic head map” model generates a classified hydraulic head shapefile from a point feature table in text format. It also generates a shapefile containing aquifer variables (chloride concentration and hydraulic head values) and aquifer parameters (storage coefficient and bottom of the aquifer). The dialog box is shown in Figure7.Figure 7
Dialog box of the “Hydraulic head map” model.This model has the same input requirements as for the “Chloride concentration map” model but focuses on hydraulic head data (m.a.s.l.). It also allows hydraulic head fields to be imported from Visual MODFLOW. The name of hydraulic head text files must be the same as the chloride concentration text files for each time period analyzed.The user has to indicate the location of the chloride shapefiles generated in the previous model (“Chloride concentration map” model). It also requires polygon shapefiles of storage coefficient and bottom (m) of the aquifer as inputs.The model generates shapefiles of hydraulic head data in an analogous way to the “Chloride concentration map” model. Moreover, it provides a polygon shapefile containing variables for each date analyzed (chloride concentration and hydraulic head values) and parameters (bottom and storage coefficient) of the aquifer. This shapefile is named “union_%name of the hydraulic head text file%_hh.shp”, where “% name of the hydraulic head text file%” is variable if different dates are analyzed.


### 3.3.3. “Summarizing SWI” Model

For the “Summarizing SWI” model, the methodology proposed in Baena-Ruiz et al. [13,15] and described in Section 3.1 has been implemented in the ArcGis environment. This tool generates Excel® tables containing statistics that summarize SWI at aquifer scale. It also generates conceptual cross sections (.shp), where the mean affected and nonaffected volumes are drawn for the aquifer (average values over a time period or instantaneous values on a specific date). If different dates are analyzed, it shows graphs representing the temporal evolution of Pa and Ta variables, percentage of affected volume, chloride concentration within the aquifer, and Ma index (or lumped vulnerability index). The dialog box for global status assessment is shown in Figure 8.Figure 8
Dialog box of the “Summarizing SWI” model.The “Summarizing SWI” model requires the folder path where the results of the “Hydraulic head model” have been previously saved to be specified. In this folder, the shapefile named “union_%name of the hydraulic head text file%_hh.shp” contains chloride concentration, hydraulic head, bottom of the aquifer, and storage coefficient fields. The user has to select from the pull-down list the corresponding column in the input shapefile for each field, as shown in Figure9.Figure 9
Selection of chloride field from the input table.The next required parameter in this tool is the “Chloride threshold.” It is defined as the chloride concentration value above which the aquifer is considered to be affected by SWI. This threshold may be set as the natural background level of the aquifer or as the relevant environmental quality standards. Hinsby et al. [14] proposed a method to calculate groundwater thresholds values based on these criteria.The shoreline length (m) is also required for subsequent calculations.X and Y axes establish the coordinate system of the conceptual cross sections. The GIS-SWIAS tool provides polyline shapefiles for X and Y axes located at (0,0), but the user can translate them to another coordinate origin or create new ones.The vertical scale factor is used to rescale the vertical magnitude (Ta) of the conceptual cross section if the factor Ta/Pa is too small. If vertical scale factor = 1, the conceptual cross section will maintain the real size ratio.Finally, two paths where the output results will be saved are required. “Output workspace statistic” will contain lumped variables reports in Excel table format for each date analyzed (Figure10) and mean statistics for the entire period. Four graphs will be also saved in this path: (1) temporal evolution of Pa and Ta variables; (2) percentage volume affected; (3) Ma index; and (4) chloride concentration within the aquifer (mean chloride concentration in the aquifer, mean chloride concentration in the affected volume, and the increment of concentration within the affected volume above the threshold).Figure 10
Summary statistics for a date.“Output workspace results” will contain the polygon shapefiles that allow the (1) mean affected and (2) nonaffected conceptual cross section within the aquifer for each date analyzed to be drawn, (3) the mean affected and (4) nonaffected conceptual cross section within a time period, and the (5) maximum affected cross section for a time period. These two paths can be the same for all results, but they have to be different from the output paths of the previous models.Figure11 and Table 1 show the graphical and statistical summaries, respectively, from the GIS-SWIAS tool.Figure 11
Graphical output from the GIS-SWIAS tool.Table 1
Lumped variables output (Excel® table format) from the GIS-SWIAS tool.
Lumped variables (Excel table)At a specific moment in timeStatistics for a time periodTotal aquifer volumeAverage aquifer volumeTotal aquifer affected volumeAverage aquifer affected volumeTotal chloride concentration∗aquifer volumeAverage chloride concentration∗aquifer volumeTotal chloride concentration∗aquifer volume in the affected volumeAverage chloride concentration∗aquifer volume in the affected volumeTotal aquifer areaAverage aquifer areaTotal aquifer affected areaAverage aquifer affected areaShoreline lengthShoreline lengthAverage aquifer thicknessAverage aquifer thickness within the periodAverage affected aquifer thicknessAverage affected aquifer thickness within the periodAverage orthogonal distance of the aquifer to the shorelineAverage orthogonal distance of the aquifer to the shoreline within the periodAverage orthogonal distance of the affected aquifer to the shorelineAverage orthogonal distance of the affected aquifer to the shoreline within the periodAverage chloride concentration within the aquiferAverage chloride concentration within the aquifer within the periodAverage chloride concentration within the affected volumeAverage chloride concentration within the affected volume within the periodIncrement of chloride concentration above the threshold valueIncrement of chloride concentration above the threshold value within the periodMa indexAverage Ma index within the periodPercentage of affected volumeAverage percentage of affected volume within the periodThe GIS-SWIAS tool can provide results for each date where information is available; these are obtained by iterative application of the described method. GIS-SWIAS allows historical [13] and future periods [15] to be analyzed if the hydraulic head and chloride maps come from a density-dependent flow model. By this means, GIS-SWIAS can be used to analyze adaptation strategies [15] in terms of reducing SWI, taking into account future potential scenarios that might include CC and/or GC, also considering projected land use change scenarios (new urbanized areas, crop transformations) [15].Moreover, this tool may be also used to summarize SWI vulnerability, for any index method applied to assess it. In this case, instead of the chloride concentration maps, generated by executing the “Chloride concentration map” model, polygon shapefiles of the vulnerability index (previously prepared by the user) would be used as inputs of the model “Summarizing SWI vulnerability” (Figure12), which also will require the “Hydraulic head map” shapefile generated by the tool, as previously described.Figure 12
Dialog box of the “Summarizing SWI vulnerability” model.The vulnerability index maps must be included as numerical fields (values obtained before defining the vulnerability classes). In order to generate the conceptual cross sections that summarize the “affected” aquifer volume, i.e., where the vulnerability to SWI is identified, the tool requires a vulnerability threshold to be input that represents the reference value chosen to distinguish between affected and nonaffected volumes. This threshold will be also used to assess the lumped vulnerability index.Just as in the definition of the Ma index, the lumped global value of vulnerability in the aquifer on a specific date is obtained by weighting the vulnerability score in each cell with its water storage. This lumped index also allows an analysis of the dynamic of SWI vulnerability of the system at aquifer scale to be performed. The lumped index can be also obtained using different threshold values [13, 15].


## 3.3.1. “Chloride Concentration map” Model

The “Chloride concentration map” model generates a classified chloride concentration shapefile from a point feature table in text format by using Inverse Distance Weighting (IDW) interpolation technique (other interpolation techniques could be implemented in this tool). It can also import chloride concentration fields from Visual MODFLOW files. The dialog box shown in Figure4 is opened by double-clicking the model tool on the ArcToolbox window.Figure 4
Dialog box of the “Chloride concentration map” model.The model requires a polygon shapefile of the aquifer and the workspace containing a text file for each date to be analyzed. The text files have to includeX and Y UTM coordinates that define the locations of the point features and chloride concentration values (mg/l) in the observation wells. Text files also have to include column header, as shown in Figure 4. The text filename cannot exceed eight characters nor include blank spaces or special characters (underscore can be used as a substitute). It is not necessary that all points contain data every date of the period to be analyzed.The user has to indicate the fields (columns) in the input table that contain theX and Y coordinates and the chloride concentration value for each point (Figure 5). Optional settings concerning IDW interpolation techniques can be changed by the user. A reclassification of the values after interpolation is also required. Finally, the user has to choose a folder where the output shapefiles will be saved.Figure 5
Setting options defined by the user.When all the required parameters are filled out, the model can be executed by clicking “OK” at the bottom of the dialog box. The execution screen (Figure6) shows the running processes and it can be closed when the execution has successfully completed.Figure 6
Execution screen of the model.The “Chloride concentration map” model provides the following shapefiles for each date analyzed: (1) point shapefile of chloride concentration data; (2) raster from data interpolation; (3) polygon shapefile from interpolation covering the default extension; and (4) polygon shapefile from interpolation clipped to the shape of the aquifer.


## 3.3.2. “Hydraulic Head Map” Model

The “Hydraulic head map” model generates a classified hydraulic head shapefile from a point feature table in text format. It also generates a shapefile containing aquifer variables (chloride concentration and hydraulic head values) and aquifer parameters (storage coefficient and bottom of the aquifer). The dialog box is shown in Figure7.Figure 7
Dialog box of the “Hydraulic head map” model.This model has the same input requirements as for the “Chloride concentration map” model but focuses on hydraulic head data (m.a.s.l.). It also allows hydraulic head fields to be imported from Visual MODFLOW. The name of hydraulic head text files must be the same as the chloride concentration text files for each time period analyzed.The user has to indicate the location of the chloride shapefiles generated in the previous model (“Chloride concentration map” model). It also requires polygon shapefiles of storage coefficient and bottom (m) of the aquifer as inputs.The model generates shapefiles of hydraulic head data in an analogous way to the “Chloride concentration map” model. Moreover, it provides a polygon shapefile containing variables for each date analyzed (chloride concentration and hydraulic head values) and parameters (bottom and storage coefficient) of the aquifer. This shapefile is named “union_%name of the hydraulic head text file%_hh.shp”, where “% name of the hydraulic head text file%” is variable if different dates are analyzed.


## 3.3.3. “Summarizing SWI” Model

For the “Summarizing SWI” model, the methodology proposed in Baena-Ruiz et al. [13,15] and described in Section 3.1 has been implemented in the ArcGis environment. This tool generates Excel® tables containing statistics that summarize SWI at aquifer scale. It also generates conceptual cross sections (.shp), where the mean affected and nonaffected volumes are drawn for the aquifer (average values over a time period or instantaneous values on a specific date). If different dates are analyzed, it shows graphs representing the temporal evolution of Pa and Ta variables, percentage of affected volume, chloride concentration within the aquifer, and Ma index (or lumped vulnerability index). The dialog box for global status assessment is shown in Figure 8.Figure 8
Dialog box of the “Summarizing SWI” model.The “Summarizing SWI” model requires the folder path where the results of the “Hydraulic head model” have been previously saved to be specified. In this folder, the shapefile named “union_%name of the hydraulic head text file%_hh.shp” contains chloride concentration, hydraulic head, bottom of the aquifer, and storage coefficient fields. The user has to select from the pull-down list the corresponding column in the input shapefile for each field, as shown in Figure9.Figure 9
Selection of chloride field from the input table.The next required parameter in this tool is the “Chloride threshold.” It is defined as the chloride concentration value above which the aquifer is considered to be affected by SWI. This threshold may be set as the natural background level of the aquifer or as the relevant environmental quality standards. Hinsby et al. [14] proposed a method to calculate groundwater thresholds values based on these criteria.The shoreline length (m) is also required for subsequent calculations.X and Y axes establish the coordinate system of the conceptual cross sections. The GIS-SWIAS tool provides polyline shapefiles for X and Y axes located at (0,0), but the user can translate them to another coordinate origin or create new ones.The vertical scale factor is used to rescale the vertical magnitude (Ta) of the conceptual cross section if the factor Ta/Pa is too small. If vertical scale factor = 1, the conceptual cross section will maintain the real size ratio.Finally, two paths where the output results will be saved are required. “Output workspace statistic” will contain lumped variables reports in Excel table format for each date analyzed (Figure10) and mean statistics for the entire period. Four graphs will be also saved in this path: (1) temporal evolution of Pa and Ta variables; (2) percentage volume affected; (3) Ma index; and (4) chloride concentration within the aquifer (mean chloride concentration in the aquifer, mean chloride concentration in the affected volume, and the increment of concentration within the affected volume above the threshold).Figure 10
Summary statistics for a date.“Output workspace results” will contain the polygon shapefiles that allow the (1) mean affected and (2) nonaffected conceptual cross section within the aquifer for each date analyzed to be drawn, (3) the mean affected and (4) nonaffected conceptual cross section within a time period, and the (5) maximum affected cross section for a time period. These two paths can be the same for all results, but they have to be different from the output paths of the previous models.Figure11 and Table 1 show the graphical and statistical summaries, respectively, from the GIS-SWIAS tool.Figure 11
Graphical output from the GIS-SWIAS tool.Table 1
Lumped variables output (Excel® table format) from the GIS-SWIAS tool.
Lumped variables (Excel table)At a specific moment in timeStatistics for a time periodTotal aquifer volumeAverage aquifer volumeTotal aquifer affected volumeAverage aquifer affected volumeTotal chloride concentration∗aquifer volumeAverage chloride concentration∗aquifer volumeTotal chloride concentration∗aquifer volume in the affected volumeAverage chloride concentration∗aquifer volume in the affected volumeTotal aquifer areaAverage aquifer areaTotal aquifer affected areaAverage aquifer affected areaShoreline lengthShoreline lengthAverage aquifer thicknessAverage aquifer thickness within the periodAverage affected aquifer thicknessAverage affected aquifer thickness within the periodAverage orthogonal distance of the aquifer to the shorelineAverage orthogonal distance of the aquifer to the shoreline within the periodAverage orthogonal distance of the affected aquifer to the shorelineAverage orthogonal distance of the affected aquifer to the shoreline within the periodAverage chloride concentration within the aquiferAverage chloride concentration within the aquifer within the periodAverage chloride concentration within the affected volumeAverage chloride concentration within the affected volume within the periodIncrement of chloride concentration above the threshold valueIncrement of chloride concentration above the threshold value within the periodMa indexAverage Ma index within the periodPercentage of affected volumeAverage percentage of affected volume within the periodThe GIS-SWIAS tool can provide results for each date where information is available; these are obtained by iterative application of the described method. GIS-SWIAS allows historical [13] and future periods [15] to be analyzed if the hydraulic head and chloride maps come from a density-dependent flow model. By this means, GIS-SWIAS can be used to analyze adaptation strategies [15] in terms of reducing SWI, taking into account future potential scenarios that might include CC and/or GC, also considering projected land use change scenarios (new urbanized areas, crop transformations) [15].Moreover, this tool may be also used to summarize SWI vulnerability, for any index method applied to assess it. In this case, instead of the chloride concentration maps, generated by executing the “Chloride concentration map” model, polygon shapefiles of the vulnerability index (previously prepared by the user) would be used as inputs of the model “Summarizing SWI vulnerability” (Figure12), which also will require the “Hydraulic head map” shapefile generated by the tool, as previously described.Figure 12
Dialog box of the “Summarizing SWI vulnerability” model.The vulnerability index maps must be included as numerical fields (values obtained before defining the vulnerability classes). In order to generate the conceptual cross sections that summarize the “affected” aquifer volume, i.e., where the vulnerability to SWI is identified, the tool requires a vulnerability threshold to be input that represents the reference value chosen to distinguish between affected and nonaffected volumes. This threshold will be also used to assess the lumped vulnerability index.Just as in the definition of the Ma index, the lumped global value of vulnerability in the aquifer on a specific date is obtained by weighting the vulnerability score in each cell with its water storage. This lumped index also allows an analysis of the dynamic of SWI vulnerability of the system at aquifer scale to be performed. The lumped index can be also obtained using different threshold values [13, 15].


## 4. Discussion

GIS-SWIAS is a user-friendly polyvalent ArcGIS tool that provides a comprehensive overview of SWI status and vulnerability at aquifer scale. It integrates three models, which are documented in order to briefly explain the tool’s description, its utility, and the data required for each item. This tool can be applied by scientists and decision makers, who may not be advanced users of GIS, to summarize SWI problems. Many GIS-based tools have demonstrated to be powerful and cost-effective to analyze groundwater issues (Criollo et al. and Perdikaki et al. [19, 32]). Moreover, GIS models as ModelBuilder models can be integrated into other platforms by using the Python script tool (Menezes and Inyang [33]).Due to the heterogeneous distribution of seawater intrusion,distributed information and assessments are required to study its impacts [8, 30]. For this reason, the methods for modelling [34, 35] SWI impacts and the user-friendly tools developed based on them [36–38] also require distributed inputs and calculations. The GIS-SWIAS is a tool that could be classified as a postprocessing tool to summarize and help in the analysis of SWI impacts at aquifer scale. This tool produces both distributed and lumped results at aquifer scale, but, logically, it also requires distributed inputs and assessments, as described in the previous sections. In this group of postprocessing tools, we find in the literature, for example, [39]. GIS-SWIAS is a new tool, in which the method proposed by [13, 15] has been implemented. A significant novelty of this method with respect to other previously developed methods is that the proposed lumped index to summarize SWI status at aquifer scale is based on a variable with physical meaning (mass of chloride that causes the concentration in some areas to exceed the natural threshold). On the other hand, a novel aspect of this tool is that, from the distributed information and calculations, GIS-SWIAS allows easy computation of the affected volume containing a chloride concentration above a threshold. This tool also helps to produce lumped SWI outputs (indices) at aquifer scale. It produces valuable information that helps draw conclusions about the dynamic at aquifer scale, in terms of affected volume and global SWI intensity. Thus, it also provides insight into aquifer resilience and trends. Therefore, it will help to identify coastal groundwater bodies that require new management strategies to be implemented to achieve a good status.The identification of SWI (the phenomenon that we want to analyze) requires a threshold value established that defines the inflection point beyond which the aquifer begins to register an impact. Previous research shows that the impact of SWI is significantly sensitive to the choice of the threshold value adopted [13]. The significant uncertainties around determining these threshold values [14] and the sensitivity of whether the aquifer is reported as being impacted by SWI or not increase the practical interest of the GIS-SWIAS tool: it is capable of performing the extensive calculations required to summarize SWI at aquifer scale, for the analyses of both historical and potential scenarios, considering different threshold values, which allow the comparison of the results.With respect to the maps employed as inputs, the tool allows two options: to generate maps from available data using different interpolation techniques integrated in the tool and to take the maps directly from SEAWAT files. This functionality—which allows maps to be generated from point data or to be loaded from other commonly employed tools—has also been implemented in other SWI assessment tools [36, 37]. However, as far as we know, it is not available in postprocessing tools. In cases where map inputs are taken from density-dependent models, a comparative assessment of different scenarios (climatic conditions and/ or management strategies) could be performed. The physical-process approach can be applied to propagate and compare various potential conditions, and so in this case maps can be obtained and compared for different scenarios (e.g., management scenarios or future potential CC scenarios); this means that the output of the tool can support the decision-making process [15]. In contrast, when the maps employed to define the indices are obtained by applying simple interpolation approaches, analysis is limited to the historical period for which the data are available.The tool also helps to analyze the vulnerability to seawater intrusion at aquifer scale. In the literature, we find different methods to assess groundwater vulnerability depending on the drivers of pollution (Aller et al., Vias et al., and Baena-Ruiz and Pulido-Velazquez [40–42]), pumping (Pulido-Velazquez et al. [43, 44]), and SWI [12, 16]. User-friendly tools have appeared to assist in this assessment; some of them were developed in a GIS environment [45]. Nevertheless, there are no tools that help in the assessment of SWI vulnerability with that focus on postprocessing. Therefore, this is the first postprocessing tool described that integrates SWI status and vulnerability assessment, which is very valuable information to help identify the significance of SWI problems in aquifers and potential sustainable solutions.The GIS-SWIAS tool has been applied to two different case studies in the Mediterranean area of Spain (Plana de Oropesa Torreblanca and Plana de Vinaroz), obtaining the results described in previous papers [13, 15]. In [13], the process automation to generate the interpolated maps enabled authors to analyze SWI status and vulnerability over an extended time period (1977–2015) and to prove the sensitivity of results to the chloride threshold value (two threshold values were analyzed: 250 mg/l and 1100 mg/l) in Plana de Oropesa-Torreblanca and Plana de Vinaroz aquifers. In [15], the impacts of future GC scenarios were analyzed in Plana de Oropesa-Torreblanca aquifer. The methodology from [13] was adapted to compare six potential future scenarios including adaptation strategies. The historical period came from 1973 to 2010 and the six future scenarios covered the period 2011–2035.The underlying methodology implemented in GIS-SWIAS was applied in [13] by interpolating chloride maps and hydraulic head from observation points, whereas the information to generate the field maps in [15] was loaded from SEAWAT model. The results of these studies for the Plana de Oropesa-Torreblanca aquifer show differences that reveal that the physical-process SWI approximations obtained using the density-dependent flow model give a more accurate representation. Despite these differences, the results are in the same order of magnitude. Other authors who have developed indices related to SWI [6, 9] have also proved that results do not differ considerably by including three-dimensional salinity data. Furthermore, the approximation obtained using interpolation will depend on the number of observations and the distribution of these points within the aquifer.Although this tool has been developed to analyze SWI problem, it could be applied to study the lumped impact of any contaminant on groundwater and/or the groundwater vulnerability by applying any vulnerability index. In this case, instead of the “Shoreline length” parameter to generate the cross section, other equivalent lengths (e.g., the aquifer length orthogonal to the groundwater flow direction) should be considered. Therefore, GIS-SWIAS fulfils the requirements of flexibility, sturdiness, easy interaction, and user-friendliness, which make it a useful tool in the decision-making process. It will allow using them as “share vision models/tools” to help in the discussion of management alternatives between stakeholders and administration representatives [46]. Many Decision Support Systems tools were not successful because they were not user-friendly [47, 48].


### 4.1. Assumptions and Limitations

In this section we summarize the main assumptions/limitations of the GIS-SWIAS tool and in the implemented methodology.


#### 4.1.1. Underlying Methodology

(i)
The methodology implemented [13] provides lumped results at aquifer scale. An analysis based exclusively on these lumped results would lose information about the heterogeneity of SWI in the system.(ii)
Although the method can be applied in cases where a distributed flow model is not available, a proper assessment for a specific date requires sufficient distributed data about hydraulic head, chloride concentration, and aquifer bottom. On the other hand, to generate the required fields from the data, it needs to apply an interpolation approach whose parameters influence on the results. It is recommended that the user tests different interpolation approaches and parameters in order to obtain results in agreement with expert opinion.(iii)
In order to assess the dynamic of SWI and the resilience of the aquifer, there is need to have data for several dates over a long period.(iv)
If the aim is to assess sustainable management strategies, appropriate hydrological models to propagate the impacts are required, not only of the climatic conditions, but also of the management strategies.


#### 4.1.2. GIS-SWIAS Tool

(i)
This tool has been implemented in ArcGIS and is licensed GIS software. Future implementations may be developed in an open environment.(ii)
The GIS-SWIAS tool had to be divided into three models due to some difficulties that arose to implement all the geoprocesses into one model. Nevertheless, in the future, we plan to combine these three steps into one, improving the usability of the model.(iii)
Graphic results are based on a predefined template. Although graphs in ArcGIS are customizable, there is not so much potential.(iv)
Output reports are generated automatically in Excel format and lack a refined style.


## 4.1. Assumptions and Limitations

In this section we summarize the main assumptions/limitations of the GIS-SWIAS tool and in the implemented methodology.


### 4.1.1. Underlying Methodology

(i)
The methodology implemented [13] provides lumped results at aquifer scale. An analysis based exclusively on these lumped results would lose information about the heterogeneity of SWI in the system.(ii)
Although the method can be applied in cases where a distributed flow model is not available, a proper assessment for a specific date requires sufficient distributed data about hydraulic head, chloride concentration, and aquifer bottom. On the other hand, to generate the required fields from the data, it needs to apply an interpolation approach whose parameters influence on the results. It is recommended that the user tests different interpolation approaches and parameters in order to obtain results in agreement with expert opinion.(iii)
In order to assess the dynamic of SWI and the resilience of the aquifer, there is need to have data for several dates over a long period.(iv)
If the aim is to assess sustainable management strategies, appropriate hydrological models to propagate the impacts are required, not only of the climatic conditions, but also of the management strategies.


### 4.1.2. GIS-SWIAS Tool

(i)
This tool has been implemented in ArcGIS and is licensed GIS software. Future implementations may be developed in an open environment.(ii)
The GIS-SWIAS tool had to be divided into three models due to some difficulties that arose to implement all the geoprocesses into one model. Nevertheless, in the future, we plan to combine these three steps into one, improving the usability of the model.(iii)
Graphic results are based on a predefined template. Although graphs in ArcGIS are customizable, there is not so much potential.(iv)
Output reports are generated automatically in Excel format and lack a refined style.


## 4.1.1. Underlying Methodology

(i)
The methodology implemented [13] provides lumped results at aquifer scale. An analysis based exclusively on these lumped results would lose information about the heterogeneity of SWI in the system.(ii)
Although the method can be applied in cases where a distributed flow model is not available, a proper assessment for a specific date requires sufficient distributed data about hydraulic head, chloride concentration, and aquifer bottom. On the other hand, to generate the required fields from the data, it needs to apply an interpolation approach whose parameters influence on the results. It is recommended that the user tests different interpolation approaches and parameters in order to obtain results in agreement with expert opinion.(iii)
In order to assess the dynamic of SWI and the resilience of the aquifer, there is need to have data for several dates over a long period.(iv)
If the aim is to assess sustainable management strategies, appropriate hydrological models to propagate the impacts are required, not only of the climatic conditions, but also of the management strategies.


## 4.1.2. GIS-SWIAS Tool

(i)
This tool has been implemented in ArcGIS and is licensed GIS software. Future implementations may be developed in an open environment.(ii)
The GIS-SWIAS tool had to be divided into three models due to some difficulties that arose to implement all the geoprocesses into one model. Nevertheless, in the future, we plan to combine these three steps into one, improving the usability of the model.(iii)
Graphic results are based on a predefined template. Although graphs in ArcGIS are customizable, there is not so much potential.(iv)
Output reports are generated automatically in Excel format and lack a refined style.


## 5. Conclusions

In this paper we describe a new general tool, GIS-SWIAS. It is an ArcGIS-based tool, designed to analyze SWI status and vulnerability at aquifer scale by applying the method presented by [13, 15]. It is a user-friendly tool that allows georeferenced information to be dealt with, and it is easy to introduce the required data (inputs) and to efficiently perform the demanding computational operations required. Its outputs are in the form of reports and images to summarize the magnitude, intensity, and temporal evolution of SWI within an aquifer.The GIS-SWIAS tool can be applied to assess historical SWI dynamic in case studies where we do not have a previous model. Nevertheless, if we want to analyze in a rational quantitative analysis of the various alternative management scenarios to manage SWI in a sustainable manner, the GIS-SWIAS tool will need to take information on hydraulic head and chloride concentration distribution generated from simulations of their impacts by a calibrated density-dependent flow model. In such cases, adaptation strategies to potential future scenarios, whose distributed impacts have to be propagated within the previously calibrated models, could usefully be analyzed and compared using this tool. GIS-SWIAS can be applied to assess not only SWI status at aquifer scale, but also vulnerability to any contaminant.Given all these ways that the GIS-SWIAS tool can be applied, it provides a valuable tool for both researcher and technician to assess SWI dynamics and aquifer resilience under different management scenarios. It can support the decision-making process in the rational selection of sustainable management strategies. The tool’s performance has been tested and confirmed in two case studies described in previous research works.It can be applied to any case study. The easy-to-use workflow and the few input data required facilitate its application to a large number of case studies in order to compare SWI.



---

*Source: 8818634-2021-05-12.xml*