Investigating Groundwater Condition and Seawater Intrusion Status in Coastal Aquifer Systems of Eastern India

Providing sustainable water supply for domestic needs and irrigated agriculture is one of the most significant challenges for the current century. This challenge is more daunting in coastal regions. Groundwater plays a pivotal role in addressing this challenge and hence, it is under growing stress in several parts of the world. To address this challenge, a proper understanding of groundwater characteristics in an area is essential. In this study, spatio-temporal analyses of pre-monsoon and post-monsoon groundwater levels of two coastal aquifer systems (upper leaky confined and underlying confined) were carried out in Purba Medinipur District, West Bengal, India. Trend analysis of seasonal groundwater levels of the two aquifers systems was also performed using Mann-Kendall test, Linear Regression test, and Innovative Trend test. Finally, the status of seawater intrusion in the two aquifers was evaluated using available groundwater-quality data of Chloride (Cl−) and Total Dissolved Solids (TDS). Considerable spatial and temporal variability was found in the seasonal groundwater levels of the two aquifers. Further, decreasing trends were spotted in the pre-monsoon and post-monsoon groundwater-level time series of the leaky confined and confined aquifers, except pre-monsoon groundwater levels in Contai-I and Deshpran blocks, and the post-monsoon groundwater level in Ramnagar-I block for the leaky confined aquifer. The leaky confined aquifer in Contai-I, Contai-III, and Deshpran blocks and the confined aquifer in Nandigram-I and Nandigram-II blocks are vulnerable to seawater intrusion. There is an urgent need for the real-time monitoring of groundwater levels and groundwater quality in both the aquifer systems, which can ensure efficient management of coastal groundwater reserves.


