The dataset contains phenological camera observation data collected at the Arou Superstation in the midstream of the Heihe integrated observatory network from June 13 to November 16, 2018. The instrument was developed with data processed by Beijing Normal University. The phenomenon camera integrates data acquisition and data transmission functions. The camera captures high-quality data with a resolution of 1280×720 by looking-downward. The calculation of the greenness index and phenology are following 3 steps: (1) calculate the relative greenness index (GCC, Green Chromatic Coordinate, calculated by GCC=G/(R+G+B)) according to the region of interest, (2) perform gap-filling for the invalid values, filtering and smoothing, and (3) determine the key phenological parameters according to the growth curve fitting (such as the growth season start date, Peak, growth season end, etc.) There are also 3 steps for coverage data processing: (1) select images with less intense illumination, (2) divide the image into vegetation and soil, and (3) calculate the proportion of vegetation pixels in each image in the calculation area. After the time series data is extracted, the original coverage data is smoothed and filtered according to the time window specified by the user, and the filtered result is the final time series coverage. This data set includes relative greenness index (Gcc). Please refer to Liu et al. (2018) for sites information in the Citation section.

The dataset contains phenological camera observation data collected at the Arou Superstation in the midstream of the Heihe integrated observatory network from June 13 to November 16, 2018. The instrument was developed with data processed by Beijing Normal University. The phenomenon camera integrates data acquisition and data transmission functions. The camera captures high-quality data with a resolution of 1280×720 by looking-downward. The calculation of the greenness index and phenology are following 3 steps: (1) calculate the relative greenness index (GCC, Green Chromatic Coordinate, calculated by GCC=G/(R+G+B)) according to the region of interest, (2) perform gap-filling for the invalid values, filtering and smoothing, and (3) determine the key phenological parameters according to the growth curve fitting (such as the growth season start date, Peak, growth season end, etc.) There are also 3 steps for coverage data processing: (1) select images with less intense illumination, (2) divide the image into vegetation and soil, and (3) calculate the proportion of vegetation pixels in each image in the calculation area. After the time series data is extracted, the original coverage data is smoothed and filtered according to the time window specified by the user, and the filtered result is the final time series coverage. This data set includes relative greenness index (GCC), phenological phase and fractional cover (FC). Please refer to Liu et al. (2018) for sites information in the Citation section.

The dataset contains phenological camera observation data collected at the Arou Superstation in the midstream of the Heihe integrated observatory network from June 13 to November 16, 2018. The instrument was developed with data processed by Beijing Normal University. The phenomenon camera integrates data acquisition and data transmission functions. The camera captures high-quality data with a resolution of 1280×720 by looking-downward. The calculation of the greenness index and phenology are following 3 steps: (1) calculate the relative greenness index (GCC, Green Chromatic Coordinate, calculated by GCC=G/(R+G+B)) according to the region of interest, (2) perform gap-filling for the invalid values, filtering and smoothing, and (3) determine the key phenological parameters according to the growth curve fitting (such as the growth season start date, Peak, growth season end, etc.) There are also 3 steps for coverage data processing: (1) select images with less intense illumination, (2) divide the image into vegetation and soil, and (3) calculate the proportion of vegetation pixels in each image in the calculation area. After the time series data is extracted, the original coverage data is smoothed and filtered according to the time window specified by the user, and the filtered result is the final time series coverage. This data set includes relative greenness index (Gcc). Please refer to Liu et al. (2018) for sites information in the Citation section.

The dataset contains phenological camera observation data collected at the Arou Superstation in the midstream of the Heihe integrated observatory network from June 13 to November 16, 2018. The instrument was developed with data processed by Beijing Normal University. The phenomenon camera integrates data acquisition and data transmission functions. The camera captures high-quality data with a resolution of 1280×720 by looking-downward. The calculation of the greenness index and phenology are following 3 steps: (1) calculate the relative greenness index (GCC, Green Chromatic Coordinate, calculated by GCC=G/(R+G+B)) according to the region of interest, (2) perform gap-filling for the invalid values, filtering and smoothing, and (3) determine the key phenological parameters according to the growth curve fitting (such as the growth season start date, Peak, growth season end, etc.) There are also 3 steps for coverage data processing: (1) select images with less intense illumination, (2) divide the image into vegetation and soil, and (3) calculate the proportion of vegetation pixels in each image in the calculation area. After the time series data is extracted, the original coverage data is smoothed and filtered according to the time window specified by the user, and the filtered result is the final time series coverage. This data set includes relative greenness index (Gcc). Please refer to Liu et al. (2018) for sites information in the Citation section.

