Abstract. An assessment of the soil and bedrock thermal structure of the Sierra de Guadarrama, in central Spain, is provided using subsurface and ground surface temperature data coming from four deep (20 m) monitoring profiles belonging to the Guadarrama Monitoring Network (GuMNet) and two shallow profiles (1 m) from the Spanish Meteorology Service (Agencia Estatal de Meteorología, AEMET) covering the time spans of 2015–2021 and 1989–2018, respectively. An evaluation of air and ground surface temperature coupling showed that soil insulation due to snow cover is the main source of seasonal decoupling, being especially relevant in winter at high-altitude sites. Temperature propagation in the subsurface was characterized by assuming a heat conductive regime by considering apparent thermal diffusivity values derived from the amplitude attenuation and phase shift of the annual cycle with depth. This methodology was further extended to consider the attenuation of all harmonics in the spectral domain, which allowed for analysis of thermal diffusivity from high-frequency changes in the soil near the surface at short timescales. For the deep profiles, the apparent thermal diffusivity ranges from 1 to 1.3×10-6 m2 s−1, which is consistent with values for gneiss and granite, the major bedrock components in the Sierra de Guadarrama. However, thermal diffusivity is lower and more heterogeneous in the soil layers close to the surface (0.4–0.8×10-6 m2 s−1). An increase in diffusivity with depth was observed that was generally larger in the soil–bedrock transition at 4–8 m depth. The outcomes are relevant for the understanding of soil thermodynamics in relation to other soil properties. Results with the spectral method suggest that changes in near-surface thermal diffusivity are related to changes in soil moisture content, which makes it a potential tool to gain information about soil drought and water resource availability from soil temperature data.