Hydrochemical appraisal of groundwater quality for drinking and irrigation: a case study in parts of southwest coast of Tamil Nadu, India

The quality of groundwater standards in Muttom–Mandaikkadu coastal stretch is the focus of the present study, whose coastal aquifers are particularly at risk due to intrusion of marine water. Thirty groundwater samples were scrutinized for the assessment of physical and chemical parameters during January and June. Hydrochemical characteristics were spatially depicted to understand the spatial variations such as (pH, EC, TDS, Na+, K+, Ca2+ and Mg2+, SO42−, HCO3−, and Cl−). Drinking water quality index based on those 11 parameters and irrigation water quality index based on EC, Na%, sodium adsorption ratio and permeability index was used to assess the water quality for drinking and irrigation, respectively. These results demonstrate that dominant hydrochemical facies for groundwater in both months are Na-K-Cl-SO4 type and Ca-Mg-Cl-SO4 type. The USSL diagram endorses that most of the water samples belong to low-medium salinity with low sodium hazards. Cl−/HCO3− ratio indicates that the majority of the samples show low to moderate seawater intrusion in the study area. Additionally, six vertical electrical sounding measuring points (Schlumberger array) were carried out in order to determine the number of the underlying layers, aquifer depths and their thicknesses as well as its influence by the marine water. The geophysical self-potential measurements suggest that the groundwater in the Manavalakurichi area is prone to contamination by seawater intrusion, confirmed by the use of Schlumberger vertical electrical sounding. For better understanding, the subsurface layers were shown in a 2D model using the constructed geoelectrical cross section.


Introduction
In India, groundwater suites the drinking and domestic needs of nearly 90% of the rural and 30% of the urban population (Bhunia 2020). Usually, groundwater is considered a safer resource than surface water. Still, several factors such as uncontrolled discharge of industrial, agricultural and domestic effluents, improper land use, certain geological formations, shortage of rainfall, and low infiltration rate affect groundwater quality (Adimalla 2019). Increasing demand for water during these days makes the problem of water contamination a major concern and has also caused water shortage in different parts of the world. According to WHO reports, 65% of rural and 36% of urban inhabitants do not consume safe drinking water in India (Sexenal and Mishra 2011).
On the other hand, deprived quality of groundwater also affects agriculture as it constitutes 53% of the country's irrigation potential (Wilcox 1948;Hem 1991;FAO 2003;WHO 2004). If the groundwater in the aquifer is contaminated, it persists for several years and it is not easy to retrieve the quality (Jakeman et al. 2016). So, there is a need to monitor and assess the quality of groundwater. To understand the suitability for drinking as well as irrigation purposes, a hydrochemical assessment of the groundwater must be carried out (Nakhaei et al. 2016). Figuring of several 1 3 53 Page 2 of 20 physico-chemical criteria like pH, electrical conductivity (EC), total dissolved solids (TDS), Ca 2+ , Mg 2+ , K + , Na + , Cl − , HCO 3 − and SO 4 2− and their comparison with existing firm standards disclose the quality of water and its suitability for drinking and irrigation needs. Several hydrochemical appraisal studies were carried out in recent days worldwide as well as India to understand the nature of the groundwater (Srinivas et al. 2014;Tiwari et al. 2017;Adimalla et al. 2018;He et al. 2018;Li et al. 2018;Chidambaram et al. 2018;Kumari and Rai 2020;Shaikh et al. 2020). Golchin and Azhdary Moghaddam (2016) in their study named hydrogeochemical characteristics and groundwater quality assessment in Iranshahr plain aquifer, Iran, obtained results such that the groundwater salinity hazard is low to high, but the Na hazard is low to medium, and in regard to irrigation water, the quality is low to medium. Ghalib (2017) had carried out a study that focused on assessing the groundwater quality of the shallow aquifer in the northeastern Wasit Governorate, Iraq. The studied groundwater samples were found to be oversaturated with carbonate minerals and undersaturated with evaporated minerals. A comparison of groundwater quality concerning drinking water standards showed that most groundwater samples were unsuitable for drinking purposes. A study was carried out by Mostafa et al (2017) to understand the hydrogeochemistry and groundwater quality in the Rajshahi City of Bangladesh. This study demonstrated that the rock-water interaction was the major geochemical process controlling groundwater chemistry in the study area. The study results also revealed that the groundwater quality in Rajshahi City area was of great concern and not suitable for human consumption without adequate treatment. The research paper by Saraswat et al (2018) provides an exploratory investigation in the hydrogeochemical characteristics of groundwater and thereby assesses the suitability of groundwater as an alternative and reliable resource for public water supply in the Indian city of Surat. The results of this study outline the unsustainability of groundwater for direct consumption, especially without any improved on-site water treatment, but it is appropriate for irrigation purposes. Mahanta et al (2019) had authored a study that aims to approach the analysis of groundwater samples of the Western part of Odisha, i.e., Maneswar Block of Sambalpur District, to assess the quality of groundwater in terms of drinking and irrigation purposes.
Nowadays, geographical information system (GIS) has been widely used as an effective tool to integrate various databases and develop solutions for assessing water quality, reducing water resources problems, preventing flooding and understanding the natural environment (Venkatramanan et al. 2019). Precautionary indicators of potential environment-related health problems include groundwater quality distribution maps and risk assessment maps (Arnous et al. 2011;Gnanachandrasamy et al. 2015;Ameen 2019;Li et al 2019).
In the Manavalakurichi and enclosing areas of India's southern apex, a total of 30 groundwater samples have been acquired for the current study. This study's primary focus is to delineate the level of spoilation in groundwater in the Muttom-Mandaikkadu coastal stretch. There is a possibility of groundwater deterioration due to saline water intrusion as the study area's southern boundary is surrounded by sea. So, the present study finds its importance in delineating the pollution in groundwater using water quality analysis.

