The data set is the global vegetation productivity data, including Gross Primary Productivity(GPP) and Net Primary Productivity (NPP). It was obtained by the CNRM-CM6-1 mode simulation of CMIP6 under the Historical scenario. The time range of the data covers from 1850 to 2014, the time resolution is a month, and the spatial resolution is about 1.406°×1.389°. For the simulated data details, please go to the following link: http://www.umr-cnrm.fr/cmip6/spip.php?article11.

Taking 2005 as the base year, the future population scenario was predicted by adopting the Logistic model of population. It not only can better describe the change pattern of population and biomass but is also widely applied in the economic field. The urbanization rate was predicted by using the urbanization Logistic model. Based on the existing urbanization horizontal sequence value, the prediction model was established by acquiring the parameters in the parametric equation by nonlinear regression. The urban population was calculated by multiplying the predicted population by the urbanization rate. The data adopted the non-agricultural population. The Logistic model was used to predict the future gross national product of each county (or city), and then, according to the economic development level of each county (or city) in each period (in terms of GDP per capita),the corresponding industrial structure scenarios in each period were set, and the output value of each industry was predicted. The trend of changes in industrial structure in China and the research area lagged behind the growth of GDP and was therefore adjusted according to the need of the future industrial structure scenarios of the research area.

This dataset includes data recorded by the Cold and Arid Research Network of Lanzhou university obtained from an observation system of Meteorological elements gradient of Xiyinghe Station from January 1 to December 31, 2018. The site (101.853E, 37.561N) was located on a alpine meadow in the Menyuan,Qinghai Province. The elevation is 3639 m. The installation heights and orientations of different sensors and measured quantities were as follows: air temperature and humidity profile (2, 4, and 8 m, towards north), wind speed and direction profile (windsonic; 2, 4, and 8 m, towards north), air pressure (1.5 m), rain gauge (4 m), four-component radiometer (4 m, towards south), infrared temperature sensors (4 m, towards south, vertically downward), photosynthetically active radiation (4 m, towards south), soil heat flux (-0.05 m and -0.1m in south of tower), soil soil temperature/ moisture/ electrical conductivity profile (-0.2 and -0.4 m in south of tower), sunshine duration sensor (4 m, towards south). The observations included the following: air temperature and humidity (Ta_2 m, Ta_4 m, and Ta_8 m; RH_2 m, RH_4 m, and RH_8 m) (℃ and %, respectively), wind speed (Ws_2 m, Ws_4 m, and Ws_8 m) (m/s), wind direction (WD_2 m, WD_4 m, and WD_8 m) (°), air pressure (press) (hpa), precipitation (rain) (mm), four-component radiation (DR, incoming shortwave radiation; UR, outgoing shortwave radiation; DLR_Cor, incoming longwave radiation; ULR_Cor, outgoing longwave radiation; Rn, net radiation) (W/m^2), infrared temperature (IRT) (℃), photosynthetically active radiation (PAR) (μmol/ (s/m^2)), soil heat flux (Gs_5 cm, Gs_10cm) (W/m^2), soil temperature (Ts_20 cm, Ts_40 cm) (℃), soil moisture (Ms_20 cm, Ms_40 cm) (%, volumetric water content), soil water potential (SWP_20cm , SWP_40cm)(kpa) , soil conductivity (Ec_20cm, Ec_40cm)(μs/cm), sun time (h). The data processing and quality control steps were as follows: (1) The AWS data were averaged over intervals of 10 min for a total of 144 records per day. The meteorological data were missing during Aug. 29 to Oct.18 because of unstable power supply due to battery box flooding; The wind speed and direction profile data were rejected because of sensor failure; The precipitation data were rejected because of program error; The air humidity data before Mar. 2 were rejected due to program error; (2) Data in duplicate records were rejected. (3) Unphysical data were rejected. (4) The data marked in red are problematic data. (5) The format of the date and time was unified, and the date and time were collected in the same column, for example, date and time: 2018-6-10 10:30.

