Human paleodiet and animal utilization strategies during the Bronze Age in northwest Yunnan Province, southwest China

Abstract

Reconstructing ancient diets and the use of animals and plants augment our understanding of how humans adapted to different environments. Yunnan Province in southwest China is ecologically and environmentally diverse. During the Neolithic and Bronze Age periods, this region was occupied by a variety of local culture groups with diverse subsistence systems and material culture. In this paper, we obtained carbon (δ¹³C) and nitrogen (δ¹⁵N) isotopic ratios from human and faunal remains in order to reconstruct human paleodiets and strategies for animal exploitation at the Bronze Age site of Shilinggang (ca. 2500 Cal BP) in northwest Yunnan Province. The δ¹³C results for human samples from Shilinggang demonstrate that people’s diets were mainly dominated by C₃-based foodstuffs, probably due to both direct consumption of C₃ food and as a result of C₃ foddering of consumed animals. Auxiliary C₄ food signals can also be detected. High δ¹⁵N values indicate that meat was an important component of the diet. Analysis of faunal samples indicates that people primarily fed pigs and dogs with human food waste, while sheep/goats and cattle were foddered with other food sources. We compare stable isotope and archaeobotanical data from Shilinggang with data from other Bronze Age sites in Yunnan to explore potential regional variation in subsistence strategies. Our work suggests that people adopted different animal utilization and subsistence strategies in different parts of Yunnan during the Bronze Age period, probably as local adaptations to the highly diversified and isolated environments in the region.

About the Authors:

Lele Ren

Affiliation: MOE Key Laboratory of Western China's Environmental System, College of Earth Environmental Sciences, Lanzhou University, Lanzhou, Gansu Province, China

Xin Li

Affiliation: School of Public Health, Lanzhou University, Lanzhou, Gansu Province, China

Lihong Kang

Affiliation: Yunnan Provincial Institute of Cultural Relics and Archaeology, Kunming, China

Katherine Brunson

Affiliation: Joukowsky Institute for Archaeology and the Ancient World, Brown University, Providence, Rhode Island, United States of America

Honggao Liu

Affiliations MOE Key Laboratory of Western China's Environmental System, College of Earth Environmental Sciences, Lanzhou University, Lanzhou, Gansu Province, China, College of Agronomy and Biotechnology, Yunnan Agricultural University, Kunming, China

Weimiao Dong

Affiliation: Department of Cultural Heritage and Museology, Fudan University, Shanghai, China

Haiming Li

Affiliation: MOE Key Laboratory of Western China's Environmental System, College of Earth Environmental Sciences, Lanzhou University, Lanzhou, Gansu Province, China

Rui Min

Affiliation: Yunnan Provincial Institute of Cultural Relics and Archaeology, Kunming, China

Xu Liu

Affiliation: Yunnan Provincial Institute of Cultural Relics and Archaeology, Kunming, China

Guanghui Dong

* E-mail: [email protected]

Affiliation: MOE Key Laboratory of Western China's Environmental System, College of Earth Environmental Sciences, Lanzhou University, Lanzhou, Gansu Province, China

ORCID http://orcid.org/0000-0001-7796-9600

Introduction

Analyzing ancient dietary signatures and the utilization of animal and plant resources can reveal how humans choose appropriate subsistence strategies in different cultural contexts, or under different environmental constraints [e.g. 1–9]. Archaeobotanical, zooarchaeological, and stable isotopic analyses of floral and faunal remains unearthed from prehistoric sites in China provide valuable datasets for exploring these issues. For example, recent research reveals that people adopted different subsistence strategies in order to adapt to environmental variations between China’s Loess Plateau and the high altitude Tibetan Plateau [1–2, 10], and between northern and southern China [6, 11–12]. The roles of plants and animals in ancient subsistence strategies have been studied intensively in recent years [e.g. 1, 13–16], especially through stable carbon and nitrogen isotopic analysis of human and animal bone collagen [e.g. 12, 17–24]. Isotopic methods have been widely used at prehistoric archaeological sites in the Yellow River valley [e.g. 17–20], where there is a significant distribution of Neolithic and Bronze Age sites with good bone preservation. However, in south China, especially in the Yunnan-Guizhou Plateau of southwest China, the application of stable isotope analysis in archaeological research remains limited. There is great potential to use isotopic techniques to expand our understanding of ancient subsistence systems in this environmentally and culturally diverse region.

In contrast to the Yellow River and Yangtze River valleys, where Neolithic cultures were present since the early Holocene, the Neolithic in Yunnan did not begin until around 5000 BP [25]. Neolithic and subsequent Bronze Age cultures in Yunnan Province are characterized by a good deal of variation in cultural features between micro-regions [25, 26]. The better studied Bronze Age cultures (ca. 3100–2000 BP) can be divided into four main groups, including the Dian culture type centered in the Dianchi Lake area (near Kunming), the Erhai culture type located around Erhai Lake, the Dianbian culture type in the mountainous western frontier of Yunnan, and the Honghe culture type in the Honghe (Red River) valley [26]. Although archaeologists assume that subsistence strategies were a crucial factor in the differentiation between these culture types [25], the similarities and differences in how groups exploited natural resources in various parts of Yunnan remain enigmatic.

Previous archaeobotanical studies have provided insight into plant utilization strategies in Yunnan and neighboring Guizhou Provinces during the late Neolithic and Bronze Age periods [e.g. 27–29]. The identification of plant remains from domestic crops indicates that there were three distinct phases of agricultural production in the region: from 4800–3900 BP when rice was first cultivated; from 3900–3400 BP when rice was planted together with foxtail millet and broomcorn millet; and from 3400–2300 BP when there was mixed rice, millet, wheat, and barley agriculture [28]. This last phase corresponds to the Bronze Age period. However, we do not know how important these various domestic crops were in ancient human diets. Previous zooarchaeological studies indicate that people utilized a variety of livestock including dogs, pigs, sheep, goats, and cattle during the Bronze Age in Yunnan [e.g. 30–31]. How people foddered and raised these domestic animals is not well understood.

Carbon and nitrogen stable isotope analysis provides a method for reconstructing the importance of various domestic crops and wild resources in human and animal diets and for exploring how humans utilized plant and animal resources in antiquity [32–34]. For Yunnan, data on the carbon and nitrogen isotopic values of human bone collagen have been reported for only two Bronze Age sites [24, 35]. In this paper, we present isotopic data from both human and animal bones excavated from a third site, Shilinggang, which greatly increases our understanding of Bronze Age subsistence systems in the region. Shilinggang is located in the middle course of the Nujiang River in northwest Yunnan. It is classified as belonging to the Dianbian culture type. Shilinggang was the first site to be excavated in the Nujiang River Valley, which makes it a critical location for understanding human subsistence strategies in this sub-region and for expanding our understanding of the Dianbian culture. Multi-disciplinary research involving archaeologists, archaeobotanists, zooarchaeologists, and physical anthropologists has been carried out at Shilinggang since excavations began in 2013, resulting in several preliminary publications [28, 36, 37]. Here we report a new dataset from the stable carbon and nitrogen isotope analysis of 16 human bones and 58 animal bones unearthed from Shilinggang. We compare these data with the results of archaeobotanical and stable isotopic analyses from the two other published Bronze Age sites in Yunnan in order to determine the nature of human and animal paleodiets and to analyze how people raised domesticated animals in Bronze Age Yunnan. Our results indicate that the human diet at Shilinggang relied mainly on C₃ food, and that people used a diverse set of strategies to raise livestock.

The study site

