The China Meteorological Forcing Dataset (CMFD) is a high spatial-temporal resolution gridded near-surface meteorological dataset that was developed specifically for studies of land surface processes in China. The dataset was made through fusion of remote sensing products, reanalysis dataset and in-situ observation data at weather stations. Its record starts from January 1979 and keeps extending (currently up to December 2018) with a temporal resolution of three hours and a spatial resolution of 0.1°. Seven near-surface meteorological elements are provided in CMFD, including 2-meter air temperature, surface pressure, specific humidity, 10-meter wind speed, downward shortwave radiation, downward longwave radiation and precipitation rate.
YANG Kun, HE Jie, WENJUN TANG , LU Hui, QIN Jun , CHEN Yingying, LI Xin
The Chinese regional surface meteorological element data set is a set of near-surface meteorological and environmental element reanalysis data set developed by the Qinghai-Tibet Plateau Research Institute of the Chinese Academy of Sciences. The data set is based on the existing Princeton reanalysis data, GLDAS data, GEWEX-SRB radiation data and TRMM precipitation data in the world, and is made by combining the conventional meteorological observation data of China Meteorological Administration. The temporal resolution is 3 hours and the horizontal spatial resolution is 0.1, including 7 factors (variables) including near-surface air temperature, near-surface air pressure, near-surface air specific humidity, near-surface full wind speed, ground downward short wave radiation, ground downward long wave radiation and ground precipitation rate. The physical meaning of each variable: | Meteorological Element || Variable Name || Unit || Physical Meaning | near-surface temperature ||temp|| K || instantaneous near-surface (2m) temperature | surface pressure || pres|| Pa || instantaneous surface pressure | specific humidity of near-surface air || shum || kg/ kg || instantaneous specific humidity of near-surface air | near ground full wind speed || wind || m /s || instantaneous near ground (anemometer height) full wind speed | downward short wave radiation || srad || W/m2 || 3-hour average (-1.5 HR ~+1.5 HR) downward short wave radiation | Downward Long Wave Radiation ||lrad ||W/m2 ||3-hour Average (-1.5 hr ~+1.5 hr) Downward Long Wave Radiation | precipitation rate ||prec||mm/hr ||3-hour average (-3.0 HR ~ 0.0 HR) precipitation rate For more information, please refer to the "User's Guide for China Meteorological Al Forcing Dataset" published with the data. The main changes in the latest version (01.06.0014) are: 1. Extend the data to December 2015 (except for short-wave and long-wave data, only until October 2015; the data from November to December 2015 are interpolated based on GLDAS data, and the error may be too large); 2. Set the minimum wind speed at 0.05 m/s; 3. Fixed a bug in the previous radiation algorithm to make our short wave and long wave data more reasonable in the morning and evening periods. 4. bug of precipitation data has been corrected, and the period involved in the change is 2011-2015.
YANG Kun, HE Jie
The dataset is a nearly 36-year (1983.7-2018.12) high-resolution (3 h, 10 km) global SSR (surface solar radiation) dataset, which can be used for hydrological modeling, land surface modeling and engineering application. The dataset was produced based on ISCCP-HXG cloud products, ERA5 reanalysis data, and MODIS aerosol and albedo products with an improved physical parameterization scheme. Validation and comparisons with other global satellite radiation products indicate that our SSR estimates were generally better than those of the ISCCP flux dataset (ISCCP-FD), the global energy and water cycle experiment surface radiation budget (GEWEX-SRB), and the Earth's Radiant Energy System (CERES). This SSR dataset will contribute to the land-surface process simulations and the photovoltaic applications in the future. The unit is W/㎡, instantaneous value.
TANG Wenjun
CMADS V1.1(The China Meteorological Assimilation Driving Datasets for the SWAT model Version 1.1) Version of the data set introduced the 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 (2m), air pressure, humidity, and wind speed data (10m) 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 (2m), average pressure, maximum and minimum temperature (2m), specific humidity, cumulative precipitation, and average wind speed (10m). The CMADS comprises other variables for any hydrological model(under 'For-other-model' folder): Daily Average Temperature (2m), Daily Maximum Temperature (2m), Daily Minimum Temperature (2m), Daily cumulative precipitation (20-20h), Daily average Relative Humidity, Daily average Specific Humidity, Daily average Solar Radiation, Daily average Wind (10m), 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 2m temperature(℃) Wind-58500 Daily average 10m 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 2m temperature(℃) Maximum-Temperature-txt Daily maximum 2m temperature(℃) Minimum-Temperature-txt Daily minimum 2m 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 10m wind speed(m/s) Data storage information: data set storage format is .dbf and .txt Other data information: Total data:45GB Occupied space: 50GB Time: From year 2008 to year 2014 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
Meng Xianyong, Wang Hao
The dataset of automatic meteorological observations was obtained at the A'rou freeze/thaw observation station from Jul. 25, 2008 to Dec. 31, 2009, in Wawangtan pasture (E100°28′/N38°03′, 3032.8), Daban, A'rou. The experimental area, situated in the valley highland of south Babaohe river, an upper stream branch of Heihe river, with a flat and open terrain slightly sloping from southeast to southeast and hills and mountains stretching for 3km is ideal for a horizontal homogeneous underlying surface. Observation items included multilayer (2m and 10m) of the wind speed, the air temperature and air humidity, the air pressure, precipitation, four components of radiation, the multilayer soil temperature (10cm, 20cm, 40cm, 80cm, 120cm and 160cm), soil moisture (10cm, 20cm, 40cm, 80cm, 120cm and 160cm), and soil heat flux (5cm & 15cm). The raw data were level0 and the data after basic processes were level1, in which ambiguous ones were marked; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate. Level2 or above were strongly recommended to domestic users. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
HU Zeyong, MA Mingguo, Wang Weizhen, HUANG Guanghui, Zhang Zhihui, TAN Junlei
1) The data set is composed of global atmospheric reanalysis data jointly produced by the National Centers for Environmental Prediction (NCEP) and the National Center for Atmospheric Research (NCAR). These grid data are generated by reanalysing the global meteorological data from 1948 to present by applying observation data, forecasting models and assimilation systems. The data variables include surface, near-surface (.995 sigma layer) and multiple meteorological variables in different barospheres, such as precipitation, temperature, relative humidity, sea level pressure, geopotential height, wind field, heat flux, etc. 2) The coverage time is from 1948 to 2018, and the data from 1948 to 1957 are non-Gaussian grid data. The data cover the whole world. The spatial resolution is a 2.5° latitude by 2.5° longitude grid. The vertical resolution is a 17-layer standard pressure barosphere, with layer boundaries at 1000, 925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, and 10 hPa, and 28 sigma levels. Some variables are calculated for 8 layers (omega) or 12 layers (humidity), with temporal resolutions of 6 hours, daily, monthly or a long-term monthly average (from 1981 to 2010). The daily data are obtained by averaging the daily values of 0Z, 6Z, 12Z and 18Z. 3) Missing values are assigned a value of -9.99691e+36f. The data are stored in the .nc format with the file name var.time.stat.nc, and each file includes data on latitude, longitude, time, and atmospheric variables. For detailed data specifications, please visit http://www.esrl.noaa.gov/pad/data.
National Oceanic and Atmospheric Administration, National Center for Atmospheric Research
The dataset of automatic meteorological observations was obtained at the Binggou cold region hydrometerological station (N38°04′/E100°13′), south of Qilian county, Qinghai province, from Sep. 25, 2007 to Dec. 31, 2009. The experimental area with paramo and riverbed gravel, situated in the upper stream valley of Heihe river, is ideal for the flat and open terrain and hills and mountains stretching outwards. The items were multilayer (2m and 10m) of the air temperature and air humidity, the wind speed, the air pressure, precipitation, four components of radiation, the multilayer soil temperature (5cm, 10cm, 20cm, 40cm, 80cm and 120cm), soil moisture (5cm, 10cm, 20cm, 40cm, 80cm and 120cm), and soil heat flux (5cm and 15cm). The raw data were level0 and the data after basic processes were level1, in which ambiguous ones were marked; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate. Level2 or above were strongly recommended to domestic users. The period from Sep. 25, 2007 to Mar. 12, 2008 was the pre-observing duration, during which hourly precipitation data (fragmented) and the soil temperature and soil moisture data were to be obtained. Stylized observations began from Mar. 12, 2008. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
WANG Jian, CHE Tao, MA Mingguo, Wang Weizhen, LI Hongyi, HAO Xiaohua, HUANG Guanghui, Zhang Zhihui, TAN Junlei
The dataset of automatic meteorological observations was obtained at the Yingke oasis station from Nov. 5, 2007 to Oct. 31, 2009. The observation site is located in an irrigation farmland in Yingke (E100°24′37.2″/N38°51′25.7″, 1519.1m), Zhangye city, Gansu province. The experimental area, situated in the middle stream Heihe river basin and with windbreaks space of 500m from east to west and 300m from south to north, is an ideal choice for its flat and open terrain. Observation items were multilayer (2m and 10m) of the wind speed and direction, air temperature and humidity, air pressure, precipitation, four components of radiation; the surface infrared temperature; the multilayer soil temperature (10cm, 20cm, 40cm, 80cm, 120cm and 160cm), the soil moisture (10cm, 20cm, 40cm, 80cm, 120cm and 160cm), and soil heat flux (5cm & 15cm). The raw data were level0 and the data after basic processes were level1, in which ambiguous ones were marked; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate. Level2 or above were strongly recommended to domestic users. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
MA Mingguo, Wang Weizhen, TAN Junlei, HUANG Guanghui, Zhang Zhihui
The dataset of eddy covariance observations was obtained at the A'rou freeze/thaw observation station from Jul. 14, 2008 to Dec. 31, 2010, in Wawangtan pasture (E100°28′/N38°03′, 3032.8m), Daban, A'rou. The experimental area with a flat and open terrain slightly sloping from southeast to northwest and hills and mountains stretching outwards is an ideal horizontal homogeneous underlying surface. The original observation items included the latitudinal wind speed Ux (m/s), the latitudinal wind speed Uy (m/s), the longitudinal wind speed Uz (m/s), the ultrasonic temperature Ts (°C), co2 consistency (mg/m^3), h2o consistency (g/m^3), air pressure (KPa) and the abnormal ultrasonic signal (diag_csat). The instrument height was 2.81m, the ultrasound direction was at an azimuth angle of 0°, the distance between Li7500 and CSAT3 was 30m and sampling frequency was 10HZ/s. The instrument mount was 3.15m, the ultrasound direction was at an azimuth angle of 86°, the distance between Li7500 and CSAT3 was 22cm and sampling frequency was 10HZ/s. The dataset was released at three levels: Level0 were the raw data acquired by instruments; Level1, including the sensible heat flux (Hs), the latent heat flux (LE_wpl), and co2 flux (Fc_wpl), were real-time eddy covariance output data and stored in .csv month by month; Level2 were processed data in a 30-minute cycle after outliers elimination, coordinates rotation, frequency response correction, WPL correction and initial quality control. The data were named as follows: station name +data level+data acquisition date. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide and Eddy Covariance Observation Manual.
Wang Weizhen, MA Mingguo, LI Xin, HUANG Guanghui, Zhang Zhihui, TAN Junlei
The dataset of eddy covariance observations was obtained at the Yingke Oasis station from 27 Dec. 2007 to 31 Dec. 2009. The observation site is located in an irrigation farmland in Yingke (E100°24′37.2″/N38°51′25.7″, 1519.1m), Zhangye city, Gansu province. The experimental area, situated in the middle stream Heihe river basin and with windbreaks space of 500m from east to west and 300m from south to north, is an ideal choice for its flat and open terrain. The original observation items included the latitudinal wind speed Ux (m/s), the latitudinal wind speed Uy (m/s), the longitudinal wind speed Uz (m/s), the ultrasonic temperature Ts (°C), co2 consistency (mg/m^3), h2o consistency (g/m^3), air pressure (KPa) and the abnormal ultrasonic signal (diag_csat). The instrument mount was 2.81m, the ultrasound direction was at an azimuth angle of 0°, the distance between Li7500 and CSAT3 was 30cm and the sampling frequency was 10HZ/s. The dataset was distributed at three levels: Level0 were the raw data acquired by instruments; Level1, including the sensible heat flux (Hs), the latent heat flux (LE_wpl), and co2 flux (Fc_wpl), were real-time eddy covariance output data and stored in .csv month by month; Level2 were processed data in a 30-minute cycle after outliers elimination, coordinates rotation, frequency response correction, WPL correction and initial quality control. The data files were named as follows: station name +data level+data acquisition date. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide and Eddy Covariance Observation Manual.
Liu Qiang, LIU Qinhuo, MA Mingguo, Wang Weizhen, HUANG Guanghui, Zhang Zhihui, TAN Junlei
The dataset of LST (land surface temperature) observed by the thermal camera (ThermaCAM SC2000 and ThermaCAM S60) at 24°×18° was obtained in the Yingke oasis, Huazhaizi desert steppe and Linze grassland foci experimental areas on May 20, 24,28 and 30, Jun. 1, 4, 16 and 29, Jul. 7, 8 and 11, 2008. Meanwhile, the optical photos were acquired in Yingke oasis maize field, Huazhaizi desert No. 1 and 2 plots, Huazhaizi desert maize field and Linze grassland. The dataset of ground truth measurement was synchronizing with WiDAS (Wide-angle Infrared Dual-mode line/area Array Scanner), OMIS-II, Landsat TM and ASTER.
HE Tao, KANG Guoting, REN Huazhong, YAN Guangkuo, WANG Haoxing, WANG Tianxing, LI Hua, Liu Qiang, XIA Chuanfu, ZHOU Chunyan, ZHOU Mengwei, CHEN Shaohui, YANG Tianfu
The dataset of eddy covariance observations was obtained at the Dayekou Guantan forest station (E100°15′/N38°32′, 2835m), south of Zhangye city, Gansu province, from Dec. 27, 2007 to Dec. 31, 2009. Guantan forest station was dominated by the spruce 15-20m high and the surface was covered by moss 10cm deep. All the vegetation was in good condition. The original observation items included the latitudinal wind speed Ux (m/s), the latitudinal wind speed Uy (m/s), the longitudinal wind speed Uz (m/s), the ultrasonic temperature Ts (°C), co2 consistency (mg/m^3), h2o consistency (g/m^3), air pressure (KPa) and the abnormal ultrasonic signal (diag_csat). The instrument mount-height was 20.02m, the ultrasound direction was at an azimuth angle of 74°, the distance between Li7500 and CSAT3 was 30cm and sampling frequency was 10HZ. The dataset was distributed at three levels: Level0 were the raw data acquired by instruments; Level1, including the sensible heat flux (Hs), the latent heat flux (LE_wpl), and co2 flux (Fc_wpl), were real-time eddy covariance output data and stored in .csv month by month; Level2 were processed data in a 30-minute cycle after outliers elimination, coordinates rotation, frequency response correction, WPL correction and initial quality control. The data were named as follows: station name +data level+data acquisition date. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide and Eddy Covariance Observation Manual.
LI Xin, MA Mingguo, Wang Weizhen, HUANG Guanghui, TAN Junlei, Zhang Zhihui
The dataset of automatic meteorological observations was obtained at the Dadongshu mountain snow observation station (E100°14′/N38°01′, 4101m) from Oct. 29, 2007 to Oct. 1, 2009. The experimental area with a flat and open terrain was slightly sloping from southeast to northwest. With alpine meadow and stones, and snow in autumn, winter and spring, the landscape was ideal. Observation items were multilayer (2m and 10m) of the wind speed, the air temperature and air humidity, the air pressure, rain and snow gauges, snow depth, four components of radiation, the multilayer soil temperature (5cm, 10cm, 20cm, 40cm, 80cm, and 120cm), soil moisture (5cm, 10cm, 20cm, 40cm, 80cm, and 120cm), and soil heat flux (5cm & 15cm). The raw data were level0 and the data after basic processes were level1, in which ambiguous ones were marked; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate. Level2 or above were strongly recommended to domestic users. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
WANG Jian, CHE Tao, LI Hongyi, HAO Xiaohua
The dataset of automatic meteorological observations was obtained at the Dayekou Guantan forest station (E100°15′/N38°32′, 2835m), south of Zhangye city, Gansu province, from Oct. 1, 2007 to Dec. 31, 2009. Guantan forest station was dominated by the 15-20m high spruce and the surface was covered by 10cm deep moss. All the vegetation was in good condition. Observation items were the multilayer (2m and 10m) wind speed and direction, the air temperature and moisture, rain and snow gauges, snow depth, photosynthetically active radiation, four components of radiation from two layers (, 1.68m and 19.75 m), stem sap flow, the surface temperature, the multi-layer soil temperature (5cm, 10cm, 20cm, 40cm, 80cm and 120cm),soil moisture (5cm, 10cm, 20cm, 40cm, 80cm and 120cm) and soil heat flux (5cm & 15cm). As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
MA Mingguo, Wang Weizhen, TAN Junlei, HUANG Guanghui, Zhang Zhihui
The GAME/Tibet project conducted a short-term pre-intensive observing period (PIOP) at the Amdo station in the summer of 1997. From May to September 1998, five consecutive IOPs were scheduled, with approximately one month per IOP. More than 80 scientific workers from China, Japan and South Korea went to the Tibetan Plateau in batches and carried out arduous and fruitful work. The observation tests and plans were successfully completed. After the completion of the IOP in September, 1998, five automatic weather stations (AWS), one Portable Atmospheric Mosonet (PAM), one boundary layer tower and integrated radiation observatory (Amdo) and nine soil temperature and moisture observation stations have been continuously observed to date and have obtained extremely valuable information for 8 years and 6 months consecutively (starting from June 1997). The experimental area is located in Nagqu, in northern Tibet, and has an area of 150 km × 200 km (Fig. 1), and observation points are also established in D66, Tuotuohe and the Tanggula Mountain Pass (D105) along the Qinghai-Tibet Highway. The following observation stations (sites) are set up on different underlying surfaces including plateau meadows, plateau lakes, and desert steppe. (1) Two multidisciplinary (atmosphere and soil) observation stations, Amdo and NaquFx, have multicomponent radiation observation systems, gradient observation towers, turbulent flux direct measurement systems, soil temperature and moisture gradient observations, radiosonde, ground soil moisture observation networks and multiangle spectrometer observations used as ground truth values for satellite data, etc. (2) There are six automatic weather stations (D66, Tuotuohe, D105, D110, Nagqu and MS3608), each of which has observations of wind, temperature, humidity, pressure, radiation, surface temperature, soil temperature and moisture, precipitation, etc. (3) PAM stations (Portable Automated Meso - net) located approximately 80 km north and south of Nagqu (MS3478 and MS3637) have major projects similar to the two integrated observation stations (Amdo and NaquFx) above and to the wind, temperature and humidity turbulence observations. (4) There are nine soil temperature and moisture observation sites (D66, Tuotuohe, D110, WADD, NODA, Amdo, MS3478, MS3478 and MS3637), each of which has soil temperature measurements of 6 layers and soil moisture measurement of 9 layers. (5) A 3D Doppler Radar Station is located in the south of Nagqu, and there are seven encrypted precipitation gauges in the adjacent (within approximately 100 km) area. The radiation observation system mainly studies the plateau cloud and precipitation system and serves as a ground true value station for the TRMM satellite. The GAME-Tibet project seeks to gain insight into the land-atmosphere interaction on the Tibetan Plateau and its impact on the Asian monsoon system through enhanced observational experiments and long-term monitoring at different spatial scales. After the end of 2000, the GAME/Tibet project joined the “Coordinated Enhanced Observing Period (CEOP)” jointly organized by two international plans, GEWEX (Global Energy and Water Cycle Experiment) and CL IVAR (Climate Change and Forecast). The Asia-Australia Monsoon Project (CAMP) on the Tibetan Plateau of the Global Coordinated Enhanced Observation Program (CEOP) has been started. The data set contains POP data for 1997 and IOP data for 1998. Ⅰ. The POP data of 1997 contain the following. 1. Precipitation Gauge Network (PGN) 2. Radiosonde Observation at Naqu 3. Analysis of Stable Isotope for Water Cycle Studies 4. Doppler radar observation 5. Large-Scale Hydrological Cycle in Tibet (Link to Numaguchi's home page) 6. Portable Automated Mesonet (PAM) [Japanese] 7. Ground Truth Data Collection (GTDC) for Satellite Remote Sensing 8. Tanggula AWS (D105 station in Tibet) 9. Syamboche AWS (GEN/GAME AWS in Nepal) Ⅱ. The IOP data of 1998 contain the following. 1. Anduo （1） PBL Tower, 2） Radiation, 3） Turbulence SMTMS 2. D66 （1） AWS （2） SMTMS （3） GTDC （4) Precipitation 3. Toutouhe （1） AWS （2） SMTMS （3 ）GTDC 4. D110 （1） AWS （2） SMTMS (3) GTDC (4) SMTMS 5. MS3608 （1） AWS （2） SMTMS （3） Precipitation 6. D105 （1） Precipitation (2) GTDC 7. MS3478(NPAM) （1） PAM （2） Precipitation 8. MS3637 （1） PAM （2） SMTMS （3） Precipitation 9. NODAA （1） SMTMS (2) Precipitation 10. WADD （1） SMTMS （2） Precipitation （3） Barometricmd 11. AQB （1） Precipitation 12. Dienpa (RS2) （1） Precipitation 13. Zuri （1） Precipitation （2） Barometricmd 14. Juze （1） Precipitation 15. Naqu hydrological station （1） Precipitation 16. MSofNaqu （1） Barometricmd 16. Naquradarsite （1）Radar system （2） Precipitation 17. Syangboche [Nepal] （1） AWS 18. Shiqu-anhe （1） AWS （2） GTDC 19. Seqin-Xiang （1） Barometricmd 20. NODA （1）Barometricmd （2） Precipitation （3) SMTMS 21. NaquHY （1） Barometricmd （2） Precipitation 22. NaquFx(BJ) （1） GTDC（2) PBLmd (3) Precipitation 23. MS3543 （1） Precipitation 24. MNofAmdo （1） Barometricmd 25. Mardi （1） Runoff 26. Gaize （1） AWS （2） GTDC （3） Sonde A CD of the data GAME-Tibet POP/IOP dataset cd （vol. 1) GAME-Tibet POP/IOP dataset cd （vol. 2)
MA Yaoming
The dataset of automatic meteorological observations was obtained from Jun. 1, 2008 to Dec. 31, 2009 at the Huazhaizi desert station which is located in Anyangtan (E100°19'06.9″/N38°45'54.7″), south of Zhangye city, Gansu province,. The experimental area, situated in the middle stream of Heihe river, with a flat and open terrain and sparse vegetation cover is an ideal desert observing field. Observation items included the multi-layer (2m and 10m) wind speed and direction, the air temperature, precipitation, the four components of radiation, the surface infrared temperature, the multi-layer soil temperature (5cm, 10cm, 20cm, 40cm, 80cm and 160cm), soil moisture (5cm, 10cm, 20cm, 40cm, 80cm and 160cm) and soil heat flux (5cm & 10cm). The raw data were level0 and the data after basic processes were level1; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate.. As for detailed information, please refer to “Meteorological and Hydrological Flux Data Guide".
LI Xin, XU Ziwei
The dataset of ground truth measurements synchronizing with ASTER was obtained in the Yingke oasis and Huazhaizi desert steppe foci experimental areas on May 28, 2008. Observation items included: (1) Atmospheric parameters in Huazhaizi desert No. 2 plot by CE318 (produced by CIMEL in France). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in .k7 format and can be opened by ASTPWin. ReadMe.txt is attached for detail. Processed data (after retrieval of the raw data) in Excel format are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (2) Photosynthesis by LI-6400. Raw data were archived in the user-defined format (by notepat.exe) and processed data were in Excel format. (3) Reflectance spectra in Yingke oasis maize field by ASD FieldSpec (350-2500nm, the vertical canopy observation and the transect observation) from Institute of Remote Sensing Applications (CAS), and in Huazhaizi desert No. 2 plot by ASD FieldSpec (350-1603nm, the vertical observation and the transect observation for reaumuria soongorica and the bare land) from Beijing Academy of Agriculture and Forestry Sciences. The grey board and the black and white cloth were also used for calibration spectrum. Raw data were binary files direct from ASD (by ViewSpecPro), and pre-processed data on reflectance were in Excel format. (4) Coverage fraction of maize and wheat by the self-made instrument and the camera (2.5m-3.5m above the ground) in Yingke oasis maize field. Based on the length of the measuring tape and the bamboo pole, the size of the photo can be decided. GPS date were also collected and the technology LAB was applied to retrieve the coverage of the green vegetation. Besides, such related information as the surrounding environment was also recorded. Data included the primarily measured image and final fraction of vegetation coverage. (5) the radiative temperature of maize, wheat and the bare land in Yingke oasis maize field by ThermaCAM SC2000 using ThermaCAM SC2000 (1.2m above the ground, FOV = 24°×18°),. The data included raw data (read by ThermaCAM Researcher 2001), recorded data and the blackbody calibrated data (archived in Excel format). (6) the radiative temperature by the automatic thermometer (FOV: 10°; emissivity: 0.95), 3 for maize canopy, the bare land and wheat canopy in Yingke oasis maize field, one for maize canopy in Huazhaizi desert maize field, and 2 for vegetation and the desert bare land in Huazhaizi desert No. 2 plot,at nadir at a time interval of one second. Raw data, blackbody calibrated data and processed data were all archived in Excel format. (7) Maize albedo by the shortwave radiometer in Yingke oasis maize field. R =10H (R for FOV radius; H for the probe height). Data were archived in Excel format. (8) LAI in Yingke oasis maize field. The maximum leaf length and width of each maize and wheat were measured. Data were archived in Excel format. (9) FPAR (Fraction of Photosynthetically Active Radiation) of maize and wheat by SUNSACN and the digital camera in Yingke oasis maize field. FPAR= (canopyPAR－surface transmissionPAR－canopy reflection PAR+surface reflectionPAR) /canopy PAR; APAR=FPAR* canopy PAR. Data were archived in the table format of Word. (10) The radiative temperature in Yingke oasis maize field (the transect observation), Yingke oasis wheat field (the transect observation), Huazhaizi desert maize field (the transect observation) and Huazhaizi desert No. 2 plot (the diagonal observation) by the handheld infrared thermometer (BNU and Institute of Remote Sensing Applications). Raw data (in Word format), blackbody calibrated data and processed data (in Excel format) were all archived.
CHAI Yuan, CHEN Ling, KANG Guoting, QIAN Yonggang, REN Huazhong, WANG Haoxing, WANG Jianhua, SHU Lele, LI Li, LIU Sihan, XIN Xiaozhou, ZHANG Yang, ZHOU Chunyan, ZHOU Mengwei, TAO Xin, WANG Dacheng, LI Xiaoyu, CHENG Zhanhui, YANG Tianfu, HUANG Bo, LI Shihua, LUO Zhen
The dataset of automatic meteorological observations was obtained at the Linze grassland station (E100 °04'/N39°15', 1394m) from Oct. 1, 2007 to Oct. 27, 2008. The landscape is dominated by wetland and saline land. Observation items were multilayer (2m, 4m and 10m) of the wind speed and direction, air temperature and humidity, air pressure, precipitation, four components of radiation, the surface temperature, the soil temperature (5cm, 10cm, 20cm and 40cm), and the multilayer soil temperature (2cm, 5cm and 10cm). The dataset was released at different levels: Level1 were transformed raw data and stored in .csv month by month; Level2 were processed data after correction and quality control. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
HU Zeyong, MA Mingguo, Wang Weizhen, TAN Junlei, HUANG Guanghui, Zhang Zhihui
The dataset of automatic meteorological observations was obtained at the Dayekou Maliantan grassland station (E100°18′/N38°33′, 2817m) from Nov. 2, 2007 to Dec. 31, 2009. The experimental area with a flat and open terrain was slightly sloping from southeast to northwest. The landscape was mainly grassland, with vegetation 0.2-0.5m high. Observation items were multilayer gradient (2m and 10m) of the wind speed, the air temperature and air humidity, the air pressure, precipitation, four components of radiation, the multilayer soil temperature (5cm, 10cm, 20cm, 40cm, 80cm, and 120cm), soil moisture (5cm, 10cm, 20cm, 40cm, 80cm, and 120cm), and soil heat flux (5cm & 15cm). The raw data were level0 and the data after basic processes were level1, in which ambiguous ones were marked; the data after strict quality control were defined as Level2. The data files were named as follows: station+datalevel+AMS+datadate. Level2 or above were strongly recommended to domestic users. As for detailed information, please refer to Meteorological and Hydrological Flux Data Guide.
MA Mingguo, Wang Weizhen, TAN Junlei, HUANG Guanghui, Zhang Zhihui
The dataset of LAS (Large Aperture Scintillometer: BLS450, made in Germany) observations was obtained at the A'rou freeze/thaw observation station from Mar. 11 to Jul. 11, 2008. The transmitter (E100°28′16.4″, N38°03′24.3″, 11.2m) and the receiver (E100°27′25.9″, N38°02′18.1″, 11.5m) were 2390m away from each other and the operating altitude was 9.5m. The observation item was the atmospheric refractive index structural parameters (Cn2). The transmitting frequency was 5HZ and the data were output per minute. The processed data were archived in a 30 minutes cycle. The data were named after WATER_LAS_A'rou_yyyymmdd-yyyymmdd.csv (yyyymmdd-yyyymmdd for observation time). The missing data were marked "None".
LIU Shaomin, LI Xin, XU Ziwei
The dataset of LAS (Large Aperture Scintillometer, made in Holland) observations was obtained in the Linze grassland station, Linze county (Gansu province), from May 19 to Aug. 31, 2008. The instrument was composed of the transmitter (100°04′10.4″E, 39°15′02.8″N, 9.25m), the receiver (100°03′36.8″E, 39°15′02.8″N, 9.1m) and the data acquisition system. The transmitter and the receiver were 1550m away from each other and the operating altitude was 9.2m. The observation item was natural logarithm of structural parameters of the refractive index (UCn2). The transmitting frequency was 0.5HZ. The data were named after WATER_LAS_Linze_yyyymmdd-yyyymmdd.csv (yyyymmdd-yyyymmdd for observation time). The missing data were marked "None". For more detailed information, please refer to Directions on LAS (Large Aperture Scintillometer) observations.
LIU Shaomin, LI Xin, XU Ziwei
The dataset of ground truth measurement synchronizing with PROBA CHRIS was obtained in the Yingke oasis and Huazhaizi desert steppe foci experimental areas on Jun. 22, 2008. Observation items included: (1) Albedo by the shortwave radiometer in Huazhaizi desert No. 2 plot. R =10H (R for FOV radius; H for the probe height). Data were archived in Excel format. (2) BRDF of maize in Yingke oasis maize field by ASD (350-2 500 nm) from Beijing University and the observation platform of BNU make. The maximum height of the platform was 5m above the ground with the azimuth 0~360° and the zenith angle -60°~60°; BRDF in Huazhaizi desert No. 2 plot by ASD from Institute of Remote Sensing Applications (CAS) and the observation platform of its own make, whose maximum height was 2m above the ground with the zenith angle -70°~70°. Raw data were binary files direct from ASD (by ViewSpecPro), and pre-processed data on reflectance were in Excel format. (3) Atmospheric parameters in Huazhaizi desert No. 2 plot by CE318 (produced by CIMEL in France). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in .k7 format and can be opened by ASTPWin. ReadMe.txt is attached for detail. Processed data (after retrieval of the raw data) in Excel format are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number.
CHEN Ling, GUO Xinping, REN Huazhong, ZOU Jie, LIU Sihan, ZHOU Chunyan, FAN Wenjie, TAO Xin
The dataset of ground truth measurement synchronizing with Envisat ASAR was obtained in the arid region hydrological experimental area on Sep. 19, 2007 during the pre-observation period. One scene of Envisat ASAR image was captured on Sep. 19. The data were in AP mode and VV/VH polarization combinations, and the overpass time was approximately at 11:29 BJT. Those provide reliable ground data for remote sensing retrieval and validation of soil moisture from Envisat ASAR image. Observation items included: (1) soil moisture measured by the cutting ring method in Linze reed land, Zhangye farmland, Zhangye gobi, Linze maize land, Linze alfalfa land, Zhangye weather station, and Linze wetland. (2) GPS measured by GARMIN GPS 76 (3) vegetation measurements including the vegetation height, the green weight, the dry weight, the sampling method, and descriptions on the land type, uniformity and dry and wet conditions (4) atmospheric parameters at Daman Water Management office measured by CE318 (produced by CIMEL in France). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in .k7 and can be opened by ASTPWin. ReadMetext files (.txt) is attached for detail. Processed data (after retrieval of the raw data) archived as Excel files are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (5) roughness measured by the roughness plate together with the digital camera. The coordinates of the sample would be got with the help of ArcView; and after geometric correction, surface height standard deviation (cm) and correlation length (cm) could be acquired based on the formula listed on pages 234-236, Microwave Remote Sensing (Vol. II). The roughness data were initialized by the sample name, which was followed by the serial number, the name of the file, standard deviation and correlation length. Each text files (.txt) file is matched with one sample photo and standard deviation and correlation length represent the roughness. In addition, the length of 101 radius is also included for further checking.
CHE Tao, LI Xin, BAI Yunjie, DING Songchuang, GAO Song, HAN Xujun, HAO Xiaohua, LI Hongyi, LI Zhe, LIANG Ji, PAN Xiaoduo, QIN Chun, RAN Youhua, WANG Xufeng, WU Yueru, YAN Qiaodi, ZHANG Lingmei, FANG Li, LI Hua, Liu Qiang, Wen Jianguang, MA Hongwei, YAN Yeqing, YUAN Xiaolong
This dataset contains the flux measurements from the Zhangye wetland station eddy covariance system (EC) in the flux observation matrix from 25 June to 26 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.00 m. The EC was installed at a height of 5.2 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&Li7500A) was 0.25 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.
LIU Shaomin, LI Xin, XU Ziwei
The dataset of TIR spectral emissivity was obtained in the arid region hydrology experiment area and A'rou foci experiment area. Observations were by: (1) Spectral emissivity obtained from 102F at 2-25um in cooperation with the handheld infrared thermometer (BNU) for the surface radiative temperature and one au-plating board for downward atmospheric radiation. The radiative transfer equation and TES methods were applied to retrieve emissivity. The grassland and the concrete floor were measured on May, 27, 2008, the wheat field and the maize field at ICBC resort on May, 29, 2008, the concrete floor (multiangle measurements) at ICBC resort on Jun. 3, 2008, the bare soil and the maize leaf in Yingke oasis maize field on Jun. 22, 2008, the maize and wheat canopy in Yingke oasis maize field on Jun. 23, 2008, the rape field in Biandukou experimental area on Jun. 24, 2008, the alfalfa, the saline land, the grassland and the barley land on Jun. 26, 2008, the wheat field and the maize field in Yingke oasis maize field on Jun. 29, 2008, the desert bare land and vegetation (Reaumuria soongorica) in No. 2 Huazhaiai desert plot on Jun. 30, 2008, the rape field and the grassland in Biandukou experimental area on Jul. 6, 2008, and the grassland and the bare land (multiangle) in A'rou experimental area on Jul. 14, 2008. The cold blackbody calibration (*.CBX/*.CBB), the warm blackbody calibration (*.WBX/*.WBB), the ground objects measurements (*.SAX), au-plating board measurements, and the downward atmospheric radiation (*.DWX) were all needed during observation. Moreover, the spectral radiance and emissivity were also archived. The response function of various bands could be acquired by 102F. And then emissivity of 2-25um could be retrieved. Two results of emissivity were developed: one was direct from 102F and the other was retrieved by ISSTES (Iterative spectrally smooth temperature-emissivity separation). Spectral resolution for raw data and proprecessed data was 4cm-1. (2) Spectral emissivity obtained from BOMAN at 2 -13μm in cooperation with the blackbody barrel and the blackbody from Institute of Remote Sensing Applications and the blackbody (BNU). The desert was measured on Jun. 30 and Jul. 1, 2008, A'rou foci experimental area on Jul. 14, 2008, indoor observations on the deep and shallow layer soil, vegetation, small stones, two maize plants from Yingke No.2 (YKYZYMD02) field and one maize plant and bare land from No. 3 (YKYZYMD03)field on on Jul. 16, 2008, Linze experimental area on Jul. 17, 2008, and gobi on Jul. 18, 2008. The sample site, coordinates, time and photos were all archived. During each observation, BOMAN was preheated and the blackbody was set at the predicted target temperature, which would be changed after the infrared radiation of the blackbody was measured by BOMAN. And then the target infrared radiation, the downward atmospheric radiation (reflected by the au-plating board) and the infrared radiation of the blackbody would be measured one by one. Raw data were archived in Igm, and after processed by FTSW500, the result was Rad (radiation). Finally, Rad would be changed into txt files by Matlab programs.
REN Huazhong, CHEN Ling, YAN Guangkuo, DU Yongming, LI Hua, LIU Yani, WANG Heshun, XIAO Qing, ZHOU Chunyan
The dataset of ground truth measurement synchronizing with the airborne imaging spectrometer (OMIS-II) mission was obtained in the Yingke oasis and Huazhaizi desert steppe foci experimental areas on Jun. 4, 2008. Observation items included: (1) ground object reflectance spectra of maize and wheat in Yingke oasis maize field by ASD FieldSpec (350~2500 nm, the vertical canopy observation and the transect observation) from Institute of Remote Sensing Applications (CAS); and of the black and white cloth, the water body, vegetation and the cement floor in the resort calibration site by ASD (350-2500nm, fixed points observation) from BNU. Raw data were binary files direct from ASD (by ViewSpecPro), and pre-processed data on reflectance were in Excel format. (2) The radiative temperature in Yingke oasis maize field (the transect observation), Yingke oasis wheat field (the transect observation), the maize field (intensive) near the resort (the transect observation) and Huazhaizi desert No. 1 plot (the diagonal and the fixed point observation) by the handheld infrared thermometer (emissivity: 1.00). As for the fixed point observation, 25 corner points were chosen in the plot of 30m×30m, and at each point, the bare land was measured twice and the vegetation once. Raw data (in Word format), blackbody calibrated data and processed data (in Excel format) were all archived. (3) Atmospheric parameters on the ICBC resort office roof by CE318 (produced by CIMEL in France) from Institute of Remote Sensing Applications. The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1640nm, 1020nm, 936nm, 870nm, 670nm, 550nm, 440nm, 380nm and 340nm were all acquired by CE318. Those data include the raw data in .k7 format and can be opened by ASTPWin. ReadMe.txt is attached for detail. Processed data (after retrieval of the raw data) in Excel format are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (4) Photosynthesis of wheat and maize by LI6400 in Yingke oasis maize field, carried out according to WATER specifications. Raw data were archived in the user-defined format (by notepat.exe) and processed data were in Excel format. (5) the radiative temperature vegetation (Reaumuria soongorica) and the bare land in Huazhaizi desert No. 1 plot by ThermaCAM SC2000 ( (1.2m above the ground, FOV = 24°×18°),. The data included raw data (read by ThermaCAM Researcher 2001), recorded data and the blackbody calibrated data (archived in Excel format). (6) the radiative temperature by the automatic thermometer at nadir in Yingke oasis maize field (2 from BNU, FOV: 10°; emissivity: 0.95, at intervals of 1s, set above the maize canopy and the bare land between ridges and the third from Institute of Remote Sensing Applications, emissivity: 1.0, at intervals of 0.05s, set above the maize canopy), Yingke wheat field (one set above the wheat canopy), Huazhaizi desert No. 1 plot (one set above the barley canopy), and in the resort calibration site (one for the cement floor). Raw data, blackbody calibrated data and processed data were all archived in Excel format. (7) Wheat albedo by the shortwave radiometer in Yingke oasis maize field. R =10H (R for FOV radius; H for the probe height). Data were archived in Excel format. (8) Wheat FPAR (Fraction of Photosynthetically Active Radiation) by SUNSACN and the digital camera in Yingke oasis maize field. FPAR= (canopyPAR－surface transmissionPAR－canopy reflection PAR+surface reflectionPAR) /canopy PAR; APAR=FPAR* canopy PAR. Data were archived in the table format of Word. (9) LAI in Yingke oasis maize field. The maximum leaf length and width of each maize and wheat were measured. Data were from Jun. 6, 2008, archived in Excel format.
CHEN Ling, REN Huazhong, ZHOU Hongmin, CAO Yongpan, SHU Lele, WU Yueru, XU Zhen, LI Li, LIU Sihan, XIA Chuanfu, XIN Xiaozhou, ZHOU Chunyan, ZHOU Mengwei, FAN Wenjie, TAO Xin, FENG Lei, LIANG Wenguang, YU Fan, WANG Dacheng, YANG Guijun, LI Xiaoyu, Liu Liangyun
This dataset contains the spectra of white cloth and black cloth obtained in the simultaneous time during the airborn remote sensing which supports the airboren data preprocessing as CASI, SASI and TASI , and the spetra of the typical targets in the middle reaches of the Heihe River Basin. Instruments: SVC-HR1024 from IRSA, ASD Field Spec 3 from CEODE, Reference board Measurement method: the spectra radiance of the targets are vertically measured by the SVC or ASD; before and after the target, the spectra radiance of the reference board is measured as the reference. This dataset contains the spectra recorded by the SVC-HR1024 ( in the format of .sig which can be opened by the SVC-HR1024 software or by the notepad ) and the ASD (in the format of .asd), the observation log (in the format of word or excel), and the photos of the measured targets. Observation time: 15-6-2012, the spectra of typical targets in the EC matrix using SVC 16-6-2012, the spectra of typical targets in the wetland by SVC 29-6-2012, the spectra of typical vegetation and soil in Daman site and Gobi site by ASD 29-6-2012, the spectra of white cloth and black cloth by ASD which is simultaneous with the airborne CASI data 30-6-2012, the spectra of vegetation and soil in the desert by ASD 5-7-2012, the spectra of white cloth and black cloth by ASD which is simultaneous with the airborne CASI data 7-7-2012, the spectra of corn in the Daman site for the research of daily speral variation. 8-7-2012, the spectra of white cloth and black cloth by ASD which is simultaneous with the airborne CASI data 8-7-2012, the spectra of corn in the Daman site by ASD for the research of daily speral variation 9-7-2012, the spectra of corn in the Daman site by ASD for the research of daily speral variation 10-7-2012, the spectra of corn in the Daman site by ASD for the research of daily speral variation 11-7-2012, the spectra of corn in the Daman site by ASD for the research of daily speral variation. The time used in this dataset is in UTC+8 Time.
XIAO Qing, MA Mingguo
The dataset of diurnal change of FPAR observations was obtained by the quantum meter in the Linze grassland foci experimental area. Incident and reflected radiation of canopy, and land surface in reed, saline grass, alfalfa, cumin and barley were measured and diurnal changes of PAR and Fpar were also acquired. Observations were carried out: In plot E (barley) and cumin field on Jun. 6, 2008; plot D (alfalfa) and plot E on Jun. 11; plot D and E on Jun. 15; plot E on Jun. 16; plot A (reed) on Jun. 20; plot B (saline) on Jun. 22; plot D and E on Jun. 23; plot B (saline) on Jun. 24; plot A and plot E on Jun. 29. 14 Excel files, one Word and one .TXT were archived. See Water: The dataset of setting of the sampling plots and stripes in the Linze grassland foci experimental area for more information.
CAO Yongpan, CHAO Zhenhua, GE Chunmei, HU Xiaoli, HUANG Chunlin, LIANG Ji, NIAN Yanyun, WANG Shuguo, WANG Xufeng, WU Yueru, LI Xiaoyu
The dataset of ground truth measurement synchronizing with the airborne WiDAS mission was obtained in the Yingke oasis and Huazhaizi desert steppe foci experimental areas on Jun. 1, 2008. WiDAS, composed of four CCD cameras, one mid-infrared thermal imager (AGEMA 550), and one infrared thermal imager (S60), can acquire CCD, MIR and TIR band data. The simultaneous ground data included: (1) The radiative temperature of maize, wheat and the bare land in Yingke oasis maize field and Huazhaizi desert No. 1 plot by ThermaCAM SC2000 (1.2m above the ground, FOV = 24°×18°). The data included raw data (read by ThermaCAM Researcher 2001), recorded data and the blackbody calibrated data (archived in Excel format). (2) The radiative temperature by the automatic thermometer (FOV: 10°; emissivity: 1.0; from Institute of Remote Sensing Applications), observing straight downwards at intervals of 1s in Yingke oasis maize field. Raw data, blackbody calibrated data and processed data were all archived in Excel format. (3) FPAR (Fraction of Photosynthetically Active Radiation) of maize and wheat by SUNSACN and the digital camera in Yingke oasis maize field. FPAR= (canopyPAR－surface transmissionPAR－canopy reflection PAR+surface reflectionPAR) /canopy PAR; APAR=FPAR* canopy PAR. Data were archived in Excel format. (4) The reflectance spectra by ASD in Yingke oasis maize field (350-2500nm , from BNU, the vertical canopy observation and the transect observation), and Huazhaizi desert No. 1 plot (350-2500nm , from Cold and Arid Regions Environmental and Engineering Research Institute, CAS, the NE-SW diagonal observation at intervals of 30m). The data included raw data (in .doc format), recorded data and the blackbody calibrated data (in Excel format). (5) Maize albedo by the shortwave radiometer in Yingke oasis maize field. R =10H (R for FOV radius; H for the probe height). Data were archived in Excel format. (6) The radiative temperature by the handheld radiometer in Yingke oasis maize field (from BNU, the vertical canopy observation, the transect observation and the diagonal observation), Yingke oasis wheat field (only for the transect temperature), and Huazhaizi desert No. 1 plot (the NE-SW diagonal observation). Besides, the maize radiative temperature and the physical temperature were also measured both by the handheld radiometer and the probe thermometer in the maize plot of 30m near the resort. The data included raw data (in .doc format), recorded data and the blackbody calibrated data (in Excel format). (7) Atmospheric parameters on the playroom roof at the resort by CE318 (produced by CIMEL in France). The underlying surface was mainly composed of crops and the forest (1526m high). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in .k7 format and can be opened by ASTPWin. ReadMe.txt is attached for detail. Processed data (after retrieval of the raw data) in Excel format are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (8) Narrow channel emissivity of the bare land and vegetation by the W-shaped determinator in Huazhaizi desert No. 1 plot. Four circumstances should be considered for emissivity, with the lid plus the au-plating board, the au-plating board only, the lid only and without both. Data were archived in Word.
CHEN Ling, HE Tao, REN Huazhong, REN Zhixing, YAN Guangkuo, ZHANG Wuming, XU Zhen, LI Xin, GE Yingchun, SHU Lele, JIANG Xi, HUANG Chunlin, GUANG Jie, LI Li, LIU Sihan, WANG Ying, XIN Xiaozhou, ZHANG Yang, ZHOU Chunyan, LIU Xiaocheng, TAO Xin, CHEN Shaohui, LIANG Wenguang, LI Xiaoyu, CHENG Zhanhui, Liu Liangyun, YANG Tianfu
The dataset of mobile meteorological station observations was obtained in the foci experiment area from March to April, 2008. To synergize the very high resolution airborne remote sensing and ground-based measurements, 11 mobile observations, including meteorological stations (for meteorological data) and GPS (for observation sites), were carried out in Binggou, A'rou and Biandukou. The items included the wind speed and direction at 3.03m (the truck height 1.84m plus the vane height 1.19m), the air temperature and humidity at 3.04m (the truck height 1.84m plus the vane height 1.2m), the surface temperature (the truck height 1.84m plus 1.06m) and the total radiation (the truck height 1.84m plus 1.39m). The observation sites and time were as follows: Dadongshu mountain pass-A'rou 15-3-2008 Biandukou-Qilian 18-3-2008 A'rou-Biandukou 19-3-2008 Qilian-Minle 20-3-2008 Mingle-Zhangye 21-3-2008 Binggou-Dadongshu mountain pass 22-3-2008 Binggou-Dadongshu mountain pass 24-3-2008 Binggou-Dadongshu mountain pass 29-3-2008 Binggou-Dadongshu mountain pass 30-3-2008 Qilian-A'rou 31-3-2008 A'rou 01-4-2008 The data were named after WATER_Mobile_ AWS_yyyymmdd (yyyymmdd for observation time).
HU Zeyong, GU Lianglei, SUN Fanglei, GAO Hongchun, MA Weiqiang, LI Maoshan, ZHOU Xiuyun, HOU Xuhong, REN Yanxia, MA Xiaowei
The dataset of ground truth measurement synchronizing with PROBA CHRIS was obtained in the Yingke oasis and Huazhaizi desert steppe foci experimental areas on Jul. 1, 2008. Observation items included: (1) FPAR (Fraction of Photosynthetically Active Radiation) of maize and wheat by SUNSACN and the digital camera in Yingke oasis maize field. FPAR= (canopyPAR－surface transmissionPAR－canopy reflection PAR+surface reflectionPAR) /canopy PAR; APAR=FPAR* canopy PAR. Data were archived in the table format of Word. (2) BRDF of maize by ASD (350～2 500 nm) from Institute of Remote Sensing Applications (CAS) and the self-made multi-angluar observation platform of BNU make in Yingke oasis maize field. The maximum height of the platform was 5m above the ground with the azimuth 0~360° and the zenith angle -60°~60°. An automatic thermometer was attached to the platform for the multiangle radiative temperature. Raw data were binary files direct from ASD (by ViewSpecPro), and pre-processed data on reflectance were in Excel. (3) The radiative temperature of the maize canopy by the automatic thermometer (emissivity: 0.95),at a hight of 50cm from the crown in Yingke oasis maize field. Raw data, blackbody calibrated data and processed data were all archived in Excel format. (4) Atmospheric parameters at the resort by CE318 (produced by CIMEL in France). The total optical depth, aerosol optical depth, Rayleigh scattering coefficient, column water vapor in 936 nm, particle size spectrum and phase function were then retrieved from these observations. The optical depth in 1020nm, 936nm, 870nm, 670nm and 440nm were all acquired by CE318. Those data include the raw data in k7 format and can be opened by ASTPWin. ReadMe.txt is attached for details. Processed data (after retrieval of the raw data) in Excel format are on optical depth, rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. (5) The multiangle radiative temperature by the automatic thermometer (emissivity: 1.0) attached on the observation platform, at an interval of 0.05s. The data were archived in .txt files (.dat format). The first seven lines were the header file, including acquisition date, time, and intervals; besides, Time (starting time), TObj (target temperature), Tint (the interior temperature of the probe), TBox (the temperature of the box) and Tact (the actual temperature calculated from the given emissivity) were also listed.
CHEN Ling, REN Huazhong, XIAO Yueting, SU Gaoli, WU Mingquan, WU Chaoyang, XIA Chuanfu, ZHOU Chunyan, ZHOU Mengwei, SHEN Xinyi, YANG Guijun
This data set contains the observation data of Zhangye National Climate Observatory from 2008 to 2009. The station is located in Zhangye, Gansu Province, with longitude and latitude of 100 ° 17 ′ e, 39 ° 05 ′ N and altitude of 1456m. The observation items include: atmospheric wind temperature and humidity gradient observation (2cm, 4cm, 10cm, 20m and 30m), wind direction, air pressure, photosynthesis effective radiation, precipitation, radiation four components, surface temperature, multi-layer soil temperature (5cm, 10cm, 15cm, 20cm and 40cm), soil moisture (10cm, 20cm, 50cm, 100cm and 180cm) and soil heat flux (5cm, 10cm and 15cm). Please refer to the instruction document published with the data for specific header and other information.
Zhangye city meteorological bureau
This data set contains the eddy related data of Zhangye National Climate Observatory from 2008 to 2009. The station is located in Zhangye, Gansu Province, with longitude and latitude of 100 ° 17 ′ e, 39 ° 05 ′ N and altitude of 1456m. For more information, see the documentation that came with the data.
Zhangye city meteorological bureau
This dataset includes the emissivity spectrum (8-14 µm) of typical ground objects in Zhangye City, Zhangye airport, desert and farmland at Wuxing experiment area. The data was measured by the BOMEM MR304 FTIR (Fourier Transform Infrared Spectrometer). A. Objective The objective of the thermal infrared (TIR) spectrum measurement lies in: Radiometric calibration for the airborne TIR sensor, land surface emissivity products validation and collecting typical surface spectrum working as priori knowledge in land surface temperature inversion and ecological and hydrological models. B. Instruments and theory Instruments: BOMEM MR304 FTIR, Mikron M340 blackbody, BODACH BDB blackbody, diffused golden plate, Fluke 50-series II thermometer Measurement theory: The target radiance is directly measured by the MR304 FTIR under clear-sky condition while the atmospheric downward radiance is obtained through a diffused golden plate, and emissivity is retrieved by the Iterative Spectrally Smooth Temperature and Emissivity Separation (ISSTES) algorithm C. Experiment site and targets 29-5-2012: Stone bricks, grassland and asphalt, etc at square of Zhangye. 20-6-2012: Roof of the building in Zhangye, water and sand sample collected from the desert, etc. 30-6-2012: Cement road at Zhangye airport, desert around the Zhangye airport. 3-7-2012: Corn leaves, soil and road in the farmland at Wuxing village, Zhangye City. 4-7-2012: Corn leaves, wheat canopy at Xiaoman town, Zhangye City. 10-7-2012: Bricks of Runquanhu park, Zhangye City. 13-7-2012: Corn leaves and other plants at Wuxing village, Zhangye City. D. Data processing The original data collected by BOMEM FTIR is firstly calibrated using the calibration data and get the radiance spectrum of the targets and sky (*.rad), then, the radiance data is converted to the easy readably text file (ASCII format). The time used in this dataset is in UTC+8 Time.
MA Mingguo, XIAO Qing
This dataset contains the flux measurements from the large aperture scintillometer (LAS) at A’rou Superstation in the hydrometeorological observation network of Heihe River Basin between 14 October, 2012, and 31 December, 2013. There were two types of LASs at A’rou Superstation: German BLS450 and China zzlas. The north tower was set up with the zzlas receiver and the BLS450 transmitter, and the south tower was equipped with the zzlas transmitter and the BLS450 receiver. Zzlas has been in use since 14 October, 2012, and the observation period of BLS450 was from 9 August to 10 December, 2013. The site (north: 100.467° E, 38.050° N; south: 100.450° E, 38.033° N) was located in Caodaban village of A’rou town in Qilian county, Qinghai Province. The underlying surface between the two towers was alpine meadow. The elevation is 3033 m. The effective height of the LASs was 9.5 m, and the path length was 2390 m. The data were sampled at 5 Hz and 1 Hz intervals for BLS450 and zzlas, respectively, and then averaged over 1 min. 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>7.25E-14, zzlas: Cn2>7.84E-14). (2) The data were rejected when the demodulation signal was small (BLS450: Average X Intensity<1000; zzlas: Demod>-20 mv). (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 zzlas, respectively. Several instructions were included with the released data. (1) The data were primarily obtained from BLS450 measurements, and missing flux measurements from the BLS450 instrument were substituted with measurements from the zzlas instrument. The missing data were denoted by -6999. Due to the drift of the zzlas signal, data from 10 November to 23 November, 2012, and 14 March to 10 April, 2013, were excluded. Due to the LAS tower’s lean, the data from 10 April to 31 May, 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.
LI Xin, CHE Tao, XU Ziwei, ZHANG Yang, TAN Junlei
This dataset contains the flux measurements from site No.10 eddy covariance system (EC) in the flux observation matrix from 4 June to 17 September, 2012. The site (100.39572° E, 38.87567° N) was located in a cropland (maize surface) in Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1534.73 m. The EC was installed at a height of 4.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&Li7500) 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/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.
LIU Shaomin, LI Xin, XU Ziwei
This dataset includes the emissivity spectrum of typical ground objects in middle researches of the Heihe river basin. This dataset was acquired in oasis, desert, Gobi and wetland of experiment area. Time range starts from 2012-05-25 to 2012-07-18 (UTC+8). Instrument: MODEL 102F PORTABLE FTIR (Fourier Transform Infrared Spectrometer), Handheld infrared thermometer. Measurement methods: at the first step, measure the thermal radiance of cold blackbody, warm blackbody, sample and gold plate (Downwelling Radiance). The radiance of cold blackbody and warm blackbody was used to calibrate the instrument, and eliminate the “noise” caused by the device itself. The retrieval of emissivity and temperature was then performed using iterative spectrally smooth temperature-emissivity separation (ISSTES) algorithm. The retrieved emissivity spectrum range from 8 to 14 μm, with spectral resolution of 4cm-1. Dataset contains the original recorded spectra (in ASCII format) and the log files (in doc format). The processed data are emissivity curves (ASCII) that ranged from 8 to 14 μm, and the temperatures of samples. Thermal photos of the sample, digital photo of the scene and the object are recorded in some cases.
MA Mingguo
The dataset of sun photometer observations was obtained in the Binggou watershed foci experimental areas (N38°04′1.4″/E100°13′15.6″, 3414.41m) from Mar. 15 to Apr. 2, 2008 (to be specific, the daytime of 15-03-2008, 16-03-2008, 17-03-2008, 18-03-2008, 19-03-2008, 21-03-2008, 22-03-2008, 23-03-2008, 24-03-2008, 25-03-2008, 26-03-2008 and 27-03-2008). Those provide reliable data for retrieval of optical depth, Rayleigh scattering, aerosol optical depth, column water vapor (through data in 936 nm) and with various parameters in 550nm, the horizontal visibility can be further developed by MODTRAN or 6S. The optical depth in 1640nm, 1020nm, 936nm, 870nm, 670nm, 550nm, 440nm, 380nm and 340nm were all acquired. Those data include the raw data in .k7 and can be opened by ASTPWin. ReadMe.txt is attached for detail. Processed data (after retrieval of the raw data) in Excel format are on optical depth, Rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. Accuracy of CE318 could be influenced by local air pressure, instrument calibration parameters, and convertion factors. (1) Most air pressure was derived from elevation-related empirical method, which was not reliable. For more accurate result, simultaneous data from the weather station are needed. (2) Errors in instrument calibration parameters need correcting. Thus field calibration based on Langly or interior instrument calibration in the standard light is required. (3) Convertion factors for retrieval of aerosol optical depth and the water vapor of the water vapor channel were also from the empirical method, and need further validation. Raw data were archived in .k7 format and can be opened by ASTPWin. ReadMe.txt is attached for detail. Preprocessed data (after retrieval of the raw data) in Excel format are on optical depth, Rayleigh scattering, aerosol optical depth, the horizontal visibility, the near surface air temperature, the solar azimuth, zenith, solar distance correlation factors, and air column mass number. Langley was used for the instrument calibration. Two subfolders including raw data and processed data (Geometric Positions and the Total Optical Depth of Each Channel and Rayleigh Scattering and Aerosol Optical Depth of Each Channel), and three data files (Directions on Data Observations, Raw Data and Proprocessed Data) were archived.
FANG Li, SU Gaoli, LIU Qinhuo
The dateset of TIR (Patent No.: ZL 02 2 37640.2) emissivity measurements was obtained in No. 3 quadrate of the A'rou foci experimental area on Mar. 14, 2008. The observation site was covered with dry pasture with height less than 5cm, in which the center point of each grid was measured twice and was named in the form of A3-9 (number 9 point in No. 3 quadrate of A'rou). Each measurement was carried out at 45° and followed strictly the order: Tsky, Tcha, Tsm and Tcm. Meanwhile, the surface temperature was also acquired by the handheld infrared thermometer and the thermal imager (FLIR ThermaCAM). [emissivity=1- (Tcm^4 – Tsm^4)/ (Tcha^4 – Tsky^4)]. Those provide reliable data for retrieval and study of the surface temperature, and energy and radiation balance.
CAO Yongpan, GU Juan, LI Hua
This dataset contains the flux measurements from the large aperture scintillometer (LAS) at site No.1 in the flux observation matrix. There were two types of LASs at site No.1: German BLS900 and China zzlas. The observation periods were from 7 June to 19 September, 2012, and 16 June to 19 September, 2012, for the BLS900 and the zzlas, respectively. The north tower is placed with the receiver of BLS900 and the transmitter of zzlas, and the south tower is placed with the transmitter of BLS900 and the receiver of zzlas. The site (north: 100.352° E, 38.884° N; south: 100.351° E, 38.855° N) was located in the Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1552.75 m. The underlying surface between the two towers contains corn, greenhouse, and village. The effective height of the LASs was 33.45 m; the path length was 3256 m. Data were sampled at 1 min intervals. Raw data acquired at 1 min intervals were processed and quality-controlled. The data were subsequently averaged over 30 min periods. The main quality control steps were as follows. (1) The data were rejected when Cn2 was beyond the saturated criterion (Cn2>3.05E-14). (2) Data were rejected when the demodulation signal was small (BLS900: Average X Intensity<1000; zzlas: Demod<-40 mv). (3) Data were rejected within 1 h of precipitation. (4) Data were rejected at night when weak turbulence occurred (u* was less than 0.1 m/s). The sensible heat flux was iteratively calculated by combining with meteorological data and based on Monin-Obukhov similarity theory. There were several instructions for the released data. (1) The data were primarily obtained from BLS900 measurements; missing flux measurements from the BLS900 were filled with measurements from the zzlas. Missing data were denoted by -6999. (2) The dataset contained the following variables: data/time (yyyy-mm-dd hh:mm:ss), the structural parameter of the air refractive index (Cn2, m-2/3), and the sensible heat flux (H_LAS, W/m^2). (3) 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. 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.
LIU Shaomin, LI Xin, XU Ziwei
This dataset contains the flux measurements from site No.14 eddy covariance system (EC) in the flux observation matrix from 30 May to 21 September, 2012. The site (100.35310° E, 38.85867° N) was located in a cropland (maize surface) in Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1570.23 m. The EC was installed at a height of 4.6 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&Li7500) 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.
LIU Shaomin, LI Xin, XU Ziwei
Based on the geostationary satellites and reanalysis data, the China Regional Atmospheric Driving Dataset is a set of atmospheric driving data sets with high spatiotemporal resolution prepared by the China Meteorological Administration, with a spatial resolution of 0.1 ° × 0.1 ° and a temporal resolution of 1 Hours, covering a range of 75 ° -135 ° east longitude and 15 ° -55 ° north latitude, include 6 elements of near-surface temperature, relative humidity, ground pressure, near-surface wind speed, incident solar radiation on the ground, and ground precipitation rate. The preparation process of precipitation products is as follows: The 6-hour cumulative precipitation estimated from the multi-channel data of the China Fengyun-2 geostationary satellite is integrated with the 6-hour cumulative precipitation from conventional ground observations to obtain 6-hour cumulative precipitation spatial distribution data, and then use the high-resolution cloud classification information retrieved from the multi-channel inversion of the geostationary satellites determines the interpolation time weight of the cumulative precipitation and obtains an estimated one-hour cumulative precipitation. The preparation process of the radiation data is as follows: The surface incident solar radiation based on FY-2C, uses the radiation transmission model DISORT (Discrete Ordinates Radiative Transfer Program for a Multi-Layered Plane-parallel Medium) to calculate the radiation transmission and obtains the data of surface incident solar radiation in China. Preparation process of other elements: The space and time interpolation method is used for the NCEP reanalysis data of 1.0 ° × 1.0 ° to obtain driving factors such as near-surface air temperature, relative humidity, ground pressure, and near-surface wind speed of 0.1 ° × 0.1 ° per hour. Physical meaning of each variable: Meteorological Elements || Variable Name || Unit || Physical Meaning | Surface temperature || TBOT || K || Surface temperature (2m) | Surface pressure || PSRF || Pa || Surface pressure | Relative humidity on the ground || RH || kg / kg || Relative humidity near the ground (2m) | Wind speed on the ground || WIND || m / s || Wind speed near the ground (anemometer height) | Surface incident solar radiation || FSDS || W / m2 || Surface incident solar radiation | Precipitation Rate || PRECTmms || mm / hr || Precipitation Rate For more information, see the data documentation published with the data.
SHI Chunxiang
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.
LIU Shaomin, LI Xin, XU Ziwei
Our project entrust the L band radiosonde sounding encrypt observations to Zhangye National Climate Observatory, and collect regular observation twice a day. The dataset contains three times one day at 8:00, 14:00, 20:00, which can support the remote sensing image atmospheric correction and atmospheric science research. Observation Site: Zhangye National Climate Observatory located in Shajing Town, west of ZhangYe. The coordinates of this site: 39°5′15.68" N, 100°16′39.11" E。 Observation Instrument: China Meteorological Administration Operational L Band radiosonde system. Observation Time: The observation date last from 1 May, 2012 to 31 September, 2012, among which: Three times observations at 7:00-8:00, 13:00-14:00 and 19:00-20:00 during 1 June, 2012 to 31 August, 2012; twice at 7:00-8:00 and 19:00-20:00 during 2012-5-1 to 5-31 and 2012-9-1 to 9-31. Accessory data： Pressure, temperature, relative humidity, wind speed and wind direction profiles data.
MA Mingguo
The dataset of the albedo measurements was obtained by the shortwave radiometer (KippZonen CMP3, 310nm-2800nm, 1m above the ground) in the Linze station foci experimental area. Sand, psammophyte and withered annual herbs in A9 of the south-north desert strip and LY07, and flax, maize and tomatoes in Linze station were measured on May 28, Jun. 5, 6, 15, 22, 25, 30 and Jul. 4, 2008. Voltage was measured manually by the digital multimeter (UNIT) at intervals of 2 minutes for albedo from May 28 to Jun. 22; self-recording Campbell CR1000 was used at intervals of 1s from Jun. 25 to Jul. 4. TIMESTAMP (observation time), SOLAR_UP_AVG (downward shortwave radiation), SOLAR_DOWN_AVG (upward shortwave radiation), SOLAR_NET_AVG (net radiation)= SOLAR_UP_AVG - SOLAR_DOWN_AVG, albedo_Avg (albedo) = SOLAR_DOWN_AVG / SOLAR_UP_AVG, batt_volt_Min (voltage), and ptemp (CR1000 temperature) were all recorded. Manual data were archived as Excel files and the self-recording data in .dat, which were processed into Excel.
BAI Yanfen, Qian Jinbo, ZHU Shijie, SONG Yi
The data set contains the observation data of the eddy covariance system of Sidaoqiao superstation which is located along the lower reaches of the Heihe Hydrometeorological observation network, and the data set covers data from January 1, 2017 to December 31, 2017. The station is located in Sidao Bridge, Ejina Banner, Inner Mongolia, and the underlying surface is Tamarix. The latitude and longitude of the observation station is 101.1374E, 42.0012N, and the altitude is 873 m. The height of the eddy covariance system is 8 meters, the sampling frequency is 10Hz, the ultrasonic orientation is positive north, and the distance between the ultrasonic wind speed and temperature monitor (CSAT3) and the CO2/H2O analyzer (Li7500) is 15cm. The original observation data of the eddy covariance system is 10 Hz, and the released data is a 30-minute data processed by Eddypro software. The main steps of the processing include: outlier eliminating, delay time correction, coordinates rotation (secondary coordinates rotation), frequency response correction, ultrasonic virtual temperature correction and density (WPL) correction, etc. Meanwhile, the quality evaluation of each flux value was performed,mainly includes atmospheric stability (Δst) test and turbulence similarity (ITC) test. The 30-min flux value output of Eddypro software was also screened: (1) Data from the instrument error was eliminated; (2) Data obtained with one hour before and after precipitation was removed; (3) Data with a deletion rate greater than 10% of the 10 Hz raw data every 30 minutes was eliminated; (4) Observation data of weak turbulence at night (u* less than 0.1 m/s) was excluded. The average period of observation data is 30 minutes, 48 data per day, and the missing data is marked as -6999. The data was missing due to Li7500 calibration of the eddy system on April 7 and 8; the suspicious data caused by instrument drift and other reasons was marked by red fonts. Published observation data include: date/time Date/Time, wind direction(°), horizontal wind speed(m/s), lateral wind speed standard deviation(m/s), ultrasonic virtual temperature (°C), water vapor density (g/m3), carbon dioxide concentration(mg/m3), friction velocity (m/s), length (m), sensible heat flux(W/m2), latent heat flux (W/m2), carbon dioxide flux (mg/(m2s)), sensible heat flux quality identification QA_Hs, latent heat flux quality identification QA_LE, carbon dioxide flux quality identification QA_Fc. The quality identification of sensible heat, latent heat, and carbon dioxide flux is divided into three levels (quality mark 0: (Δst <30, ITC<30); 1: (Δst <100, ITC<100); the rest is 2). The meaning of the data time, such as 0:30 represents an average data of 0:00-0:30; the data is stored in *.xls format. For hydrometeorological network or station information, please refer to Li et al. (2013). For observation data processing, please refer to Liu et al. (2011).
LIU Shaomin, LI Xin, CHE Tao, XU Ziwei, REN Zhiguo, TAN Junlei
This dataset contains the flux measurements from site No.11 eddy covariance system (EC) in the flux observation matrix from May 29 to September 18, 2012. The site (100.34197° E, 38.86991° N) was located in a cropland (maize surface) in Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1575.65 m. The EC was installed at a height of 3.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&Li7500) 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.
LIU Shaomin, LI Xin, XU Ziwei
This dataset contains the flux measurements from the Daman superstation eddy covariance system (EC) at the highest layer in the flux observation matrix from 30 May to 15 September, 2012. The site (100.37223° E, 38.85551° N) was located in a cropland (maize surface) in Daman irrigation district, which is near Zhangye, Gansu Province. The elevation is 1556.06 m. The EC was installed at a height of 34 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/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.
LIU Shaomin, LI Xin, XU Ziwei
This dataset contains the flux measurements from site No.12 eddy covariance system (EC) in the flux observation matrix from 28 May to 21 September, 2012. The site (100.36631° E, 38.86515° N) was located in a cropland (maize surface) in Daman irrigation district, which is near Zhangye, Gansu Province. The elevation is 1559.25 m. The EC was installed at a height of 3.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&Li7500) 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.
LIU Shaomin, LI Xin, XU Ziwei
This dataset contains the flux measurements from site No.1 eddy covariance system (EC) in the flux observation matrix from 4 June to 17 September 2012. The site (100.35813° E, 38.89322° N) was located in a cropland (vegetable surface) in the Yingke irrigation district, which is near Zhangye, Gansu Province. The elevation is 1552.75 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 (Gill&Li7500A) 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.
LIU Shaomin, LI Xin, XU Ziwei
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS
0931-4967287
poles@itpcas.ac.cn
Links
National Tibetan Plateau Data CenterFollow Us
A Big Earth Data Platform for Three Poles © 2018-2020 No.05000491 | All Rights Reserved
｜ 
No.11010502040845
Tech Support: westdc.cn
