The initial colonization of the Americas remains a highly debated topic1, and the exact timing of the first arrivals is unknown. The earliest archaeological record of Mexico—which holds a key geographical position in the Americas—is poorly known and understudied. Historically, the region has remained on the periphery of research focused on the first American populations2. However, recent investigations provide reliable evidence of a human presence in the northwest region of Mexico3,4, the Chiapas Highlands5, Central Mexico6 and the Caribbean coast7,8,9 during the Late Pleistocene and Early Holocene epochs. Here we present results of recent excavations at Chiquihuite Cave—a high-altitude site in central-northern Mexico—that corroborate previous findings in the Americas10,11,12,13,14,15,16,17of cultural evidence that dates to the Last Glacial Maximum (26,500–19,000 years ago)18, and which push back dates for human dispersal to the region possibly as early as 33,000–31,000 years ago. The site yielded about 1,900 stone artefacts within a 3-m-deep stratified sequence, revealing a previously unknown lithic industry that underwent only minor changes over millennia. More than 50 radiocarbon and luminescence dates provide chronological control, and genetic, palaeoenvironmental and chemical data document the changing environments in which the occupants lived. Our results provide new evidence for the antiquity of humans in the Americas, illustrate the cultural diversity of the earliest dispersal groups (which predate those of the Clovis culture) and open new directions of research.
The data that support the findings of this study are available in the Article and its Supplementary Information. Raw data and sequence alignments are available from the European Nucleotide Archive under accession number PRJEB37914. The exact coordinates of Chiquihuite Cave are available from C.F.A. on reasonable request. C.F.A. can also be contacted at firstname.lastname@example.org.
Code for R and OxCal is noted within the Supplementary Information. Code for environmental DNA data analysis can be found at https://github.com/miwipe/ngsLCA.
Meltzer, D. J. The Great Paleolithic War: How Science Forged an Understanding of America’s Ice Age Past (Univ. Chicago Press, 2015).
Ardelean, C. F. et al. in People and Culture in Ice Age Americas: New Dimensions in Paleoamerican Archaeology (eds Suarez, R. & Ardelean, C. F.) 108–133 (Univ. Utah Press, 2019).
Sanchez, G. et al. Human (Clovis)–gomphothere (Cuvieronius sp.) association ~13,390 calibrated yBP in Sonora, Mexico. Proc. Natl Acad. Sci. USA 111, 10972–10977 (2014).
Des Lauriers, M. R., Davis, L. G., Turnbull, J., Southon, J. R. & Taylor, R. E. The earliest fish hooks from the Americas reveal fishing technology of Pleistocene maritime foragers. Am. Antiq. 82, 498–516 (2017).
Acosta, G. et al. Climate change and peopling of the Neotropics during the Pleistocene–Holocene transition. Bol. Soc. Geol. Mex. 70, 1–19 (2018).
González, S., Jiménez López, C., Hedges, R., Pompa y Padilla, J. A. & Huddart, D. Early humans in Mexico: new chronological data. In El Hombre Temprano en América y sus Implicaciones en el Poblamiento de la Cuenca de México: Primer Simposio Internacional (eds Jiménez López, C. et al.) 67–77 (Instituto Nacional de Anthrología e Historia, 2006).
González, A. et al. in Paleoamerican Odyssey (eds Graf, K. E. et al.) 323–337 (Center for the Study of the First Americans, 2014).
Chatters, J. C. et al. Late Pleistocene human skeleton and mtDNA link Paleoamericans and modern Native Americans. Science 344, 750–754 (2014).
Stinnesbeck, W. et al. The earliest settlers of Mesoamerica date back to the late Pleistocene. PLoS ONE 12, e0183345 (2017).
Williams, T. J. et al. Evidence of an early projectile point technology in North America at the Gault Site, Texas, USA. Sci. Adv. 4, eaar5954 (2018).
Waters, M. R. et al. Pre-Clovis projectile points at the Debra L. Friedkin site, Texas–implications for the Late Pleistocene peopling of the Americas. Sci. Adv. 4, eaat4505 (2018).
Jenkins, D. L. et al. Clovis age Western Stemmed projectile points and human coprolites at the Paisley Caves. Science 337, 223–228 (2012).
Adovasio, J. M., Gunn, J. D., Donahue, J. & Stuckenrath, R. Meadowcroft Rockshelter, 1977: an overview. Am. Antiq. 43, 632–651 (1978).
Halligan, J. J. et al. Pre-Clovis occupation 14,550 years ago at the Page–Ladson site, Florida, and the peopling of the Americas. Sci. Adv. 2, e1600375 (2016).
Dillehay, T. D. Monte Verde, a Late Pleistocene Settlement in Chile: The Archaeological Context and Interpretation (Smithsonian Institution Press, 1997).
Bourgeon, L., Burke, A. & Higham, T. Earliest human presence in North America dated to the Last Glacial Maximum: new radiocarbon dates from Bluefish Caves, Canada. PLoS ONE 12, e0169486 (2017).
Davis, L. G. et al. Late Upper Paleolithic occupation at Cooper’s Ferry, Idaho, USA, ~16,000 years ago. Science 365, 891–897 (2019).
Clark, P. U. et al. The Last Glacial Maximum. Science 325, 710–714 (2009).
Ardelean, C. F. Archaeology of Early Human Occupations and the Pleistocene–Holocene Transition in the Zacatecas Desert, Northern Mexico. PhD thesis, Univ. Exeter (2013).
Brock, F., Higham, T., Ditchfield, P. & Bronk Ramsey, C. Current pretreatment methods for AMS radiocarbon dating at the Oxford Radiocarbon Accelerator Unit (ORAU). Radiocarbon 52, 103–112 (2010).
Wang, Y., Amundson, R. & Trumbore, S. Radiocarbon dating of soil organic matter. Quat. Res. 45, 282–288 (1996).
Pessenda, L. C. R., Gouveia, S. E. M. & Aravena, R. Radiocarbon dating of total soil organic matter and humin fraction and its comparison with 14C ages of fossil charcoal. Radiocarbon 43, 595–601 (2001).
Rasmussen, S. O. et al. A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. Quat. Sci. Rev. 106, 14–28 (2014).
Dunham, R. J. Classification of carbonate rocks according to depositional texture. In Classification of Carbonate Rocks—A Symposium (ed. Ham, W. E.) 108–121 (The American Association of Petroleum Geologists, 1962).
Folk, R. L. Practical petrographic classification of limestones. Bull. Am. Assoc. Petrol. Geol. 43, 1–38 (1959).
Smallwood, A. M. & Jennings, T. A. Clovis: On the Edge of a New Understanding (Texas A&M Univ. Press, 2014).
Goebel, T. & Keene, J. L. in Archaeology in the Great Basin and Southwest: Papers in Honor of Don D. Fowler (eds Parezo, N. J. & Janetski, J. C.) 35–60 (Univ. Utah Press, 2014).
Graf, K. E. & Buvit, I. Human dispersal from Siberia to Beringia: assessing a Beringian standstill in light of the archaeological evidence. Curr. Anthropol. 58, S583–S603 (2017).
Willerslev, E. et al. Diverse plant and animal genetic records from Holocene and Pleistocene sediments. Science 300, 791–795 (2003).
Pedersen, M. W. et al. Postglacial viability and colonization in North America’s ice-free corridor. Nature 537, 45–49 (2016).
Villaseñor, J. L. Checklist of the native vascular plants of Mexico. Rev. Mex. Biodivers. 87, 559–902 (2016).
Henderson, A., Bernal, R. & Galeano-Garces, G. Field Guide to the Palms of the Americas (Princeton Univ. Press, 1997).
Slon, V. et al. Neandertal and Denisovan DNA from Pleistocene sediments. Science 356, 605–608 (2017).
Middleton, W. D. et al. The study of archaeological floors: methodological proposal for the analysis of anthropogenic residues by spot tests, ICP-OES, and GC-MS. J. Archaeol. Method Theory 17, 183–208 (2010).
Rasmussen, S. O. et al. A new Greenland ice core chronology for the last glacial termination. J. Geophys. Res. 111, D06102 (2006).
Sedlock, R. L., Ortega-Gutiérrez, F. & Speed, R. C. Tectonostratigraphic Terranes and Tectonic Evolution of Mexico (GSA Special Papers Volume 278) (Geological Society of America, 1993).
Padilla y Sánchez, R. J. Geological Map of the Curvature of Monterrey, Mexico (GSA, 2006).
Ramsey, C. B., Higham, T. & Leach, P. Towards high-precision AMS: progress and limitations. Radiocarbon 46, 17–24 (2004).
Hajdas, I. Radiocarbon dating and its applications in Quaternary studies. E&G Quat. Sci.J. 57, 2–24 (2008).
Abbott, M. B. & Stafford, T. W. Radiocarbon geochemistry of modern and ancient Arctic lake systems, Baffin Island, Canada. Quat. Res. 45, 300–311 (1996).
Longin, R. New method of collagen extraction for radiocarbon dating. Nature 230, 241–242 (1971).
Solís, C. et al. A new AMS facility in Mexico. Nucl. Instrum. Methods Phys. Res. B 331, 233–237 (2014).
International Chemical Analysis. International Chemical Analysis. https://www.radiocdating.com/ (accessed 18 July 2018) (2017).
Bronk Ramsey, C. Bayesian analysis of radiocarbon dates. Radiocarbon 51, 337–360 (2009).
Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55, 1869–1887 (2013).
Bronk Ramsey, C. Dealing with outliers and offsets in radiocarbon dating. Radiocarbon 51, 1023–1045 (2009).
Aitken, M. J. Introduction to Optical Dating: The Dating of Quaternary Sediments by the Use of Photon-stimulated Luminescence (Clarendon, 1998).
Mirazón Lahr, M. et al. Inter-group violence among early Holocene hunter-gatherers of West Turkana, Kenya. Nature 529, 394–398 (2016).
Bøtter-Jensen, L., Bulur, E., Duller, G. A. T. & Murray, A. S. Advances in luminescence instrument systems. Radiat. Meas. 32, 523–528 (2000).
Richter, D., Richter, A. & Dornich, K. Lexsyg smart—a luminescence detection system for dosimetry, material research and dating application. Geochronometria 42, 202–209 (2015).
Murray, A. S. & Wintle, A. G. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiat. Meas. 32, 57–73 (2000).
Wintle, A. G. & Murray, A. S. A review of quartz optically stimulated luminescence characteristics and their relevance in single-aliquot regeneration dating protocols. Radiat. Meas. 41, 369–391 (2006).
Banerjee, D., Murray, A. S., Bøtter-Jensen, L. & Lang, A. Equivalent dose estimation using a single aliquot of polymineral fine grains. Radiat. Meas. 33, 73–94 (2001).
Wallinga, J., Murray, A. S. & Bøtter-Jensen, L. Measurement of the dose in quartz in the presence of feldspar contamination. Radiat. Prot. Dosimetry 101, 367–370 (2002).
Durcan, J. A., King, G. E. & Duller, G. A. T. DRAC: dose rate and age calculator for trapped charge dating. Quat. Geochronol. 28, 54–61 (2015).
Barba, L. Chemical residues in lime-plastered archaeological floors. Geoarchaeology 22, 439–452 (2007).
Piperno, D. R. Phytoliths: A Comprehensive Guide for Archaeologists and Paleoecologists (Rowman Altamira, 2006).
Piperno, D. R. & Pearsall, D. M. The Silica Bodies of Tropical American Grasses: Morphology, Taxonomy, and Implications for Grass Systematics and Fossil Phytolith Identification (Smithsonian Institution, 1998).
Gallego, L. & Distel, R. A. Phytolith assemblages in grasses native to central Argentina. Ann. Bot. 94, 865–874 (2004).
Fredlund, G. G. & Tieszen, L. T. Modern phytolith assemblages from the North American Great Plains. J. Biogeogr. 21, 321–335 (1994).
Colinvaux, P., De Olieira, P. E. & Moreno Patino, J. E. Amazon Pollen Manual and Atlas (Harwood Academic, 1999).
Roubik, D. W. & Moreno Patiño, J. E. Pollen and Spores of Barro Colorado Island (Missouri Botanical Garden, 1991).
Markgraf, V. & d’Antoni, H. L. Pollen Flora of Argentina (Univ. Arizona Press, 1978).
Johnston, I. M. Plants of Coahuila, eastern Chihuahua, and adjoining Zacatecas and Durango, V. J. Arnold Arbor. 25, 133–182 (1944).
González-Tagle, M. A., Schwendenmann, L., Pérez, J. J. & Schulz, R. Forest structure and woody plant species composition along a fire chronosequence in mixed pine–oak forest in the Sierra Madre Oriental, Northeast Mexico. For. Ecol. Manage. 256, 161–167 (2008).
Ludlow Wiechers, B., Almeida Leñero, L. & Sugiura, Y. Palinomorfos del Holoceno en la cuenca alta del Río Lerma, Estado de México, México. Bol. Sociedad Botánica de México 72, 59–105 (2003).
Meyer, M. & Kircher, M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010, db.prot5448 (2010).
Jónsson, H., Ginolhac, A., Schubert, M., Johnson, P. L. F. & Orlando, L. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics 29, 1682–1684 (2013).
The field explorations and part of the laboratory studies were made possible by the special sponsorship from the Government of the State of Zacatecas, through the Consejo Zacatecano de Ciencia, Tecnología e Innovación (COZCyT); we thank its consecutive directors, G. A. Mercado-Sánchez and A. Enciso-Muñoz, and the Governor of Zacatecas, A. Tello-Cristerna. Seed money for cave exploration came from the Center for American Paleolithic Research (CAPR); we thank S. Holen, K. Holen and the members of the board. Fieldwork, laboratory analyses and publication expenses were partially covered by CONACYT grant CB-2016-286130. The Concepción del Oro municipality, and A. Maldonado-Falcón, offered administrative and occasional financial assistance. Radiocarbon work at Oxford was supported by the NERC Radiocarbon Facility (NRCF), Merton College, Santander and the Clarendon Fund. We thank all colleagues at the ORAU. D. Peat contributed greatly to laboratory preparation of OSL samples. A. Ocaña and I. Alarcón participated in the identification of animal bones. The environmental DNA work was supported by the Lundbeck Foundation, the Novo Nordic Foundation, the Wellcome Trust Foundation, the Carlsberg Foundation and the Danish National Research Foundation. We thank INAH’s Archaeology Council for authorization and legal permits; the inhabitants of Guadalupe Garzarón for accepting this project in their territory and participating in the caravans; and J. Martínez-Ledezma for his constant support as an on-site administrator and liaison with the community.
The authors declare no competing interests.
Peer review information Nature thanks Deborah M. Pearsall, Fiona Petchey and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
a, Digital elevation model (DEM) map of Mexico with the study area (rectangle) in relationship to relevant modern and prehistoric localities. b, The study region (Concepción del Oro endorheic basin), with Chiquihuite Cave (green dot) on the Astillero Mountains. DEM mosaic generated from ortophotographs and elevation data from National Institute of Statistics and Geography (INEGI) of Mexico. c, Chiquihuite Peak, seen from the south; the arrow indicates the entrance of the cave. d, The main chamber, looking west towards the double-eyed entrance. e, Contact between the limestone lintel of the ancient entrance (25 m west of X-12), and the debris that sealed it, probably at the end of the Pleistocene epoch. f, The south profile of the central–eastern squares, as in Fig. 1b. g, Western profiles of the central sector, showing the exposed interface 1210 tilted in the centre of the photograph. h, East–west cross-section of the northern sector of the main chamber of the cave, along the southern profile of the dig. i, X-12 cross-section, showing the inclination of the cave floor and the stepped approach. The DEMs in a, b were created with ArcMap/ArcScene (by ESRI) using open-access topographical data provided by INEGI.
Extended Data Fig. 2 Three-dimensional photogrammetric model of excavation X-12 and the location of samples extracted and analysed for ancient environmental DNA.
a, North-facing profile. b, East-facing profile. Triangles indicate the respective locations of environmental DNA samples extracted in 2019, corresponding to different stratigraphic units (‘UE’, as shown in the key). The circled samples are the initial ones, from the 2016–2017 excavation. c, Complete 3D photogrammetric view of the excavation towards the southern profile. The reddish-brown layers visible along the entire southern profile represent the stratigraphic component C (SC-C), starting downwards with stratum 1212, the terminal LGM chronostratigraphic marker that defines the separation between SC-B and SC-C (Fig. 1b, Extended Data Fig. 7).
a, Horizontal (x and y axes) spatial distribution of lithic artefacts, radiocarbon and OSL samples on the excavation grid. b, Side view (from the south) of the vertical distribution (x and z axes) of finds and samples. Illustrated dating samples: (1) LEMA-635.1.1 and LEMA-635.1.2; (2) OxA-36608 and LEMA-978.1.1; (3) OxA-36610; (4) LEMA-574.1.1; (5) LEMA-640.1.1; (6) OxA-36496; (7) OxA-36613; (8) X-7227; (9) OxA-36614; (10) LEMA-636.1.1 and LEMA-636.1.2; (11) LEMA-576.1.1; (12) LEMA-577.1.1; (13) X-7228; (14) X-7231 and X-7232; (15) X-7229; (16) OxA-36530; (17) OxA-36633; (18) OxA-36609; (19) OxA-36612; (20) OxA-36360; (21) BETA-345055; (22) X-4135; (23) OxA-34965; (24) ICA-16OS/0510; (25) PRI-5414; (26) OxA-36616; (27) OxA-36620; (28) OxA-36618; (29) OxA-36615; (30) OxA-36617; (31) OxA-36621; (32) OxA-36619; (33) OxA-36753; (34) OxA-36611; (35) X-7233; (36) LEMA-573.1.1; (37) OxA-36359; (38) LEMA-892.1.2; (39) OxA-36622; (40) OxA-36623; (41) OxA-36634; (42) LEMA-893.1.1; (43) LEMA-977.1.1; (44) OxA-36624; and (45) OxA-36625.
Extended Data Fig. 4 Taxonomic profiles of animals (Amniota) and plants (Viridiplantae) identified by ancient environmental DNA.
a, b, Animals presented as the proportion of reads found of each taxa and plotted as a bar plot (b) and stratigraphic plot (a). c, d, All plants are presented as the proportion of reads found of each taxa and plotted as a bar plot (c) and stratigraphic plot (d). *Taxa also found by pollen, phytoliths or faunal remains.
a–c, Cores. d, e, Bifacial preforms on ovoid nodules. f–n, Flakes. o–t, Blades. u–x, Microliths. y, w, Burins. z, a′, Scrapers. b′–l′, Points and point-like shapes. m′–p′, Geometric items made by fracturing calcite laminae. Artefacts i, n, o, b′, c′, f′, and k′ are from SC-C; all others are from SC-B. Scale bars, 1 cm.
a, Core. b, Flake with isolated platform. c, Flake with lipped platform. d–i, Blades and microlith blade segments. j, Circular scraper on trimmed flake. k, Possible preform. l, Point on plaquette. m, Bifacial point preform. n, o, s–a′, d′–g′, Transversal points (obtained by slightly modifying transversal flake blanks). p–r, b′, c′, Other points. h′, A geometric, point-like shape on calcite sheet. Most artefacts were discovered in SC-B. Specimens x–z, b′, g′, h′ are from SC-C. Scale bars, 1 cm.
a, b, Excavation grid diagram (a) and maximum-depth diagram (b) showing the position of the profiles depicted below. For correlation reasons, all profiles show the upper contact of stratum 1210 with 1212 marked with a white contour. The position of OSL samples is indicated (labels beginning with X-). c, South profile, squares K2, K3, L2 and L3. d, South profile, square H2. e, South profile, square G2. f, South profile, square F2. g, South profile, squares M3, M4, N3 and N4. h, West profile, square J4. i, West profile, squares I5 and I6. j, Eastern profile, squares M4, M5 and M6 (profile removed during the 2016–2017 winter excavations). k, Eastern profile, squares N4, N5 and N6 (new eastern profile after the excavation of the one shown in j).
This file contains the concentration of data from all the different scientific approaches included in the paper. Each proxy-related section is numbered consecutively, from 1 to 10, such as: 1.1. Geology; 1.2. Radiocarbon dating and Bayesian modelling. 1.3. Lithic artefact metrics. 1.4. Chemical residues. 1.5. Faunal remains. 1.6. Phytoliths and pollen. 1.7. Thin-section and micromorphology. 1.8. Commercial radiocarbon dating methods. 1.9. Optically stimulated luminescence (OSL) dating. 1.10. Environmental DNA. Each section number is referred as such in the manuscript.
This file contains plots used in the assessment of ancient DNA authentication for the animal taxa found. The first column contains the read length distribution (in Bp) of the reads aligning to the reference (see supplementary metadata file), the second columns plots the number of mismatches (edit distance) of each read to the same reference. In the third column, the fraction of C-T transitions due to DNA damage, as calculated by mapDamage2.0, is plotted. Each new row contains a new species from a specific layer.
This file contains plots used in the assessment of ancient DNA authentication for the plant taxa found. The first column contains the read length distribution (in Bp) of the reads aligning to the reference (see supplementary metadata file), the second columns plots the number of mismatches (edit distance) of each read to the same reference. In the third column, the fraction of C-T transitions due to DNA damage, as calculated by mapDamage2.0, is plotted. Each new row contains a new species from a specific layer.
A link to the 3D digital photogrammetric model of excavation X-12. This is an on-line dynamic model that shows the tridimensional photographic modelling of the excavation, same as in Fig. 1c and Extended Data Fig. 2. It includes the display of the entire spectrum of eDNA samples extracted in 2019 for aDNA reassessment, as well as the northern wall of the cave. The viewer can rotate the model around the x,y,z axis.
The eDNA metadata file contains information relating to layer in excavation, excavation coordinate, total sequences sequenced, reads after quality control, sequencing batch, extraction batch, number of replicated samples and total counts of total taxonomically classified reads and total classified reads at the different kingdoms; Bacteria, Viridiplantae, and Amniota. Including reads checked for DNA damage at different mapping qualities (MQ) and the proportion of damage found.
: A dynamic tridimensional view of the artefact-and-sample plot depicted in Extended Data Fig. 3. This video allows the viewer to appreciate the 3D relationship between the items without the need to run the model in an ArcScene (by ESRI) software. Motion starts at 00: 08.
About this article
Cite this article
Ardelean, C.F., Becerra-Valdivia, L., Pedersen, M.W. et al. Evidence of human occupation in Mexico around the Last Glacial Maximum.
Nature 584, 87–92 (2020). https://doi.org/10.1038/s41586-020-2509-0