Study area
The study area Muttom to Mandaikkadu coastal zone located within the Kanyakumari district in the southern part of Tamil Nadu. Figure 1 shows the study area map with geology. The climate is humid with extremely hot summers; the southern part of the study area is bounded by the Indian Ocean and in the west by the Arabian Sea. The area being part of the west coast enjoys a typical monsoonal climate and receives an average annual rainfall of about 1465 mm. The maximum rain is received during April, May, June, and July (southwest monsoon) and in October, November (northeast monsoon). Rainfall forms the only source of aquifer replenishment of these aquifers. Physiographically, the region of study is almost outstretched with a gentle slope, generally facing south. Valliyar is the only river leading in the study area, and the drainage design is sub dendritic, a seasonal river. Manavalakurichi is one of the familiar coastal sectors for its mineral deposits. Several mining zones under the authority of Indian Rare Earth Limited (IREL) are present here (Srinivas et al. 2015). The typical temperature within the study area ranges between 22.8 °C and 33.6 °C. From the geological perspective, most of the part in the study area underlies Quaternary fluviomarine deposits. The Western Ghats Mountain occupies the entire north and northwestern region of the district with a maximum elevation of 1658 m. Sand dunes (Teri sands) and marshy swamps typically cover the coastal tract of the area. The dominant soil types of the study area are coastal sand, red soil, brown soil, and red lateritic soil.
The area is underlain by an Archaean Metamorphic Complex rock, comprising mainly gneisses and charnockites Pegmatite intruded and laterite cappings and quaternary deposits. In a geological perspective, the coastal part in the study area underlies Quaternary fluviomarine deposits. Borehole lithology reveals that the aquifer material consists of fine to coarse grain sand, sandy soil, clay sand and small lens patches of clay (Thamarai 2014). In the fluvial and shallow marine environments, the clay that occurs as thin lenses or small patches surrounded the charnockites and these charnokites that exist under this sandy formation act as impermeable strata. The charnockite group of rocks is well exposed around Muttom, Manavalakurichi and the Eastern part of the study area. The rock shows varying composition from acid to ultramafic. The topography is mild at places in the coastal belt, not exceeding 60 m height from sea level, where crystalline rocks are highly weathered to laterite (GSI 2005). Weathering is relatively higher in granite gneissic rock rather than charnockite. The entire Kanyakumari to Colachel area is covered by hard rock formations like charnockite and gneisses. The groundwater occurrence is limited to only the weathered mantle of the hard rock. The weathered layer thickness generally ranges from 10 to 35 m below ground level. The groundwater occurrence is also limited from 1 to 10 m above mean sea level (MSL) (Perumal 2008). The study area's general water level varies from 3 to 25 m below ground level (IWS 2000). The shallow alluvial aquifers along Thamirabarani and Valliyar rivers serve as significant drinking water irrigation development for the Kanyakumari district. Alluvium occurs as the upper layer and is characterized by sand, gravel, and sandy clay with thickness ranging from 1 to 20 m in the study area. The alluvium and weathered charnockites function as an unconfined aquifer system (Perumal et al. 2010). The discharge from the aquifer is mainly by the withdrawal of water by pumping for domestic supply, mostly for drinking purposes, and for industries, both private and government agencies.

