Groundwater Quality Evolution Model in the Ring of Cenotes, Yucatan, Mexico

Karst aquifers show dissolution/precipitation processes of the minerals present in the carbonate rocks. The Ring of Cenotes (RC) extends along the edge of the Chicxulub crater, in the limestone platform of the Yucatan Peninsula (YP), where groundwater shows preferential flow paths toward the coast near Celestun and Dzilam Bravo towns. This study aimed to describe the regional hydrogeochemical evolution of groundwater of the RC, and its association with the dissolution/precipitation of the minerals present along its path to the ocean. To achieve this aim, we: a) characterized groundwater's hydrogeochemistry; b) determined the calcite, dolomite, and gypsum saturation indexes (reaction phases with the groundwater) in the study area; c) proposed a hydrogeochemical model developed through PHREEQC using an inverse modelling approach. The model predictions confirmed that there are two evolution pathways of the groundwater consistent with the preferential flow paths suggested in a previous regionalization of the RC. On the western path, where groundwater flows towards Celestun, an important marine intrusion influences the hydrogeochemical processes and represents a risk for the prevalence of freshwater. On the eastern path, where groundwater flows toward Dzilam Bravo, the hydrogeochemistry in the sinkholes correlates well with rainfall, suggesting a higher vulnerability during droughts than during rainy periods.