The Antarctic and Arctic bacterial distribution data set provides distribution characteristics of bacteria in the Arctic and Antarctic. The collection period of the samples was from December 13,2005, to December 8,2006; 52 samples were obtained from 3 Arctic regions (Spitsbergen Slijeringa, Spitsbergen Vestpynten, and Alexandra Fjord_Highlands), and 171 samples were obtained from 5 Antarctic regions (the Mitchell Peninsula, Casey station main Power house, Robinsons Ridge, Herring Island, and Browning Peninsula). The soil surface samples were stored in liquid nitrogen after collection, shipped to a Sydney laboratory, and extracted using the FastPrep DNA kit. The extracted DNA samples were processed by 27F (5'-GAGTTTGATCNTGGCTCA-3' and 519R (5'-GTNTTACNGCGGCKGCTG-3') to amplify the 16S rRNA gene fragments. The amplified fragments were sequenced by the 454 method, and the raw data were analyzed by Mothur software. First, the sequences with poor sequencing quality were removed, the sequences were then sorted, and the chimera sequences were removed. The similarities between the sequences were calculated, the sequences with similarities above 97% were clustered into one OTU, and the OTU representative sequence was defined. By comparison with the Silva database, the OTU sequences with reliabilities greater than 80% were identified as level one. This data system compared the diversity of microorganisms in the eastern Antarctic with that in the Arctic and is of great significance for the study of the distributions of microorganisms in the Antarctic and Arctic.

As the main parameter in the land surface energy balance, surface temperature indicates the degree of land-atmosphere energy and water transfer and is widely used in research on climatology, hydrology and ecology. In the study of frozen soil, climate is one of the decisive factors for the existence and development of frozen soil. The surface temperature is the main climatic factor affecting the distribution of frozen soil and affects the occurrence, development and distribution of frozen soil. It is the upper boundary condition for modelling frozen soil and is significant to the study of hydrological processes in cold regions. The data set was based on the DEM and observation station data of the Tibetan Plateau Engineering Corridor and analysed the changing trend of surface temperature on the Tibetan Plateau from 2000 to 2014. Using the surface temperature data products MOD11A1/A2 and MYD11A1/A2 of MODIS aboard Terra and Aqua, the surface temperature information under cloud cover was reconstructed based on the spatio-temporal information of the images. The reconstruction information and surface temperature representativeness problems were analysed using information obtained from 8 sites, including the Kunlun Mountains (wetland, grassland), Beiluhe (grassland, meadow), Kaixinling (meadow, grassland), and Tanggula Mountain (meadow, wetland). According to the correlation coefficient (R2), root-mean-square error (RMSE), mean absolute error (MAE) and mean deviation (MBE), the following results were obtained: (1) the reconstruction accuracy of MODIS surface temperature under cloud cover is higher when it is based on spatio-temporal information; (2) the weighted average representation is the best when generalizing four observations of Terra and Aqua. By analysing the reconstruction of MODIS surface temperature information and representativeness problems, the average annual MODIS surface temperature data of the Tibetan Plateau and the engineering corridor from 2000 to 2010 were obtained. According to the data set, the surface temperature from 2000 to 2010 also experienced volatile rising trends from 2000 to 2010, which is basically consistent with the changing trend of the climate change in the permafrost regions of the Tibetan Plateau and the Qinghai-Tibet Engineering Corridor.

The variation in the duration of snow on the Tibetan Plateau is relatively great, and the high mountainous areas around the plateau are rich in snow and ice resources. Taking full account of the terrain of the Tibetan Plateau and the snow characteristics in the mountains, the data set adopted AVHRR data to gradually realize generating data products for daily, ten-day, and monthly snow cover areas while maintaining the snow classification accuracy. These data included the daily/10-day/monthly snow cover area data for the Tibetan Plateau from 2007 to 2015, the average accuracy of which is 0.92. It can provide reliable data for snow changes during the historical periods of the Tibetan Plateau.

The long-time series data set of snow cover area on the qinghai-tibet plateau is derived from the fusion of MODIS 005 version and IMS data set, andThe cloud-free products of daily snow cover area were obtained by using interpolation de-cloud algorithm.The projection is latitude and longitude, the spatial resolution is 0.005 degrees (about 500m), and the time is a long time series from January 1, 2003 to December 31, 2014. Each file is the result of the proportion of snow cover area on that day, and the value is 0-100 (%). It is the ENVI standard file, The naming convention： ims_mts_yyyyddd.tif, where YYYY stands for year and DDD stands for Julian day (001-365/366).Files can be directly used ENVI or ARCMAP software open view. Document description: 200 snow, 100 lake ice, 25 land, 37 sea

The past frozen soil map of the Tibetan Plateau was based on a small number of temperature station observations and used a classification system based on continuity. This data set used the geographically weighted regression model (GWR) to synthesize MODIS surface temperature, leaf area index, snow cover ratio and multimodel soil moisture forecast products of the National Meteorological Information Center through spatiotemporal reconstruction. In addition, precipitation observations of more than 40 meteorological stations, the precipitation products of FY2 satellite observations and the multiyear average temperature observation data of 152 meteorological stations from 2000 to 2010 were integrated to simulate the average temperature data of the Tibetan Plateau, and the permafrost thermal condition classification system was used to classify permafrost into several types: Very cold, Cold, Cool, Warm, Very warm, and Likely thawing. The map shows that, after deducting lakes and glaciers, the total area of permafrost on the Tibetan Plateau is approximately 1,071,900 square kilometers. Verification shows that this map has higher accuracy. It can provide support for future planning and design of frozen soil projects and environmental management.

Taking 2005 as the base year, the future population scenario prediction adopted the Logistic model of population, and it not only can better describe the change pattern of population and biomass but is also widely applied in the economic field. The urbanization rate was predicted using the urbanization Logistic model. Based on the existing urbanization horizontal sequence value, the prediction model was established by acquiring the parameters in the parametric equation applying nonlinear regression. The urban population was calculated by multiplying the predicted population by the urbanization rate. The Logistic model was used to predict the future gross national product of each county (or city), and then, according to the economic development level of each county (or city) in each period (in terms of real GDP per capita),the corresponding industrial structure scenarios in each period were set, and each industry’s output value was predicted. The trend of changes in industrial structure in China and the research area lagged behind the growth of GDP, and, therefore, it was adjusted according to the need of the future industrial structure scenarios of the research area.

The Karuola Glacier of Tibet is located at the junction of Langkazi County, the Shannan Area of the Tibetan Autonomous Region and Jiangzi County of the Shigatse Region. Latitude: 28°54'23.30′′～28°56'50.95′′N, Longitude 90°11′42.21′′～90°09′26.23′′E. It is a continental glacier with an average elevation of 5042 meters. It is the north-south spreading part of the Ningjingangsang peak. Based on the integration of the first glacier inventory data of China from the Cold and Arid Regions Environmental and Engineering Research Institute，Chinese Academy of Sciences, the 1:100,000 inventory data of the Yarlu Zangbu River Basin Glacier from the Sharing Platform for the Earth Systematic Science Data, and Google Earth remote sensing image and field survey data, the dataset was obtained with the help of ArcGIS, ENVI and other software by the following steps: first, the research and development of the data was achieved by band combination, research area clipping, manual visual interpretation and other techniques, and the accuracy of the obtained data was then verified. This dataset includes a total of 25 statistics of vector and area data of Tibet’s Karuola Glacier. It recorded the changes at the borders of Karuola Glacier in the past 45 years and could be used as reference data for the study of glacier and climate changes on the Tibetan Plateau.

By applying Supply-demand Balance Analysis, the water resource supply and demand of the whole river basin and each county or district were calculated, based on which the vulnerability of the water resources system of the basin was evaluated. The IPAT equation was used to set a future water resource demand scenario, setting variables such as future population growth rate, economic growth rate, and unit GDP water consumption to establish the scenario. By taking 2005 as the base year and using assorted forecasting data of population size and economic scale, the future water demand scenarios of various counties and cities from 2010 to 2050 were forecast. By applying the basic structure of the HBV conceptual hydrological model of the Swedish Hydrometeorological Institute, a model of the variation tendency of the basin under climate change was designed. The glacial melting scenario was used as the model input to construct the runoff scenario under climate change. According to the national regulations of the water resources allocation of the basin, a water distribution plan was set up to calculate the water supply comprehensively. Considering of the supply and demand situation, the water resource system vulnerability was evaluated by the water shortage rate. By calculating the (grain production) land pressure index of the major counties and cities in the basin, the balance of supply and demand of land resources under the climate change, glacial melt and population growth scenarios was analyzed, and the vulnerability of the agricultural system was evaluated. The Miami formula and HANPP model were used to calculate the human appropriation of net primary biomass and primary biomass in the major counties and cities for the future, and the vulnerability of ecosystems from the perspective of supply and demand balance was assessed.

Taking 2000 as the base year, the future population scenario prediction adopted the Logistic model of population, and it not only can better describe the change pattern of population and biomass but also is widely applied in the economic field. The urbanization rate was predicted using the urbanization Logistic model. Based on the existing urbanization horizontal sequence value, the prediction model was established by acquiring the parameters in the parametric equation applying nonlinear regression. The urban population was calculated by multiplying the predicted population by the urbanization rate. The Logistic model was used to predict the future gross national product of each county (or city), and then, according to the economic development level of each county (or city) in each period (in terms of real GDP per capita), the corresponding industrial structure scenarios in each period were set, and the output value of each industry was predicted. The trend of industrial structure changing in China and the research area lagged behind the growth of GDP, so it was adjusted according to the need of the future industrial structure scenarios of the research area.

Based on the Global 1,000,000 Basic Geographic Data (2010) of the Resource and Environment Science Data Center of the Chinese Academy of Sciences, the railway and highway networks of Arctic countries (USA, Canada, Russia, Norway (including Greenland and the Faroe Islands), Denmark, Sweden, Finland, and Iceland) are extracted via ArcGIS. The data are stored separately by nation. The data format is the .shp format of ArcGIS, and the projection mode is GCS_WGS_1984. The railway network data are from http://www.resdc.cn/data.aspx?DATAID=208, and the highway network data are from http://www.resdc.cn/data.aspx?DATAID=207

Arctic administrative boundary data sets include Arctic_National, Arctic_Provincial, and Arctic_Prefecture vector spatial data sets of arcti-bound countries and Its corresponding name, TYPE related attribute data :(LOCAL_NAME), (ENG_NAME), (CNTRY_NAME), (TYPE), (CNTRY_CODE), (CONTINENT) The data comes from the 1:1,000,000 ADC_WorldMap global data set, which is a comprehensive, up-to-date and seamless geographic digital data. The world map coordinate system is latitude and longitude, WGS84 datum surface, and the arctic data set is the special projection parameter for the arctic (North_Pole_Stereographic).

This data set is output from WRF model. The data include ‘LU_INDEX’ (land use category), ‘ZNU’(eta values on half (mass) levels), ‘ZNW’(eta values on full (w) levels)，’ZS’(depths of centers of soil layers), ‘DZS’ (thicknesses of soil layers), ‘VAR_SSO’ (variance of subgrid-scale orography), ‘U’(x-wind component), ‘V’(y-wind component),’W’(z-wind component),’T’(perturbation potential temperature (theta-t0)), ‘Q2’ ('QV at 2 M), ‘T2’ (TEMP at 2 M), ‘TH2’ ('POT TEMP at 2 M), ‘PSFC’ (SFC pressure), ‘U10’ (U at 10 M), ‘V10’ (V at 10 M), ‘QVAPOR’ (Water vapor mixing ratio), ‘QLOUD’ (Cloud water mixing ratio),’QRAIN’ (Rain water mixing ratio), ‘QICE’ (Ice mixing ratio), ‘QSNOW’ (Snow mixing ratio), ‘SHDMAX’ (annual max veg fraction), ‘SHDMIN’ (annual min veg fraction), ‘SNOALB’ (annual max snow albedo in fraction), ‘TSLB’ (soil temperature), ‘SMOIS’ (soil moisture), ‘GRDFLX’ (ground heat flux), ‘LAI’ (Leaf area index),’ HGT’ (Terrain Height), ‘TSK’ (surface skin temperature), ‘SWDOWN’ (downward short wave flux at ground surface), ‘GLW’ (downward long wave flux at ground surface), ‘HFX’ (upward heat flux at the surface), ‘QFX’ (upward moisture flux at the surface), ‘LH’ (latent heat flux at the surface), ‘SNOWC’ (flag indicating snow coverage (1 for snow cover)), and so on. The data is in netCDF format with a spatial resolution of 10 km.