Mapping of Salinity Ingress Using Galdit Model for Sirkali Coastal Region : A Case Study

A model is traced to evaluate and enumerate the significance of vulnerability to seawater intrusion due to excessive ground water withdrawals and some anthropogenic activities at coastal aquifers. So taking these issues into account few thematic maps which were influencing the saline water intrusion were prepared and overlaid using Geographical Information Systems (GIS). Based on GALDIT method, the groundwater vulnerability cartography has been assessed. To reckon the GALDIT index it requires six parameters like aquifer type, aquifer hydraulic conductivity, depth to groundwater level (AMSL), distance from the shore, impact of existing status of seawater intrusion and thickness of the aquifer. This GALDIT is the indicator scores and summing them and dividing by the total weight for determining the relative role of each one. Apart from this an identification of saltwater intruded area is done by using indicators of saltwater intrusion like Cl/(HCO3 + CO3) ratio and Na/Cl ratio. The vulnerability areas are classified as moderate with an area of 147.31 sq. km and low covering an area of 168.72 sq. km respectively based on the thematic maps. The final thematic map can be used for management of the coastal ground water resources.


Introduction
During the last decades, the coastal aquifer systems are under to several pressures, due to urbanization and intensive agriculture.Approximately 70% of the population of the earth life is in coastal areas and the majority of these people depend on coastal aquifers for freshwater.The mis-management of water resources in an area causes negative effects including depletion of the aquifer storage, groundwater level decline, and seawater intrusion in coastal areas, land subsidence, quality deterioration and environmental problems in other water bodies.As this coastal aquifers constitute an important source of fresh water supply but are often confronted with the problem of seawater intrusion.In coastal plains, due to inadequate storage facilities, most of the rainwater flows towards sea as runoff.Population growth, agricultural, industrial and domestic requirements exploit the available ground water with decreasing recharge areas.Regulation of ground water extraction is essential to optimize the overdraft, reducing the risk of seawater intrusion.The entire seawater intrusion phenomenon is governed by Ghyben-Herzberg relation [1].The lighter fresh water lies over the seawater and the boundary surface between them is known as the freshwater-seawater interface.This distribution was attributed to a hydrostatic equilibrium existing between the two fluids of different densities.For each meter of fresh ground water found above sea level, 40 m of freshwater exists below sea level at that point.When drawdown occurs, the base of the freshwater lens is adjusted at a rate of 40 m for each meter of drawdown in the well through upcoming effect.
Rainfall plays a major role in recharge of an oceanic island.When this precipitation is sufficient it gives a sufficient need of fresh water and this freshwater exists as a lens above seawater due to the variation in there densities.Freshwater lenses have been described as some of the most vulnerable aquifer systems in the world, with seawater intrusion (SWI) (i.e., the encroachment of seawater into fresh coastal aquifers) being a critical management issue [2].This interface layer is weaken by various factors like low topography of oceanic island, less freshwater availability and heavy pumping, heavy reliance by local communities with few alternative water supplies [3] and [2].
Generally, the term of vulnerability refers to the potential degree of harm that can be expected depending on the characteristics of an element at risk with respect to a certain hazard [4] [5].Relating to groundwater, the vulnerability is defined by [6] as "the sensitivity of groundwater quality to an imposed contaminant load, which is determined by the intrinsic characteristics of the aquifer".Thus the vulnerability of groundwater to different pollutants or to seawater intrusion constitutes a subject of analysis in several studies [7]- [9].Also the vulnerability of soil to salinization is demonstrated in many studies [10]- [14].
In this study an attempt is made to outline the area of scope of saline water intrusion by preparation of vulnerability map in the Sirkali coastal region.This study will be helpful for the groundwater management and plans to preserve the freshwater.

