Here we track the water-table position and temperature of the Mount San Giovanni aquifer (Iglesiente-Sulcis mining district, SW Sardinia, Italy) during the past 600 ka by determining the ages (U/Th dating) and stable isotope compositions (δ18O, δ13C and Δ47) of a variety of subaqueous carbonate speleothems (e.g. calcite spars, dogtooth calcite crystals and calcite coatings). Clumped isotopes (Δ47) provide quantitative estimates of carbonate formation temperatures (and thus water temperatures) that are independent of the oxygen isotope composition of water (δ18Ow). Then, the δ18Ow of the paleo-water has been reconstructed from the clumped isotope temperature (TΔ47) and the δ18O of the carbonate (δ18Oc). We find that some high-temperature calcite spars formed prior to 600 ka at temperatures above ~120 °C. Lower-temperature spars (~70 °C) precipitated at ~400 ka, and cold-water subaqueous speleothems (~10–20 °C) formed in perched ponds at different levels of the karst systems between 410 ka and 110 ka, while coeval precipitation of subaerial flowstones occurred in the upper level of the shallower caves until 82 ka. We infer that the groundwater level dropped by ~120 m from ~400 to ~250 ka, with a relatively rapid rate of ~0.8 mm/yr. Considering the tectonic stability of Sardinia during the Quaternary, this high rate derives from climate driven geomorphological processes at the surface rather than tectonic uplift. The δ18O values of the paleo-aquifer water range from −6.0 ± 0.7‰ during MIS 5c, and similar to modern cave water values (−5.1‰), to −7.7± 0.4‰ during the colder MIS 8. These values indicate that the groundwater reflected the δ18O signal of meteoric water, with no significant contributions from metasomatic and metamorphic waters. The observed δ18Ow variability can be explained by glacial/interglacial paleoclimate changes affecting rainfall δ18O (and thus groundwater δ18O). We conclude that clumped isotope thermometry on subaqueous carbonate speleothems is a useful tool for tracking paleo-aquifer temperatures and δ18Ow reconstructions.