Advanced Search
Article Contents

Understanding the Soil Temperature Variability at Different Depths: Effects of Surface Air Temperature, Snow Cover, and the Soil Memory

Fund Project:

This work was sponsored by the National Key R&D Program of China (Grant No. 2016YFA0600404) and the National Natural Science Foundation of China (Grant Nos. 41530532 and 41675088). N. Y. also thanks the support from the Chinese Academy of Sciences Pioneer Hundred Talents Program


doi: 10.1007/s00376-020-0074-y

  • The soil temperature (ST) is closely related to the surface air temperature (AT), but their coupling may be affected by other factors. In this study, significant effects of the AT on the underlying ST were found, and the time taken to propagate downward to 320 cm can be up to 10 months. Besides the AT, the ST is also affected by memory effects—namely, its prior thermal conditions. At deeper depth (i.e., 320 cm), the effects of the AT from a particular season may be exceeded by the soil memory effects from the last season. At shallower layers (i.e., < 80 cm), the effects of the AT may be blocked by the snow cover, resulting in a poorly synchronous correlation between the AT and the ST. In northeastern China, this snow cover blockage mainly occurs in winter and then vanishes in the subsequent spring. Due to the thermal insulation effect of the snow cover, the winter ST at layers above 80 cm in northeastern China were found to continue to increase even during the recent global warming hiatus period. These findings may be instructive for better understanding ST variations, as well as land−atmosphere interactions.
    摘要: 土壤温度与地表气温密切相关,但两者之间的耦合会受到多种因素的影响。本文基于中国地表气温、积雪及多层土壤温度资料,利用线性相关及非线性因果分析方法研究了不同因素对土壤温度的影响。研究发现,地表气温可以显著的影响多层土壤温度,但随土壤深度的增加,地表气温的影响逐渐变弱且向下传递其影响所需的时间变长。地表气温的影响最长需要10个月才能传递到320厘米深的土壤。除了受地表气温的影响,土壤温度本身的记忆性特征也可以对未来的土壤温度产生显著的影响。在深层土壤(如320厘米),前一季度土壤温度记忆性对当季的影响甚至可以超过同期来自地表气温的影响。在浅层土壤(如80厘米以上),地表气温与土壤温度的耦合可能会被积雪阻挡。这种阻挡作用在我国东北地区的冬季尤其显著。受积雪的保温作用影响,即便是在本世纪初以来的全球变暖停滞期,我国东北地区冬季浅层土壤(80厘米以上)温度依旧保持升温趋势。这些研究结果有助于更好的理解多层土壤温度变率以及开展地-气相互作用的研究。
  • 加载中
  • Figure 1.  Cross mapping skills of the CCM between the AT and the ST at different layers (ST40, ST80, ST160, ST320) (left-hand column) and the cross mapping lag for the response of the STs to the variations in the AT (right-hand column; units: months). Cross mapping skills of the CCM at all stations are statistically significant at the 99% confidence level.

    Figure 2.  Cross mapping skills of the CCM between the STs and (i) the AT (black dashed line) and (ii) the STs in the previous season (blue dashed line). The right-hand column shows the results with all stations included, while the left-hand column excludes the stations in northeastern China. Above the red dashed line indicates the greater than 95% confidence level. MAM, March−May; JJA, June−August; SON, September−October; DJF, December−February.

    Figure 3.  Linear trends of the AT and STs (ST0, ST40, ST80, ST160, ST320) in winter (first column) and spring (second column) from 1998 to 2013. The spatial averages over northeastern China are shown in the third (winter) and fourth (spring) columns. “trend1” indicates the temperature changes from 1960 to 2013 and “trend2” shows the results from 1998 to 2013. “*” indicates the trends are statistically significant at the 95% confidence level. The curves are the time series of the anomalous AT and STs in northeastern China (right-hand vertical axis of the panels in the third and fourth columns; units: °C). The bars show the trends calculated from the sliding 15-year windows [left-hand vertical axis of panels in the third and fourth columns; units: °C (10 yr)−1].

    Figure 4.  The curves show the 15-year sliding correlation coefficients between the winter ST and four variables in northeastern China: the AT (first column), SSD (second column), SD (third column) and autumn ST (fourth column). Above the upper red line or below the bottom red line indicates the correlation coefficients are statistically significant at the 95% confidence level. The bars show the correlation coefficients from 1960 to 2013. The solid bars indicate the correlation coefficients are statistically significant at the 95% confidence level. All correlation coefficients were calculated after removing the trends of the time series.

  • Beltrami, H., and L. Kellman, 2003: An examination of short- and long-term air-ground temperature coupling. Global and Planetary Change, 38(3−4), 291−303, https://doi.org/10.1016/S0921-8181(03)00112-7.
    Beltrami, H., G. Ferguson, and R. N. Harris, 2005: Long-term tracking of climate change by underground temperatures. Geophys. Res. Lett., 32(19), L19707, https://doi.org/10.1029/2005GL023714.
    Chudinova, S. M., O. W. Frauenfeld, R. G. Barry, T. J. Zhang, and V. A. Sorokovikov, 2006: Relationship between air and soil temperature trends and periodicities in the permafrost regions of Russia. J. Geophys. Res., 111, F02008, https://doi.org/10.1029/2005JF000342.
    Cohen, J. L., J. C. Furtado, M. Barlow, V. A. Alexeev, and J. E. Cherry, 2012: Asymmetric seasonal temperature trends. Geophys. Res. Lett., 39, L04705, https://doi.org/10.1029/2011GL050582.
    Fan, X. G., 2009: Impacts of soil heating condition on precipitation simulations in the weather research and forecasting model. Mon. Wea. Rev., 137(7), 2263−2285, https://doi.org/10.1175/2009MWR2684.1.
    Fyfe, J. C., and Coauthors, 2016: Making sense of the early-2000s warming slowdown. Nature Climate Change, 6(3), 224−228, https://doi.org/10.1038/nclimate2938.
    Hu, Q., and S. Feng, 2004a: A role of the soil enthalpy in land memory. J. Climate, 17(18), 3633−3643, https://doi.org/10.1175/1520-0442(2004)017<3633:AROTSE>2.0.CO;2.
    Hu, Q., and S. Feng, 2004b: Why has the land memory changed? J. Climate, 17(16), 3236−3243, https://doi.org/10.1175/1520-0442(2004)017<3236:WHTLMC>2.0.CO;2.
    Hu, Q., and S. Feng, 2005: How have soil temperatures been affected by the surface temperature and precipitation in the Eurasian continent? Geophys. Res. Lett., 32, L14711, https://doi.org/10.1029/2005GL023469.
    Huang, Y., C. L. E. Franzke, N. M. Yuan, and Z. T. Fu, 2020: Systematic identification of causal relations in high-dimensional chaotic systems: Application to stratosphere-troposphere coupling. Climate Dyn., 55, 2469−2481, https://doi.org/10.1007/s00382-020-05394-0.
    Kosaka, Y., and S. P. Xie, 2013: Recent global-warming hiatus tied to equatorial Pacific surface cooling. Nature, 501(7467), 403−407, https://doi.org/10.1038/nature12534.
    Luo, D. H., Y. A. Chen, A. Dai, M. Mu, R. H. Zhang, and S. Ian, 2017: Winter Eurasian cooling linked with the Atlantic multidecadal oscillation. Environmental Research Letters, 12(12), 125002, https://doi.org/10.1088/1748-9326/aa8de8.
    Ma, H., S. Y. Leng, and L. N. Chen, 2018: Data-based prediction and causality inference of nonlinear dynamics. Science China Mathematics, 61(3), 403−420, https://doi.org/10.1007/s11425-017-9177-0.
    Mahanama, S. P. P., R. D. Koster, R. H. Reichle, and M. J. Suarez, 2008: Impact of subsurface temperature variability on surface air temperature variability: An AGCM study. Journal of Hydrometeorology, 9(4), 804−815, https://doi.org/10.1175/2008JHM949.1.
    McCracken, J. M., and R. S. Weigel, 2014: Convergent cross-mapping and pairwise asymmetric inference. Physical Review E, 90(6), 062903, https://doi.org/10.1103/PhysRevE.90.062903.
    Medhaug, I., M. B. Stolpe, E. M. Fischer, and R. Knutti, 2017: Reconciling controversies about the ‘global warming hiatus’. Nature, 545(7652), 41−47, https://doi.org/10.1038/nature22315.
    Mønster, D., R. Fusaroli, K. Tylén, A. Roepstorff, and J. F. Sherson, 2016: Inferring causality from noisy time series data-a test of convergent cross-mapping. Proc. 1st Int. Conf. on Complex Information Systems, Rome, Italy, 48−56, https://doi.org/10.5220/0005932600480056.
    Mønster, D., R. Fusaroli, K. Tylén, A. Roepstorff, and J. F. Sherson, 2017: Causal inference from noisy time-series data−Testing the convergent cross-mapping algorithm in the presence of noise and external influence. Future Generation Computer Systems, 73, 52−62, https://doi.org/10.1016/j.future.2016.12.009.
    Risbey, J. S., and S. Lewandowsky, 2017: Climate science: The 'pause' unpacked. Nature, 545(7652), 37−39, https://doi.org/10.1038/545037a.
    Sugihara, G., and R. M. May, 1990: Nonlinear forecasting as a way of distinguishing chaos from measurement error in time series. Nature, 344(6268), 734−741, https://doi.org/10.1038/344734a0.
    Sugihara, G., 1994: Nonlinear forecasting for the classification of natural time series. Philosophical Transactions of the Royal Society A: Mathematical Physical and Engineering Sciences, 348(1688), 477−495, https://doi.org/10.1098/rsta.1994.0106.
    Sugihara, G., R. May, H. Ye, C. H. Hsieh, E. Deyle, M. Fogarty, and S. Munch, 2012: Detecting causality in complex ecosystems. Science, 338(6106), 496−500, https://doi.org/10.1126/science.1227079.
    Takens, F., 1981: Detecting strange attractors in turbulence. Dynamical Systems and Turbulence, Warwick 1980, D. Rand and L. S. Young, Eds., Springer, 366−381, https://doi.org/10.1007/BFb0091924.
    Tang, M. C., and E. R. Reiter, 1986: The similarity between the maps of soil temperature and precipitation anomaly of the subsequent season. Plateau Meteorology, 5, 293−307. (in Chinese)
    Tang, M. C., J. X. Wang, and J. Zhang, 1987: A primary method for predicting the spring rainfall by the winter soil temperature depth 80cm. Plateau Meteorology, 6(3), 244−255. (in Chinese)
    Tang, M. C., S. H. Sun, Q. Zhong, and S. J. Wu, 1982: The energy variation of the underlying surface and the changes of the weather and climate. Plateau Meteorology, 1(1), 24−34. (in Chinese)
    Tang, Q. H., X. J. Zhang, X. H. Yang, and J. A. Francis, 2013: Cold winter extremes in northern continents linked to Arctic sea ice loss. Environmental Research Letters, 8(1), 014036, https://doi.org/10.1088/1748-9326/8/1/014036.
    Trenberth, K. E., J. T. Fasullo, G. Branstator, and A. S. Phillips, 2014: Seasonal aspects of the recent pause in surface warming. Nature Climate Change, 4(10), 911−916, https://doi.org/10.1038/NCLIMATE2341.
    Tsonis, A. A., E. R. Deyle, H. Ye, and G. Sugihara, 2018: Convergent cross mapping: Theory and an example. Advances in Nonlinear Geosciences, A. A. Tsonis, Ed., Springer, 587−600, https://doi.org/10.1007/978-3-319-58895-7_27.
    Wang, C. Y., S. P. Xie, Y. Kosaka, Q. Y. Liu, and X. T. Zheng, 2017: Global influence of tropical Pacific variability with implications for global warming slowdown. J. Climate, 30(7), 2679−2695, https://doi.org/10.1175/JCLI-D-15-0496.1.
    Wang, L., M. Henderson, B. H. Liu, X. J. Shen, X. W. Chen, L. Y. Lian, and D. W. Zhou, 2018a: Maximum and minimum soil surface temperature trends over China, 1965−2014. J. Geophys. Res., 123(4), 2004−2016, https://doi.org/10.1002/2017JD027283.
    Wang, Y. Q., J. Yang, Y. N. Chen, P. De Maeyer, Z. Li, and W. L. Duan, 2018b: Detecting the causal effect of soil moisture on precipitation using convergent cross mapping. Scientific Reports, 8, 12171, https://doi.org/10.1038/s41598-018-30669-2.
    Wu, L. Y., and J. Y. Zhang, 2014: Strong subsurface soil temperature feedbacks on summer climate variability over the arid/semi-arid regions of East Asia. Atmospheric Science Letters, 15(4), 307−313, https://doi.org/10.1002/asl2.504.
    Xu, W. H., C. H. Sun, J. Q. Zuo, Z. G. Ma, W. J. Li, and S. Yang, 2019: Homogenization of monthly ground surface temperature in China during 1961−2016 and performances of GLDAS reanalysis products. J. Climate, 32(4), 1121−1135, https://doi.org/10.1175/JCLI-D-18-0275.1.
    Xue, Y. K., R. Vasic, Z. Janjic, Y. M. Liu, and P. C. Chu, 2012: The impact of spring subsurface soil temperature anomaly in the western U.S. on North American summer precipitation: A case study using regional climate model downscaling. J. Geophys. Res., 117, D11103, https://doi.org/10.1029/2012JD017692.
    Yang, K., and J. Y. Zhang, 2015: Spatiotemporal characteristics of soil temperature memory in China from observation. Theor. Appl. Climatol., 126(3−4), 739−749, https://doi.org/10.1007/s00704-015-1613-9.
    Ye, H., E. R. Deyle, L. J. Gilarranz, and G. Sugihara, 2015: Distinguishing time-delayed causal interactions using convergent cross mapping. Scientific Reports, 5, 14750, https://doi.org/10.1038/srep14750.
    Zhang, N. N., G. L. Wang, and A. A. Tsonis, 2019: Dynamical evidence for causality between Northern Hemisphere annular mode and winter surface air temperature over Northeast Asia. Climate Dyn., 52, 3175−3182, https://doi.org/10.1007/s00382-018-4317-x.
    Zhang, T. J., 2005: Influence of the seasonal snow cover on the ground thermal regime: An overview. Rev. Geophys., 43, RG4002, https://doi.org/10.1029/2004RG000157.
    Zhang, T., R. G. Barry, D. Gilichinsky, S. S. Bykhovets, V. A. Sorokovikov, and J. P. Ye, 2001: An amplified signal of climatic change in soil temperatures during the last century at Irkutsk, Russia. Climatic Change, 49(1−2), 41−76, https://doi.org/10.1023/A:1010790203146.
  • [1] LIU Yonghe, FENG Jinming, MA Zhuguo, 2014: An Analysis of Historical and Future Temperature Fluctuations over China Based on CMIP5 Simulations, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 457-467.  doi: 10.1007/s00376-013-3093-0
    [2] Changyu ZHAO, Haishan CHEN, Shanlei SUN, 2018: Evaluating the Capabilities of Soil Enthalpy, Soil Moisture and Soil Temperature in Predicting Seasonal Precipitation, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 445-456.  doi: 10.1007/s00376-017-7006-5
    [3] CHEN Wei, and LU Riyu, 2014: A Decadal Shift of Summer Surface Air Temperature over Northeast Asia around the Mid-1990s, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 735-742.  doi: 10.1007/s00376-013-3154-4
    [4] MAN Wenmin, and ZHOU Tianjun, 2014: Regional-scale Surface Air Temperature and East Asian Summer Monsoon Changes during the Last Millennium Simulated by the FGOALS-gl Climate System Model, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 765-778.  doi: 10.1007/s00376-013-3123-y
    [5] Siqiong LUO, Zihang CHEN, Jingyuan WANG, Tonghua WU, Yao XIAO, Yongping QIAO, 2024: Impact of Initial Soil Conditions on Soil Hydrothermal and Surface Energy Fluxes in the Permafrost Region of the Tibetan Plateau, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 717-736.  doi: 10.1007/s00376-023-3100-z
    [6] CHEN Wei, LU Riyu, 2014: The Interannual Variation in Monthly Temperature over Northeast China during Summer, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 515-524.  doi: 10.1007/s00376-013-3102-3
    [7] Jing Peng, Li Dan, xiba tang, 2023: Spatial variation in CO2 concentration improves simulated surface air temperature increase in the Northern Hemisphere, ADVANCES IN ATMOSPHERIC SCIENCES.  doi: 10.1007/s00376-023-3249-5
    [8] Li Guo ping, Lu Jinghua, Jin Bingling, Bu Nima, 2001: The Effects of Anomalous Snow Cover of the Tibetan Plateau on the Surface Heating, ADVANCES IN ATMOSPHERIC SCIENCES, 18, 1207-1214.  doi: 10.1007/s00376-001-0034-0
    [9] LIU Ge, WU Renguang, ZHANG Yuanzhi, and NAN Sulan, 2014: The Summer Snow Cover Anomaly over the Tibetan Plateau and Its Association with Simultaneous Precipitation over the Mei-yu-Baiu region, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 755-764.  doi: 10.1007/s00376-013-3183-z
    [10] LI Qingxiang, LIU Xiaoning, ZHANG Hongzheng, Thomas C. PETERSON, David R. EASTERLING, 2004: Detecting and Adjusting Temporal Inhomogeneity in Chinese Mean Surface Air Temperature Data, ADVANCES IN ATMOSPHERIC SCIENCES, 21, 260-268.  doi: 10.1007/BF02915712
    [11] Jeong-Hyeong LEE, Keon-Tae SOHN, 2007: Prediction of Monthly Mean Surface Air Temperature in a Region of China, ADVANCES IN ATMOSPHERIC SCIENCES, 24, 503-508.  doi: 10.1007/s00376-007-0503-1
    [12] Wei CHEN, Xiaowei HONG, Riyu LU, Aifen JIN, Shizhu JIN, Jae-Cheol NAM, Jin-Ho SHIN, Tae-Young GOO, Baek-Jo KIM, 2016: Variation in Summer Surface Air Temperature over Northeast Asia and Its Associated Circulation Anomalies, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 1-9.  doi: 10.1007/s00376-015-5056-0
    [13] Chaofan LI, Riyu LU, Philip E. BETT, Adam A. SCAIFE, Nicola MARTIN, 2018: Skillful Seasonal Forecasts of Summer Surface Air Temperature in Western China by Global Seasonal Forecast System Version 5, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 955-964.  doi: 10.1007/s00376-018-7291-7
    [14] Wu Renguang, Chen Lieting, 1995: Interannual Fluctuations of Surface Air Temperature over North America and Its Relationship to the North Pacific SST Anomaly, ADVANCES IN ATMOSPHERIC SCIENCES, 12, 20-28.  doi: 10.1007/BF02661284
    [15] WANG Jia, ZHI Xiefei, and CHEN Yuwen, 2013: Probabilistic multimodel ensemble prediction of decadal variability of East Asian surface air temperature based on IPCC-AR5 near-term climate simulations, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1129-1142.  doi: 10.1007/s00376-012-2182-9
    [16] Shangfeng CHEN, Linye SONG, 2018: Impact of the Winter North Pacific Oscillation on the Surface Air Temperature over Eurasia and North America: Sensitivity to the Index Definition, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 702-712.  doi: 10.1007/s00376-017-7111-5
    [17] YUAN Wei, SUN Jianqi, 2009: Enhancement of the Summer North Atlantic Oscillation Influence on Northern Hemisphere Air Temperature, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 1209-1214.  doi: 10.1007/s00376-009-8148-x
    [18] Qian Weihong, Zhang Henian, Zhu Yafen, Dong-Kyou Lee, 2001: lnterannual and lnterdecadal Variability of East Asian Acas and Their Impact on Temperature of China in Winter Season for the Last Century, ADVANCES IN ATMOSPHERIC SCIENCES, 18, 511-523.  doi: 10.1007/s00376-001-0041-1
    [19] TIAN Di, GUO Yan*, DONG Wenjie, 2015: Future Changes and Uncertainties in Temperature and Precipitation over China Based on CMIP5 Models, ADVANCES IN ATMOSPHERIC SCIENCES, 32, 487-496.  doi: 10.1007/s00376-014-4102-7
    [20] Fuqiang YANG, Li DAN, Jing PENG, Xiujing YANG, Yueyue LI, Dongdong GAO, 2019: Subdaily to Seasonal Change of Surface Energy and Water Flux of the Haihe River Basin in China: Noah and Noah-MP Assessment, ADVANCES IN ATMOSPHERIC SCIENCES, 36, 79-92.  doi: 10.1007/s00376-018-8035-4
  • ESM-200074.pdf

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 18 March 2020
Manuscript revised: 20 September 2020
Manuscript accepted: 24 September 2020
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Understanding the Soil Temperature Variability at Different Depths: Effects of Surface Air Temperature, Snow Cover, and the Soil Memory

    Corresponding author: Naiming YUAN, naimingyuan@hotmail.com
  • 1. Key Laboratory of Regional−Climate Environment Research for Temperate East Asia, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
  • 2. University of Chinese Academy of Sciences, Beijing 100049, China
  • 3. Laboratory for Climate and Ocean−Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871, China