Study Area
The study area lies between longitudes: 79˚05'00" to 11˚15'30"N with an area of 316.2 sq.km (Figure 1).Sirkazhi is located in the Uppanar river basin with an area of 168.8 sq.km and Mannampandal in Cauvery river basin with an area of 147.4 sq.km.Topographically the area is flat, except beach ridges tending towards the Bay of Bengal.The elevation of the Uppanar River Basin which is the upper part of the study area varies from 0 to 14.33 m above mean sea level (amsl) and in low areas was <3 m (amsl).In Uppanar river basin backwater, proceed to a distance of above 20 km towards inland.In Cauvery river basin which is the lower part of the study area the relief varies from 0.0 to 17.3 m (amsl) and the backwater inundation is within 5 km from the coast where the topography is around <3 m amsl.The annual average rainfall ranges from 1200 -1500 mm/yr with a major contribution from NE monsoon during October-December.Economically this area depends on agriculture, prawn cultivation and brick industries.
Geologically the area known as younger (quaternary) deposits [15]; the sediments classified as alluvial plain deposits (Cauvery formation) of the Cauvery River and its distributaries, fluviomarine deltaic plain deposits (Nagapattinam formation), marine coastal plain deposits (East Coast formation).Paleo-channels noticed with admixtures of sand, silt clay and gravels.The deltaic plain includes paleo tidal flats with clays and sand, sand ridges of grey brown sand.The marine coastal plains include tidal flats, deposits of sand, clay and tidal clays.The study area underlined by sedimentary formations and geomorphologically it is divided into the three major units (i.e., fluvial, fluviomarine, and marine units).Groundwater occurs under unconfined, semi-confined, and represented by sands, gravels, sandy clays and variegated clays, and its thickness ranges from 10 -35 m [16].Groundwater extracted through the shallow tube wells and dug-cum-bore wells.In the coastal area, fresh water is available only in beach ridges and in sand dunes.

Methodology
In this study we carry out the application of the GALDIT method to aquifer of Sirkali coastal region.This method acts to distinguish the vulnerability of the coastal aquifers, based on six parameters existing (Table 1).The parameters of the intrusion of sea water were described by [17].
A numerical grade system has been framed using the parameters to assess saline water intrusion extent.The system contains three significant parts: weights, ranges and ratings.Each parameter has been assessed with respect to the other parameter to determine the relation of each factor by assigning a relative weight.The other parameter using GALDIT method indicates an identification of saltwater intruded area using the ration like Cl/(HCO 3 + CO 3 ).Parameters that are normally used for seawater intrusion analysis are Ca/Mg, Total alkalinity/ Total hardness, Cl/(CO 3 + HCO 3 ) ratios, Chloride content, Electrical conductivity (EC) etc.In this study, the Cl/(CO 3 + HCO 3 ) ratio is used to delineate the interface (Table 2).Chloride is most dominant in ocean water and normally occurs in small amounts in ground water, while HCO 3 is usually the most abundant negative ion in ground water but it occurs in minor amounts in the seawater.The Cl/(CO 3 + HCO 3 ) ratio is calculated for all the wells and these point values are taken to prepare the spatial variation map using the kriging and IDW interpolation technique.Similarly, the other factors are evaluated and used.From these, the GALDIT (Combining the first letters of the above significant factors) Index is calculated as: where, W 1 to W 6 are the relative weights assigned to the six factors.The ratings of for all the factors were assigned by reclassifying the influencing factors with different ranges through a trial and error method.In this study, the ranges and weight ages were modified to suit to the study area [18] [19].

Mapping of GALDIT Index
Mapping of GALDIT index gives its variation in the study area, which is done by Arc GIS 10.3 software.For mapping of this GALDIT index it requires the location of all samplimg sites and there parameters data.The important steps involved for mappings are georeferencing, digitization and spatial interpolation.Georeferencing defines the location of a dataset using known map coordinates and assigns it a coordinate system.This allows for the dataset to be viewed, queried, and analyzed with other geographic data.Digitzing is the process of making features which can be editable and these features have an additional spatial and non-spatial attributes that can be assigned.By digitizing these features, you make them available for mapping once you have added the tabular data to the attribute table.Spatial interpretation is the process of using points with known values to estimate values at other points.Spatial interpretation is therefore a means of creating surface data from sample points so that the surface data can be used for analysis.

Groundwater Occurrence, G
The parameter G (aquifer type) affects the degree of advancement of the marine water into the groundwater.
From the pump test data and TDS contours it is evident that the aquifer is leaky and unconfined in nature with rich groundwater potential and hence a rating of 5.5 to 7.5 is adopted as per the specifications.Also, from the field observations it is evident that the aquifer is shallow and confined in nature.An unconfined aquifer, in natural conditions, is more affected by marine water intrusion than a confined one.In the study area, the aquifer is unconfined and corresponds to class 7.5 as in Figure 2.

Aquifer Hydraulic Conductivity, A
Hydraulic conductivity or the permeability is the aptitude of a soil or rock to let itself cross by water under the effect of a hydraulic gradient.A layer of containment is a geological unit of low or very low hydraulic conductivity (< to 10 −7 m/s) whereas the formations considered as aquifers.It consists of materials whose hydraulic conductivity exceeds 10 −4 m/s.Finally from all the studies conducted to evaluate the hydraulic conductivity in the study area, we can conclude that the hydraulic conductivity varies up to 5 m/sec and hence a common GALDIT rating of 2.5 can be assigned for the entire study area and as in Figure 3.

