#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Simulating the route of the Tang-Tibet Ancient Road for one branch of the Silk Road across the Qinghai-Tibet Plateau


Authors: Zhuoma Lancuo aff001;  Guangliang Hou aff001;  Changjun Xu aff002;  Yuying Liu aff003;  Yan Zhu aff001;  Wen Wang aff004;  Yongkun Zhang aff004
Authors place of work: Key Laboratory of Physical Geography and Environmental Process, College of Geography, Qinghai Normal University, Xining, Qinghai Province, China aff001;  Key Laboratory of Geomantic Technology and Application of Qinghai Province, Provincial geomantic Center of Qinghai, Xining, Qinghai Province, China aff002;  Department of computer technology and application, Qinghai University, Xining, Qinghai Province, China aff003;  State Key Laboratories of Plateau Ecology and Agriculture, Qinghai University, Xining, Qinghai Province, China aff004
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0226970

Summary

As the only route formed in the inner Qinghai-Tibet plateau, the Tang-Tibet Ancient Road promoted the extension of the Overland Silk Roads to the inner Qinghai-Tibet plateau. Considering the Complex geographical and environmental factors of inner Qinghai-Tibet Plateau, we constructed a weighted trade route network based on geographical integration factors, and then adopted the principle of minimum cost and the shortest path on the network to simulate the ancient Tang-Tibet Ancient Road. We then compared the locations of known key points documented in the literature, and found a significant correspondence in the Qinghai section. However, there was a certain deviation between the key points recorded in Tibetan section and the simulated route; we found that the reason is the relative oxygen content (ROC) became a limited factor of the choice of the Tibetan section road. Moreover, we argue that the warm and humid climate and the human migration to the hinterland of the Qinghai-Tibet plateau were the fundamental driving forces for the formation of the Tang-Tibet Ancient Road.

Keywords:

Roads – Forests – Temperate forests – Population density – Rivers – Deserts – Plateaus – Tibetan people

Introduction

The Belt and Road Initiative is designed to bean open, international network of cooperation and is becoming a multi-nation platform for exploring new mechanisms of international economic governance [1]. More specifically, the Overland Silk Roads (OSR) is a complex trade network that was formed by uniting large trading centers along with ancient transportation routes. Through this network, China has conducted extensive cultural and economic exchanges with Central, West, and South Asia [26]. The Qinghai-Tibet Plateau (QTP) is a key region connecting all of East, Central, West, and South Asia; owing to this geographical fact, this region has become the core of the Belt and Road Initiative [78]. Due to the harsh natural environment of the QTP, early trade routes existed around the margins of the QTP and historical and archaeological evidence (e.g., human activities, wheat and barley crops) from along these trade routes can be traced back to prehistoric times[914]. Using archeological evidence like bronze-handled mirrors (Archaeologists determined that such bronze-handled mirrors were introduced into the Tibetan Plateau in the early metal ages from West Asia and Central Asia. This is the evidence of communication between China and the west.), silk fabrics, tea leaves and Eurasian steppe-style rock painting, researchers have inferred that during the Han Dynasty (202 BC), people from the Southwest QTP reached the Silk Road through Xinjiang province (Yutian county) [1518]. People from the Northeast QTP reached the Silk Roads via other roads, such as those going through Gansu province (Zhangye county and Tunhuang county) and Xinjiang province (Ruoqiang county) (Fig 1) [1921]. These routes were all located on the geographical edge of the QTP. Up until the Tubo Dynasty (TBD) (618 AD-842 AD), these routes were important trade routes constituting OSR in the QTP [2226].

Trade routes on the edge of the Qinghai-Tibet Plateau connected with the Overland Silk Roads before the existence of the Tang-Tibet Ancient Road.
Fig. 1. Trade routes on the edge of the Qinghai-Tibet Plateau connected with the Overland Silk Roads before the existence of the Tang-Tibet Ancient Road.
Note: A, Evidence of Eurasian steppe style rock painting (Ritu Rock Art Sites in Ngari); B, Evidence of bronze mirrors with handles (Quta Cemetery in Ngari); C, Evidence of silk and tea (Gurgyam Cemetery in Ngari); D, Evidence of early metal (Qugong site in Lhasa).

