This data is precipitation data, which is the monthly precipitation product of tropical rainfall measurement mission TRMM 3b43. It integrates the main area of the Qinghai Tibet Plateau (25 ~ 40 ° n; 25 ~ 40 ° n); The precipitation data of 332 meteorological stations are from the National Meteorological Information Center of China Meteorological Administration. The reanalysis data set is obtained by the station 3 ° interpolation optimization variational correction method. For the monthly sample data from January 1998 to December 2018, the spatial coverage is 25 ~ 40 ° n; 73 ~ 105 ° e, the spatial resolution is 1 ° * 1 °.
XU Xiangde, SUN Chan
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
The 0.25 Degree climate data set in the northeastern part of the Tibetan Plateau from 1957 to 2009 contains four meteorological elements, which are precipitation, maximum and minimum temperatures, and wind speed. The time resolution is daily. The data set contains 2400 text files, each with precipitation (the 1st column), highest (the 2nd column) and lowest (the 3rd column) temperatures and wind speed (the 4th column). Each file name contains latitude and longitude. Each file represents the four meteorological element values for the corresponding grid point (0.25*0.25 degrees). These data are formed by gridding the observation data of 81 meteorological stations in the northeast of the plateau, considering the change of meteorological conditions with the elevation. The gridding methods and steps are as follows. Download the original daily maximum and minimum temperatures, precipitation, and wind speed from the China Meteorological Data Network (http://data.cma.cn). Then, perform quality control on the data. The principle used is 1) to remove daily precipitations below 0 and greater than 150 mm, daily temperatures below -50 °C and greater than 50 °C and wind speeds below 0 m / s, 2) draw annual sequence precipitation, temperature and wind speed, check for abnormal year-to-year changes, and conduct quality control through station migration records. For data with abnormal changes but with station migration records, the data are segmented by modifying the station name. For example, at Xining Station (52866), abnormal temperature changes occurred in 1996, which was found through records that Xining Station migrated after 1996. Therefore, the records before 1996 are recorded as virtual station 52867 data, and after 1996, the data are still recorded as 52866 stations. If the data change abnormally but there is no station migration record, the abnormally changed data are eliminated, for example, the data from Delingha Station before 1975. Some stations have migration records, but the data do not change abnormally; then, it is assumed that the stations before and after the migration are still in the same climate environment, so there is no change in station name and data record. Interpolation begins after quality control. The method begins with (1) calculating the changes in daily average temperature, precipitation and wind speed as the altitude changes. It is concluded that the temperature decreases with altitude by 4.3 °C/km, and the coefficient of determination R2 is 0.65. In the warm and humid season (from May to September), the average daily precipitation has an insignificant increase with altitude (0.5 mm/km, R2 is 0.1). The average daily precipitation in the cold dry season (from October to April) does not change with altitude. The wind speed also has an insignificant increase with altitude, with an increase rate of 0.4 m/s/km and R2 of 0.1. (2) The spatial interpolation is performed using the Synographic Mapping System (SYMAP, Shepard, 1984) method. In this method, the distance between stations and the angle between surrounding stations are taken into account in interpolation to indicate the density of stations. The distance and angle are integrated into a weight. In addition, the stations that are close and have a large angle between each other are given a large weight. (3) The latitude and longitude of the station, the meteorological element value, the altitude, the rate of change with the altitude, and the weight are considered simultaneously, and the value of the target grid is interpolated. The maximum search range for interpolation is 55 stations around, and the smallest search range is 4 stations around. (4) Integrate the precipitation in the warm and dry seasons to form a precipitation sequence throughout the period. (5) During the method test period, some stations are set aside to check the gridded data. (6) After the verification is passed, all 81 stations are used in the final gridding process and form this set of data sets. Shepard, D. S., 1984: Computer Mapping: The SYMAP interpolation algorithm. Spatial Statistics and Models, G.Gaile and C. Willmot, Eds., Reidel 133-145.
LAN Cuo
The atmospheric forcing dataset for along the Belt and Road from 2000 to 2015 comes from CRUNCEP. CRUNCEP is an atmospheric forcing dataset used forcing the land surface models. Specifically, this long time series data set (including temperature, precipitation, temperature, etc.) is used to drive the Community Land Model (CLM) Land Model in the long term. The CRUNCEP is a combination of two existing datasets; the CRU TS3.2 0.5 X 0.5 monthly data covering the period 1901 to 2002 and the NCEP reanalysis 2.5 X 2.5 degree 6-hourly data covering the period 1948 to 2016. The CRUNCEP dataset has been used to force CLM for studies of vegetation growth, evapotranspiration, and gross primary production and for the TRENDY (trends in net land-atmosphere carbon exchange over the period 1980-2010) project, among many other use cases. The CRUNCEP data archived in this dataset is Version 7.
The National Center for Atmospheric Research, CAO Wei
The long-time series data set of extreme precipitation index in the arid region of Central Asia contains 10 extreme precipitation index long-time series data of 49 stations. Based on the daily precipitation data of the global daily climate historical data network (ghcn-d), the data quality control and outlier elimination were used to select the stations that meet the extreme precipitation index calculation. Ten extreme precipitation indexes (prcptot, SDII, rx1day, rx5day, r95ptot, r99ptot, R10, R20) defined by the joint expert group on climate change detection and index (etccdi) were calculated 、CWD、CDD)。 Among them, there are 15 time series from 1925 to 2005. This data set can be used to detect and analyze the frequency and trend of extreme precipitation events in the arid region of Central Asia under global climate change, and can also be used as basic data to explore the impact of extreme precipitation events on agricultural production and life and property losses.
YAO Junqiang, CHEN Jing, LI Jiangang
Data description: This dataset includes the grid data of annual temperature and annual precipitation on the Tibetan Plateau from 1998 to 2017. It is the basic data for study of climate change and its impact on the ecological environment. Data source and processing: The meta data was aquired from the temperature and precipitation daily data of China's ground high-density stations (above 2,400 national meteorological stations) based on the latest compilation of the National Meteorological Information Center's basic data. After removing the missing stations, the software's thin plate spline method in ANUSPLIN was used to perform spatial interpolation, in order to generate grid data with spactial resolution of 1 km on the Tibetan Plateau . Data application: This data can be used to indentify the impact of climate change on the ecological environment.
DING mingjun
This data is the aridity index (AI) under the rcp4.5 scenario. AI data is the ratio of precipitation to potential evapotranspiration. This data is calculated by the average of 14 models. These 14 modes are canesm2; ccsm4; cnrm-cm5; csiro-mk3-6-0; giss-e2-r; hadgem2-cc; hadgem2-es; inmcm4; ipsl-cm5a-lr; miroc5; miroc-esm-chem; miroc-esm; mpi-esm-lr; mri-cgcm3. The spatial resolution is 2 * 2 degrees, and the temporal resolution is from January 2020 to December 2099. This data set can be used to analyze the future dry and wet change scenarios in the Great Lakes region of Central Asia, as well as the dry and wet past and pattern in other regions of the world under the future scenarios.
HUA Lijuan
Coupled Model Intercomparison Project Phase 5 (CMIP5) provides a multiple climate model environment, which can be used to predict the future climate change in the key nodes in the Belts and Road to deal with the environmental and climate problems. Key nodes in the Belt and Road are taken as the study regions of this dataset. The ability of 43 climate models in CMIP5 to predict the future climate change in the study regions was assessed and the optimal models under different scenarios were selected according to the RMSE between the prediction results and real observations. This dataset is composed of the prediciton results of precipitation and near-surface air temperature between 2006 and 2065 using the optimal models in monthly temporal frequncy. The spatial resolution of the dataset has been downscaled to 10 km using statistical downscaling method. Data of each period has three bands, namely maximum near-surface air temperature, minimum near-surface air temperature and precipitation. In this data set, the unit of precipitation is kg / (m ^ 2 * s), and the unit of near-surface air temperature is K. This dataset provides data basis for solving environmental and climate problems of the key nodes in the Belts and Road.
LI Xinyan, LING Feng
Precipitation estimates with fine quality and spatio-temporal resolutions play significant roles in understanding the global and regional cycles of water, carbon, and energy. Satellite-based precipitation products are capable of detecting spatial patterns and temporal variations of precipitation at fine resolutions, which is particularly useful over poorly gauged regions. However, satellite-based precipitation products are the indirect estimates of precipitation, inherently containing regional and seasonal systematic biases and random errors. Focusing on the potential drawbacks in generating Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) and its recently updated retrospective IMERG in the Tropical Rainfall Measuring Mission (TRMM) era (finished in July 2019), which were only calibrated at a monthly scale using ground observations, Global Precipitation Climatology Centre (GPCC, 1.0◦/monthly), we aim to propose a new calibration algorithm for IMERG at a daily scale and to provide a new AIMERG precipitation dataset (0.1◦/half-hourly, 2000–2015, Asia) with better quality, calibrated by Asian Precipitation – Highly Resolved Observational Data Integration Towards Evaluation of Water Resources (APHRODITE, 0.25◦/daily) at the daily scale for the Asian applications. Considering the advantages from both satellite-based precipitation estimates and the ground observations, AIMERG performs better than IMERG at different spatio-temporal scales, in terms of both systematic biases and random errors, over mainland China.
MA Ziqiang
Qiangyong glacier: 90.23 °E, 28.88° N, 4898 m asl. The surface is bedrock. The record contains data of 1.5 m temperature, 1.5 m humidity, 2 m wind speed, 2 m wind orientation, surface temperature, etc. Data from the automated weather station was collected using USB equipment at 19:10 on August 6, 2019, with a recording interval of 10 minutes, and data was downloaded on December 20, 2020. There is no missing data but a problem with the wind speed data after 9:30 on July 14, 2020 (most likely due to damage to the wind vane). Jiagang glacier: 88.69°E, 30.82°N, 5362 m asl. The surface is rubble and weeds. The records include 1.5 meters of temperature, 1.5 meters of humidity, 2 meters of wind speed, 2 meters of wind direction, surface temperature, etc. The initial recording time is 15:00 on August 9, 2019, and the recording interval is 1 minute. The power supply is mainly maintained by batteries and solar panels. The automatic weather station has no internal storage. The data is uploaded to the Hobo website via GPRS every hour and downloaded regularly. At 23:34 on January 5, 2020, the 1.5 meter temperature and humidity sensor was abnormal, and the temperature and humidity data were lost. The data acquisition instrument will be retrieved on December 19, 2020 and downloaded to 19:43 on June 23, 2020 and 3:36 on September 25, 2020. Then the temperature and humidity sensors were replaced, and the observations resumed at 12:27 on December 21. The current data consists of three segments (2019.8.9-2020.6.30; 2020.6.23-2020.9.25; 2020.12.19-2020.12.29), Some data are missing after inspection. Some data are duplicated in time due to recording battery voltage, which needs to be checked. The meteorological observation data at the front end of Jiagang mountain glacier are collected by the automatic weather station Hobo rx3004-00-01 of onset company. The model of temperature and humidity probe is s-thb-m002, the model of wind speed and direction sensor is s-wset-b, and the model of ground temperature sensor is s-tmb-m006. The meteorological observation data at the front end of Jianyong glacier are collected by the US onset Hobo u21-usb automatic weather station. The temperature and humidity probe model is s-thb-m002, the wind speed and direction sensor model is s-wset-b, and the ground temperature sensor model is s-tmb-m006.
ZHANG Dongqi
The RCM employed is the International Center for Theoretical Physics (ICTP) Regional Climate Model version 4 (RegCM4, Giorgi et al., 2012). The domain used is the Coordinated Regional Climate Downscaling Experiment (CORDEX) Phase II East Asia domain, covering whole of China and its surrounding East Asia areas. The model is run at 25 km gird spacing, with its standard configuration of 18 vertical sigma layers with a model top at 10 hPa. Configuration of the model follows Gao et al. (2016, 2017), with land cover data over China was updated as reported by Han et al. (2015) to better represent the realistic vegetation. The initial and lateral boundary conditions needed to drive RegCM4 are derived from the CMIP5 models of HadGEM2-ES (RCP4.5 pathways), and the data set include temperature and precipitation.
GAO Xuejie
Kara batkak glacier meteorological station in Western Tianshan Mountains, Kyrgyzstan (42 ° 9'46 ″ n, 78 ° 16'21 ″ e, 3280m). The observation data include hourly meteorological elements (temperature (℃), maximum temperature (℃), time of maximum temperature occurrence, minimum temperature (℃), minimum temperature occurrence time, 0.1mm hourly rainfall (mm), 0.5mm hourly rainfall (mm), instantaneous wind direction (°), instantaneous wind speed (M / S), 2-minute wind direction (°), 2-minute wind speed (M / s), 10 minute wind speed (M / s), maximum wind direction (°), maximum wind speed (°), maximum wind speed (°) Major wind speed (M / s), maximum wind speed time, wind direction (°), maximum wind speed (M / s), maximum wind speed time, maximum instantaneous wind speed and direction (°), maximum instantaneous wind speed (M / s), relative humidity (%), minimum relative humidity (%), occurrence time of minimum relative humidity, water pressure (HPA), dew point temperature (℃), air pressure (HPA), sea level pressure (HPA), maximum pressure (HPA) The time of the highest air pressure, the lowest pressure (HPA) and the lowest air pressure (time). Meteorological observation elements are processed into climatic data after accumulation and statistics, providing important data for planning, design and research of agriculture, forestry, industry, transportation, military, hydrology, medical and health care and environmental protection departments.
HUO Wen
This data set includes the temperature, precipitation, relative humidity, wind speed, wind direction and other daily values in the observation point of Kunsha Glacier. The data is observed from October 3, 2015 to September 19, 2017. It is measured by automatic meteorological station (Onset Company) and a piece of data is recorded every 2 hours. The original data forms a continuous time series after quality control, and the daily mean index data is obtained through calculation. The original data meets the accuracy requirements of China Meteorological Administration (CMA) and the World Meteorological Organization (WMO) for meteorological observation. Quality control includes eliminating the systematic error caused by the missing point data and sensor failure. The data is stored as an excel file.
ZHANG Yinsheng
This data set includes meteorological data observed by the carbon flux station in the Guoluo Army Ranch in Qinghai. The temporal coverage is from 2005 to 2009, and the temporal resolution is 1 day. Meteorological and carbon flux data observation methods: vorticity-related observation instruments were used for automatic recording; biomass observation method: harvest method, weighing in a 60-degree oven for 48 hours. Both carbon flux and meteorological data were automatically recorded by the instruments and manually checked. During the data observation process, the operation of the instrument and the selection of the observation objects were in strict accordance with professional requirements, and the data could be applied to plant leaf photosynthetic parameter simulation and productivity estimation. This data contains observation items as follows: Temperature °C Precipitation mm Wind speed m/s Soil temperature at 5 cm depth °C Photosynthetically active radiation µmol/m²s Total radiation W/m²
ZHAO Xinquan
This data set describes the temporal and spatial distribution of precipitation in the Upper Brahmaputra River Basin. We integrate (CMA, GLDAS, ITP-Forcing, MERRA2, TRMM) five sets of reanalysis precipitation products and satellite precipitation products, and combine the observation precipitation of 9 national meteorological stations from China Meteorological Administration (CMA) and 166 rain gauges of the Ministry of Water Resources (MWR) in the basin. The time range is 1981-2016, the time resolution is 3 hours, the spatial resolution is 5 km, and the unit is mm/h. The data will provide better data support for the study of Upper Brahmaputra River Basin, and can be used to study the response of hydrological process to climate change. Please refer to the instruction document uploaded with the data for specific usage information.
WANG Yuanwei, WANG Lei, LI Xiuping, ZHOU Jing
The Central Asia Reanalysis (CAR) dataset is generated based on the Weather Research and Forecast (WRF) model version 4.1.2 and WRF Data Assimilation (WRFDA) Version 4.1.2. Variables include temperature,, pressure, wind speed, precipitation and radiation. The reanalysis is established through cyclic assimilation, which performs data assimilation every 6 hours by 3DVAR. The assimilated data include conventional atmospheric observation and satellite radiation data. The main source of conventional data is Global Teleconnection System (GTS), including surface station, automatic station, radiosonde and aircraft report, and the observation elements include temperature, air pressure, wind speed and humidity. Satellite observations include retrievals and radiation data, The retrievals are mainly atmospheric motion vectors from polar orbiting meteorological satellites (NOAA-18, NOAA-19, MetOP-A and MetOP-B) and resampled to a horizontal resolution of 54km; the radiation data includes microwave radiation from MSU, AMSU and MHS and HIRS infrared radiation data. The simulation applies nesting with a horizontal resolution of 27km and 9km respectively, a total of 38 layers in the vertical direction and a top of the model layer of 10hPa. The lateral boundary conditions of the model are provided by ERA-Interim every 6 hours. The physical schemes used in the model are Thompson microphysics scheme, CAM radiation scheme, MYJ boundary layer scheme, Grell convection scheme and Noah land surface model. The data covers five countries in Central Asia, including Kazakhstan, Tajikistan, Kyrgyzstan, Turkmenistan and Uzbekistan, as well as lakes in Central Asia, such as Caspian Sea, Aral Sea, Balkash lake and Isaac lake, which can be used for the study of climate, ecology and hydrology in the region. Compared with gauge-based precipitation in Central Asia, the simulation by CAR shows similar performance with MSWEP ( a merged product) and outperforms ERA5 and ERA-Interim.
YAO Yao
Precipitation over the Tibetan Plateau (TP) known as Asia's water tower plays a critical role in regional water and energy cycles, largely affecting water availability for downstream countries. Rain gauges are indispensable in precipitation measurement, but are quite limited in the TP that features complex terrain and the harsh environment. Satellite and reanalysis precipitation products can provide complementary information for ground-based measurements, particularly over large poorly gauged areas. Here we optimally merged gauge, satellite, and reanalysis data by determining weights of various data sources using artificial neural networks (ANNs) and environmental variables including elevation, surface pressure, and wind speed. A Multi-Source Precipitation (MSP) data set was generated at a daily timescale and a spatial resolution of 0.1° across the TP for the 1998‒2017 period. The correlation coefficient (CC) of daily precipitation between the MSP and gauge observations was highest (0.74) and the root mean squared error was the second lowest compared with four other satellite products, indicating the quality of the MSP and the effectiveness of the data merging approach. We further evaluated the hydrological utility of different precipitation products using a distributed hydrological model for the poorly gauged headwaters of the Yangtze and Yellow rivers in the TP. The MSP achieved the best Nash-Sutcliffe efficiency coefficient (over 0.8) and CC (over 0.9) for daily streamflow simulations during 2004‒2014. In addition, the MSP performed best over the ungauged western TP based on multiple collocation evaluation. The merging method could be applicable to other data-scarce regions globally to provide high quality precipitation data for hydrological research. The latitude and longitude of the left bottom corner across the TP, the number of rows and columns, and grid cells information are all included in each ASCII file.
HONG Zhongkun , LONG Di
The precipitation dataset of the Third Pole region mainly contains two EXCEL files: (1) Daily precipitation data in China in the Third Pole region, named as China_daily.xlsx. The precipitation data in China were obtained from the China Meteorological Administration-National Meteorological Information Center (http://data.cma.gov.cn/site/index.html). (2) Daily precipitation data in other countries in the Third Pole region, named as Foreign_daily.xlsx. The precipitation data in other countries were obtained from NCDC International Climatic Data Center - NOAA Satellite Information Service Center (http://www7.ncdc.noaa.gov/CDO/country), Pakistan Meteorological Administration, Nepal Meteorological Administration, etc. There are seven variables in these two EXCEL data files: precipitation, corrected precipitation, correction factor, wind-induced loss, evaporation loss, wet loss, and trace precipitation. The detail characteristics of TPE stations were described in an EXCEL file either, named as "TPE station and gauge type.xls". The raw data has been strictly quality controlled by the relevant meteorological departments and has been applied in relevant academic papers.
ZHANG Yinsheng
Agricultural irrigation consumes a large amount of available freshwater resources and is the most immediate human disturbance to the natural water cycle process, with accelerated regional water cycles accompanied by cooling effects. Therefore, estimating irrigation water use (IWU) is important for exploring the impact of human activities on the natural water cycle, quantifying water resources budget, and optimizing agricultural water management. However, the current irrigation data are mainly based on the survey statistics, which is scattered and lacks uniformity, and cannot meet the demand for estimating the spatial and temporal changes of IWU. The Global Irrigation Water Use Estimation Dataset (2011-2018) is calculated by the satellite soil moisture, precipitation, vegetation index, and meteorological data (such as incoming radiation and temperature) based on the principle of soil water balance. The framework of IWU estimation in this study coupled the remotely sensed evapotranspiration process module and the data-model fusion algorithm based on differential evolution. The IWU estimates provided from this dataset have small bias at different spatial scales (e.g., regional, state/province and national) compared to traditional discrete survey statistics, such as at Chinese provinces for 2015 (bias = −3.10 km^3), at U.S. states for 2013 (bias = −0.42 km^3), and at various FAO countries (bias = −10.84 km^3). Also, the ensemble IWU estimates show lower uncertainty compared to the results derived from individual precipitation and soil moisture satellite products. The dataset is unified using a global geographic latitude and longitude grid, with associated metadata stored in corresponding NetCDF file. The spatial resolution is about 25 km, the time resolution is monthly, and the time span is 2011-2018. This dataset will help to quantitatively assess the spatial and temporal patterns of agricultural irrigation water use during the historical period and support scientific agricultural water management.
ZHANG Kun, LI Xin, ZHENG Donghai, ZHANG Ling, ZHU Gaofeng
Central Asia (referred to as CA) is among the most vulnerable regions to climate change due to the fragile ecosystems, frequent natural hazards, strained water resources, and accelerated glacier melting, which underscores the need of high-resolution climate projection datasets for application to vulnerability, impacts, and adaption assessments. We applied three bias-corrected global climate models (GCMs) to conduct 9-km resolution dynamical downscaling in CA. A high-resolution climate projection dataset over CA (the HCPD-CA dataset) is derived from the downscaled results, which contains four static variables and ten meteorological elements that are widely used to drive ecological and hydrological models. The static variables are terrain height (HGT, m), land use category (LU_INDEX, 21 categories), land mask (LANDMASK, 1 for land and 0 for water), and soil category (ISLTYP, 16 categories). The meteorological elements are daily precipitation (PREC, mm/day), daily mean/maximum/minimum temperature at 2m (T2MEAN/T2MAX/T2MIN, K), daily mean relative humidity at 2m (RH2MEAN, %), daily mean eastward and northward wind at 10m (U10MEAN/V10MEAN, m/s), daily mean downward shortwave/longwave flux at surface (SWD/LWD, W/m2), and daily mean surface pressure (PSFC, Pa). The reference and future periods are 1986-2005 and 2031-2050, respectively. The carbon emission scenario is RCP4.5. The results show the data product has good quality in describing the climatology of all the elements in CA, which ensures the suitability of the dataset for future research. The main feature of projected climate changes in CA in the near-term future is strong warming (annual mean temperature increasing by 1.62-2.02℃) and significant increase in downward shortwave and longwave flux at surface, with minor changes in other elements. The HCPD-CA dataset presented here serves as a scientific basis for assessing the impacts of climate change over CA on many sectors, especially on ecological and hydrological systems.
QIU Yuan QIU Yuan
Contact Support
Northwest Institute of Eco-Environment and Resources, CAS 0931-4967287 poles@itpcas.ac.cnLinks
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