Under the sunshine in summer, the snow covered on the ice surface melts, forming a pool of different shapes and sizes on the ice surface. The melting pool caused by the melting of sea ice surface will reduce the albedo of sea ice, which will have a significant impact on the energy balance of the polar region, increase the absorption and accelerate the melting process of sea ice. Among the factors that affect the albedo of sea ice, the melting pool is one of the most important and dramatic. With the change of climate, the melting rate of ice is faster and faster in summer. It has an important impact on the energy balance of the earth's surface. The acceleration of ice melting may also make the melting pool, an important natural phenomenon, become one of the most significant ice surface features in the Arctic sea ice melting season. The albedo of the melting pool is between sea water and sea ice. The study of the ice melting pool is also an important part of the study of the rapid change mechanism of the Arctic sea ice. Because of the similar microwave signal characteristics between sea ice melting pool and sea surface, and the uncertainty of mapping the coverage of melting pool with microwave data influenced by wind speed and other factors, the most reliable remote sensing method for the coverage of melting pool is to map the coverage of sub-pixel melting pool with MODIS and MERIS data. This data set includes the Arctic sea ice melting pool coverage retrieved by subpixel decomposition using MODIS observation data.

This map was compiled by Li Xin and others in 2008 in order to re-count the permafrost area in China and based on the analysis of the existing permafrost map in China. It consists of three parts, of which the Qinghai-Tibet Plateau part uses the simulated permafrost map of the Qinghai-Tibet Plateau (Nanzhuo Copper, 2002), the northeast part comes from the "14 million map of China's Glacier, Frozen Soil and Desert" (Institute of Environment and Engineering in Cold and Arid Regions, Chinese Academy of Sciences, 2006), and the other part uses the map of China's permafrost zoning and types (1: 10 million) (Zhou Youwu and others, 2000). More Information References (Institute of Environment and Engineering in Cold and Arid Regions, Chinese Academy of Sciences, 2006; Nanzhuo Copper, 2002; Zhou Youwu et al., 2000; Li et al, 2008）。

The time coverage of this data is (1961-1990). The station data set includes 222 stations of precipitation data and 202 stations of temperature data. In order to fill the meteorological data in the surrounding area of ​​Xinjiang in the study area, this data set uses the Central Asia Temperature and Precipitation Data (1879-2003), and some site data of Pakistan, Afghanistan, Mongolia (Global Historical Climate Network) and CRU dataset, in addition to the Xinjiang Meteorological Data Set, Qinghai, and Gansu Daily Data. There is a large amount of missing data in the used dataset, which will affect the accuracy of the grid data generated by the extrapolation method. Therefore, this article deletes sites with consecutive missing years, and uses sites adjacent to the site to replace missing sites with fewer years (less than 3 years). For sites where the spatial distribution of sites is too sparse, BP neural network is used to fit and reconstruct sites with severely missing data, such as Tazhong (51747), Andi Township (51848), and Hangya (51915). Based on the pre-processed data, the interpolation method of this data set is the Cressman objective analysis method. The monthly average temperature and monthly precipitation are extrapolated to the study area, and the grid period observation data with a horizontal resolution of 0.5 ° is obtained. This data contains two files: temperature data of xinjiangtemp.nc and precipitation data of xinjiangpre2.nc.

Based on the WRF model, using ERA5 reanalysis data as the initial and boundary fields, the high-resolution low-level atmospheric structure and the earth atmosphere exchange data set of the Qinghai Tibet Plateau are preliminarily obtained by the method of dynamic downscaling. The time range of this data set is from August 1 to August 31, 2014, with a time resolution of 1 hour, a horizontal range of 25 °N-40 °N, 70oE-105oE, and a horizontal resolution of 0.05 °. The data format is NetCDF, and one file is output every hour. The file is named after the date. The lower atmospheric structure data includes temperature, relative humidity, water vapor mixing ratio, potential height, meridional wind and latitudinal wind meteorological elements, with 34 isobaric surfaces in the vertical direction; the surface air exchange data set includes the upward / downward short wave radiation, upward / downward long wave radiation, surface sensible heat and flux, 2m air temperature and water vapor mixing ratio, 10m wind, etc. The data set can provide data support for the study of weather process and climate environment in the Tibetan Plateau.

