Search

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 x 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 x 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.

The anthropogenically intensified greenhouse effect has caused a radiative imbalance at the top of the atmosphere during the industrial period. This, in turn, has led to an energy surplus in various components of the Earth system, with the ocean storing the largest part. The land contribution ranks second with the latest observational estimates based on borehole temperature profiles, which quantify the terrestrial energy surplus to be 6 % in the last 5 decades, whereas studies based on state-of-the-art climate models scale it down to 2 %. This underestimation stems from land surface models (LSMs) having a subsurface that is too shallow, which severely constrains the land heat uptake simulated by Earth system models (ESMs). A forced simulation of the last 2000 years with the Max Planck Institute ESM (MPI-ESM) using a deep LSM captures 4 times more heat than the standard shallow MPI-ESM simulations in the historical period, well above the estimates provided by other ESMs. However, deepening the LSM does not remarkably affect the simulated surface temperature. It is shown that the heat stored during the historical period by an ESM using a deep LSM component can be accurately estimated by considering the surface temperatures simulated by the ESM using a shallow LSM and propagating them with a standalone forward model. This result is used to derive estimates of land heat uptake using all available observational datasets, reanalysis products, and state-of-the-art ESM experiments. This approach yields values of 10.5–16.0 ZJ for 1971–2018, which are 12 %–42 % smaller than the latest borehole-based estimates (18.2 ZJ).

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 x 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 x 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.

Understanding the near-surface soil thermal regime and its connection to the atmospheric state is important for the assessment of several climate-related processes. However, the lack of in situ soil temperatures measurements limits the analysis of such processes. In this study, we have developed a quality-controlled soil temperature database for Spain that consists of 39 sites spanning from 1987 to 2018. We have used this database to assess the near-surface soil thermal regime. Likewise, we evaluate at seasonal to multidecadal timescales the land–air temperature coupling over Spain by analysing the structure of the surface air temperature (SAT) and the ground surface temperature (GST) covariance and also their long-term evolution. In addition, we have employed the ERA5-Land reanalysis to test the consistence between observations and reanalysis. The results show that the near-surface soil thermal structure is dominated by conduction despite some influence of hydrology-related processes. Regarding the land–air temperature coupling, we have found a strong connection between SAT and GST. However, in the summer months there is an offset in SAT–GST at some sites due to limited evaporation and enhanced sensible heat fluxes. Furthermore, multidecadal SAT–GST decoupling may exist over some sites as a response to decreasing precipitation. The ERA5-Land represents the observations' climatology well, but it underestimates the summer soil temperature observations and the long-term trends at some sites.

This work was supported by the projects IlModelS, project no. CGL2014-726 59644-R and GReatModelS, project no. RTI2018-102305-B-C21. The work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under project ID bm1026. Vladimir Alexeev was supported by the Interdisciplinary Research for Arctic Coastal Environments (InteRFACE) project through the Department of Energy, Office of Science, Biological and Environmental Research Program's Regional and Global Model Analysis program, and by NOAA project NA18OAR4590417. We also wish to thank Veronika Gayler for technical support on JSBACH and Christian Reick for helpful comments and discussion., Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. This study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models. Plain Language Summary Many current-generation climate models have land components that are too shallow. Under climate change conditions, the long-term warming trend at the surface propagates deeper into the ground than the commonly used 3-10 m. Shallow models alter the terrestrial heat storage and distribution of temperatures in the subsurface, influencing the simulated land-atmosphere interactions. Previous studies focusing on annual timescales suggest that deeper models are required to match subsurface-temperature observations and the classic analytical heat conduction solution. However, for a systematic investigation of land-model deepening in the frame of anthropogenic climate change, the classic analytical solution is inaccurate because it does not mimic the timescale and amplitude of the simulated warming trend. This study intends to bridge the gap between analytical and simulatio, Ministerio de Economía y Competitividad (MINECO), Scientific Steering Committee (WLA), Interdisciplinary Research for Arctic Coastal Environments (InteRFACE) project through the Department of Energy, Office of Science, Biological and Environmental Research Program's Regional and Global Model Analysis program, NOAA project, Depto. de Física de la Tierra y Astrofísica, Fac. de Ciencias Físicas, TRUE, pub

© 2021 American Meteorological Society. L Thanks to R. Schnur and V. Gayler from MPI Hamburg. We gratefully acknowledge the IlModels (CGL2014-59644-R) and GreatModelS (RTI2018102305-B-C21 and RTI2018-102305-A-C22) projects funded by the Spanish MINECO. SH contributed in the frame of the ERANET-plus-Russia project SODEEP (Study Of the Development of Extreme Events over Permafrost areas) supported by BMBF (Grant 01DJ18016A). This work used resources of the Deutsches Klimarechenzentrum (DKRZ) granted by its Scientific Steering Committee (WLA) under project ID bm1026., The representation of the thermal and hydrological states in land surface models is important for a realistic simulation of land-atmosphere coupling processes. The available evidence indicates that the simulation of subsurface thermodynamics in Earth system models is inaccurate due to a zero-heat-flux bottom boundary condition being imposed too close to the surface. To assess the influence of soil model depth on the simulated terrestrial energy and subsurface thermal state, sensitivity experiments have been carried out in piControl, historical, and RCP scenarios. A deeper bottom boundary condition placement has been introduced into the JSBACH land surface model by enlarging the vertical stratification from 5 to 12 layers, thereby expanding its depth from 9.83 to 1416.84 m. The model takes several hundred years to reach an equilibrium state in stand-alone piControl simulations. A depth of 100 m is necessary, and 300 m recommendable, to handle the warming trends in historical and scenario simulations. Using a deep bottom boundary, warming of the soil column is reduced by 0.5 to 1.5 K in scenario simulations over most land areas, with the largest changes occurring in northern high latitudes, consistent with polar amplification. Energy storage is 3-5 times larger in the deep than in the shallow model and increases progressively with additional soil layers until the model depth reaches about 200 m. While the contents of Part I focus on the sensitivity of subsurface thermodynamics to enlarging the space for energy, Part II addresses the sensitivity to changing the space for water and improving hydrological and phase-change interactions., Ministerio de Economía y Competitividad (MINECO), BMBF Federal Ministry of Education & Research, Scientific Steering Committee (WLA), Depto. de Física de la Tierra y Astrofísica, Fac. de Ciencias Físicas, TRUE, pub