Height of Water above the Sea Level, L
The level of groundwater compared to the average altitude of the sea is a very significant factor in the evaluation of the sea water intrusion in any area.By this it determines the possibility of the water pressure to move back the sea front [7].In general the minimal values of groundwater level below the sea level remain most significant, because they provide the strongest possible vulnerability to this marine water intrusion.The groundwater levels were measured at the monitoring wells identified.The parameter required for the present study, "L" is obtained by reducing the water level with respect to mean sea level and is shown in Figure 4.

Distance from the Shore, D
The impact of the intrusion of sea water generally decreases when moving perpendicularly to the shore towards the interior.This parameter was estimated according to three distances (500 m, 750 m and 1000 m) perpendicular to the line of coast and the rivers of the uppanar and cauvery River.The maximum estimate of 10 is adopted for the distance lower than 500 m of the coast, whereas the minimal one (2.5) is allotted for all those higher than 1000 m.The values of 7.5 and 5 are given, respectively, with the distances from 500 to 750 m and from 750 to 1000 m.The distribution of the parameter D of GALDIT of the aquifer of coast is represented in Figure 5.

Impact of Existing Status of Saltwater Intrusion, I
The chloride concentration of groundwater defines the extent of saltwater intrusion.The chloride concentration in the range of 120 to 3260 mg/l is indicative of saltwater intrusion [20].The present results infer that the aquifer water is contaminated with saltwater as Cl − , the most abundant ion in saltwater, is in higher proportions.The 3 HCO − ion which is the most dominant ion in fresh groundwater occurs generally in small amounts in saltwater.Chloride concentration greater than 250 mg/l is considered unfit for drinking purpose.The existing imbalance in the seawater-freshwater interface should be considered while mapping the aquifer vulnerability to seawater intrusion [21].Chloride is the dominant ion in the seawater and it is only available in small quantities in groundwater while bicarbonate, which is available in large quantities in groundwater, occurs only in very small quantities in  CO − ] is a criterion to identify the extent of seawater intrusion into the coastal aquifers [22] and can be used if the chemical analysis data is available.The evaluation of this parameter was given starting from an ionic analysis of chromatography at the laboratory of the Centre for Water Resources Development and Management (CWRDM) Calicut, Kerala.The distribution of this parameter GALDIT is presented in Figure 6(a).

Thickness of the Aquifer, Z
The thickness of the aquifer is obtained from the electrical resistivity survey conducted in the study area at 39 locations Figure 6(b).The resistivity survey indicated that the area consists of shallow unconfined aquifer with the thickness ranging from 2.6 m to 75.2 m.From the resistivity survey we can observe that all the values are ranging from 2.6 m to 75.2 m and hence the GALDIT rating of 2.5 to 10 is adopted in the study area as shown in Figure 7.

Computing the GALDIT Index
The calculation of GALDIT Index and the superposition of the various layers of the parameters in a Geographical Information System (GIS) make it possible to identify the significant and susceptible areas that could be affected by a seawater intrusion.The map derived for this study area is shown in Figure 8.The figure clearly states that the area closer to the coast is in moderate vulnerable and away from the coast its under low vulnerable state and in the upper part of the study area it is observed that NNW part is in the moderate vulnerable condition that is due to the surface phenomena like aquaculture and the backwaters as seen in the Figure 1.This moderate vulnerability class is very well matched with the TDS of groundwater in the study area Figure 9.

Conclusions
Geologically the area known as younger (quaternary) deposits; the sediments classified as alluvial plain deposits    (Cauvery formation) of the Cauvery River and its distributaries, fluviomarine deltaic plain deposits (Nagapattinam formation), marine coastal plain deposits (East Coast formation).Paleo-channels noticed with admixtures of sand, silt clay and gravels.
The application of GALDIT method to groundwater of Nagapattinam coast made it possible to evaluate the impact of the increase of saline level.The groundwater is characterized by a low vulnerability in the lower part of area i.e.Cauvery basin and a moderate vulnerability in upper part i.e.Uppanar basin with a strong contamination of the marine intrusion in the coast zone and in proximity of Uppanar River which is seen fresh from Figure 9.According to GALDIT method the aquifer is classified as moderate with an area of 147.31 sq.km and low covering an area of 168.72 sq.km in the entire coastal area taken into consideration.

Figure 1 .
Figure 1.Key map of the study area showing different geomorphological units.

Figure 4 .
Figure 4. Parameter L (height of water above the sea level).

Figure 6 .
Figure 6.(a) Parameter I (impact of existing status of saltwater intrusion); (b) Vertical electrical sounding location map.

Figure 9 .
Figure 9.Total Dissolved Solids (TDS) values observed in Uppanar and Cauvery River Basin along cross sections.

Table 2 .
Impact of existing status of seawater intrusion.