Shilinggang (N 25°38′57″, E 98°53′16″) is situated on a hill in Lushui County, Nujiang Lisu Autonomous Prefecture, northwest Yunnan, to the west of the Nujiang River (Fig 1). The site has a low-latitude plateau monsoonal climate, with an altitude of 842 m above sea level, a mean annual temperature of 20.2°C, and a mean annual precipitation level of 977.2 mm [38]. The region has high biodiversity due to its location in a low latitude area with abundant heat and light. In addition, there are numerous mountains and rivers in the area, leading to distinctively complicated topographic features that form many microclimates with intrazonal vegetation. Modern vegetation types and animal communities tend to vary with elevation [39]. Today, staple crops in the region consist of rice, maize, buckwheat, wheat, roots, and tubers, while meat is provided by livestock and occasional hunting [40]. In the past, abundant wild and domestic floral and faunal resources would have provided a wide range of options to fulfill human dietary needs.

[Figure omitted. See PDF.]

Fig 1. The location of Shilinggang and other sites mentioned in the text.

DEM (digital elevation model) data was downloaded from Geospatial Data Cloud (http://www.gscloud.cn/), and map features in the figure were modified from Li et al. [28].

https://doi.org/10.1371/journal.pone.0177867.g001

The site of Shilinggang was first discovered and investigated in 2003 by the Nujiang Prefectural Administration of Cultural Relics. It was further investigated in 2007 by the Yunnan Provincial Institute of Cultural Relics and Archaeology and the Nujiang Prefectural Administration of Cultural Relics. The first formal excavations of the site, which were also the first archaeological excavations in the Nujiang River Valley, were conducted in 2013 and 2014 by the Yunnan Provincial Institute of Cultural Relics and Archaeology. The site has a higher elevation in the north than the south. It covers a total area of nearly 100,000 m² with 2 m deep accumulations of cultural materials. The northern part of the site has thicker cultural layers that were excavated according to rectangular excavation units (Fig 2a). All of the materials analyzed in this study were excavated from the northern part of the site. Additional excavations in the central and southern parts of the site were conducted by excavating long narrow trenches. A total area of 500 m² has been excavated to date.

[Figure omitted. See PDF.]

Fig 2. Plan map of northern Shilinggang and stratigraphic column from unit TN6W3.

a) Plan map of the northern part of Shilinggang. Rectangular excavation units are labeled starting with the letter T. Features include burials (labeled with the letter “M”), pits (labeled with the letter “H”), and building foundations (labeled with the letter “F”). The features shown in Fig 2a are modified from Li et al. [28]. b) Photograph of stratigraphic profile in unit TN6W3 (center right of plan map in Fig 2a). Circles around numbers indicate stratigraphic levels.

https://doi.org/10.1371/journal.pone.0177867.g002

Stratigraphic divisions separate the site into seven cultural layers from top to bottom (Fig 2b). In total, excavators uncovered 42 graves, 4 ash pits, 2 activity floors, 2 house foundations, and the remains of burnt soil. Excavated artifacts include pottery, stone tools, bronze vessels, plant remains, and large quantities of human and animal bones. Graves account for the overwhelming majority of excavated features. Skeletons in these graves are relatively complete even though some bones are poorly preserved. There are no obvious traces of groundwater flow or other transportation of bones out of the graves. We believe that remains in grave features were not significantly disturbed through taphonomic processes.

All of the tombs are rectangular with a vertical earthen shaft. Most graves are single burials while a few are joint burials containing both adults and children buried together. In terms of size, graves are relatively small (110–260 cm long, 36–68 cm wide and 10–48 cm deep) and the vast majority contain no burial goods. Only about 10% of graves contain one or two bronzes, stone tools, or ceramic objects. There is no obvious evidence from mortuary remains to indicate the existence of strict social class differences. While house foundations and ash pits are few in number, the distribution of graves among the houses and other features suggests that there were no spatial divisions between areas of the site used for artifact production, living, burial, and other activities.

Artifacts excavated from all seven cultural layers at Shilinggang are consistent with the material culture of the Spring and Autumn and Warring States periods (770–221 BC) [36]. Two charred rice seeds from the fourth cultural layer of square TN6W2 and fifth cultural layer of square TN5W3 were directly radiocarbon dated to between 2723 Cal BP and 2339 Cal BP at 2σ, confirming that the site dates to the late Bronze Age [28]. Flat land in this area of Yunnan is scarce, and there are many large canyons that cut through the landscape. It is possible that the gentle slope of Shilinggang was a preferred place to live and grow crops, which may explain the site’s deep accumulation of Bronze Age cultural materials [36].

Material and method

Material selection

A total of 16 human bone samples (all from graves in the northern part of Shilinggang, including 2 child samples and 14 adult samples) and 58 animal bone samples (all from cultural layers in northern Shilinggang, including one sheep (Ovis aries), 4 goats (Capra hircus), 5 samples that could be either sheep or goats (Ovis aries/Capra hircus), 10 cattle (Bos sp.), 4 dogs (Canis familiaris), 13 pigs (Sus sp.), 2 Old World monkeys (Cercopithecidae), 12 muntjacs (Muntiacus sp.), and 7 unidentified deer (Cervidae)) were selected for stable isotope analysis. Bones from human graves were assigned to age and sex, if possible, based on standard morphological traits [41], by Dongyue Zhao of Jilin University. Animal bones were identified by Lele Ren of Lanzhou University. All samples date to the Spring and Autumn and Warring States periods. All human and animal bones are currently stored in Yunnan Provincial Institute of Cultural Relics and Archaeology, Kunming. Tables 1 and 2 list each sample including information on species, skeletal element, sample code, sample excavation context, and bone collagen quality control indicators.

[Figure omitted. See PDF.]

Table 1. Isotopic composition and quality indicators of animal samples from Shilinggang, Yunnan.

Samples marked in bold italics were found to be contaminated, and were not included in further statistical analysis. Context locations are shown in Fig 2a. Among the sample context codes, TNnWn (where n is an Arabic numeral) refers to excavation unit; circles with a number inside refer to the stratigraphic layer.

https://doi.org/10.1371/journal.pone.0177867.t001

[Figure omitted. See PDF.]

Table 2. Isotopic composition and quality indicators of human samples from Shilinggang, Yunnan.

The samples marked in bold italics were found to be contaminated, and were not included in further statistical analysis. Context locations are shown in Fig 2a. Among the sample context codes, TNnWn (where n is an Arabic numeral) refers to excavation unit; Mn (where n is an Arabic numeral) refers to the grave number. “Ind.” stands for indeterminate sex.

https://doi.org/10.1371/journal.pone.0177867.t002

Collagen extraction

Bone collagen was extracted from the human and animal samples following the protocol outlined by Jay and Richards [42], with some modification described in Ambrose et al. [43]. For all human and animal bones, a dense bone fragment of approximately 3 g was removed from each whole bone sample. We used a large 3 g sample size, which is more than the usual 200 mg used in similar studies, in order to ensure that enough bone collagen could be extracted from each sample since we predicted that the acidic soil at Shilinggang would lead to poor collagen preservation. An electric grinder was used to grind the outer and inner surface of the bone to remove surface contamination before further processing. After cleaning, the bone fragments were demineralized by soaking in 0.5 M HCl at 4°C, refreshing the solution every two days, until the bone became soft and no bubbles were emitted. The residue was washed repeatedly with deionized water to neutral pH, and then rinsed in 0.125 M NaOH for 20 h at 4°C to remove soil humic acids. Samples were washed again with deionized water to neutral pH, dissolved in 0.001 M HCl at 70°C, and put in a dry oven for two days to gelatinize. The gelatinized solution was filtered when hot and then freeze-dried using a Labconco* FreeZone Lypholizer for 48 h to retain the collagen. The final collagen was saved and weighed for later analysis.

Measurement of organic carbon and organic nitrogen content and stable isotope ratios