Abstract: The soil temperature (ST) is closely related to the surface air temperature (AT), but their coupling may be affected by other factors. In this study, significant effects of the AT on the underlying ST were found, and the time taken to propagate downward to 320 cm can be up to 10 months. Besides the AT, the ST is also affected by memory effects—namely, its prior thermal conditions. At deeper depth (i.e., 320 cm), the effects of the AT from a particular season may be exceeded by the soil memory effects from the last season. At shallower layers (i.e., < 80 cm), the effects of the AT may be blocked by the snow cover, resulting in a poorly synchronous correlation between the AT and the ST. In northeastern China, this snow cover blockage mainly occurs in winter and then vanishes in the subsequent spring. Due to the thermal insulation effect of the snow cover, the winter ST at layers above 80 cm in northeastern China were found to continue to increase even during the recent global warming hiatus period. These findings may be instructive for better understanding ST variations, as well as land−atmosphere interactions.

摘要: 土壤温度与地表气温密切相关,但两者之间的耦合会受到多种因素的影响。本文基于中国地表气温、积雪及多层土壤温度资料,利用线性相关及非线性因果分析方法研究了不同因素对土壤温度的影响。研究发现,地表气温可以显著的影响多层土壤温度,但随土壤深度的增加,地表气温的影响逐渐变弱且向下传递其影响所需的时间变长。地表气温的影响最长需要10个月才能传递到320厘米深的土壤。除了受地表气温的影响,土壤温度本身的记忆性特征也可以对未来的土壤温度产生显著的影响。在深层土壤(如320厘米),前一季度土壤温度记忆性对当季的影响甚至可以超过同期来自地表气温的影响。在浅层土壤(如80厘米以上),地表气温与土壤温度的耦合可能会被积雪阻挡。这种阻挡作用在我国东北地区的冬季尤其显著。受积雪的保温作用影响,即便是在本世纪初以来的全球变暖停滞期,我国东北地区冬季浅层土壤(80厘米以上)温度依旧保持升温趋势。这些研究结果有助于更好的理解多层土壤温度变率以及开展地-气相互作用的研究。