Before the formation of Tang-Tibet Ancient Road (TTAR), both Qinghai and Tibet areas (today’s Tibet and Qinghai province) had their own separate routes to the OSR. However, there was no direct trade route in the inner plateau that could connect these two marginal areas. Given this, the TTAR is the only historically recorded Silk Road across the hinterland of the plateau. The road is approximately 2,240 km long, starting from Chang’an(today’s Xi'an province) and ending in Lhasa. Chen (1987) reported that the establishment of the TTAR connected the Xi'an-Qinghai and Tibet regions, which then linked OSR with the inner Plateau [27]. Hence, studies of the TTAR are crucial for a better understanding of the formative drivers and evolutionary processes of trade routes in the inner Plateau.

Currently, spatial data and geographic information systems (GIS) are widely applied to restore ancient trade routes [2829]. For instance, Ma et al. (2017) used 30 m resolution Shuttle Radar Topography Mission (SRTM) data and adopted an optimal path analysis method to reconstruct the ancient OSR from Chang’an to Constantinople across the Eurasian continent [30]. Frachettiet al. (2017) used a flow accumulation model to reveal a high-resolution flow network that simulated how centuries of seasonal nomadic herding shaped discrete routes of connectivity across the OSR [31]. When combining these previous studies, these renderings relied primarily on terrain-based least cost to predict likely routes. This approach was effective in low altitudes OSR areas where economic networks and mobility between sites were consistent with human ease of travel [3233]. However, this method may not be effective for the study of TTAR located in the inner QTP areas.

Previous studies on TTAR in the QTP areas relied primarily on historical literature and archaeological discoveries to determine node locations [3438]. However, the path from node to node remains relatively unknown. Unfortunately, a better understanding of these paths is difficult, as the names of the nodes are quite different across historical records and modern literals due to linguistic changes. In particular, the lack of historical records has made it difficult to ascertain the route of the Tibetan section of the TTAR. Moreover, the complex and changeable environments (e.g., anoxia, low temperature, vast rivers, low vegetation primary production, and sparse population) of the plateau have increased the difficulty of obtaining high-resolution images of ancient Tibetan road paths. Therefore, it is difficult to obtain complete, high-resolution road maps of the TTAR.

In this study, we constructed a cost-response data set based on 90 m DEM (Digital Elevation Model) data (including slope, altitude, and relief data)[39], as well as data regarding river classification [40], accumulated temperature [41], vegetation [42], and population density [43]. We then established a weighted trade route network based node set. Based on this, a complete and high-resolution route of the TTAR was simulated using an optimal path analysis method. Finally, the accuracy and reliability of the simulation results were verified using historical records. This study sought to reveal the formative factors of ancient key trade routes inside the QTP, broaden our understanding of the relationship between the TTAR and the OSR, and strengthen human knowledge of extreme environments such as high altitudes and the resulting hypoxic conditions.

Materials and methods

Study area

The QTP has an average altitude of 4,400 m and covers an estimated area of 250 ×104 km2. The latitude is from 26°00′12″N to 39°46′50″N and its longitude is from 73°18′52″E to 104°46′59″E [44]. The altitude rises rapidly from approximately 2,000 m in the northeast to above 4,000 m in the hinterland of the plateau, forming a natural barrier for human entry. The primary landforms of the plateau are mountains, plains, and basins. The plateau includes the source of the Three-River areas—the Yangtze, Yellow, and the Lancang Rivers. The average temperatures of the coldest and hottest months are -10°C ~ -15°C and < 10°C, respectively. There are also dramatic diurnal and seasonal temperature variations [4546]. The highest value of total solar radiation exceeds 8500 MJ/ (m2·a), and most parts of the plateau exceed 5000 MJ/ (m2·a) [47]. Along the TTAR in the QTP, vegetation zone communities mainly consist of subalpine forests, alpine meadows, alpine steppes, and alpine desert [48]. Low temperatures lead to insufficient soil nutrients, weak erosion resistance, slow plant growth, and low natural productivity [49]. In summary, our study area is characterized as an extreme natural environment that in ancient times presented a geographical barrier for human communication.

Methods

In this study, we assumed that Xining and Lhasa were two separate trade centers, with the nodes along the route between the two centers functioning as courier stations. In accordance with these rules, we constructed a network Eq (1), which defined a node (courier stations) as the node of the network, and the route connecting the nodes as the side of the network. We then abstracted the problem into a path search algorithm on the weighted network. Based on the principle of minimum cost, the shortest path of the network was used to simulate the TTAR (starting node Vs [Xining] and the ending node Ve [Lhasa]) (Fig 2).


Where V is the node set, E is the network connection, and W is the weight. The method of this paper sought to solve the following four main questions: (1) collection of a network node set; (2) generation of network connection; (3) creation of a weighted model; (4) search for the shortest path.

Flowchart for simulating the route of the Tang-Tibet Ancient Road.
Fig. 2. Flowchart for simulating the route of the Tang-Tibet Ancient Road.

Collection of a network node set

In this study, we assumed that nodes were courier stations along the aforementioned route. Therefore, we asserted the following types as nodes.

  • The site data showed direct evidence of the existence of ancient human activities. It is possible there was communication within and between the sites. The condition set was all sites data in the study area, which totaled to 512 [50] (“512” refers to the total number of archaeological sites found in the study area during the Tubo Dynasty (618 AD-842 AD) Available from: Atlas of Chinese Cultural Relics: Tibet Volume and Qinghai Volume.)

  • The criteria used for acceptable literature records regarding the TTAR selection used in this paper included the following: a chronology that spanned from ancient to modern literature documentation, literature types with official historical record, field investigation reports, and private travel notes. A total of 43 documented nodes were sorted; these node GPS data sets are shown in S1 Table. (Available from: documents [5155])

  • Communication between inhabited settlements was the most frequent communication type in ancient times. Both modern and ancient settlements have a degree of historical continuity, but the scale of large and medium-sized modern settlements is quite different from those during ancient times. The number of modern small settlements may not be completely consistent with that of ancient times. After all, the increase of modern population is unsurpassed in any historical period, and the rapid increase of population requires human to develop more areas suitable for human survival, Because of this the number of small settlements must be far more than that of historical period. But now the small settlements must contain all the plateau settlements of the historical period, which also include the fixed central camps (Nomadic peoples settlements). Therefore, we used modern, township-level administrative settlements as node data. We totally collected 1,792 data points regarding the geographic location data of the township-level administrative centers(National Earth System Science Data Sharing Platform. The geographical location data of administrative centers at the township level 2012, Database: [Internet]. Available from: http://www.tianditu.gov.cn/[56])

Generation of network connection

This study was based on the following conditions: all nodes of a node set were connected, communication routes could be formed between each node, and that there was no intersection in the communication route. Accordingly, a spatial data model of the Triangulate Irregular Network (TIN) was used to create the communication route network [57]. The connected nodes formed a triangle in line with the Delaunay standard, i.e. the collection conformed to a series of connected—but non-overlapping—triangles, and the outer circle of those triangles did not contain any other points in the polygon. We required that the maximum TIN margin length of each node did not exceed the calculated value of 22 km. Of note, the sighting of the TBD courier stations was usually 22 km (Note: During the Tang Dynasty and the Tubo Dynasty, the communication between the central government and the local governments was completed through special postal caravanserais system. The post caravanserais established at that time had special mileage regulations, which roughly translate into the present mileage of 22 km. In other words, the government stipulated that the distance between the post stations was equivalent to 22 km [58]). This study combines the actual mileage of post stations in the Tang Dynasty and Tubo Dynasty when setting the maximum side length of the model. Therefore, we required that the maximum TIN margin length of each node did not exceed the calculated value of 22 km). The step adapts the 3D analyst tool ‘data management-TIN’ in ArcGIS10.5version (ESRI Corporation, Redlands, CA, USA).Because the connection route was on the ground plane, the field with altitude value (z) was defined as 0 when the node parameter was set and the final output was a spatial network that consisted of edges.

Creation of a weighted model

In this study, the creation of a weighted model sought to solve the following three main questions: (1) calculate the weight of each edge; (2) calculate the unit cost of a raster; (3) calculate the length of edge through raster.

1. Calculation of the weight of each edge: The weight of each edge was determined as follows: if: p1p2V and there are edges (existing routes) between p1 and p2, the edge is ep1p2=(p1,p2)∈E, the cost function at coordinate (x,y) is f(x,y), and then the total cost of edges can be expressed as follows Eq (2):


If ep1p2 passes through m rasters, and the inside length of each raster is L(i)(i = 1,2,3,⋯m), then the unit cost of edges (Wep1p2) and the unit cost of raster(S(i)) and can be expressed as follows Eqs (3 and 4):



Where: S(i) is the cost estimate of unit raster i, cj is the weight of the parameter j and fj(i) is the weight of the parameter j in raster i.

Calculation of the unit cost of a raster:The unit cost of a grid was calculated using formula 4 with a cost response data set (W) of a grid as an influencing factor. The parameters were calculated using Application of Analytical Hierarchy Process (AHP)[59].

Creating the cost-response data set: The cost-response data include: geographical environment factors composed of Slope (A), river (R), altitude (H), relief (U), and accumulated temperature (T), socio-economic factors composed of vegetation net primary production (NPP) (N) and population density (P). The function of the resulting cost-response data set is as expressed in Eq (5).


Parameters in the calculation: Human road choice is related to accumulated, long-term experience. The AHP was used to make quantitative experience judgments (Table 1). The more accurate square root method was used to calculate the standard layer single factor weight, which was also used to test consistency. The result of the consistency test was CR = 0.05, CI = 0.05, with a smaller CI value indicating a better consistency. The overall CI and CR were less than 0.1 (detailed description in S1 Text.).

Tab. 1. Weight calculation result.
Weight calculation result.

Calculation of the length of the edge through raster: If: the coordinates of p1 and p2 are (x1y1) and (x2y2), and p1 and p2 constitute an equation of a straight line, then: the intersection coordinates of p1p2¯ and each raster boundary can be calculated as follows Eq (6): If: p1p2¯ has m intersection points with the rasters, and the intersection points are (p(1),p(2),⋯p(m)), Then: the coordinate of the intersection point i (p(i)) satisfy Eq (6).


Eq (7) shows the calculation of the distance between p1 and p2 actually passing through the raster (L(i) shown in Fig 3): the data set (p1p(1),⋯p(m)p2) was sorted by x-coordinate size (p(0),p(m)p(m+1)) and the Euclidean distance (D) between two adjacent points was calculated.


Where D(⋅,⋅)is the Euclidean distance and i is the intersection i after sorting.

Calculation of the length (<i>L</i><sup>(<i>i</i>)</sup>) of the edge through raster.
Fig. 3. Calculation of the length (L(i)) of the edge through raster.

Detailed description of data standardization

Table 2 provides the results of our data standardization grading value assignments. The cost-response data include: slope, river, altitude, relief, accumulated temperature, vegetation NPP, and population density. The detailed descriptions for data standardization for each are as follows:

  • Slope: Slope affects regional material flow and energy conversions, which are important factors impacting traffic lines [60]. First, this dataset considered the impact of the regional gradient on the road. Past studies have shown that humans cannot walk on a gradient of more than 70°; moreover, human activity has been shown to be mainly concentrated on terrain where the gradient is less than 40° [60]. The interval in the areas with a gradient of 0–40° was 5°, and the interval with gradient >40° was 10°. Values were assigned according to the following principle: The lower the slope, the lower the assigned value (Geospatial Data Cloud. Qinghai-Tibet Plateau DEM data product with 90m×90m spatial resolution. 2012, Database: [Internet]. Available from: http://www.gscloud.cn/)

  • River: Rivers are important supply routes. Within the study area there were a number of large rivers flowing from west to east, effectively blocking north-south traffic. The larger the river, the harder the crossing; comparatively, a small tributary river is more easily crossed and may also provide a source to replenish water supplies. The classification used was based on river order (1, 3, 4, and 5) and river buffer size (5 km, 7.5 km, and 10 km). The buffer zone design of rivers is that we think the ancient human walking passages mainly traveled on the river terraces, and did not sail against the water directly on the river. We divided the buffer zones into 5km, 7.5km and 10km for the convenience of water intake. Take level 1 rivers as an example. The terraces 5km away from the level 1 river are mainly V-shaped valleys with large slope, and water intake is not easy, therefore, the cost is relatively large. The terraces 7.5 km away from the river are relatively flat and easy to fetch water (In Level 1 Rivers, the flow is high; therefore an additional buffer value of 7.5km was added.). It is very difficult to get water 10km away from the river. The principle governing river assignment was as follows: the lower the river grade, the lower the assignment; similarly, the lower the river buffer from the river, the lower the assignment (National Earth System Science Data Sharing Platform. China 1:250000 first, third, fourth and fifth river grading data sets. 2002, Database: [Internet]. Available from: http://www.geodata.cn/)

  • Altitude: Medical research has shown that humans cannot survive at altitudes above 5,500 m for long periods of time; comparatively, when altitude is lower than 1,600 m, humans have no hypoxic response [61]. Past work has shown that when people enter the QTP for the first time, the incidence of acute high altitude diseases (AHAD) in the areas with altitudes ranging from 2,500–3,000 m is generally not more than 10%. When the altitude ranges from 3,000–4,000 m, incidence can exceed 30%, and from 4,000–5,000 m the incidence of AHAD is more than 50%. When the elevation is greater than 5,000 m, death can occur[62]. Therefore, elevation is a natural obstacle to human plateau traffic. In areas with altitude ≤ 4,000 m, the interval was 500 m, and in areas with altitude > 4,000 m, the interval was 250 m. These values were assigned according to the principle that the lower the altitude, the lower the value (Geospatial Data Cloud. Qinghai-Tibet Plateau DEM data product with 90m×90m spatial resolution. 2012, Database: [Internet]. Available from: http://www.gscloud.cn/).

  • Relief: Relief is an important index to quantitatively describe geomorphological form and can be used to divide the landscape into different geomorphologic types [63]. The level of fluctuation determines the role of the terrain in the choice of roads. This is because people will choose flatter plains for roads, to avoid as much as possible the high, undulating mountainous regions. Nine levels were used in this study to define relief: 0–30 m, 30–50 m, 50–70 m, 70–150 m, 150–300 m, 300–450 m, 450–600 m, 600–1000 m, 1000–2500 m, and >2500 m. The values were assigned according to the principle that the lower the relief, the lower the assigned value (Geospatial Data Cloud. Qinghai-Tibet Plateau DEM data product with 90m×90m spatial resolution. 2012, Database: [Internet]. Available from: http://www.gscloud.cn/).

  • ≥ 0°C accumulated temperature: ≥ 0°C accumulated temperature is an important index to sub-divide the QTP into climatic zones [64]. Travelers entering the plateau had a poor adaptive response to the colder climate and preferentially chose areas with higher temperatures. Therefore, the accumulated temperature used was from high to low: 1200°C, 1100°C, 700°C, 500°C, 300°C, and 200°C for interval classification; the higher the accumulated temperature, the higher the value (Database Available from: ≥0°C accumulated temperature is obtained by the Chinese Academy of Agricultural Sciences, 1981–1990; Project sources: Database of the temperature and humidity data level in the background layer of China's ecological environment completed by the Agricultural Regional Planning Institute Database: [Internet]. Available from: http://www.caas.cn).

  • Vegetation net primary production (NPP): The NPP of vegetation refers to the total amount of dry, organic matter produced by plants in a unit of time and per unit area [65]. Given this, NPP can be used as an indicator of production for different vegetation types. The traditional economic model in this area is animal husbandry, so grassland has become an important seasonal productivity measure for QTP. For ancient travelers who had to carry a large amount of supplies, the supply of animal feed in long journeys was extremely important. Notably, animal feed could be provided by plants. Given the methodological complexity of accurately rendering vegetation changes in spatial distribution of grasslands, vegetation NPP data were most generalizable for reflecting the productivity level of vegetation cohorts under natural environment conditions in our model. The NPP can be used to effectively reflect the geographic distribution of QTP pasture resources associated with the productivity level in the animal husbandry economic system. Furthermore, the satellite image archive is limited to vegetation indices (NDVI) mapped from more than 30 years of data from MODND1D or AVHRR, so the longer-term averages are not necessarily better for modeling scalar changes in vegetation geography over the past 2000 years or more [66]. Thus, the NPP value we chose was not fixed per time and per unit area, but the fixed value with general significance calculated by combining vegetation types. Although not representative of every year over the past 2000 or more years, this vegetation NPP was reclassified according to the size associated with the productivity level of vegetation cohorts, for use as a generic modeling dataset. Our classification of vegetation types by NPP value was based on the calculation results of Chen et al. (2012), Zhou et al. (2004), and He et al. (2004) [6769]. We predicted that the higher the NPP value, the higher the productivity level (Data Center of the Chinese Academy of Sciences by the Resource and Environmental Science. Spatial distribution data of 1 million vegetation types in China. 2012, Database: [Internet]. Available from: http://www.resdc.cn/data).

  • Population Density: Population density reflects the regional distribution of the population and is closely related to natural environment and level of socio-economic development [70]. Population density is not only an important indicator for gauging the level of human activities, but also an important factor in determining settlement location, which played a vital role in the development of ancient trade routes[71]. The areas with the highest population densities on the QTP are the Yellow River and the Huang Shui River Valley area in the eastern (Xining, 1000 people/km2) and the Lhasa River Valley area in the western (Lhasa, 1000 people/km2), which is consistent with the distribution pattern of population density in ancient times. The reason for this result is that the extreme environment limits the location of the settlement on the QTP, and the environment background in which human beings can settle down and live has relatively fixed[72]. Although the modern and ancient population densities cannot be exactly the same, the population density is a historical continuity, and the distribution of modern populations is a result of the past evolution of the regional population density. Therefore, modern population density data was used in this study. Our population density factor was classified based on various ranges of density rankings[73], with the higher population density giving a lower value. Briefly, more human traffic would be required in areas with higher population density, and thus these areas would be more possible to form roads. (Data cloud of Chinese academy of sciences. Population data of 1 km2 in China. 2000, Database: [Internet]. Available from: http://www.data.ac.cn/).

Tab. 2. Grading value assignments for all data standardization.
Grading value assignments for all data standardization.

Search for the shortest path

The following process was implemented using the Arc GIS10.5 version (ESRI Corporation, Redlands, CA, USA). Step 1: Conversion of cost data set to raster data. Output raster dataset was based on the results shown in Table 2. Step 2: Calculation of weight using raster calculator. Based on the results in shown in Table 1 and Eq (8), the step 2 adapts the spatial analyst ‘map algebra-raster calculator’ in ArcGIS. A larger cost value indicates higher cost, and vice versa. (Fig 4) Step 3: Calculation of weight of each edge in the TIN. (Fig 5)


Raster data set after reclassification and cost weight.
Fig. 4. Raster data set after reclassification and cost weight.
Weight result of each edge in the TIN.
Fig. 5. Weight result of each edge in the TIN.

Step 4: Search of the shortest route in the network. The Dijkstra algorithm was used to calculate the shortest path from a source node to all other nodes. It is a classic algorithm that produces the shortest path in increasing order of path length [74]. Since the Dijkstra algorithm was used to search for the shortest path between two specified nodes in an existing network, the optimal route was calculated by setting Xining (Vs) as the starting node and Lhasa (Ve) as the target node. The step 3 and step 4 were executed in Python using ArcGIS’s native Python geoprocessing tools. The full programming codes necessary for executing the steps are available as a unified Python script in S2 Text.

Step 5: Processing route. We originally considered the network for the modeling (macro scale) to be a triangular network of straight lines, but the actual route (micro scale) was not straight lines. The path through each raster (micro-scale) cannot be a straight line, and these routes will also pass through the grid according to the actual cost. So we treated the routes that go through each grid. In order to produce a more accurate route, in this study, the model’s output result (macro-scale) was processed as close as possible to match the actual route (micro-scale). On the premise of transforming the raster into a matrix, the order of the matrix was equal to the number of rows of the raster multiplied by the number of columns of the raster. If it was a route, the value of the raster was 1, and the value of the element in the corresponding matrix was also 1; if it was not a route, the raster value was 0, and the value of the element in the corresponding matrix was also 0. Under this premise, a matrix was constructed, starting from the first route (V1) to the end of (Vn) (n = 52) Eq (9). The route formed by the cumulative sum was the simulation result (R).The step 5 was designed using the R Programming Language (Alcatel-Lucent Bell Labs, USA) script as shown in S2 Text.


Where: 3070 is columns count of the raster, 1975 is rows count of the raster, and n is the route count.

Results

Simulation results of the Tang-Tibet Ancient Road

We used the comprehensive cost data as a weight to calculate the shortest path. Based on final results from our online search, the shortest path from Xining (starting node) to Lhasa (ending node) was calculated with a total weight of 50 nodes. Data from the collecting node set included 29 settlement nodes, 15 document nodes, 7 site points, and a total weight of 5391 (Table 3). The route was the one with the smallest sum of weights and the shortest distance connecting Xining to Lhasa with a total length of 1356 km. The maximum and minimum distances between these 50 nodes were 108 and 1 km, respectively, with an average distance of 24.2 km between these 50 nodes. We confirmed that 14% of these 50 nodes had a distance between 50 and 79 km, while 41% had a distance between 20 and 50 km, the remaining 45% had a distance between 0 and 20 km.

Tab. 3. Output simulation result.
Output simulation result.

Simulation results of the comprehensive cost

The calculation results of the weights showed that slope (25%) was the most important factor in determining the route of the South Asian corridor across the QTP. Other important factors in determining the corridor included: river (18%), NPP factors (18%), altitude (14%), population density (12%), relief (11%), and accumulated temperature (2%). The highest and lowest slopes through which the simulated route passed were 32.34° and 0.21°, respectively. The highest relief through which the simulated route passed was 551 m, with a minimum of 6 m; The maximum and minimum altitudes through which the simulated route passed were 5,400 m and 2,285 m, The highest and lowest accumulated temperatures through which the simulated route passed were 23,810°C and 2,096°C, respectively(Table 4).

Tab. 4. Statistical results of simulated route comprehensive cost value.
Statistical results of simulated route comprehensive cost value.

As shown in Table 5, 11% of the simulated route crossed 10 km buffer of grade 5 rivers, and 10% crossed 5 km buffer. This is in line with the convenience of travelling on the plateau, mainly relying on small rivers. With population density as a route choice determinant, 93% of the simulated route passed through an area with a population density of 10 people/km2. Population density is closely related to route choice. Temperate grass and meadows were the main vegetation determinants and accounted for more than 45% of total vegetation through which the simulated route passed. The vegetation type meets the needs of animals to obtain feed, and can supply the feeding needs along the route.

Tab. 5. Statistical results of simulated route comprehensive cost value.
Statistical results of simulated route comprehensive cost value.

As shown in Table 6, with slope as a determining factor, a simulated route with a length of 180 km and passing over an approximate 20° slope accounted for 13% of the total route. The remaining 856 km of the simulated route passed over terrain with slope ≤ 10° and accounted for 63% of the total simulated route. As shown in Table 6 and Fig 6, the simulated route passed through relief areas between 100–200 m, which accounted for 320 km or 23% of the total simulated route. The simulated route passed through relief areas < 100 m, which accounted for a total length of 526 km or 39% of the total simulated route. A larger length (970 km) of the simulated route passed through areas with altitudes >4,000 m, which accounted for 71% of the total simulated route. This is very consistent with the topography of the study area. The slope and relief in over cross the high mountains is the highest. The section with the lowest terrain cost is concentrated in Qinghai province section. (Description of accumulated temperature of simulated route is explained in the discussion section: Comparison between simulated route and literature records).The Pearson correlation coefficient was determined using the weight value of the simulated route and the route distance. The correlation coefficient R2 = 0.9782, which was significant (two-sided) at the 0.01 test level.

Relationship between district, distance, and comprehensive cost.
Fig. 6. Relationship between district, distance, and comprehensive cost.
Tab. 6. Relationship between comprehensive cost, distance, and through the area in our simulated route.
Relationship between comprehensive cost, distance, and through the area in our simulated route.

Discussion

Comparison between simulated route and literature records

According to the literature [5355], the length of the TTAR is known to be 2,240 km. Table 7 shows the comparative results of the simulated TTAR distance with that found in the literature. These comparisons were composed of six sections of distance, including Xining-Xinghai County (X-X), Xinghai-MaduoCounty (X-M), Maduo-Chengduo County (M-C), Chenduo-YushuCounty (C-Y), Yushu-Zaduo County (Y-Z), and ZaduoCounty-Lhasa (Z-L). The distances of the first two sections (X-X and X-M) were almost identical to those from the literature. Moreover, in the first two sections, the key station nodes (Daheba, Kuhaitan, and Huashixia) that were documented in the literature also appeared in our simulated routes (Fig 7). The upper Yellow River lies in the M-C section. According to the literature, there were two ferries (Fig 7A and 7B) that crossed the Yellow River. The simulated route crossed the Yellow River at ferry b, owing to the shorter distance and lower cost relative to ferry a. In the M-C section, the distance of our simulated route was 79 km less than that recorded in the literature. This may be due to the choice of different routes over the Bayan Har Mountains. Across the entire route, the simulated mountain crossing routes were always shorter than those obtained from the literature. This was because the simulated route followed the principle of shortest distance and lowest cost. This was also the reason why the distance of the simulated route crossing the Tanggula Mountains was much shorter than that recorded in the literature. The distance of the simulated route from ChengduotoYushu County (C-Y section) was 134 km less than that recorded in the literature. This was because the route documented in the literature passed through Yushu, which was an important economic and trade center in this region. However, the simulated route only passed through Anchong Township located in the north of Yushu, 120 km away from Yushu (show in Fig 7). Further, passing through Anchong Township is in line with the principle of the shortest distance and lowest cost.

Comparison between simulated routes and key points in the literature.
Fig. 7. Comparison between simulated routes and key points in the literature.
Tab. 7. Comparison of distance between simulated routes and literature documentation.
Comparison of distance between simulated routes and literature documentation.

In the Z-L section crossing over the Tanggula Mountains, the distance between the simulated and the literature routes were different. Comparing the key points between the route recorded in the literature and our simulated route, we found that the Qinghai section had 38 key nodes while the Tibetan section had five key nodes recorded in the literature. In the Qinghai section, 55.26% (21/38) of the key points documented in the literature appeared in our simulated routes. In the Tibetan section, the five nodes (Baixiong, Nierong, Naqu, Dangxiong, and Yangbajing) from the literature showed deviations from our simulated route(Fig 7).

When comparing the five variables in the region model where the five nodes were located with the region where the simulated routes were located, only cumulative temperature was significantly different in the Z-L section (Table 8). The accumulated temperature found in the literature route was 5723.69°C higher than that of the simulated route. The simulated route passed through 370 km of areas with accumulated temperatures between 23810–10000°C, which accounted for 28% of the total simulated route. The simulated route passed through 850 km of areas with accumulated temperatures between 4000–10000°C, which accounted for 63% of the total and went from Xinghai County to Lhasa. Previous studies have suggested that the relative oxygen concentration (ROC) in the near-ground air shows no obvious changes at different altitudes, but that atmospheric pressure and oxygen partial pressure significantly decrease with increasing altitude [7577]. Further research has also shown that aside from altitude, surface vegetation coverage and weather conditions may also have an impact on ROC. These previous results indicated [78] that altitude, vegetation coverage, and 500 h Pa-T accounted for a total of 65.5% of the total variance in ROC based on principal component analysis. Moreover, the individual variance interpretation rate for atmospheric temperature (500 h Pa-T) variance, vegetation coverage variance, and altitude variance was 33.1%, 28.5%, and 3.9%, respectively. Owing to the same altitude of the simulated route and the literature route, there was no significant difference in absolute oxygen concentration. However, since the accumulated temperature increased, ROC also increased.

Tab. 8. Comparison of Z-L cost factors between simulated results and literature documentation.
Comparison of Z-L cost factors between simulated results and literature documentation.

In Tibet area, the ability of humans to adapt to high altitude hypoxic conditions decreased with an average altitude of 4,419 m. Given this, it was reasonable to forego the route with shorter distance—but low oxygen concentration(simulation routes)—and choose the route with longer distance—but higher accumulated temperature and oxygen concentration(Baixiong to Nierong and Yangbajing to Lhasa). Since the simulated route focused primarily on distance and cost, it was 508 km shorter than the route recorded in the literature. As the route to Tibet was closely related to difference in human adaptations and needs, there may be multiple routes for crossing the Tanggula Mountains to Lhasa.Our simulation revealed one such route.

Based on our findings, the constraint factors on people's choice of roads in complex, natural environments are complex. These factors also demonstrated that the model used in our research needs to be modified in order to adapt to the simulation of ancient routes at high altitudes. The simulation of the Qinghai section showed that the model was applicable for altitudes ranging from 2,200 m to 4,400 m. In the areas with an altitude higher than 4,400 m, additional influential factors need to be considered.

The formation and evolution of the Tang-Tibet Ancient Road

A few studies have indicated that climatic changes during the Common Era were linked with the human history of China [7980]. It is generally assumed that TTAR was formed during the Tang Dynasty (618 AD). This period also coincided with the rise of the Tubo Dynasty (618 AD-842 AD) (TBD). High-resolution climate reconstruction sequences have shown that the climatic characteristics of this period were warm and humid[8182]. The response of the plateau vegetation to this climatic change was sensitive and rapid [83]. With an increase in temperature, the vegetation zones of the plateau expanded from the Southeast to the Northwest, from the edge to the hinterland, and from low to high altitudes. The increase in vegetation provided resources for nomadism and animal husbandry, which were the primary tribal economic activities in the heartland of the QTP[5152]. Of the main vegetation types on the QTP, the range of alpine grasslands and meadows expanded during warm periods. This trend extended to high-altitude areas, which created the necessary conditions for the formation of trade routes. When this warm period ended, the vegetation belt moved wholly southward and the high-altitude routes were gradually abandoned. Given this, it is clear that during the cold period, the main route into Tibet was formed in the southern part of the QTP. Previous work has shown that the Buddhist cultural transmission route formed in the southern part of Tibet was also the main route to Tibet [55]. The formation time of this road coincides with the introduction of Buddhism into China, corresponding to the cold period of the Western Jin Dynasty (266AD-316AD)(Fig 8).

Distributions of simulated route, Buddhist cultural communication route, and sites.
Fig. 8. Distributions of simulated route, Buddhist cultural communication route, and sites.

In addition, there was a correlation between climatic change, site spatial distribution, and site number. Corresponding to the cold period of the Tubo tribe (TBT, period before 618 AD) and the warm period of the TBD, the number of sites in the TBT period (261) was significantly less than the number of sites in the TBD period (512). The TBT sites were mainly distributed along the Southeast edge of the plateau, while the distribution of sites during the TBD expanded to the inner QTP (Fig 8). The reason for this trend in site expansion was that the warm and humid climate provided good climatic conditions for human activities. Therefore, the ancient road into Tibet gradually transitioned from the southeastern edge to the hinterland of the plateau. This trend was the result of plateau settlement and the expansion of human activities to the plateau hinterland given the warm climate of the medieval period.

Conclusions

Here, we constructed a weighted trade route network, with data on slope, river, altitude, relief, accumulated temperature, vegetation NPP, and population density. We then used the principle of minimum cost and the shortest path on the network to simulate the Tang-Tibet Ancient Road. We found that slope, river, and vegetation NPP were the most important geographical determinants of the ancient road. Choosing a gentle slope, and routes along rivers and river valleys, and providing vegetation for feeding also corroborated that these three, major controlling factors must be present when traveling on the plateau. When comparing the simulated route with that obtained from the literature, the Qinghai section had a higher degree of fit in terms of both distance and direction; however, the Tibetan section deviated from this close fit. This deviation was due to ROC becoming a limiting factor when choosing the Tibetan section of the road. This finding further illustrated that avoiding hypoxia in high-altitude areas (4,400 m) was a fundamental factor for ancient human road choice. This also suggested that oxygen content in high-altitude areas needs to be fully considered in future models. Moreover, we argue that warm and humid climates as well as human settlement migration to the hinterland of the QTP were the fundamental driving forces for the formation of the Tang-Tibet Ancient Road. This evolution formed a stable settlement node in the center of the QTP. These nodes were like an island chain, running through the route connecting Xining to Lhasa and making the Tang-Tibet Ancient Road an important route for crossing the inner plateau.

Supporting information

S1 Table [docx]
GPS data set of the key points in the literature.

S1 Text [docx]
Supplementary materials.

S2 Text [r]
Script programming codes.


Zdroje

1. Liu WD, Michael D, Gao BY. Discursive construction of the Belt and Road Initiative: From neoliberal-al to inclusive globalization. Prog Geogr. 2017; 36(11): 1321–1331.

2. Yan GW. A Textual Study of Traffic Maps in Tang Dynasty.1st ed. Shang Hai: Shanghai Classics Publishing House; 2007. pp. 421–585.

3. Su HY, Yong JC. The comprehensive conclusion obtained from the traffic investigation of the Silk Road in Qinghai. Silk Route. 2009; (6): 39–42.

4. Liu XR. The Silk Road in World History.1st ed. New York: Oxford university press; 2010. pp. 30–65.

5. Valerie H.The Silk Road: A New History. 1st ed.New York: Oxford university press; 2012. pp. 46–57.

6. Zhou J. TheSilk Road Traffic Routes (China section) Historical Geography Study. 1st ed.Nan Jing: Jiangsu People’s Press; 2012. pp. 7–280.

7. Yao T. D, Chen FH, Cui P, Ma YM, Xu BQ, Zhu LP, et al. From the Qinghai-Tibet Plateau to the Third and Pan-third poles. Chin Sci Bull. 2017; 32(09): 924–931.

8. Chen FH, An CB, Dong GH, Zhang DJ. Human Activities, Environmental Changes, and Rise and Decline of Silk Road Civilization in Pan-Third Pole Region. Chin Sci Bull.2017; 32(9): 967–975.

9. Han JY. Sino-Western exchange culture South Road in 5,000 years ago. Soc Sci Front. 2012; (06): 102–106.

10. Lv HL. Understanding the Neolithic Age in Western Tibet from the perspective of crossing the Himalayas. Archaeology. 2014;12:77–89.

11. Zhang DJ, Dong GH, Wang H, Ren XY, Ha BB, Qiang MR, et al. History and possible mechanisms of prehistoric human migration to the Tibetan Plateau. Sci China Earth Sci. 2016;59: 1765–1778.

12. Liu XY, Lister DL, Zhao ZJ, Staff RA, Jones PJ, Zhou LP, et al. The virtues of small grain size: Potential pathways to a distinguishing feature of Asian wheat. Quatern Int. 2016 Dec 28. pii: http://dx.doi.org/10.1016/j.quaint.2016.02.059.

13. Hou GL, Xu CJ, Cao GC, Er CY. The spatial-temporal simulation Of mankind's expansion on the Tibetan Plateau during last DE glaciation- middle Holocene. Quatern Sci. 2017; 37(4):709–720.

14. Zeng XQ, Guo Y, Xu QJ, Mascher M, Guo GG, Li SC, et al. Origin and evolution of qingke barley in Tibet. Nature. 2018 Dec 21.pii:https://doi.org/10.1038/s41467-018-07920-5.

15. Hou W. Viewing Prehistoric Traffic and Trade in Tibet from Archeological Discoveries. China Tibet ology. 2013; 107(02): 5–24.

16. Huo W. On the early metallic wares and the early metallic Age in Tibet. Acta ArchaeologicaSinica. 2014; (3):327–350.

17. Tong T, Li LH. Himalayan golden mask in the Eurasian view.Archaeology.2015; (02): 92–102.

18. Lu HY, Zhang JP, Yang YM, Yang XY, Xu BQ, Yang WZ, et al. Earliest tea as evidence for one branch of the Silk Road across the Tibetan Plateau.Sci Rep-UK. 2016; 6(1): 18955.

19. Zhou WZ. Ancient Road investigates in Qinghai Area. J N W Univ (Phil Soc Sci). 1982; (01): 65–72.

20. Cui YH. Changes in the development of the Qinghai part of the Silk Road. Qinghai Soc Sci. 2016; (01):9–16:31.

21. Jia X, Dong GH, Wang L, Ma MM, Lee HF., Zhang ZC, et al. How Humans inhabited the Northeastern Tibetan Plateau during the Little Ice Age: a case study at Hualong County, Qinghai Province, China. Journal of Archaeological Sci Reports 2016;7:27–36.

22. Wu Z. Research on Qinghai road. N W J Ethnol. 1992; (2): 123–140.

23. Hou W. Viewing Tubo and Ancient Transport in Central Asia and West Asia from Archeological Material- The role of western Tibet in the spread of Tibetan Buddhism is also discussed. China Tibet ology. 1995; (04): 48–63.

24. Sun XS. Four questions on China-Federal Democratic Republic of Nepal traffic in the early Tang Dynasty.China Tibet ology. 2000; (04):64–74.

25. Li ZJ. The Southern Route to Tibet along the Hexi Corridor at the Tang Dynasty. Dunhuang Res.2007; (3):44–49.

26. Su HY. On the Evolution of Qinghai Silk Road from the International Perspective. Qinghai J Ethnol.2012; (3):136–139.

27. Chen XP. The Tang-Tibet Ancient Road. Qinghai Soc Sci. 1987;(03):70–76.

28. Bi JT, Wang XX. Researches and application of geo-spatial Technology to the history of the Silk Road. Arid Land Geogr. 2007; 30(6): 954–961.

29. Yue L, Zhang P. A GIS-based study of the real conditions of Ling Zhou Road as a part of the Silk Roads in the early Northern Song dynasty. J Yunnan Univ (Soc Sci). 2017 (16):62.

30. MaJ Li FY, PangGW Li CR, Liu W. The Restoration of the Ancient Silk Road on Land and Analysis of Geographical Features along the Route. Geogr Geo-Inform Sci. 2017 (4):123–128.

31. Frachetti MD, Smith CE, Traub CM, Williams T. Nomadic ecology shaped the highland geography of Asia’s Silk Roads. Nature. 2017;543(7644):193–198. doi: 10.1038/nature21696 28277506

32. CiolkTM. Old Word Trade Routes (OWTRAD) Project. 2012 Jul15.pii: http://www.ciolek.com/owtrad.html.

33. WilliamsT. The Silk Roads: an ICOMOS thematic study. Icomos. 2014Sep 08.pii: http://openarchive.icomos.org/id/eprint/1487

34. Wu J. From the Jiezhi Bridge to XinuoluoDak of the Trend of the Tang-Tibet Ancient Road—A Discussion on Some Issues in Sato's "Research on the Historical Geography of Tibet". China Tibet ology. 1988; (02):116–123.

35. Qing YJ, Zhang HS. Analysis of the several stations and the Dafei location in the Tang-Tibet Ancient Road in Hainan Prefecture. Qinghai J Ethnol. 1994; (04): 57–63.

36. Su HY, Yong JC, Yan B, You XN, et al. A new test on the trend of the Tang-Tibet Ancient Road from Da Zhen Pass to Shan Cheng section. J Qinghai Natl Univ (Soc Sci). 2011;37(03): 62–67.

37. Luo Y. A brief talk on the formation of Ji Long Roads used by the government during Tubo Period. Theor Platform Tibetan Dev. 2015; (01): 69–74.

38. Yu XH, Xi L. Archeological discoveries and reconstruction of road network in the Tang-Tibet Ancient Road. J Xizang Minus Univ (PhilosSoc Sci Edit). 2017; 38(06):53–59+155.

39. Chinese Academy of Sciences Computer Network Information Center. Tibetan Plateau DEM data (90m×90m); 2000 [cited 2018 Nov 10].Database:Open source data [Internet]. Available from:http://www.gscloud.cn

40. National Earth System Science Data Sharing Platform.China 1:250,000 data set of grade 1, 3, 4 and 5 rivers; 2002[cited 2018 Nov 10].Database: Open source data [Internet]. Available from:http://www.geodata.cn

41. CAAS. Establishment of National Basic Resource and Environmental Remote Sensing Dynamic Information Service System (96-B02-01): Construction of Temperature and Humidity Data Level at Ecological Environment Background (96-B02-01-02); 1995[cited 2018 Nov 10]. Database: Open source data[Internet]. Available from: http://www.caas.cn

42. China by the Resource and Environmental Science Data Center of the Chinese Academy of Sciences. Spatial distribution data of 1:100,000 vegetation types in China; 2001[cited 2018 Nov 10].Database:Open source data [Internet]. Available from:http://www.resdc.cn/

43. Thematic Database for Human-Earth System. China population density data; 2000[cited 2018 Nov 10]. Database: Open source data[Internet]. Available from: http://www.data.ac.cn.

44. Zhang YL, Li BY, Zheng D. On the scope and area of the Qinghai-Tibet Plateau.GeogrRes. 2002; (01):1–8.

45. Ma YM, Hu ZY, Tian LD, Zhang F, Duan AM, Yang K, et al. Study progresses of the Tibet Plateau climate system change and mechanism of its impact on East Asia. Prog Geogr.2014; 29(2): 207–215.

46. ZhengD,ZhaoDS. Natural Environmental Characteristics of Qinghai—Tibet Plateau. Sci &Technol Rev. 2017;35(06): 13–22.

47. Wu RS, Ma YM. Comparative analyses on radiation characteristics in different areas over the Tibetan Plateau. PlateauMeteoro.2010; 29(2): 251–259.

48. Zheng D, Zhang RZ, Yang QY. On the natural zonation in the Qinghai-Tibetan Plateau. Acta Geogr Sin. 1979; 34(1): 1–11.

49. Zhao ZP, Liu JY, Shao QQ. Characteristic analysis of land cover change in nature reserve of Three River's source regions. Sci Geogr Sin. 2010; 30(3): 415–420.

50. National cultural heritage administration. Atlas of Chinese Cultural Relics: Tibet Volume and Qinghai Volume. 1st ed. Bei Jing: Cultural Relics Press; 2010. pp. 207–404(Tibet); 36-92(Qinghai).

51. Zhou XW. Yu Shu Survey. Xi Ning: Qinghai People’s Publishing House; 1986. pp. 22–106.

52. Ma HT. Gan Su-Qinghai-Tibetan Border Area Investigation. 1st ed. Lanzhou: Gan Su People’s Publishing House; 2003. pp. 137–276.

53. Chen XP. The Tang-Tibet Ancient Road. 1st ed. Xi’an: San Qin Publishing House; 1989. pp. 58–66.

54. Edited by the editorial department of Zhong Hua Book Company. Twenty-Four Books of History: [The Song Dynasty] New Tang Book • Geography • Long You Road • Xi Ping County to Shan Zhou.1st ed. Bei Jing: Zhong Hua Book Company Press; 2017. pp. 1040–1120.

55. Shaanxi Institute of Archaeology et al. From Chang'an to Lhasa: the Tang-Tibet Ancient Road Investigation Tour in 2014. Shang Hai: Shanghai Classics Publishing House; 2017. pp. 73–169.

56. National Earth System Science Data Sharing Platform. The geographical location data of administrative centers at the township level; 2012[cited 2018 Nov 10]. Database: Open source data [Internet]. Available from: http://www.tianditu.gov.cn/

57. Tsai Victor JD. Delaunay triangulations in TIN creation: an overview and a linear-time algorithm. Inter J GeogrInformSyst.1993; 7(6): 501–524.

58. Chen MJ. The difference between a mile and an acre. Archaeology. 1966; (01): 36–45.

59. Saaty ThomasL. Decision making—the Analytic Hierarchy and Network Processes (AHP/ANP). J Syst Sci Syst Eng. 2004;13(1):1–35.

60. Groucutt HS, White TS, Clark-Balzan L, Parton A, Crassard R, Shipton C, et al. Human occupation of the Arabian Empty Quarter during MIS 5: Evidence from Mundafan Al-Buhayrah, Saudi Arabia. Quaternary Sci Rev.2015; 119: 116–135.

61. Wu TY. Challenges to Human Being under Hypoxia Environment on Plateau. J Med Res. 2006; 35(10): 1–3.

62. Xue HJ, Fang CJ. Theory of human adaptation to plateau and its application to traveling to the Tibet. J Arid LandResourEnviron.2011; 25(02): 183–189.

63. Feng ZM, Tang Y, Yang YZ, Zhang Dan. The ReliefDegree of Land Surface in China and Its Correlation with Population Distribution. Acta Geogr Sin. 2007; 62(10): 1073–1082.

64. Dai JX. The Climate of Qinghai-Tibet Plateau. Bei Jing: China Meteorological Press; 1990. pp.4–5.

65. Zhang YL, Qi W, Zhou CP, Ding MJ, Liu LS, Gao JG, et al. Spatial and temporal variability in the net primary production of alpine grassland on the Tibetan Plateau since 1982. Acta Geogr Sin.2013; 68(09): 1197–1211.

66. Frachetti MD, Smith CE, Traub CM, Williams T. Nomadic ecology shaped the highland geography of Asia’s Silk Roads(Supplementary Discussion). Nature. 2017; 543(7644):8–11

67. Chen ZQ, Shao QQ, Liu JY; Wang JB. Analysis of net primary productivity of terrestrial vegetation on the Qinghai-Tibetan Plateau based on MODIS remote sensing data. Sci Sin (Terrae). 2012;42(03): 402–410.

68. Zhou CP, Ou YH, Wang QX, Du BZX, Sun QQ. Estimation of Net Primary Productivity in Tibetan Plateau ActaGeogr Sin. 2004; (01): 74–79.

69. He JS,FangJY, Piao SL, Xiao Y. Spatial Distribution of Grassland Biomass in China. Chinese J Plant Ecol.2004; (04):491–498.

70. Lu DD, Wang Z, Feng ZM. Academic debates on Hu Huan Yong population line. Geogr Res. 2016; 35(5): 805–824.

71. Jia X, Lee HF, Cui MC, Cheng GQ, Zhao Y, Ding H, et al. Differentiations of geographic distribution and subsistence strategiesbetween Tibetan and other major ethnic groups are determined by the physical environment in Hehuang Valley. Sci China Earth Sci. 2019;62:412–422.

72. Jia X, Lee HF, Cui MC, Liu C, Zeng L, Yue RPH, et al.Habitat Variability and Ethnic Diversity in Northern Tibetan Plateau. Sci Reports. 2017;7: 918.

73. Ge ML, Feng ZM. GIS——based Analysis of Population Distribution in 2000 in China. Popul Res. 2008; (32): 51–57.

74. Zhang FH, Liu JP, Li QY. A New Way of Network Analysis Based on Dijkstra. Remote Sens Inform. 2004; (02): 38–41.

75. West JB. Prediction of barometric pressures at high altitudes with the use of model atmospheres. J Appl Physiol. 1996; 81(4): 1850–1854. doi: 10.1152/jappl.1996.81.4.1850 8904608

76. Zha RB, Sun GN, Dong ZB, Yu ZK. Assessment of Atmospheric Oxygen Practical Pressure and Plateau Reaction of Tourists in the Qinghai-Tibet Plateau.Ecol and Environ Sci.2016; (1): 92–98.

77. Beall CM. Two routes to functional adaptation: Tibetan and Andean high-altitude natives. P Natl Acad Sci USA. 2007; 104 (Supplement 1): 8655–8600.

78. Shi PJ, Chen YQ, Zhang AY, He Y, Gao MN, Yang J, et al. Factors contribution to oxygen concentration in Qinghai-Tibetan Plateau (in Chinese). Chin Sci Bull. 2018Nov 23. Pii:http://kns.cnki.net/kcms/detail/11.1784.N.20181122.1426.010.html

79. Zhang ZB, Tian HD, Cazelles B, Kausrud KL, Bräuning A, Guo F, et al. Periodic climate cooling enhanced natural disasters and wars in China during AD 10–1900. P Roy Soc B-Biol Sci. 2010; 277(1701): 3745–3753.

80. Feng XP,ZhaoC, D’Andrea WJ, Liang J Zhou AF, Shen J. Temperature fluctuations during the Common Era in subtropical southwestern China inferred from brGDGTs in a remote alpine lake. Earth Planet Sc Lett. 2019; 510: 26–36.

81. Ge QS, Zheng JY, Hao ZX. PAGES synthesis study on climate changes in Asia over the last 2000 years: Progresses and perspectives. Acta Geogr Sin. 2015; 70(03): 355–363.

82. Yang B, Qin C, Wang J, He M, Melvin TM, Osborn TJ, Briffa KR. A 3500-year tree-ring record of annual precipitation on the north-eastern Tibetan Plateau. P Natl Acad Sci USA. 2014;111:2903–2908.

83. Hou GL, Yang P, Cao GC, E CY, Wang QB. Vegetation evolution and human expansion on the Qinghai–Tibet Plateau since the Last DE glaciation. Quatern Int.2017 Feb12. pii: https://doi.org/10.1016/j.quaint.2015.03.035


Článek vyšel v časopise

PLOS One


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#