The organic C and N content and the stable isotope ratios of the collagen samples were measured using a multi-flow isotope-rationing mass-spectrometer (Isoprime 100 IRMS) combined with an elemental analyzer (Vario PYRO cube) in the Archaeological Stable Isotope Lab, Department of Scientific History and Archaeometry, University of the Chinese Academy of Sciences, Beijing, China. Before the measurement, sulfanilamide was added as a reference standard to calibrate the C and N content. Caffeine (IAEA-600, δ¹³C: -27.8±0.0‰, δ¹⁵N: 1.0±0.2‰), sucrose (IAEA-cH-6, δ¹³C: -10.5±0.1‰) and (NH₄)₂SO₄ (IAEA-N-2, δ¹⁵N: 20.3±0.2‰) were used to normalize both N₂ (AIR as standard) and CO₂ (VPDB as standard). During the testing process, a sulfanilamide sample and a laboratory standard sample (CAAS, δ¹³C: -14.7±0.1‰, δ¹⁵N: 7.0±0.1‰) were inserted into the collagen list every ten samples for real-time monitoring and calibration. The analytical precision on the results of samples repeated in triplicate was better than 0.2‰ for both carbon and nitrogen isotopic ratios.

Data analysis

All statistical analyses were carried out with IBM SPSS Statistics 19 for Windows. Hierarchical cluster analysis was used to analyze differences in collagen δ¹³C and δ¹⁵N values among samples (S1 Text, S1 Fig). The resulting dendrograms (constructed using average linkage between groups) were used to separate samples into groups with significantly different isotopic signatures that reflect groups of individuals with different diets or foddering strategies (S2 and S3 Figs). Hierarchical cluster analysis is a preferred statistical method over other tests such as K-means because it auto-generates clusters of similar data points without requiring the analyst to pre-establish the number of output groups. We also tested the dietary groups identified through hierarchical cluster analysis using the K-means method and found that these groups are apparent in either method. A One-sample Kolmogorov-Smirnov Test was used to test whether the data was normally distributed, and if they were not, a non-parametric Mann-Whitney U test was used. The significance level was set at p < 0.05.

Results and discussion

Bone collagen preservation and contamination

The chemical composition and biological characteristics of bone are influenced and changed by temperature, pH levels, humidity, and microorganisms in the burial environment, processes which together are referred to as bone diagenesis [44]. A prerequisite of dietary reconstruction is that bone must retain its original chemical and isotopic composition even after a long period of deposition. In Tables 1 and 2, we include three reliable measurements of collagen preservation quality, including the collagen yield (% yield), the content of C (% C) and N (% N), and the atomic ratio of C/N. Out of the total 86 samples, 74 samples with collagen representing more than 1% of the original bone mass were used for analysis (mean collagen yield = 4.1%, SD = 3.1). The other 12 samples with less than 1% collagen yield were discarded because it is impossible to eliminate potential contaminants in samples with very low collagen yield [45–46]. In addition, the percentages of C (mean = 42.8±0.2%) and N (mean = 15.1±0.1%) in the 74 samples included in the analysis are in close proximity to those of modern bone collagen (C: 41%, N: 15%) [45], indicating that the ancient collagen is well preserved. The atomic C/N of these samples (mean = 3.3), ranged from 3.2 to 3.6, which falls within the 2.9–3.6 range considered to be a quality indicator for excellent collagen preservation [47]. The 16 human and 58 animal samples which meet the above criteria can be considered well preserved for stable isotope analysis.

Animal stable isotope data and paleodietary analysis at Shilinggang

The isotopic composition of ingested food is reflected in body tissues [48]. Stable isotope ratios in bone collagen and apatite generally indicate the isotopic composition of an individual’s diet over a period of years [49, 50]. The average carbon stable isotopic value (δ¹³C) of an organism can be related to C₃ vs. C₄ ecosystem dominance [51] or marine vs. terrestrial food webs [52]. C₃-plants have more negative δ¹³C values (-27.1±2.0‰) compared to C₄-plants (-13.1±1.2‰) because of different mechanisms of CO₂ fixation during photosynthesis [53]. In China, foxtail millet, a common C₄ crop found at archaeological sites, has a much less negative δ¹³C value compared to C₃ crops such as rice or other wild vegetation [54, 55]. Macrobotanical remains and phytoliths of both foxtail millet and rice have been identified at Shilinggang [28]. Trophic level also influences stable isotope ratios. Bone collagen δ¹³C values are enriched about 0–2‰ between prey animals and their predators [32]. A significant enrichment in δ¹⁵N also occurs between an organism’s diet and its body tissues, leading to δ¹⁵N values 3–5‰ higher in the body than in the average diet [56], or up to 6‰ in some cases [57]. However, nitrogen isotopic ratios are affected by many other factors as well, such as manuring, breastfeeding, aridity, salinity, and nutritional stress [58–63].

The δ¹³C and δ¹⁵N values from the Shilinggang bone samples (Tables 1 and 2) are presented in the scatter diagram in Fig 3. Stable isotope analysis of the animal bones provides a baseline for understanding the local food web [64]. For convenience in further discussion, the animal data were classified into the following groups: wild herbivores (muntjac and deer); omnivores (dogs and pigs); domesticated herbivores (sheep/goats and cattle); and primates (monkeys).

[Figure omitted. See PDF.]

Fig 3. Scatter diagram of human and animal collagen carbon and nitrogen values from Shilinggang.

Analytical error is so small that it is contained within the symbols.

https://doi.org/10.1371/journal.pone.0177867.g003

The δ¹³C values of wild herbivores (n = 19, 12 muntjacs and 7 deer) range from -25.9‰ to -19.9‰ (mean = -22.5±1.4‰) and the δ¹⁵N values range from 4.1‰ to 7.0‰ (mean = 5.2±1.0‰), indicating a C₃ plant-based herbivorous diet. The narrow range of δ¹³C values suggests that these wild herbivores ate similar types of plants, and is consistent with the natural ecological background of C₃ plant-dominated vegetation in the region [39].

For the omnivore group (n = 17, 4 dogs and 13 pigs), all dogs were identified as fully domesticated in accordance with their bone morphological characteristics while pigs were identified as either domestic or wild. The δ¹³C values of dogs range between -20.5‰ and -14.8‰ (mean = -18.3±2.4‰). It appears that three of the dogs were fed predominantly with C₃ food, which possibly includes rice, root and tuber crops, and/or animal protein based on C₃ food. The one other dog with the highest δ¹³C value (-14.8‰) consumed a mixed C₃ and C₄ diet, suggesting that its diet likely contained millet and/or animal protein from animals that ate millet. The δ¹⁵N values of the dog bones have a narrow range from 7.4‰ to 8.4‰ (mean = 7.8±0.5‰). The offset of nitrogen stable isotope between wild herbivores and dogs is 2.6‰, which approaches the range of the N isotopic enrichment (3‰-5‰) of one trophic level [32, 65], indicating that Shilinggang dogs consumed animal protein in quantities possibly through association with humans [17]. Dogs likely lived in close proximity to their human owners. These dietary signatures may reflect intentional provisioning of dogs, but could also indicate that dogs ate human waste and food scraps available around the site.

