Assessing water quality in aquifers has become increasingly important as water demand and pollution concerns rise. In the Yucatan Peninsula, sinkholes, locally known as cenotes, are karst formations that intercept the water table. Cenotes are distributed across the peninsula, but are particularly dense and aligned along a semicircular formation called the Ring of Cenotes. This area exhibits particular hydrogeological properties because it concentrates, channels, and discharges fresh water toward the coasts. In this study, we identify spatial and temporal variations in chemical and physical variables at twenty-two cenotes to identify groups that share similar characteristics. Water samples from each cenotes were taken at three depths (0.5, 5.5, and 10.5 m) and during three seasons (dry, rainy, and cold-fronts season). Field measurements of pH, temperature, electrical conductivity, and dissolved oxygen were taken, and the concentrations of major ions (K + , Na + , Mg 2+ , Ca 2+ , HCO 3 ) were quantified. Identifying regions of the cenotes were done by applying multivariate statistical techniques (PCA, PERMANOVA, CAP). The chemical variables revealed spatial trends among the cenotes. We identified three main regions. Region 1 is associated with sea-water encroachment and high levels of sulfate that travel through preferential groundwater flowpaths from evaporites in the southern Yucatan Peninsula; Region 2 is a recharge zone, and Region 3 is characterized by sea water encroachment and by the high chemical and physical variability associated with groundwater flow from the east.