This dataset contains the automatic weather station (AWS) measurements from Zhangye wetland station in the flux observation matrix from 25 June to 21 September, 2012. The site (100.44640° E, 38.97514° N) was located in a wetland surface, which is near Zhangye city, Gansu Province. The elevation is 1460 m. The installation heights and orientations of different sensors and measured quantities were as follows: air temperature and humidity (HMP45AC; 5 m and 10 m, towards north), air pressure (CS100; 2 m), rain gauge (TE525M; 10 m), wind speed (03002; 5 m and 10 m, towards north), wind direction (03002; 10 m, towards north), a four-component radiometer (NR01; 6 m, towards south), two infrared temperature sensors (SI-111; 6 m, vertically downward), soil temperature profile (109ss-L; 0, -0.02, -0.04, -0.1, -0.2, and -0.4 m), and soil heat flux (HFP01; 3 duplicates, 0.06 m). The observations included the following: air temperature and humidity (Ta_5 m and Ta_10 m, RH_5 m and RH_10 m) (℃ and %, respectively), air pressure (press, hpa), precipitation (rain, mm), wind speed (Ws_5 m and Ws_10 m, m/s), wind direction (WD_10 m, °), four-component radiation (DR, incoming shortwave radiation; UR, outgoing shortwave radiation; DLR_Cor, incoming longwave radiation; ULR_Cor, outgoing longwave radiation; Rn, net radiation; W/m^2), infrared temperature (IRT_1 and IR_2, ℃), soil heat flux (Gs_1, Gs_2 and Gs_3, W/m^2), and soil temperature profile (Ts_0 cm, Ts_2 cm, Ts_4 cm, Ts_10 cm, Ts_20 cm, Ts_40 cm, ℃). The data processing and quality control steps were as follows. (1) The AWS data were averaged over intervals of 10 min; therefore, there were 144 records per day. The missing data were filled with -6999. (2) Data in duplicate records were rejected. (3) Unphysical data were rejected. (4) In this dataset, the time of 0:10 corresponds to the average data for the period between 0:00 and 0:10; the data were stored in *.xlsx format. (5) Finally, the naming convention was AWS+ site no. Moreover, suspicious data were marked in red. For more information, please refer to Liu et al. (2016) (for multi-scale observation experiment or sites information), Xu et al. (2013) (for data processing) in the Citation section.