This dataset contains the flux measurements from site No.16 eddy covariance system (EC) in the flux observation matrix from 6 June to 17 September, 2012. The site (100.36411° E, 38.84931° N) was located in a cropland (maize surface) in Daman irrigation district, which is near Zhangye, Gansu Province. The elevation is 1564.31 m. The EC was installed at a height of 4.9 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 (Gill&Li7500) was 0.2 m. Raw data acquired at 10 Hz were processed using the Eddypro post-processing software (Li-Cor Company, http://www.licor.com/env/products/ eddy_covariance/software.html), including spike detection, lag correction of H2O/CO2 relative to the vertical wind component, sonic virtual temperature correction, angle of attack 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/m^2), 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.

This dataset contains the flux measurements from site No.7 eddy covariance system (EC) in the flux observation matrix from 29 May to 18 September, 2012. The site (100.36521° E, 38.87676° N) was located in a cropland (maize surface) in Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1556.39 m. The EC was installed at a height of 3.8 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.15 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/m^2), 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 data set includes meteorological data from the Ngari Desert Observation and Research Station from 2009 to 2017. It includes the following basic meteorological parameters: temperature (1.5 m from the ground, once every half hour, unit: Celsius), relative humidity (1.5 m from the ground, once every half hour, unit: %), wind speed (1.5 m from the ground, once every half hour, unit: m/s), wind direction (1.5 m from the ground, once every half hour, unit: degrees), atmospheric pressure (1.5 m from the ground, once every half hour, unit: hPa), precipitation (once every 24 hours, unit: mm), water vapour pressure (unit: kPa), evaporation (unit: mm), downward shortwave radiation (unit: W/m2), upward shortwave radiation (unit: W/m2), downward longwave radiation (unit: W/m2), upward longwave radiation (unit: W/m2), net radiation (unit: W/m2), surface albedo (unit: %). The temporal resolution of the data is one day. The data were directly downloaded from the Ngari automatic weather station. The precipitation data represent daily precipitation measured by the automatic rain and snow gauge and corrected based on manual observations. The other observation data are the daily mean value of the measurements taken every half hour. Instrument models of different observations: temperature and humidity: HMP45C air temperature and humidity probe; precipitation: T200-B rain and snow gauge sensor; wind speed and direction: Vaisala 05013 wind speed and direction sensor; net radiation: Kipp Zonen NR01 net radiation sensor; atmospheric pressure: Vaisala PTB210 atmospheric pressure sensor; collector model: CR 1000; acquisition interval: 30 minutes. The data table is processed and quality controlled by a particular person based on observation records. Observations and data acquisition are carried out in strict accordance with the instrument operating specifications, and some data with obvious errors are removed when processing the data table.

