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
The data is based on the Harmonized World Soil Database version 1.1 (HWSD) constructed by the Food and Agriculture Organization of the United Nations (FAO) and the Vienna International Institute for Applied Systems (IIASA). The data source of China is 1: 1 million soil data in the second national land survey provided by the Nanjing Soil Research Institute. The data can provide model input parameters for modelers, in agricultural perspective, it can be used to study eco-agricultural zoning, food security and climate change. The data format is grid and the projection is WGS84. The soil classification system used is mainly FAO-90. The main fields of the soil property table include: SU_SYM90 (the soil name in the FAO90 soil classification system); SU_SYM85 (FAO85 classification); T_TEXTURE (top soil texture); DRAINAGE (19.5); REF_DEPTH (soil reference depth); AWC_CLASS (19.5); AWC_CLASS (soil effective water content); PHASE1: Real (soil phase); PHASE2: String (soil phase); ROOTS: String (depth classification with obstacles to the bottom of the soil); SWR: String (soil moisture characteristics); ADD_PROP: Real (a specific soil type related to agricultural use in the soil unit); T_GRAVEL: Real (gravel volume percentage); T_SAND: Real (sand content); T_SILT: Real (silt content); T_CLAY: Real (clay content); T_USDA_TEX: Real (USDA soil texture classification); T_REF_BULK: Real (soil bulk density); T_OC: Real (organic carbon content); T_PH_H2O: Real (pH) T_CEC_CLAY: Real (cation exchange capacity of clay soil); T_CEC_SOIL: Real (cation exchange capacity of soil) T_BS: Real (basic saturation); T_TEB: Real (exchangeable base); T_CACO3: Real (carbonate or lime content) T_CASO4: Real (sulfate content); T_ESP: Real (exchangeable sodium salt); T_ECE: Real (conductivity). The attribute field beginning with T_ indicates the upper soil attribute (0-30cm), and the attribute field beginning with S_ indicates the lower soil attribute (30-100cm). For the meaning of specific attribute values, please refer to the documentation * .pdf and database * .mdb in the folder.
Food and Agriculture Organization of the United Nations（FAO） International Institute for Applied Systems Analysis
This data set is an upgraded version of the “Long-term series of daily snow depth dataset in China". This dataset provides daily data of snow depth distribution in China from January 1, 1979, to December 31, 2019, with a spatial resolution of 0.25 degrees. The original data used to derive the snow depth dataset are the daily passive microwave brightness temperature data (EASE-Grid) from SMMR (1979-1987), SSM/I (1987-2007) and SSMI/S (2008-2019) which were archived in the National Snow and Ice Data Center (NSIDC). Because the brightness temperatures come from different sensors, there is a certain system inconsistency among them. Therefore, before the derivation of snow depth, the inter-sensor calibration were performed to improve the temporal consistency of the brightness temperature data. Based on the calibrated brightness temperatures, the modified Chang algorithm developed by Dr. Tao Che, was used to retrieve daily snow depth. The algorithm details were introduced in the data specification document- “Long-term Sequence Data Set of China Snow Depth (1979-2019) Introduction. doc". The projection of the data set is latitude and longitude. The data of each day was stored in a file, and the naming convention of which is year + day; for example, 1990001 represents the first day of 1990, and 1990207 represents the 207th day of 1990. For a detailed data description, please refer to the data specification document.
CHE Tao DAI Liyun
This dataset includes the monthly precipitation data with 0.0083333 arc degree (~1km) for China from Jan 1901 to Dec 2017. The data form belongs to NETCDF, namely .nc file. The unit of the data is 0.1 mm. The dataset was spatially downscaled from CRU TS v4.02 with WorldClim datasets based on Delta downscaling method. The dataset was evaluated by 496 national weather stations across China, and the evaluation indicated that the downscaled dataset is reliable for the investigations related to climate change across China. The dataset covers the main land area of China, including Hong Kong, Macao and Taiwan regions, and excluding islands and reefs in South China Sea.
The dataset is a 34-year (1983.7-2017.6) 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.
China's second glacier inventory uses the high-resolution Landsat TM/ETM+ remote sensing satellite data as the main glacier boundary data source and extracts the data source with the latest global digital elevation model, SRTM V4, as the glacier attribute, using the current international ratio threshold segmentation method to extract the glacier boundary in bare ice areas. The ice ridge extraction algorithm is developed to extract the glacier ice ridge, and it is used for the segmentation of a single glacier. At the same time, the international general algorithm is used to calculate the glacier attributes, so that the vector data and attribute data that contain the glacier information of the main glacier regions in west China are obtained. Compared with some field GPS field measurement data and higher resolution remote sensing images (such as from QuickBird and WorldView), the glacial vector data in the second glacier inventory data set of China have higher positioning accuracy and can meet the requirements for glacial data in national land, water conservancy, transportation, environment and other fields. Glacier inventory attributes: Glc_Name, Drng_Code, FCGI_ID, GLIMS_ID, Mtn_Name, Pref_Name, Glc_Long, Glc_Lati, Glc_Area, Abs_Accu, Rel_Accu, Deb_Area, Deb_A_Accu, Deb_R_Accu, Glc_Vol_A, Glc_Vol_B, Max_Elev, Min_Elev, Mean_Elev, MA_Elev, Mean_Slp, Mean_Asp, Prm_Image, Aux_Image, Rep_Date, Elev_Src, Elev_Date, Compiler, Verifier. For a detailed data description, please refer to the second glacier inventory data description.
LIU Shiyin GUO Wanqin XU Junli
This dataset (version 1.5) is derived from the complementary-relationship method, with inputs of CMFD downward short- and long-wave radiation, air temperature, air pressure, GLASS albedo and broadband longwave emissivity, ERA5-land land surface temperature and humidity, and NCEP diffuse skylight ratio, etc. This dataset covers the period of 1982-2017, and the spatial coverage is Chinese land area. This dataset would be helpful for long-term hydrological cycle and climate change research. Land surface actual evapotranspiration (Ea)，unit: mm month-1. The spatial resolution is 0.1-degree; The temporal resolution is monthly; The data type is NetCDF; This evapotranspiration dataset is only for land surface.
MA Ning Jozsef Szilagyi ZHANG Yinsheng LIU Wenbin
The multi-decadal lake number and area changes in China during 1960s–2015 are derived from historical topographic maps and >3831 Landsat satellite images, including lakes as fine as ≥1 km2 in size. The total area of lakes in China has increased by 5858.06 km2 (9%) between 1960s and 2015, and with heterogeneous spatial variations. Lake area changes in the Tibetan Plateau, Xinjiang, and Northeast Plain and Mountain regions reveal significant increases of 5676.75, 1417.15, 1134.87 km2 (≥15%), respectively, but the Inner-Mongolian Plateau shows an obvious decrease of 1223.76 km2 (22%). We find that 141 new lakes have appeared predominantly in the arid western China; but 333 lakes, mainly located in the humid eastern China, have disappeared over the past five decades.
The dataset includes soil physical and chemical attributes: pH value, organic matter fraction, cation exchange capacity, root abundance, total nitrogen (N), total phosphorus (P), total potassium (K), alkali-hydrolysable N, available P, available K, exchangeable H+, Al3+, Ca2+, Mg2+, K+ , Na+, horizon thickness, soil profile depth, sand, silt and clay fractions, rock fragment, bulk density, porosity, structure, consistency and soil color. Quality control information (QC) was provided. The resolution is 30 arc-seconds (about 1 km at the equator). The vertical variation of soil property was captured by eight layers to the depth of 2.3 m (i.e. 0- 0.045, 0.045- 0.091, 0.091- 0.166, 0.166- 0.289, 0.289- 0.493, 0.493- 0.829, 0.829- 1.383 and 1.383- 2.296 m) for convenience of use in the Common Land Model and the Community Land Model (CLM). 1.THSCH.nc: Saturated water content of FCH 2.PSI_S.nc: Saturated capillary potential of FCH 3.LAMBDA.nc: Pore size distribution index of FCH 4.K_SCH.nc: Saturate hydraulic conductivity of FCH 5.THR.nc: Residual moisture content of FGM 6.THSGM.nc: Saturated water content of FGM 7.ALPHA.nc: The inverse of the air-entry value of FGM 8.N.nc: The shape parameter of FGM 9.L.nc: The pore-connectivity parameter of FGM 10.K_SVG.nc: Saturated hydraulic conductivity of FGM 11.TH33.nc: Water content at -33 kPa of suction pressure, or field capacity 12.TH1500.nc: Water content at -1500 kPa of suction pressure, or permanent wilting point
DAI Yongjiu SHANGGUAN Wei
Soil data is important both on a global scale and on a local scale, and due to the lack of reliable soil data, land degradation assessments, environmental impact studies, and sustainable land management interventions have received significant bottlenecks . Affected by the urgent need for soil information data around the world, especially in the context of the Climate Change Convention, the International Institute for Applied Systems Analysis (IIASA) and the Food and Agriculture Organization of the United Nations (FAO) and the Kyoto Protocol for Soil Carbon Measurement and FAO/International The Global Agroecological Assessment Study (GAEZ v3.0) jointly established the Harmonized World Soil Database version 1.2 (HWSD V1.2). Among them, the data source in China is the second national land in 1995. Investigate 1:1,000,000 soil data provided by Nanjing Soil. The resolution is 30 seconds (about 0.083 degrees, 1km). The soil classification system used is mainly FAO-90. The core soil system unit unique verification identifier: MU_GLOBAL-HWSD database soil mapping unit identifier, connected to the GIS layer. MU_SOURCE1 and MU_SOURCE2 source database drawing unit identifiers SEQ-soil unit sequence in the composition of the soil mapping unit; The soil classification system utilizes the FAO-7 classification system or the FAO-90 classification system (SU_SYM74 resp. SU_SYM90) or FAO-85 (SU_SYM85). The main fields of the soil property sheet include: ID (database ID) MU_GLOBAL (Soil Unit Identifier) (Global) SU_SYMBOL soil drawing unit SU_SYM74 (FAO74 classification); SU_SYM85 (FAO85 classification); SU_SYM90 (name of soil in the FAO90 soil classification system); SU_CODE soil charting unit code SU_CODE74 soil unit name SU_CODE85 soil unit name SU_CODE90 soil unit name DRAINAGE (19.5); REF_DEPTH (soil reference depth); AWC_CLASS(19.5); AWC_CLASS (effective soil water content); PHASE1: Real (soil phase); PHASE2: String (soil phase); ROOTS: String (depth classification to the bottom of the soil); SWR: String (soil moisture content); ADD_PROP: Real (specific soil type in the soil unit related to agricultural use); T_TEXTURE (top soil texture); T_GRAVEL: Real (top gravel volume percentage); (unit: %vol.) T_SAND: Real (top sand content); (unit: % wt.) T_SILT: Real (surface layer sand content); (unit: % wt.) T_CLAY: Real (top clay content); (unit: % wt.) T_USDA_TEX: Real (top layer USDA soil texture classification); (unit: name) T_REF_BULK: Real (top soil bulk density); (unit: kg/dm3.) T_OC: Real (top organic carbon content); (unit: % weight) T_PH_H2O: Real (top pH) (unit: -log(H+)) T_CEC_CLAY: Real (cation exchange capacity of the top adhesive layer soil); (unit: cmol/kg) T_CEC_SOIL: Real (cation exchange capacity of top soil) (unit: cmol/kg) T_BS: Real (top level basic saturation); (unit: %) T_TEB: Real (top exchangeable base); (unit: cmol/kg) T_CACO3: Real (top carbonate or lime content) (unit: % weight) T_CASO4: Real (top sulfate content); (unit: % weight) T_ESP: Real (top exchangeable sodium salt); (unit: %) T_ECE: Real (top conductivity). (Unit: dS/m) S_GRAVEL: Real (bottom crushed stone volume percentage); (unit: %vol.) S_SAND: Real (bottom sand content); (unit: % wt.) S_SILT: Real (bottom sludge content); (unit: % wt.) S_CLAY: Real (bottom clay content); (unit: % wt.) S_USDA_TEX: Real (bottom USDA soil texture classification); (unit: name) S_REF_BULK: Real (bottom soil bulk density); (unit: kg/dm3.) S_OC: Real (underlying organic carbon content); (unit: % weight) S_PH_H2O: Real (bottom pH) (unit: -log(H+)) S_CEC_CLAY: Real (cation exchange capacity of the underlying adhesive layer soil); (unit: cmol/kg) S_CEC_SOIL: Real (cation exchange capacity of the bottom soil) (unit: cmol/kg) S_BS: Real (underlying basic saturation); (unit: %) S_TEB: Real (underlying exchangeable base); (unit: cmol/kg) S_CACO3: Real (bottom carbonate or lime content) (unit: % weight) S_CASO4: Real (bottom sulfate content); (unit: % weight) S_ESP: Real (underlying exchangeable sodium salt); (unit: %) S_ECE: Real (underlying conductivity). (Unit: dS/m) The database is divided into two layers, with the top layer (T) soil thickness (0-30 cm) and the bottom layer (S) soil thickness (30-100 cm). For other attribute values, please refer to the HWSD1.2_documentation documentation.pdf, The Harmonized World Soil Database (HWSD V1.2) Viewer-Chinese description and HWSD.mdb.
Meng Xianyong Wang Hao
Gridded climatic datasets with fine spatial resolution can potentially be used to depict the climatic characteristics across the complex topography of China. In this study we collected records of monthly temperature at 1153 stations and precipitation at 1202 stations in China and neighboring countries to construct a monthly climate dataset in China with a 0.025° resolution (~2.5 km). The dataset, named LZU0025, was designed by Lanzhou University and used a partial thin plate smoothing method embedded in the ANUSPLIN software. The accuracy of LZU0025 was evaluated based on three aspects: (1) Diagnostic statistics from the surface fitting model during 1951–2011. The results indicate a low mean square root of generalized cross validation (RTGCV) for the monthly air temperature surface (1.06 °C) and monthly precipitation surface (1.97 mm1/2). (2) Error statistics of comparisons between interpolated monthly LZU0025 with the withholding of climatic data from 265 stations during 1951–2011. The results show that the predicted values closely tracked the real true values with values of mean absolute error (MAE) of 0.59 °C and 70.5 mm, and standard deviation of the mean error (STD) of 1.27 °C and 122.6 mm. In addition, the monthly STDs exhibited a consistent pattern of variation with RTGCV. (3) Comparison with other datasets. This was done in two ways. The first was via comparison of standard deviation, mean and time trend derived from all datasets to a reference dataset released by the China Meteorological Administration (CMA), using Taylor diagrams. The second was to compare LZU0025 with the station dataset in the Tibetan Plateau. Taylor diagrams show that the standard deviation, mean and time trend derived from LZU had a higher correlation with that produced by the CMA, and the centered normalized root-mean-square difference for this index derived from LZU and CMA was lower. LZU0025 had high correlation with the Coordinated Energy and Water Cycle Observation Project (CEOP) - Asian Monsoon Project, (CAMP) Tibet surface meteorology station dataset for air temperature, despite a non-significant correlation for precipitation at a few stations. Based on this comprehensive analysis, we conclude that LZU0025 is a reliable dataset. LZU0025, which has a fine resolution, can be used to identify a greater number of climate types, such as tundra and subpolar continental, along the Himalayan Mountain. We anticipate that LZU0025 can be used for the monitoring of regional climate change and precision agriculture modulation under global climate change.
HUANG Wei ZHAO Hong
This dataset includes the monthly minimum temperature data with 0.0083333 arc degree (~1km) for China from Jan 1901 to Dec 2017. The data form belongs to NETCDF, namely .nc file. The unit of the data is 0.1 ℃. The dataset was spatially downscaled from CRU TS v4.02 with WorldClim datasets based on Delta downscaling method. The dataset was evaluated by 496 national weather stations across China, and the evaluation indicated that the downscaled dataset is reliable for the investigations related to climate change across China. The dataset covers the main land area of China, including Hong Kong, Macao and Taiwan regions, and excluding islands and reefs in South China Sea.
The data set was produced based on the SRTM DEM data collected by Space Shuttle Radar terrain mission in 2016, the reference data such as river, lake and other water system auxiliary data , using the arcgis hydrological model to analyze and extract the river network. There are 12 sub-basins over the Tibet Plateau, including AmuDayra、Brahmaputra、Ganges、Hexi、Indus、Inner、Mekong、Qaidam、Salween、Tarim、Yangtze、Yellow. The outer boundary is based on the 2500-metre contour line and national boundaries.
PML_V2 terrestrial evapotranspiration and total primary productivity dataset, including gross primary product (GPP), vegetation transpiration (Ec), soil evaporation (Es), vaporization of intercepted rainfall , Ei) and water body, ice and snow evaporation (ET_water), a total of 5 elements. The data format is tiff, the space-time resolution is 8 days, 0.05°, and the time span is 2002.07-2019.08. Based on the Penman-Monteith-Leuning (PML) model, PML_V2 is coupled to the GPP process based on stomatal conductance theory. GPP and ET mutually restrict and restrict each other, which makes PML_V2 in ET simulation accuracy, which is greatly improved compared with the previous model. The parameters of PML_V2 are divided into different vegetation types and are determined on 95 vorticity-related flux stations around the world. The parameters were then migrated globally according to the MODIS MCD12Q2.006 IGBP classification. PML_V2 uses GLDAS 2.1 meteorological drive and MODIS leaf area index (LAI), reflectivity (Albedo), emissivity (Emissivity) as inputs, and finally obtains PML_V2 terrestrial evapotranspiration and total primary productivity data sets.
This dataset includes the monthly maximum temperature data with 0.0083333 arc degree (~1km) for China from Jan 1901 to Dec 2017. The data form belongs to NETCDF, namely .nc file. The unit of the data is 0.1 ℃. The dataset was spatially downscaled from CRU TS v4.02 with WorldClim datasets based on Delta downscaling method. The dataset was evaluated by 496 national weather stations across China, and the evaluation indicated that the downscaled dataset is reliable for the investigations related to climate change across China. The dataset covers the main land area of China, including Hong Kong, Macao and Taiwan regions, and excluding islands and reefs in South China Sea.
The integration dataset of Tibetan Plateau boundary includes: TPBoundary_2500m：Based on ETOPO5 Global Surface Relief, ENVI+IDL is used to extract the longitude of the Tibetan Plateau (65~105) and the altitude of 2500 meters above the latitude (20~45); TPBoundary_3000m：Based on ETOPO5 Global Surface Relief, ENVI+IDL is used to extract the longitude of the Tibetan Plateau (65~105) and the altitude of 3000 meters above the latitude (20~45); TPBoundary_HF (High Frequency):Li Bingyuan (1987) has conducted a systematic discussion on the principle and specific boundary of determining the scope of the Qinghai-Tibet Plateau. From the perspective of the formation and basic characteristics of plateau geomorphology, Based on the geomorphological features, the plateau surface and its altitude, and considering the integrity of the mountain as the basic principle for determining the plateau range.Zhang Yili (2002) according to the results of new research in related fields and years of field practice, demonstration principles to determine the scope and boundaries of the Tibetan Plateau, Based on the information technology method, the location and boundary position of the Qinghai-Tibet Plateau are accurately located and quantitatively analyzed. It is concluded that the Qinghai-Tibet Plateau is partly in the Pamir Plateau in the west, the Hengduan Mountains in the east, the southern margin of the Himalayas in the south, and the Kunlun Mountains in the north. Mountain - north side of Qilian Mountain. On April 14, 2017, the Ministry of Civil Affairs of the People's Republic of China issued the "Announcement on Supplementing the Public Use of Place Names in the Southern Region of Tibet (First Batch)", adding Wujianling, Mirage, Qu Dengbu, and Mechuca 6 places in southern Tibet such as Baimingla Mountain Pass and Namkam;. TPBoundary_rectangle：According to the range Lon (63~105E) & Lat (20~45N), The data is projected using latitude and longitude WGS84.. Project source: national natural science foundation of China (41571068,41301063) Spatial range and projection mode of data: elevation greater than 2500m, WGS84 projection As the basic data, the boundary of qinghai-tibet plateau can be used as a reference for all kinds of geoscientific research on Qinghai-Tibet Plateau.
ZHANG Yili REN Huixia PAN Xiaoduo
This dataset includes the monthly mean temperature data with 0.0083333 arc degree (~1km) for China from Jan 1901 to Dec 2017. The data form belongs to NETCDF, namely .nc file. The unit of the data is 0.1 ℃. The dataset was spatially downscaled from CRU TS v4.02 with WorldClim datasets based on Delta downscaling method. The dataset was evaluated by 496 national weather stations across China, and the evaluation indicated that the downscaled dataset is reliable for the investigations related to climate change across China. The dataset covers the main land area of China, including Hong Kong, Macao and Taiwan regions, and excluding islands and reefs in South China Sea.
This data is originated from the 1:100,000 national basic geographic database, which was open freely for public by the National Basic Geographic Information Center in November 2017. The boundary of the Qinghai-Tibet Plateau was spliced and clipped as a whole, so as to facilitate the study on the Qinghai-Tibet plateau. This data set is the 1:100,000 administrative boundaries of the qinghai-tibet plateau, including National_Tibet_line、 Province_Tibet、City_Tibet、County_Tibet_poly and County_Tibet_line. Administrative boundary layer (County_Tibet_poly) property name and definition: Item Properties Describe Example PAC Administrative division code 513230 NAME The name of the County line name Administrative boundary layer (BOUL) attribute name and definition: Item Properties Describe Example GB classification code 630200 Administrative boundary layer (County_Tibet_line) attribute item meaning: Item Properties Describe Example GB 630200 Provincial boundary GB 640200 Prefectural, municipal and state administrative boundaries GB 650201 county administrative boundaries (determined)
National Basic Geographic Information Center
Based on 2015 ESA global land cover data (ESA GlobCover), combined with the Tsinghua university global land cover data (FROM GLC)、NASA MODIS global land cover data (MCD12Q1)、University of Maryland global land data (UMD)、USGS global land data (IGBP DISCover),we build the LUC classification system in the Tibet Plateau and the rest of the data transformation rules of the classification system. We also build the land cover classification confidence function and the rules of fusing land classification to finish the Integration and modification of land cover products and finally complet the land use data in the Tibet Plateau V1.0.
A multi-layer soil particle-size distribution dataset (sand, silt and clay content), based on USDA (United States Department of Agriculture) standard for regional land and climate modelling in China. was developed The 1:1,000,000 scale soil map of China and 8595 soil profiles from the Second National Soil Survey served as the starting point for this work. We reclassified the inconsistent soil profiles into the proper soil type of the map as much as possible because the soil classification names of the map units and profiles were not quite the same. The sand, silt and clay maps were derived using the polygon linkage method, which linked soil profiles and map polygons considering the distance between them, the sample sizes of the profiles, and soil classification information. For comparison, a soil type linkage was also generated by linking the map units and soil profiles with the same soil type. The quality of the derived soil fractions was reliable. Overall, the map polygon linkage offered better results than the soil type linkage or the Harmonized World Soil Database. The dataset, with a 1-km resolution, can be applied to land and climate modelling at a regional scale. Data characteristics： projection：projection Coverage: China Resolution: 0.00833 (about 1 km) Data format: FLT, TIFF Value range: 0%-100% Document describing： Floating point raster files include: Sand1. FLT, clay1. FLT -- surface (0-30cm) sand, clay content. Sand2. FLT, clay2. FLT -- content of sand and clay in the bottom layer (30-100cm). PSD. HDR -- header file: Ncols - the number of columns Nrows- rows Xllcorner - latitude in the lower left corner Yllcorner - longitude of the lower left corner Cellsize - cellsize NODATA_value - a null value byteorder - LSBFIRST, Least Significant Bit First TIFF raster files include: Sand1. Tif, clay1. Tif - surface (0-30cm) sand, clay content. Sand2. Tif, clay2. Tif - bottom layer (30-100cm) sand, clay content.
SHANGGUAN Wei DAI Yongjiu