Abstract:The ground surface temperature (GST) serves as a crucial indicator for understanding land-atmosphere mass and energy exchange. The shift from manual measurement to automated station for GST in China after 2002 introduced inconsistencies at certain stations, potentially distorting research findings. Here, daily automatedly observed GST from 2003 to 2017 at 615 selected meteorological stations were updated by constructing linear regression model based on manually observed air temperature (AT) and GST from 1960 to 2002. Then, the spatiotemporal variations of GST from 1960 to 2017 and its driving factors were investigated. Results indicated that: (1) The AT-GST linear regression model could effectively mitigate the inconsistency caused by the change of GST observation methods, enhancing data reliability. (2) GST in China showed little change from 1960-1980, but increased significantly across all regions from 1980 to 2000, with the increase rate slowed down except in the Qinghai-Tibet Plateau (QTP) and southwest China after 2000. Notable GST increase is concentrated in colder regions, including the QTP, northeast (NEC), and northwest China (NWC). (3) Evapotranspiration (ET) and vapor pressure deficit (VPD) were the primary drivers of annual GST variations at the regional scale, while their contributions to GST variations exhibited notable seasonal variability. Our findings could offer valuable scientific insights for addressing climate change, enhancing surface environmental models, and safeguarding ecological environments.
What problem does this paper attempt to address?
This paper attempts to address the issues of the spatio - temporal variation of the ground surface temperature (GST) in China and its driving factors. Specifically, the paper focuses on the following two main issues:
1. **How to construct a simple yet reliable GST correction model to eliminate the uncertainties caused by changes in GST observation methods?**
- The paper corrects the inconsistencies in the automatic observation data after 2003 by establishing a linear regression model based on the manual observation data from 1960 to 2002. Specifically, the correction is carried out by using the linear relationship between air temperature (AT) and GST to obtain continuous and reliable historical GST data.
2. **What are the spatio - temporal variation patterns of GST in China and its main driving factors?**
- The paper analyzes the spatio - temporal variation patterns of GST in China from 1960 to 2017 and explores the main driving factors affecting these changes. The study finds that there are significant differences in the spatio - temporal variation of GST in different regions of China. Especially after 1980, GST has increased significantly in most regions, but the change trends in the Qinghai - Tibet Plateau (QTP) and the southwestern region (SWC) are different. In addition, evapotranspiration (ET) and vapor pressure deficit (VPD) are the main driving factors for the annual GST change, but their contributions show significant seasonal differences in different seasons.
### Main research contents and methods
1. **Data sources**:
- The study uses the daily GST and AT data of 615 meteorological stations from 1960 to 2017 provided by the China Meteorological Administration.
- Other auxiliary data include precipitation (Prep), soil moisture (SM), evapotranspiration (ET), surface solar radiation (Rs), snow water equivalent (SWE), vapor pressure deficit (VPD), leaf area index (LAI), as well as irrigated area (Irri) and land cover types (VA, UA, WA), etc.
2. **GST correction model**:
- Using the manual observation data from 1960 to 2002, a monthly linear regression model is established to correct the automatic observation data from 2003 to 2017.
- The linear relationship of the model is expressed as:
\[
GST_{m,n}=K_{m,n}\times AT_{m,n}+B_{m,n}
\]
where \( GST_{m,n} \) is the corrected GST of the \( n \) - th station in the \( m \) - th month, \( AT_{m,n} \) is the air temperature of the \( n \) - th station in the \( m \) - th month, and \( K_{m,n} \) and \( B_{m,n} \) are the slope and intercept between GST and AT, respectively.
3. **Spatio - temporal variation analysis**:
- The Kriging method is used to perform spatial interpolation on the updated monthly GST data to generate a GST distribution map with a resolution of 0.1°×0.1°.
- The Theil - Sen median method and the Mann - Kendall trend test method are used to analyze the change trends of GST at each grid point from 1960 to 2017.
- The period from 1960 to 2017 is divided into three stages (1960 - 1979, 1980 - 1999, 2000 - 2017), and the change trends of GST in each region are analyzed respectively.
4. **Driving factor analysis**:
- Canonical analysis is used to evaluate the main driving factors of GST change, including precipitation, evapotranspiration, surface solar radiation, soil moisture, snow water equivalent, vapor pressure deficit, leaf area index, irrigated area and land cover types, etc.
### Main conclusions
- **Effectiveness of the GST correction model**:
- The correction model can effectively reduce the inconsistencies caused by changes in GST observation methods and improve the reliability of the data.
- The evaluation indicators (MAE, RMSE, NSE, R²) of the model all perform well, which verifies the accuracy of the model.
- **Spatio - temporal variation of GST**:
- From 1960 to 1980, the GST in China did not change much; from 1980 to 2000, GST increased significantly in all regions; after 2000, the increase rate of GST in other regions except QTP and SWC has slowed down.
- Spring and winter are...