Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2017; 12(8): e0183345.
Published online 2017 August 30. doi:  10.1371/journal.pone.0183345
PMCID: PMC5576649

The earliest settlers of Mesoamerica date back to the late Pleistocene

Wolfgang Stinnesbeck, Funding acquisition, Investigation, Supervision, Writing – original draft, Writing – review & editing,1,* Julia Becker, Formal analysis, Methodology, Writing – original draft, Writing – review & editing,2 Fabio Hering, Data curation, Formal analysis, Methodology, Visualization, Writing – original draft,1 Eberhard Frey, Conceptualization, Funding acquisition, Writing – original draft,3 Arturo González González, Conceptualization, Funding acquisition, Investigation,4 Jens Fohlmeister, Formal analysis, Funding acquisition, Investigation, Methodology,5 Sarah Stinnesbeck, Data curation, Investigation, Visualization, Writing – original draft,3 Norbert Frank, Data curation, Formal analysis, Methodology, Supervision,6 Alejandro Terrazas Mata, Data curation, Investigation, Methodology,7 Martha Elena Benavente, Data curation, Investigation, Methodology,7 Jerónimo Avilés Olguín, Data curation, Investigation,8 Eugenio Aceves Núñez, Data curation, Investigation, Methodology,8 Patrick Zell, Visualization, Writing – review & editing,9 and Michael Deininger, Formal analysis, Methodology, Writing – original draft, Writing – review & editing10
Michael D. Petraglia, Editor


Preceramic human skeletal remains preserved in submerged caves near Tulum in the Mexican state of Quintana Roo, Mexico, reveal conflicting results regarding 14C dating. Here we use U-series techniques for dating a stalagmite overgrowing the pelvis of a human skeleton discovered in the submerged Chan Hol cave. The oldest closed system U/Th age comes from around 21 mm above the pelvis defining the terminus ante quem for the pelvis to 11311±370 y BP. However, the skeleton might be considerable older, probably as old as 13 ky BP as indicated by the speleothem stable isotope data. The Chan Hol individual confirms a late Pleistocene settling of Mesoamerica and represents one of the oldest human osteological remains in America.