Materials and methods
Thirty groundwater samples were acquired from dug and bore wells in the study area during January and June 2018. The samples were scrutinized to interpret the chemistry variations of water quality and interpret observed groundwater anomalies in the region. Specimens were obtained in high-density 1L polyethylene bottles. The acquired samples were subjected to evaluating water quality criteria such as pH, electrical conductivity (EC), major cations, and anions. In situ measurements were taken for pH, electrical conductivity (EC), and total dissolved solids (TDS) with portable water quality analysis instruments HANNA (HI-9828, USA). Distinct titration procedures valued the total hardness (TH), major cations, and anions. Simultaneously, the sodium and potassium were analyzed using a Flame photometer instrument (DEEP VISION, Model -381).
The amount of sulfate ions was found using the UV-Visible photometer. The spatial analyst tool is used to draw maps that alleviated the spatial analysis of hydrological parameters to describe the water quality of the study area. The water quality index analysis has been done with the help of the Canadian Water Quality index (CWQI) programmed excel software (John-Mark 2006). ArcGIS (version10.1) software was used for spatial analysis of various physico-chemical factors. An inverse distance weighted (IDW) technique was used to interpolate the data spatially and enumerated the value for each grid node by inspecting the encompassing data points that lie within a user-defined search area (Burrough 1998).
The self-potential signals are naturally occurring electric fields measured at the ground surface (Naudet et al. 2003) using gradient electrode configuration with a 4-m interval. The non-polarizing electrodes P1, P2 (Fig. 2) are filled by copper sulfate crystals and copper sulfate solution. The equipment used to measure the self-potential was a WDDS-2 resistivity meter, accompanied by two non-polarizable electrodes and connecting cables. A total of twenty-four points was collected throughout the study area using the Selfpotential gradient method.
The vertical electrical sounding (VES) is carried out at 4 locations using the Schlumberger method by spacing the current electrode (AB) up to a maximum of 200 m in the study area. As the electrode spacing increases with unearthing the depth of investigation, it will always be less than the electrode spacing. The final results are produced by interpreting the observed data using one-dimensional software IPI2win.