The 13 pig samples can be divided into three groups (labeled as pig A, pig B and pig C, respectively in Fig 4) based on differences in δ¹³C and δ¹⁵N values revealed through Hierarchical Cluster analysis (S2 Fig). The pig A group includes 5 samples with the lowest δ¹³C values (mean = -21.2±0.9‰) and δ¹⁵N values (mean = 4.8±0.8‰). These values are similar to those of the wild herbivores. These individuals likely represent wild boars consuming C₃-based plants and little animal protein. The morphological characteristics of the bones that provided samples in the pig A group are also more similar to wild boar than domestic pigs. However, we can’t rule out the possibility that individuals in the pig A group could be domesticated, but free-range. Three pig samples with the highest δ¹³C values (mean = -17.3±1.3‰) and middle-range δ¹⁵N values (mean = 6.5±0.2‰) fall into the pig B group. These animals had a mixed C₃ and C₄ diet. C₄ plant species comprise a minority of terrestrial plants (less than 4%) and are maladapted to most shady temperate environments [66], such as the environment of western Yunnan. Pollen cores from Lake Erhai indicate that the most common plants in Yunnan during the late Bronze Age were C₃ plants, especially Pinus, and deciduous arboreal trees [67]. Although the distribution of wild C₄ plants in China is not well known, we assume that most of the natural flora in the region has low δ¹³C values. Given that the only known C₄ plant at Shilinggang is foxtail millet, we argue that the C₄ signal in the pig B group represents the consumption of domestic millet, its by-products, and/or animals that fed on millet.

[Figure omitted. See PDF.]

Fig 4. Isotopic value distributions for domesticated animals showing groups with different husbandry strategies.

https://doi.org/10.1371/journal.pone.0177867.g004

The remaining five pig samples correspond to the pig C group, with a middle-range mean δ¹³C value of -20.1±0.4‰ and the highest mean δ¹⁵N value of 8.6±0.3‰. These animals appear to have been at a similar trophic level to dogs and humans (Fig 3). It is likely that the pigs in group C were domesticated and, like most dogs, scavenged on mainly C₃ food scraps or were provisioned by people with C₃ foods. The diet also included plenty of animal protein. Stable carbon and nitrogen isotope analysis has been applied to track pig domestication at many other archaeological sites in China [18, 68, 69]. Our results indicate that δ¹³C values are also useful for distinguishing between remains of domesticated pigs (pig groups B and C) and wild boars (pig group A) at Shilinggang.

Hierarchical Cluster Analysis revealed that domesticated herbivores (10 sheep/goats and 10 cattle) could also be divided into three groups according to their δ¹³C values (S3 Fig and Fig 4). Based on previous isotopic studies of ancient Chinese fauna [18–19, 55], these groups represent three distinct dietary characteristics: group A represents individuals with diets of mainly C₄ plants (> -12‰), group B represents individuals with diets of mixed C₃ and C₄ plants (-18‰ to -12‰), and group C represents individuals with diets of predominantly C₃ plants (< -18‰). Group A is composed of two sheep/goats samples and one cattle sample with the highest δ¹³C values (mean = -10.4±0.8‰) and lowest δ¹⁵N values (mean = 3.4±0.3‰). These animals largely fed on C₄ plants, very likely millet by-products. One sheep/goats sample and six cattle samples in group B are identified as mixed C₃ and C₄ eaters, with mean δ¹³C and δ¹⁵N values of -15.5±1.9‰ and 4.5±0.6‰, respectively. The δ¹³C values likely represent a mixed C₄ signal from consumption of millet and a C₃ signal from consumption of rice and/or wild forage [28, 67]. Group C contains seven sheep/goats samples and three cattle samples with a mean δ¹³C value of -22.5±0.7‰ and a mean δ¹⁵N value of 5.4±1.2‰. This group consumed predominantly C₃ plants, likely rice and/or the natural vegetation surrounding the site [28, 67]. The mean of both δ¹³C values (Mann-Whitney test: U = 91.500, Z = -0.161, P = 0.872) and δ¹⁵N values (Mann-Whitney test: U = 84.000, Z = -0.505, P = 0.614) between group C and the wild herbivores are not significantly different, demonstrating that the domestic sheep/goats and cattle in group C were likely free-range animals that ate wild C₃ vegetation [67]. Another possibility is that some of these animals actually represent wild bovids native to the region such as gaur, takin, serow, or goral, rather than domestic bovids [70]. The diverse isotopic signatures for the domesticated herbivores, representing three distinct dietary groups, may reflect three different husbandry strategies used by people at Shilinggang. The subtropical natural environment in northwestern Yunnan may have provided abundant year-round forage for most herd animals. However, people also provided supplemental fodder to some sheep, goats, and cattle using millet agricultural products. A few individuals were foddered very intensively with millet or may have grazed in millet fields.

