This study investigated the groundwater dynamics within the Cuitzeo groundwater flow system, an area encompassing central Mexico's second and third largest lakes. We combined water level observations, structural geology, hydrogeochemical data, stable isotopes, with statistical analysis and mixing models to investigate groundwater flow dynamics and geochemical evolution in the system.