Reconstruction of cave air temperature based on surface atmosphere temperature and vegetation changes: Implications for speleothem palaeoclimate records

Cave temperature is a parameter of major importance for palaeoclimate studies based on speleothems since most proxies used are temperature-dependent. The general assumption that cave temperature reflects the mean annual surface atmosphere temperature (SAT) over the cave is based on a limited number of specific studies, and detailed mechanisms supporting this link are not properly described. Eagle Cave, in Spain, was used to understand the detailed mechanisms connecting SAT and cave temperature. A monitoring programme conducted from 2009 to 2011 allows the characterization of the thermal dynamics in the cave air. In the studied karstic environment thermal conduction is the main mechanism responsible for transferring the ground temperature signal to the cave. The calculated thermal diffusion coefficient is 0.756×10^-6±0.013m²s^-1. The SAT signal is recorded at the cave with a delay ranging from ~3 to 11yr, depending on the thickness of the bedrock over the ceiling (5-18m). A cave cooling of 2°C has been recorded since mid-1970s, when cave temperature measurements were made. The application of a model exclusively based on thermal conduction of SAT does not reproduce the observed thermal changes. However, vegetation cover over the cave was modified dramatically during recent decades. This impacts ground temperature due to the changes in exposure to insolation and the modification of soil properties, which is eventually recorded in the cave. The incorporation of changes in vegetation cover in the model provides a realistic thermal reconstruction. The impact of vegetation cover was much more significant than SAT in order to explain the large thermal change in the cave over past decades. Thus, the cooling in Eagle Cave during past decades is mainly controlled by a variable decoupling from mean annual SAT as a result of the afforestation over the cave. The effect of this thermal decoupling on speleothem δ¹⁸O records is potentially significant and should be considered when performing accurate palaeoclimate reconstructions.