Changes in the freeze-thaw cycles of shallow soil have important consequences for surface and subsurface hydrology, land-atmosphere energy and moisture interaction, carbon exchange, and ecosystem diversity and productivity. This work examines the shallow soil freeze-thaw cycle on the Tibetan Plateau (TP) using in-situ soil temperature observations in 0-20 cm soil layer during July 1982-June 2017. The domain and layer averaged beginning frozen day is November 18 and delays by 2.2 days per decade; the ending frozen day is March 9 and advances by 3.2 days per decade; the number of frozen days is 109 and shortens by 5.2 days per decade. Altitude and latitude combined could explain the spatial patterns of annual mean freeze-thaw status well. Stations located near 0 degrees C contour line experienced dramatic changes in freeze-thaw cycles as seen from subtropical mountain coniferous forest in the southern TP. Soil completely freezes from surface to 20-cm depth in 15 days while completely thaws in 10 days on average. Near-surface soil displays more pronounced changes than deeper soil. Surface air temperature strongly influences the shallow soil freeze-thaw status but snow exerts limited effects. Different thresholds in freeze-thaw status definition lead to differences in the shallow soil freeze-thaw status and multiple-consecutive-day approach appears to be more robust and reliable. Gridded soil temperature products could resolve the spatial pattern of the observed shallow soil freeze-thaw status to some extent but further improvement is needed.

A change in soil temperature (ST) is a significant indicator of climate change, so understanding the variations in ST is required for studying the changes of the Qinghai-Tibet Plateau (QTP) permafrost. We investigated the performance of three reanalysis ST products at three soil depths (0-10 cm, 10-40 cm, and 40-100 cm) on the permafrost regions of the QTP: the European Centre for Medium-Range Weather Forecasts interim reanalysis (ERA-Interim), the second version of the National Centers for Environmental Prediction Climate Forecast System (CFSv2), and the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). Our results indicate that all three reanalysis ST products underestimate observations with negative mean bias error values at all three soil layers. The MERRA-2 product performed best in the first and second soil layers, and the ERA-Interim product performed best in the third soil layer. The spatiotemporal changes of annual and seasonal STs on the QTP from 1980 to 2017 were investigated using Sen's slope estimator and the Mann-Kendall test. There was an increasing trend of ST in the deeper soil layer, which was less than that of the shallow soil layers in the spring and summer as well as annually. In contrast, the first-layer ST warming rate was significantly lower than that of the deeper soil layers in the autumn and winter. The significantly (P < 0.01) increasing trend of the annual ST indicates that the QTP has experienced climate warming during the past 38 years, which is one of the factors promoting permafrost degradation of the QTP.