1.   Introduction
  • Soil temperature (ST) is a primary component of the thermal regimes of land and plays an important role in land−atmosphere interaction. Similar to the oceans, land areas can record their early anomalous conditions, which subsequently affect the heat exchange between the land surface and the overlying atmosphere, making ST an important climate indicator at seasonal time scales (Hu and Feng, 2004a, b). As early as in the 1980s, the winter ST was used to predict the precipitation during the flood season (Tang et al., 1982, 1987; Tang and Reiter, 1986). From then on, understanding the variability of ST has received increasing attention (Mahanama et al., 2008; Fan, 2009; Xue et al., 2012; Wu and Zhang, 2014; Yang and Zhang, 2015). Many different factors, such as the variations of solar radiation, precipitation, vegetation, and surface air temperature (AT), have been reported to be vital in modifying ST (Zhang et al., 2001; Beltrami et al., 2005).

    It is well recognized that the ST is normally closely coupled with the AT, and their temporal variations often exhibit similar characteristics (Beltrami and Kellman, 2003; Hu and Feng, 2005; Chudinova et al., 2006). This similarity, however, largely depends on the magnitude of the soil heat flux and the surface sensible heat flux, which may be affected by meteorological factors and physical properties of the underlying surfaces, such as vegetation cover, snow cover, etc. As a result, divergent changes between the AT and the ST may arise when the heat conduction between the near-surface atmosphere and the soil becomes relative weak or disappears (Wang et al., 2018a). Furthermore, in the vertical downward direction till deeper soil depths, the effects of the AT on the ST decreases, and the deeper ST may deviate from the shallower one due to the stronger persistence effects accumulated from preceding conditions (Yang and Zhang, 2015). Accordingly, although the AT plays an important role in modulating the variability of the ST, the detailed processes are complex, and sometimes controlled by multiple factors.

    In recent years, many researchers have evidenced a slowdown in the global mean surface temperature during the first decade of the 21st century, known as the global warming hiatus (Fyfe et al., 2016; Medhaug et al., 2017; Risbey and Lewandowsky, 2017). Modulation of the decadal variability of ocean heat seems to be the explanation behind this phenomenon (Kosaka and Xie, 2013; Trenberth at al., 2014; Wang at al., 2017). One of the strongest pauses occurred over northern Eurasia in winter, and it seems to be associated with the large-scale wintertime cooling over Eurasia that started in the late 1980s (Cohen et al., 2012; Trenberth at al., 2014). As the largest continent, Eurasia is a broad interface where the heat and water vapor transfer between land and atmosphere, but the role of the land in this hiatus is not clear. Did the ST also experience a warming hiatus? If yes, were there any regional differences? More generally speaking, what are the physical processes responsible for modulating the ST? How deep can the ST be significantly affected by the signals from the AT? What are the roles of other factors, such as the snow cover or the ST memory (i.e., the interseasonal persistence of the anomalous ST that reflects the effects of the ST from previous seasons on the current ST variations)? These questions have not been fully addressed owing to the lack of observational ST records.

    In this study, we collected in-situ ST records at multiple depths over China for a relatively long and continuous period (1960−2013). To investigate the relations between the ST records with other factors, such as AT, snow cover etc., both linear (time-lagged) correlation analysis and nonlinear causality analysis (convergent cross mapping and its time-lagged version) were conducted. It was found that the AT may significantly affect the ST at all the different layers considered in this study, but at the deepest depth (320 cm) the accumulated memory effects from the previous seasons may take over and become a more important factor. Affected by the AT, the warming hiatus was revealed in the winter ST over most regions of China during the recent global warming hiatus period. In the northeast, however, the winter ST at shallower layers (< 80 cm) continued to increase, which may be attributable to the thermal insulation effect of the snow cover. These findings may be helpful for a better understanding of ST variabilities, as well as land−air interactions.

    The rest of the paper is organized as follows: In section 2, we briefly introduce the data and methods used in the study. Causal relationships between the AT and the ST are reported in section 3, with the memory effects from previous seasons also considered. After a comparison of recent trends (1998−2013) in the winter AT and ST, the effects of snow cover and soil memory in modulating the ST are demonstrated in section 4. Finally, in section 5, we discuss the results and conclude the paper.