This dataset contains the flux measurements from the Daman superstation eddy covariance system (EC) at the lowest layer in the flux observation matrix from 25 May to 15 September, 2012. The site (100.37223° E, 38.85551° N) was located in a cropland (maize surface) in the Daman irrigation district, which is near Zhangye, Gansu Province. The elevation is 1556.06 m. The EC was installed at a height of 4.5 m; the sampling rate was 10 Hz. The sonic anemometer faced north, and the separation distance between the sonic anemometer and the CO2/H2O gas analyzer (CSAT3&Li7500A) was 0.17 m. Raw data acquired at 10 Hz were processed using the Edire post-processing software (University of Edinburgh, http://www.geos.ed.ac.uk/abs/research/micromet/EdiRe/), including spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, coordinate rotation (2-D rotation), corrections for density fluctuation (Webb-Pearman-Leuning correction), and frequency response correction. The EC data were subsequently averaged over 30 min periods. Moreover, the observation data quality was divided into three classes according to the quality assessment method of stationarity (Δst) and the integral turbulent characteristics test (ITC), which was proposed by Foken and Wichura [1996]: class 1 (level 0: Δst<30 and ITC<30), class 2 (level 1: Δst<100 and ITC<100), and class 3 (level 2: Δst>100 and ITC>100), representing high-, medium-, and low-quality data, respectively. In addition to the above processing steps, the half-hourly flux data were screened in a four-step procedure: (1) data from periods of sensor malfunction were rejected; (2) data before or after 1 h of precipitation were rejected; (3) incomplete 30 min data were rejected when the missing data constituted more than 3% of the 30 min raw record; and (4) data were rejected at night when the friction velocity (u*) was less than 0.1 m/s. There were 48 records per day; the missing data were replaced with -6999. Moreover, suspicious data were marked in red. The released data contained the following variables: data/time, wind direction (Wdir, °), wind speed (Wnd, m/s), the standard deviation of the lateral wind (Std_Uy, m/s), virtual temperature (Tv, ℃), H2O mass density (H2O, g/m^3), CO2 mass density (CO2, mg/m^3), friction velocity (ustar, m/s), stability (z/L), sensible heat flux (Hs, W/m2), latent heat flux (LE, W/m^2), carbon dioxide flux (Fc, mg/ (m^2s)), quality assessment of the sensible heat flux (QA_Hs), quality assessment of the latent heat flux (QA_LE), and quality assessment of the carbon flux (QA_Fc). In this dataset, the time of 0:30 corresponds to the average data for the period between 0:00 and 0:30; the data were stored in *.xlsx format. For more information, please refer to Liu et al. (2016) (for multi-scale observation experiment or sites information), Xu et al. (2013) (for data processing) in the Citation section.

The High Asia region is an area sensitive to global changes in mid-latitude regions and is a hotspot for research. The lakes in the territory are scattered, and the lake freeze-thaw process is one of the key factors sensitive to global change. Due to the large difference in the dielectric constant between ice and water, satellite-borne passive microwave remote sensing is weather insensitive and has a high revisiting rate; thus, it can achieve rapid monitoring of the freeze-thaw state of lakes. According to the area ratio of the lake and the land surface in the sub-pixels of passive microwave radiometer data, this data set represents the lake brightness temperature information of the pixel (sub-pixel level) by applying the hybrid pixel decomposition method in order to monitor the lake freeze-thaw process in the High Asia region. Thus, by adopting a variety of passive microwave data, time series of lake brightness temperature and freeze-thaw status were obtained for a total of 51 medium to large lakes from 2002 to 2016 in the High Asia region. Using cloudless MODIS optical products as validation data, three lakes of different sizes in different regions of High Asia, i.e., Hoh Xil Lake, Dagze Co Lake, and Kusai Lake, were selected for freeze-thaw detection validation. The results indicated that the lake freeze-thaw parameters obtained by microwave and optical remote sensing were highly consistent, and the correlation coefficients reached 0.968 and 0.987. This data set contained the time series brightness temperature of lakes and the freeze-thaw parameters of lake ice, which could be used to further invert the characteristic parameters of lakes and enhance the understanding of lake ice freezing and thawing in the High Asia region. This database will be useful in the assessment of climatic and environmental changes in the High Asia region and in global climatic change response models. The data set consists of two parts: the passive microwave remote sensing brightness temperature data set of 51 lakes in the High Asia region from 2002 to 2016, with an observation interval of 1 to 2 days, and the lake ice freeze-thaw data set obtained by estimation of the lake brightness temperature. The files are the lake brightness temperature data via the nearest neighbour method and pixel decomposition in the form of a .zip file (12 MB) and the lake freeze-thaw data set for 51 lakes in the High Asia region from 2002 to 2016 in the form of an .xls file (0.1 MB).

This dataset contains the flux measurements from the large aperture scintillometer (LAS) at Daman Superstation in the hydrometeorological observation network of Heihe River Basin between 15 September, 2012, and 31 December, 2013. There were two types of LASs at Daman Superstation: German BLS450 (labeled as NQ and AR) and Netherlands Kipp&zonen. The north tower was set up with the Kipp&zonen/BLS450_AR receiver and the BLS450_NQ transmitter, and the south tower was equipped with the Kipp&zonen/BLS450_AR transmitter and the BLS450_NQ receiver. BLS450_NQ has been in use since 26 September, 2012, Kipp&zonen has been in use since 23 September, 2013, and the observation period of BLS450_AR was from 15 September, 2012, to 25 July, 2013. The site (north: 100.379° E, 38.861° N; south: 100.369° E, 38.847° N) was located in Daman irrigation district, which is near Zhangye, Gansu Province. The underlying surfaces between the two towers were corn, orchard, and greenhouse. The elevation is 1556 m. The effective height of the LASs was 22.45 m, and the path length was 1854 m. The data were sampled at 5 Hz and 1 Hz intervals for BLS450 and Kipp&zonen, respectively, then averaged to 1 minute. The raw data acquired at 1 min intervals were processed and quality controlled. The data were subsequently averaged over 30 min periods, in which sensible heat flux was iteratively calculated by combining Cn2 with meteorological data according to the Monin-Obukhov similarity theory. The main quality control steps were as follows: (1) The data were rejected when Cn2 exceeded the saturated criterion (BLS450: Cn2>1.43E-13, Kipp&zonen: Cn2>1.54E-13). (2) The data were rejected when the demodulation signal was small (BLS450: Average X Intensity<1000; Kipp&zonen: Demod>-20mv). (3) The data were rejected when collected during precipitation. (4) The data were rejected if collected at night when weak turbulence occurred (u* was less than 0.1 m/s). In the iteration process, the universal functions of Thiermann and Grassl, 1992 and Andreas, 1988 were selected for BLS450 and Kipp&zonen, respectively. Several instructions were included with the released data. (1) The data were primarily obtained from BLS450_NQ measurements, and missing flux measurements from the BLS450_NQ instrument were substituted with measurements from the BLS450_AR and Kipp&zonen instrument. The missing data were denoted by -6999. Due to the problems of BLS450_NQ SPU storing and wireless transmission, large amount of data from 11 August to 17 August, 18 August to 20 August, 22 August to 24 August, 27 August to 30 August, 2013, and 1 September to 3 September, 19 September to 23 September, 2013, were not collected. (2) The dataset contained the following variables: data/time (yyyy-m-d h:mm), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). In this dataset, a time of 0:30 corresponds to the average data for the period between 0:00 and 0:30, and the data were stored in *.xls format. Moreover, suspicious data were marked in red. For more information, please refer to Li et al. (2013) (for hydrometeorological observation network or sites information), Liu et al. (2011) (for data processing) in the Citation section.