The early settlement of the Americas is a controversial subject. While genetic evidence suggests a Beringian origin of the earliest inhabitants of the continent [15], migration routes used for the southward spread of these humans and the timing of human arrival across the Americas are presently reevaluated [611]. The hypothesis of a routing across the exposed Bering land bridge through an ice-free corridor between retreating North American glaciers, at about 12.6 thousand years (ky) ago [12], is increasingly challenged by the discovery of evidence predating the earliest North American widespread archaeological complex, the Clovis culture [1315]. Based on new sites in both North and South America this emerging consensus suggests that people must have arrived in North America as early as 22 ky ago (e.g. [2, 6, 10, 1618].

Osteological evidence for early American settlers is scarce and majorly fragmentary, with at present only a few individuals, from North-, Central- and South America, securely predating 11 ky BP [1924]. Recently, Chatters et al. [25] documented a well-preserved prehistoric skull of a young girl from the submerged Hoyo Negro (Black Hole) sinkhole of the Tulum area, southern Mexico. The individual was 14C-dated to 10976±20 y BP (12910–11750 cal y BP) based on bioapatite from tooth enamel [25]. Previously, González González et al. [20, 21] published a similar 14C age for a human skeleton from the nearby Naharon cave, also located close to Tulum (11670±60 14C y BP; 13721–13354 y cal BP). These two skeletons belong to the oldest 14C-dated New World Homo sapiens. They emphasize the enormous potential of the Tulum system of submerged caves as an archive for the human settlement history in America.

Nevertheless, 14C data from the Tulum submerged caves must be considered with extreme caution because the amount of collagen in both bones and teeth of these individuals is extremely low. This lack of collagen is the result of the exposure of the osteological remains for thousands of years to alternating salt- and fresh water environments [21, 2527]. In addition, bioapatite is highly susceptible to contamination with fossil carbon resulting in false, mostly older ages [25], apart from a general problem with radiocarbon ages in the tropics [26]. Here we use U-series techniques to date a human skeleton that was discovered in the Chan Hol cave near Tulum (here referred to as Chan Hol II, because there are other prehistoric human skeletons in this cave [21]. The pelvis of this skeleton, previously documented briefly by [21], is overgrown by a stalagmite (CH-7; Fig 1).

Fig 1
Geographical position and site of the Chan Hol skeleton.

Regional setting

The Tulum submerged cave system developed in almost horizontally layered, thick-bedded, shallow-water carbonate bedrock of Neogene ages [28, 29] and is the result of intensive karstification during the Pleistocene [30], caused by a series of sea-level oscillations and changes in the overall hydrology of the area. During Pleistocene stadials, sea-level was more than 100 m below the recent level [31, 32], thus leaving large parts of the cave system dry and accessible. In contrast, during interstadials of the Pleistocene and the early Holocene, between 13 and 7.6 ky BP [33], the karst labyrinth was flooded preserving both archaeological and climate archives. Recent water levels were reached at approximately 4500 y BP [34], although oscillations of up to a few meters are known to have occurred during Maya times [3537]. The Tulum cave system contains a coastal density stratified aquifer with a freshwater layer overlying penetrating seawater. The depth of the halocline depends on the global sea-level as well as on the thickness of the superimposed freshwater layer and is controlled by the distance to the coastline as well as the amount of precipitation, with a hydraulic gradient across the Yucatan Peninsula (YP) of between 0.5 and 100 mm km-1 [3840]. Water level in the Tulum area is thus approximately equivalent to mean sea-level. Sea-level rise on the YP was predominantly controlled by eustasy, because the peninsula is tectonically stable and glacial isostatic adjustments in this tropical area were minor (e.g. [32, 40, 41]).

The Chan Hol II site and the skeleton

The Chan Hol II skeleton was located at 20°9.467' N, 87°34.165' W, 15 km southwest of Tulum, Quintana Roo, southern Mexico, and about 11.5 km from the coast line (Fig 1). It was discovered in a low cave tunnel approximately 1240 m southwest of the Chan Hol sinkhole, at about 8.5 m water depth. The maximum depth of the Chan Hol cave is about 13 m below present day sea-level. The halocline is at a depth of about 9 m. Contrary to most other preceramic human skeletons discovered so far in the submerged caves of Tulum, which have been located at water depths of 20 to >30 m [20, 21, 25], the shallow Chan Hol cave must have been accessible up to the middle Holocene [33, 34]. This interpretation is supported by an U/Th age of 5700 y BP of a stalagmite tip from the 8.5 m depth level at Chan Hol. For the nearby Outland cave, [42] postulate a flooding initiating from 8100 cal y BP to a complete inundation at around 6000 cal y BP, which agrees with the results provided here.

Fossil remains were also discovered by us in the extended Chan Hol cave system, though not close to the Chan Hol II skeleton. They include isolated bones of a megalonychid ground sloth, and of extant pacas (Dasyprocta punctata), spider monkeys (Ateles geoffroyi), peccaries (cf. Tayassu tajacu) and white-tailed deer (Odocoileus virginianus).

The Chan Hol II skeleton was brought to our knowledge (JAO) in February 2012 through photos in social media. Soon after, the site was vandalized between the 16th and 23rd of March 2012 and all easily collectable bones were stolen. Photographs of the Chan Hol II skeleton prior to this vandalism provide strong evidence that the skeleton must have been more than 80% complete with the skull excellently preserved (Figs (Figs1B,1B, ,1C1C and and22).

Fig 2
The original Chan Hol II skeleton.

The photographs also allow us to reconstruct the original position of the skeleton and indicate that it was preserved nearly articulated (Fig 1B), with the corpse lying on its back. This is concluded from the ribs covering the vertebral column (Figs (Figs1B1B and and2A)2A) and the position of the left angled femur showing its caudal face (Fig 2B1 and 2C1). The head was inclined slightly to the right. The right leg was fully extended, while the left leg was flexed at the knee at an angle of 20° (Fig 1C). The right femur was still in an articulated position with the pelvis. Based on these data we speculate that the Chan Hol II human died in the cave and that it was not intentionally buried, but there is no positive evidence for this interpretation. Also, no artifacts were identified close to the skeleton on the photographs of the original site or during our collection.

After the looting of Chan Hol II only about 10% of the skeleton remained on site (Fig 2). Among the 155 bone fragments collected there are two auditive labyrinths, an incus, fragments of the temporal, the occipital condyle, four teeth (two incisiors, two molars), a mandibular fragment, the hyoid, numerous ribs, carpals and metacarpals, the right pelvis, a patella, tarsals and metatarsals. Embedding of the pelvis in a stalagmite (CH-7) likely prevented this bone from being stolen (Fig 3). Interpretation of sex and age of the Chan Hol II human is speculative, given that our collection only consists of highly fragmentary bones and a few photographs from the original site. We suggest that the individual was a young adult based on the osteophytes in the vertebral bodies, eruption of the third molar in the right half of the mandible, and an epiphysis that was completely fused. Based on the sciatic notch, it was likely a male.

Fig 3
The Chan Hol II pelvis.


The CH-7 stalagmite

Stalagmite CH-7 is 107 mm high with a mean diameter of ~70 mm (Figs (Figs33 and and4).4). It encloses the human pelvis of the Chan Hol II skeleton at ~95 mm below the top of the stalagmite. The bone is under- and overlain by brown-colored calcareous stalagmite layers, which are each between 1 and 3 mm thick. The internal section of the stalagmite exhibits a succession of milky white with less frequent dark brown calcitic laminae along its long axis. Layers underlying the 3–5 mm thick solid layers below the bone show a wide range of porosity, resembling lime tufa, and they are thus distinct from the dense overlying layers. In addition, the underlying layers are irregularly flexed and bent, which is not seen directly below and above the pelvic fragment, where the lime layers are substantially more homogenous.

Fig 4
Cross-section of CH-7 and results of U/Th age measurements.

U/Th analysis

17 samples were taken along the growth axis of stalagmite CH-7 and were dated using mass spectrometric measurements of the natural disequilibrium isotope ratios of U and Th [43]. U concentrations of these samples are variable with values ranging between 0.2605 and 12.78 ppm. Initial δ234U values are close to the secular equilibrium with values ranging between -33.7‰ and +23.4‰ (Table 1).

Table 1
U/Th measurements of the CH-7 stalagmite from Chan Hol cave.

Contamination of the samples with non-carbonate particles is mostly insignificant. This is indicated by 232Th concentrations of <10 ng g-1, introducing only minor age corrections when using the bulk Earth 230Th/232Th ratio as a correlate [43]. Thus, corrected U/Th ages for CH-7 are variable, ranging between 7.82 ky BP and 12.09 ky BP with most data, however, clustering between 9.8 and 12.1 ky (Fig 4A, Table 1).

Only the calculated U-Th ages of samples taken from the stalagmite growth axis above 72 mm can be interpreted as in stratigraphic order and reflecting a closed system behavior. The U-Th ages of samples below the 72 mm cluster reveal much younger ages than those above. This contradicts a quasi linear growth model for the stalagmite. Furthermore, the U-Th age of a sample at the left flank of CH-7 at 45 mm distance from top (dft) is much older compared to this linear age-depth relationship (Fig 4B). The relationship between the U concentration and the initial δ234U isotopy is asymptotic between these two parameters, indicating that δ234U decreases for higher U concentrations and proximity to the pelvis. The highest U concentrations, which reach 12.8 ppm, and smallest δ234U values are measured for samples that are adjacent to the pelvis (Fig 5).

Fig 5
Uranium concentrations and initial δ234U values of CH-7.

It is a matter of fact that the U concentration increases in bones post mortem to as high as 100 ppm (e.g. [44]). This is due to the “soaking up effect” of U into the bone due to diffusion (e.g. [44]). One approach to date a bone or tooth by U-series is to use this effect and to determine the U-isotopes of subsamples across the bone by application of an adsorption-diffusion (D-A) model (e.g. [45, 46]). Because of the porous structure of the pelvis, which makes it likely that pore water in the spongiosa disturbed the original D-A relationship, this approach was not adopted here. Instead, we identified the termini ante quem and post quem, by using the reliable closed system U/Th ages of carbonate above and below the pelvis from the overgrowing stalagmite, respectively. While the stalagmite base around the pelvic bone fragment is characterized by a complex morphological texture, as well as age which complicate the determination of the terminus post quem and the use of these ages, determination of the terminus ante quem is possible from ages above the pelvis. These show a linear age-depth model and likely no influence on the U- and Th-system (Fig 5).

The exceedingly porous carbonate texture adjacent to and below the pelvis is likely influenced by U diffusion, as indicated by extremely high U concentrations. These high U concentrations result from an opening of the uptake series system, for which the source of U is unlikely seawater because the δ234U values stay in close ranges. Hence, it could be that U is redistributed from the very local environment resulting in variable ages, such as the very young age of the sample in direct contact with the bone. Most likely, this effect is restricted to the lowermost section of the stalagmite below 72 mm (dft). We cannot exclude the possibility that samples above 72 mm (dft) are also influenced by U-exchange, but the low variability and quasi linear age-depth relationship suggests that this effect is likely negligible and within ranges of standard age-uncertainty. Based on U/Th ages above 72 mm of the CH-7 stalagmite the minimum terminus ante quem of the pelvis is 11311±370 y BP.

Regarding the significantly older age of the sample at the flank of CH-7, one must assume U loss possibly due to the mineralogical change from presumably aragonite to calcite, as is indicated by the needle like texture and present day calcite configuration.

Stable isotope analysis

Stable oxygen and carbon isotope values, δ18O and δ13C (expressed in the δ-notation relative to V-PDB), vary from -3.59‰ to -7.47‰ and from -5.22‰ to -11.61‰, respectively (Fig 6; S1 Table).

Fig 6
Stable isotope profile.

A significant change in the δ18O profile of CH-7 occurs in the lower part of CH-7 between 100 mm and 70 mm dft of CH-7. The δ18O values vary by as much as 3.87‰ and show a w-shape in this section. Above 70 mm dft, the δ18O profile of CH-7 is more variable but a trend towards higher values is identified towards the top of CH-7 being -5.41‰ on average, with a 1-sigma standard deviation of 0.69‰. Internal δ13C variability is less compared to that of the δ18O. Heaviest δ13C values occur at the bottom and top-most isotope tracks, with lighter δ13C values in between.

Age assessment

Our Chan Hol (CH-7) isotope record indicates a very pronounced 4‰ shift across tracks c and d (Fig 6), from the most positive (-3.5‰) δ18O values at 95 mm to the most negative (-7.5‰) δ18O values at 80 mm and back to -4‰ at around 11.3 ky BP, the oldest open system U-Th date. This near 4‰ excursion presents a significant shift in the isotopic signal.

Our oldest open system U/Th date of 11311±370 y BP coincides with the end of the Younger Dryas (YD), a time episode characterized by a brief return to near glacial conditions interrupting the general amelioration of climate conditions at the last deglaciation [4750]. The ~ 4‰ shift in our CH-7 δ18O record occurs across a 21 mm interval below this last U/Th date, and hence potentially falls into the time interval of the YD. Seen the amplitude and signature of the CH-7 δ18O signal, it seems likely that the YD, or part of the YD time interval, is displayed in our Chan Hol speleothem. As we do not have an absolute (U/Th) date in this lower part of the speleothem we can only, as a first guess, assume a continuous growth rate for the CH-7 stalagmite in the interval below our last independent U/Th date and linearly extrapolate our age model back in time. We can then compare our stable isotope record of CH-7 to other well-known independently dated climate archives from a similar time interval and area.

Unfortunately, speleothem records covering the YD interval are rare from the wider Caribbean region. The closest records come from New Mexico and Arizona (Fig 7A) and have been discussed in detail by [5153]. These precipitation sensitive speleothem records have been U/Th dated and interpreted to reflect changes in the intensity of the North American monsoon region during the YD time interval, in concert with global climate as recorded in Greenland Ice cores [54, 55], Asian speleothems (Fig 7B; [56]), or Cariaco Basin (off Venezuela) Ti % (Fig 7C; [57]). In terms of signature, absolute values and amplitude our Yucatan δ18O record resembles the American speleothem records and hence reflects a similar climate signal. It further resembles the global climate signal during the YD time interval as reflected in the Yamen speleothem (China) δ18O and Caricao Basin Ti % records (Fig 7B and 7C). Therefore, we are confident that the Chan Hol stalagmite has indeed grown throughout the YD time interval.

Fig 7
Comparison of oxygen isotope and other environmental data across the Pleistocene-Holocene boundary.

Confirmation for this age assessment comes from the comparison of the CH-7 carbon isotope record with radio carbon dated Lake Peten Itza (Guatemala) magnetic susceptibility [58] (Fig 7D). Variations in the Lake Peten Itza magnetic susceptibility record reflect changes in the sediment lithology, with high values associated with clay-rich horizons and low values associated with gypsum deposits, thereby representing episodes of high and low lake levels, respectively [61].

Speleothem carbon isotope values depend on the amount of local rainfall and infiltration of vegetation cover. During the late Pleistocene, the Tulum area was dominated by ‘steppe’ and possibly looked like the Irish Burren today [62]. The dense tropical forest that dominates the present-day landscape only developed around 9000 years ago [6365]. Due to the low sea-level, local rainfall at the time of CH-7 growth would have immediately infiltrated the epikarst, with water level in the karst caves dependent on the amount of rainfall. Consequently, the CH-7 δ13C record should reflect the availability of CO213C≈-10 per mil) during the dissolution of limestone in the epikarst (δ13C≈0 per mil), with decreasing (increasing) δ13C of CH-7 reflecting lower (higher) infiltration, respectively. In this way, both the Peten Itza and CH-7 records reflect changes in water runoff into the lake and cave, respectively, with lowest run off signals during the early YD. This is also similar to the Cariaco Basin Ti % record, where low Ti percentages reflect low river runoff into the basin at that time (Fig 7C; [57]).

Although the proxy signals of both Peten Itza and CH-7 contain a certain portion of nonlinear response to vegetation cover retaining water and CO2 dissolution levels in the epikarst and lake water, respectively, the agreement between the two different records is good. Both records show coeval episodes of low lake stand (Peten Itza) and infiltration (CH-7), and high lakes stand and infiltration, respectively, across the growth interval of CH-7.

A last confirmation for our age assessment comes from comparison of the CH-7 oxygen isotope record with the detrended residual Δ14C data [5960]. This comparison, presented in Fig 7E, is inspired by [51] who showed a correlation between the δ18O variability of Pink Panther cave and the detrended residual Δ14C data, postulating a linkage of North American monsoonal precipitation and solar forcing through modulation in the Walker circulation, and the tropical Pacific Decadal Oscillation and El Niño–Southern Oscillation systems [51]. From the visual inspection, CH-7 δ18O and detrended residual Δ14C data show a remarkable similarity, with episodes of increased solar activity expressed as negative Δ14C, well correlated with positive δ18O anomalies in CH-7. Correlation is reasonably high (0.48 for inversed δ18O) considering the different nature of the proxies and the independence of the age models.

To summarize, we are confident that our age assessment is reasonably correct and that the YD time interval is indeed recorded in our CH-7 speleothem. This rises the age of the pelvis from the U/Th derived terminus ante quem of 11311 y BP to an age as old as 13 kyr BP (Fig 7F).


The Chan Hol II individual was discovered at about 1240 meters away from the nearest modern entrance, the Chan Hol sinkhole. The skeleton disarticulated slightly during the final stages of decomposition and, probably again, during the early to middle Holocene flooding of the cave, but most bones still lie close to their original anatomical position. Even small bones, like auditory ossicles, hyoid, or ungual phalanges are present. The person must thus have died in the cave at a time when the cave floor was dry [21]. The decay of the carcass occurred in-situ. Growth of the CH-7 stalagmite began after the decay of the Chan Hol II individual was complete. This interpretation is consistent with the macroscopic sequence of basal-most stalagmite laminae. At that time, the drip point was located near the margin of the pelvis (Fig 5). The calcite layer precipitated by the lateral run-off dripping water embraced the pelvic bone from laterally to ventrally, with its ventral surface coalescent with the stalagmite, due to flow extension below the pelvis. This close overgrowth could not have taken place with bones covered by soft tissue. The porous tufa-like layer conforming the basis of the stalagmite must therefore have formed at that time, when the cave was dry and the pelvis completely exposed on the cave floor for an unknown amount of time.

No data are at hand to define the amount of time that elapsed between the death of the individual and initial growth of the CH-7 stalagmite, nor the lapse needed for maceration and decay of this individual under the environmental conditions prevailing in the cave during the YD. Corpses decaying in caves are mostly decomposed by fungi and insects. Both are not able to move bones [66]. According to [66], a 25 kg kangaroo carcass, deposited in a cave in southern Australia, is completely decayed after a little more than 1000 days. For the complete decay of a human carcass with a body mass of 60 kg one would expect a minimum decay time of 2000 days as an estimate. Our U/Th datum of 11311±370 y BP at 72 mm of the CH-7 stalagmite and even the 13 ky BP age assignment resulting from the speleothem stable isotope record must therefore be regarded as minimum ages of the Chan Hol II skeleton.

Validation of our age assessment has been done by comparison of our CH-7 stable isotope data with different independently dated climate records (Fig 7). We stress, that no age correlation or wiggle matching has been carried out to any of the records. All we did was to apply the linear age model from the closed system U/Th dates to the lower part of the stalagmite.

The comparison of our CH-7 stable isotope record with other climate records indicates that the Chan Hol speleothem indeed covers the Younger Dryas time interval (Fig 7). In terms of amplitude and absolute value, the Yucatan δ18O record fits well the δ18O signal of speleothem records of New Mexico [5152] and Arizona [53] that have been demonstrated to record the global climate signal of the YD. The American speleothem δ18O records have further been interpreted to reflect changes in the contribution, intensity and source of winter versus summer precipitation, the latter being fed from the Caribbean, and these changes have been linked to changes in the positioning of the polar jet stream related to the still northerly expansion of ice sheets causing modulation of winter storm tracks across the continent [5253]. On the other hand, the relationship between the Pink Panther cave oxygen isotope record and solar forcing has been explained through changes in the Walker circulation and the Pacific Decadal Oscillation and El Niño–Southern Oscillation systems [51] but shows a significant similarity to Northern Hemisphere records. Finally, the climate signals at Peten Itza and Cariaco Basin have been discussed to reflect swings in the position of the ITCZ [57, 58].

In the end, these climate components are all linked to a complex system [67] and likely influenced our Chan Hol record. However, it is beyond the scope of this paper to disentangle the different components of this complex climate system. This deserves a thorough discussion in a separate paper. The focus of the current paper is on the dating of the Chan Hol II skeleton and we can confidentially state that with the U/Th dates and the stable isotope record at hand we can approximate the age of the Chan Hol II individual to ~13 ky BP.


Speleothem (U/Th) age data indicate that the Chan Hol underwater cave south of Tulum, state of Quintana Roo, Mexico, was accessed by humans during the Younger Dryas period, i.e. during the late Pleistocene. This is indicated by a minimum speleothem age of 11311±370 y BP of a stalagmite encrusting and overgrowing the pelvic bone of an almost articulated human individual in this cave. The age was measured at 72 mm depth from the top of the CH-7 stalagmite, at about 21 mm above the pelvis and 33 mm above the base of the stalagmite, while ages in the immediate bone vicinity are altered due to uranium dissolution. 11311±370 y BP is thus a minimum age for the skeleton. Based on a linear growth model and extension of growth rates from the well-dated upper part of the CH-7 stalagmite to its lower portion and base, the age of the Chan Hol II human rises to ~13 ky BP. The Chan Hol II skeleton thus represents one of the oldest directly dated osteological heritage of a human from the American continent. Age of the Chan Hol II human equals that of other skeletons in the Tulum cave system (e.g. Naia, Najaron), thus emphasizing the importance of these caves for early human settlement in the Americas [20, 21, 25].


The CH-7 stalagmite consists of only calcite (no aragonite was detected), as was confirmed by 25 measurements using Raman spectroscopy techniques at the Institute of Earth Sciences at Mainz University, Germany (Figs (Figs55 and and8).8). For Raman spectroscopy a Horiba Jobin Yvon was used that was connected to an Olympus BX41 microscope using a Nd-YAK laser at a wavelength of 532.12 nm (hole = 400 μm; slit = 100 μm).

Fig 8
Raman spectrum of CH-7.

U/Th dating

Samples for U/Th-dating were drilled from stalagmite CH-7 using a Dremel Fortifex precision tool, 1 mm in diameter. The samples were taken along laminar growth layers to minimize mixing of material of different age and thus age uncertainties. Individual sample thickness is typically 2 mm (in growth direction), with an individual sample weight of between 100 and 150 μg. All samples (carbonate powder) were prepared for measurements in the clean laboratories at the Institute of Earth Sciences and Institute of Environmental Physics, both Heidelberg University, by wet-column chemistry using UETVA® resin. All samples were spiked using a Th-U multi-spike. U- and Th-isotopes were analyzed using ICP-MS (Thermo Finnigan Neptune Plus and iCAP (RQ), respectively) at the Institute of Environmental Physics at Heidelberg University. Details for sample preparation and U- and Th-isotope analysis are documented in [69]. Ages were calculated using the half-lives of both elements as determined by [70]. Detrital correction was performed using a bulk Earth value of 3.8 ± 1.9. Age uncertainties are quoted at the 2-σ level and do not include half-life uncertainties. The reference year for all ages given in the study is 1950 AD (i.e. 0 BP).

U/Th ages and the growth of the CH-7 stalagmite

U/Th ages in the upper 72 mm of the CH-7 stalagmite are approximately consistent with the macroscopic sequence of individual laminae growing onto each other (Fig 4). The basal-most layer was dated to 11363±304 y BP. In this layer, the drip point (the highest point of each stalagmite layer) is located at about 20 mm lateral to the pelvis (Fig 4). The layer embraces the bone from lateral to ventral and its ventral surface is coalescent with the stalagmite. The dripping water accumulated lateral of the bone and ran off diffusing laterally, with the flow expanding below the pelvis. To do so, the body must then have been completely decayed already. In a second step, calcium carbonate-rich water, dripping from the ceiling, accumulated next to the pelvis and enclosed the bone completely. The porous tufa-like layer conforming the base of the stalagmite and dated to >11363 y BP, must therefore have formed when the cave was dry and the carcass already completely skeletonized. This is concluded from the spongy carbonate crust that formed beneath the pelvis, at a time when this bone blank of soft tissue.

C and O isotope analysis of CH-7 stalagmite

Stable oxygen and carbon isotope samples were micro-milled and measured at the Institute for Earth Sciences (GeoZentrum Nordbayern), Friedrich-Albert-Universität Erlangen, Germany. A total of 117 data points was sampled along five transects, each along the growth axis of stalagmite CH-7 (Figs (Figs44 and and6).6). A minimum of 0.05 to 0.1 mg CaCO3 was analyzed to ensure precise measurement. Carbonate powders were reacted with 100% phosphoric acid at 70°C, using a Gasbench II connected to a ThermoFisher Delta V Plus mass spectrometer All values are reported in per mil relative to V-PDB through international standard NBS19. Reproducibility was monitored by international and in house laboratory standards and was 0.5‰ and 0.8‰ for δ13C and δ18O, respectively.

Supporting information

S1 Table

Stable isotope (δ13O, δ18O) measurements of the CH-7 stalagmite from Chan Hol cave.



We acknowledge Valentina Cucchiara and Nick Poole (Liquid Jungle) and Thomas Spamberg for the use of underwater photographs of the ChanHol II site prior to the robbery of the skeleton and Ben McGeever (DiveXtras) for the permission to use diver propulsion vehicles, and we thank Carmen Rojas Sandoval and Adriana Velásquez (INAH, Tulum) for helpful discussions. Réne Eichstaedter is gratefully acknowledged for his expertise on MC-ICPMS and Th/U dating and Michael Joachimski for enabling stable isotope analysis. We gratefully acknowledge support of the project “Atlas Arqueológico Subacuático para el Registro, Estudio y Protección de los Cenotes en la Península de Yucatán” and “Estudio de los grupos humanos preceramicos de la costa oriental de Quintana Roo, México, a través de los contextos actualmente inundados” by the Instituto Nacional de Antropología e Historia (INAH). Michael Waters and one anonymous reviewer, as well as journal editor Michael Petraglia, are gratefully acknowledged for their many helpful comments and corrections to this manuscript.


The Chan Hol II osteological material, including the CH-7 stalagmite encrusting pelvic remains, is housed at the Área de Prehistoria y Evolución of the Instituto de Investigaciones Antropológicas at the Universidad Nacional Autónoma de México (UNAM), in Mexico City. Here the material is publicly deposited and accessible by others in a permanent repository. Specimen numbers: Chan Hol II human: PQR2012-CHAN HOL-2; the Chan Hol stalagmite is housed under the number CH-7. The material was recovered from the Chan Hol cave sytem at 20°9.467' N, 87°34.165' W, 15 km southwest of Tulum, Quintana Roo, southern Mexico, and about 11.5 km from the coast line (Fig 1).


All necessary permits were obtained for the described study, which complied with all relevant regulations from Instituto Nacional de Antropología e Historia, Mexico (INAH permit number: C.A. 401.B (4)19.2011/36/1723).

Funding Statement

All the funding or sources of support received during this specific study have been presented. This financial support was granted to us by the Internationales Büro of the German Bundesministerium für Bildung und Forschung (BMBF project 01DN119) and the Deutsche Forschungsgemeinschaft (DFG project STI 128/28-1). MD acknowledges support by the Irish Research Council (IRC) by a Government of Ireland Postdoctoral Fellowship (GOIPD/2015/789). BMBF and DFG financed our field work in Mexico and provided funds for laboratory work. Author Michael Deininger received a salary from the IRC during part of this study. Other than that the funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

All relevant data are within the paper and its Supporting Information files.


1. Kemp BM, Malhi RS, McDonough J, Bolnick DA, Eshleman JA, Rickards O, et al. Genetic Analysis of Early Holocene Skeletal Remains from Alaska and its Implications for the Settlement of the Americas. Am J Phys Anthropology.2007;32: 605–621. [PubMed]
2. Jenkins DL, Davis LG, Stafford TW, Campos PF, Connolly TJ, Cummings LS, et al. Geochronology, archaeological context, and DNA at the Paisley Caves In: Graf KE, Ketron CV, Waters M, editors. Paleoamerican Odyssey, Texas, A&M University, Center for the Study of the First Americans; 2013. pp. 485–510.
3. Rasmussen M, Anzick SL, Waters MR, Skoglund P, DeGiorgio M, et al. The genome of a late Pleistocene human from a Clovis burial site in western Montana. Nature. 2014;506: 225–229. doi: 10.1038/nature13025 [PMC free article] [PubMed]
4. Raghavan M, Steinrücken M, Harris K, Schiffels S, Rasmussen S, DeGiogio M, et al. Genomic evidence for the Pleistocene and recent population history of Native Americans. Science.2015;349: aab3884 doi: 10.1126/science.aab3884 [PMC free article] [PubMed]
5. Hoffecker JF, Elias SA, O'Rourke DH, Scott GR, Bigelow NH. Beringia and the global dispersal of modern humans. Evolutionary Anthropology: Issues, News, and Reviews.2016;25(2): 64–78. [PubMed]
6. Goebel T, Waters MR, O’Rourke DH. The Late Pleistocene Dispersal of Modern Humans in the Americas. Science.2008;319: 1497–1502. doi: 10.1126/science.1153569 [PubMed]
7. Stanford DJ, Bradley DA. Across Atlantic Ice: The Origin of America’s Clovis Culture. U California Press, Berkeley; 2012. 336 pp.
8. Erlandson JM, Rick TC, Braje TJ. Paleoindian seafaring, maritime technologies, and coastal foraging on California’s Channel Islands. Science.2011;33: 1181–1185. [PubMed]
9. Erlandson JM, Braje TJ, McGill KM, Graham MH. Ecology of the Kelp Highway: Did Marine Resources Facilitate Human Dispersal From Northeast Asia to the Americas? J I Coast Archaeol.2015;10(3): 392–411.
10. Dillehay TD, Ocampo C, Saavedra J, Oliveira Sawakuchi A, Vega RM, Pino M, et al. New archaeological evidence for an early human presence at Monte Verde, Chile. PLoSOne.2015;10, e0141923 (2015). [PMC free article] [PubMed]
11. Pedersen MW, Ruter A, Schweger C, Friebe H, Staff RA, Kjeldsen KK, et al. Postglacial viability and colonization in North America’s ice-free corridor. Nature.2016;537 doi: 10.1038/nature19085 [PubMed]
12. Heintzman PD, Froese D, Ives JW, Soares AE, Zazula GD, Letts B, et al. Bison phylogeography constrains dispersal and viability of the Ice Free Corridor in western Canada. Proceedings of the National Academy of Sciences.2016;113(29), 8057–8063. [PubMed]
13. Fiedel SJ. The Peopling of the New World: Present Evidence, New Theories, and Future Directions. J Archaeological Res.2000;8(1): 39–103.
14. Ferring CR. The Archaeology and Paleoecology of the Aubrey Clovis Site (41DN479) Denton County, Texas. University of North Texas; 2001; 1–276.
15. Sanchez G, Holliday VT, Gaines EP, Arroyo Cabrales J, Martinez-Taqüeña N, Kowler A, et al. Human (Clovis)-gomphothere (Cuvieronius sp.) association ~13,390 calibrated y BP in Sonora, Mexico. Proc Nac Acad Sci.2014;111: 10972–10977. [PubMed]
16. Waters MR, Stafford TW. The first Americans: a review of the evidence for the Late-Pleistocene peopling of the Americas In: Graf KE, Ketron CV, Waters M, editors. Paleoamerican Odyssey, Texas, A&M University, Center for the Study of the First Americans; 2014. pp. 541–560.
17. Joyce DJ. Pre-Clovis megafauna butchery sites in the Western Great Lakes Region In: Graf KE, Ketron CV, Waters M, editors. Paleoamerican Odyssey, Texas, A&M University, Center for the Study of the First Americans;2013. pp. 467–483.
18. Halligan JJ, Waters MR, Perrotti A, Owens IJ, Feinberg JM, Bourne MD, et al. Pre-Clovis occupation 14,550 years ago at the Page-Ladson site, Florida, and the peopling of the Americas. Sci Adv. 2016;2: e1600375 doi: 10.1126/sciadv.1600375 [PMC free article] [PubMed]
19. Lessa A, Guidon N. Osteobiographic analysis of skeleton I, Sítio Toca dos Coqueiros, Serra da Capivara National Park, Brazil, 11,060 BP: First results. American Journal of Physical Anthropology.2002;118(2): 99–110. doi: 10.1002/ajpa.10084 [PubMed]
20. González González AH, Sandoval CR, Terrazas Mata A, Benavente Sanvicente M, Stinnesbeck W, Aviles OJ, et al. The Arrival of Humans on the Yucatan Peninsula: Evidence from Submerged Caves in the State of Quintana Roo, Mexico. Current Res in the Pleistocene.2008;25: 1–24.
21. González González AH, Terrazas Mata A, Stinnesbeck W, Benavente Sanvicente M, Aviles OJ, Rojas C, et al. The First Human Settlers on the Yucatan Peninsula: Evidence from Drowned Caves in the State of Quintana Roo (South Mexico) In: Graf KE, Ketron CV, Waters M, editors. Paleoamerican Odyssey, Texas, A&M University, Center for the Study of the First Americans; 2013. pp 399–413.
22. Bueno L, Schmidt Dias A, Steele J. The late Pleistocene/early Holocene archaeological record in Brazil: A geo-referenced database. Quaternary International.2013;301: 74–93.
23. Neves WA, Hubbe M, Bernardo D, Strauss A, Araujo A, Kipnis R. Early Human Occupation of Lagoa Santa Eastern Central Brazil: Craniometric Variation of the Initial Settlers of South America In: Graf KE, Ketron CV, Waters M, editors. Paleoamerican Odyssey, Texas, A&M University, Center for the Study of the First Americans; 2013. pp. 397–412.
24. Jackson D, Méndez C, de Saint Pierre M, Aspillaga E, Politis G. Direct Dates and mtDNA of Late Pleistocene Human Skeletons from South America: A Comment on Chatters et al. (2014). PaleoAmerica.2015;1(3): 213–216.
25. Chatters JC, Kennett DJ, Asmerom Y, Kemp BM, Polyak V, Blank AN, et al. Late Pleistocene human skeleton and mtDNA link Paleoamericans and modern Native Americans. Science.2014;344(6185): 750–754. doi: 10.1126/science.1252619 [PubMed]
26. Taylor RE. Six decades of radiocarbon dating in New World Archaeology. Radiocarbon.2009;51(1): 173–212.
27. Prüfer K, Meyer M. Anthropology. Comment on "Late Pleistocene human skeleton and mtDNA link Paleoamericans and modern Native Americans". Science.2015;347: 835–835. [PubMed]
28. Weidie AE. Geology of the Yucatan Platform In: Ward WC, Weidie AE, Back W, editors. Geology and Hydrogeology of the Yucatan and Quaternary Geology of Northeastern Yucatan Peninsula, New Orleans, NOGS Publication; 1985. pp. 1–19.
29. Lefticariu M, Perry EC, Ward WC, Lefticariu L. Post-Chicxulub depositional and diagenetic history of the northwestern Yucatan Peninsula, Mexico. Sed Geol.2006;183(1–2): 51–69.
30. Ward WC. Quaternary Geology of Northeastern Yucatan Peninsula In: Ward WC, Weidie AE, Back W, editors. Geology and Hydrogeology of the Yucatan and Quaternary Geology of Northeastern Yucatan Peninsula, New Orleans, NOGS Publication; 1985. pp. 23–91.
31. Blanchon P, Shaw J. Reef drowning during the last deglaciation: Evidence for catastrophic sea-level rise and ice-sheet collapse. Geology.1995;23(1): 4–8.
32. Moseley GE, Richards DA, Smart PL, Standish CD, Hoffmann DL, ten Hove H, et al. Early-middle Holocene relative sea-level oscillation events recorded in a submerged speleothem from the Yucatan Peninsula, Mexico. The Holocene.2015;25(9): 1–11.
33. Smart PL, Beddows P, Coke J, Doerr S, Whitaker FF. Cave Development on the Caribbean coast of the Yucatan Peninsula, Quintana Roo, Mexico. Geol Soc America Spec Pap.2006;404: 105–128.
34. Grant KM, Rohling EJ, Bar-Matthews M, Ayalon A, Medina-Elizalde M, Ramsey CB, et al. Rapid coupling between ice volume and polar temperature over the past 150,000 years. Nature.2012;491(7246):744–747. [PubMed]
35. Curtis J, Hodell D, Brenner M. Climate Variability on the Yucatan Peninsula (Mexico) during the Past 3500 Years, and implications for Maya Cultural Evolution. Quat Res.1996;46: 37–47.
36. Kennett DJ, Breitenbach S, Aquino V, Asmerom Y, Awe J, Baldini J, et al. Development and Disintegration of Maya Political Systems in Response to Climate Change. Science.2012;338: 788–791. doi: 10.1126/science.1226299 [PubMed]
37. Douglas PM, Pagani M, Canuto MA, Brenner M, Hodell DA, Eglinton TI, Curtis JH. Drought, agricultural adaptation, and sociopolitical collapse in the Maya Lowlands. Proc Nac Acad Sci.2015;112(18): 5607–5612. [PubMed]
38. Marín LE. Field investigations and numerical simulation of the karstic aquifer of northwest Yucatan, Mexico. Ph.D. thesis, Northern Illinois University, DeKalb, IL; 1990. 183 p.
39. Moore YH, Stoessell RK, Easley DH. Fresh-water/sea-water relationship within a ground-water flow system, northeastern coast of the Yucatan Peninsula. Ground Water.1992,30(3): 343–350.
40. Beddows PA. Groundwater hydrology of a coastal Conduit carbonate aquifer: Caribbean Coast of the Yucatán Peninsula México. Ph.D. thesis, University of Bristol, UK; 2004.
41. Blanchon P, Eisenhauer A, Fietzke J, Liebetrau V. Rapid sea-level rise and reef back-stepping at the close of the last interglacial highstand. Nature. 2009;458: 881–884. doi: 10.1038/nature07933 [PubMed]
42. Collins SV, Reinhardt EG, Rissolo D, Chatters JC, Nava Blank A, Luna Erreguerena P. Reconstructing water level in Hoyo Negro, Quintana Roo, Mexico, implications for early Paleoamerican and faunal access. Quat Sci Rev 2015;124: 68–83.
43. Fontugne M, Shao Q, Frank N, Thil F, Guidon N, Boeda E. Cross-dating (Th/U-C-14) of calcite covering prehistoric paintings at Sierra Da Capivara National Park, Piaui, Brazil. Radiocarbon. 2013;55: 1191–1198.
44. Pike AWG, Pettitt PB. U-series Dating and Human Evolution. Rev Mineral Geochem.2003;52: 607–629.
45. Millard AR, Hedges REM. A diffusion–adsorption model of uranium uptake by archaeological bone. Geochim Cosmochim Acta.1996;60: 2139–2152.
46. Sambridge M, Grün R, Eggins S. U-series dating of bone in an open system: The diffusion-adsorption-decay model. Quat Geochron.2012;9: 42–53.
47. Haug GH, Gunther D, Peterson LC, Sigman DM, Hughen KA, Aeschlimann B. Climate and the collapse of Maya civilization. Science.2003;299: 1731–1735. doi: 10.1126/science.1080444 [PubMed]
48. Dansgaard W. Stable isotopes in precipitation. Tellus.1964;16: 436–468.
49. Stuiver M, Grootes PM, Braziunas TF. The GISP2 δ18O climate record of the past 16,500 years and the role of the sun, ocean and volcanoes. Quat Res.1995;44: 341–354
50. Carlson AE. The Younger Dryas Climate Event In: Elias SA, editor. The Encyclopedia of Quaternary Science, Amsterdam, Elsevier; 2013. 3, pp. 126–134.
51. Asmerom Y, Polyak V, Burns S, Rasmussen J. Solar forcing of Holocene climate: New insights from a speleothem record, Southwestern United States. Geology.2007;35: 1–4.
52. Asmerom Y, Polyak VJ, Burns SJ. Variable winter moisture in the southwestern United States linked to rapid glacial climate shifts. Nature Geoscience.2010;3: 114–117.
53. Wagner JDM, Cole JE, Beck JW, Patchett PJ, Henderson GM, Barnett HR. Moisture variability in the southwestern United States linked to abrupt glacial climate change. Nature Geoscience.2010;3: 110–113.
54. Alley RB. The Younger Dryas cold interval as viewed from central Greenland. Quat Sci Rev.2000;19: 213–226.
55. Rasmussen SO, Andersen KK, Svensson AM, Steffensen JP, Vinther BM, Clausen HB, et al. A new Greenland ice core chronology for the last glacial termination. Jl Geophys Res.2006;111: D6.
56. Yang Y, Yuan DX, Cheng H, Zhang ML, Qin JM, Lin YS, et al. Precise dating of abrupt shifts in the Asian Monsoon during the last deglaciation based on stalagmite data from Yamen Cave, Guizhou Province, Science China. Earth Sci.2010;53: 633–641.
57. Haug GH, Hughen KA, Sigman DM, Peterson LC, Röhl U. Southward migration of the Intertropical Convergence Zone through the Holocene. Science.2001;293: 1304–1308. doi: 10.1126/science.1059725 [PubMed]
58. Escobar J, Hodell DA, Mark Brenner M, Curtis JH, Gilli A, Mueller AD, et al. 43-ka record of paleoenvironmental change in the Central American lowlands inferred from stable isotopes of lacustrine ostracods. Quat Sci Rev.2012;37: 92–104.
59. Reimer PJ, Baillie MGL, Bard E, Bayliss A, Beck JW, Bertrand C, et al. Radiocarbon.2004;46:1029–1058.
60. IntCal04. Calibration Issue of Radiocarbon (Volume 46, 3, 2004)
61. Hodell DA, Anselmetti FS, Arizegui D, Brenner M, Curtis JH, Gilli A, et al. An 85-ka record of climate change in lowland Central America. Quat Sc Rev.2008;27: 1152–1165.
62. Feehan J. The Rocks and Landforms of the Burren In: O'Connell JW, Korff A, editors. The Book of the Burren, Tir Eolas; 2001. pp. 14–23.
63. Leyden BW, Brenner M, Hodell DA, Curtis JH. Late Pleistocene climate in the Central American lowlands In: Climate Change in Continental Isotopic Records, Swart PK, Lohmann KC, McKenzie J, Savin S, editors. Geophysical Monograph 78, American Geophysical Union, Washington, DC; 1993. pp. 165–178.
64. Leyden BW. Evidence of the Younger Dryas in Central America? Quat Sci Rev.1995;14: 833–839.
65. Brenner M, Rosenmeier MF, Hodell DA, Curtis JH. Paleolimnology of the Maya Lowlands. Long-term perspectives on the interactions among climate, environment, and humans. Ancient Mesoamerica.2002;13: 141–157.
66. Reed E. Decomposition and Disarticulation of Kangaroo Carcasses in Caves at Naracoorte, South Australia. J Taph.2009;7(4): 265–284.
67. Felis T, Scholz D, Lohmann G, Giry C, Fensterer C, Wei W, Mangini A. Control of Seasonality and Interannual to Centennial Climate Variability in the Caribbean During the Holocene Combining Coral Records, Stalagmite Records and Climate Models In: Schulz M, Paul A, editors. Integrated Analysis of Interglacial Climate Dynamics (INTERDYNAMIC), SpringerBriefs in Earth System Sciences;2015. pp. 69–75.
68. Lafuente B, Downs RT, Yang H, Stone N. The power of databases: the RRUFF project In: Armbruster T, Danisi RM, editors. Highlights in Mineralogical Crystallography, Berlin, Germany, W. De Gruyter;2015. pp. 1–30.
69. Douville E, Sallé E, Frank N, Eisele M, Pons-Branchu E, Ayrault S. Rapid and accurate U-Th dating of ancient carbonates using inductively coupled plasma-quadrupole mass spectrometry. Chem Geol.2010;272: 1–11.
70. Cheng H, Edwards RL, Hoff J, Gallup CD, Richards DA, Asmerom Y. The half-lives of uranium-234 and thorium-230. Chem Geol.2000;169: 17–33.

Articles from PLoS ONE are provided here courtesy of Public Library of Science