2.   Data and methods
  • The data used in this study were obtained from the China Meteorological Administration. They include the monthly mean surface AT records (recorded at a height of 1.5 m), the monthly mean ST records at the five layers of 0 cm, 40 cm, 80 cm, 160 cm and 320 cm (ST0, ST40, ST80, ST160, ST320, respectively), the sunshine duration (SSD), and the snow depth (SD). All the data cover the period from 1960 to 2013.

    Due to the serious lack of long-term ST data, we barely found enough stations with all meteorological elements recorded. Hence, we used a criterion that, for each element, only those stations with no more than 5.5% missing data during 1960−2013 were selected. The numbers of selected stations for the AT, ST0, ST40, ST80, ST160, ST320, SD and SSD were 2038, 1549, 239, 150, 23, 19, 550 and 1932, respectively. Northeastern China was defined as the region covering (40°−55°N, 122°−135°E). The number of stations with ST data decreased sharply from ST0 to ST320, but there were still eight stations for ST320 in northeastern China. Before analysis, the gaps in the time series caused by the missing data were filled using a simple linear interpolation algorithm.

  • Pearson correlation coefficients were calculated to identify the linear relations between different variables. The temporal changes of the correlations between different pairs of variables were also revealed by linear correlation in a sliding 15-year window, and the time steps of each sliding were 1 year. In addition, within each 15-year window, least-squares linear trends were also fitted (denoted as 15-year trends) to detect turning points of the trends.

  • We used convergent cross mapping (CCM) and its time-lagged version (time-lagged CCM) to identify variables that cause the variations in ST. This analysis method is effective in detecting whether variables are coupled in a nonlinear system, and their causality (Sugihara et al., 2012; Wang et al., 2018b; Zhang et al., 2019). Moreover, it can reveal the time delay if the variables are delay-coupled. Brief introductions to CCM and time-lagged CCM are given in the following two subsections.

  • Suppose we have two observational variables from a complex system, denoted as X and Y. CCM considers that if Y causes X, then the information of Y can be encoded in the nonlinear system attractor of X. Thus, by using an algorithm called cross mapping, the attractor of X is able to reconstruct the time series of Y (Sugihara et al., 2012). The cross mapping from X to Y can be carried out according to the following steps (Tsonis et al., 2018):

    (i) Based on Takens’ theorem (Takens, 1981), we first use the time series of X (with data length L) to reconstruct the nonlinear system attractor of X, which is represented by the vector ${{M}_X}({t_i}) = \left\{ {{x_{{t_i}}},{x_{{t_i} - {\tau _0}}}, \ldots,{x_{{t_i} - (E - 1){\tau _0}}}} \right\}$, where “ti” denotes any historical time point in the observations, E is the embedding dimension, and τ0 stands for the time delay. As we know, the attractor of X is acting as the dynamical state of the underlying nonlinear system, and it is encoded with much information on the temporal evolution of X. The two parameters E and τ0 can be determined through the simplex projection (Sugihara and May, 1990), more details of which are provided in the electronic supplementary material (ESM).

    (ii) Before using the attractor of $ {X} $ to estimate Y at time t, it is required to find the most similar dynamical states of $ {X} $ around the time t, and then Y at time t can be estimated through these historical similar states. Thus, the state of $ {X} $ at time t is denoted as vector ${{M}_X}(t)$ (t denotes the expected time in this cross mapping) and the state of $ {X} $ at any past time is denoted as vector ${{M}_X}({t_i})$ (${t_i}$ denotes any past time except t). The weight parameter ${w_i}$ that represents the associated weight between the states of $ {X} $ at time t and ${t_i}$ is calculated as:

    where ${\rm{d}}[{{M}_X}(t),{{M}_X}({t_i})]$ denotes the Euler distance between vectors ${{M}_X}(t)$ and ${{M}_X}({t_i})$. The most similar state to the time t will occupy the largest weight.

    (iii) According to Tsonis et al. (2018), Y at time t can be estimated by

    Then, we can calculate the cross mapping skill (Sugihara et al., 2012) that quantifies how close the estimated value $\widehat Y\left(t \right)$ is to the real value Y(t). Sugihara et al. (2012) and Tsonis et al. (2018) suggested to use the Pearson correlation ${\rho _{XY}}$ between Y(t) and $\widehat Y\left(t \right)$ to quantify the cross mapping skill, as shown in Eq. (4):

    Likewise, the cross mapping skill from Y to X can also be measured referring to the above steps, labeled as ${\rho _{YX}}$. ${\rho _{XY}}$ represents how well the attractor of X can be used to reconstruct the time series of Y, and if ${\rho _{XY}}$ is of high magnitude, then it means that the information of Y is encoded in the attractor of X through the causal influence from Y to X (Sugihara et al., 2012). Sugihara et al. (2012) and Tsonis et al. (2018) defined the causal inference from ${\rho _{XY}}$ and ${\rho _{YX}}$ as follows: (i) if ${\rho _{XY}}$ is convergent when L is increased, and ${\rho _{XY}}$ is statistically significant in the surrogate test, then Y is suggested to be a causation of X. (ii) Moreover, if ${\rho _{YX}}$ is also convergent when L is increased, and statistically significant in the surrogate test, then the causal relationship between X and Y is bidirectional (X and Y cause each other).

    CCM is an effective method for detecting the causality in a coupled nonlinear system; however, one should note that there are limitations, especially under the circumstance of strong coupling. As previously reported (Sugihara et al., 2012; McCracken and Weigel, 2014; Mønster et al., 2016, 2017; Ma et al., 2018), the non-zero cross mapping skill may not reflect the real coupling direction in the presence of strong coupling or even intermediate coupling. Using CCM, it may not be easy to distinguish the real causality relations from the symmetric coupling. Furthermore, the existence of noise may also reduce the cross mapping skill. Therefore, the detected causality from the CCM needs to be considered with caution, and its time-lagged version, the time-lagged CCM, is sometimes required.

  • Time-lagged CCM performs better in detecting causality relations in the presence of a strong coupling effect and random noise and can further estimate the delay effect of causal interaction, according to Ye et al. (2015). When detecting a causal influence from Y(t + tp) to X(t) (tp denotes the time lag and accepts both 0, positive and negative values) in a delayed-coupled system, it is assumed that the information of Y(t + tp) has been encoded to X(t), so that the value of Y(t + tp) can be reconstructed from X(t) (Takens, 1981; Ye et al., 2015). Following the computation algorithm of the CCM, one can use X(t) to cross map Y(t + tp) with different values of tp, and estimate the cross mapping skill ${\rho _{XY}}$(tp). Previous studies (Ye et al., 2015; Wang et al., 2018b) suggest that an optimal lag will exist where the cross mapping skill reaches its maximum. This optimal lag is determined by the delay time of the dynamical system (Ye et al., 2015), and we can use this optimal lag to identify the causality direction and the dynamical delay time between two coupled variables.

    As suggested in Ye et al. (2015), the causal relationship between two variables can be identified according to the following new criterion of causal inference: (i) if the maximum of ${\rho _{XY}}$(tp) occurs at a negative lag (i.e., tp < 0), this means that the information of Y at the past time is encoded in the current state of X. In other words, the current state of X can be used to reconstruct Y at the past time. In this case, Y causes X; (ii) if the maximum of ${\rho _{XY}}$(tp) occurs at a positive lag, this means that the information of Y at the future time is determined by the current state of X. In this case, the current state of X can be used to reconstruct Y at the future time, and X causes Y. Based on this new criterion, the time-lagged CCM can produce robust results of the causality directions and avoid the weakness of the CCM. Therefore, compared with the causal inference of CCM (without time lag), this new causal inference from time-lagged CCM analysis has been found to be more accurate in the presence of a strong coupling effect and random noise (Ye et al., 2015; Huang et al., 2020); see also Figs. S1−S3 in the ESM. In this study, we apply the time-lagged CCM to study the causal relations between the AT and the ST at multiple layers.