On 29 June 2012, CASI sensor carried by the Harbin Y-12 aircraft was used in a visible near Infrared hyperspectral airborne remote sensing experiment, which is located in the observation experimental area (30×30 km). The land cover pattern product in the middle reaches of the Heihe River Basin were obtained at a spatial resolution of 1 m, using CASI aerial data with high spatial and spectral resolution．A hierarchical classification structure integrated by pixel-based classification and object-based classification is used to obtain production．According to surveyed reference data about land cover and visual interpretation from high resolution imagery，the accuracy of the classification result of land cover was evaluated，and the result showed that overall accuracy was 84.61 %，Kappa coefficient was 0.8262．

The dataset includes the fractional vegetation cover data generated from the stations of crop land, wetland, Gebi desert and desert steppe in Yingke Oasis and biomass data generated from the stations of crop land (corn) and wetland. The observations lasted for a vegetation growth cycle from 19 May, 2012 to 15 September, 2012. 1. Fractional vegetation cover observation 1.1 Observation time 1.1.1 Station of the crop land: The observations lasted from 20 May, 2012 to 15 September, 2012, and in five-day periods for each observation before 31 July and in ten-day periods for each observation after 31 July. The observation time for the station of crop land (corn) are 2013-5-20, 2013-5-25, 2013-5-30, 2013-6-5, 2013-6-10, 2013-6-16, 2013-6-22, 2013-6-27, 2013-7-2, 2013-7-7, 2013-7-12, 2013-7-17, 2013-7-27, 2013-8-3, 2013-8-13, 2013-8-25, 2013-9-5 and 2013-9-15. 1.1.2 The other four stations: The observations lasted from 20 May, 2012 to 15 September, 2012 and in ten-day periods for each observation. The observation time for the crop land are 2013-5-20, 2013-6-5, 2013-6-16, 2013-6-27, 2013-7-7, 2013-7-17, 2013-7-27, 2013-8-3, 2013-8-13, 2013-8-25, 2013-9-5 and 2013-9-15. 1.2 method 1.2.1 Instruments and measurement method Digital photography measurement is implemented to measure the FVC. Plot positions, photographic method and data processing method are dedicatedly designed. In field measurements, a long stick with the camera mounted on one end is beneficial to conveniently measure various species of vegetation, enabling a larger area to be photographed with a smaller field of view. The stick can be used to change the camera height; a fixed-focus camera can be placed at the end of the instrument platform at the front end of the support bar, and the camera can be operated by remote control. 1.2.2 Design of the samples Three and two plots with the area of 10×10 m^2 were measured for the station of the crop land and wetland, respectively. One plot with the area of 10×10 m^2 was measured for the other three stations. Shoot 9 times along two perpendicularly crossed rectangular-belt transects. The picture generated of each time is used to calculate a FVC value. “True FVC” of the plot is then acquired as the average of these 9 FVC values. 1.2.3 Photographic method The photographic method used depends on the species of vegetation and planting pattern. A long stick with the camera mounted on one end is used for the stations of crop land and wetland. For the station of the crop land, rows of more than two cycles should be included in the field of view (<30), and the side length of the image should be parallel to the row. If there are no more than two complete cycles, then information regarding row spacing and plant spacing are required. The FVC of the entire cycle, that is, the FVC of the quadrat, can be obtained from the number of rows included in the field of view. For other three stations, the photos of FVC were obtained by directly photographing for the lower heights of the vegetation. 1.2.4 Method for calculating the FVC The FVC calculation was implemented by the Beijing Normal University. The detail method can be found in the reference below. Many methods are available to extract the FVC from digital images, and the degree of automation and the precision of identification are important factors that affect the efficiency of field measurements. This method, which is proposed by the authors, has the advantages of a simple algorithm, a high degree of automation and high precision, as well as ease of operation (see the reference). 2. Biomass observation 2.1. Observation time 2.1.1 Station of the crop land: The observations lasted from 20 May 2012 to 15 September 2012, and in five-day periods for each observation before 31 July and in ten-day periods for each observation after 31 July. The observation time for the crop land are 2013-5-25, 2013-5-30, 2013-6-5, 2013-6-10, 2013-6-16, 2013-6-22, 2013-6-27, 2013-7-2, 2013-7-7, 2013-7-12, 2013-7-17, 2013-7-27, 2013-8-3, 2013-8-13, 2013-8-25, 2013-9-5 and 2013-9-15. 2.1.2 The station of wetland: The observations lasted from 20 May 2012 to 15 September 2012, and in ten-day periods for each observation. The observation time for the crop land are 2013-6-5, 2013-6-16, 2013-6-27, 2013-7-7, 2013-7-17, 2013-7-27, 2013-8-3, 2013-8-13, 2013-8-25, 2013-9-5 and 2013-9-15. 2.2. Method Station of the crop land: Three plots were selected and three strains of corn for each observation were random selected for each plot to measure the fresh weight (the aboveground biomass and underground biomass) and dry weight. Per unit biomass can be obtained according to the planting structure. Station of the wetland: Two plots of reed with the area of 0.5 m × 0.5 m were random selected for each observation. The reed of the two plots was cut to measure the fresh weight (the aboveground biomass) and dry weight. 2.3. Instruments Balance (accuracy 0.01 g); drying oven 3. Data storage All observation data were stored in excel. Other data including plant spacing, row spacing, seeding time, irrigation time, the time of cutting male parent and the harvest time of the corn for the station of cropland were also stored in the excel.

