Varying cultivation strategies in eastern Tianshan corresponded to growing pastoral lifeways between 1300 BCE and 300 CE
Duo Tian1† 
Yufeng Sun2† 
Melissa M. Ritchey2 Tongyuan Xi1 Meng Ren1 Jian Ma1 Jianxin Wang3 
Zhijun Zhao1,3 Xue Ling1* 
Xinyi Liu2*- 1Key Laboratory of Cultural Heritage Research and Conservation, China-Central Asia “the Belt and Road” Joint Laboratory on Human and Environment Research, School of Culture Heritage, Northwest University, Xi’an, China
- 2Department of Anthropology, Washington University in St. Louis, St. Louis, MO, United States
- 3Institute of Archaeology, Chinese Academy of Social Sciences, Beijing, China
This study combines plant stable isotope and archaeobotanical analyses to explore how ancient pastoral communities in varying landscapes of eastern Tianshan managed their barley fields. The question is less archaeologically investigated, as recent discussions have focused on pastoral and nomadic activities. Results show that diversified cultivation strategies were employed in barley cultivation at different locations in eastern Tianshan. We also observed a diachronic transition toward less labour-intensive crop management corresponding to a growing pastoral lifeway from the late Bronze Age (1300–800 BCE) to historical periods (400 BCE–300 CE). These results inform us about the mechanism by which southwest Asian originated domesticates were adapted to the Inner Asian environments in the context of the early food globalisation.
Introduction
In this article, we explore the cultivation strategies employed by ancient communities in the eastern Tianshan foothills between 1,300 BCE and 300 CE. Did the herders carry out intensive barley (Hordeum vulgare) cultivation with irrigation and manuring in pastoral-optimal environments? One of the recent developments in Old World archaeology has been the realisation of the movements of cultigens (and livestock) East and West across the Eurasian landmass thousands of years before the historic Silk Route (e.g., Jones et al., 2011; Liu et al., 2019). Accumulating evidence suggests an eastern dispersal of the Fertile Crescent cereals, notably free threshing wheat (Triticum aestivum/durum) and barley, was paralleled by the western movement of broomcorn and foxtail millet (Panicum miliaceum and Setaria italica) across the continental interior during the second millennium BCE. Recent analyses ascertained the chronology of the cereal movement (e.g., Liu et al., 2017; Zhou et al., 2020) and established the importance of Inner Asian Mountain Corridors (Frachetti, 2012; Spengler, 2015). However, the mechanism by which novel crops were adapted to the existing agropastoral landscapes remains unclear and is subjected to cultural-specific investigations (e.g., Vaiglova et al., 2020; Li et al., in press).
There has been a considerable recent development in the eastern Tianshan mountains with intensive field investigation and specialist research, resulting in a better understanding of the settlement pattern and subsistence strategies of a series of agropastoral sites commonly referred to as the eastern Tianshan complexity (Ma, 2014; Ma et al., 2014a,b; Xi et al., 2020a,b). These communities were herders with sheep, goats, cattle, and horses who also carried out barley and millet cultivation (Ling et al., 2013, 2016; Li et al., 2020, 2021; Tian et al., 2021a,b). Central to our inquiry is the conceptual conflict between pastoral and cultivation activities and the labour arrangement to facilitate both. Pastoral economies employ extensive labour approaches, in which productivity is driven by the size of herds and the quality and expansiveness of the pastures. Cultivation in traditional farming society, on the other hand, is often limited by land availability, and within this context, productivity is driven by intensive labour input. Recent excavations at eastern Tianshan and burgeoning plant stable isotope analysis provide a unique opportunity to measure the historically less investigated cultivation strategies in pastoral landscapes.
Stable carbon and nitrogen isotope compositions measured in archaeobotanical remains are informative to the plants’ growing conditions. Particularly, δ13C values of C3 plants can help infer crop water status (Farquhar et al., 1989; Araus et al., 1997; Flohr et al., 2011; Wallace et al., 2013). Plant δ15N values enable the assessment of soil 15N enrichment caused by aridity, soil denitrification, and anthropogenic activities such as manuring (Bogaard et al., 2007; Fraser et al., 2011). Macrofossil analysis focusing on weed taxonomy enables insights into field conditions and crop processing (e.g., Hillman, 1984; Stevens, 2003; Harvey and Fuller, 2005). These approaches allow direct assessment of past crop water status and growing conditions, which have shed light on early crop management in Europe, the Middle East, and East Asia (Bogaard et al., 2016, 2018; Styring et al., 2017; Alagich et al., 2018; Li et al., in press). Here we apply a similar approach to the study of agricultural practices in the eastern Tianshan Mountains.
Eastern Tianshan complexity and the study area
Situated at the easternmost part of the Inner Asian mountain ranges, the eastern Tianshan has been a cultural crossroad since antiquity, connecting the Hexi Corridor and Gobi Desert to the eastern Eurasian Steppe and Inner Asian mountains. In the past decade, intensive field investigations led by a team of archaeologists from Northwest University, Xi’an, revealed abundant archaeological remains from numerous sites dating between 2000 BCE and 1000 CE and transformed our understanding of the region’s early history (Ma, 2014). Most sites are found in three counties, Barkol, Hami, and Yiwu, located along the northern slopes of the Tianshan Mountains with altitudes ranging between 1,600 and 2,200 metres above sea level (m.a.s.l.). This region has highly diverse ecosystems, consisting of mountain steppes, coniferous forests, moraine hills, alluvial fans, and saline lakes. Nowadays, the climate is characterised by a relatively short-cool summer (17°C in July on average) and a long-frozen period (−17.5°C in January on average) and is generally considered a semi-arid environment (100–200 mm for annual precipitation) (Zhang, 1993; Liu, 1994). Modern residents practice spring-sown agriculture and seasonal transhumant pastoralism between foothill and mountainous environments.
Isotopic values measured on skeletal remains suggest a C3 plant-based human diet at Shirenzigou during the first millennium BCE (Ling et al., 2013, 2016), although other research has documented millet consumption at other locations in the region (Liu and Reid, 2020). Lacustrine sediments from Lake Tuolekule (also known as Yanchi) and Barkol revealed a gradual contraction of lake areas during the late Holocene (An et al., 2011; Xue and Zhong, 2011). The lake area of Barkol, for example, decreased significantly around 4000 years ago with subsequent a shrinkage occurred around 2000 years ago (Wang et al., 2014).
The study sites in this research can be divided into two phases according to radiocarbon measurements (detailed below): The Bronze and Iron Age (Phase I) occupations date to c. 1300–800 BCE, and the later historical period (Phase II) dates to between c. 400 BCE and 300 CE. Most sites are dated to the earlier Phase I, including sites located on moraine hills (2,000–2,200 m.a.s.l.), such as Shirenzigou IV (Wang et al., 2009), Hongshankou (Ma et al., 2014a), and Kuola (Xi et al., 2020a), and lakeside sites such as Haiziyan at the shore of Lake Barkol (1,600 m.a.s.l.) (Ren et al., 2021) and Yanchigucheng near Lake Tuolekeule (1,950 m.a.s.l.) (Xi et al., 2020b). Those sites consist of stone-built dwellings with compartments enclosed by walls and roofed by wooden and mudbrick structures. Inside the dwellings, fireplaces, and associated ash dumps and trash pits were recovered, together with saddle-shaped grinding stones, stone pestles and hoes, and pottery jars, among other artefacts. Notably, storages with charred barley grains were recovered from Shirenzigou IV and Haiziyan, and archaeobotanical research established the importance of barley cultivation to subsistence (Ma et al., 2021; Ren et al., 2021; Tian et al., 2021b).
Only Shirenzigou III in Barkol is dated to Phase II. The site is fully excavated, including four dwellings (F1, F2, F3, and F4) located on hill slopes between 2,000 and 2,100 m.a.s.l. F1 is the largest stone-enclosed structure (85.5 m × 50 m) and could be the central building of the site (Ren, 2012; Tian et al., 2018). By contrast, other stone structures are smaller semi-subterranean buildings (<100 m2), surrounded by lower and thinner stone walls. This resonates with the hypothesis of growing nomadism during phase II, marked by changing burial customs and horse-riding evidence and unambiguous cultural exchange with the Altai region that has been discussed extensively (Ma, 2014; Ning et al., 2019; Li et al., 2021).
Materials and methods
Sampling and archaeobotanical analysis
During the excavations of Shirenzigou IV, Hongshankou, Haiziyan, and Shirenzigou III in Barkol, flotation samples (n = 350, 2,290 L) were collected from different features, including pits, ash dumps, brunt floors, and fireplaces within the dwellings (Supplementary Tables 2, 3). Additional soil samples from Yanchigucheng (n = 3, 17 L) and Kuola (n = 3, 23 L) were acquired during a field survey and collected from cultural layers on exposed profiles. All samples were floated using the bucket flotation method described in Pearsall (2016) and Zhao (2010). Light fractions were collected and suspended through a 0.2 mm mesh. Archaeobotanical analysis was subsequently conducted in the Paleoethnobotany Laboratory at Northwest University, Xi’an, China. Analysis followed the laboratory procedure described in Zhao (2010), using seed atlases and taxonomic keys for identification (Wei, 1993; Zhang, 2000; Qiang, 2002; Liu et al., 2008; Ullrich, 2011), as well as the Flora of China online database1 for taxonomic nomenclature. Considering the varying quantities of samples and differences in contexts at each site, we use calculated density (the number of seed remains per litre of samples) as well as portions in our analysis (Marston, 2015).
Plant stable isotope analysis
We selected 36 charred barley grains from Shirenzigou VI (n = 10), Haiziyan (n = 10), Yanchigucheng (n = 4), Kuola (n = 2), and Shirenzigou III (n = 10) (each sample consistes of a single charred barley grain). We did not include millet in the isotope analysis as the very small number of millet remains (two broomcorn and one foxtail millet from Shirenzigou III only) precluded the possibility of sampling. After removing visible surface contaminants by gentle scrapings, plant samples were pre-treated with 0.5 M HCl acid for 30 min at 80°C as described by Vaiglova et al. (2014). After freeze-drying, samples were weighed and transferred into tin capsules. Subsequent stable carbon and nitrogen isotopic values were measured using an elemental analyser coupled to a Delta V Plus continuous flow isotope ratio mass spectrometer located at the Biogeochemistry Lab, Washington University in St. Louis. δ13C and δ15N values were calibrated relative to VPDB and AIR, respectively, using international standard of USGS 40 (δ13C = −26.4‰, δ15N = −4.5‰) and USGS 41a (δ13C = +36.6‰, δ15N = −4.5‰) and two in-house standards acetanilide (δ13C = −29.5‰, δ15N = +47.6‰) and BR millet (Bob’s Red Mill millet flour; δ13C = −13.2‰, δ15N = +3.3‰). Following Szpak et al. (2017)’s procedure, the variability of the calibration standards (Ssrm) was determined to be 0.09‰ for δ13C and 0.17‰ for δ15N. According to the repeated measurements of international calibration standards, check standards, and sample replicates, the precision [u(Rw)] was determined to be 0.18‰ for δ13C and 0.43‰ for δ15N, respectively; accuracy [u(bias)] was determined to be 0.11‰ for δ13C and 0.24‰ for δ15N based on the difference between the observed and defaulted δ values of the check standards and standard deviations of them. Overall, the standard uncertainty for δ13C is 0.20‰ and for δ15N is 0.49‰.
Considering the fluctuation of atmospheric δ13C values through times, we calculated Δ13Cplant–air values to infer differences in the water status of crops from the study sites (Farquhar et al., 1989). The Δ13C values were calculated from the determined δ13Cplant values and a δ13Cair value approximated by the AIRCO2_LOESS system (Ferrio et al., 2005), using the equation defined by Farquhar et al. (1989):
Accelerator mass spectrometry radiocarbon measurements
We selected fourteen samples, including charred grains, animal bones, and wood charcoal, for radiocarbon (14C) analyses at Beta-Analytic based in Miami, USA (n = 11) and Radiocarbon Accelerator Laboratory, Peking University, China (n = 3). The pretreatment methods were similar, with a standard acid-base-acid (ABA) chemical pretreatment method followed by combustion and graphitisation prior to accelerator mass spectrometry (AMS) analyses. All radiocarbon dates were calibrated in OxCal 4.4 (Ramsey, 2017) using the IntCal20 calibration curve (Reimer et al., 2020).
Results and analytical frameworks
In this section, we shall announce the results of the three methodological strains: AMS radiocarbon measurement, archaeobotany, and stable isotope analysis. We shall also inform the readers about the analytical framework applied in evaluating the results.
Radiocarbon results and chronology
Accelerator mass spectrometry radiocarbon results are listed in Figure 1 and Supplementary Table 1. Most samples are dated (calibrated) from approximately 1300 BCE to the first centuries CE (Supplementary Table 1). Accordingly, study sites can be assigned into two phases, confirming observations based on uncovered archaeological remains. The early phase (Phase I) associated with the Bronze and Iron Ages are dated between ca. 1300 to 800 BCE, and the late phase (Phase II) is dated between ca. 400 BCE and 300 CE.
Figure 1. Accelerator mass spectrometry (AMS) radiocarbon results generated by this study together with previously published 14C measurements as detailed in Supplementary Table 1. These results confirm the two occupational phases of the eastern Tianshan complexity (ETC) sites.
Archaeobotanical assemblages
Results of archaeobotanical analyses show that barley (likely six row naked form, Hordeum vulgare subsp. vulgare L.) was the dominant crop (10,594 grains out of 17,064 total seeds remain count) of all sites in this study. Only two charred broomcorn and one foxtail millet grain are recovered from Shirenzigou III (Tian et al., 2018). Seeds of weedy taxa were also recovered in considerable quantities (n = 4,798) (Supplementary Tables 2, 3 and Figures 2, 3). Wood charcoal is abundant in the assemblages. Weed taxa can be grouped into two categories. First, arable weeds common for agricultural fields such as Setaria viridis, Avena fatua, Chenopodium album, Medicago sp., Silene conoidea, Lepyrodiclis holosteoides, Vaccaria segetalis, Thlaspi arvense, and Convolvulus sp. (n = 2,779). The second group includes weed taxa that are adaptive to well-watered soils, notably, Galium sp., Polygonum sp., Medicago sp., and Carex sp. (n = 2,019) may reflect well-watered soils for barley cultivation (Wei, 1993; Qiang, 2002; Ullrich, 2011; Motuzaite Matuzeviciute et al., 2020). Other plant taxa, for example, Astragalus sp., Melilotus sp., Corispermum sp., Artemisia sp., Nitraria sp., etc., have also been recovered in these assemblages.
Figure 2. Charred grains recovered from excavated eastern Tianshan complexity (ETC) sites. (A) Hordeum vulgare (ventral side) and (B) Hordeum vulgare (dorsal side) from barley cache at Shirenzigou IV F7; (C) Setaria italica from Shirenzigou III F1; (D) Panicum miliaceum from Shirenzigou III F1; (E) Setaria viridis from Shirenzigou III F1; (F) Lepyrodiclis holosteoides from Shirenzigou III F1; (G) Chenopodium album from Shirenzigou III F1; (H) Thlaspi arvense from Shirenzigou III F4; (I) Silene conoidea from Shirenzigou III F1; (J) Avena fatua from Haiziyan; (K) Convolvulus sp. from Shirenzigou III F1; (L) Vaccaria segetalis from Haiziyan; (M) Carex sp. from Shirenzigou IV F7; (N) Galium sp. from Shirenzigou III F1; (O) Polygonum sp. from Haiziyan; (P) Medicago sp. from Shirenzigou III F1.
Figure 3. The density (A) and proportion (B) of charred barley (purple), weeds (including arable weeds and water-demanding weeds, green), and water-demanding weeds (yellow) recovered from eastern Tianshan complexity (ETC) study sites.
Plant isotope compositions
Carbon and nitrogen isotope compositions (δ13C and δ15N values) measured as well as converted Δ13Cplant–air values are reported in Supplementary Table 4. We follow Wallace et al. (2013)’s interpretative framework of the barley isotope composition based on modern experiments in the semi-arid Eastern Mediterranean to infer the water status of past crops. The original thresholds between “poorly watered” and “moderately watered” (17‰ for barley), and between “moderately watered” and “well-watered” (18.5‰) are refer as “Optimal Watering Threshold” (OWT) and “Superfluous Watering Threshold” (SWT), respectively. In order to assess the manuring condition of archaeological barley, we established a local vegetation baseline of c. 2.8‰ by using average δ15N values of archaeological herbivores’ collagen and subtract 4.5‰, the average offset between contiguous trophic levels (DeNiro and Epstein, 1981; Minagawa and Wada, 1984; Ling et al., 2016).
Figure 4A displays the variances of Δ13C values of different sites. A considerable number of barley samples (n = 24) are below OWT—Shirenzigou III (n = 6) and Haiziyan (n = 8)—signifying water stress during plant growth. Shirenzigou IV and Kuola samples display a much better water condition as most samples are above the OWT threshold (n = 8 for Shirenzigou IV and n = 2 for Kuola). The altitude effect on the C3 plants may apply due to changes in atmospheric pressure and high carboxylation efficiency at high altitudes leading to foliar 13C depletion (Körner and Diemer, 1987; Morecroft and Woodward, 1990). However, the pattern is largely driven by taxonomic combinations rather than foliar δ13C variations within single species (Szpak et al., 2013). Currently, there are no available data to evaluate the altitude effect of barley. Given the ETC sites are located around 2,000 m.a.s.l. We estimate the OWT and SWT thresholds would be slightly higher than what Wallace and colleagues established originally, making the observed “water stress” more pronounced. Since our interpretations are based on large relative differences between sites at similar evaluations, the altitude effect is not deemed meaningful.
Figure 4. (A) Δ13C and (B) δ15N values measured in charred barley grains from eastern Tianshan complexity (ETC) sites. Dashed lines represent the threshold of “Optimal Water Status” (OWT) and “Superfluous Water Status” (SWT) for barley (Wallace et al., 2013; Vaiglova et al., 2020), and the estimated δ15N baseline of local vegetation using average δ15N values of archaeological herbivores (Ling et al., 2016).
Turning to δ15N values, Figure 4B demonstrates that all samples yield significantly higher δ15N values than the local vegetation baseline. The large variations of δ15N values, both between and within sites, suggests a variably but universally favorable nutrient supply of barley cultivation in this region.
Discussion
We will next consider the temporal and spatial patterns of barley cultivation in eastern Tianshan. We will first focus on changes in cultivation strategies through time. This is most clear at Shirenzigou. With its long-term occupational history, different sites at Shirenzigou offer insights into changes in subsistence strategies between the prehistoric and historical periods (Phase I and II). Secondly, as study sites are situated in markedly different environments, we shall explore the potential distinctions at the hill sites (Shirenzigou, Hongshankou, and Kuola) compared to areas near the lakes (Haiziyan and Yanchigucheng).
Changing cultivation strategies between 1300 BCE–300 CE
The isotopic results show a clear difference in water management and potentially soil preparation between the early and late occupational phases of Shirenzigou. The site has been central to archaeological investigations as one of the largest settlements in eastern Tianshan. Previous research revealed an increase in nomadic lifeways from the late Bronze Age to the historical period (Ma, 2014; Ma et al., 2014b). The two localities under investigation, Shirenzigou IV (SRZIV) and Shirenzigou III (SRZIII) were occupied during Phase I (1300–800 BCE) and Phase II (400 BCE–300 CE), respectively (Figure 1). This provides an opportunity to explore changes through time.
Most of the barley grains from the Bronze Age SRZIV display Δ13C values above the OWT (Figure 5A), indicating that their growth was not limited by water availability. Barley is adapted to arid environments since it is drought-resistant and is the pivot crop of rainfed agriculture in southwestern Asia (e.g., Riehl, 2012). In the eastern Tianshan, where the climate is dry, no evidence for enhanced precipitation is attested in palaeoclimatic records for the investigated period (Xue and Zhong, 2011; Wang et al., 2014). For this reason, the high-water status of barley crop observed isotopically at Shirenzigou could not be achieved without artificial water management. This doesn’t necessarily mean channel irrigation (such evidence is absent archaeologically). Barley could be strategically sown in areas with high water retention close to runoffs of snowmelts perhaps with small ditches to diffuse water or more sophisticated management of snowmelts at higher elevations with human or animal labour. In contrast, barley remains from the historical SRZIII show unambiguous water stress, with the plant development being limited by water availability (as Δ13C values below the OWT). This suggests that barley crops in the later occupational phase were either less watered or sown at drier locations than in the Bronze Age.
Figure 5. (A) Δ13C and (B) δ15N values of archaeological barley recovered from Shirenzigou IV of Phase I and Shirenzigou III of Phase II; (C) proportion of barley, arable weeds (excluding water-demanding weeds) and water-demanding weeds in the archaeobotanical assemblages of Shirenzigou IV (Phase I) and Shirenzigou III (Phase II).
A similar diachronic change, though less clear, is seen in barley δ15N values at SRZIV and SRZIII (Figure 5B). Crop δ15N values reflect nitrogen isotope compositions of the soil in which they are grown. The system is driven by nitrogen cycle openness and mycorrhizal associations that appear to be the dominant factors controlling nitrogen isotope variation in plants and soils (Szpak, 2014). In this context, soil 15N enrichment can be caused by a host of both environmental factors (e.g., nutrient status, salinity, topography, water saturation, etc.) and anthropogenic factors (e.g., Vaiglova et al., 2020). The latter has received much archaeological attention concerning application of animal manure has been found to increase the δ15N values of soil and crops significantly (e.g., Bogaard et al., 2013). Barley δ15N values from SRIV and SRZIII display wide ranges. This observation suggests that crops were cultivated in variable soil conditions in both periods, which could result from several natural causes and the manuring effect combined (Groffman and Hanson, 1997; Finlay et al., 2008; Fraser et al., 2011; Styring et al., 2016). As barley grains from both sites display δ15N values significantly higher than the local vegetation baseline, anthropogenic activities such as the application of animal manure or organic waste likely contributed to the soil 15N enrichment. Otherwise, we expect the nitrogen isotope values to be closer to the baseline. Although natural causes (e.g., soil denitrification or aridity) cannot be excluded, the difference between SRZIV and SRZIII might indicate that barley fields were less manured during Phase II compared to Phase I. This resonates with the plant water status discussed above, suggesting less labour input/care in agricultural fields in a later phase.
The archaeobotanical results seem to support such a conclusion further. Tian et al. (2021a) observed a change in weed abundance through time, with a lower weed/barley ratio in the Bronze Age and a higher weed/barley ratio in the historical period. While the density of weeds decreases over time (Figure 3), their proportional input increases from the early to the late period (Figure 5C). The proportional input of water-demanding weeds decreased over time echoing the observation of less irrigation as suggested by carbon isotope compositions. It should be noted that this observed pattern could be driven by crop processing with increases in post-harvest and pre-storage labour biasing toward higher weed/barley ratios and higher weed proportion as the grains were removed for consumption (Stevens, 2003) or due to the scale of processing and/or producer/consumer status of the site (Jones, 1985; Van der Veen and Jones, 2006). Future research considering assemblage formation processes in contexts will guide our understanding. A definitive conclusion cannot be reached with the current available evidence. However, an alternative interpretation of the high weed proportion and weed/barley ratio at historical SRZIII is less pre-harvest labour, such that fewer weeding activities led to higher weed ubiquity in the assemblage. This possibility supports the isotopic results, suggesting less labour investment in field conditions before the harvest with regards to irrigation, manuring, and weeding. The decrease of water-demanding weed taxa (Figure 5C) over time could also be understood in the context of reduced irrigation.
In modern-day eastern Tianshan, spring-sown wheat is an important crop cultivated on the alluvial landscape with channel irrigation. This is an extensive system. According to local residences, each household (2–4 family members) is responsible for more than 10,000 square metres of arable land. As such, labour deficiency restrains intensive field activities such as weeding, and many weeds can be observed in the fields during the growing season.
Field management system between lakeside and foothill site
As discussed, the study sites are situated in markedly different environments, including areas on the piedmont hills such as the Shirenzigou group, Hongshankou, and Kuola and sites near the lakeshores Haiziyan and Yanchigucheng (Figures 6, 7). Barley grains from lakeside sites display significantly lower Δ13C values (below OWT) than barley from the hill sites (Figure 6). The lower Δ13C values at the lakeside sites could be the results of high soil salinity near the saline Lake Barkol and Tuolekule. Controlled experiments show that increases in soil salinity correspond to reducing the leaf stomatal conductance, resulting in lower Δ13C values and yield (Isla et al., 1998). In addition to the high soil salinity, lakeside sites are located on flat terrains at distances to the foothill runoffs or springs, making irrigation more difficult. Combined, agricultural fields near the lakes were likely more water-stressed.
Figure 6. (A) The typical landscapes of Bronze Age lakeside sites (the excavation area of Haiziyan in 2017) and (B) piedmont sites (the piedmont landscape of Shirenzigou) in the eastern Tianshan region; (C) Δ13C and (D) δ15N values of archaeological barley from lakeside (Haiziyan, Yanchigucheng) and piedmont sites (Shirenzigou IV, Hongshankou, Kuola); (E) proportion of barley, arable weeds (excluding water-demanding weeds) and water-demanding weeds uncovered from lakeside and piedmont sites.
Figure 7. The locations of Shirenzigou IV (SRZ IV, 2,140 m.a.s.l.), Shirenzigou III (SRZ III, 2,200 m.a.s.l.), Haiziyan (HZY, 1,600 m.a.s.l.), Yanchigucheng (YCGC, 2,020 m.a.s.l.), Kuola (KL, 2,200 m.a.s.l.) and Hongshankou (HSK, 2,200 m.a.s.l.) in the eastern Tianshan region. The lowest and highest elevations of Barkol are at Lake Barkol (1,581 m.a.s.l.) and Yueyashan (4,308 m.a.s.l.), respectively. The lowest and highest elevations of Yiwu are at Lake Tuolekule (1,890 m.a.s.l.) and Tuomuerti (4,886 m.a.s.l.), respectively. Maps were generated using ArcMap v. 10.2 with data credit to National Park Service Natural Earth physical map (public domain).
In contrast, barley samples from piedmont hill sites yield Δ13C values that are mostly above OWT (with two exceptions). Some samples are even above SWT (Figure 6). This suggests barley crops from those sites were well-watered. These results enable an observation of distinct cultivation strategies concerning water regimes (including both irrigation and strategic planting in soils with high water retention) at the piedmonts and lakeshores.
This observation resonates with the archaeobotanical results, suggesting a difference between the environmental conditions of barley fields by the lakeside and on the foothills. The densities of barley grains, weeds, and water-demanding taxa in the piedmont sites (Shirenzigou IV F7, Kuola, and Hongshankou) all outnumber their counterparts on the lakeside (Haiziyan and Yanchigucheng, Figure 3). The proportions show a different pattern with arable and water-demanding weeds are lower at the lakeside in comparison to the piedmont sites. The percentage data are subject to further scrutiny considering assemblage formation as grain/weed ratio typically reflect crop processing activities rather than field condition (e.g., Stevens, 2003; Van der Veen and Jones, 2006). That issue aside, some arable weeds are informative to crop growing conditions and can provide complementary information (e.g., Bogaard et al., 2016). It is noticeable that water-demanding taxa such as Carex, Galium, Polygonum, and Medicago are more ubiquitous in the piedmont assemblages compared to the lakeside assemblages (Figures 3, 6E). Carex, in particular, is sometimes seen as an indicator of field irrigation due to its adaptation to open and well-watered soils (Rühl et al., 2015; Motuzaite Matuzeviciute et al., 2020).
The isotopic and archaeobotanical evidence, therefore, suggests two distinct barley cultivation strategies co-existing in eastern Tianshan. Barley crops were grown in stressful conditions near the saline lakes, while crops were relatively well-looked after in the piedmont zone. As mentioned, archaeologically visible irrigation channels are absent in this region, but water management can be achieved by small-scale modifications or strategic planting without leaving traits in projected archaeology. Ethnographic records in Tuva, for example, demonstrated the application of simple channels for irrigating fields among southern Siberian pastoralists (Vainshtein, 1980). Some piedmont fields in modern-day Barkol are irrigated by simple ditches bringing glacial water from the top of the mountain.
Conclusion
The data presented in this paper allows two inferences. First, our results indicate that different cultivation strategies were employed at different ecosystems. Plant isotope results show, unequivocally, water management at hill sites to meet barley’s growing demand in an otherwise arid environment. However, barley fields near the lake shores were stressed due to soil salinity.
The second inference concerns changing cultivation strategies over time. Stable carbon and nitrogen isotope compositions and archaeobotanical results combined suggest a shift in barley cultivation from the Bronze Age to the historical period manifested by the difference between SRZIV and SRZIII. The former likely represents a labour-intensive system with water management, application of animal manure or midden for soil improvement, and potentially frequent weeding activities during 1300 BCE and 800 BCE. The latter is associated with an extensive cultivation strategy with diminished input in irrigation, manuring, and weeding between 400 BCE and 300 CE, corresponding to the development of nomadic lifeways during this time that has been discussed elsewhere (Ma, 2014; Ning et al., 2019; Li et al., 2021).
These findings shed light on an extensification trend of barley cultivation from the late Bronze Age to the historical period. During this later period, with communities increasingly relying on pastoral activities as the main subsistent strategies, intensive cultivation with high labour input per area seen in the Bronze Age was partially diminished. The latter was given way to an extensive approach in which productivity is driven by the size and expansiveness of herds, pastures and cultivation, a trend that has been speculated by scholars but we provided tangible evidence here. Of course, there are limitations. We began our interpretation with carbon isotope evidence that is quite clear; the nitrogen isotope compositions also present clear trends but are more enigmatic in distinguishing natural and anthropogenetic drivers; the archaeobotanical data, however, has more than one interpretation. This is where the ambiguity begins, and we draw attention to the value of context-based assemblage formation process in future research.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
XLiu, DT, and YS conceptualised the research with MMR’s contribution and wrote the draft which was read and further edited by all authors. YS, XLing, and XLiu conducted the isotopic analysis. DT and ZZ carried out macrofossil analysis and sample selection. TX, MRen, JM, and JW led excavations and archaeological surveys. All authors contributed to the article and approved the submitted version.
Funding
DT, TX, MRen, and JM are grateful to the program of “Youth Innovation Team of Shaanxi Universities”. DT and the archaeobotanical analysis were supported by the National Social Science Foundation of China (19CKC031). DT would like to acknowledge the “111 Project” of Northwest University (D18004) in supporting the open-access publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.966366/full#supplementary-material
References
Alagich, R., Gardeisen, A., Alonso, N., Rovira, N., and Bogaard, A. (2018). Using stable isotopes and functional weed ecology to explore social differences in early urban contexts: The case of Lattara in Mediterranean France. J. Archaeol. Sci. 93, 135–149. doi: 10.1016/j.jas.2018.03.006
An, C. B., Zhao, J., Tao, S., Lv, Y., Dong, W., Li, H., et al. (2011). Dust variation recorded by lacustrine sediments from arid Central Asia since~ 15 cal ka BP and its implication for atmospheric circulation. Quat. Res. 75, 566–573. doi: 10.1016/j.yqres.2010.12.015
Araus, J. L., Febrero, A., Buxó, R., Rodrıguez-Ariza, M. O., Molina, F., Camalich, M. D., et al. (1997). Identification of ancient irrigation practices based on the carbon isotope discrimination of plant seeds: A case study from the South-East Iberian Peninsula. J. Archaeol. Sci. 24, 729–740. doi: 10.1006/jasc.1997.0154
Bogaard, A., Fraser, R., Heaton, T. H., Wallace, M., Vaiglova, P., Charles, M., et al. (2013). Crop manuring and intensive land management by Europe’s first farmers. Proc. Natl. Acad. Sci. U.S.A. 110, 12589–12594. doi: 10.1073/pnas.1305918110
Bogaard, A., Heaton, T. H., Poulton, P., and Merbach, I. (2007). The impact of manuring on nitrogen isotope ratios in cereals: Archaeological implications for reconstruction of diet and crop management practices. J. Archaeol. Sci. 34, 335–343. doi: 10.1016/j.jas.2006.04.009
Bogaard, A., Hodgson, J., Nitsch, E., Jones, G., Styring, A., Diffey, C., et al. (2016). Combining functional weed ecology and crop stable isotope ratios to identify cultivation intensity: A comparison of cereal production regimes in Haute Provence, France and Asturias, Spain. Veg. Hist. Archaeobot. 25, 57–73. doi: 10.1007/s00334-015-0524-0
Bogaard, A., Styring, A., Ater, M., Hmimsa, Y., Green, L., Stroud, E., et al. (2018). From traditional farming in Morocco to Early Urban Agroecology in Northern Mesopotamia: Combining present-day arable weed surveys and crop isotope analysis to reconstruct past Agrosystems in (Semi-)arid regions. Environ. Archaeol. 23, 303–322. doi: 10.1080/14614103.2016.1261217
DeNiro, M. J., and Epstein, S. (1981). Influence of diet on the distribution of nitrogen isotopes in animals. Geochim. Cosmochim. Acta 45, 341–351. doi: 10.1016/0016-7037(81)90244-1
Farquhar, G. D., Ehleringer, J. R., and Hubick, K. T. (1989). Carbon isotope discrimination and photosynthesis. Annu. Rev. Plant Biol. 40, 503–537. doi: 10.1146/annurev.pp.40.060189.002443
Ferrio, J. P., Araus, J. L., Buxó, R., Voltas, J., and Bort, J. (2005). Water management practices and climate in ancient agriculture: Inferences from the stable isotope composition of archaeobotanical remains. Veg. Hist. Archaeobot. 14, 510–517. doi: 10.1007/s00334-005-0062-2
Finlay, J. C., Kendall, C., Michener, R., and Lajtha, K. (2008). “Stable isotope tracing of temporal and spatial variability in organic matter sources to freshwater ecosystems,” in Stable isotopes in ecology and environmental science, eds R. Michener and K. Lajtha (Hoboken, NJ: Wiley), 283–324. doi: 10.1002/9780470691854.ch10
Flohr, P., Müldner, G., and Jenkins, E. (2011). Carbon stable isotope analysis of cereal remains as a way to reconstruct water availability: Preliminary results. Wat. Hist. 3:121. doi: 10.1007/s12685-011-0036-5
Frachetti, M. D. (2012). Multiregional emergence of mobile pastoralism and nonuniform institutional complexity across Eurasia. Curr. Anthropo. 53, 2–38. doi: 10.1086/663692
Fraser, R. A., Bogaard, A., Heaton, T., Charles, M., Jones, G., Christensen, B. T., et al. (2011). Manuring and stable nitrogen isotope ratios in cereals and pulses: Towards a new archaeobotanical approach to the inference of land use and dietary practices. J. Archaeol. Sci. 38, 2790–2804. doi: 10.1016/j.jas.2011.06.024
Groffman, P. M., and Hanson, G. C. (1997). Wetland denitrification: Influence of site quality and relationships with wetland delineation protocols. Soil Sci. Soc. Am. J. 61, 323–329. doi: 10.2136/sssaj1997.03615995006100010047x
Harvey, E. L., and Fuller, D. Q. (2005). Investigating crop processing using phytolith analysis: The example of rice and millets. J. Archaeol. Sci. 32, 739–752. doi: 10.1016/j.jas.2004.12.010
Hillman, G. C. (1984). “Reconstructing crop husbandry practices from charred remains of crops,” in Farming practice in British prehistory, ed. R. Mercer (Edinburgh: Edinburgh University Press), 123–162.
Isla, R., Aragüés, R., and Royo, A. (1998). Validity of various physiological traits as screening criteria for salt tolerance in barley. Field Crops Res. 58, 97–107. doi: 10.1016/S0378-4290(98)00088-4
Jones, M., Hunt, H., Lightfoot, E., Lister, D., Liu, X., and Motuzaite-Matuzeviciute, G. (2011). Food globalization in prehistory. World Archaeol. 43, 665–675. doi: 10.1080/00438243.2011.624764
Jones, M. K. (1985). “Archaeobotany beyond subsistence reconstruction,” in Beyond Domestication in Prehistoric Europe, eds G. W. Barker and C. Gamble (New York, NY: Academic Press), 107–128.
Körner, C., and Diemer, M. (1987). In situ photosynthetic responses to light, temperature and carbon dioxide in herbaceous plants from low and high altitude. Funct. Ecol. 1, 179–194. doi: 10.2307/2389420
Li, H., Sun, Y., Yang, Y., Cui, Y., Ren, L., Li, H., et al. (in press). Distinct crop management strategies of first wheat and barley cultivators in north China. Antiquity
Li, Y., Ma, J., Zhang, C., Liu, H., Zong, T., Chen, T., et al. (2020). Worked bone production in ancient pastoral societies in China: A case study of the Shirenzigou site in Balikun, Xinjiang. Quaternary Sci. 40, 331–342.
Li, Y., Zhang, C., Taylor, W. T. T., Chen, L., Flad, R. K., Boivin, N., et al. (2021). Early evidence for mounted horseback riding in northwest China. Proc. Natl. Acad. Sci. U.S.A. 117, 29569–29576. doi: 10.1073/pnas.2004360117
Ling, X., Chen, X., Wang, J., Chen, L., Ma, J., Ren, M., et al. (2013). Carbon and nitrogen isotopic analysis of human bones excavated from the Dongheigou site, Ba Likun County, Xinjiang Autonomous Region. Acta Anthropo. Sin. 32, 219–225. In Chinese with English Abstract.
Ling, X., Lan, D., Chen, X., Ma, J., Wang, J., and You, Y. (2016). Carbon and nitrogen isotopic analysis on animal bones excavated from Dongheigou site in Barkol County, Xinjiang Autonomous Region. West. Archaeol. 11, 289–299. In Chinese with English Abstract.
Liu, C., Jin, G., and Kong, Z. (eds) (2008). Archaeobotany: Research on seeds and fruit. Beijing: Science Press. In Chinese.
Liu, X., Jones, P. J., Matuzeviciute, G. M., Hunt, H. V., Lister, D. L., An, T., et al. (2019). From ecological opportunism to multi-cropping: Mapping food globalisation in prehistory. Quat. Sci. Rev. 206, 21–28. doi: 10.1016/j.quascirev.2018.12.017
Liu, X., Lister, D. L., Zhao, Z., Petrie, C. A., Zeng, X., Jones, P. J., et al. (2017). Journey to the east: Diverse routes and variable flowering times for wheat and barley en route to prehistoric China. PLoS One 12:e0187405. doi: 10.1371/journal.pone.0187405
Liu, X., and Reid, R. E. B. (2020). The prehistoric roots of Chinese cuisines: mapping staple food systems of China, 6000 BC–220 AD. PLoS One 15:e0240930. doi: 10.1371/journal.pone.0240930
Ma, J. (2014). “The rise of nomads on the Barkol Steppe and their cultural interflow with the Altai Region during the 1st millennium BCE,” in Eurasian studies II, eds T. Yu and J. Li (Sydney: Asia Publishing Nexus), 29–37.
Ma, J., Xi, T., Ren, M., and Wang, J. (2014a). Field survey report of Hongshankou Site, Barkol, Xinjiang in 2008. Cult. Relics. 7, 17–30. In Chinese.
Ma, J., Wang, J., Zhao, H., Xi, T., Yahepu, J., and Jiang, X. (2014b). Excavation report on F2 of Shirenzigou Site, Barkol, Xinjiang in 2009. Archaeol. Cult. Relics. 5, 25–36. In Chinese with English abstract.
Ma, Z., Liu, S., Ren, M., Huan, X., Xi, T., Wang, J., et al. (2021). Analysis of grain storage pit of the Shirenzigou site in Balikun County, east Tianshan area, Xinjiang. Quat. Sci. 41, 214–223. (In Chinese with English abstract). doi: 10.11928/j.issn.1001-7410.2021.01.19
Marston, J. M. (2015). “Ratios and simple statistics in paleoethnobotanical analysis: Data exploration and hypothesis testing,” in Method and theory in paleoethnobotany, eds J. M. Marston, J. d’Alpoim Guedes, and C. Warinner (Boulder, CO: University Press of Colorado), 163–179. doi: 10.5876/9781607323167.c009
Minagawa, M., and Wada, E. (1984). Stepwise enrichment of 15N along food chains: Further evidence and the relation between δ15N and animal age. Geochim. Cosmochim. Acta 48, 1135–1140.
Morecroft, M. D., and Woodward, F. I. (1990). Experimental investigations on the environmental determination of δ13C at different altitudes. J. Exp. Bot. 41, 1303–1308.
Motuzaite Matuzeviciute, G., Mir-Makhamad, B., and Tabaldiev, K. R. (2020). The first comprehensive archaeobotanical analysis of prehistoric agriculture in Kyrgyzstan. Veg. Hist. Archaeobot. 30, 743–758. doi: 10.1007/s00334-021-00827-0
Ning, C., Wang, C. C., Gao, S., Yang, Y., Zhang, X., Wu, X., et al. (2019). Ancient genomes reveal Yamnaya-related ancestry and a potential source of Indo-European speakers in Iron Age Tianshan. Curr. Biol. 29, 2526–2532. doi: 10.1016/j.cub.2019.06.044
Qiang, S. (2002). Weed diversity of arable land in China. Korean J. Weed Sci. 22, 87–198. doi: 10.7717/peerj.13583
Ramsey, C. B. (2017). Methods for summarizing radiocarbon datasets. Radiocarbon 59, 1809–1833. doi: 10.1017/RDC.2017.108
Reimer, P. J., Austin, W. E., Bard, E., Bayliss, A., Blackwell, P. G., Ramsey, C. B., et al. (2020). The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). Radiocarbon 62, 725–757. doi: 10.1017/RDC.2020.41
Ren, M. (2012). Study on archaeological remains around the eastern Tianshan Mountain in first millennium BC. Dissertation. Xi’an: Northwest University. In Chinese with English abstract.
Ren, M., Ma, J., Xi, T., Wang, J., Li, W., Tian, Y., et al. (2021). Achievements and reflection of archaeological excavation at Hanziyan Site in Balikun, Xinjiang. West. Reg. Stud. 4, 141–148. (In Chinese with English abstract). doi: 10.16363/j.cnki.xyyj.2021.04.015
Riehl, S. (2012). Variability in ancient Near Eastern environmental and agricultural development. J. Arid Environ. 86, 113–121.
Rühl, L., Herbig, C., and Stobbe, A. (2015). Archaeobotanical analysis of plant use at Kamennyi Ambar, a Bronze Age fortified settlement of the Sintashta culture in the southern Trans-Urals steppe, Russia. Veg. Hist. Archaeobot. 24, 413–426. doi: 10.1007/s00334-014-0506-7
Spengler, R. N. (2015). Agriculture in the central Asian bronze age. J. World Prehist. 28, 215–253. doi: 10.1007/s10963-015-9087-3
Stevens, C. J. (2003). An investigation of agricultural consumption and production models for prehistoric and Roman Britain. Environ. Archaeol. 8, 61–76. doi: 10.1179/env.2003.8.1.61
Styring, A. K., Ater, M., Hmimsa, Y., Fraser, R., Miller, H., Neef, R., et al. (2016). Disentangling the effect of farming practice from aridity on crop stable isotope values: A present-day model from Morocco and its application to early farming sites in the eastern Mediterranean. Anthr. Rev. 3, 2–22. doi: 10.1177/2053019616630762
Styring, A. K., Charles, M., Fantone, F., Hald, M. M., McMahon, A., Meadow, R. H., et al. (2017). Isotope evidence for agricultural extensification reveals how the world’s first cities were fed. Nat. Plants 3, 1–11. doi: 10.1038/nplants.2017.76
Szpak, P. (2014). Complexities of nitrogen isotope biogeochemistry in plant-soil systems: Implications for the study of ancient agricultural and animal management practices. Front. Plant Sci. 5:288. doi: 10.3389/fpls.2014.00288
Szpak, P., Metcalfe, J. Z., and Macdonald, R. A. (2017). Best practices for calibrating and reporting stable isotope measurements in archaeology. J. Archaeol. Sci. Rep. 13, 609–616. doi: 10.1016/j.jasrep.2017.05.007
Szpak, P., White, C. D., Longstaffe, F. J., Millaire, J. F., and Vásquez Sánchez, V. F. (2013). Carbon and nitrogen isotopic survey of northern Peruvian plants: Baselines for paleodietary and paleoecological studies. PLoS One 8:e53763. doi: 10.1371/journal.pone.0053763
Tian, D., Ma, J., Wang, J., Pilgram, T., Zhao, Z., and Liu, X. (2018). Cultivation of naked barley by early Iron Age agro-pastoralists in Xinjiang, China. Environ. Archaeol. 23, 416–425. doi: 10.1080/14614103.2017.1415121
Tian, D., Xi, T., Ren, M., Ma, J., Wang, J., and Zhao, Z. (2021a). Agricultural production in Barkol, Xinjiang in the first millennium BC: From the perspective of farmland weeds. West. Reg. Stud. 3, 72–82. In Chinese with English abstract.
Tian, D., Xi, T., Ren, M., Ma, J., Wang, J., and Zhao, Z. (2021b). Barley Production in Early Xinjiang:Evidence from archaeobotanical remains in the Shirenzigou Site, Northern Tianshan Mountains. Agricul. Hist. Chin. 40, 44–55. In Chinese with English abstract.
Vaiglova, P., Gardeisen, A., Buckley, M., Cavanagh, W., Renard, J., Lee-Thorp, J., et al. (2020). Further insight into Neolithic agricultural management at Kouphovouno, southern Greece: Expanding the isotopic approach. Archaeol. Anthropol. Sci. 12, 1–17. doi: 10.1007/s12520-019-00960-y
Vaiglova, P., Snoeck, C., Nitsch, E., Bogaard, A., and Lee-Thorp, J. (2014). Impact of contamination and pretreatment on stable carbon and nitrogen isotopic composition of charred plant remains. Rapid Commun. Mass Spectr. 28, 2497–2510. doi: 10.1002/rcm.7044
Vainshtein, S. (1980). Nomads of South Siberia: The pastoral economies of Tuva. Cambridge: Cambridge University Press.
Van der Veen, M., and Jones, G. (2006). A re-analysis of agricultural production and consumption: Implications for understanding the British Iron Age. Veg. Hist. Archaeobot. 15, 217–228. doi: 10.1007/s00334-006-0040-3
Wallace, M., Jones, G., Charles, M., Fraser, R., Halstead, P., Heaton, T. H., et al. (2013). Stable carbon isotope analysis as a direct means of inferring crop water status and water management practices. World Archaeol. 45, 388–409. doi: 10.1080/00438243.2013.821671
Wang, H., Yue, L., Li, J., and Yang, L. (2014). Changing of the Lake Area and records of climate and environment of Barkol Lake During Holocene. Acta Sedimentol. Sin. 32, 93–100. In Chinese. doi: 10.14027/j.cnki.cjxb.2014.01.011
Wang, J., Zhang, F., Ren, M., Yahepu, J., and Yu, J. (2009). Excavation report on Dongheigou Site in Barkol of Xinjiang, Year 2006–2007. Archaeology 1, 3–27. In Chinese with English abstract.
Wei, Q. (ed.) (1993). Xinjiang farmland weeds and control. Urumqi: Xinjiang Science and Health Publishing House. In Chinese.
Xi, T., Fu, Y., Maihemaiti, R., Ma, J., and Ren, M. (2020a). Field survey report of Kuola Site, Yiwu, Xinjiang, 2017-2018. Cult. Relics. 8, 17–30. In Chinese. doi: 10.13619/j.cnki.cn11-1532/k.2020.08.002
Xi, T., Wang, W., Maihemaiti, R., Yiliyazi, S., Ma, J., and Ren, M. (2020b). Field survey report of Yanchigucheng Site, Yiwu, Xinjiang in 2017. West. Archaeol. 19, 26–38. In Chinese with English abstract.
Xue, J., and Zhong, W. (2011). Holocene climate variation denoted by Barkol Lake sediments in northeastern Xinjiang and its possible linkage to the high and low latitude climates. Sci. China Earth Sci. 54, 603–614. doi: 10.1007/s11430-010-4111-z
Zhang, J. (ed.) (1993). Gazetteer of Barkol Kazakh county. Urumqi: Xinjiang University Press. (In Chinese).
Zhang, Z. (2000). Chinese colored weed illustrated book. Beijing: Institute for the Control of Agrochemicals, Ministry of Agriculture, P. R. China. In Chinese.
Zhao, Z. (2010). Paleoethnobotany: Theories, methods and practice. Beijing: Science Press. (In Chinese with English abstract).
Keywords: isotope, barley, eastern Tianshan, irrigation, manuring, archaeobotany
Citation: Tian D, Sun Y, Ritchey MM, Xi T, Ren M, Ma J, Wang J, Zhao Z, Ling X and Liu X (2022) Varying cultivation strategies in eastern Tianshan corresponded to growing pastoral lifeways between 1300 BCE and 300 CE. Front. Ecol. Evol. 10:966366. doi: 10.3389/fevo.2022.966366
Received: 10 June 2022; Accepted: 21 September 2022;
Published: 12 October 2022.
Edited by:
Anna Maria Mercuri, University of Modena and Reggio Emilia, ItalyReviewed by:
Gabriel Servera-Vives, University of Modena and Reggio Emilia, ItalyCristiano Vignola, Sapienza University of Rome, Italy
Copyright © 2022 Tian, Sun, Ritchey, Xi, Ren, Ma, Wang, Zhao, Ling and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinyi Liu, liuxinyi@wustl.edu; Xue Ling, lxue@nwu.edu.cn
†These authors have contributed equally to this work