Salinity Distribution as a Hydrogeological Limit in a Karstic Watershed in Yucatan

In coastal regions, the interaction between freshwater and seawater creates a dynamic system in which the spatial distribution of salinity critically constrains the use of freshwater for human consumption. Although saline intrusion is a globally widespread phenomenon, its inland extent varies significantly with hydrological conditions, posing a persistent threat to groundwater quality and sustainability. This study aimed to characterize salinity distribution using an integrated karst-watershed approach, thereby enabling the identification of both lateral and vertical salinity gradients. The study area is in the northwestern Yucatan Peninsula. Available hydrogeological data were analyzed to determine aquifer type, soil texture, evidence of saline intrusion, seawater fraction, vadose zone thickness, and field measurements. These included sampling from 42 groundwater sites (open sinkholes and dug wells), which indicated a fringe zone approximately 5 km in size influenced by seawater interaction, in mangrove areas and in three key zones of salinity patterns: west of Mérida (Celestun and Chunchumil), and northern Yucatan (Sierra Papacal, Motul, San Felipe). Vertical Electrical Sounding (VES) and conductivity profiling in two piezometers indicated an apparent seawater influence. The interface was detected at a depth of 28 m in Celestun and 18 m in Chunchumil. These depths may serve as hydrogeological thresholds for freshwater abstraction. Results indicate that saltwater can extend several kilometers inland, a factor to consider when evaluating freshwater availability. This issue is particularly critical within the first 20 km from the coastline, where increasing tourism exerts substantial pressure on groundwater reserves. A coastal-to-inland salinity was identified, and an empirical equation was proposed to estimate the seawater fraction (fsea%) as a function of distance from the shoreline in the Cenote Ring trajectory. Vertically, a four-layer model was identified in this study through VES in the western watershed: an unsaturated zone approximately 2.6 m thick, a confined layer in the coastal Celestun profile about 9 m thick, a freshwater lens floating above a brackish layer between 8 and 25 m, and a saline interface at 37 m depth. The novelty of this study, in analyzing all karstic water surfaces together as a system, including the vadose zone and the aquifer, and considering the interactions with the surface, is highlighted by the strength of this approach. This analysis provides a better understanding and more precise insight into the integrated system than analyzing each component separately. These findings have significant implications for water resource management in karst regions such as Yucatan, underscoring the urgent need for sustainable groundwater management practices to address seawater intrusion.