On 25 July 2012, a Leica ALS70 airborne laser scanner boarded on the Y-12 aircraft was used to obtain the point cloud data. Leica ALS70 airborne laser scanner has unlimited numbers of returns intensities measurements including the first, second, third return intensities. The wavelength of laser light is 1064 nm. The absolute flight altitude is 4800 m with the point cloud density 1 points per square meter. Aerial LiDAR- DSM was obtained through parameter calibration, automatic classification of point cloud density and manual editing.

This dataset includes data recorded by the Cold and Arid Research Network of Lanzhou university obtained from an observation system of Meteorological elements gradient of Dunhuang Station from January 1 to December 31, 2018. The site (93.708° E, 40.348° N) was located on a wetland in the Dunhuang west lake, Gansu Province. The elevation is 990 m. The installation heights and orientations of different sensors and measured quantities were as follows: air temperature and humidity profile (4m and 8 m, towards north), wind speed and direction profile (windsonic; 4m and 8 m, towards north), air pressure (1 m), rain gauge (4 m), infrared temperature sensors (4 m, towards south, vertically downward), soil heat flux (-0.05 and -0.1m ), soil soil temperature/ moisture/ electrical conductivity profile (below the vegetation in the south of tower, -0.05 and -0.2 m), photosynthetically active radiation (4 m, towards south), four-component radiometer (4 m, towards south), sunshine duration sensor(4 m, towards south). The observations included the following: air temperature and humidity (Ta_4 m, Ta_8 m; RH_2 m, RH_4 m, RH_8 m) (℃ and %, respectively), wind speed (Ws_4 m, Ws_8 m) (m/s), wind direction (WD_4 m, WD_8 m) (°), air pressure (press) (hpa), precipitation (rain) (mm), four-component radiation (DR, incoming shortwave radiation; UR, outgoing shortwave radiation; DLR_Cor, incoming longwave radiation; ULR_Cor, outgoing longwave radiation; Rn, net radiation) (W/m^2), infrared temperature (IRT) (℃), photosynthetically active radiation (PAR) (μmol/ (s m-2)), soil heat flux (Gs_0.05m, Gs_0.1m) (W/m^2), soil temperature (Ts_0.05m, Ts_0.2m) (℃), soil moisture (Ms_0.05m, Ms_0.2m) (%, volumetric water content), soil conductivity (Ec_0.05m, Ec_0.2m)(μs/cm), sun time(h). The data processing and quality control steps were as follows: (1) The AWS data were averaged over intervals of 10 min for a total of 144 records per day. The data were missing during Jan. 23 to Jan. 24 because of collector failure; the data during Mar. 17 and May 24 were wrong because of the tower body tilt; The air humidity data were rejected due to program error. (2) Data in duplicate records were rejected. (3) Unphysical data were rejected. (4) The data marked in red are problematic data. (5) The format of the date and time was unified, and the date and time were collected in the same column, for example, date and time: 2018-6-10 10:30.

