In 2018, the Second Tibetan Plateau Scientific Expedition and Research Programme tasked a research team to conduct an "Investigation of the water vapor channel of the Yarlung Zangbo Grand Canyon" in the southeastern Tibetan Plateau. This team subsequently established a three-dimensional comprehensive observation system of land-air interaction, water vapor transport, cloud cover, and rainfall activity along different altitude from south to north along the YGC canyon. Motuo meteorological station is an integrated observation base of a variety of large observation equipment. The whole comprehensive observation network includes 2 cloud radars, 2 micro rain radars, seven sets of eddy covariance flux sites for measuring land-air interaction, 3 sets of microwave radiometers , 6 sets of GPS water vapor observatories, 2 sets of automatic weather stations and 19 sets of rainfall buckets.

The data set contains nearly 15 years of eddy covariance data from an alpine steppe ecosystem on the central Tibetan Plateau.The data was processed following standardized quality control methods to allow for comparability between the different years of our record and with other data sets. To ensure meaningful estimates of ecosystem atmosphere exchange, careful application of the following correction procedures and analyses was necessary: (1) Due to the remote location, continuous maintenance of the eddy covariance (EC) system was not always possible, so that cleaning and calibration of the sensors was performed irregularly. Furthermore, the high proportion of bare soil and high wind speeds led to accumulation of dirt in the measurement path of the infrared gas analyzer (IRGA). The installation of the sensor in such a challenging environment resulted in a considerable drift in CO2 and H2O gas density measurements. If not accounted for, this concentration bias may distort the estimation of the carbon uptake. We applied a modified drift correction procedure following Fratini et al. (2014) which, instead of a linear interpolation between calibration dates, uses the CO2 concentration measurements from the Mt. Waliguan atmospheric observatory as reference time series. (2) We applied rigorous quality filtering of the calculated fluxes to retain only fluxes which represent actual physical processes. (3) During the long measurement period, there were several buildings constructed in the near vicinity of the EC system. We investigated the influence of these obstacles on the turbulent flow regime to identify fluxes with uncertain land cover contribution and exclude them from subsequent computations. (4) We calculated the de-facto standard correction for instrument surface heating during cold conditions (hereafter called sensor self heating correction) following Burba et al. (2008) and a revision of the original method following Frank and Massman (2020). (5)Subsequently, we applied the traditional and widely used gap filling procedure following Reichstein et al. (2005) to provide a more complete overview of the annual net ecosystem CO2 exchange.(6) We estimated the flux uncertainty by calculating the random flux error (RE) following Finkelstein and Sims (2001) and by using the standard deviation of the fluxes used for gap filling(NEE_fsd) as a measure for spatial and temporal variation.

Data set contains tree age of trees growing at different glacier moraines in the central Himalayas. The data were obtained using tree ring samples. Cores samples were collected (almost near to the ground level to estimate the minimum age of the related moraine) using an increment borer. Samples were processed by using standard dendrochronological techniques.

We use waveform cross-correlation to analyze the recordings of eight earthquakes (2009-2018) beneath the Indian Ocean at stations from the Chinese Digital Seismic Network. We obtain 929 high quality residual traveltime differences between the phases ScS and S (Differential traveltimes.dat). We interpret variations of δt up to 10 seconds as due to horizontal shear-velocity variations in D” beneath northern India, Nepal, and southwestern China. The shear velocity can vary by as much as 7% over distances shorter than 300 km. Our observations provide additional observational evidence that compositional heterogeneity and possibly melt contribute to the seismic structure of the lower mantle characterized by long-term subduction and mantle downwelling.

The land cover dataset of Pan third pole major cities contains 14 cities (Urumqi, Xining, Lanzhou, Dhaka, Kathmandu, Lucknow, Delhi, Lahore, Islamabad, Kabul, Dushanbe, Tashkent, Bishkek and Almaty) in 2000 / 2010 / 2017, the spatial resolution of this dataset is 30 m. It includes vegetation, cultivated land, artificial surface, water body and others. Based on globeland30, mcd12q1 and globcover2009, the consistent regions were identified and retained. The inconsistent regions were reclassified by deep learning method, and the final classification results were obtained by fusing the above regions. The data has been verified by visual interpretation. The data are applied to the study of construction land dynamics and anthropogenic influence in Pan-Third Pole cities. Data type: grid. Projection mode: UTM projection.

This data comes from the result of teleseismic data, mainly including the velocity and radial anisotropic structures beneath western Tibet. In the process of processing, bandwidth filtering is adopted, and the filtering range is 0.05-2 Hz. Due to the use of teleseismic data, the cross-correlation method is used in the acquisition process to "align" the waveform. The data quality is good, because the extracted data are all from the earthquakes with magnitude greater than 5.0 located in the global seismic catalog, and each event has an obvious take-off point. The data can be used by other seismologists to reconstruct and analyze the underground structures in this area.

Data content: Standard ring-width chronology derived from Wilson juniper shrub around the northern shore of the Nam Co Lake; May-June SZI (Standardized Moisture Anomaly Index) drought reconstruction for the Nam Co region. Time span: 1605 to 2010. Temporal resolution: Yearly. Application and prospects: Hydroclimate study on the south-central Tibetan Plateau.

