The data set mainly includes the ice observation frequency (ICO) of north temperate lakes in four periods from 1985 to 2020, as well as the location, area and elevation of the lakes. Among them, the four time periods are 1985-1998 (P1), 1999-2006 (P2), 2007-2014 (P3) and 2015-2020 (P4) respectively, in order to improve the "valid observation" times in the calculation period and improve the accuracy. The ICO of the four periods is calculated by the ratio of "icing" times and "valid observation" times counted by all Landsat images in each period. Other lake information corresponds to the HydroLAKEs data set through the "hylak_id" column in the table. In addition, the data only retains about 30000 lakes with an area of more than 1 square kilometer, which are valid for P1-P4 observation. The data set can reflect the response of Lake icing to climate change in recent decades.
WANG Xinchi
This biophysical permafrost zonation map was produced using a rule-based GIS model that integrated a new permafrost extent, climate conditions, vegetation structure, soil and topographic conditions, as well as a yedoma map. Different from the previous maps, permafrost in this map is classified into five types: climate-driven, climate-driven/ecosystem-modified, climate-driven/ecosystem protected, ecosystem-driven, and ecosystem-protected. Excluding glaciers and lakes, the areas of these five types in the Northern Hemisphere are 3.66×106 km2, 8.06×106 km2, 0.62×106 km2, 5.79×106 km2, and 1.63×106 km2, respectively. 81% of the permafrost regions in the Northern Hemisphere are modified, driven, or protected by ecosystems, indicating the dominant role of ecosystems in permafrost stability in the Northern Hemisphere. Permafrost driven solely by climate occupies 19% of permafrost regions, mainly in High Arctic and high mountains areas, such as the Qinghai-Tibet Plateau.
RAN Youhua, M. Torre Jorgenson, LI Xin, JIN Huijun, Wu Tonghua, Li Ren, CHENG Guodong
Based on the analysis of brgdgts and hydrogen isotopes of leaf wax in lake sediments from Tengchong Qinghai (tcqh) in Yunnan Province, this study shows for the first time the high-resolution annual average temperature change history of low latitude land since the last glacial period (since the last 88000 years). According to the annual average temperature of South Asia established by tcqh core, there are two warm periods of 88000-71000 years and 45000-22000 years in this region, and the temperature range is about 2-3 ° C. Since the Holocene, the temperature has been increasing for about 1-2 years ° C。
ZHAO Cheng
A comprehensive understanding of the permafrost changes in the Qinghai Tibet Plateau, including the changes of annual mean ground temperature (Magt) and active layer thickness (ALT), is of great significance to the implementation of the permafrost change project caused by climate change. Based on the CMFD reanalysis data from 2000 to 2015, meteorological observation data of China Meteorological Administration, 1 km digital elevation model, geo spatial environment prediction factors, glacier and ice lake data, drilling data and so on, this paper uses statistics and machine learning (ML) method to simulate the current changes of permafrost flux and magnetic flux in Qinghai Tibet Plateau The range data of mean ground temperature (Magt) and active layer thickness (ALT) from 2000 to 2015 and 2061 to 2080 under rcp2.6, rcp4.5 and rcp8.5 concentration scenarios were obtained, with the resolution of 0.1 * 0.1 degree. The simulation results show that the combination of statistics and ML method needs less parameters and input variables to simulate the thermal state of frozen soil, which can effectively understand the response of frozen soil on the Qinghai Tibet Plateau to climate change.
Ni Jie, Wu Tonghua
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.
CAO Xianyong, TIAN Fang, LI Kai, NI Jian
These datasets include mean annual ground temperature (MAGT) at the depth of zero annual amplitude (approximately 3 m to 25 m), active layer thickness (ALT), the probability of the permafrost occurrence, and the new permafrost zonation based on hydrothermal condition for the period of 2000-2016 in the Northern Hemisphere with an 1-km resolution by integrate unprecedentedly large amounts of field data (1,002 boreholes for MAGT and 452 sites for ALT) and multisource geospatial data, especially remote sensing data, using statistical learning modelling with an ensemble strategy, and thus more accurate than previous circumpolar maps.
RAN Youhua, LI Xin, CHENG Guodong, CHE Jinxing, Juha Aalto, Olli Karjalainen, Jan Hjort, Miska Luoto, JIN Huijun, Jaroslav Obu, Masahiro Hori, YU Qihao, CHANG Xiaoli
We comprehensively estimated water volume changes for 1132 lakes larger than 1 km2. Overall, the water mass stored in the lakes increased by 169.7±15.1 Gt (3.9±0.4 Gt yr-1) between 1976 and 2019, mainly in the Inner-TP (157.6±11.6 or 3.7±0.3 Gt yr-1). A substantial increase in mass occurred between 1995 and 2019 (214.9±12.7 Gt or 9.0±0.5 Gt yr-1), following a period of decrease (-45.2±8.2 Gt or -2.4±0.4 Gt yr-1) prior to 1995. A slowdown in the rate of water mass increase occurred between 2010 and 2015 (23.1±6.5 Gt or 4.6±1.3 Gt yr-1), followed again by a high value between 2015 and 2019 (65.7±6.7 Gt or 16.4±1.7 Gt yr-1). The increased lake-water mass occurred predominately in glacier-fed lakes (127.1±14.3 Gt) in contrast to non-glacier-fed lakes (42.6±4.9 Gt), and in endorheic lakes (161.9±14.0 Gt) against exorheic lakes (7.8±5.8 Gt) over 1976−2019.
ZHANG Guoqing
Effective evaluation of future climate change, especially prediction of future precipitation, is an important basis for formulating adaptation strategies. This data is based on the RegCM4.6 model, which is compatible with multi-model and different carbon emission scenarios: CanEMS2 (RCP 45 and RCP85), GFDL-ESM2M (RCP2.6, RCP4.5, RCP6.0 and RCP8.5), HadGEM2-ES (RCP2.6, RCP4.5 And RCP8.5), IPSL-CM5A-LR (RCP2.6, RCP4.5, RCP6.0 and RCP8.5), MIROC5 (RCP2.6, RCP4.5, RCP6.0 and RCP8.5). The future climate data (2007-2099) has 21 sets, with a spatial resolution at 0.25 degrees and the temporal resolution at 3 hours (or 6 hours), daily and yearly scales.
PAN Xiaoduo, ZHANG Lei
This dataset is the monthly precipitation data of China, with a spatial resolution of 0.0083333 ° (about 1km) and a time range of 1901.1-2021.12. The data format is NETCDF, i.e.. Nc format. This dataset is generated in China through the Delta spatial downscaling scheme based on the global 0.5 ° climate dataset released by CRU and the global high-resolution climate dataset released by WorldClim. In addition, 496 independent meteorological observation point data are used for verification, and the verification results are reliable. This data set covers the main land areas in China (including Hong Kong, Macao and Taiwan), excluding islands and reefs in the South China Sea. In order to facilitate storage, the data are all int16 type and stored in nc files, with precipitation units of 0.1mm. NC data can be mapped using ArcMAP software; Matlab software can also be used for extraction processing. Matlab has released the function to read and store nc files. The read function is ncread, and switch to the nc file storage folder. The statement is expressed as: ncread ('XXX.nc ',' var ', [i j t], [leni lenj lent]), where XXX.nc is the file name, and is the string required' '; Var is from XXX The variable name read in NC. If it is a string, '' is required; i. J and t are the starting row, column and time of the read data respectively, and leni, lenj and lent i are the length of the read data in the row, column and time dimensions respectively. In this way, this function can be used to read in any region and any time period in the study area. There are many commands about NC data in the help of Matlab, which can be viewed. WGS84 is recommended for data coordinate system.
PENG Shouzhang
The field observation platform of the Tibetan Plateau is the forefront of scientific observation and research on the Tibetan Plateau. The land surface processes and environmental changes based comprehensive observation of the land-boundary layer in the Tibetan Plateau provides valuable data for the study of the mechanism of the land-atmosphere interaction on the Tibetan Plateau and its effects. This dataset integrates the 2005-2016 hourly atmospheric, soil hydrothermal and turbulent fluxes observations of Qomolangma Atmospheric and Environmental Observation and Research Station, Chinese Academy of Sciences (QOMS/CAS), Southeast Tibet Observation and Research Station for the Alpine Environment, CAS (SETORS), the BJ site of Nagqu Station of Plateau Climate and Environment, CAS (NPCE-BJ), Nam Co Monitoring and Research Station for Multisphere Interactions, CAS (NAMORS), Ngari Desert Observation and Research Station, CAS (NADORS), Muztagh Ata Westerly Observation and Research Station, CAS (MAWORS). It contains gradient observation data composed of multi-layer wind speed and direction, temperature, humidity, air pressure and precipitation data, four-component radiation data, multi-layer soil temperature and humidity and soil heat flux data, and turbulence data composed of sensible heat flux, latent heat flux and carbon dioxide flux. These data can be widely used in the analysis of the characteristics of meteorological elements on the Tibetan Plaetau, the evaluation of remote sensing products and development of the remote sensing retrieval algorithms, and the evaluation and development of numerical models.
MA Yaoming
Mean annual ground temperature (MAGT) at a depth of zero annual amplitude and permafrost thermal stability type are fundamental importance for engineering planning and design, ecosystem management in permafrost region. This dataset is produced by integrating remotely sensed freezing degree-days and thawing degree-days, snow cover days, leaf area index, soil bulk density, high-accuracy soil moisture data, and in situ MAGT measurements from 237 boreholes for the 2010s (2005-2015) on the Tibetan Plateau (TP) by using an ensemble learning method that employs a support vector regression (SVR) model based on distance-blocked resampling training data with 200 repetitions. Validation of the new permafrost map indicates that it is probably the most accurate of all available maps at present. The RMSE of MAGT is approximately 0.75 °C and the bias is approximately 0.01 °C. This map shows that the total area of permafrost on the TP is approximately 115.02 (105.47-129.59) *104 km2. The areas corresponding to the very stable, stable, semi-stable, transitional, and unstable types are 0.86*104 km2, 9.62*104 km2, 38.45*104 km2, 42.29*104 km2, and 23.80*104 km2, respectively. This new dataset is available for evaluate the permafrost change in the future on the TP as a baseline. More details can be found in Ran et al., (2020) that published at Science China Earth Sciences.
RAN Youhua, LI Xin
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
This data set comprises the plateau soil moisture and soil temperature observational data based on the Tibetan Plateau, and it is used to quantify the uncertainty of model products of coarse-resolution satellites, soil moisture and soil temperature. The observation data of soil temperature and moisture on the Tibetan Plateau (Tibet-Obs) are from in situ reference networks at four regional scales, which are the Nagqu network of cold and semiarid climate, the Maqu network of cold and humid climate, and the Ali network of cold and arid climate，and Pali network. These networks provided representative coverage of different climates and surface hydrometeorological conditions on the Tibetan Plateau. - Temporal resolution: 1hour - Spatial resolution: point measurement - Measurement accuracy: soil moisture, 0.00001; soil temperature, 0.1 °C; data set size: soil moisture and temperature measurements at nominal depths of 5, 10, 20, 40 - Unit: soil moisture, cm ^ 3 cm ^ -3; soil temperature, °C
BOB Su, YANG Kun
Soil bulk density, porosity, water content, water characteristic curve, saturated hydraulic conductivity, particle analysis, infiltration rate, and sampling point location information in the upper reaches of the Heihe River Basin. 1. The data is for 2014 supplementary sampling for 2012, using the ring knife to take the original soil; 2. The soil bulk density is the dry bulk density of the soil and is measured by the drying method. The original ring-shaped soil sample collected in the field was thermostated at 105 ° C for 24 hours in an oven, and the soil dry weight was divided by the soil volume (100 cubic centimeters) , unit: g/cm 3 . 3. Soil porosity is obtained according to the relationship between soil bulk density and soil porosity; 4. Soil infiltration analysis data set, the data is the field experimental measurement data from 2013 to 2014. 5. The infiltration data is measured by “MINI DISK PORTABLE TENSION INFILTROMETER”, and the approximate saturated hydraulic conductivity under a certain negative pressure is obtained. 6. Soil particle size data was measured at the Grain Granulation Laboratory of the Key Laboratory of the Ministry of Education of Lanzhou University. The measuring instrument is a Malvern laser particle size analyzer MS2000. 7. The saturated hydraulic conductivity is measured according to the enamel hair self-made instrument of Yi Yanli (2009). The Marioot bottle was used to maintain the head during the experiment; at the same time, the Ks measured at the time was converted to the Ks value at 10 °C for analysis and calculation. 8. Soil water content data is measured using ECH2O, including 5 layers of soil water content and soil temperature. 9. The water characteristic curve is measured by the centrifuge method: the undisturbed soil of the ring cutter collected in the field is placed in a centrifuge, and each of the speeds is measured at 0, 310, 980, 1700, 2190, 2770, 3100, 5370, 6930, 8200, 11600. The secondary rotor weight is obtained.
HE Chansheng
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
This data set is collected from the supplementary information part of the paper: Yao, T. , Thompson, L. , & Yang, W. . (2012). Different glacier status with atmospheric circulations in tibetan plateau and surroundings. Nature Climate Change, 1580, 1-5. This paper report on the glacier status over the past 30 years by investigating the glacial retreat of 82 glaciers, area reductionof 7,090 glaciers and mass-balance change of 15 glaciers. This data set contains 8 tables, the names and content are as follows: Data list: The data name list of the rest tables; t1: Distribution of Glaciers in the TP and surroundings; t2: Data and method for analyzing glacial area reduction in each basin; t3: Glacial area reduction during the past three decades from remote sensing images in the TP and surroundings; t4: Glacial length fluctuationin the TP and surroundings in the past three decades; t5: Detailed information on the glaciers for recent mass balance measurement in the TP and surroundings; t6: Recent annual mass balances in different regions in the TP; t7: Mass balance of Long-time series for the Qiyi, Xiaodongkemadi and Kangwure Glaciers in the TP. See attachments for data details: Supplementary information.pdf, Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings.pdf.
YAO Tandong
DEM is the English abbreviation of Digital Elevation Model, which is the important original data of watershed topography and feature recognition.DEM is based on the principle that the watershed is divided into cells of m rows and n columns, the average elevation of each quadrilateral is calculated, and then the elevation is stored in a two-dimensional matrix.Since DEM data can reflect local topographic features with a certain resolution, a large amount of surface morphology information can be extracted through DEM, which includes slope, slope direction and relationship between cells of watershed grid cells, etc..At the same time, the surface flow path, river network and watershed boundary can be determined according to certain algorithm.Therefore, to extract watershed features from DEM, a good watershed structure pattern is the premise and key of the design algorithm. Elevation data map 1km data formed according to 1:250,000 contour lines and elevation points in China, including DEM, hillshade, Slope and Aspect maps. Data set projection: Two projection methods: Equal Area projection Albers Conical Equal Area (105, 25, 47) Geodetic coordinates WGS84 coordinate system
TANG Guoan
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
Qinghai Tibet Plateau is the largest permafrost area in the world. At present, some permafrost distribution maps have been compiled. However, due to the limited data sources, unclear standards, insufficient verification and lack of high-quality spatial data sets, there is great uncertainty in drawing Permafrost Distribution Maps on TP. Based on the improved medium resolution imaging spectrometer (MODIS) surface temperature (LSTS) model of 1 km clear sky mod11a2 (Terra MODIS) and myd11a2 (Aqua MODIS) product (reprocessing version 5) in 2003-2012, the data set simulates the distribution of permafrost and generates the permafrost map of Qinghai Tibet Plateau. The map was verified by field observation, soil moisture content and bulk density. Permafrost attributes mainly include: seasonally frozen ground, permafrost and unfrozen ground. The data set provides more detailed data of Permafrost Distribution and basic data for the study of permafrost in the Qinghai Tibet Plateau.
ZHAO Lin
This dataset contains five types of boundaries. 1. TPBoundary_ 2500m: Based on ETOPO5 Global Surface Relief, ENVI+IDL was used to extract data at an elevation of 2500m within the longitude (65~105E) and latitude (20~45N) range in the Tibetan Plateau. 2. TPBoundary_ 3000m: Based on ETOPO5 Global Surface Relief, ENVI+IDL was used to extract data at an elevation of 3000m within the longitude (65~105E) and latitude (20~45N) range in the Tibetan Plateau. 3. TPBoundary_ HF (high_frequency): This boundary is defined according to 2 previous studies. Bingyuan Li (1987) had a systematic discussion on the principles for determining the extent of the Tibetan Plateau and the specific boundaries. From the perspective of the formation and basic characteristics of the Tibetan Plateau, he proposed the basic principles for determining the extent of the Tibetan Plateau based on the geomorphological features, the plateau surface and its altitude, while considering the integrity of the mountain. Yili Zhang (2002) determined the extent and boundaries of the Tibetan Plateau based on the new results of research in related fields and years of field practice. He combined information technology methods to precisely locate and quantitatively analyze the extent and boundary location of the Tibetan Plateau, and concluded that the Tibetan Plateau in China extends from the Pamir Plateau in the west to the Hengduan Mountains in the east, from the southern edge of the Himalayas in the south to the northern side of the Kunlun-Qilian Mountains in the north. On April 14, 2017, the Ministry of Civil Affairs of the People's Republic of China issued the Announcement on Adding Geographical Names for Public Use in the Southern Tibetan Region (First Batch), adding six geographical names in the southern Tibetan region, including Wo’gyainling, Mila Ri, Qoidêngarbo Ri, Mainquka, Bümo La, and Namkapub Ri. ４. TPBoundary_ New (2021): Along with the in-depth research on the Tibetan Plateau, the improvement of multidisciplinary research and understanding inside and outside the plateau, and the progress of geographic big data and Earth observation science and technology, the development of the 2021 version of the Tibetan Plateau boundary data by Yili Zhang and et al. was completed based on the comprehensive analysis of ASTER GDEM and Google Earth remote sensing images. The range boundary starts from the northern foot of the West Kunlun Mountain-Qilian Mountain Range in the north and reaches the southern foot of the Himalayas and other mountain ranges in the south, with a maximum width of 1,560 km from north to south; from the western edge of the Hindu Kush Mountains and the Pamir Plateau in the west to the eastern edge of the Hengduan Mountains and other mountain ranges in the east, with a maximum length of about 3,360 km from east to west; the latitude and longitude range is 25°59′30″N~40°1′0″N, 67°40′37″E~104°40′57″E, with a total area of 3,083,400km2 and an average altitude of about 4,320m. Administratively, the Tibetan Plateau is distributed in nine countries, including China, India, Pakistan, Tajikistan, Afghanistan, Nepal, Bhutan, Myanmar, and Kyrgyzstan. ５. TPBoundary_ Rectangle: The rectangle was drawn according to the range of Lon (63~105E) and Lat (20~45N). The data are in latitude and longitude projection WGS84. As the basic data, the boundary of the Tibetan Plateau can be used as a reference basis for various geological data and scientific research on the Tibetan Plateau.
ZHANG Yili
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