3.   Effects of AT on the underlying ST of different layers
  • The AT variability can propagate downward to the deep soil layers (Beltrami and Kellman, 2003; Hu and Feng, 2005; Chudinova et al., 2006), but the signal decays during the propagation, and the significance of the AT effects on the ST in deep layers is doubtful. Compared to the time-lagged linear correlation analysis (Fig. S4 in the ESM), time-lagged CCM can distinguish the true causal relationship from the phenomenon of synchrony even if the variables are nonlinearly coupled (Sugihara et al., 2012; Ye et al., 2015; Wang et al., 2018b). Before the analysis of the causality from the AT to the STs, we first tested whether nonlinear processes were reflected in the time series of the variables of interest. As shown in Fig. S5 and Fig. S6 in the ESM, using simplex projection (Sugihara and May, 1990) the optimal embedding dimension (E) was calculated. Then, using the S-map method (Sugihara, 1994) based on this E, we found the time series of both the AT and the STs have nonlinear signals, suggesting that nonlinear causal inference analysis (CCM or its time-lagged version) is necessary. Accordingly, based on the monthly anomaly time series (data length: 648 months), time-lagged CCM was thus employed to detect the strength (the maximum of the cross mapping skills for different time lags in time-lagged CCM) and the time lag of the impacts of the monthly AT on the ST at different layers. As shown in Fig. 1, the AT impacts on the ST are regionally different. With the increase in depth, the impact intensity decays rapidly, but the cross mapping skills of the CCM at all layers are statistically significant at the 99% confidence level (first column in Fig. 1). It takes time for the signals from the AT to propagate downward to the soil layers. At the depth of 40 cm, the responses of the ST to the AT variability at most stations are either “synchronous” or have one-month lags (right-hand column in Fig. 1). Only at several stations in northeastern China and northern Xinjiang province is the time taken to propagate downward to 40 cm longer than one month. It is worth noting that the detected “synchronization” is actually limited by the temporal resolution (i.e., monthly) of the data, as the time-lagged CCM cannot detect time lags shorter than one month. For more precise information, datasets with higher temporal resolution (e.g., daily) are required. With the depth increasing, the time lags increase from around one month at 40 cm to even ten months at some stations at 320 cm. Note that the AT and ST are weakly coupled in northeastern China. The cross mapping skills in this region are obviously smaller than in other regions. This may be a result of the fact that the influences of the AT take more time to propagate downward to the deep soil layers, and thus the intensity decays more. However, the mechanisms for this phenomenon are still unclear. It may be related to the soil type, but further studies are needed.

    Figure 1.  Cross mapping skills of the CCM between the AT and the ST at different layers (ST40, ST80, ST160, ST320) (left-hand column) and the cross mapping lag for the response of the STs to the variations in the AT (right-hand column; units: months). Cross mapping skills of the CCM at all stations are statistically significant at the 99% confidence level.

  • Since the impacts of AT decay with the increase in depth, one may ask whether its impacts still dominate the variations of the ST at deeper layers. Particularly, in view of the strong soil memory in deep layers, is there a depth where the effects of soil memory become more important? To answer this question, the impact intensities of the AT and the soil memory on the ST in different seasons were compared in terms of the CCM coefficient ${\rho _{XY}}$, where Y represents the ST in the preceding seasons or the AT, while X represents the ST in the considered season. Note that, due to the data resolution (seasonal data, for example, spring series from 1960 to 2013) required in this analysis, we could not apply time-lagged CCM here. Instead, we employed CCM and only compare the AT effects of a considered season with the soil memory effects accumulated from one season before. Although the CCM coefficients may be affected by noise or the strength of the coupling, since the true causality from the AT to the STs (AT effects) has been detected by the time-lagged CCM (Fig. 1) and the pairwise asymmetric inference (PAI) method (Fig. S7 in the ESM) suggested by McCracken and Weigel (2014), to some extent the CCM coefficients may reflect the relative intensities of the AT effects versus the soil memory effects. Figure 2 shows that, in all seasons, the AT effects are always statistically significant at different layers, while the soil memory effects become consistently significant only from deep layers (≥ 160 cm) (right-hand column in Fig. 2). With the increase in depth, the AT effects decrease but the soil memory effects increase. At the depth of 160 cm, one can already see that effects of the soil memory can exceed those of the AT in summer (June−August), autumn (September−November) and winter (December−February). With further depth, i.e., at 320 cm, the soil memory effects eventually exceed the AT effects in all seasons and become more important. It is worth noting that the soil memory effects in northeastern China may be stronger, because when excluding the stations in northeastern China (left-hand column in Fig. 2), the soil memory effects become weaker. At the 160 cm layer, the soil memory effects are weaker than the AT effects (except for the effects in autumn), and only at the 320 cm layer are the soil memory effects significant in all seasons.

    Figure 2.  Cross mapping skills of the CCM between the STs and (i) the AT (black dashed line) and (ii) the STs in the previous season (blue dashed line). The right-hand column shows the results with all stations included, while the left-hand column excludes the stations in northeastern China. Above the red dashed line indicates the greater than 95% confidence level. MAM, March−May; JJA, June−August; SON, September−October; DJF, December−February.