Drinking water quality
The statistical parameters (minimum, maximum and average) that help to understand and relate the concentration of different physico-chemical parameters such as pH, electrical conductivity (EC), total hardness (TH) and total dissolved solids (TDS) and major ions (Ca 2+ , Mg 2+ , Na + , K + , Cl − , HCO 3 − , SO 4 2− ) are given in Table 1. The number of samples from both months, which exceeds the standards for the  Table 2.
In the study area, the pH values ranged between 5.5 and 7.8 with an average of 6 and from 5.4 to 7.2 with an average of 6.1 during June and January, respectively. Figure 3a shows the spatial distributions of pH for both the months in the study area. This indicates that the groundwater of the Mandaikkadu area is acidic (pH value of samples < 6.0) in nature. The presence of low-pH laterite soil is responsible for the acidity of groundwater in the Mandaikadu area. The electrical conductivity (EC value) varies from 93 to 2022 (µS/cm) in January, and it varies from 82 to 2163 (µS/ cm) during June as shown in Fig. 3b. As per WHO 2011 standard, few water samples exceed the limit for drinking during January and June. Higher EC content in the samples endorses the enrichment of salts in the groundwater ). In the study area, the TDS value ranged from 60 to 1294 mg/l during January and varied from 53 to 1384 mg/l during June (Fig. 3c). As a result of saltwater intrusion and anthropogenic activities, the TDS value is also higher than the permitted value. It is noted that the complete Manavalakurichi area has recorded a high value of EC and TDS.
Following the value of total hardness, groundwater is classified according to Sawyer and McCarthy DL 1967 in Table 3 and the distribution is shown in Fig. 3d. During January, 57% of samples and 47% of samples fall in the hard to very hard water category during June. The presence of calcium-bearing minerals, dissolution of sulfate and carbonate minerals in sediments and anthropogenic activities adds hardness in the study area (Ayman and Mohamed 2011). Calcium and magnesium are the most abundant components in the surface and groundwater that prevail mostly as  (2011) BIS (2000) WHO (2011) BIS (2000) Jan (2011) Jun (2011) pH 6.5-8.5 6.5-8. Bitter in taste bicarbonates and in the form of sulfate and chloride in small amounts (Kumar et al 2015). The concentration of calcium varied between 2 to 72 mg/l with an average of 27 mg/l in January (Fig. 4a). Manavalakurichi samples show that the Calcium level is higher due to the presence of shell-rich marine mudstone and some carbonate rocks. In January, the concentration of Mg 2+ in the study area ranges between 2 and 178 mg/l with an average of 55 mg/l and in June, it ranges between 3 and 126 mg/l with an average value of 45 mg/l as shown in Fig. 4b. Reverse cationic exchange, i.e., the replacement of calcium and magnesium ions by sodium ions in the groundwater, results in a higher concentration of Ca 2+ and Mg 2+ in some areas (Srinivas et al. 2015).
The chloride ion concentration in groundwater originates from various sources like leaching of sedimentary rocks and soils, weathering, domestic and industrial waste discharges, and saltwater intrusion (Krishnakumar et al. 2009). In the study area, Clconcentration ranges between 20-426 mg/l during January with an average of 161 mg/l, whereas in June, it ranges from 17 mg/l to 646 mg/l with an average of 141 mg/l. Clion spatial distribution is shown in Fig. 4c. The sodium concentration in January ranges from 16 to 420 mg/l, with an average value of 105 mg/l. In June, it ranges from 23 to 470 mg/l, with an average of 125 mg/l and the distribution is shown in Fig. 4d, where few samples exceed the permissible limit. In the study area, about 40% of samples in January and 30% of samples in June exceed the permissible limit for potassium. High potassium ion concentration is the result of agricultural runoff and saline water intrusion (Alfarrah and Walraevens 2018) (Fig. 5a). In the case of bicarbonate (HCO 3 − ) and sulfate (SO 4 2− ), no samples exceed the permissible limit of WHO. The spatial distribution of HCO 3 − and SO 4 2− ions in the study area is depicted in Fig. 5b, c, respectively.
In our study, spatial distribution maps of Figs .3, 4 and 5 show good groundwater quality in the Mandaikadu region, while Manavalakurichi regions appear to pose contaminated groundwater. Generally, groundwater contamination is reasoned by both natural and human-induced chemicals. In our study area, probable causes for natural contamination include nearby seawater intrusion, backwater from Valliyar river estuary and geomorphic features. The main human activities in the study area contribute to contamination are agriculture and coir industries. The above analysis results will be used to calculate a weighted average to categorize the groundwater quality for suitability of drinking purposes.

Water quality index (WQI)
WQIs have been applied worldwide and are used to assess the overall water quality within a particular region quickly and effectively. The water quality index (WQI) was calculated to enumerate the impact of natural and anthropogenic activities (Bilgin 2018). Based on the WQI of British Columbia Canadian council, the Canadian Water Quality index (CWQI) was developed. To measure CWQI, a set of 11 physical and chemical parameters such as pH, EC, TDS, TH, Ca 2+ , Mg 2+ , Na + , K + , HCO 3 − , Cl − and SO 4 2− were resolved. The water quality objective depends on three attributes of CWQI (Khan et al. 2004;Hurley et al. 2012;Gupta et al. 2017) and it can be calibrated using Eq. (1) (Lumb et al. 2006) According to CWQI, the water can be classified into five types, namely poor, marginal, fair, good and excellent, given in Table 4. If the CWQI ranking stands between poor and marginal, then the water is dangerous and needs to be purified. In the case of good and excellent ranking, the water does not need purification. (Khan et al. 2004). Finally, the fair ranking of CWQI needs to be improved for the water quality by purification. In the study area, 63% of the water samples are categorized as good to excellent in January and June. Figure 6 shows the spatial distribution map and most of the water samples in and around the Manavalakurichi area show marginal to fair condition.

Irrigation water quality
The quality criteria for determining groundwater viability for agricultural purposes include salinity indices, comprising Na%, SAR and PI.
(1) CWQI = 100 − √ (F 2 1 + F 2 2 + F 2 3 )∕1.732 Fig. 3 a Spatial distribution of pH in the study area, b spatial distribution of EC in the study area, c spatial distribution of TDS in the study area, d spatial distribution of TH in the study area ◂

Sodium Absorption ratio (SAR)
The SAR value represents the capability of soils to adsorb Na + from irrigation water and, therefore, the risk of sodium hazard to soils posed by irrigation with water containing massive amounts of sodium (Davraz and Ozdemir 2014). The estimation of sodium hazard value by the relative concentration of major cations such as Na + , Ca 2+ and Mg 2+ is expressed in meq/l, and SAR is calculated by Eq.
(2) (Hem 1985;Wen et al. 2020) The calculated values of SAR in the study area range between 0.64 and 16.86 meq/l during January and 1.45 to 9 meq/l during June with averages of 3.3 and 3.6 meq/l. The classification of groundwater samples based on SAR values is given in Table. 5. The correlation plot between SAR and EC in (USSL) salinity diagram (Richards 1954) gives a perfect idea about the salinity hazard in irrigation water (Ravikumar et al. 2010). Utmost, the samples collected in the study area belong to the C3-S1 category deciphering high salinity and low sodium type. Two samples in January and three samples in June indicate high salinity (C3) and medium sodium hazard (S2), which falls under C3-S2 classification (Fig. 7). On soil with restricted drainage, this class cannot be used (Ravikumar et al. 2011). The samples plotted in C1S1 and C2S1 indicate good irrigation quality (Singh et al. 2020).

Sodium percentage
The sodium concentration in irrigation waters is generally denoted as sodium percentage. Na% is measured in meq/l, and it can be calculated using Eq. (3).
High sodium concentration in irrigation water leads to Na + absorption by clay particles and displacing Mg 2+ and Ca 2+ that reduces the soil permeability, eventually resulting in poor internal drainage, which indirectly affects plant growth (Ravikumar 2011). Wilcox diagram (Wilcox 1955) for irrigation water classification correlates Na% and EC (Fig. 8). It shows that 47% of samples come under the grade of very good to good, 40% of samples belong to good to permissible, and the remaining 13% of samples the permissible to the doubtful cadre for irrigation purposes in January. During June, 60% of samples fall in very good to good, 30% of samples fall in the good to permissible, 7% of samples fall in permissible to doubtful, and the remaining 3% fall in doubtful to unsuitable for irrigation.

Permeability index
Permeability index (PI) is essential for irrigation water that affects soil permeability, influenced by sodium, calcium, magnesium and bicarbonate contents in the soil (Davraz and Ozdemir 2014). The PI values replicate the groundwater quality for irrigation calculated using Eq. (4) (Doneen 1964;Raghunath 1987) where all the ions were given in milliequivalents per liter.
The PI values greater than 75 indicate the excellent quality of water for irrigation. The PI values ranging between 25 and 75 fall in the insensible water quality, and ranging below 25 is unfit for irrigation (Al-Amry 2008). In January, 47% of samples fall in excellent quality, and 47% of samples under sensible quality and the remaining 6% of samples in the unsuitable category for irrigation. During June, 50% of samples were excellent, and the remaining 50% of samples with sensible quality for irrigation.