CMADS V1.0(The China Meteorological Assimilation Driving Datasets for the SWAT model Version 1.0)Version of the data set introduces the technology of STMAS assimilation algorithm . It was constructed using multiple technologies and scientific methods, including loop nesting of data, projection of resampling models, and bilinear interpolation. The CMADS series of datasets can be used to drive various hydrological models, such as SWAT, the Variable Infiltration Capacity (VIC) model, and the Storm Water Management model (SWMM). It also allows users to conveniently extract a wide range of meteorological elements for detailed climatic analyses. Data sources for the CMADS series include nearly 40,000 regional automatic stations under China’s 2,421 national automatic and business assessment centres. This ensures that the CMADS datasets have wide applicability within the country, and that data accuracy was vastly improved. The CMADS series of datasets has undergone finishing and correction to match the specific format of input and driving data of SWAT models. This reduces the volume of complex work that model builders have to deal with. An index table of the various elements encompassing all of East Asia was also established for SWAT models. This allows the models to utilize the datasets directly, thus eliminating the need for any format conversion or calculations using weather generators. Consequently, significant improvements to the modelling speed and output accuracy of SWAT models were achieved. Most of the source data in the CMADS datasets are derived from CLDAS in China and other reanalysis data in the world. The integration of air temperature, air pressure, humidity, and wind velocity data was mainly achieved through the LAPS/STMAS system. Precipitation data were stitched using CMORPH’s global precipitation products and the National Meteorological Information Center’s data of China (which is based on CMORPH’s integrated precipitation products). The latter contains daily precipitation records observed at 2,400 national meteorological stations and the CMORPH satellite’s inversion precipitation products.The inversion algorithm for incoming solar radiation at the ground surface makes use of the discrete longitudinal method by Stamnes et al.(1988)to calculate radiation transmission. The resolutions for CMADS V1.0, V1.1, V1.2, and V1.3 were 1/3°, 1/4°, 1/8°, and 1/16°, respectively. In CMADS V1.0 (at a spatial resolution of 1/3°), East Asia was spatially divided into 195 × 300 grid points containing 58,500 stations. Despite being at the same spatial resolution as CMADS V1.0, CMADS V1.1 contains more data, with 260 × 400 grid points containing 104,000 stations. For both versions, the stations’ daily data include average solar radiation, average temperature, average pressure, maximum and minimum temperature, specific humidity, cumulative precipitation, and average wind velocity. The CMADS comprises other variables for any hydrological model（under 'For-other-model' folder ）: Daily Average Temperature, Daily Maximum Temperature, Daily Minimum Temperature, Daily cumulative precipitation (20-20h), Daily average Relative Humidity, Daily average Specific Humidity, Daily average Solar Radiation, Daily average Wind, and Daily average Atmospheric Pressure. Introduction to metadata of CMADS CMADS storage path description:(CMADS was divided into two datesets) 1.CMADS-V1.0\For-swat\ --specifically driving the SWAT model 2.CMADS-V1.0\For-other-model\ --specifically driving the other hydrological model(VIC,SWMM,etc.) CMADS--\For-swat-2009\ folder contain:(Station\ and Fork\) 1).Station\ Relative-Humidity-58500\ Daily average relative humidity(fraction) Precipitation-58500\ Daily accumulated 24-hour precipitation(mm) Solar radiation-58500\ Daily average solar radiation(MJ/m2) Tmperature-58500\ Daily maximum and minimum temperature(℃) Wind-58500\ Daily average wind speed(m/s) Where R, P, S, T, W+ dimensional grid number - the number of longitude grid is the station in the above five folders respectively.(Where R,P,S,T,W respective Daily average relative humidity,Daily cumulative precipitation(24h),Daily mean solar radiation(MJ/m2),Daily maximum and minimum temperature(℃) and Daily mean wind speed (m/s)) respectively.Data format is (.dbf) 2).Fork\ (Station index table over East Asia) PCPFORK.txt (Precipitation index table) RHFORK.txt (Relative humidity index table) SORFORK.txt (Solar radiation index table) TMPFORK.txt (Temperature index table) WINDFORK.txt (Wind speed index) CMADS--\For-swat-2012\ folder contain:(Station\ and Fork\) Storage structure is consistency with \For-swat- 2009\.However, all the data in this directory are only available in TXT format and can be readed by SWAT2012. 3)\For-other-model\ (Includes all weather input data required by the any hydrologic model (daily).) Atmospheric-Pressure-txt\ Daily average atmospheric pressure(hPa) Average-Temperature-txt\ Daily average temperature(℃) Maximum-Temperature-txt\ Daily maximum temperature(℃) Minimum-Temperature-txt\ Daily minimum temperature(℃) Precipitation-txt\ Daily accumulated 24-hour precipitation (mm) Relative-Humidity-txt\ Daily average relative humidity(fraction) Solar-Radiation-txt\ Daily average solar radiation(MJ/m2) Specific-Humidity-txt\ Daily average Specific Humidity(g/kg) Wind-txt\ Daily average wind speed(m/s) Data storage information: data set storage format is .dbf and .txt Other data information: Total data: 33.6GB Occupied space: 35.2GB Time: From year 2008 to year 2016 Time resolution: Daily Geographical scope description: East Asia Longitude: 60°E The most east longitude: 160°E North latitude: 65°N Most southern latitude: 0°N Number of stations: 58500 stations Spatial resolution: 1/3 * 1/3 * grid points Vertical range: None