4.   Seasonal AT and ST trends during the recent global warming hiatus
  • During the first decade of the 21st century there was a warming slowdown in the global mean surface temperature. When different seasons are considered, however, the hiatus did not occur everywhere. In the last section, significant AT impacts on the subsurface ST were found at all stations and all soil layers. Therefore, the ST trends are supposed to be similar to the AT trends, regardless of warming or cooling. As shown in Fig. S8 in the ESM, the spatial distribution of the AT and ST trends are similar in spring, summer and autumn. As for winter, however, the cooling trends of the AT in northeastern China have not been replicated by the STs in this region. On the contrary, continuous warming trends in the STs from 0 cm to 320 cm can be seen during the warming hiatus period (Fig. 3).

    Figure 3.  Linear trends of the AT and STs (ST0, ST40, ST80, ST160, ST320) in winter (first column) and spring (second column) from 1998 to 2013. The spatial averages over northeastern China are shown in the third (winter) and fourth (spring) columns. “trend1” indicates the temperature changes from 1960 to 2013 and “trend2” shows the results from 1998 to 2013. “*” indicates the trends are statistically significant at the 95% confidence level. The curves are the time series of the anomalous AT and STs in northeastern China (right-hand vertical axis of the panels in the third and fourth columns; units: °C). The bars show the trends calculated from the sliding 15-year windows [left-hand vertical axis of panels in the third and fourth columns; units: °C (10 yr)−1].

    For a closer look of the trends, the third column of Fig. 3 shows the time series of the winter AT and winter ST at different depths averaged over northeastern China. It is clear that they all have similar trends from 1960, but since 1998, the trends between the AT and the ST display opposite signs (Fig. 3). A dramatic increase in ST0 is apparent from 1998 to 2013, and this is mainly attributable to the change in the observational infrastructure in 2005. Before 2005, the winter surface ST was usually measured manually on the bare ground. After the use of automatic measurement instruments in 2005, however, the ST0 was measured automatically even when the ground was snow covered. Thus, higher temperatures may be measured after 2005 (Xu et al., 2019). In spite of this, however, the ST0 still increases after 2005, and this warming cannot be explained by the change in the observational infrastructure. With increased depth, the ST trend becomes weaker. For the ST40, ST80, ST160 and ST320 layers, the linear trends from 1998 to 2013 are 0.544°C (10 yr)−1, 0.472°C (10 yr)−1, 0.071°C (10 yr)−1 and 0.15°C (10 yr)−1, respectively. To identify the start time of the warming trend, 15-year trends of the ST were calculated. It was found that the trends of ST40 and ST80 experienced a negative-to-positive change at the beginning of the 21st century, while those of ST160 and ST320 have been positive since the mid-1970s. This implies that the warming of the soil below 160 cm from 1998 to 2013 may be part of a persistent warming at a longer time scale. In other words, the physical processes related to the wintertime ST variability seem to vary with depth.

    When spring comes, the situation becomes completely different. The trend deviations between the AT and the ST almost disappear in northeastern China. Except for the ST at 320 cm, decreasing trends in ST are found at most stations in northeastern China, as well as in the spatial averages during the warming hiatus period (Fig. 3). The 15-year trends further reveal a positive-to-negative change at the beginning of the 21st century.

  • To understand the observed warming trend of the ST at different layers in winter, potential factors including the SD (the SD exceeds 4 cm in northeastern China, as shown in Fig. S9 in the ESM), the AT, the autumn ST, and the SSD were first studied using linear correlation analysis (Fig. 4). Considering the inhomogeneity of the ST0 (see previous section), in this section we only focus on the subsurface ST from 40 cm to 320 cm. The weak coupling between the AT and the ST in northeastern China was reconfirmed by the nonsignificant linear correlation coefficients (first column in Fig. 4). The SSD was found to be closely related with the ST40, but the correlation is negative, reflecting indirect relations between the SSD and the ST40. The effects of AT and SSD seem to be blocked by the effects of the snow cover, as significant positive correlations between the ST and the SD were found at the depths of 40 cm and 80 cm. The 15-year sliding correlations further indicate that the positive relations between the SD and the ST were stable in the past four decades and became more significant after the mid-1990s. This finding may be related to the insulation effect of the thicker snow cover. As reported in Zhang (2005), an increase in SD by 5 to 15 cm can lead to a 1°C increase in ground temperature. In the past few decades, large-scale wintertime cooling and increased snowfall over Eurasia have been found (Cohen et al., 2012; Tang et al., 2013; Luo et al., 2017). In this context, the SD in northeastern China has become thicker, and the spatial pattern of the SD increasing trend (Fig. S9) was found to be similar to the ST trends at the layers of 40 cm and 80 cm. Moreover, the 15-year trends of the SD in northeastern China also reveal a negative-to-positive change at the beginning of the 21st century (figure not shown). These facts support the inference that the thicker snow cover may contribute to the warmer ST at these two layers.

    Figure 4.  The curves show the 15-year sliding correlation coefficients between the winter ST and four variables in northeastern China: the AT (first column), SSD (second column), SD (third column) and autumn ST (fourth column). Above the upper red line or below the bottom red line indicates the correlation coefficients are statistically significant at the 95% confidence level. The bars show the correlation coefficients from 1960 to 2013. The solid bars indicate the correlation coefficients are statistically significant at the 95% confidence level. All correlation coefficients were calculated after removing the trends of the time series.

    With the increase in depth, the correlations between the ST and the SD become nonsignificant, but the ST memory gradually starts to dominate. At the depth of 160 cm, the ST in autumn was found to be closely related to the ST in winter, but the 15-year sliding correlations reveal significant correlations only before 1980. For the ST at 320 cm, the memory from the previous season seems to have bigger impacts, as the correlation coefficients between the ST in autumn and the ST in winter are consistently significant over the past few decades. Accordingly, the memory effects from autumn may be important for the winter ST at deep layers. It is worth noting, however, that the persistent impact from autumn seems to mainly contribute to the interannual variations of the winter ST. It may not be possible to explain the slightly increasing trends of the winter ST at 160 cm and 320 cm by the memory effect from autumn, as the autumn ST at 160 cm and 320 cm does not exhibit increasing trends in northeastern China during the warming hiatus (Fig. S8).

    Due to the time series of the SD (seasonal data) having nonlinear signals (Fig. S10 in the ESM), to further confirm the effects of the SD and the soil memory on the variations in ST at different layers, in the rest of this section the CCM method was applied (Fig. S11 in the ESM). In line with the results of linear correlation analysis (Fig. 4), the causal effects of the SD on the ST40 and ST80 in winter are remarkable, with the cross mapping skills nearly converging to 0.4 and 0.3, respectively. As for the ST at deeper layers, however, the cross mapping skills between the SD and the ST160 (ST320) almost vanish, indicating no causal effects. To rule out the potential biases induced by the drawbacks of the CCM (see Methods section), we further applied the PAI method (McCracken and Weigel, 2014) to check the causal relations detected by the CCM. As shown in Fig. S12 in the ESM, the causal effects of the SD on the winter ST at shallower layers (i.e., 40 cm and 80 cm) are confirmed well, without doubt. Regarding the effect of the soil memory on the winter ST, the CCM analysis indicates that the autumn ST has significant impacts on the winter ST at 320 cm, with the cross mapping skill nearly converging to 0.75. As for the ST at 160 cm, the CCM analysis does not infer a certain causality between the autumn ST and the winter ST, as the cross mapping skill is not convergent, indicating this depth might be a transition layer where the ST is influenced by the combined effects of both the snow cover and the soil memory.

