Bayesian Morphometric Analysis for Archaeological Seed Identification: Phoenix (Arecaceae) Palms from the Canary Islands (Spain)

The taxonomic complexity of Phoenix palms in the Canary Islands, where multiple morphotypes representing at least four taxa currently exist, presents significant challenges for archaeobotanical identification. We developed a Bayesian probabilistic framework to identify archaeological Phoenix seeds within the context of genus-wide morphological diversity. Our analysis incorporated thousands of specimens including modern reference collections, archaeological materials from pre-Hispanic sites in Gran Canaria and La Gomera (3–16th centuries CE), and fossil remains. We recorded quantitative measurements and qualitative characteristics for each specimen. To understand taphonomic effects, we conducted experimental carbonization of modern P. canariensis seeds and documented the resulting morphological alterations. We performed a hierarchical cluster analysis using Ward’s minimum variance method and calculated taxonomic assignment probabilities for archaeological specimens using Bayesian inference, where likelihood was derived from taxon proportions within assigned clusters. The results indicated a high probability (0.69–1.00) that the archaeological specimens belong to P. canariensis var. canariensis, with no evidence for P. dactylifera presence. These findings provide critical insights into pre-Hispanic exploitation of Phoenix palms, particularly the endemic P. canariensis, which served as a vital resource, providing food, fiber, and construction materials. Our methodological approach offers a robust framework for addressing taxonomic uncertainty in archaeobotanical research while enhancing understanding of historical palm biogeography and resource use patterns in the Canary Islands.