Abstract
The critical role of groundwater in meeting diverse needs, including drinking, industrial and agricultural, highlights the urgency of effective resource management. Excessive groundwater extraction, especially in coastal regions including Urmia Plain in NW Iran, disrupts the equilibrium between freshwater and saline boundaries within aquifers. Influential parameters governing seawater intrusion – groundwater occurrence (G), aquifer hydraulic conductivity (A), the height of groundwater level above the mean sea level (L), distance from the shore (D), impact of the existing status of seawater intrusion (I), and thickness of the saturated aquifer (T) – merge to shape the GALDIT vulnerability index for coastal aquifers. This study enriches the GALDIT framework by incorporating two additional hydrogeological variables: hydraulic gradient (i) and pumping rate (P). This expansion produces seven distinct vulnerability maps (GALDIT, GAiDIT, GAiDIT-P, GALDIT-i, GALDIT-iP, GALDIT-P, GAPDIT). In the Urmia Plain, the traditional GALDIT index reveals vulnerability values ranging from 2 to 8.1, categorized into six classes from negligible to very high vulnerability. However, the modified indices, GAiDIT and GAiDIT-P, yield a three-class categorization, ranging from low to high vulnerability. The introduction of the 'i' and 'P' parameters in GALDIT-i and GALDIT-iP enhances the precision of vulnerability mapping, altering class distribution and intensifying vulnerability ratings. The eastern, central, and coastal areas of the Urmia Plain demonstrate high to very high vulnerability levels, in contrast to the lower vulnerability observed in the western regions. Both the GALDIT-P (r = 0.82) and GALDIT-iP (r = 0.81) indices show strong correlations with chloride concentration, thereby improving mapping accuracy over the traditional GALDIT index (r = 0.72). A sensitivity analysis highlights the critical influence of the 'i' parameter, suggesting its weighting should be revised. Parameter recalibration serves to amplify the significance of 'G', 'L', 'D', and 'i' parameters, while diminishing others. The integration of multiple hydrogeological variables considerably enhances the precision of groundwater vulnerability assessments.