Next Article in Journal
Dictionary Learning- and Total Variation-Based High-Light-Efficiency Snapshot Multi-Aperture Spectral Imaging
Next Article in Special Issue
Spatial Trend and Impact of Snowmelt Rate in Spring across China’s Three Main Stable Snow Cover Regions over the Past 40 Years Based on Remote Sensing
Previous Article in Journal
Assessing Spatiotemporal Dynamics of Land Use and Cover Change and Carbon Storage in China’s Ecological Conservation Pilot Zone: A Case Study in Fujian Province
Previous Article in Special Issue
Forest Fire Effects on Landscape Snow Albedo Recovery and Decay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating the Effects of Snow Cover and Vegetation on Soil Temperature Using Remote Sensing Indicators in the Three River Source Region, China

1
Key Laboratory of Land Surface Process and Climate Change in Cold and Arid Regions, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Qinghai Climate Center, Xining 810001, China
4
Key Laboratory of Remote Sensing of Gansu Province, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 4114; https://doi.org/10.3390/rs14164114
Submission received: 20 July 2022 / Revised: 18 August 2022 / Accepted: 19 August 2022 / Published: 22 August 2022
(This article belongs to the Special Issue Remote Sensing for Mountain Vegetation and Snow Cover)

Abstract

:
Soil temperature is an important physical variable that characterises geothermal conditions and influences geophysical, biological and chemical processes in the earth sciences. Soil temperature is not only affected by climatic and geographical factors; it is also modulated by local factors such as snow cover and vegetation. This paper investigates the relationship between snow cover and vegetation and soil temperature with the help of two classical remote sensing indicators, the Snow Cover Days (SCD) based Advanced Very High Resolution Radiometer and the Normalized Difference Vegetation Index (NDVI)-based Global Inventory Modelling and Mapping Studies, to analyse the influence of local factors on soil temperature in the Three River Source Region (TRSR). Combing multi-layer geothermal observations from 23 stations in the TRSR with meteorological dataset, soil properties datasets, snow cover and vegetation indices, a non-linear model, the Random Forest model, is used to establish a multi-layer soil temperature dataset to analyse the influence of surface cover factors in each depth. The results showed that the annual SCD had a decreasing trend during 1982–2015 and was negatively correlated with the annual mean soil temperature; the annual NDVI had no significant trend, but it was positively correlated with the annual mean soil temperature. Regionally, there was a significant decrease in SCD in the mountainous areas bordering the source areas of the three rivers, and there was a trend of increasing NDVI in the northwest and decreasing vegetation in the southwest in the TRSR. The stronger the correlation with soil temperature in areas with a larger SCD, the more the snow has a cooling effect on the shallower soil temperatures due to the high albedo of the accumulated snow and the repeated melting and heat absorption of the snow in the area. The snow has an insulating effect on the 40 cm soil layer by impeding the cooling effect of the atmosphere in winter. In sparsely vegetated areas, vegetation lowers ground albedo and warms the soil, but in July and August, in areas with more vegetation, NDVI is negatively correlated with soil temperature, with heavy vegetation intercepting summer radiant energy and having a cooling effect on the soil.

Graphical Abstract

1. Introduction

The Three River Source Region (TRSR), located at the headwaters of three important rivers in China, is an area sensitive to climate change. The warming trend in this region is more than twice the global average warming rate and much higher than the global average for the same latitudes [1,2,3]. In addition, precipitation in the TRSR has also shown a significant increase, and there has been a tendency for increased precipitation to become more pronounced since the beginning of the 21st century [1,4,5]. As a result of climate change, soil temperatures in this region have warmed to varying degrees in all layers [6,7], and the stability of permafrost and the freeze-thaw process of seasonal permafrost have changed. Studies have shown that the area of permafrost is shrinking, the thickness of the active layer is increasing, and the temperature of the active layer is increasing; the maximum freeze depth of seasonal permafrost is decreasing, the freeze start date is delayed, the thaw end date is earlier, and the freeze-thaw duration is shorter [6,8,9,10].
There are a variety of factors that contribute to the changes in soil temperature and permafrost. At the global scale, the soil temperature is controlled by latitudinal zonation patterns; at the regional scale, it is influenced by elevation and climatic wetness, while at the local scale, it is influenced by the surface cover, soil properties and topographic slope orientation. Climatic and geographic factors have received more scholarly attention in previous studies [10,11,12,13,14]. However, as a mediator between the lithosphere and the atmosphere, the influence of vegetation and snow cover on soil temperature should not be neglected, even within a small area with different superficial thermal conditions of snow cover and vegetation [15,16,17].
Snow cover, with properties such as high albedo and low thermal conductivity, alters the surface radiation balance and energy distribution processes, leading to differences in soil temperature and soil freeze-thaw processes [18,19]. Increasing snow depth keeps soil temperatures relatively stable at higher levels and increases the temperature difference between the ground and the atmosphere, while snow cover weakens the correlation between meteorological factors and soil temperatures [20,21]. The effect on soil temperature varies according to the characteristics of the snow cover, which can have both a cooling and a holding effect. Thicker snow cover can prevent the soil from radiating long-wave radiation to the atmosphere [15,22]. Conversely, in areas with low snowfall where the snowpack is thin and of short duration, the melting of snow absorbs heat and sublimation dissipates it, leading to a drop in soil temperature [19].
Similarly, differences in vegetation cover have a significant impact on annual differences in soil heat flux, soil freeze-thaw and surface temperature [16,23]. The effect of vegetation on soil temperature is highly complex and varies across seasons, types and percentages of cover. In warm seasons, vegetation can prevent some heat from entering the soil and reduce land surface temperatures; in cool seasons, vegetation can reduce heat loss from the soil and help maintain soil temperatures [7,24,25,26]. Different vegetation types have different abilities to reflect solar radiation and root insulation, which have significant effects on mean monthly surface temperatures and soil freeze-thaw processes [27]. It has been found that as the height and density of vegetation decreases, the mean daily minimum soil temperature increases significantly and the timing of the start and end of freeze-thaw and the duration of freeze-thaw change [28,29].
Observational studies have been carried out in and around the TRSR. In Li′s study, observations from three stations were used to investigate the effects of different snow cover conditions on soil water heat exchange processes [19]. Luo et al. used site monitoring and vegetation survey data as a basis to analyse changes in temperature and moisture in the active layer at three sites in the Yellow River source area with similar local conditions but significant differences in vegetation [28]. Jin′s study was conducted along the highway, and the results showed that the effect of snow and vegetation on ground temperature was twofold. Ground temperature increases gradually as vegetation cover decreases, but decreases when vegetation cover is reduced to 0–20%. The insulation effect is generally weak when the snow cover is thin for a short period of time, and begins to increase when the snow cover exceeds 20 cm in thickness [15].
However, in situ observations are subject to high spatial limitations. Compared to field observations, satellite-based remote sensing is an important tool for detecting large-scale vegetation and snow cover dynamics in space and time [24]. To better explore the influence of snow cover and vegetation on soil temperature over the region, this study uses machine learning methods to calculate multi-layer soil temperature data in the TRSR region with the help of snow cover and vegetation remote sensing indicators, combined with multi-layer geothermal observations.