Hydrogeochemical facies
The hydrogeochemical facies interpretation is useful to determine the geochemical histories and flow patterns of groundwater (Rajkumar et al. 2019). Piper trilinear plot shows the changes in water quality within an aquifer using the two triangle-shaped fields. Each represents the composition of cations and anions, and both compositions present in the groundwater represented by the diamond-shaped field (Piper 1944). The classification of hydrogeochemical facies for groundwater in January and June plotted by piper trilinear diagram is shown in Fig. 9.
During January and June, most water samples fall in the NaCl segment, followed by mixed CaMgCl > CaNaHCO 3 > CaHCO 3 . The dominance of NaCl and CaMgCl in the study area indicates saltwater intrusion to the aquifer and anthropogenic contamination sources like domestic wastewater. However, the presence of CaNaHCO 3 and CaHCO 3 suggests recharge of freshwater and mineral dissolution (Senthilkumar et al. 2014).

Gibbs plot
Gibbs plot provides a clear idea about the close relationship between aquifer lithology and the chemical composition of groundwater. It is based on the mechanisms that control . 4 a Spatial distribution of Ca 2+ in the study area, b spatial distribution of Mg 2+ in the study area, c spatial distribution of Cl − in the study area, d spatial distribution of Na + in the study area ◂ Fig. 5 a Spatial distribution of K + in the study area, b spatial distribution of HCO 3 − in the study area, c spatial distribution of SO 4 2in the study area groundwater chemistry like rock weathering, atmospheric precipitation and evaporation (Gibbs 1970). Two plots proposed by Gibbs are, i) TDS versus Cl -/ (Cl -+HCO 3 -) ii) TDS versus Na + / (Na + + Ca 2+ ) Both the plots drawn for January and June samples are shown in Fig. 10. The Gibbs plots endorse that most of the samples acquired in the study area fall in the rock dominant region. The curve relating the TDS and cation indicates the presence of cations in a few wells, maybe because of evaporation or crystallization processes. The major ion chemistry of groundwater in this sector is primarily controlled by the weathering mechanism of rocks, which has been theorized by the Gibbs diagrams (Senthilkumar and Elango 2013).

Cl − /HCO 3 − ratio
Saltwater intrusion can be identified by the ratio of Cl − to HCO 3 ions and is studied to find the source of salinity in groundwater (Shin et al. 2020). In our study area, the ratio of Cl − to HCO 3 − ions ranges from 0.3 to 9.4 in January. During June, it ranges from 0.23 to 8.5, which shows a strong positive correlation with Cl − ion concentration. This relation (ratio of Cl − /HCO 3 − vs. Cl) proves the occurrence of saltwater intrusion, and Fig. 13 shows that the groundwater is slightly affected by seawater intrusion.