Introduction
Water is a vital substance on earth and forms the principal constituent of all living things. Therefore, it is called the lifeblood of the biosphere. The unabated population growth and implacable rise of water demand in different sectors have huge repercussions for a growing freshwater scarcity in several parts of the world, including India. Providing sustainable water for domestic use and irrigated agriculture is one of the big challenges for the 21st century. To address this challenge, efficient planning and management of water resources are of utmost importance. Groundwater is a renewable but finite resource. Generally, it is preferred for drinking water supply because of its good quality, pleasant taste, and from a safety point of view. Due to rapid population growth, water resources potential is decreasing day by day in general, and groundwater resources potential in particular. Furthermore, mismanagement of water resources and climate change creates an imbalance between water supply and demand. For proper planning and management of groundwater, assessment of groundwater is very important. Nowadays, the management of coastal aquifer systems is becoming more essential due to dramatic climate change. Unfortunately, groundwater management has been neglected in coastal areas, though coastal regions provide one-third of the world's ecosystem services and natural capital [1]. Anthropogenic activities such as excessive pumping of coastal aquifers and increased impervious surfaces in urbanized areas are the major causes of saltwater intrusion. Moreover, the increasing sea level due to climate change could aggravate this problem. Human water use is a key driver in the hydrology of coastal aquifers, but water management in the coastal areas is misguided [2].
For efficient planning of water resources, statistical analyses of hydrological time series play a vital role by identifying and describing quantitatively each of the hydrological processes [3] and forecasting future values of the time series variable. The trend is one of the essential characteristics of hydrological time series and it exists in a dataset if there is a significant correlation between observation and time. Different parametric and nonparametric methods have been used for detecting monotonic trends. However, efficient management of water resources requires both monotonic trends as well as separate trends for low, medium, and high values of the concerned variable. Such an approach helps to find drought and flood occurrences in their increasing or decreasing frequencies [4]. A novel method has been proposed by Sen (2012) [4], which is capable of determining monotonic trends as well as trends for low, medium, and high values. In the past, many studies have been focused on the detection of trends in rainfall [5][6][7][8][9][10][11][12][13][14]. The application of a set of statistical tests for the same objective in a time series analysis increases the chance of rejecting a true null hypothesis [5]. Therefore, it is important to analyze the results of at least more than two statistical tests while deciding on rejecting the null hypothesis.
India, with a vast coastline in order of about 7500 km, has a network of major and minor delta systems and estuaries all along the coastal regions [15]. Coastal aquifers serve as major freshwater reservoirs in the coastal areas. Due to rapid population growth, stresses on coastal aquifers are increasing day by day. Nowadays, the management of coastal aquifer systems is becoming more important due to global climate change. Serious problems of saltwater intrusion exist in many coastal regions of the country including West Bengal, Odisha, Andhra Pradesh, and Tamil Nadu. In India, seawater intrusion is observed along with the coastal areas of Gujarat and Tamil Nadu (http://cgwb.gov.in/documents/ papers/INCID.html, accessed on 15 April 2016). A better understanding of the patterns of movement and mixing of freshwater and seawater, as well as the factors influencing these processes, is necessary to manage coastal groundwater resources efficiently. Bhosale and Kumar (2002) [16] used SUTRA software package for the simulation of seawater intrusion in Ernakulam coast, Kerala (India) and found that the most sensitive zone for seawater intrusion laid between 400 and 2000 m from the high tide line in Ernakulam coast. Elango and Sivakumar (2008) [17] carried out a modeling study using MODFLOW on seawater intrusion in a coastal aquifer of Chennai (India) and the model showed that the total abstraction from this aquifer should be restricted to 4.25 MGD to prevent seawater intrusion. Datta et al. (2009) [18] used FEMWATER for modeling and control of seawater intrusion in a coastal aquifer in two Blocks of Nellore District, Andhra Pradesh (India). Ganesan and Thayumanavan (2009) [19] identified management strategies by the simulation of groundwater flow and solute transport in the seawater intruded coastal aquifer system of Chennai (India) using MODFLOW and MT3D. Simulation-optimization modelling studies were carried out by Sindhu et al. (2012) [20], Lathashri and Mahesha (2015) [21], and Mohan and Pramada (2015) [22] in Trivandrum, Kerala (India), sea-coast of Karnataka, and southern Chennai respectively.
In Eastern India, Rejani et al. (2007) [23] carried out simulation-optimization modelling in Balasore coastal basin, Odisha and established the most promising management strategy of reducing pumpage from the second aquifer by 50% and an increase in pumpage up to 150% from the first and second aquifer. Goswami (1968) [24] applied a geochemical approach and Shivanna et al. (1993) [25] employed environmental isotopes D, 18O, 34S, 3H, and 14C to examine seawater salinity in the coastal part of Purba Medinipur District, West Bengal. Bhattacharya and Basack (2012) [26] found out the extent of seawater intrusion by analyzing EC, Chloride, and TDS concentration of groundwater in coastal areas of Contai-Digha, West Bengal.
The literature review reveals that only a couple of preliminary geochemical investigations have been conducted in the coastal portion of West Bengal. No comprehensive scientific studies on seawater intrusion have been conducted to date in general and in West Bengal in particular. West Bengal has its own share of water problems in India. The erratic nature of the southwest monsoon has rendered surface water supply unreliable and inadequate, which in turn has compelled the people to depend more and more on groundwater for meeting the growing water demands of different sectors. Therefore, there is an urgent need for integrated management of water resources in this state using modern tools/techniques to ensure their long-term sustainability. Keeping this in mind, the present study has been undertaken with an overall objective to investigate groundwater conditions and evaluate seawater intrusion status in the coastal aquifer systems of Purba Medinipur district, West Bengal.

Overview of Study Area
The study area selected for this study lies in the coastal area of Purba Medinipur District, West Bengal, India ( Figure 1). It is located in the lower deltaic part of the Ganga basin and lies between latitude 21 • [26] found out the extent of seawater intrusion by analyzing EC, Chloride, and TDS concentration of groundwater in coastal areas of Contai-Digha, West Bengal. The literature review reveals that only a couple of preliminary geochemical investigations have been conducted in the coastal portion of West Bengal. No comprehensive scientific studies on seawater intrusion have been conducted to date in general and in West Bengal in particular. West Bengal has its own share of water problems in India. The erratic nature of the southwest monsoon has rendered surface water supply unreliable and inadequate, which in turn has compelled the people to depend more and more on groundwater for meeting the growing water demands of different sectors. Therefore, there is an urgent need for integrated management of water resources in this state using modern tools/techniques to ensure their long-term sustainability. Keeping this in mind, the present study has been undertaken with an overall objective to investigate groundwater conditions and evaluate seawater intrusion status in the coastal aquifer systems of Purba Medinipur district, West Bengal.

Overview of Study Area
The study area selected for this study lies in the coastal area of Purba Medinipur District, West Bengal, India ( Figure 1). It is located in the lower deltaic part of the Ganga basin and lies between latitude 21°34′41″ N and 22° 00′45″ N and longitude 87° 18′32″ E and 88° 03′22″ E. It covers an area of about 1398.7 km 2 bounded by Haldi River in the north, Hugli River outfall in the northeast, the Bay of Bengal in the southeast, Subarnarekha River in the south, Orissa Coast Canal (OCC) in the south-west and Hugli Tidal Canal (HTC) second range in the northwest (Figure 1). The area is characterized by gently sloping flat alluvial terrain, which gradually merges to the deltaic plain towards the south. The average annual rainfall in the area is 1700 mm, out of which 73% occurs during the monsoon season lasting from June to September. The study area has a humid subtropical type of climate with maximum and minimum temperatures of 39 °C in summer and 12 °C in winter, respectively. The study area is underlain by unconsolidated alluvial sediments of quaternary age, and is overlain by linear continuous beach ridges (sand dunes) along the coastline around Digha-Ramnagar-Contai deposited by the wave action of the sea. The quaternary sediments are underlain by tertiary sediments of the Mio-Pliocene age (http://cgwb.gov.in/District_Profile/WestBangal/Purba%20Medinipur.pdf, accessed on 10 May 2016). It covers an area of about 1398.7 km 2 bounded by Haldi River in the north, Hugli River outfall in the northeast, the Bay of Bengal in the southeast, Subarnarekha River in the south, Orissa Coast Canal (OCC) in the south-west and Hugli Tidal Canal (HTC) second range in the northwest (Figure 1). The area is characterized by gently sloping flat alluvial terrain, which gradually merges to the deltaic plain towards the south. The average annual rainfall in the area is 1700 mm, out of which 73% occurs during the monsoon season lasting from June to September. The study area has a humid subtropical type of climate with maximum and minimum temperatures of 39 • C in summer and 12 • C in winter, respectively. The study area is underlain by unconsolidated alluvial sediments of quaternary age, and is overlain by linear continuous beach ridges (sand dunes) along the coastline around Digha-Ramnagar-Contai deposited by the wave action of the sea. The quaternary sediments are underlain by tertiary sediments of the Mio-Pliocene age (http://cgwb.gov.in/District_Profile/WestBangal/Purba%20Medinipur.pdf, accessed on 10 May 2016).

Datasets Used
The monthly total rainfall (1976-2015), pre-monsoon and post-monsoon average groundwater-level (2001-2015) and groundwater quality (2011-2014), groundwater irrigation (2006-2007 and 2013-2014), and lithological data required for the study were collected from various state and central government organizations. The collected hydrological and meteorological monthly data were processed in block-level and analyzed in time series using different statistical methods. Moreover, the mapping and spatial distribution of hydrological data were done on the ArcGIS platform.

Detection of Trend in Groundwater-Level Time Series
In the present study, powerful statistical tests, namely the 'Linear Regression test', 'Mann-Kendall test' and 'Innovative Trend test' proposed by Sen (2012) [4] were used to find trends in the pre-monsoon and post-monsoon groundwater-level time series (2001-2015) of the leaky confined and confined aquifers in the study area for the available data. To carry out trend analysis, spreadsheet programs were developed using MS Excel software. Brief descriptions of these tests are given below.

Linear Regression Test
The least-square linear regression test is a parametric test. Before the application of this test, a normality check should be done as hydrologic and climatic data are skewed most of the time. This test used to describe the presence of linear trends in time series [27].
where m is the estimated slope of the regression line between the observed value and time and S m stands for the standard error of the estimated slope. The test statistic (t) follows a student's t-distribution with n − 2 degree of freedom, where n is the number of samples. The null hypothesis of slope zero will be rejected when the test statistic (t) value is greater than the critical value t α/2 with α significance level.

Mann-Kendall Test
The Mann-Kendall test is a non-parametric test for exploring trends in a time series without specifying the type of trend (i.e., linear or non-linear). This test has been found to be an excellent tool for trend detection [28]. Using Mann-Kendall's test, it is possible to find the existence of a rising or falling trend, but it is hard to quantify its magnitude. The Mann-Kendall test statistic (S) is calculated as: where, n is the number of data points, x i and x j are the data values in time series i and j (j > i), respectively, and sgn (x j − x i ) is the sign function defined as: Moreover, for n > 10, the standard normal test statistic Z S is computed as: where, In Equations (4) and (5), p = 1 for S < 0 and p = −1 for S > 0, k is the number of tied groups (a set of sample data having the same value), and t i is the numbers of data in the ith tied group. The value of the test-statistic Z S is considered as zero for S = 0. Positive values of Z S indicate rising trends while negative Z S values show falling trends. If the computed absolute value of Z S is greater than the critical value of the standard normal distribution, the hypothesis of a rising or falling trend cannot be rejected at the α significance level.

Innovative Trend Test
The innovative trend test proposed by Sen (2012) [4] depends on a 1:1 line on the Cartesian coordinate system, irrespective of the assumptions of distribution, sample length, and serial correlation. In this test, the hydrologic time series is divided into two equal halves and is sorted in ascending order. Thereafter, the first half of the series is plotted on the horizontal axis and the second half on the vertical axis leading to a graph with a 1:1 (45 • ) straight-line on the Cartesian coordinate system. If all the data points lie exactly or close to the 1:1 line, it can be said that there exists no trend in the time series. The presence of data points above or below the 1:1 line indicates increasing or decreasing monotonic trends in the time series, respectively.
In this test, the null hypothesis, H 0 , implies that there is no significant trend if the calculated slope value, s, remains below a critical value, s cri . Here, s cri , represents the critical standard deviation for standardized time series at ±1.96 (1.65) for 95% (90%) significance levels (α).
The slope of the trend (s) is expressed as: In Equation (6), n is the numbers of data in whole time series and y 1 , y 2 are the first and second half time series average.
Otherwise, an alternative hypothesis, Ha, is valid when s > s cri . If at α percent significance level, the confidence limits of a standard normal probability density function (PDF) with zero mean and standard deviation s cri , then the confidence limits (CL) of the trend slope can be represented by the following expression: where, In Equations (7) and (8), σ s = standard deviation of trend slope, σ = standard deviation of the whole time series, ρ y 1 y 2 = cross-correlation coefficient between the sorted two halves arithmetic averages in ascending order, and n = numbers of data in whole time series.
If the slope (s) value falls outside the lower and upper confidence limits, the null hypothesis of no trend is rejected.

Evaluation of Pre-Monsoon and Post-Monsoon Groundwater-Levels
In India, the monsoon season spans from July to September, during which maximum groundwater recharge takes place. Generally, well hydrograph follows a similar trend like stream hydrograph. Groundwater-levels during April and November are taken as the representatives of pre-monsoon and post-monsoon groundwater-levels, respectively. In this study, the pre-monsoon (April) and post-monsoon (November) groundwater-levels were plotted in the same graphs for the years 2001-2015 for both leaky confined and confined aquifers. Block was considered as the spatial unit for this analysis. The analysis was performed for 9 blocks and 6 blocks underlying confined and leaky confined aquifers, respectively.

Evaluation of Seawater Intrusion Status
The spatial variation of Cl − and TDS concentrations of groundwater of the leaky confined and confined aquifers during pre-monsoon and post-monsoon seasons were analyzed for the years 2011-2014. The Inverse Distance Weighted (IDW) technique available in ArcGIS was used to interpolate Cl − and TDS concentrations in groundwater. Furthermore, the interpolated raster files were classified based on 250 mg/L interval for Cl − and TDS concentrations up to 1000 mg/L, and then 1000 mg/L interval for the same beyond 1000 mg/L concentration.

Trend in Groundwater-Levels of the Leaky Confined Aquifer
The Mann-Kendall (MK) test, Linear Regression (LR) test, and Innovative Trend (IT) test were used for trend detection at a 5% significance level (α = 0.05). The results of the trend test for pre-monsoon and post-monsoon groundwater-levels are summarized in Tables 1 and 2, respectively.  MK and LR test results show that pre-monsoon groundwater-levels (elevations) of leaky confined aquifers for Nandigram-I, Nandigram-II, and Contai-III blocks have a decreasing trend. However, other blocks have no trend in the pre-monsoon groundwater level of the leaky confined aquifer. The MK test and LR test show no trend in the premonsoon groundwater levels of Ramnagar-I and Deshpran, whereas the Innovative trend test shows a significant decreasing trend as the Innovative trend test is more sensitive and capable to capture mild trends in the hydrologic time series. No trend of pre-monsoon groundwater-level (elevation) is shown from MK, LR and IT test results for Contai-I block, most probably due to availability of short-time series data (Table 1). It can be observed from Table 2 that all the blocks reveal a decreasing trend in the post-monsoon groundwater level. In the case of Contai-I and Ramnagar-I blocks, the results of LR and MK tests do not agree with those of the Innovative trend test.
The majority among the three test results is taken into consideration to decide the presence or absence of time series data of groundwater-level levels.

Trend in Groundwater-Levels of the Confined Aquifer
There are decreasing trends for pre-monsoon and post-monsoon groundwater elevations of confined aquifers in all the blocks. The calculated test statistic values are greater than the critical test statistic values for all blocks based on these three tests. For the confined aquifer, the results of the Innovative Trend test support those of the two traditional tests.

Seasonal Groundwater Fluctuation
Average groundwater elevations over piezometers (monitoring wells) present in the study area were taken into account for the study of pre-monsoon and post-monsoon groundwater fluctuation at a block level for both the aquifers. For ease of analysis, the study area was divided into three parts, i.e., North-East (Nandigram-I, Nandigram-II, Khejuri-I, and Khejuri-II), Central Part (Contai-I, Deshpran, and Contai-III), and Southwest (Ramnagar-I, Ramnagar-II, and Bhograi).

Groundwater Fluctuations in the Leaky Confined Aquifer
The pre-monsoon and post-monsoon groundwater-level-time series data for 6 blocks are shown in Figure 2a-f along with the annual rainfall time series data. Among the 6 blocks, pre-monsoon and post-monsoon groundwater-levels (GWL) are below the mean sea level (MSL) in Nandigram-I, Nandigram-II, and Contai-I blocks (Table 3).   Overall, the pre-monsoon GWL varies from 0.18 to −9.07 m, whereas the postmonsoon GWL varies from 2.06 to −7.82 m. This should be because of enhancement of groundwater extraction due to an increased number of irrigation Shallow Tube Wells (STW) and additional installation of a good number of irrigation Medium Deep Tube Wells (MDTW) after 2006-2007 ( Figure 3). It can be observed from Figure 2a-f that the North-East part shows a more decreasing trend in both pre-monsoon and post-monsoon seasons. Nandigram-I block has a less decreasing trend of groundwater-level compared to Nandigram-II in the post-monsoon groundwater-level, but Nandigram-I is more vulnerable to seawater intrusion as compared to Nandigram-II as it is located near to the sea. This might be due to withdrawal of a larger quantity of groundwater through greater numbers of MDTWs in Nandigram-II as compared to Nandigram-I. In the Central part of the study area, Contai-I and Contain-III have a more declining trend compared to Deshpran, as a result of heavy withdrawal of groundwater by a large number of MDTWs in Contai-I and Contai-III as compared to Deshpran (Figure 3). In the South-West part (Ramnagar-I), pre-monsoon groundwater-level shows a declining trend, but it recovers during the post-monsoon season compared to other blocks.

Groundwater Fluctuations in the Confined Aquifer
Pre-monsoon and post-monsoon groundwater-level data for 9 blocks are shown in Figure 4a-i along with annual rainfall. Among the 9 blocks, pre-monsoon and post-monsoon groundwater-levels are below the mean sea level (MSL) in Nandigram-I, Nandigram-II, Khejuri-I, Khejuri-II, Deshpran, and Contai-III blocks (Table 4). Table 4. Temporal variation of groundwater-level in confined aquifer.

Block
Groundwater

Seawater Intrusion Status
The study areapart of the coastal areas of Purba Midnapur District possesses high geologic heterogeneity. The pattern of variation of salinity measured as chloride (Cl -) and Total Dissolved Solids (TDS) concentration is heterogeneous in three-dimensional spaces. For better understanding of seawater intrusion, the study area was divided into three parts, i.e., North-East, Central, and South-West parts.
Chloride (Cl -) and TDS concentrations of groundwater in the leaky confined aquifer show remarkable spatial and temporal variations during 2011-2014. In Nandigram-I and Ramnagar-I, both the Chloride (Cl -) and TDS concentration in groundwater remains within the WHO guideline value i.e., <250 mg/L and <1000 mg/L, respectively (https://www.researchgate.net/figure/World-Health-Organization-WHO-Drinking-Wate r-Qualy-Guidelines_tbl2_255700290, accessed on 10 July 2016). During the pre-monsoon of 2014, Clconcentration greater than 250 mg/L was found in Nandigram-II and Contai-I. In Deshpran, both the Cland TDS concentrations fluctuate seasonally and annually, and high Cland TDS concentrations were noticed during the pre-monsoon of 2014. Furthermore, in Contai-III block, Cland TDS concentrations are above the WHO guideline value.
Significant spatial variation and comparatively higher concentration of Cland TDS in pre-monsoon groundwater were found in the confined aquifer [ Figure 6(a-h)]. The central part of the study area shows a higher concentration of Cl -(> 250 mg/L) and TDS (> 1000 mg/L) i.e., beyond the WHO guideline value. This is because of high tidal saline water ingress into the coastal aquifers through leakage from the beds of Rasulpur River (passing through Deshpran), Pichaboni Khal (passing through Contai-I), and Mirzapur Khal (connector of Rasulpur River and Pichaboni Khal; passing through both Deshpran and Contai-I). Furthermore, the confined aquifer in the North-East part (particularly in Nandigram-I (part) and Nandigram II) shows high Cl -(>250 mg/L) and moderate TDS (750-1000 mg/L). This may be due to the saline water ingress through leakage from the Haldi River bed, and declining trend of groundwater-level during both pre-monsoon and post-monsoon seasons.
Based on distribution of Chloride (Cl -) and TDS concentration in groundwater, it can be inferred that the effect of seawater intrusion is persistent in the leaky confined aquifer in Contai-III block (Central part of the study area). However, in Deshpran block, both the chloride and TDS concentrations fluctuate seasonally and annually. Mild effect of sea-

Seawater Intrusion Status
The study area-part of the coastal areas of Purba Midnapur District possesses high geologic heterogeneity. The pattern of variation of salinity measured as chloride (Cl − ) and Total Dissolved Solids (TDS) concentration is heterogeneous in three-dimensional spaces. For better understanding of seawater intrusion, the study area was divided into three parts, i.e., North-East, Central, and South-West parts.
Chloride (Cl − ) and TDS concentrations of groundwater in the leaky confined aquifer show remarkable spatial and temporal variations during 2011-2014. In Nandigram-I and Ramnagar-I, both the Chloride (Cl − ) and TDS concentration in groundwater remains within the WHO guideline value i.e., <250 mg/L and <1000 mg/L, respectively (https://www. researchgate.net/figure/World-Health-Organization-WHO-Drinking-Water-Qualy-Guidelines_ tbl2_255700290, accessed on 10 July 2016). During the pre-monsoon of 2014, Cl − concentration greater than 250 mg/L was found in Nandigram-II and Contai-I. In Deshpran, both the Cl − and TDS concentrations fluctuate seasonally and annually, and high Cl − and TDS concentrations were noticed during the pre-monsoon of 2014. Furthermore, in Contai-III block, Cl − and TDS concentrations are above the WHO guideline value.
Significant spatial variation and comparatively higher concentration of Cl − and TDS in pre-monsoon groundwater were found in the confined aquifer (Figure 6a-h). The central part of the study area shows a higher concentration of Cl − (>250 mg/L) and TDS (>1000 mg/L) i.e., beyond the WHO guideline value. This is because of high tidal saline water ingress into the coastal aquifers through leakage from the beds of Rasulpur River (passing through Deshpran), Pichaboni Khal (passing through Contai-I), and Mirzapur Khal (connector of Rasulpur River and Pichaboni Khal; passing through both Deshpran and Contai-I). Furthermore, the confined aquifer in the North-East part (particularly in Nandigram-I (part) and Nandigram II) shows high Cl − (>250 mg/L) and moderate TDS (750-1000 mg/L). This may be due to the saline water ingress through leakage from the Haldi River bed, and declining trend of groundwater-level during both pre-monsoon and post-monsoon seasons.  Based on distribution of Chloride (Cl − ) and TDS concentration in groundwater, it can be inferred that the effect of seawater intrusion is persistent in the leaky confined aquifer in Contai-III block (Central part of the study area). However, in Deshpran block, both the chloride and TDS concentrations fluctuate seasonally and annually. Mild effect of seawater intrusion is found in Nandigram-I, Nandigram-II, and Contai-I blocks. However, in Ramnagar-I, the effect of seawater intrusion is not markedly noticeable. This might be due to the presence of soil of low hydraulic conductivity of the aquifer in the area [26]. On the other hand, the effect of seawater intrusion is persistent in the confined aquifer in Nandigram-II block showing greater salinity than corresponding Nandigram-I block adjacent to the seacoast. This situation occurs as the filaments of high hydraulic conductivity interspersed with soil of low hydraulic conductivity protrude deep inland from the sea causing ingress of sea water inland through this preferential path (Maity et al., 2017) [29]. However, in the Central part of the study area (i.e., Contai-I, Deshpran, and Contai-III blocks), the concentrations of both Cl − and TDS in groundwater of confined aquifer vary considerably from season to season as well as from year to year.

Conclusions
The goal of the present study was to investigate the spatio-temporal characteristics of groundwater-level and evaluate the seawater intrusion status of the leaky confined and confined aquifers. For this, trend analysis of seasonal groundwater-levels (i.e., pre-monsoon and post-monsoon seasons) was carried out by the Mann-Kendall test, Linear Regression test, and Innovative Trend test for both leaky-confined and confined aquifers. Thereafter, the fluctuations of pre-monsoon and post-monsoon groundwater-levels with time and space were analyzed. As the groundwater-levels during pre-monsoon and post-monsoon season decline below the mean sea level leading to a reversal of hydraulic gradient, it causes inland movement of seawater in aquifers of some coastal blocks and hence, they are vulnerable to seawater intrusion. By analyzing the chloride (Cl − ) and TDS concentrations of groundwater, the status of seawater intrusion in the study area was evaluated. The methodology of this study can be useful for management of the groundwater resources in other coastal regions of India as well as other coastal parts of the world. Statistical analysis of available datasets supported with simulation-optimization modelling in this study area can be suggested for future research to address the seawater intrusion problem more efficiently.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.