5.   Conclusions
  • In this study, responses of ST to different factors, including the surface AT, the snow cover, as well as the memory effects, were analyzed using ST data of different depths (0 cm to 320 cm) over China. Beyond previous studies that were mainly based on linear correlation analyses, this work investigated their relations using the CCM method and its time-lagged version, which can reveal the causality relations between the factors and the ST. As one of the most important factors, the AT was found to have significant impacts on the ST at all the considered depths, but the influence becomes weaker with the increase in depth. The memory effects, on the contrary, become more important with increased depth. At the depth of 320 cm, the thermal conditions from previous seasons were found to have stronger impacts than those from the AT of the considered season, indicating a possible change in the main controlling factor of the subsurface ST.

    At shallower layers, the effects of AT on the ST are strong, but they may still be blocked by the change in land cover, e.g., the snow cover, in which case it is not easy for the signals from the AT to propagate downward into the soil and the corresponding effects are weakened. For instance, in northeastern China, the wintertime AT exhibited a decreasing trend during the recent global warming hiatus period, while the underlying ST showed an increasing trend. This unexpected increase was attributed to the increasing SD, which prevented the loss of heat from the soil to the atmosphere in winter. A detailed CCM analysis further revealed that the thermal insulation effect of the snow cover significantly affected the ST till the layer at 80 cm. It is worth noting that the snow cover in northeastern China becomes very thin in spring, and the thermal isolation effect very weak. The springtime ST showed a decreasing trend from 1998 to 2013, which was consistent with the AT.

    These findings describe a physical picture of how the ST at different depths is influenced by different factors. They may be helpful for better understanding ST variabilities and land−atmosphere interactions. However, there are still some issues to be discussed. Firstly, limited by the observational ST dataset, we only focused here on stations in China. Whether the conclusions drawn in this study are applicable in other regions needs further study in the future. Secondly, one should note that the CCM analysis only revealed the causality relations from the AT to the ST. In the opposite direction, the influences of the ST on the AT were not detected. Recently, the mystery of the warming hiatus has been unpacked, and the heat redistribution processes in the oceans have been reported as a key reason. The fact that more energy is transferred and restored in the deep ocean is suggested as the main explanation for the slowdown in the warming trend. Similarly, did land-related processes also contribute to the warming hiatus? More specifically, since the ST has been considered as an important predictor for future climate change on seasonal time scales, does the ST affect the AT? And if so, how? To address these questions and further improve our understanding of land−atmosphere interactions, further research efforts in the future remain an urgent necessity.

    Acknowledgements. This work was sponsored by the National Key R&D Program of China (Grant No. 2016YFA0600404) and the National Natural Science Foundation of China (Grant Nos. 41530532 and 41675088). N. Y. also thanks the support from the Chinese Academy of Sciences Pioneer Hundred Talents Program.

    Electronic supplementary material: Supplementary material is available in the online version of this article at https://doi.org/10.1007/s00376-020-0074-y.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return