2. Materials and Methods

2.1. Description of the Study Area

The study area is TRSR and surrounding areas, with a coordinating range of 89°24′–103°30′E, 31°–37°12′N (Figure 1). The TRSR is located in the interior of the Tibet Plateau (TP), China, and contains the headwaters of the Yellow, Yangtze and Mekong River, also known as the ‘Chinese water tower’. The TRSR is a sensitive area for climate change, where the temperature increase trend far exceeds the global average warming, and also significantly exceeds the average for the same latitudes and the whole of China. In addition, the TRSR is also at the transition zone between permafrost and seasonally frozen ground, and the permafrost degradation caused by climate change is more apparent in this region.

2.2. Data

2.2.1. Site Data

In this paper, observation data of 23 sites during 1960–2020 (Figure 1) are collected as training data and validation data to derive the soil temperature regression model. The observations used in this study were obtained from the Climate Center of Qinghai Province, China, and included precipitation, wind speed, and soil temperature at 0 cm, 5 cm, 10 cm, 20 cm, 40 cm, 80 cm, 160 cm, and 320 cm depths. The air temperature data for these stations are obtained from the China Meteorological Data website (https://data.cma.cn/data/cdcdetail/dataCode/SURF_CLI_CHN_MUL_MMON_19712000.html (accessed on 6 July 2021)), and the original data are daily data, which are calculated to obtain monthly averages. The elevation of these 23 sites is distributed between 2491.4 m and 4612.2 m, including 4 sites at 2000–3000 m, 11 sites at 3000–4000 m, and 8 sites at 4000–5000 m.

2.2.2. Reanalysis Data

Monthly air temperature, precipitation, and wind speed data at 0.1° × 0.1° spatial resolution for the period 1982–2015 were obtained from the China Meteorological Forcing Dataset (CMFD) (http://poles.tpdc.ac.cn/zh-hans/data/8028b944-daaa-4511-8769-965612652c49/ (accessed on 30 July 2020)), a high spatial and temporal resolution gridded near-surface meteorological dataset developed by the Institute of Tibetan Plateau Research, Chinese Academy of Sciences (ITP/CAS), specifically for the study of land surface processes in China [30,31,32]. The CMFD was made through the fusion of remote sensing products, a reanalysis dataset and in situ observation data at weather stations.

2.2.3. Remote Sensing Data

The monthly snow cover days (SCD) and the number of days with snow on the ground were calculated from the China 5 km Cloud-gap-filled Advanced Very High Resolution Radiometer Snow Cover Extent product dataset (CGF-AVHRR SCE) which was provided by the National Cryosphere Desert Data Center (http://www.ncdc.ac.cn (accessed on 6 April 2022)) [33]. The CGF-AVHRR SCE product is prepared for the snow distribution area in China, based on the AVHRR CDR SR product, using a multi-level decision tree classification method combined with a vacancy-filling algorithm such as the Spatio-temporal interpolation algorithm of the Hidden Markov Random Field Model, removing cloud effects from products.
The Normalized Difference Vegetation Index (NDVI) was used to characterize the land cover condition. The NDVI product of the National Oceanic and Atmospheric Administration (NOAA) Global Inventory Monitoring and Modelling System (GIMMS NDVI3g v1) was selected since it is the classic dataset used to detect vegetation change [34]. The GIMMS NDVI data were downloaded from National Tibetan Plateau Data Center (https://data.tpdc.ac.cn/en/data/9775f2b4-7370-4e5e-a537-3482c9a83d88/ (accessed on 9 November 2020)), and the temporal resolution is twice a month, while the spatial resolution is 1/12 of a degree. The monthly data are calculated by the maximum value composite method.

2.2.4. Auxiliary Data

Digital elevation model (DEM) data with a spatial resolution of 0.008° × 0.008° was downloaded from the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn (accessed on 24 September 2020)).
Soil properties data were obtained from the Dataset of Soil Properties for Land Surface Modelling over China (http://globalchange.bnu.edu.cn/research/soil2 (accessed on 24 July 2020)) [35], and six variables were selected, bulk density (BD), sand content (SA), silt content (SI), clay content (CL), gravel content (GRAV), and soil organic matter (SOM). The resolution is 30 arc-seconds (about 1 km at the equator). The vertical variation of soil property was captured by eight layers to the depth of 2.3 m (i.e., 0–0.045, 0.045–0.091, 0.091–0.166, 0.166–0.289, 0.289–0.493, 0.493–0.829, 0.829–1.383 and 1.383–2.296 m). The data used for the computations were calculated by weighted averaging to a specific soil depth close to the observation level.
In this study, the period from 1982 to 2015 was chosen as the study period, and all the grid data were resampled and interpolated to 0.01°×0.01° by bilinear interpolation. For variables lacking observations, such as SCD and NDVI, site values are obtained by interpolating on the corresponding grid data.

2.3. Method

2.3.1. Model Evaluation

Four metrics were used to assess the performance of the model’s simulation of soil temperature: root-mean-square error (RMSE), bias, Pearson correlation coefficient (r) and the model determination coefficient (R2).
Suppose n as the sample size, yi as the ith true value and y ι ^ as the corresponding simulated value. Then, RMSE is calculated as follows:
R M S E = 1 n i = 1 n y i y ι ^ 2  
Bias is calculated as follows:
b i a s = E y ι ^ y i  
R2 is calculated as follows:
R 2 = 1 i = i n y i y ι ^ 2 i = 1 n y i y ¯ 2  
where y ¯ is expected value of y.

2.3.2. Correlation Analysis

The relationship between the ground temperature and its influencing factors was analysed by Pearson′s correlation method. The Pearson′s correlation coefficient (r) is calculated as follows:
r = C o v x , y s x s y  
where the numerator is the sample covariance of the variables x and y and the denominator is the product of the standard deviations of the two variables.

2.3.3. Trend Analysis

The Theil–Sen Median is a robust non-parametric statistics method applied to fit a straight line that estimates the overall trend by searching for the median slope of connecting lines between every two data points [36]. For asymmetrical distribution or heteroscedastic data, the accuracy of the method is much superior to that of the simple linear regression, while for normally distributed data it has comparable statistical efficacy to the non-robust model. Sen′s slope is calculated as follows:
s l o p e = M e d i a n x i x j i j f o r     1 j < i n
for a data set containing n data, there are n n 1 2 data pairs.
The Theil–Sen estimator is generally combined with the Mann–Kendall method to test trends in the data. The test statistic for the MK method is:
S = i = 1 n 1 j = i + 1 n s g n x i + 1 x i = i < j s g n x j x i  
s g n Δ x =       1 ,     i f   Δ x > 0       0 ,     i f   Δ x = 0 1 ,     i f   Δ x < 0
Then, if the sample size n >10, the variance of S is:
V a r S = 1 18 n n 1 2 n + 5 j = 1 J t i t j 1 2 t j + 5  
where j is the number of groups of repeated values and tj is the number of repeated values in the j-th group.
However, for sequences with autocorrelation, the modified Mann–Kendall (MMK) method is more commonly used [36]. The MMK method takes into account the effect of autocorrelation on the variance of the Mann–Kendall trend test statistic. The MMK method modifies the variance calculation method for S based on the MK method as follows:
V a r * S = η × V a r S  
η is the variance correction factor, which is calculated as follows:
η = 1 + 2 n n 1 n 2 × i = 1 n 1 n i n i 1 n i 2 r i  
where r(i) is the autocorrelation function of the rank of the observations:
r i = k = 1 n i ( R k R ) R k + 1 R k = 1 n R k R 2  
Ri is the rank of yi, and R is the average order of yi:
y i = x i β × i  
β = M e d i a n x j x i j i , f o r   1 j < i N
where β is the trend estimator based on the rank of the series. The corrected test statistic Z is as follows:
Z = S 1 V a r * S ,     i f   S > 0           0       ,     i f   S = 0 S + 1 V a r * S ,     i f   S < 0 r i = k = 1 n i ( R k R ) R k + 1 R k = 1 n R k R 2  
Positive values of Z indicate increasing trends while negative values show decreasing trends. Absolute values of Z greater than or equal to 1.645, 1.96 and 2.576 indicate that they pass the significance test at 90%, 95% and 99% confidence levels respectively.

2.3.4. Regression Method

The study compared three different methods for the computation of regression models, respectively Multiple Linear Regression (MLR) in statistical methods, Support Vector Machine (SVM) in machine learning, and Random Forests (RF) in ensemble learning.
MLR is the most common method to investigate the effect of multiple characters on a physical variable, assuming a linear relationship between soil temperature and multiple influencing factors. For a multiple linear regression equation with K arguments, the equation is:
y = β 0 + β 1 x 1 + β 2 x 2 + β 3 x 3 + + β i x i + + β K x K + μ  
where β0 is the regression constant, βi is the partial regression coefficient of the i-th variable and μ is the residual. The least-squares estimation of the parameters is used to evaluate the parameter βi.
Soil temperature is influenced by many factors and often requires the use of high-dimensional data. SVM has advantages in computing high-dimensional data by mapping high-dimensional data to a low-dimensional space to find the regression hyperplane [37]. SVM allows for some deviation between predicted and true values, which will reduce the effect of anomalous data on the model and reduce overfitting and make the model more generalizable [38,39]. The choice of different nonlinear kernel functions within the SVM model can capture different nonlinear relationships; the Radial Basis Function was chosen in the study.
The process of training multiple machine learning models and combining their outputs is called ensemble learning. There are two different types of integrated learning: bagging and boosting. RF is one of the bagging methods, which builds multiple independent simple models by repeated random sampling, and each mini-model can consider the influence of impact factors on soil temperature from a different perspective [40]. Thus, RF prevents overfitting and makes the model more generalizable. RF also has the characteristics of fast learning speed and insensitivity to noise, which can achieve good results without extensive parameter adjustment [41].

2.3.5. Data Pre-Processing

Periodic variation between soil temperature and surface factors are not always synchronous, especially in deep layers, and there is a lag correlation between them. The response of deeper soils to upper layer changes has a sluggish response [6,42]. To investigate the lagging relationship between soil temperature and surface factors at different depths, Pearson correlation coefficients r were calculated for soil temperature and surface factors in the current month, one month earlier, two months earlier, three months earlier, and four months earlier.
The results indicate that the lagged correlations of soil temperature with air temperature, precipitation, snow cover and vegetation were all significant, and the correlations with wind speed were poor. Soil temperature at 0–40 cm depth had the highest correlation with surface factors in the current month; soil temperature at 80 cm depth had the highest correlation with air temperature, precipitation, and snow cover one month ago, and higher correlation with vegetation both in the current month and one month ago; soil temperature at 160 cm had the highest correlation with air temperature two months ago, high correlation with precipitation both in the current month and one month ago, and highest correlation with vegetation and snow cover one month ago; soil temperature at 320 cm had a higher correlation with air temperature, precipitation, and snow cover two months ago and three months ago, and higher correlation with vegetation two months before ago.
Due to the limited effect of wind speed on soil temperature [12], wind speed is not involved in the calculation of soil temperature below 40 cm depth. Considering the lagged response of each layer, the surface influence factor of the current month was used for the regression of 0–40 cm soil temperature, and the lags of 1, 2 and 3 months were set for 80 cm, 160 cm and 320 cm soil temperature, respectively.
Additionally, based on the findings of previous studies, the inclusion of a periodicity factor can increase the accuracy of the model when modelling periodic changes in soil temperature. Therefore, this study adds the month in which each data sample is located to the model to indicate the phase in which the data is situated in the annual cycle variation.
Of the 23 site observations, 6 sites were retained as validation sites for the final regional soil temperature dataset that were not involved in the model building process. The remaining 17 sites were used for the model building process, for which 70% of the data was applied as training data for model training and parameter tuning by the five-fold cross validation method, and the rest as testing data for testing the regression models. Soil temperature at 0 cm was observed during 1960–2020 and soil temperature at 5–320 cm was observed during 1960–2017, and the amount of data used for model building in these nine layers were 6491, 4985, 4985, 4984, 4983, 2722, 2717, 2729, and 2726, respectively.
This study is based on python software for calculations; the scikit-learn package is used to calculate statistical model and machine learning models and the pymannkendall package is used for trend calculations and MK tests [43].

3. Results

3.1. Models Selection

In this study, three regression models (MLR, SVM, RF) were chosen to evaluate the simulation ability for soil temperature in each layer by r, RMSE, and bias (Table 1). The R2 of all three regression models were almost greater than 0.95 for the testing data at 0–40 cm, the RMSEs were above 1 °C for both MLR and SVM, and the RMSEs were below 1 °C for RF. At 80 cm, 160 cm and 320 cm depth, the R2 of the MLR simulation results decreased sharply with depth. The R2 of the SVM model is around 0.95, while that of the RF model is above 0.97.
The observation–simulation plots for each layer (Figure 2) show that MLR and SVM are poorly simulated for low temperatures at 0–40 cm depth, with MLR being under-simulated for low temperatures and SVM being over-simulated for low temperatures, with individual data simulated significantly under-simulated in the higher temperature of the SVM.

3.2. Result Validation

Based on collected grid data, RF was chosen to calculate the soil temperature at 0–320 cm in TRSR, and a 34-year regional soil temperature dataset from 1982 to 2015 was obtained.
Site verification was carried out on the TRSR regional soil temperature dataset. The six sites used for verification were Qumalai, Gande, Yushu, Jiuzhi, Tongde and Tonren, with their elevations of 4175, 4050, 3716.9, 3628.5, 3148.2 and 2491.4, all in meters, respectively.
Due to the difference between the site data and the grid data, elevation correction for soil temperature in each layer is required when validating the region dataset with the site data. Linear regression using multi-year mean ground temperature versus elevation at 23 stations calculated the soil temperature lapse rate for each layer to be 0.43, 0.38, 0.38, 0.38, 0.39, 0.43, 0.42, 0.41, 0.41 all in °C/100m. The simulated data were corrected based on the lapse rate and compared with the observation data (Figure 3). The validation results show that the R2 between simulated and observed values are above 0.92 at 0–20 cm, and decrease with depth at 40 cm, 80 cm, 160 cm and 320 cm, with 0.86, 0.72, 0.70 and 0.50 respectively.

3.3. Results Analysis

3.3.1. Interannual Variation

Based on the TRSR soil temperature dataset, the anomalies of regional annual soil temperature, NDVI and SCE from 1983 to 2015 were calculated. What can be clearly seen in Figure 4a,b is the gradual rise in soil temperature in each layer; from a negative to a positive distance level after 1998, the regional average slopes for the nine layers were 0.49, 0.39, 0.36, 0.36, 0.33, 0.27, 0.24, 0.16 and 0.08, all in °C/decade, respectively. There is a decreasing trend in SCD, with an average decrease of 9.58 days per decade. Before 2000, the SCD was mostly at a positive anomaly, and after 2000 it was almost at a negative anomaly. Alpine meadows and alpine steppes are the main vegetation types in TRSR [44]. The NDVI in the region is low overall, and the MK test results showed no significant trend in NDVI (Figure 4c,d).
On an interannual scale, the variation in soil temperature was significantly negatively correlated with SCD, with the highest correlation at 0 cm with a correlation coefficient of −0.7, and the correlation coefficient between soil temperature and SCD decreased with increasing depth. The correlation between soil temperature and NDVI was significantly positive and decreased with depth, with correlation coefficients above 0.5 at 0–40 cm and above 0.4 at 80–320 cm. The correlation between soil temperature and vegetation was relatively small compared to the correlation between soil temperature and SCD.
As the most essential land cover of the TRSR during the cold season (from November to March) and warm seasons (from May to September) respectively, snow and vegetation affect the energy exchange processes between air and soil, and thus soil temperature. To explore the effect of snow and vegetation on the energy transfer between the soil and atmosphere, the temperature difference between soil temperature and air temperature was calculated for each layer. The results show that the difference between soil and air temperatures tends to decrease significantly in all layers except 0 cm during the study period, decreasing by 0.18, 0.20, 0.21, 0.23, 0.30, 0.32, 0.41 and 0.48 °C, all in °C/decade, in 5–320 cm, respectively. The deeper the depth, the more the difference in ground temperature decreases. The difference in ground temperature at 0 cm is much smaller than at other levels of the soil and its interannual fluctuations are less pronounced than at other levels.
A decrease in the number of snow days always corresponds to a decrease in the difference in ground temperature, and a decrease in NDVI always corresponds to an increase in ground temperature difference (Figure 4c,d). The temperature difference in the soil and atmosphere is positively correlated with the SCD with a correlation coefficient of around 0.7, and negatively correlated with the NDVI value with a correlation coefficient of around −0.6.

3.3.2. Seasonal Changes

(1)
Regional Average
From November to April, the SCD (Figure 5a) is greater than 5 days in all months except June to September, with January and February averaging greater than 10 days. The NDVI (Figure 5b) is greater than 0.2 from May to October, reaching a maximum in July and August.
In the annual variation, only the situation at the depth of 0–40 cm is calculated for the monthly mean ground temperature difference because of the phase difference between the soil temperature below 40 cm and the shallow climate factor.
The temperature difference at 0 cm is smaller in the cold season and larger in the warm season. The difference in ground temperature at 0 cm is more variable in the cold season and decreases with the number of snow days. Overall, however, the annual variation in the 0 cm ground temperature difference is small, ranging from 3.35 to 5.24 °C.
The depth of the ground temperature difference at 5–40 cm depth has greater intra-annual amplitude, ranging from 4.31–7.70 °C, 3.37–9.36 °C and 2.42–11.20 °C for 5 cm, 20 cm and 40 cm depths respectively. It is smaller in the warm season, larger in the cold season, positively correlated with changes in the number of days of snowfall and negatively correlated with the vegetation index in the warm season, and, with increasing depth, the correlation decreases.
(2)
Spatial Distribution
Snow cover in the TRSR mainly occurs from October to March, with a regional average of 6 to 10 days per month (Figure 6a). The central part of the Three Rivers Source, especially the junction of the Yangtze River source and Lancang River source, the south-western border of the Lancang River source, and the central part of the Yellow River source and its south-western border are the main areas with large values of SCD, with a monthly average of more than 15 days in the cold season.
The regional trend in SCD shows (Figure 6b) a significant decrease in the number of snow days from January to March and in December, especially in the central areas of the Yangtze and Lancang River sources and the western and eastern areas of the Yellow River source, with January and December decreasing by more than 1.5 days per decade on average. In contrast, May, October and November saw an increase in the number of snow days in some areas.
The NDVI of the region (Figure 7a) decreases from southeast to northwest, with vegetation growing from May until October when it gradually degrades. The vegetation is lush in July, August and September, with a regional average NDVI for 34 years greater than 0.4. In July and August, the NDVI was above 0.8 in the eastern part of the Yellow River source. The western part of the Yangtze River source and the north-western part of the Yellow River source is sparsely vegetated, with an NDVI no greater than 0.4 throughout the year.
Overall, the NDVI index is increasing in the TRSR, especially from May to October, with a clear trend of increasing NDVI of 0.01 per decade or more in most areas. The trend of increasing NDVI was greater in May and June in the southeast and in July, August and September in the northwest.
As a whole, SCD was negatively correlated with soil temperature from 0–40 cm, with a stronger negative correlation for shallow soils (Figure 8). At 0 cm, the number of snow days in December and January had the greatest effect on soil temperature, with more snow days resulting in lower soil temperatures and snow acting as a cooling agent in this layer. However, at a depth of 40 cm, soil temperature and the SCD in the north-western part of the TRSR are positively correlated, with more snow days leading to a higher soil temperature and the insulating effect of snow on this layer.
Different from SCD, the NDVI index was generally positively correlated with soil temperature, with a significantly stronger positive correlation in the cold season than in the warm season (Figure 9). In July, August and September the NDVI showed a non-significant negative correlation with soil temperature in the southeast, with vegetation playing a cooling role. In November, December, January and February, NDVI and soil temperature in the western part of the Yangtze River source were significantly positively correlated, with lower NDVI correlating with lower shallow soil temperature; however, they were negatively correlated with soil temperature at 40 cm, with vegetation playing a cooling role.
At 0 cm (Figure 10a), during the warm season, the difference between ground temperature and air temperature does not differ significantly across the TRSR. Relatively speaking, the difference in ground temperature is slightly lower in the south-east. In the cold season, the temperature difference gradually increases in the north-western part of the TRSR, while the difference in ground temperature gradually decreases in other regions, with the overall temperature difference being smaller than in the warm season. In some areas in the east and south, the difference in ground temperature even decreases to negative values. Spatial variability in the regional distribution of winter temperature differences is high.
However, the opposite is true for soil layers below 0 cm, and at a depth of 20 cm (Figure 10b), the difference in ground temperature is significantly smaller in the warm season than in the cold season. At 40 cm soil, the difference in ground temperature is significantly higher in the cold season than in shallower soil layers, where the soil temperature is much warmer than the air temperature. The warm season ground temperature difference was negative in the southwest and east, which implies that the soil temperature was lower than the air temperature.

4. Discussion

4.1. Soil Temperature in Correlation with Snow Cover

The lower the number of annual snow days during the study period, the higher the average annual soil temperature and the cooling effect of snow. This is because the snow on the ground has a high albedo, which reduces net shortwave radiation, causing the net radiative flux to decrease, which in turn reduces sensible, latent and soil heat fluxes [19,21]. In addition, Li et al. point out that the repeated melting of snow in all seasons causes the snow to have a cooling effect on the soil. In a warming climate, a significant reduction in the SCD will reduce its cooling effect on the soil, which will make soil temperatures more sensitive to climate change [19]. The correlation coefficients between the ground temperature and the number of snow days in each layer show that the effect of snow on soil temperature diminishes with this depth, but can affect up to 320 cm on an annual scale.
The correlation coefficient between SCD and soil temperature is predominantly negative every month; the more days of snow cover, the lower the soil temperature, and the snow cover mainly plays a cooling role. For 0 cm soil temperature, the cooling effect of snow accumulation is strongest in January and December in the north-western part of the TRSR. This may be since snow accumulation increases the surface albedo, making the ground receive less short-wave radiation during the day. At the same time, the snow covering the soil repeatedly melts and absorbs heat, cooling the 0 cm soil temperature. At 40 cm, in areas with more snow in the north-western part of the Yangtze River source, the soil temperature is positively correlated with the number of snow days, with the snow acting as insulation for the deeper soil layers. This is mainly because in winter, the soil radiates long-wave radiation to the atmosphere and the presence of snow impedes the long-wave energy from the soil to the atmosphere, making the soil temperature warmer at 40 cm depth in the snow-covered area. In contrast, in the mountainous areas of the central TRSR, where the snow is thicker, the insulating effect of the snow can be seen at 5–20 cm. The decreasing in SCD in the context of climate change reduces its cooling effect on the shallow ground temperature and contributes to the rise in shallow surface temperatures. As for soil temperatures in deep layers, the insulating effect of the snow cover is reduced in winter, making the energy exchange between the deep soil and the atmosphere more direct and the permafrost more susceptible to temperature changes.

4.2. Soil Temperature in Correlation with Vegetation

There was no clear trend in the annual mean NDVI variation in the study area, but soil temperatures were generally higher when the vegetation was in good condition, and there was a clear positive correlation between soil temperature and NDVI variation, with NDVI correlating better with shallow soil temperatures than with deeper layers. Previous studies have suggested that vegetation intercepts radiant energy reaching the surface, but this study found that vegetation may have a warming effect on soil. This is probably because the vegetation type of the TRSR is dominated by alpine grasslands and alpine meadows, where low vegetation has a limited role in blocking radiant energy, and its low albedo compared to bare soil raises the near surface temperature and leads to soil warming.
In winter, NDVI and soil temperature are positively correlated mainly at 0 cm. Especially in sparsely vegetated areas, such as the western part of the Yangtze River source, lower NDVI is associated with lower soil temperatures. This may be due to the fact that the less vegetation there is, the more the surface is exposed in winter when the wind is strong near the surface, the increased sensible heat exchange, which facilitates the outward emission of long-wave radiation from the soil and soil cooling. In addition, dry, bare soil has a higher albedo and the sparser the vegetation, the more radiant energy is reflected from the soil. In the south-eastern region, the vegetation has a more developed root system and the soil has a higher organic matter content, which favours soil insulation. In July and August, vegetation is abundant in the south-eastern part of the TRSR, with NDVI greater than 0.5. The vegetation played a cooling role, probably due to the shading effect of the lush vegetation on the high solar energy in summer and the reduction of evaporation of surface soil water, which slows down the summer soil warming.
There is a tendency for vegetation to increase during the warm season in the north-western part of the TRSR, which will lead to an increased cooling effect of vegetation in summer, slowing down soil warming in summer; while the accumulation of soil organic matter continues to increase the insulating effect of vegetation roots in winter.

5. Conclusions

From 1982 to 2015, SCD decreased significantly in the TRSR, and was negatively correlated with each soil temperature and positively correlated with the difference in ground temperature, with the correlation weakening with increasing depth. NDVI showed no significant trend during this period, but was positively correlated with soil temperature and negatively correlated with the difference in ground temperature.
Regionally, the confluence areas of the three river source areas and the north-western part of the three river sources are areas of large values of snow days, with the former showing a significant trend of decreasing snow days from December to March. Vegetation decreases from southeast to northwest, with NDVI below 0.3 throughout the year in northwest areas, which are sparsely vegetated, but vegetation in these areas has a significant increasing trend over 34 years, but not in densely vegetated areas.
In general, soil temperature is negatively correlated with SCD, snow has a cooling effect, soil temperature is positively correlated with NDVI and vegetation plays a warming role.
The negative correlation between soil temperature and SCD is greatest at 0 cm, where the high albedo of snow and the repeated melting of thinner thicknesses of snow cools the soil. However, at 40 cm, especially in the north-western part of the Yangtze River source, snow is positively correlated with soil temperature in winter and its insulating effect on the deeper soils, where the presence of snow slows down the cooling of the deeper soils.
NDVI is positively correlated with soil temperature in winter, and in less vegetated areas not only is albedo high, but high surface winds carry heat away from the bare soil, allowing it to cool. In areas with heavy vegetation in summer, the vegetation intercepts the strong radiant energy, which slows down soil warming.
In the background of climate change, a reduction in the number of snow days will reduce its cooling effect on shallow soils as well, but make deep soil temperatures more susceptible to fluctuations with temperature. In contrast, the trend towards increased vegetation in the north-west of the TRSR in summer will result in greater interception of radiation by vegetation in the area, as well as greater insulation of the soil in winter with the accumulation of organic matter in the root system.
In this study, only 23 sites were used to build the soil temperature model, and the sites were unevenly distributed, with more in the east and less in the west, making regional representation limited. In addition, the lack of site observations for snow and vegetation makes it difficult to verify the accuracy of remote sensing information in the TRSR area. The reliability of the study results will be further improved if the observation stations and observation variables can be supplemented.

Author Contributions

Conceptualization, X.T. and S.L.; methodology, X.T.; software, X.T.; analysis, X.T. and S.L.; investigation, X.T., J.W., Z.C. and Q.D.; data curation, X.T., H.L., X.H. and J.W.; writing—original draft preparation, X.T.; writing—review and editing, X.T., S.L., J.W. and Q.D.; visualization, X.T. and Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported jointly by the National Natural Science Foundation of China (U20A2081, 41975096, 41971325), and the West Light Foundation of the Chinese Academy of Sciences (xbzg-zdsys-202102).

Data Availability Statement

The air temperature site data can be obtained from the China Meteorological Data website (https://data.cma.cn (accessed on 6 July 2021)). The China Meteorological Forcing Dataset can be accessed online (http://poles.tpdc.ac.cn/zh-hans/data/8028b944-daaa-4511-8769-965612652c49/ (accessed on 30 July 2020)). The China 5 km Cloud-gap-filled AVHRR Snow Cover Extent product dataset from National Cryosphere Desert Data Center can be accessed online (http://www.ncdc.ac.cn (accessed on 6 April 2022)). The NDVI product of NOAA Global Inventory Monitoring and Modeling System (GIMMS NDVI3g v1) was obtained online (https://data.tpdc.ac.cn/en/data/9775f2b4-7370-4e5e-a537-3482c9a83d88/ (accessed on 9 November 2020)). The Dataset of Soil Properties for Land Surface Modeling over China can be accessed online (http://globalchange.bnu.edu.cn/research/soil2 (accessed on 24 July 2020)).

Acknowledgments

We thank the Qinghai Provincial Climate Center for providing the observed soil temperature data at the site.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liang, L.; Li, L.; Liu, C.; Cuo, L. Climate Change in the Tibetan Plateau Three Rivers Source Region: 1960–2009. Int. J. Climatol. 2013, 33, 2900–2916. [Google Scholar] [CrossRef]
  2. Bai, Y.; Guo, C.; Degen, A.A.; Ahmad, A.A.; Wang, W.; Zhang, T.; Li, W.; Ma, L.; Huang, M.; Zeng, H. Climate Warming Benefits Alpine Vegetation Growth in Three-River Headwater Region. China Sci. Total Environ. 2020, 742, 140574. [Google Scholar] [CrossRef] [PubMed]
  3. Jin, Z.; You, Q.; Wu, F.; Sun, B.; Cai, Z. Changes of climate and climate extremes in the Three-Rivers Headwaters’ Region over the Tibetan Plateau during the past 60 years. Trans. Atmos. Sci. 2020, 43, 1042–1055. (In Chinese) [Google Scholar]
  4. Piao, S.; Ciais, P.; Huang, Y.; Shen, Z.; Peng, S.; Li, J.; Zhou, L.; Liu, H.; Ma, Y.; Ding, Y. The Impacts of Climate Change on Water Resources and Agriculture in China. Nature 2010, 467, 43–51. [Google Scholar] [CrossRef] [PubMed]
  5. Li, H.; Yan, L.; Wen, T.; Feng, X. Characteristics of Climate Change and Its Impact Assessment in the Three-River Regions. Plateau Meteorol. 2022, 41, 306–316. (In Chinese) [Google Scholar]
  6. Luo, S.; Fang, X.; Lyu, S.; Ma, D.; Chang, Y.; Song, M.; Chen, H. Frozen Ground Temperature Trends Associated with Climate Change in the Tibetan Plateau Three River Source Region from 1980 to 2014. Clim. Res. 2016, 67, 241–255. [Google Scholar] [CrossRef] [Green Version]
  7. Cheng, G.; Zhao, L.; Li, R.; Wu, X.; Sheng, Y.; Hu, G.; Zou, D.; Jin, H.; Li, X.; Wu, Q. Characteristic, changes and impacts of permafrost on Qinghai-Tibet Plateau. Chin. Sci. Bull. 2019, 64, 2783–2795. (In Chinese) [Google Scholar]
  8. Zou, D.; Zhao, L.; Sheng, Y.; Chen, J.; Hu, G.; Wu, T.; Wu, J.; Xie, C.; Wu, X.; Pang, Q.; et al. A New Map of Permafrost Distribution on the Tibetan Plateau. Cryosphere 2017, 11, 2527–2542. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, T.; Yang, D.; Fang, B.; Yang, W.; Qin, Y.; Wang, Y. Data-Driven Mapping of the Spatial Distribution and Potential Changes of Frozen Ground over the Tibetan Plateau. Sci. Total Environ. 2019, 649, 515–525. [Google Scholar] [CrossRef]
  10. Luo, S.; Wang, J.; Pomeroy, J.W.; Lyu, S. Freeze–Thaw Changes of Seasonally Frozen Ground on the Tibetan Plateau from 1960 to 2014. J. Clim. 2020, 33, 9427–9446. [Google Scholar] [CrossRef]
  11. Deluigi, N.; Lambiel, C.; Kanevski, M. Data-Driven Mapping of the Potential Mountain Permafrost Distribution. Sci. Total Environ. 2017, 590–591, 370–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gao, S.; Wu, Q.; Zhang, Z.; Xu, X. Impact of Climatic Factors on Permafrost of the Qinghai–Xizang Plateau in the Time-Frequency. Domain Quat. Int. 2015, 374, 110–117. [Google Scholar] [CrossRef]
  13. Xing, L.; Li, L.; Gong, J.; Ren, C.; Liu, J.; Chen, H. Daily Soil Temperatures Predictions for Various Climates in United States Using Data-Driven Model. Energy 2018, 160, 430–440. [Google Scholar] [CrossRef]
  14. Ran, Y.; Li, X.; Cheng, G.; Nan, Z.; Che, J.; Sheng, Y.; Wu, Q.; Jin, H.; Luo, D.; Tang, Z.; et al. Mapping the Permafrost Stability on the Tibetan Plateau for 2005–2015. Sci. China Earth Sci. 2020, 64, 62–79. [Google Scholar] [CrossRef]
  15. Jin, H.; Sun, K.; Wang, S.; He, R.; Lv, L.; Yu, S. Dual Influences of Local Environmental Variables on Ground Temperatures on the Interior-Eastern Qingha-i Tibet Plateau (Ⅰ): Vegetation and Snow Cover. J. Glaciol. Geocryol. 2008, 30, 535–545. (In Chinese) [Google Scholar]
  16. Smith, M.W. Microclimatic Influences on Ground Temperatures and Permafrost Distribution, Mackenzie Delta, Northwest Territories. Can. J. Earth Sci. 1975, 12, 1421–1438. [Google Scholar] [CrossRef]
  17. Ni, J.; Cheng, Y.; Wang, Q.; Ng, C.W.W.; Garg, A. Effects of Vegetation on Soil Temperature and Water Content: Field Monitoring and Numerical Modelling. J. Hydrol. 2019, 571, 494–502. [Google Scholar] [CrossRef]
  18. Flanner, M.G.; Shell, K.M.; Barlage, M.; Perovich, D.K.; Tschudi, M. Radiative Forcing and Albedo Feedback from the Northern Hemisphere Cryosphere between 1979 and 2008. Nat. Geosci. 2011, 4, 151–155. [Google Scholar] [CrossRef]
  19. Li, W.; Luo, S.; Hao, X.; Wang, J.; Wang, Y. Observations of East Qinghai-Xizang Plateau Snow Cover Effects on Surface Energy and Water Exchange in Different Seasons. Plateau Meteorol. 2021, 40, 455–471. (In Chinese) [Google Scholar]
  20. Fu, Q.; Hou, R.; Li, T.; Wang, M.; Yan, J. The Functions of Soil Water and Heat Transfer to the Environment and Associated Response Mechanisms under Different Snow Cover Conditions. Geoderma 2018, 325, 9–17. [Google Scholar] [CrossRef]
  21. Zhang, T. Influence of the Seasonal Snow Cover on the Ground Thermal Regime: An Overview. Rev. Geophys. 2005, 43, RG4002. [Google Scholar] [CrossRef]
  22. Qian, B.; Gregorich, E.G.; Gameda, S.; Hopkins, D.W.; Wang, X.L. Observed Soil Temperature Trends Associated with Climate Change in Canada. J. Geophys. Res. Atmos. 2011, 116, D02106. [Google Scholar] [CrossRef]
  23. Aalto, J.; le Roux, P.C.; Luoto, M. Vegetation Mediates Soil Temperature and Moisture in Arctic-Alpine Environments. Arct. Antarct. Alp. Res. 2018, 45, 429–439. [Google Scholar] [CrossRef] [Green Version]
  24. Feng, H.; Ye, S.; Zou, B. Contribution of Vegetation Change to the Surface Radiation Budget: A Satellite Perspective. Glob. Planet. Change 2020, 192, 103225. [Google Scholar] [CrossRef]
  25. Lu, Y.; Yu, W.; Guo, M.; Liu, W. Spatiotemporal variation characteristics of land cover and land surface temperature in Mohe County, Helongjiang Province. J. Glaciol. Geocryol. 2017, 39, 1137–1149. (In Chinese) [Google Scholar]
  26. Way, R.G.; Lapalme, C.M. Does Tall Vegetation Warm or Cool the Ground Surface? Constraining the Ground Thermal Impacts of Upright Vegetation in Northern Environments. Environ. Res. Lett. 2021, 16, 054077. [Google Scholar] [CrossRef]
  27. Chang, X.; Jin, H.; Wang, Y.; Zhang, Y.; Zhou, G.; Che, F.; Zhao, Y. Influences of vegetation on permafrost: A review. Acta Ecol. Sin. 2012, 32, 7981–7990. (In Chinese) [Google Scholar] [CrossRef] [Green Version]
  28. Luo, D.; Jin, H.; He, R.; Yang, S. Responses of Surface Vegetation on Soil Temperature and Moisture of the Active Layer in the Source Area of the Yellow River. J. Earth Sci. J. China Univ. Geosci. 2014, 39, 421–430. (In Chinese) [Google Scholar]
  29. Song, Y.; Zhou, D.; Zhang, H.; Li, G.; Jin, Y.; Li, Q. Effects of Vegetation Height and Density on Soil Temperature Variations. Chin. Sci. Bull. 2013, 58, 907–912. [Google Scholar] [CrossRef] [Green Version]
  30. Kun, Y.; Jie, H.; China Meteorological Forcing Dataset (1979–2018). A Big Earth Data Platf. Three Poles 2019. Available online: http://data.tpdc.ac.cn/en/data/8028b944-daaa-4511-8769-965612652c49/ (accessed on 19 July 2022).
  31. Yang, K.; He, J.; Tang, W.; Qin, J.; Cheng, C.C. On Downward Shortwave and Longwave Radiations over High Altitude Regions: Observation and Modeling in the Tibetan Plateau. Agric. For. Meteorol. 2010, 150, 38–46. [Google Scholar] [CrossRef]
  32. He, J.; Yang, K.; Tang, W.; Lu, H.; Qin, J.; Chen, Y.; Li, X. The First High-Resolution Meteorological Forcing Dataset for Land. Process Stud. Over China Sci. Data 2020, 7, 1–11. [Google Scholar]
  33. Hao, X.; Huang, G.; Che, T.; Ji, W.; Sun, X.; Zhao, Q.; Zhao, H.; Wang, J.; Li, H.; Yang, Q. The Nieer Avhrr Snow Cover Extent Product over China–a Long-Term Daily Snow Record for Regional Climate Research. Earth Syst. Sci. Data 2021, 13, 4711–4726. [Google Scholar] [CrossRef]
  34. Wu, D.; Wu, H.; Zhao, X.; Zhou, T.; Tang, B.; Zhao, W.; Jia, K. Evaluation of Spatiotemporal Variations of Global Fractional Vegetation Cover Based on Gimms Ndvi Data from 1982 to 2011. Remote Sens. 2014, 6, 4217–4239. [Google Scholar] [CrossRef] [Green Version]
  35. Shangguan, W.; Dai, Y.; Liu, B.; Zhu, A.; Duan, Q.; Wu, L.; Ji, D.; Ye, A.; Yuan, H.; Zhang, Q. A China Data Set of Soil Properties for Land Surface Modeling Journal of Advances in Modeling. Earth Syst. 2013, 5, 212–224. [Google Scholar]
  36. Hamed, K.H.; Rao, A.R. A Modified Mann-Kendall Trend Test for Autocorrelated Data. J. Hydrol. 1998, 204, 182–196. [Google Scholar] [CrossRef]
  37. Deluigi, N.; Lambiel, C. Permal: A Machine Learning Approach for Alpine Permafrost Distribution Modeling; Eidg. Forschungsanstalt WSL: Zürcherstrasse, Switzerland, 2013. [Google Scholar]
  38. Gunn, S.R. Support Vector Machines for Classification and Regression. ISIS Tech. Rep. 1998, 14, 5–16. [Google Scholar]
  39. Vapnik, V.N. An Overview of Statistical Learning Theory. IEEE Trans. Neural Netw. 1999, 10, 988–999. [Google Scholar] [CrossRef] [Green Version]
  40. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  41. Hutengs, C.; Vohland, M. Downscaling Land Surface Temperatures at Regional Scales with Random Forest Regression. Remote Sens. Environ. 2016, 178, 127–141. [Google Scholar] [CrossRef]
  42. Wang, J.; Luo, S.; Li, Z.; Wang, S.; Li, Z. The Freeze/Thaw Process and the Surface Energy Budget of the Seasonally Frozen Ground in the Source Region of the Yellow River. Theor. Appl. Climatol. 2019, 138, 1631–1646. [Google Scholar] [CrossRef] [Green Version]
  43. Hussain, M.M.; Mahmud, I. Pymannkendall: A Python Package for Non Parametric Mann Kendall Family of Trend. Tests J. Open Source Softw. 2019, 4, 1556. [Google Scholar] [CrossRef]
  44. Luo, S.; Li, H.; Ma, D.; Li, W.; Wang, J.; Tan, X.; Dong, Q. Review and Prospects of Frozen Soil-Vegetation Interaction and Climate Effects in the Three Rivers Source Region. Plateau Meteorol. 2022, 41, 255–267. (In Chinese) [Google Scholar]
Figure 1. Three River Source Region and observation station.
Figure 1. Three River Source Region and observation station.
Remotesensing 14 04114 g001
Figure 2. Observation-simulation comparison of MLR, SVM, and RF models.
Figure 2. Observation-simulation comparison of MLR, SVM, and RF models.
Remotesensing 14 04114 g002aRemotesensing 14 04114 g002b
Figure 3. Comparison of simulated values and observed values.
Figure 3. Comparison of simulated values and observed values.
Remotesensing 14 04114 g003aRemotesensing 14 04114 g003b
Figure 4. The series of SCD anomaly (a,c), NDVI anomaly (b,d), soil temperature anomaly and the temperature difference between the atmosphere and the soil.
Figure 4. The series of SCD anomaly (a,c), NDVI anomaly (b,d), soil temperature anomaly and the temperature difference between the atmosphere and the soil.
Remotesensing 14 04114 g004aRemotesensing 14 04114 g004b
Figure 5. Multi-year average monthly SCD (a), NDVI (b) and temperature difference between atmosphere and soil from 0 cm to 40 cm.
Figure 5. Multi-year average monthly SCD (a), NDVI (b) and temperature difference between atmosphere and soil from 0 cm to 40 cm.
Remotesensing 14 04114 g005
Figure 6. SCD (a) and its trends (b) (Pass significance test at 90% confidence level).
Figure 6. SCD (a) and its trends (b) (Pass significance test at 90% confidence level).
Remotesensing 14 04114 g006
Figure 7. NDVI (a) and its trends (b) (Pass significance test at 90% confidence level).
Figure 7. NDVI (a) and its trends (b) (Pass significance test at 90% confidence level).
Remotesensing 14 04114 g007
Figure 8. The correlation coefficient of SCD and soil temperature (slashed areas indicate passing the 95% significance test).
Figure 8. The correlation coefficient of SCD and soil temperature (slashed areas indicate passing the 95% significance test).
Remotesensing 14 04114 g008
Figure 9. The correlation coefficient of NDVI and soil temperature (slashed areas indicate passing the 95% significance test).
Figure 9. The correlation coefficient of NDVI and soil temperature (slashed areas indicate passing the 95% significance test).
Remotesensing 14 04114 g009
Figure 10. The difference between soil temperature and air temperature.
Figure 10. The difference between soil temperature and air temperature.
Remotesensing 14 04114 g010
Table 1. Soil temperature simulation testing results for MLR, SVM, and RF models.
Table 1. Soil temperature simulation testing results for MLR, SVM, and RF models.
R2RMSE (°C)Bias (°C)
MLRSVMRFMLRSVMRFMLRSVMRF
0 cm0.9868 0.9578 0.9913 1.17092.07410.95810.0250−0.17260.0168
5 cm0.9841 0.9720 0.9935 1.11711.48710.7172−0.0676−0.1495−0.0524
10 cm0.9821 0.9727 0.9933 1.14511.41710.7032−0.0550−0.1138−0.0369
15 cm0.9811 0.9727 0.9931 1.15571.40050.6974−0.0621−0.1150−0.0320
20 cm0.9779 0.9720 0.9926 1.20551.35840.7063−0.0630−0.1140−0.0181
40 cm0.9495 0.9588 0.9862 1.65901.48040.8570−0.0520−0.12500.0102
80 cm0.9289 0.9437 0.9811 1.60371.43260.8267−0.1547−0.1561−0.0433
160 cm0.8786 0.9517 0.9769 1.63671.03030.7104−0.0813−0.0401−0.0446
320 cm0.8542 0.9400 0.9735 1.34200.85270.5704−0.0472−0.0283−0.0275
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, X.; Luo, S.; Li, H.; Hao, X.; Wang, J.; Dong, Q.; Chen, Z. Investigating the Effects of Snow Cover and Vegetation on Soil Temperature Using Remote Sensing Indicators in the Three River Source Region, China. Remote Sens. 2022, 14, 4114. https://doi.org/10.3390/rs14164114

AMA Style

Tan X, Luo S, Li H, Hao X, Wang J, Dong Q, Chen Z. Investigating the Effects of Snow Cover and Vegetation on Soil Temperature Using Remote Sensing Indicators in the Three River Source Region, China. Remote Sensing. 2022; 14(16):4114. https://doi.org/10.3390/rs14164114

Chicago/Turabian Style

Tan, Xiaoqing, Siqiong Luo, Hongmei Li, Xiaohua Hao, Jingyuan Wang, Qingxue Dong, and Zihang Chen. 2022. "Investigating the Effects of Snow Cover and Vegetation on Soil Temperature Using Remote Sensing Indicators in the Three River Source Region, China" Remote Sensing 14, no. 16: 4114. https://doi.org/10.3390/rs14164114

APA Style

Tan, X., Luo, S., Li, H., Hao, X., Wang, J., Dong, Q., & Chen, Z. (2022). Investigating the Effects of Snow Cover and Vegetation on Soil Temperature Using Remote Sensing Indicators in the Three River Source Region, China. Remote Sensing, 14(16), 4114. https://doi.org/10.3390/rs14164114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop