Integrating magnetic susceptibility, hydrogeochemical, and isotopic data to assess the seawater invasion in coastal aquifers of Digha, West Bengal, India

Seawater intrusion in coastal aquifers is a major concern due to geogenic and anthropogenic activities leading to declining groundwater quality. The present study focuses on deciphering the sea water intruded zones and its extent in the Quaternary alluvial aquifer system in the coastal belt of Digha, West Bengal, India. In this study, 36 groundwater samples were collected during pre-monsoon (2020). Subsequently, an integrated approach of hydrogeological, hydrogeochemistry, bulk magnetic susceptibility, isotopic, multivariate statistical, and geochemical modeling is adopted. Spatial distribution maps of hydrological parameters (salinity, conductivity, TDS) and major ion concentration (Na⁺, K⁺, Ca²⁺, Mg²⁺, Cl⁻, SO₄²⁻, F⁻, and Br⁻) suggest that the northern, south-west, and eastern parts of the study area are largely affected by saltwater intrusion and are corroborated with seawater mixing index (SMI). Based on sodium adsorption ratio (SAR), sodium percentage (Na%), and Permeability index (PI) distribution maps, the same locations are identified under critical condition for the suitability of groundwater for irrigation. The order of concentration of cations and anions in the water samples are Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ and HCO₃⁻ > SO₄⁻ > Cl⁻ > Br⁻ > F⁻ respectively. Piper diagram shows three principal hydrochemical water types with water composition changes from fresh (86%) to saline water mix (14%). The hydrochemical facies evolution diagram depicts 81% of water samples are in the freshening phase, and 19% are in the intrusion phase. The various bivariate plots revealed that ion exchange, reverse ion exchange, silicate weathering, seawater mixing, and anthropogenic inputs are the governing factors that control groundwater evolution. R-mode factor analysis, principal component analysis (PCA), and agglomerative hierarchical cluster (AHC) also indicate the influence on groundwater from seawater mixing and/or seawater intrusion. The superlativeness of bulk magnetic susceptibility (χ) analysis of water samples in delineating seawater intruded zones is elaborated. Saturation index (SI) shows that groundwater is saturated (> 0) with calcite, dolomite, and aragonite, plausibly due to seawater ingression. Stable isotopic analysis of δ²H (− 53.979 to − 16.9578‰) and δ¹⁸O (− 7.00183 to − 1.37 ‰) suggests precipitation recharge/paleo-water at some locations and evaporation enrichment of groundwater. It is recommended to increase groundwater recharge, reduce groundwater extraction at critically affected locations, and have regular monitoring and management to control seawater intrusion.