Abstract
The origin of millet from Neolithic China has generally been accepted, but it remains unknown whether common millet (Panicum miliaceum) or foxtail millet (Setaria italica) was the first species domesticated. Nor do we know the timing of their domestication and their routes of dispersal. Here, we report the discovery of husk phytoliths and biomolecular components identifiable solely as common millet from newly excavated storage pits at the Neolithic Cishan site, China, dated to between ca. 10,300 and ca. 8,700 calibrated years before present (cal yr BP). After ca. 8,700 cal yr BP, the grain crops began to contain a small quantity of foxtail millet. Our research reveals that the common millet was the earliest dry farming crop in East Asia, which is probably attributed to its excellent resistance to drought.
Keywords: Holocene, origins of agriculture, phytoliths, Neolithic, Cishan
Foxtail millet (Setaria italica) and common millet (or broomcorn millet; Panicum miliaceum) were among the world's most important and ancient domesticated crops. They were staple foods in the semiarid regions of East Asia (China, Japan, Russia, India, and Korea) and even in the entire Eurasian continent before the popularity of rice and wheat (1–4), and are still important foods in these regions today (5, 6).
Thirty years ago, the world's oldest millet remains, dating to ca. 8,200 calibrated years before present (cal yr BP), were discovered at the Early Neolithic site of Cishan, northern China. The site contained >50,000 kg of grain crops stored in the storage pits (7–9). Until now, the importance of these findings has been constrained by limited taxonomic identification with regard to whether they are from foxtail millet (S. italica) or common millet (P. miliaceum), because the early reported S. italica identifications are not all accepted (4, 9–12). This article presents the phytoliths, biomolecular records, and new radiocarbon dating from newly excavated grain crop storage pits at the Cishan site. Large modern reference collections are used to compare and contrast microfossil morphology and biomolecular components in different millets and related grass species (13). The renewed investigations show that common millet agriculture arose independently in the semiarid regions of China by 10,000 cal yr BP. Our findings contribute to our knowledge of agricultural origins across the globe and have broader implications for understanding the development of human societies.
The Cishan site (36°34.511′ N, 114°06.720′ E) is located near the junction between the Loess Plateau and the North China Plain at an elevation of 260–270 m above sea level (Fig. 1). The archaeological site, containing a total of 88 storage pits with significant quantities (≈109 m3) of grain crop remains, was excavated from 1976 to 1978 (7, 8). Each storage pit included 0.3- to 2-m-thick grain crops, which were well preserved and found in situ in the 3- to 5-m-deep loess layer (9). All grain remains have been oxidized to ashes soon after they were exposed to air. Archaeological excavations also revealed the remains of houses and numerous millstones (Fig. S1), stone shovels, grind rollers, potteries, rich faunal remains, and plant assemblages including charred fruits of walnut (Juglans regia), hazel (Corylus heterophylla), and hackberry (Celtis bungeana) (7–9). Only 2 14C dates of charcoal from previously excavated H145 and H48 storage pits yielded uncalibrated ages of 7355 ± 100 yr BP and 7235 ± 105 yr BP, respectively (8). These remains represent the earliest evidence for the significant use of dry-farming crop plants in the human diet in East Asia. They also suggest that by this time agriculture had already been relatively well developed here.
Early identification assumed these grain crop remains to be foxtail millet (S. italica). This preliminary identification was mainly based on a characteristic of very small sizes of grain crop ash (no charred grains)—often <2–3 mm in length—that resembled the foxtail millets (Fig. S2), but without any spodogram evidence for the grain crop remains in the Cishan site (7–10). These have been reported as the world's oldest foxtail millet in the literature (6, 9, 10, 14). However, the millet identification has been questionable (4, 11, 12, 15), because the macro (ash) remains were too friable to be observed under the microscope. Furthermore, very little study has been conducted on the spodograms or phytoliths of modern millets, so no clear diagnostic feature has been used to distinguish foxtail millet from common millet (12, 16). Previous study has also considered at some length how the charred dehusked grains of various native millet species might have been systematically misidentified (3). Thus, questions remain regarding whether the Cishan grain crops are from foxtail millet or common millet, or both, because they cannot be distinguished from the ash vestige (11, 12).
Results
Phytoliths and biomolecular components have provided substantial empirical evidence demonstrating the considerable antiquity of food production and crop dispersals in many regions of the world (17–20). Based on our observation and statistics of the variation of anatomy and silicon structure patterns in the glumes, lemmas, and paleas occurring in modern cultivated millets and related grass species collected from different regions in China (13), we found that 5 key diagnostic characteristics in phytolith morphology could be used to distinguish between foxtail millet and common millet (13): (i) a cross-shaped type phytolith is formed in the lower lemma and glume of S. italica, whereas a bilobe-shaped type is formed in those of P. miliaceum; (ii) regularly arranged papillae on the surface of the upper lemma and paleas are peculiar to S. italica; (iii) the epidermal long cell walls are Ω-undulated in S. italica (Fig. 2H), and η-undulated in P. miliaceum (Fig. 2 B, D, and F); (iv) the endings structures of epidermal long cells are cross-wavy type in S. italica and cross-finger type in P. miliaceum; the R value (ratio of the width of endings interdigitation to the amplitude of undulations) is higher (0.79 ± 0.12; n = 3,303) in P. miliaceum than in S. italica (0.33 ± 0.11; n = 2,774); and (v) surface ridgy line sculpture of the upper lemma is peculiar to S. italica. These 5 diagnostic characteristics used together give the only reliable way of distinguishing foxtail millet from common millet when only powder remains are available. In addition, a species-specific identification of phytoliths is possible for S. italica and P. miliaceum because they have typically well-defined silica skeletons that are distinguishable from those in Panicum bisulcatum, Setaria viridis, and Setaria plicata, which have no such demonstrable patterns (13).
To determine the taxa of foxtail millet and common millet, we analyzed the phytoliths of 46 grain crop samples stored in 5 newly excavated storage pits (CS-I to CS-V) at the Cishan site, and 1 sample (BWG) preserved in a storage bottle from the Culture Museum of Cishan. These grain crop samples are dated between ca. 10,300 and ca. 7,500 cal yr BP based on new 14C dating measurements (Fig. 3, Table S1).
All 47 archaeological samples we analyzed contained abundant diagnostic husk phytoliths that can be divided into 2 groups according to their phytolith assemblages and 14C dating results. The first group, including 27 samples from CS-I, III, and IV, is dated between ca. 10,300 and ca. 8,700 cal yr BP. All of the husk phytoliths present are diagnostic of common millet based on their shapes, η patterns (Fig. 2 A, C, and E), and average R value (>0.7) (Fig. 2I). The second group, including 20 samples from CS-II, V, and BWG, is dated between ca. 8700 and ca. 7,500 cal yr BP. More than 97% of the husk phytoliths are also diagnostic of the common millet (97.2% for CS- II, n = 1,273; 97.5% for CS-V, n = 1,000; 99.6% for BWG, n = 1,000), but a small quantity (0.4–2.8%) of the husk phytoliths in the second group can be attributed to S. italica (Fig. 2 G and I).
Fig. 4 shows that the CS-I storage pit contains ≈1.5-m-thick grain crop remains composed of 3 prominent layers of lemma and palea from common millets alternating with 3 glume plus reed layers. This storage manner indicates that prehistoric humans had known how to preserve large volumes of grains in secure storage. They did so by digging deep storage pits in the dry loess strata and by covering the floors with thick mats of millet glumes and reed (Phragmites australis) leaves.
Identification of biomolecular components was mainly based on our modern reference collection. A recent study has used the biomolecular components of P. miliaceum, the only miliacin-exclusive producer reported, as a basis for identification (20). However, because of the lack of comparable data derived from S. italica, previous investigators are still not clear how to use biomarkers to distinguish P. miliaceum from S. italica. In this study, we examined biomolecular components in 6 samples of modern Paniceae (see Materials and Methods). The results show that biomolecular components can be used to distinguish between P. miliaceum and S. italica based on the presence or absence of 5 biomarkers—miliacin, α-amyrin methyl ether (ME), and 3 pentacyclic triterpene methyl ethers (PTMEs), although the structures of the 3 PTMEs have yet to be confirmed (Fig. 5, Fig. S3). The total ion current trace shows the whole distribution of aromatic hydrocarbons and ethers extracted from P. miliaceum (Fig. 5A) and S. italica (Fig. 5B) to be in the 54- to 62-min analysis time range. The significant relative abundance of miliacin (compound 2) (Fig. 5D) was found in both modern species (89.0 ± 1.64% for P. miliaceum and 33.8 ± 22.2% for S. italica); the miliacin relative abundances in the aromatic hydrocarbons and ethers were estimated by measuring the area of the miliacin peak on the m/z 189 + 204 + 218 ion-specific chromatogram. However, compounds 1, 4, and 5 are peculiar to P. miliaceum, and compound 3 is peculiar to S. italica.
Fig. 5 shows that the prominent compound products of the archaeological samples (BWG, CS-V-03) are compounds 1, 2, 4, and 5, which are also present in P. miliaceum. However, compounds 1, 4, and 5 are notably absent from S. italica, which contain compounds 2 and 3 only. Thus, the distribution pattern of major biomarkers of BWG and CS-V-03 grain crop remains is similar to that of modern common millet (Fig. 5A). In 2 archeological samples, including BWG and CS-V-03, the relative abundance of miliacin reaches 88.5% and 88.2%, respectively, which is also very similar to that of common millet. These results provide further support to our conclusion and suggest that common millet is an important source of grain crops stored in the storage pits.
Discussion
According to archeobotanical research, the early charred grains of common millet occurred during the initial stages of various Early Neolithic sites (Fig. 1), including Dadiwan (ca. 7.8–7.35 cal kyr BP) (21), Xinglonggou (ca. 8.0–7.5 cal kyr BP) (22), and Yuezhuang (ca. 7.87 cal kyr BP) (23) in North China, but foxtail millet was barely present during these stages. Lee et al. (24) have speculated that the Early Neolithic predominance of broomcorn over foxtail millet at Xinglonggou and Yuezhuang ca. 6000 cal B.C. might be a regional phenomenon, implying that broomcorn millet might have been domesticated earlier than foxtail millet. Our analytical results of both phytoliths and biomolecular components have established that the earliest cereal remains stored in the Cishan Neolithic sites, during ca. 10,300–8,700 cal yr BP, are not foxtail millet, but only common millet. After 8,700 cal yr BP, the grain crops gradually contained 0.4–2.8% foxtail millet. Our study also suggests that common millet was used as a staple food significantly earlier than foxtail millet in northern China. It provides direct evidence to show that, by 10,000 cal yr BP, the early people in northern China had developed various methods of maintenance and multiplication of millet seeds for the next generation, and had known how to store crops of staple food in secure, dry places of storage pits during the Early Neolithic epoch.
Common millet has the lowest water requirement among all grain crops; it is also a relatively short-season crop, and could grow well in poor soils (5, 6, 25). The geographical distribution of both foxtail millet and common millet in China (Fig. S4) shows that foxtail millet is more common in the semiwet eastern areas, and its optimal growth occurs at mean annual temperature (MAT) from 8 to 10 °C and mean annual precipitation (MAP) from 450 to 550 mm. However, common millet is more adapted to the drier interior areas, and its optimal growing conditions occur at MAT from 6 to 8 °C and MAP from 350 to 450 mm (5, 6). The origin and dispersal of millet agriculture is a key problem closely related to the history of human impact on the environment and transformation of natural vegetation.
Paleoenvironmental data from the Weinan section (26–29) (Fig. 1) in the southern part of the Loess Plateau between the Cishan and Dadiwan sites are crucial for understanding the early stage of the forager–cultivator transition. The early Holocene was a period of significant environmental change marked by dry climate conditions as inferred from sediment texture (26, 28), magnetic susceptibility (26, 28), pollen (27), phytoliths (28), and mollusk assemblages (29). These proxy records show an environmental transition from cold–dry (ca. 11,000–8,700 cal yr BP) to warm–wet (ca. 8,700–5,500 cal yr BP) conditions. Many lacustrine and loess records from the Chinese Loess Plateau to Central Asia also support the scenario of a dry climate during the early Holocene (30–34). Under the drier climate conditions, soil development was slowed, and the soil developed on the underlying older and coarser loess of the glacial period was poor in nutrients (28). This raises the possibility that common millet was more significant than foxtail millet in the early stages of food production in North China because it was more adaptable than foxtail millet to the dry condition prevailing during the early Holocene. The common millet cultivation may involve complex selection by natural forces and human activities, although no clear evidence has been documented in this region for the transitions from gathering to cultivation and/or from a wild ancestor to domesticated common millet (1, 2, 5).
Conclusions
Our research indicates that the earliest significant common millet cultivation system was established in the semiarid regions of China by 10,000 cal yr BP, and that the relatively dry condition in the early Holocene may have been favorable for the domestication of common millet over foxtail millet. Our study shows that common millet appeared as a staple crop in northern China ≈10,000 years ago, suggesting that common millet might have been domesticated independently in this area and later spread to Russia, India, the Middle East, and Europe. Nevertheless, like Mesopotamia, where the spread of wheat and barley to the fertile floodplains of the Lower Tigris and Euphrates was a key factor in the emergence of civilization, the spread of common millet to the more productive regions of the Yellow River and its tributaries provided the essential food surplus that later permitted the development of social complexity in the Chinese civilization.
Materials and Methods
Phytolith Extraction from Archaeological Material.
The Cishan site is located on the north terrace of the Ming River (Fig. 1). The geography of this region consists of alluvial terraces covered by loess deposits. Five storage pits (CS-I to CS-V) are situated on the edge of the terrace. The terrace is composed of thick sediments (4–5 m thick for CS-I to CS-IV, 11 m thick for CS-V) that consist of alternating layers of loess and soil. Each of the storage pits (CS-I to CS-V) contains ≈1- to 2-m-thick grain crop remains, being well preserved in the 3- to 5-m-deep loess layer (Table S1).
A total of 46 samples from 5 storage pits (CS-I to CS-V) at the Cishan site and 1 sample (BWG) from the Culture Museum of Cishan were analyzed. Of these, 16 continuous samples were taken at 10-cm intervals from CS-I storage pit. The samples were prepared according to the procedure slightly modified from Piperno (35, 36) and Lu et al. (37). It consists of sodium pyrophosphate (Na4P2O7) deflocculation, treatment with 30% hydrogen peroxide (H2O2) and cold 15% hydrochloric acid (HCl), zinc bromide (ZnBr2; density, 2.35 g/cm3) heavy liquid separation, and mounting on a slide with Canada balsam. Phytolith counting and identification were performed by using a Leica microscope with phase-contrast and microscopic interferometer at 400× magnification.
Millet Molecule Extractions from Modern Plants and Archaeological Material.
We examined modern millet molecules from 6 domesticated samples of Paniceae, including 3 samples of Setaria italica L. Beauv. (foxtail millet) and 3 samples of Panicum miliaceum L. (common millet). Two archaeological samples of grain crop remains, including 1 sample from CS-V storage pit and 1 sample (BWG) in a storage bottle from the Culture Museum of Cishan, were analyzed.
One gram of each of dried sample (smaller than 100 mesh) was ultrasonically extracted for 3 times with a mixture of acetone and pentane at 1:1 (vol/vol). After filtration, the mixture was eluted with n-hexane to give aliphatic hydrocarbons (N). Using dichloromethane in n-hexane passing through the column resulted in satisfactory separation of aromatic hydrocarbons and ethers (A). The extracted subfractions (A) were then analyzed by gas chromatography coupled with mass spectrometry (GC/MS) .
GC-MS was conducted by using an Agilent 6890 N GC-5973 N MSD mass selective detector system (EI 70 eV; ion source temperature, 230 °C) fitted with a fused silica chemically bonded capillary column (J&W DB-5; 0.25 mm in diameter, 30 m long, 0.25 μm film thickness). Each sample was injected onto the column at 280 °C in the splitless mode. After a 1-min isothermal hold at 50 °C, the column temperature was increased by 30 °C/min to 120 °C; and then 3 °C/min to 290 °C, successively, with a 30-min isothermal hold at 290 °C. The flow rate of the helium carrier gas was 1.2 mL/min (40 cm/sec).
Compounds were identified by their retention time within the gas chromatograph, their fragmentation pattern within the mass spectrometer, and by matching their mass spectra to reference spectra from the literature or from libraries (NIST02L). The miliacin relative abundances were estimated by measuring the area of the miliacin peak on the m/z 189 + 204 + 218 ion-specific chromatogram. Note that this value may vary under different GC/MS conditions, and there will be a small change in the abundance of this ion.
Supplementary Material
Acknowledgments.
We are grateful to Dolores R. Piperno for critically reading the original manuscript and her helpful comments in improving this manuscript. We also greatly appreciate the valuable comments from two anonymous reviewers. We are grateful to Jacob Jérémy and A. Gerasimenko for their help on biomolecular analyses, and D. Q. Fuller and L. Qin for their discussions on phytolith analyses. We thank C. D. Shen, Z. Y. Gu, and B. Xu for their help in radiocarbon age measurement. This work was supported by the Chinese Academy of Sciences (100 Talents Program, kzcx2-yw-117), the National Natural Science Foundation of China (40771216; 40325002), Chinese Civilization Origin projects (2006BAK21B02), and the U.S. National Science Foundation (BCS-0623514, ATM-0402475).
Footnotes
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
See Commentary on page 7271.
This article contains supporting information online at www.pnas.org/cgi/content/full/0900158106/DCSupplemental.
References
- 1.Bellwood P. First Farmers: The Origins of Agricultural Societies. Malden, MA: Blackwell; 2005. [Google Scholar]
- 2.Zohary D, Hopf M. Domestication of Plants in the Old Word: The Origin and Spread of Cultivated Plants in West Asia, Europe and the Nile Valley. Oxford: Oxford Univ Press; 2000. [Google Scholar]
- 3.Fuller DQ. Agricultural origins and frontiers in South Asia: A working synthesis. J World Prehist. 2006;20:1–86. [Google Scholar]
- 4.Crawford G. East Asian plant domestication. In: Stark MT, editor. Archaeology of Asia. Malden, MA: Blackwell; 2005. pp. 77–95. [Google Scholar]
- 5.You XL. The question for origin and spread in both foxtail millet and common millet. Agric Hist China. 1993;12:1–13. [Google Scholar]
- 6.Li Y, Wu SZ. Traditional maintenance and multiplication of foxtail millet (Setaria italica. (L.) P. Beauv.) landraces in China. Euphytica. 1996;87:33–38. [Google Scholar]
- 7.Handan Relics Preservation Station, Education Team of Cishan Archaeology. Preliminary excavation of Neolithic sites in Cishan, Hebei Province. Kaogu (Archaeology) 1977;6:361–372. [Google Scholar]
- 8.CPAM Hebei Province, Handan Relics Preservation Station. The Cishan site in Wu'an, Hebei Province. Acta Archaeologica Sinica. 1981;3:303–338. [Google Scholar]
- 9.Tong WH. Original agricultural relics of Cishan sites and related problems. Agric Archaeology. 1984;1:194–207. [Google Scholar]
- 10.Huang QX. Spodograms study and its application in archaeology. Kaogu (Archaeology) 1982;4:418–420. [Google Scholar]
- 11.Zhao ZJ Institute of Archaeology Chinese Academy of Social Sciences and the Institute of Archaeology Swedish National Heritage Board. Archaeology in China and Sweden. Beijing: Science Press; 2006. Domestication of millet—paleoethnobotanic data and ecological perspective; pp. 97–104. [Google Scholar]
- 12.Harvey EL, Fuller DQ. Investigating crop processing using phytolith analysis: The example of rice and millets. J Archaeol Sci. 2005;32:739–752. [Google Scholar]
- 13.Lu HY, et al. Phytoliths analysis for the discrimination of foxtail millet (Setaria italica) and common millet (Panicum miliaceum) PLoS ONE. 2009;4:e4448. doi: 10.1371/journal.pone.0004448. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14.Fukunaga K, Ichitani K, Kawase M. Phylogenetic analysis of the rDNA intergenic spacer subrepeats and its implication for the domestication history of foxtail millet, Setaria italica. Theor Appl Gene. 2006;113:261–269. doi: 10.1007/s00122-006-0291-5. [DOI] [PubMed] [Google Scholar]
- 15.Hunt HV, et al. Millets across Eurasia: Chronology and context of early records of the genera Panicum and Setaria from archaeological sites in the Old World. Veget Hist Archaeobot. 2008;17(Suppl 1):S5–S18. doi: 10.1007/s00334-008-0187-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Nasu H, Momohara A, Yasuda Y, He JJ. The occurrence and identification of Setaria italica (L.) P. Beauv. (foxtail millet) grains from the Chengtoushan site (ca. 5800 cal B.P.) in central China, with reference to the domestication centre in Asia. Veget Hist Archaeobot. 2007;16:481–494. [Google Scholar]
- 17.Pearsall DM, et al. Distinguishing rice (Oryza sativa Poaceae) from wild Oryza species through phytolith analysis: Results of preliminary research. Econ Bot. 1995;49:183–196. [Google Scholar]
- 18.Piperno DR, Weiss E, Holst I, Nadel D. Processing of wild cereal grains in the Upper Palaeolithic revealed by starch grain analysis. Nature. 2004;430:670–673. doi: 10.1038/nature02734. [DOI] [PubMed] [Google Scholar]
- 19.Lu HY, et al. Millet noodles in Late Neolithic China. Nature. 2005;437:967–968. doi: 10.1038/437967a. [DOI] [PubMed] [Google Scholar]
- 20.Jacob J, et al. Millet cultivation history in the French Alps as evidenced by a sedimentary molecule. J Archaeol Sci. 2008;35:814–820. [Google Scholar]
- 21.Liu CJ, Kong ZC. The morphological comparison of grains between foxtail and broomcorn millet and its application for identification in archaeology remains. Kaogu (Archaeology) 2004;3:76–83. [Google Scholar]
- 22.Zhao ZJ. Palaeoethnobotany and its new achievements in China. Kaogu (Archaeology) 2005;7:42–49. [Google Scholar]
- 23.Crawford GW, Chen XX, Wang JH. Houli Culture rice from the Yuezhuang site, Jinan. Orient Archaeol. 2006;3:247–251. [Google Scholar]
- 24.Lee GA, Crawford GW, Li L, Chen XC. Plants and people from the Early Neolithic to Shang periods in North China. Proc Natl Acad Sci USA. 2007;104:1087–1092. doi: 10.1073/pnas.0609763104. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 25.Nesbitt M. Grains. In: Prance G, Nesbitt M, editors. Cultural History of Plants. London: Routledge; 2005. pp. 45–60. [Google Scholar]
- 26.Guo ZT, et al. High frequency pulses of East Asian monsoon climate in the last two glaciations: Link with the North Atlantic. Clim Dyn. 1996;12:701–709. [Google Scholar]
- 27.Sun XJ, Song CQ, Wang FY, Sun MR. Vegetation history of the loess plateau of China during the last 100,000 years based on pollen data. Quat Int. 1997;37:25–36. [Google Scholar]
- 28.Lu HY, Wu NQ, Liu KB, Jiang H, Liu TS. Phytoliths as quantitative indicators for the reconstruction of past environmental conditions in China II: Palaeo-environmental reconstruction in the Loess Plateau. Quat Sci Rev. 2007;26:759–772. [Google Scholar]
- 29.Lu HY, Zhang JP. Neolithic cultural evolution and Holocene climate in the Guanzhong basin, Shaanxi, China. Quat Sci. 2008;28:1050–1060. [Google Scholar]
- 30.Feng ZD, Tang LY, Wang HB, Ma YZ, Liu K-b. Holocene vegetation variations and the associated environmental changes in the western part of the Chinese Loess Plateau. Palaeogeogr Palaeoclimatol Palaeoecol. 2006;241:440–456. [Google Scholar]
- 31.Feng ZD, An CB, Tang LY, Jull AJT. Stratigraphic evidence of a megahumid climate between 10,000 and 4000 years B.P. in the western part of the Chinese Loess Plateau. Glob Planet Change. 2004;43:145–155. [Google Scholar]
- 32.Boomer I, Aladin N, Plotnikov I, Whatley R. The palaeolimnology of the Aral Sea: A review. Quat Sci Rev. 2000;19:1259–1278. [Google Scholar]
- 33.Fowell SJ, Hansen BCS, Peck JA, Khosbayar P, Ganbold E. Mid to late Holocene climate evolution of the Lake Telmen Basin, North Central Mongolia, based on palynological data. Quat Res. 2003;59:353–363. [Google Scholar]
- 34.Feng ZD, An CB, Wang HB. Holocene climatic and environmental changes in the arid and semi-arid areas of China: A review. Holocene. 2006;16:119–130. [Google Scholar]
- 35.Piperno DR. Phytolith Analysis. An Archaeological and Geological Perspective. San Diego: Academic; 1988. [DOI] [PubMed] [Google Scholar]
- 36.Piperno DR. Phytoliths: A Comprehensive Guide for Archaeologists and Paleoecologists. Oxford: Alta Mira Press; 2006. [Google Scholar]
- 37.Lu HY, et al. Phytoliths as quantitative indicators for the reconstruction of past environmental conditions in China I: Phytolith-based transfer functions. Quat Sci Rev. 2006;25:945–959. [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.