The sand drift potential (DP, in vector units (VU)) is calculated by DPi=∑U^2 [U-Ut]*fu where i represents 16 directions: N, NNE, NE, NEE, E, EES, ES, ESS, S, SSW, WS, WWS, W, WWN, NW and NNW; U is the effective sand-moving wind speed at the standard height of 10 m; Ut is the threshold wind velocity, i.e., the minimum wind velocity at the standard height to cause sand particle rolling; and fu is the fraction of time when the wind speed is higher than Ut. The 2 m s-1 bin is adopted in the effective sand-moving wind (sand-moving wind >6 m s-1 at the height of 10 m) directions, corresponding to the mean wind speeds of 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33 and 34 m s-1, to sum all the above results to obtain the final DP in the wind direction. The divisor used in calculating the frequency of effective sand-moving winds from different directions is the total hour number of Julian years (8760 hours for common years or 8784 hours for leap years). The wind speed and wind direction data from 2000 to 2008 were hourly estimates of 10 m u-component of wind and 10 m v-component wind with a horizontal resolution of 0.25°×0.25° generated with the fifth generation of ECMWF atmospheric ReAnalysis of the global climate (ERA5).

Relationship between modern pollen and climate, and its representative to vegetation are the important references in explaining and reconstructing past climate and vegetation qualitatively or quantitatively. To extrct past climate and vegetation signals from fossil pollen spectrum of a lacustrine sediment, a corresponding modern pollen dataset collected from lake-sediment surface is necessary. At present, there are a few modern pollen datasets extracted from lake sediment-surface established on the Tibetan Plateau, however, the geographic gaps (e.g. the central and east Tibetan Plateau) of available sampled lakes influence the correct understanding. To ensure the even distribution of the representative lakes, we collected lake sediment-surface samples (n=117) covering the alpine meadow evenly on the east and central Tibetan Plateau, in July and August 2018. For pollen extraction, approximately 10 g (wet original sediment) per sample were sub-sampled. Pollen sample was processed by the standard acid-alkali-acid procedures followed by 7-μm-mesh sieving. More than 500 terrestrial pollen grains were counted for each sample. Pollen assemblages of the dataset from alpine meadow are dominated by Cyperaceae (mean is 68.4%, maximum is 95.9%), with other herbaceous pollen taxa as commen taxa including Poaceae (mean is 10.3%, maximum is 87.7%), Ranunculaceae (mean is 4.8%, maximum is 33.6%), Artemisia (mean is 3.7%, maximum is 24.5%), Asteraceae (mean is 2.1%, maximum is 33.6%), etc. Salix (mean is 0.4%, maximum is 5.3%) is the major shrub taxon in these pollen assemblages, while arboreal taxa occur with low percentages generally (mean of total arboreal percentages is 0.9% (maximum is 5.8%), including mainly Pinus (mean is 0.3%, maximum is 1.8%), Betula (mean is 0.1%, maximum is 0.9%) and Alnus (mean is 0.1%, maximum is 0.7%). These pollen assemblages represent the plant components well in the alpine meadow communities, although they are influenced slightly by long-distance pollen grain transported by wind or river (such as these arboreal pollen taxa). Together with pollen counts and percentages, we also provided the modern climatic data for the sampled lakes. The China Meteorological Forcing Dataset (CMFD; gridded near-surface meteorological dataset) with a temporal resolution of three hours and a spatial resolution of 0.1° was employed, and the climatic data of the nearest pixel of one sampled lake was defined to represent climatic conditions of the lake. Finally, the mean annual precipitation (Pann), mean annual temperature (Tann) and mean temperature of the coldest month (Mtco) and warmest month (Mtwa) are calculated for each sampled lake.

The data set includes the reconstructed long-term annual, ablation-season, and cold-season glacier-wide mass balance and its components for Guliya Ice Cap and the reconstructed daily meteorological data on the glacier from 1969 to 2019. The reconstructed meteorological data includes air temperature (℃), relative humidity (%), wind speed (m s-1), air pressure (hPa) and downward shortwave radiation (W m-2) at an elevation of 6004 m a.s.l. and precipitation (mm) at an elevation of 5491 m a.s.l. ERA5 data from the grid point (35°N, 78°75′E) around AWS2 were calibrated by the measured meteorological data. The exact method has been described in the reference. The long-term mass balance of Guliya Ice Cap during 1970-2019 was reconstructed using an energy and mass balance model and calibrated ERA5 data, which was calibrated and validated by in-situ measurements and geodetic mass balances. Please see the reference. The data is stored in an Excel file. It can be used by researchers for studying the changes in climate, hydrology, glaciers, etc.

The data set includes annual mass balance of Naimona’nyi glacier (northern branch) from 2008 to 2018, daily meteorological data at two automatic meteorological stations (AWSs) near the glacier from 2011 to 2018 and monthly air temperature and relative humidity on the glacier from 2018 to 2019. In the end of September or early October for each year , the stake heights and snow-pit features (snow layer density and stratigraphy) are manually measured to derive the annual point mass balance. Then the glacier-wide mass balance was then calculated （Please to see the reference). Two automatic weather stations (AWSs, Campbell company) were installed near the Naimona’nyi Glacier. AWS1, at 5543 m a. s.l., recorded meteorological variables from October 2011 at half hourly resolution, including air temperature (℃), relative humidity (%), and downward shortwave radiation (W m-2) . AWS2 was installed at 5950 m a.s.l. in October 2010 at hourly resolution and recorded wind speed (m/s), air pressure (hPa), precipitation (mm). Data quality: the quality of the original data is better, less missing. Firstly, the abnormal data in the original records are removed, and then the daily values of these parameters are calculated. Two probes (Hobo MX2301) which record air temperature and relative humidity was installed on the glacier at half hour resolution since October 2018. The observed meteorological data was calculated as monthly values. The data is stored in Excel file. It can be used by researchers for studying the changes in climate, hydrology, glaciers, etc.