Data validation using the geophysical methods
An integrated approach of using geophysical and hydrogeochemical investigations was utilized (Abou Heleika et al. 2018;Alfaifi et al. 2019). The self-potential (SP) method is an important geophysical technique that responds directly to fluid flow. It involves the passive measurement of electric potential distribution at the surface of the earth with non-polarizable electrodes. The origin of SP has two main components if the groundwater is contaminated. That is, (1) electrokinetic contribution associated with groundwater flow through the permeable soil and (2) electrochemical component associated with oxide reduction phenomena. Many researchers used the SP method to study the hydrological properties and proved it as a suitable way to identify the direction of groundwater flow (Fournier 1989;Titov et al. 2005;Arsène et al. 2018).
Twenty-four measuring points were collected throughout the study area to delineate the self-potential in millivolts. The electrical potentials resulting from underground conditions of fluid pressure, temperature and salt concentration are the primary source for self-potential distribution in the study area (Arsène et al. 2018). SP anomaly increases with increasing ion concentration in the groundwater (Giampaolo et al. 2016). Figure 14 shows the spatial distribution of self-potential, and it ranges from -243 to 330 mV. The area with freshwater shows a clear negative SP signal at the surface, whereas a higher positive SP signal indicates groundwater contamination due to seawater intrusion . The contour map suggests that contaminated water spread is from Manavalakurichi toward the north and west direction in the study area.
The geophysical resistivity method is a useful technique to differentiate the subsurface resistivity to identify the seawater intrusion. In coastal regions, the seawater intrusion is determined by various researchers using the geophysical resistivity method to distinguish the freshwater and saltwater zones (Senthilkumar et al. 2019;Maiti and Gupta 2020). Especially to reveal the saline water intrusion, the one-dimensional array method in VES is useful (Song et al. 2006). The interpreted VES curves are shown in Fig. 15, which reveals the layer's resistivity  Table 6).
The resistivity pseudo-cross sections were constructed along five profiles based on the resistivity values (Fig. 16.). The pseudo-cross sections are prepared using IPI2win software, and it gives a simple two-dimensional image that provides an idea about the subsurface based on resistivity. The resistivity cross sections enable us to decipher subsurface resistivity distribution in both horizontal and vertical directions. For the generation of pseudo-cross section, a maximum number of VES locations have been considered along with the selected profiles. The pseudo-section profile is prepared by joining the VES stations VES 4, VES 3, VES 2, VES 1 and VES 5 along the west to east direction. On comparison of resistivity pseudo-cross-section map with Public Works Department (PWD) bore hole data and the geological information of the study area, depthwise subsurface layers in the form of sand, sandy clay and clayey sand are identified.
In the coastal region of the study area, the aquifer material consists of fine to coarse grain sand, sandy clay, clayey sand and small patches of clay occurring as lenses. These are identified using resistivity pseudo-cross section and available borehole lithology map of the study area. According to Perumal (2008), the clay that exists as thin lenses or small patches is mostly deposited in fluvial and shallow marine environments in the study area. Geologically, we see the presence of clay in a shallow depth, which has the characteristics to control the seawater intrusion. Thereby, VES aquifer resistivity of the coastal line of the study area whose value is low but not less than 5 is indicated as moderately affected by coastal salinity, whereas the other areas of western and southern parts (VES-4, VES-5 and VES-6) are noted as high resistivity values indicating the presence of fresh groundwater for measured depth.

Conclusion
In the current study, the hydrogeochemical analysis reveals that the groundwater in the Mandaikkadu area is acidic. The concentrations of TDS and EC were observed to be high due to saltwater intrusion in the coastal region of Manavalakurichi. From this study, in and around the area of Manavalakurichi, looking at the WQI value, it is clear that the water needs to be purified to improve the quality before use. Major cations and anions in the study are Na + and Cl − . Therefore, Fig. 15 The Schlumberger responses of interpreted VES curve the greater number of the samples comes under Na-Cl and mixed CaMgCl type in piper plot, which is mainly due to saltwater intrusion. Cl − /HCO 3 − ratio indicates that most of the samples show low to moderate seawater intrusion in the study area. No major differences in physico-chemical parameters of groundwater quality during January and June were noticed. It may be due to the prevailing hydrogeological condition in the study area. The geophysical self-potential method reveals a high positive value in Manavalakurichi, indicating saltwater intrusion in a particular location. Moreover, the resistivity studies confirm the presence of contaminated water due to saline water intrusion with a higher ionic concentration near the coastal region. From this result, marine water influencing the chemistry of groundwater can be proved. We recommend the use of water from the northwest area, where the water quality is best suitable for drinking. There is a need to properly monitor surrounding wells and Valliyar river backwater salinity to prevent the increase in saltwater contamination in the Manavalakurichi region.
Funding This research did not receive any specific grand from funding in the public, commercial, or not-for-profit sectors.

Data availability
The data were collected and analyzed primarily by us.

Compliance with ethical standard
Conflict of interest All authors have participated in (a) conception and design or analysis and interpretation of the data; (b) drafting the article or revising it critically for important intellectual content; and (c) approval of the final version. This manuscript has not been submitted to nor is under review at, another journal or other publishing venue. The authors have no affiliation with any organization with a direct or indirect financial interest in the subject matter discussed in the manuscript. On behalf of all authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.