Two Old World monkey samples (Cercopithecidae) have δ¹³C and δ¹⁵N values of -18.9‰ and 8.0‰ (#84) and -18.8‰ and 10.6‰ (#85), respectively, showing that they are characterized by a predominantly C₃ omnivorous diet with high-levels of protein. The stable isotope values are consistent with the diet of most monkeys in the region, which is composed of C₃ plants (including tender leaves, sprouts, flowers, fruit, and seeds), insects which feed on C₃ plants, and bird eggs [71–72]. The δ¹⁵N value of sample #84 approaches those in omnivorous nonhuman primates, which suggests that its diet included insects or vertebrate fauna [73–74]. Sample #85 has approximately 3‰ higher δ¹⁵N indicating an offset of one trophic level from sample #84. This sample was from an unfused distal ulna, indicating that it came from a sub-adult animal. The higher trophic level may be due to the consumption of maternal milk [e.g. 73–75], with this sample representing a homogenized signal from the first few years of life due to slow bone turnover rates.

Human stable isotope data and paleodietary analysis at Shilinggang

The carbon and nitrogen data for the 16 human samples from Shilinggang are listed in Tables 2 and 3, and plotted in Figs 3 and 5. The δ¹³C values for the human bones range from -19.5‰ to -16.3‰ (mean = -18.8±0.8‰), and the δ¹⁵N values vary from 7.5‰ to 11.3‰ (mean = 10.0±1.0‰). In accordance with previous studies [18–19, 55], human δ¹³C signatures can be divided into three main dietary groups: mainly C₃ consumers (< -18‰); mixed C₃ and C₄ consumers (-18‰ to -12‰); and highly C₄ dependent consumers (> -12‰). The population of Shilinggang appears to have been predominantly C₃ consumers. It is likely that people’s diet consisted primarily of rice [28], as well as roots, tubers, fruits, and/or animals fed on a C₃ diet. However, there is one individual (#130) with a δ¹³C of -16.3‰ who consumed a small amount of C₄ foodstuffs, likely foxtail millet and/or animals fed on millet [28]. The mean δ¹⁵N value for humans is 4.8‰ higher than that of wild-herbivores and 3.1‰ higher than that of omnivores, suggesting that people consumed abundant animal protein. People at Shilinggang might have consumed both wild animals available in the region or domesticated livestock kept near the site. Consumption of freshwater fish could also lead to a high δ¹⁵N value [76], which may be a complementary reason for the high δ¹⁵N values in the human samples. However, very few fish bones have been found in the excavation of Shilinggang. Zooarchaeological analysis of animal remains from Shilinggang indicates that domesticates were the most common animals present, making up about 60% of identified specimens. However, a diverse array of wild fauna including several types of deer, boar, badger, porcupine, bear, monkeys, birds, and fish were also exploited in small numbers. Together the isotopic, paleobotanical, and zooarchaeological data indicate that people at Bronze Age Shilinggang developed a broad-spectrum economic strategy that included agriculture (planting rice, millet, roots, and tubers), animal husbandry (raising pigs, dogs, sheep/goats, and cattle), hunting wild animals, and freshwater fishing.

[Figure omitted. See PDF.]

Table 3. Summary of human and animal isotopic data.

Adults include samples from individuals over 25 years old. Children include two bone samples from individuals under 25 years old. SD: Standard deviation.

https://doi.org/10.1371/journal.pone.0177867.t003

[Figure omitted. See PDF.]

Fig 5. Comparison of human isotopic signatures for adult vs. child samples.

https://doi.org/10.1371/journal.pone.0177867.g005

We divided human samples into two groups based on age (adult vs. child) to see if there are demographic differences in dietary signatures. Two bone samples from children were estimated to be around 4 to 5 years old at the time of death, which provide some early childhood dietary information for the Shilinggang population. The child samples have δ¹³C values ranging from -19.1‰ to -18.5‰ (mean = -18.9±0.2‰) and δ¹⁵N values ranging from 11.1‰ to 11.3‰ (mean = 11.2±0.1‰). The other 14 adult samples have δ¹³C values ranging from -19.5‰ to -16.3‰ (mean = -18.8±0.2‰) and are not significantly different from the two samples representing childhood diets (Mann-Whitney test: U = 13.000, Z = -0.159, P = 0.933). This indicates that people at Shilinggang consumed C₃ foodstuffs during both childhood and adulthood. The weak signal of difference between adults and children (0.1‰) may simply result from the small sample size for children. Alternately, if children consumed more C₃ plant food after being weaned this could also weaken the signal. The δ¹⁵N values for the 14 adults range from 7.5‰ to 11.1‰ (mean = 9.8±0.2‰) and display a significant difference from the two samples representing children’s diets (Mann-Whitney test, U = 1.000, Z = -2.064, P = 0.033). The samples representing children’s diets have enriched δ¹⁵N values by 1.4‰, which could likely be attributed to the consumption of breast milk [e.g. 74–75], but could also be affected by the small sample size.

As Fig 3 and Table 2 show, there are two outliers (#128 and #130) among the human samples. Sample #128 has a δ¹³C value of -18.9‰ and a δ¹⁵N value of 7.5‰, which shows that this individual had a predominately C₃ diet similar to most other people at the site. However, the δ¹⁵N value for this individual is 2.5‰ less than the average value of all human samples (10.0±1.0‰), indicating that this person consumed less animal protein. There are several possible explanations. This person may have had low social status (e.g. a slave, prisoner of war, criminal, etc.), and her diet may have included less protein than other members of the population [e.g. 35, 77–78]. The diet difference could also represent personal preference or a cultural choice. Perhaps she did not like eating animal products or was a migrant to Shilinggang from other regions where people ate less meat. Sample #130 has a δ¹⁵N value of 9.3‰, which is similar to other individuals at the site, but the δ¹³C value of -16.3‰ is relatively high, suggesting a mixed C₃ and C₄ diet. This person was not only consuming rice, roots, and/or tubers and wild plants (C₃ plants), but also some foxtail millet (C₄ plant). The mixed paleodietary signature from this sample is accordant with archaeobotanical evidence showing that people made use of diverse plant food resources at Shilingang [28].

Human diets and animal utilization strategies during the Bronze Age in Yunnan

In addition to our current study, two additional stable isotope analyses of human remains from Yunnan have been published for the Jinlianshan cemetery (2700–2300 cal BP) and Mayutian site (2500–2400 cal BP) (Fig 1) [24, 35]. No stable isotopic data for faunal remains have been reported from these two sites. Therefore, our analysis of human and animal bone samples from Shilinggang provides an important isotopic baseline for understanding human and animal dietary structures in Yunnan during the Bronze Age. Stable isotope results from our current study as well as those from Jinlianshan and Mayutian indicate that people may have had varied diets in different regions of Yunnan. For example, at Jinlianshan bone collagen δ¹³C values have a limited range from -19.3‰ to -18.2‰ (mean = -18.8±0.4‰, n = 9) [35, Table 4]. This is quite different from the large range of bone collagen δ¹³C values from -19.5‰ to -16.3‰ (mean = -18.8±0.8‰, n = 16) that we found at Shilinggang. Despite the greater range of δ¹³C values at Shilinggang, the mean δ¹³C values (Mann-Whitney test: U = 68.000, Z = -0.226, P = 0.846) and δ¹⁵N values (Mann-Whitney test: U = 59.000, Z = -0.736, P = 0.487) between Shilinggang and Jinlianshan have no significant difference. When compared to Jinlianshan and Shilinggang, the population at the site of Mayutian appears to have had a very different diet made up of both C₃ and C₄ resources, but statistical comparison is not possible because isotope data from Mayutian is derived from apatite rather than collagen [24]. More data is still needed, but the current evidence suggests that Bronze Age groups at Shilinggang and in the Jinlianshan cemetery consumed C₃-based foods, while those at Mayutian mainly had a mixed diet of both C₃ and C₄ resources [24, 35]. The high proportion of C₃ plants in people’s diets in this part of China may also indicate that in addition to growing domestic C₃ plants such as rice, people also took advantage of the wild plant resources available locally.

[Figure omitted. See PDF.]

Table 4. Human collagen carbon and nitrogen isotope results from sites in Yunnan.

https://doi.org/10.1371/journal.pone.0177867.t004

We do not yet have a complete understanding of the regional variation in agricultural practices and animal husbandry among Yunnan’s different culture groups. At Shilinggang, rice and foxtail millet are the only domestic crops that have been identified in the paleobotanical record [28]. At other Bronze Age sites such as Haimenkou and Xueshan (Fig 1), rice, wheat, barley, foxtail millet, and broomcorn millet have all been identified [79–80]. This indicates that people in Bronze Age Yunnan cultivated a wide variety of domestic crops. Because carbon and nitrogen isotopic data from animal bone collagen have not been reported from Jinlianshan and Mayutian, estimating the contribution of animal foods to ancient human diets at these two sites is impossible. Our comparison of nitrogen isotopic values for human and animal bone collagen at Shilinggang indicates that people consumed a large amount of animal protein during the late Bronze Age, and that it came primarily from C₃ fed domestic or wild animals. Our analysis also identified dietary groups such as pig group A and sheep/goats and cattle group C that may represent the additional exploitation of wild boar and wild bovids.

Bronze Age animal exploitation in Yunnan seems to have focused on domestic animal husbandry supplemented by hunting of wild animals. Zooarchaeological analyses of other archaeological sites in the region have identified similar domestic livestock and wild animals as we have found at Shilinggang [30–31]. However, the ways that humans raised domestic animals in this region have not been discussed in detail. The isotopic signatures for faunal remains from Shilinggang suggest that humans might have raised livestock using multiple herd management strategies. We found that the δ¹³C values for domesticated sheep/goats and cattle suggest they had three different kinds of diets containing variable amounts of C₄ plants (Fig 4). Pigs can also be split into three similar dietary groups that may reflect wild vs. domestic populations and different amounts of millet by-products in the diet.

The animal dietary signatures at Shilinggang may also reflect seasonal changes in herding practices or differences in the types of grazing present in different micro-regions. The landscape of the Nujiang River Valley changes dramatically with increasing altitude. During the Bronze Age, rain-fed millet crops may have been preferentially cultivated in the highlands, while rice may have been planted in the lowlands [28]. Accordingly, people might have herded animals in different ways at different points along the slopes or at different times of year, feeding animals with the by-products of the plants grown at various elevations. Transhumance may explain some of the variation in carbon and nitrogen isotopic values between the sheep/goats and cattle groups at Shilinggang, and future isotopic studies of the region should consider this possibility in more detail.

Although more data from additional sites is needed, we propose that regional cultures in different areas of Yunnan developed subsistence strategies according to local environmental conditions. Yunnan has highly diverse and isolated ecological environments split by numerous large rivers, such as the Nujiang, Lancang, and Jinsha Rivers and the Hengduan Mountains (Fig 1) [81]. Mountainous terrain with limited flat land for planting crops likely promoted the development of independent cultural systems within subregions, and the exploitation of abundant wild plant and animal resources. These geographic barriers may have also hindered cultural and technological exchange between regions [81]. Research in other parts of the world has shown that geographic variation results in diversification of subsistence strategies [e.g. 82–86]. Our research provides a starting point for examining the relationships between environmental diversity and subsistence in Yunnan. In the coming years, additional isotopic, paleobotanic, and zooarchaeological studies will further clarify the subsistence systems and cultural variation within this unique region.

Conclusion

Stable isotope data from Shilinggang reveals the diverse subsistence strategies and animal husbandry techniques that were used during the late Bronze Age in southwest China. Dogs and pigs had diets that consisted of mainly C₃ foodstuffs and were high in animal protein. Their diets were similar to those of humans at the site, likely resulting from the consumption of human food waste. However, pigs can be divided into three separate groups, including pig A representing wild boar and pig B and C two types of domestic individuals based on their isotopic signatures. Sheep/goats and cattle specimens can also be divided into three separate dietary groups. These dietary differences may reveal the presence of wild vs. domestic populations, but also indicate that some herd animals were more intensively foddered with millet by-products. Humans generally consumed a large amount of C₃ foodstuffs and a large amount of animal protein, supplemented with C₄ foods. Together the data suggest that people at Shilinggang developed diversified economic strategies that focused on agriculture and domestic animal husbandry, supplemented by collecting wild plants, hunting, and fishing. We propose that the subsistence systems used by different cultural groups in the various regions of Yunnan may have been adapted to local environmental conditions, resulting in considerable spatial differentiation of Bronze Age cultures in this region.

Supporting information

[Figure omitted. See PDF.]

S1 Text. Hierarchical cluster analysis.

https://doi.org/10.1371/journal.pone.0177867.s001

(DOCX)

S1 Fig. Operation process and parameter setting of hierarchical cluster analysis.

https://doi.org/10.1371/journal.pone.0177867.s002

(TIF)

S2 Fig. Hierarchical cluster analysis of collagen d¹⁵N and d¹³C values from pig samples.

The three main groups are labeled A-C. The dendrogram was constructed using average linkage (between groups).

https://doi.org/10.1371/journal.pone.0177867.s003

(TIF)

S3 Fig. Hierarchical cluster analysis of collagen d¹⁵N and d¹³C values from sheep/goats and cattle samples.

The three main groups are labeled A-C. The dendrogram was constructed using average linkage (between groups).

https://doi.org/10.1371/journal.pone.0177867.s004

(TIF)

Acknowledgments

We thank Professor Yaowu Hu from the University of the Chinese Academy of Sciences for constructive comments and suggestions. We also thank Editage (http://online.editage.cn/) for English editing services. We are greatly indebted to the editor and the two anonymous reviewers for their constructive suggestions and revisions.

Author Contributions

1. Formal analysis: LLR WMD GHD.

2. Funding acquisition: LLR GHD.

3. Investigation: LLR X Li LHK HGL HML RM X Liu.

4. Methodology: LLR X Li WMD GHD.

5. Project administration: GHD.

6. Resources: LHK RM X Liu.

7. Software: LLR WMD.

8. Supervision: GHD.

9. Visualization: GHD.

10. Writing – original draft: LLR WMD GHD.

11. Writing – review & editing: LLR GHD KB.

Citation: Ren L, Li X, Kang L, Brunson K, Liu H, Dong W, et al. (2017) Human paleodiet and animal utilization strategies during the Bronze Age in northwest Yunnan Province, southwest China. PLoS ONE 12(5): e0177867. https://doi.org/10.1371/journal.pone.0177867

References

1. Chen FH, Dong GH, Zhang DJ, Liu XY, Jia X, An CB, et al. Agriculture facilitated permanent human occupation of the Tibetan Plateau after 3600 BP. Science. 2015; 347(6219): 248–250. pmid:25593179

2. d’Alpoim Guedes J, Lu HL, Hein AM, Schmidt AH. Early evidence for the use of wheat and barley as staple crops on the margins of the Tibetan Plateau. Proc Natl Acad Sci U S A. 2015; 112(18): 5625–5630. pmid:25902511

3. Farley NW, Moll CJ, Blacksher LM. Environment, Diet, and Craniofacial Development: A Study of Mixed Subsistence Strategies in the Great Lakes Watershed, AD 900–1600. Int J Humanit Soc Sci. 2015; 5(9): 29–45.

4. Roberts P, Henshilwood CS, van Niekerk KL, Keene P, Gledhill A, Reynard J, et al. Climate, Environment and Early Human Innovation: Stable Isotope and Faunal Proxy Evidence from Archaeological Sites (98-59ka) in the Southern Cape, South Africa. PLoS One. 2016; 11(7): e0157408. pmid:27383620

5. Jia X, Dong GH, Li H, Brunson K, Chen FH, Ma MM, et al. The development of agriculture and its impact on cultural expansion during the mid-late Neolithic in the Western Loess Plateau, China. The Holocene. 2013; 23(1): 85–92.

6. Jia X, Sun YG, Wang L, Sun WF, Zhao ZJ, Lee HF, et al. The transition of human subsistence strategies in relation to climate change during the Bronze Age in the West Liao River Basin, Northeast China. The Holocene. 2016;

7. Li XQ, Shang X, Dodson J, Zhou XY. Holocene agriculture in the Guanzhong Basin in NW China indicated by pollen and charcoal evidence. The Holocene. 2009; 19(5): 1213–1220.

8. Somerville AD, Sugiyama N, Manzanilla LR, Schoeninger MJ. Animal Management at the Ancient Metropolis of Teotihuacan, Mexico: Stable Isotope Analysis of Leporid (Cottontail and Jackrabbit) Bone Mineral. PLoS One. 2016; 11(8): e0159982. pmid:27532515

9. Dong GH, Yang YS, Han JY, Wang H, Chen FH. Exploring the history of cultural exchange in prehistoric Eurasia from the perspectives of crop diffusion and consumption. Sci China Earth Sci. 2017;

10. Dong GH, Ren LL, Jia X, Liu XY, Dong SM, Li HM, et al. Chronology and subsistence strategy of Nuomuhong Culture in the Tibetan Plateau. Quat Int. 2016;

11. Li G. Subsistence of Neolithic Pearl River Estuary Area, South China. M. A. Thesis, The University of British Columbia. 2002. https://open.library.ubc.ca/cIRcle/collections/ubctheses/831/items/1.0058381.

12. Wang TT, Fuller BT, Wei D, Chang XE, Hu YW. Investigating Dietary Patterns with Stable Isotope Ratios of Collagen and Starch Grain Analysis of Dental Calculus at the Iron Age Cemetery Site of Heigouliang, Xinjiang, China. Int J Osteoarchaeol. 2015;

13. Erickson CL. Domestication and Subsistence Implications of Plant and Animal Utilization in the Titicaca Basin. Biochem J. 2006; 396(1): 117–126.

14. Hou LL, Hu YW, Zhao XP, Li ST, Wei D, Hou YF, et al. Human subsistence strategy at Liuzhuang site, Henan, China during the proto-Shang culture (~2000–1600 BC) by stable isotopic analysis. J Archaeol Sci. 2013; 40(5): 2344–2351.

15. Liu L, Ge W, Bestel S, Jonesd D, Shi JM, Song YH, et al. Plant exploitation of the last foragers at Shizitan in the Middle Yellow River Valley China: evidence from grinding stones. J Archaeol Sci. 2011; 38(12): 3524–3532.

16. Wang C, Lu H, Zhang J, He K, Huan X. Macro-Process of Past Plant Subsistence from the Upper Paleolithic to Middle Neolithic in China: A Quantitative Analysis of Multi-Archaeobotanical Data. PLoS One. 2016; pmid:26840560

17. Atahan P, Dodson J, Li XQ, Zhou XY, Hu SM, Fiona B, et al. Subsistence and the isotopic signature of herding in the Bronze Age Hexi Corridor, NW Gansu, China. Quat Int. 2011; 38(7): 1747–1753.

18. Barton L, Newsome SD, Chen FH, Wang H, Guiderson TP, Bettinger RL. Agricultural origins and the isotopic identity of domestication in northern China. Proc Natl Acad Sci U S A. 2009; 106(14): 5523–5528. pmid:19307567

19. Ma MM, Dong GH, Lightfoot E, Wang H, Liu XY, Jia X, et al. Stable isotope analysis of human and faunal remains in the western Loess Plateau, approximately 2000 cal BC. Archaeometry. 2014; 56(S1): 237–255.

20. Ma MM, Dong GH, Liu XY, Lightfoot E, Chen FH, Wang H, et al. Stable isotope analysis of human and animal remains at the Qijiaping site in middle Gansu, China. Int J Osteoarchaeol. 2015; 25(6): 923–934.

21. Hu Y, Ambrose SH, Wang C. Stable isotopic analysis of human bones from Jiahu site, Henan, China: implications for the transition to agriculture. J Archaeol Sci. 2006; 33(9): 1319–1330.

22. Hu YW, Wang GF, Cui YP, Dong Y, Guan L, Wang CS. A study on the dietary of the ancestors of Sanxingcun site in Jintan, Jiangsu province. Chinese Sci Bull. 2007; 52(1): 85–88 (In Chinese).

23. Hu YW, Li FJ, Wang CS, Richard MP. Carbon and nitrogen stable isotope analysis of the human bones from the Liyutun site, Zhanjiang, Guangdong: A preliminary exploration of the Neolithic human lifestyle in south China. Acta Anthropologica Sinica. 2010; 29(3): 264–269 (In Chinese).

24. Zhang X, Burton J, Jin Z, Xiao M, Fan A, Xu JF, et al. Isotope studies of human remains from Mayutian, Yunnan Province, China. J Archaeol Sci. 2014; 50(1): 414–419.

25. Yunnan Provincial Institute of Cultural Relics and Archaeology (YNICRA). Archaeological papers of Yunnan Province. Kunming: Yunnan nationalities publishing house; 1998 (In Chinese).

26. Fan Y. Study on the category of Bronze culture in Yunnan. Sichuan Cultural Relics. 2007; 2: 64–73 (In Chinese).

27. Jin HT, Liu X, Min R, Li XR, Wu XH. Results and analysis of flotation in Dadunzi site in Yuanmou county, Yunnan Province. Jianghan Archaeology. 2014; 3: 109–114 (In Chinese).

28. Li HM, Zuo XX, Kang LH, Ren LL, Liu FW, Liu HG, et al. Prehistoric agriculture development in the Yunnan-Guizhou Plateau, southwest China: archaeobotanical evidence. Sci China Earth Sci. 2016; 1–12.

29. Zhao Z J. The flotation results report of Mopandi site, Yunnan Province. Yunnan Province Institute of Cultural Relics and Archaeology. The 2001 excavation report appendix of Caiyuanzi and Mopandi sites, Yunnan Province. Acta Archaeologica Sinica. 2003; 2: 28–39 (In Chinese).

30. Zhao Y. A research of faunal remains excavated from Yinsuodao in Yunnan province. M.Sc. Thesis, Jilin University, Changchun. 2011. paperuri: (a0f92250bdc734f1ff8d7638c0d38b4a) (In Chinese).

31. Yunnan Provincial Institute of Cultural Relics and Archaeology (YNICRA). Shifodong in Gengma. Beijing: Cultural relics press; 2010 (In Chinese).

32. Ambrose SH, Norr L. Experimental evidence for the relationship of the carbon isotope ratios of whole diet and dietary protein to those of bone collagen and carbonate. In: Joseph B, Lambert GG, editors. Prehistoric human bone. Berlin: Springer Berlin Heidelberg; 1993. pp. 1–37.

33. Privat KL, O'Connell TC, Richards MP. Stable isotope analysis of human and faunal remains from the Anglo-Saxon cemetery at Berinsfield, Oxfordshire: dietary and social implications. J Archaeol Sci. 2002; 29: 779–790.

34. Pearson JA, Bogaard A, Charles M, Hillson SW, Larsen CS, Russell N, et al. Stable carbon and nitrogen isotope analysis at Neolithic Çatalhöyük: evidence for human and animal diet and their relationship to households. J Archaeol Sci. 2015; 57: 69–79.

35. Zhang QC. Stable isotope preliminary analysis of human bones of Jinlianshan tombs in Chengjiang county, Yunnan province. Archaeology. 2011; 1: 007 (In Chinese).

36. Kang LH. Bronze civilization on the bank of Nujiang River-excavation in Shilinggang site in Lushui county, Yuannan province. Popular archaeology. 2015; 6:12–13 (In Chinese).

37. Zhao DY, Zhu H, Kang LH. Paleopathological study of human remains from the Bronze age Shilinggang site in Nujiang, Yunnan, southwest China. Asian archaeology. 2015; 3: 129–144.

38. National Meteorological Information Center (CMA Meteorological Data Center); 2005 [cited 2016 Sep 1]. Database: dataset of monthly surface observation values in individual years (1981–2010) in China [Internet]. http://data.cma.cn/en/?r=data/detail&dataCode=A.0029.0004.

39. Wang CY, Ou XK, He ZR, Zhu WM, Zhang ZP, Su WH, et al. Analysis on vegetation diversity in northwestern Yunnan, China. Chinese Journal of Ecology. 2001; 20: 21–25 (In Chinese).

40. Fang T. Yunnan's diet culture and Yunnan's historical development. Dietetic Culture Research. 2007; 3: 37–48 (In Chinese).

41. Buikstra JE, Ubelaker DH, Glassman DM. Standards for data collection from human skeletal remains. Hum Biol. 1996; 68(1): 155–158.

42. Richards MP, Hedges RE. Stable isotope evidence for similarities in the types of marine foods used by Late Mesolithic humans at sites along the Atlantic coast of Europe. J Archaeol Sci. 1999; 26(6): 717–722.

43. Ambrose SH, Butler BM, Hanson DB, Hunter-Anderson RL, Krueger HW. Stable isotopic analysis of human diet in the Marianas Archipelago, Western Pacific. Am J Phys Anthropol. 1997; 104(3): 343–361. pmid:9408540

44. Price TD, Blitz J, Burton J, Ezzoa JA. Diagenesis in prehistoric bone: problems and solutions. J Archaeol Sci. 1992; 19(5): 513–529.

45. Ambrose SH. Preparation and characterization of bone and tooth collagen for isotopic analysis. J Archaeol Sci. 1990; 17(4): 431–451.

46. Van Klinken GJ. Bone collagen quality indicators for palaeodietary and radiocarbon measurements. J Archaeol Sci. 1999; 26(6): 687–695.

47. DeNiro MJ. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature. 1985; 317: 806–809.

48. Kohn MJ. "You are what you eat." Science.1999; 283(5400): 335–336. pmid:9925492

49. Schwarcz HP, Schoeninger MJ. Stable isotope analyses in human nutritional ecology. Am J Phys Anthropol. 1991; 34(S13): 283–321.

50. Ambrose SH, Norr L. Experimental evidence for the relationship of the carbon isotope ratios of whole diet and dietary protein to those of bone collagen and carbonate. In: Lambert JB, Grupe B, editors. Prehistoric human bone. New York: Spring Berlin Heidelberg; 1993. pp. 1–37.

51. Van der Merwe NJ, Vogel JC. 13C content of human collagen as a measure of prehistoric diet in woodland North America. Nature. 1978; 276: 815–816. pmid:364321

52. Schoeninger MJ, DeNiro MJ. Nitrogen and carbon isotopic composition of bone collagen from marine and terrestrial animals. Geochim Cosmochim Acta. 1984; 48(4): 625–639.

53. O'Leary MH. Carbon isotope fractionation in plants. Phytochemistry. 1981; 20(4): 553–567.

54. McGovern PE, Zhang J, Tang J, Zhang Z, Hall GR, et al. Fermented beverages of pre-and proto-historic China. Proc Natl Acad Sci U S A. 2004; 101(51): 17593–17598. pmid:15590771

55. Pechenkina EA, Ambrose SH, Ma XL, Benfer RA. Reconstructing northern Chinese Neolithic subsistence practices by isotopic analysis. J Archaeol Sci. 2005; 32(8): 1176–1189.

56. Bocherens H, Drucker D. Trophic level isotopic enrichment of carbon and nitrogen in bone collagen: case studies from recent and ancient terrestrial ecosystems. Int J Osteoarchaeol. 2003; 13(1–2): 46–53.

57. O'Connell TC, Kneale CJ, Tasevska N, Kuhnle GG. The diet-body offset in human nitrogen isotopic values: A controlled dietary study. Am J Phys Anthropol. 2012; 149(3): 426–434. pmid:23042579

58. Ambrose SH, DeNiro MJ. Bone nitrogen isotope composition and climate. Nature. 1987; 325: 201–201.

59. Ambrose SH. Effects of diet, climate and physiology on nitrogen isotope abundances in terrestrial foodwebs. J Archaeol Sci. 1991; 18(3): 293–317.

60. Hobson KA, Alisauskas RT, Clark RG. Stable-nitrogen isotope enrichment in avian tissues due to fasting and nutritional stress: implications for isotopic analyses of diet. The Condor. 1993; 95: 388–394.

61. Mays SA, Richards MP, Fuller BT. Bone stable isotope evidence for infant feeding in Mediaeval England. Antiquity. 2002; 76(293): 654–656.

62. Fuller BT, Fuller JL, Sage NE, Harris DA, O'Connell TC, Hedges RE. Nitrogen balance and δ15N: why you're not what you eat during nutritional stress. Rapid Commun Mass Spectrom. 2005; 19(18): 2497–2506. pmid:16106342

63. Bogaard A, Heaton TH, Poulton P, Merbach I. The impact of manuring on nitrogen isotope ratios in cereals: archaeological implications for reconstruction of diet and crop management practices. J Archaeol Sci. 2007; 34(3): 335–343.

64. Layman CA, Araujo MS, Boucek R, Hammerschlag-Peyer CM, Harrison E, Jud ZR, et al. Applying stable isotopes to examine food-web structure: an overview of analytical tools. Biol Rev. 2012; 87(3): 545–562. pmid:22051097

65. Hedges REM, Reynard LM. Nitrogen isotopes and the trophic level of humans in archaeology. J Archaeol Sci. 2007; 34(8): 1240–1251.

66. Sage RF, Wedin DA, Li M. The biogeography of C4 photosynthesis: patterns and controlling factors. Physiological Ecology. 1999; 313–373.

67. Shen J, Jones RT, Yang X, Dearing JA, Wang S. The Holocene vegetation history of Lake Erhai, Yunnan province southwestern China: the role of climate and human forcings. The Holocene. 2006; 16(2): 265–276.

68. Hu YW, Luan FS, Wang SG, Wang CS, Richards MP. The initial attempt to identify domesticated pigs and wild boar using stable carbon and nitrogen isotope analysis. Science China Series D. 2008; 38(6): 693–700 (In Chinese).

69. Cucchi T, Dai L, Balasse M, Zhao C, Gao J, Hu Y, et al. Social Complexification and Pig (Sus scrofa) Husbandry in Ancient China: A Combined Geometric Morphometric and Isotopic Approach. PloS one. 2016; 11(7): e0158523. pmid:27384523

70. Smith AT, Xie Y, Hoffmann RS, Lunde D, MacKinnon J, Wilson DE, Wozencraft WC. A guide to the mammals of China. Princeton: Princeton University Press; 2010.

71. He XR, Yang DH. Preliminary studies on the biology of Presbytis phayrei. Zoological research. 1982; S2: 349–354 (In Chinese).

72. Zhang Y, Chen L, Qu W, Coggins C. The Primates of China: Biogeography and Conservation Status-Past, Present, and Future. Beijing: China Forestry Publishing House; 2002. pp. 20–22.

73. Schoeninger MJ, Iwaniec UT, Glander KE. Stable isotope ratios indicate diet and habitat use in New World monkeys. Am J Phys Anthropol. 1997; 103(1): 69–83. pmid:9185952

74. Schoeninger MJ, Jim M, Sept JM. Subsistence strategies of two "savanna" chimpanzee populations: the stable isotope evidence. Am J Primatol. 1999; 49(49): 297–314.

75. Fuller BT, Fuller JL, Harris DA, Hedges REM. Detection of breastfeeding and weaning in modern human infants with carbon and nitrogen stable isotope ratios. Am J Phys Anthropol. 2006; 129: 279–293. pmid:16261548

76. Zanden MJV, Cabana G, Rasmussen JB. Comparing trophic position of freshwater fish calculated using stable nitrogen isotope ratios (δ15N) and literature dietary data. Can J Fish Aquat Sci. 1997; 54(5): 1142–1158.

77. Shelach G. The Qiang and the question of human sacrifice in the late Shang period. Asian Perspectives. 1996; 1–26.

78. Yuan J, Flad R. New zooarchaeological evidence for changes in Shang Dynasty animal sacrifice, J Anthropol Archaeol. 2005; 24: 252–270.

79. Xue YN. Pilot study of plant remains about Haimenkou site in Jianchuan County, Yunnan Province. M.A. Thesis, Shandong University, Jinan. 2010 (in Chinese). Available: paperuri: (c3d57aa2372d9453f71a2f4a2218b1d0).

80. Wang Q. Plant remains analysis at Xueshan site in Chengjiang county, Yunnan Province. M. A. Thesis, Shandong University, Jinan. 2014 (in Chinese). Available: paperuri: (c3d57aa2372d9453f71a2f4a2218b1d0).

81. Li SM. Research on the evolution of early civilization in Yunnan Province. Kunming: The Peoples Press of Yunnan. 2013; pp. 54–56 (In Chinese).

82. Mansur ME, Piqué R. Between the Forest and the Sea: Hunter-Gatherer Occupations in the Subantarctic Forests in Tierra del Fuego, Argentina. Arctic Anthropol. 2009; 82(1–2): 342–343.

83. Drass RR, Flynn P. Temporal and geographic variations in subsistence practices for plains villagers in the southern Plains. Plains Anthropol. 1990; 175–190.

84. Collard IF, Foley RA. Latitudinal patterns and environmental determinants of recent human cultural diversity: do humans follow biogeographical rules? Evol Ecol Res. 2002; 4(3): 371–383.

85. Nardoto GB, Silva S, Kendall C, Ehleringer JR, Chesson LA, Ferraz ES, et al. Geographical patterns of human diet derived from stable-isotope analysis of fingernails. Am J Phys Anthropol. 2006; 131(1): 137–146. pmid:16552735

86. Ma MM, Dong GH, Jia X, Wang H, Cui YF, Chen FH. Dietary shift after 3600 cal yr BP and its influencing factors in northwestern China: Evidence from stable isotopes. Quat Sci Rev. 2016; 145: 57–70.