This dataset contains the data of the meteorological element gradient observation system of the Sidaoqiao superstation downstream of the Heihe Hydrometeorological Observation Network from January 1, 2014 to December 31, 2014. The site is located in Sidaoqiao, Dalaihu Town, Ejin Banner, Inner Mongolia. The underlying surface is Tamarix. The latitude and longitude of the observation point is 101.1374E, 42.0012N, and the altitude is 873m. The air temperature, relative humidity and wind speed sensors are respectively set at 5m, 7m, 10m, 15m, 20m and 28m, with 6 layers facing the north; the wind direction sensor is set at 15m, facing the north; the barometer is installed in the waterproof box. The tipping bucket rain gauge is installed at 28m; the four-component radiometer is installed at 10m, facing south; two infrared thermometers are installed at 10m, facing south, the probe orientation is vertically downward; two photosynthetically active radiometers are installed At 10m, facing south, and the probe is vertically upward and downward respectively; the soil moisture sensor is installed 2m on the south side of the tower body, and the soil heat flow plates (self-correcting type) (3 pieces) are buried in turn in the ground 6cm deep; The average soil temperature sensor TCAV is buried in the ground 2cm, 4cm; the soil temperature probe is buried in the ground surface 0cm and underground 2cm, 4cm, 10cm, 20cm, 40cm, 80cm, 120cm and 160cm; soil moisture sensors are buried in the underground 2cm, 4cm, 10cm, 20cm, 40cm, 80cm, 120cm and 160cm. Observed items include: wind speed (WS_5m, WS_7m, WS_10m, WS_15m, WS_20m, WS_28m) (unit: m/s), wind direction (WD_15m) (unit: degree), air temperature and humidity (Ta_5m, Ta_7m, Ta_10m, Ta_15m, Ta_20m, Ta_28m and RH_5m, RH_7m, RH_10m, RH_15m, RH_20m, RH_28m) (unit: centigrade, percentage), pressure (unit: hectopascal), precipitation (Rain) (unit: mm), four-component radiation (DR, UR, DLR_Cor, ULR_Cor, Rn) (unit: watts/square meter), surface radiation temperature (IRT_1, IRT_2) (unit: centigrade), up and down photosynthetically active radiation (PAR_U_up, PAR_U_down) (unit: micromol/square Msec), average soil temperature (TCAV) (unit: centigrade), soil heat flux (Gs_1, Gs_2, Gs_3) (unit: watts/square meter), soil moisture (Ms_2cm, Ms_4cm, Ms_10cm, Ms_20cm, Ms_40cm, Ms_80cm) , Ms_120cm, Ms_160cm) (unit: volumetric water content, percentage), soil temperature (Ts_0cm, Ts_2cm, Ts_4cm, Ts_10cm, Ts_20cm, Ts_40cm, Ts_80cm, Ts_120cm, Ts_160cm) (unit: centigrade). Processing and quality control of the observation data: (1) ensure 144 data per day (every 10 minutes), when there is missing data, it is marked by -6999; From September 8, 2014 to November 8, due to the sensor problems, the data is missing; on May 9, 2014, the soil moisture probe was re-buried, and the data before and after is inconsistent; (2) eliminate the moment with duplicate records; (3) delete the data that is obviously beyond the physical meaning or the range of the instrument; (5) the format of date and time is uniform, and the date and time are in the same column. For example, the time is: 2014-9-10 10:30; (6) the naming rules are: AWS+ site name. For hydrometeorological network or site information, please refer to Li et al. (2013). For observation data processing, please refer to Liu et al. (2011).