# Understanding anatomical plasticity of Argan wood features at local geographical scale in ecological and archaeobotanical perspectives

Jérôme Ros, Jean-Frédéric Terral, Marie-Pierre Ruas, Sarah Ivorra, Bertrand Limier, Mohammed Ater, Laure Paradis, Ahmed Ettahiri, Abdallah Fili, Jean-Pierre van Staëvel

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Jérôme Ros, Jean-Frédéric Terral, Marie-Pierre Ruas, Sarah Ivorra, Bertrand Limier, et al. Understanding anatomical plasticity of Argan wood features at local geographical scale in ecological and archaeobotanical perspectives. Scientific Reports, 2021, 11 (1), pp.10830. 10.1038/s41598-021-90286-4. http://hal-03233492

HAL Id: hal-03233492

https://hal.science/hal-03233492

Submitted on 16 Nov 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Understanding anatomical plasticity of Argan wood features at local geographical scale in ecological and archaeobotanical perspectives

Jérôme Ros<sup>1,2,9</sup>, Jean-Frédéric Terral<sup>1,2,9✉</sup>, Marie-Pierre Ruas<sup>3</sup>, Sarah Ivorra<sup>1,2</sup>, Bertrand Limier<sup>1,2,4</sup>, Mohammed Ater<sup>2,5</sup>, Laure Paradis<sup>1,2</sup>, Ahmed S. Ettahiri<sup>6</sup>, Abdallah Fili<sup>7</sup> & Jean-Pierre Van Staëvel<sup>8</sup>

The emergence of the Argan tree as an agricultural, pastoral, cultural, economic and ecological keystone species in Southern Morocco is considered to be linked to the settlement of agropastoral communities that favored its expansion. Nevertheless, the use and exploitation of Argan tree is documented by both few medieval written sources and archaeobotanical studies, from a single location, Îgîlîz (Toughmart, Morocco), a famous medieval site of the Anti-Atlas Mountains. Therefore, data remain scarce regarding the type of Argan communities exploited at this period. In order to document this question, a quantitative eco-anatomical approach aiming to understand variations of wood characters involved in sap conduction and reserve storage, is developed from modern samples collected in the area of Îgîlîz. Results show that diameter of branches and environmental factors are the major parameters explaining plasticity of wood anatomical characters. Quantitative eco-anatomical features of Argan archaeological charcoal confronted to two statistical models, allow assessing both the diameter of the branches from which it derives and the agro-ecological conditions of tree growth and development. This preliminary study may be considered as a relevant and pioneering work for the understanding of the eco-history of the Argan tree, and of its use and exploitation during the past.

The Argan tree, *Argania spinosa* (L.) Skeels, Sapotaceae ('Argân' in Amazigh Berber language), is, along with the Atlas cedar (*Cedrus atlantica* (Manetti ex Endl.) Carrière, Pinaceae) and the two olive tree subspecies (*Olea europaea* L. subsp. *europaea* and *O. e. subsp. maroccana* (Greuter & Burdet) Vargas et al., Oleaceae), the most iconic woody plants of Southwest Morocco. It is a xerophilous endemic species and the only representative of the tropical Sapotaceae family in Morocco. The Argan tree is distributed from sea level up to 1300-1500 m, from Safi (North) to the Drâa River (South) and isolated populations extend as far as Tindouf, well inside the Western Sahara (Fig. 1).

<footnote n="1">¹Institut des Sciences de l'Evolution – Montpellier, CNRS, IRD, EPHE, Univ. Montpellier, Place Eugène Bataillon, 34095 Montpellier Cedex 5, France. ²International Associated Laboratory / International Research Program (LIA/IRP, INEE-CNRS) EVOLEA, France-Morocco, Montpellier, France. ³Archéozoologie, Archéobotanique: Sociétés, Pratiques Et Environnements (AASPE), CNRS, MNHN, CP56, 43 Rue Buffon, 75005 Paris, France. ⁴INRAE, Centre Occitanie-Montpellier, 2 Place Pierre Viala, 34000 Montpellier, France. ⁵Laboratoire Botanique Appliquée, Equipe Bio-Agrodiversité, Faculté Des Sciences, Université Abdelmalek Essaâdi, BP 2060, 93 030 Tétouan, Morocco. ⁶Institut National Des Sciences de L'Archéologie Et du Patrimoine, Département d'archéologie islamique, Madinat Al-Irfane, Hay ar-Riyad, Angle Rues 5 et 7, Rabat-Instituts, 10100 Rabat, Morocco. ⁷Faculté Des Lettres Et Des Sciences Humaines, Université Chouaib Doukkali - El Jadida, 2, Avenue Med Ben Larbi Alaoui, B.P. 299, 24000 El Jadida, Morocco. ⁸UMR Orient Et Méditerranée, CNRS, EPHE, Collège de France, Université Paris, 1-Panthéon-Sorbonne, 27 Rue Paul Bert, 94200 Ivry-sur-Seine, France. ⁹These authors contributed equally: Jérôme Ros and Jean-Frédéric Terral✉email: jean-frederic.terral@umontpellier.fr</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Location of the study area and sampled Argan trees. (a, b): Maps and Argan distribution area was drawn with the help of resources in the public domain, ASTER GDEM (https://asterweb.jpl.nasa.gov/gdem.asp) and Natural Earth (https://www.naturalearthdata.com). (c, d): Geographical location of the sampled Argan trees is presented using © Esri softwares (ArcMap and ArcGlobe 10.7, Esri World Imagery). Maps were designed and drawn by Laure Paradis (ISEM 2020).</figcaption></figure>

It forms particularly well-developed woodlands on the Essaouira and Souss plains, near the Atlantic coast, and rather open vegetation in the Anti-Atlas and the High Atlas area. The Argan woodlands may be considered as a transitional macro-ecotonal zone among two main bioclimatic and phytogeographic areas: (1) the thermo-mediterranean vegetations, further north, characterized by wild olives (Olea europaea ssp. europaea var. sylvestris (Mill) Lehr), the Mediterranean fan palm (Chamaerops humilis L., Arecaceae) and the Berber Thuya (Tetraclinis articulata (Vahl) Mast., Cupressaceae) communities; (2) Sahelo-Saharan shrublands reaching in the Argan area their septentrional boundary, dominated by Acacias and succulent species such as flat crown (Vachellia gummifera (Willd.) Kyal. & Boatwr., Fabaceae), resin spurge (Euphorbia resinifera O.Berg, Euphorbiaceae) and Swizzle Sticks (Kleinia anteuphorbium (L.) Haw., Asteraceae), respectively. The climate of the Argan tree zone benefits from temperate oceanic influences with annual precipitations between 150 and 400 mm and frequent fog throughout the year.

Its origin and biogeographical history remain enigmatic. The Argan tree is considered to have a Tertiary origin, period during which the coastal areas of the Southwestern Morocco were under tropical to subtropical influence. Its Pliocene range is thus supposed to be larger than the present distribution area¹,². However, this hypothesis is not supported by any tangible proofs. At the contrary, the palynological and palaeobotanical data
of the late Pliocene and early Quaternary seem to show that the Argan tree occupied a more limited area than
today. The Argan sporadic Quaternary pollen record pleads in favor of an ancient establishment, probably prior
to the Quaternary, of its restricted endemism³⁻⁶.

A multiproxy analysis of marine sediment cores covering the last 2000 years showed that the Argan tree expanded very recently (around 150 BP) in Southeastern Morocco, while all woodland and shrubland communities continued their decline correlatively with the rapidly increasing population, pastoralism (mainly goat farming), deforestation and agriculture (olive cultivation)⁷. In his recent synthesis, Ballouche⁸ insisted on the importance of the goats that played a significant role in vegetation degradation while being an efficient agent of dispersion of the Argan tree by endozoochory. The dominance of the Argan woodland in Southwestern Morocco is thus mainly anthropogenic. Therefore, the Argan tree is the main element of a cultural landscape that may be considered as a socio-ecosystem. The tree, neither wild nor domesticated, is the mainstay of a production system for Argan oil and constitutes the major source of vegetable fat in the local traditional diet. Its wood is used as fuel and for construction. Its foliage and fruit pulp are an important source of feed for livestock⁹⁻¹². In areas close to houses and in cultivated fields, Argan trees are protected and branches are cut only when hindering the passage of inhabitants, to develop the height of the tree, to increase fruit production, or to enable ploughing with mules or facilitate fallen fruit gathering¹³.

The growth form of the Argan tree seems to be very variable according to local conditions. When managed and protected, for shade or fruit production, trees can become very large with well-developed crowns. On the contrary, when heavily browsed or pruned for fodder and fuel, they are reduced to stunted small trees and shrubs spreading on the rocky soil (Fig. 2).

Nonetheless, the history of Argan exploitation during historical times was, till recently, poorly documented but by very few medieval written sources. In the southern part of the country, excavations have been carried out at Îgiliz for the past 15 years, as a part of a French-Moroccan cooperation research program directed by A. S. Etta-hiri, A. Fili and J.-P. Van Staëvel. This program focuses on the origins of the tribal and religious movement that will found the Almohad empire (1147-1269), and aims to reconstruct the life of a peripheral Muslim mountain society in medieval and pre-modern times¹⁴⁻¹⁶. The archaeobotanical approach led on the site has revealed the first and oldest Argan remains in Morocco, resulting from the daily activities of medieval populations (exploita-tion of the Argan tree, fuel and fodder, and oil extraction). Numerous charred wood and seed shells of *Argania spinosa* predominated in a wide range of contexts, indicating the major role of this species in the tenth-thirteenth century economy. With the remains of the Argan tree, a diversity of crops (cereals, legumes and fruit trees) are attested. This rich and abundant carpological record have raised the question of techniques that were managed to maintain them in a landscape currently constrained by the semi-arid slopes¹³,¹⁷. For the present inhabitants of the neighboring village of Tifigit, an exceptional conservatory of the history and ethnobotanical observatory of pre-Saharan agrarian human communities located at 1200 m a.s.l. on the slope opposite the archaeological site, the Argan tree is an important resource as well, in domestic and subsistence activities. It is mixed with family staple crops growing in the irrigated and manured terrace soils. These practices allowed us to directly examine the management of the Argan grove and to observe its conditions of growth following different types of exploitation: irrigated manured trees in garden, non-irrigated manured trees exploited in association with crops, principally barley, unfertilized wild trees outside the village and the cultivated areas, all of them being grazed by livestock or not, depending on the will of the inhabitants (Fig. 2).

The first part of the study deals with a quantitative analytical approach on the influence of conditions of growth (irrigation, manuring, grazing) on the Argan wood anatomy, based on wood samples collected in the area of the Tifigit village and the Îgiliz archaeological site (Table S1). In a second time, comparisons are made between the modern reference model and data derived from archaeological charcoal analyses (Table S2).

The objective of this work is threefold. First, it is intended to provide a detailed guideline of all stages of the methodological and analytical procedure of quantitative eco-anatomical analysis including preliminary tests such as the determination of the number of measurements required for an optimal assessment of each anatomi-cal character and the reproducibility of measurements. Secondly, taking into consideration the age of the wood as an important factor of anatomical variability, this study attempts to assess, at the local geographical scale, the anatomical plasticity of Argan wood, in order to distinguish on the basis of quantitative anatomical features different Argan trees growing under distinct ecological conditions. Finally, the study aims to reconstruct the growth conditions and management of the Argan tree during medieval times at Îgiliz, by confronting quantita-tive anatomical data from analysis of archaeological charcoal fragments to the eco-anatomical statistical model previously obtained.

## Results

Preliminary approaches and tests (determination of the number of required measurements for an optimal assess-ment of anatomical characters and measurement reproducibility) to quantitative eco-anatomical analyses carried out on current and archaeological Argan charcoal are presented in supplementary information (Fig. S1).

**Variation of wood anatomical features according to the diameter of samples.** Quantitative data acquired and used for multivariate statistical analysis (PCA) are summarized in Table S3. The first principal components (i.e. coordinates of modern wood samples in the first axis) of the PCA are extracted from the results of the PCA. As the first principal components (PC1) of the PCA express the greatest variance while the diameter (DIA) is supposed to be the major parameter of wood character variability, the correlation among PC1 and DIA is tested.

Diameter values were previously Log-transformed (Napierian logarithm) to satisfy the requirement of nor-mality and homoscedasticity (Shapiro-Wilk normality test: W = 0.98, P value = 0.51). Correlation between PC1
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Different growing conditions and forms of the Argan tree in the study area: (1) Argan shrub-steppe; (2) Stunted small Argan tree frequently grazed by goats; (3) Argan tree in irrigated gardens; (4) Argan tree in a non-irrigated barley field; (5) Argan tree growing in a temporary stream bed ('oued') – All photographs were taken by J.-F. Terral and J. Ros.</figcaption></figure>

and sample diameter appears to be significant (Kendall correlation test: $r = -0.761$, *P* value < 0.0001). Then, variation of PC1 in relation to wood sample caliber is modelled with a least squares linear regression model (Fig. 3): PC1 = 5.73–2.27Ln(DIA). After examining the plot of residuals, this model satisfies the assumptions made prior to its establishment: non-correlation between residuals and adjusted values, sum of residuals equal to 0 and normality of residuals (*w* = 0.97, *p* < 0.24).

Since we know that PC1 varies in relation to wood sample diameter, we have developed a predictive linear regression model aiming to estimate the diameter of a branch (Fig. 4) ($r = -0.69$, *P* value < 0.0001). Then, the coordinate in the PC1 of the archaeological specimens have been integrated in the model to infer the diameter of a branch they came from. However, as it is impossible to know the position of the charred sample from which it originates, the diameter predicted by the model corresponds to a minimum diameter. This least squared linear regression model satisfies the suppositions made prior to its establishment: homoscedasticity, non-correlation
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Linear regression model explaining the variation of PC1 according to the diameter of branch (caliber) (PC1 = f(Ln(DIA))) and diagram of the standardized residuals of the model.</figcaption></figure>

between residuals and adjusted values, sum of residuals equal to 0 and normality of residuals (Shapiro-Wilk normality test: W = 0.97, P value = 0.24) (Fig. 4): Ln(DIA) = 2.53-0.24 PC1.

**Variation of wood anatomical features according to environmental parameters.** No influence of environmental parameters on wood anatomy has been demonstrated from the examination of distribution of individuals in the plane 1-2 of the PCA. As noticed above, the first axis and, at a lesser extent, the second axis, show variations related to the diameter of samples. But a trend from samples collected on trees growing in uncultivated (wild) conditions to samples from trees located within the local agroecosystem of the village of Tifigit seems to draw on the plane 1-3 (Fig. 5).

This trend appears to be explained by the axis 3 (10.1% of the variance) which seems to express variations in density of axial parenchyma (DPA) and wood ray density (DRA). Samples from trees growing in uncultivated
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Linear predictive model of the branch diameter in relation to PC1 (Ln(DIA) = f(CP1)) and diagram of the standardized residuals of the model.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 5.** PCA analysis biplot 1–3 (66.3% of variance) showing a trend from samples of Argan trees growing in wild conditions (in red) to samples from trees in cultivated conditions (in blue). The grey dots represent the samples that do not match to the expected pattern of distinction. Test samples used to test measurement errors and repeatability are shown in green. Archaeological charcoal fragments used as additional statistical individuals are represented in black.</figcaption></figure>

(wild) zones are characterized by lower values of DPA and a higher DRA than samples collected on trees growing in cultivated areas. Nevertheless, 8 individuals do not match to the pattern shown by the PCA. The possible causes of their surprising positioning will be discussed later.

**Modeling of the minimum diameter of branches from which archaeological charred samples were produced.** Coordinates in the first PCA axis (CP1) of the archaeological samples considered as additional statistical samples were applied to the predictive regression model (Fig. 4). Considering the confidence interval (at 95%), reconstruction of the minimum branch diameter was obtained by direct estimation from the predictive model (Table 1). The estimates appear relatively accurate when they are obtained by interpolation. But, it is not the case with higher values obtained by extrapolation. However, the results show that the majority of reconstituted diameters are below 50 mm and may be attributed to secondary small branches and twigs used for firewood (Gate 1-entry 1, Latrine and Fireplace). The highest modeled diameter values are those of the plank. They seem to correspond to the function and the size of the archaeological object that has been cut in one piece from the trunk of a tree. It can also be noted that the concentrated charcoal assemblages (Gate 1—entry 1 and Fireplace) gave equivalent diameter values. These charcoal fragments could come from the same branch, rather small in diameter.

**Assessment of the ecological situation of exploited medieval Argan trees.** Measurements of anatomical characters from archaeological charcoal were confronted as additional individuals to the reference model provided by the PCA. Distribution of archaeological samples in the principal plane 1–3 showed analogies with modern Argan samples (Fig. 5). The following samples (41408-7, 56613-4, 12137-3 and 52113-3) seem to possess the features of Argan growing in wild conditions. The sample number 54613-5 located in the overlapped of the two scatter plots could not be affiliated to a precise group. The other archaeological samples, characterized by a high value of axial parenchyma density (DPA) and a low wood ray density (DRA), originate from Argan trees growing in agricultural conditions.

**Discussion and conclusion**

Previous studies have shown the major importance of age (or diameter of stem, branch or trunk) on the size of wood characters¹⁸⁻²¹. Through the negative correlation between branch diameter and coordinates (PC1) of samples from the PCA carried out on the basis of eco-anatomical data explaining mostly by ‘Vessel density’
<table><caption>**Table 1.** Predicted value of branch diameter (mm) assessed using linear regression analysis.</caption><thead><tr><th rowspan="2">Stratigraphic unit</th><th rowspan="2">Archaeological context or object</th><th rowspan="2">Charcoal accession number</th><th colspan="2">Predicted value of branch diameter (mm)</th></tr><tr><th>Mean</th><th>Confidence interval (at 95%)</th></tr></thead><tbody><tr><td rowspan="5">12137</td><td rowspan="5">Gate 1-Entry 1</td><td>12137-1</td><td>27.20</td><td>[20.12-36.80]</td></tr><tr><td>12137-2</td><td>20.32</td><td>[16.83-24.54]</td></tr><tr><td>12137-3</td><td>24.43</td><td>[18.84-31.68]</td></tr><tr><td>12137-4</td><td>24.84</td><td>[19.03-32.43]</td></tr><tr><td>12137-5</td><td>22.30</td><td>[17.81-27.92]</td></tr><tr><td>41005</td><td>Plank</td><td>41005</td><td>80.05</td><td>[38.95-164.49]</td></tr><tr><td rowspan="5">52113</td><td rowspan="5">Latrine</td><td>52113-1</td><td>53.66</td><td>[30.49-94.43]</td></tr><tr><td>52113-2</td><td>43.52</td><td>[26.83-70.62]</td></tr><tr><td>52113-3</td><td>20.75</td><td>[17.04-25.25]</td></tr><tr><td>52113-4</td><td>30.21</td><td>[21.45-42.54]</td></tr><tr><td>52113-5</td><td>26.59</td><td>[19.84-35.64]</td></tr><tr><td rowspan="4">41405</td><td rowspan="4">Beam</td><td>41405-4</td><td>55.40</td><td>[31.10-98.71]</td></tr><tr><td>41405-7</td><td>30.56</td><td>[21.60-43.23]</td></tr><tr><td>41405-11</td><td>34.33</td><td>[23.20-50.81]</td></tr><tr><td>41405-13</td><td>38.82</td><td>[25.01-60.26]</td></tr><tr><td rowspan="5">54613</td><td rowspan="5">Fireplace</td><td>54613-1</td><td>24.53</td><td>[18.89-31.87]</td></tr><tr><td>54613-2</td><td>23.20</td><td>[18.25-29.49]</td></tr><tr><td>54613-3</td><td>15.64</td><td>[14.34-17.06]</td></tr><tr><td>54613-4</td><td>25.63</td><td>[19.40-33.86]</td></tr><tr><td>54613-5</td><td>13.16</td><td>[12.90-13.43]</td></tr></tbody></table>

(DVS) and 'Density of wood fenestrated zones' (DWF), our study pinpoints the relationship between anatomy and the branch age (Fig. 3). Young branches and twigs are thus characterized by a high density of vessels and fenestrated zones that, in an eco-functional perspective, may be interpreted as a higher level of sap conduction security²². In that case, the cambium functioning favors the establishment of wood with numerous and rather narrow vessels (low-porosity). In the case of a breakdown of sap conduction (cavitation), the contiguous vessels will be able to take over and allow the sap to reach the leaves at the apex of the last growth units. Moreover, the fenestrated zones are small and correlative, the cells implied in the reserve storage (horizontal and axial parenchyma) are numerous. Antagonistically, the fibers filling the fenestrated zones are fewer, thus diminishing the support function of the stem and young branch wood. In the case of branches with increasing diameter and perimeter, characterized by a conduction limited to the last growth rings, the decrease in density of vessels and fenestrated zones linked to an increase of vessel surface area, shows that the sap conduction function of wood relies more on the efficiency than on the safety. The increase of sap efficiency allows conducting the sap, more efficiently to the upper sections of the tree.

In semi-arid and arid regions, as the Igiliz and Tifgit area, water resources in pastoral land are scarce and difficult to be detected. A careful examination of the landscape before sampling the trees allowed us to detect the possible presence of water: in the areas closest to villages and associated to irrigated gardens, most often built around or on travertine edifices (calcareous formations edified by precipitation of carbonates from hard water, mainly spring water or highly mineral rivers) and the beds or banks of generally intermittent streams that grow during severe thunderstorm episodes. However, the existence of a dry watercourse does not guarantee the presence of water in the soil, since it is rapidly drained and evaporated. Unfortunately, a small part of the rainfalls is thus absorbed by the soil. In the areas further away from villages, the water resources available to plants are very uncertain. Although adapted to conditions of intense xericity, the Argan tree can have an extremely developed root system that allows it to draw water from the deepest layers of the soil, as we have seen on some individuals where a large part of the root system was uncovered by a landslide.

In our study, a variation pattern of two anatomical characters ('Density of axial parenchyma' (DPA) and 'Ray density' (DRA)) contributes to distinguish samples issued from Argan trees growing under wild conditions from those from cultivated areas. The latter are characterized by a highly developed parenchyma vertical system (parallel to the conduction system) of cells, whose exclusive function is the reserve storage, while the horizontal system (rays) ensuring not only reserve storage but also the horizontal conduction of sap among vessel groups, seems to be weaker. Therefore, it results that growing conditions (elimination of competitors, addition of organic fertilizer and irrigation) unbalance the ratio reserve storage / use of nutrient in favor of the storage. Conversely, under more stressful, wilder conditions, Argan trees rather optimize horizontal conduction to the detriment of nutrient storage, which may be attributed to a strategy of sap conduction safety. Seven samples out of 53 analyzed samples (13.2%) seem not to match with this variation model. Trees growing in the bed or banks of intermittent streams (ARG14 and ARG27 samples) seem to benefit from a substantial water supply. For this reason, anatomical features of these samples, convergent to cultivated Argan trees such as the ARG7 and ARG10 samples, are driven by unsuspected water resources. On the other hand, trees from non-irrigated cultivated
areas (barley field) on which the ARG17, ARG18 and ARG35 samples were collected, show anatomical features of wild individuals. In these cases, trees deeply rooted do not seem to take advantage of cultural practices and of the organic fertilizers (animal dung) spread on the soil of cultivated terraces. Finally, the possible outlying anatomical features of the ARG7 sample may probably be due to measurement errors, undetected during the implementation of preliminary tests.

From a bioarchaeological point of view, the minimal diameter of the branches, reconstructed on the basis of the linear prediction model applied to quantitative eco-anatomical features of archaeological samples varies from 1 to 9 cm, with a majority ranging between 1.2 and 5 cm (75% of the individuals). Although this result is to be considered with caution, since the branches may have been of higher diameter, it seems that the Argan wood used as fuel rather comes from the cutting of branches for firewood, the gathering of dead wood, or from the gathering of branches whose leaves have been previously cut for fodder.

75% of archaeological samples (N = 15) may be affiliated to Argan trees under cultivated conditions, whereas 4 samples seem deriving from individuals growing in wild conditions. Thus, it appears that Îgilîz inhabitants have mostly used wood issued from Argan trees growing and exploited in the agro-horticultural environment located in the Îgilîz mountain itself, and/or a few hundred meters below the settlement. It is very likely that the Tifigit garden and terraces set up around a spring allowing to irrigate the village, already existed during the Middle Ages. Nevertheless, in the current stage of research, we have no tangible archaeological evidence of the existence of a settlement at this location. While it was shown that the Argan tree was the basis of the daily economy and agro-sylvo-pastoral system, our results thus bring new data about the past exploitation of this emblematic fruit tree.

Tested with success on the olive<sup>23,24</sup> and the grapevine<sup>20,25</sup>, the quantitative eco-anatomical method brings for the first time results on anatomical plasticity of Argan wood, on the ecological and functional signature of anatomical features and on the medieval exploitation of Argan wood as fuel. The preliminary application of this method to 20 charred wood fragments from Îgilîz (eleventh-thirteenth century AD—Toughmart, Anti-Atlas, Morocco) allows to infer the diameter of branches used as firewood, mostly less than or equal to 5 cm. It also shows that the exploited Argan trees derive for the most part from cultivated areas (garden or cultivated terraces).

Finally, this study opens new and fascinating ecological, palaeoecological, bioarchaeological, historical, bio-geographic and palaeoclimatic perspectives, that will be examined on the basis of a larger reference collection including samples from different diameters and various ecological contexts, from the ocean coast to the continental and mountainous areas.

## Methods

**Sampling, preparation and treatment of modern reference material.** A total of 53 modern wood samples were analyzed. The modern reference samples were collected in 2014 during the annual archaeological field mission, from 36 individuals (Table S1). For some trees, two wood samples of different diameters were collected in order to take into account anatomical variability within individual.

The collected individuals showed different conditions of growth described in the introduction section and detailed in the Table 1. With the agreement of the Tifigit inhabitants and local authorities, wood sampling was achieved but samples were not collected from trunks, to avoid injuring trees of major symbolic, ecological and economic importance. Only section samples with perfect axial symmetry were retained to avoid any impact of biomechanical constraints (formation of reaction wood) on wood characters.

Once collected, the samples were air-dried during a month at the laboratory. The samples were separately wrapped in tin foil and buried in the sand and then charred without oxygen, at 450 °C for 15 to 20 min depending on the size of the sample. As a result, samples were enriched in carbon (content > 90%)<sup>20,26</sup>, reached their maximal shrinkage<sup>27</sup>, and thus are considered to become morphologically comparable to charcoal produced in medieval fires<sup>27–31</sup>. The minimum and the maximum diameter of wood samples were measured (mm) using a digital measuring calliper before and after carbonization. The diameter used in the following analyses is the mean of the two measurements carried out before carbonization.

**Archaeological material.** Twenty archaeological charcoal fragments of Argan tree identified during a previous analysis session<sup>13</sup> were included in this study (Table S2). All the Argan charcoal fragments were collected in the medieval archaeological deposits of Îgilîz<sup>15</sup>. They come from various contexts, for the most part from living units, and belong to the period of highest human activity at the site, between the late 11th and early thirteenth centuries.

**Quantitative eco-anatomical analysis of wood applied to the Argan tree.** The approach consists in measuring constitutive elements of wood and aims to understand variations according to intrinsic (inferred by the branch diameter mainly age of tree<sup>18</sup>, linked to the existence of growth rings that are often difficult to distinguish) and environmental parameters affecting the cambial activity and thus, rate of growth and wood development<sup>20,28–30</sup>. This high resolution analysis of wood structure, particularly of conductive elements, allows addressing numerous questions that have been successfully solved in the case of the olive tree and the grapevine, such as phenology, ecology, climate, impact of human activities and agricultural practices<sup>20,24,25,31–33</sup>.

*Argania spinosa* wood is diffuse-porous with a dendritic and diagonal arrangement of vessel elements in transversal section<sup>34</sup>. The axial parenchyma bands are in tangential alignment and composed of multicellular strands. In radial alignment, woody rays are 1–3 cells wide and of heterocellular composition (Fig. 6).

To apply a quantitative eco-anatomy approach to the Argan tree, both modern charred samples and archaeological charcoal are broken manually in transverse anatomical section. The following wood constitutive elements and anatomical characters related to sap conduction and reserve storage are observed and measured under a reflected-light microscope connected with an image analysis system (DFC300 FX Leica camera and LAS Leica
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Wood anatomical features of the Argan tree (in blue) and measured anatomical characters (in red).</figcaption></figure>

software) (Fig. 6): (1) vessel density (DVS—number of vessels / mm²), (2) vessel surface area (SVS, µm²), (3) ray density (DRA—number of rays / mm²), (4) axial parenchyma density (DPA, number of bands / mm²), (5) Density of wood fenestrated zones bordered on one side by the radial alignment of axial parenchyma cells and on the other side, tangentially, by rays (DWF—number of fenestrated zones / mm²).

These anatomical features were measured several times (see ‘Statistical analyses’ section) following radial lines from the cambium inwards the sample and crossing a small number of growth rings (i.e. functional rings from a sap conduction point of view). Moreover, the hydraulic conductivity or vascular conductivity (CD) was assessed using the following formula: $CD = (SVS/π)^2/DVS$ (after³²,³⁵⁻³⁷). Finally, the ratio ‘Conductive surface / total wood area’ (SC) was calculated.

**Statistical analyses.** In order to determine the number of measurements required for an optimal assessment of anatomical features, a rarefaction method was carried out from the analysis of test wood samples. For each one, repeated measurements of anatomical characters (Surface vessel area (SVS), Density of vessels (DVS), Ray density (DRA), Axial parenchyma density (DPA) and Density of wood fenestrated zones (DWF)) were performed following the aforementioned method and the cumulative mean value was then calculated for each character²⁰,²⁹. For each test sample and anatomical character, the number of measurements of a character required for an optimal assessment was quantified as the minimum number of measurements required to stabilize the mean value (rarefaction curve or cumulative mean curve).

Furthermore, different measurement sessions were carried out with the aim of testing possible errors and reproducibility of measurements taken by one or various observers, respectively. The data sets produced were tested using the PCA performed to evaluate the Argan anatomical variability. The ARG8-2 sample was used as test sample. In addition to the initial measurements. The ARG8-2 sample was analyzed 4 times: twice by one operator (ARG8-2 (1-OP1) and ARG8-2 (2-OP1)) and twice by another (ARG8-2 (3-OP2) and ARG8-2 (3-OP2)) at different times. The additional data were incorporated into the PCA as additional individuals for comparison with initial anatomical features of ARG8-2.

After showing that measurement errors have no impact on the validity of results and the measurements are reproducible, quantitative eco-anatomical data were processed using a principal component analysis (PCA) in
order to evaluate anatomical plasticity in the reference modern material, to appreciate relationships between char-
acters and wood sample caliber and to confront archaeological data to the reference model. PCA was applied on
53 reference modern samples and 7 quantitative variables (anatomical characters) to: (1) validate the hypothesis
that there is a significant relationship between the size of the branch and anatomy, as previously demonstrated by
analyses of wood development and structure¹⁸,²⁰,³⁸ and dendrochronology³⁹; (2) identify the anatomical characters
most affected by the age of the branch and, in that case, model the 'anatomical characters—caliber of the branch'
relationship; (3) develop predictive model that might estimate the minimum branch caliber from eco-anatomical
data of archaeological charcoal.

Finally, data from analysis of the 20 archaeological charcoal fragments were included in PCA as additional
statistical samples. They do not contribute to the development of the reference model, but are compared to the
modern reference samples in order to infer the ecological conditions under which Argan trees grew during the
Middle Ages.

Received: 21 October 2020; Accepted: 10 May 2021

Published online: 24 May 2021

References

1. Quézel, P. & Barbero, M. Variations climatiques du Sahara et en Afrique sèche depuis le Pliocène: enseignement de la flore et de la végétation actuelles. *Bulletin d'Ecologie* **24**(2), 191-202 (1993).

2. M'hirit, O. et al. *L'arganier, une espèce fruitière-forestière à usages multiples* (Mardaga, 1998)

3. Ballouche, A. & Carruesco, C. Evolution holocène d'un écosystème lagunaire: la lagune de Oualidia (Maroc atlantique). *Rev. Géol. Dyn. Géog. Phys.* **27**, 113-118 (1986).

4. Ballouche, A. & Damblon, F. Nouvelles données palynologiques sur la végétation holocène du Maroc. *Travaux de la section scientifique et technique. Institut français de Pondichéry* **25**, 83-90 (1988).

5. Marret, F. & Turon, J.-L. Paleohydrology and paleoclimatology off Northwest Africa during the last glacial-interglacial transition and the Holocene: palynological evidence. *Mar. Geol.* **118**(1), 107-117 (1994).

6. Leroy, S.A.G. & Dupont, L.M. Development of vegetation and continental aridity in northwestern Africa during the Late Pliocene: the pollen record of ODP 658. *Palaeogeogr., Palaeoclim., Palaeoecol.* **109**(2-4), 295-316 (1994).

7. McGregor, et al. Vegetation change, goats, and religion: a 2000-year history of land use in southern Morocco. *Quatern. Sci. Rev.* **28**(15-16), 1434-1448 (2009).

8. Ballouche, A. (2015). Le palimpseste du paysage ou comment lire le temps dans l'espace. L'exemple de l'arganeraie du sud-ouest marocain in *Traces du végétal. Nouvelles recherches sur l'imaginaire* (ed. Taibi, A. & Oghina-Pavie, C.) 23-42 (PUR, 2015).

9. Benchekroun, F. & Buttoud, G. L'arganeraie dans l'économie rurale du sud-ouest marocain. *Forêt méditerranéenne XI*(2), 127–136 (1989).

10. Charrouf, Z. & Guillaume, D. Ethnoeconomical, ethnomedical, and phytochemical Study of *Argania spinosa* (L.) Skeels. *Journal of Ethnopharmacology* **67**, 7–14 (1999).

11. Msanda, F., El Aboudi, A. & Peltier, J.-P. Biodiversité et biogéographie de l'Arganeraie marocaine. *Cahiers Agric.* **14**, 357-364 (2005).

12. Simenel, R. et al. L'argan: l'huile qui cache la forêt domestique. De la valorisation du produit à la naturalisation de l'écosystème in Les produits de terroir au service de la diversité biologique et culturelle (ed. Cormier-Salem, M.C. & B. Roussel B.) 51-74 (Armand Colin-Autrepart, 2009).

13. Ruas, M.-P. et al. History an archaeology of the emblematic Argan tree in the medieval Anti-Atlas Moutains (Morocco). *Quat. Int.* https://doi.org/10.1016/j.quaint.2015.09.030 (2015).

14. Ettahiri, A.S., Fili, A., & Van Staëvel, J.-P. La montagne d'Igiliz et le pays des Arghen (Maroc). Enquête archéologique sur une société de montagne, de la révolution almohade à la constitution des terroirs précoloniaux. *Les Nouvelles de l'archéologie* **124**, 49-53 (2011).

15. Ettahiri, A.S., Fili, A. & Van Staëvel, J.-P. Nouvelles recherches archéologiques sur la période islamique au Maroc: Fès, Aghmat et Igiliz in *Histoire et Archéologie de l'Occident musulman* (VIle-XVe siècles), *Al-Andalus, Maghreb, Sicile* (ed. P. Sénac, P.) 157-181 (Collection Méridiennes, Maison de la Recherche, 2013a).

16. Ettahiri, A.S., Fili, A. & Van Staëvel, J.-P. Nouvelles recherches archéologiques sur les origines de l'empire almohade au Maroc : les fouilles d'Igiliz. *Comptes rendus des Séances de l'Académie des Inscriptions et Belles Lettres*, 1109-1143 (2013b).

17. Ruas, M.-P. et al. Archaeobotanical research at the medieval fortified site of Igiliz (Anti-Atlas, SW Morocco) with special attention to the exploitation of the Argan tree. *Veg. Hist. Archaeobot.* https://doi.org/10.1007/s00334-011-0306-2 (2011).

18. Nicault, A., Bégin, Y. & Guiot, J. Standardisation des séries dendrochronologiques in *La Dendroécologie: Principes, méthodes et applications* (ed. Payette, S. & Filion, L.) 199-227 (Presses de l'Université Laval - Québec, 2010)

19. Achim, A. & Cloutier, A. (2010). Caractéristiques anatomiques du bois des Gymnospermes et des Angiospermes in *La Dendroécologie: Principes, méthodes et applications* (ed. Payette, S. & Filion, L.) 7-32 (Presses de l'Université Laval - Québec, 2010).

20. Limier, B. et al. Documenting the history of the grapevine and viticulture: a quantitative eco-anatomical perspective applied to modern and archaeological charcoal. *J. Archaeol. Sc.* https://doi.org/10.1016/j.jas.2018.10.001 (2018).

21. Kassout, J. et al. Resisting to aridity and aridification: adaptation of sap conduction performance in *Olea europaea* L. in a heterogeneous and changing environment (submitted).

22. Carlquist, S. & Debhur, L. Wood anatomy of Penaceae (Myrtales): comparative, phylogenetic and ecological implications. *Bot J Linn Soc* **75**, 211-227 (1977).

23. Chabal, L., Fabre, L., Terral, J.-F. & Théry-Parisot, I. L'anthracologie in *La botanique* (ed. Ferdière A.). 43-104 (Errance, 1999).

24. Durand, A. & Terral, J.-F. Regarder autrement le charbon de bois archéologique: l'exemple de l'irrigation des plantations d'oliviers en France méridionale et en Catalogne (IXe-XVe siècle). *Archéol. Mid. Médiév.* **23-24**, 75-92 (2005-2006).

25. Terral, J.-F. Quantitative anatomical criteria for discriminating wild grapevine (*Vitis vinifera* ssp. *sylvestris*) from cultivated vines (*Vitis vinifera* ssp. *vinifera). British Archaeological Reports (International Series) **1063**, 59-64 (2002).

26. Briane, D. & Doat, J. Guide technique de la carbonisation. La fabrication de charbon de bois (Edisud, 1985)

27. Beall, F. C., Blankenhorn, P. R. & Moore, G. R. Carbonized wood-physical properties and use as an SEM preparation. *Wood Sci. Technol.* **6**(3), 212-219 (1974).

28. Terral, J.-F. & Arnold-Simard, G. Beginnings of olive cultivation in Eastern Spain in relation to Holocene bioclimatic changes. *Quat. Res.* **46**, 176-185 (1996).

29. Terral, J.-F. Débuts de la domestication de l'olivier (*Olea europaea* L.) en Méditerranée nord-occidentale, mise en évidence par l'analyse morphométrique appliquée à du matériel anthracologique. *Comptes Rendus de l'Académie des Sciences de Paris* **324**(5), 417-425 (1997).
30. Marguerie, D., Bernard, V., Bégin, Y. & Terral, J.-F. Dendroanthracologie in *La Dendroécologie: Principes, méthodes et applications* (ed. Payette, S. & Filion, L.) 311–347 (Presses de l'Université Laval – Québec, 2010).

31. Terral, J.-F. Exploitation and management of the olive tree during prehistoric times in Mediterranean France and Spain. *J. Archaeol. Sc.* https://doi.org/10.1006/jasc.1999.0444 (2000).

32. Terral, J.-F. & Mengüal, X. Reconstruction of Holocene climate in southern France and eastern Spain using quantitative anatomy of olive wood and archaeological charcoal. *Palaeogeogr. Palaeoclim. Palaeoecol.* https://doi.org/10.1016/S0031-0182(99)00079-6 (1999).

33. Terral, J.-F. et al. A hydraulic conductivity model points to post-Neogene survival of the Mediterranean Olive in riparian habitat. *Ecology* https://doi.org/10.1890/03-3081 (2004).

34. Neumann, K. et al. *Woods of the Sahara and the Sahel: an anatomical atlas* (Haupt, 2000).

35. Zimmermann, M.H. *Xylem Structure and the Ascent of Sap* (Springer, 1983).

36. Carlquist, S. *Comparative Wood Anatomy. Systemic, Ecological, and Evolutionary Aspects of Dicotyledon Wood* (Springer, 1988).

37. Terral, J.-F. & Durand, A. Bio-archaeological evidence of olive tree (Olea europaea L.) irrigation during the Middle Ages in Southern France and North Eastern Spain. *J. Archaeol. Sc.* https://doi.org/10.1016/j.jas.2005.10.004 (2006).

38. Larson, P.R. *The vascular cambium. Development and structure* (Springer-Verlag, 1994).

39. Payette, S. & Filion L. *La Dendroécologie. Principes, Méthodes et Applications* (Presses de l'Université Laval, 2010).

## Acknowledgements

This research is a part of the international archaeological studies co-directed by Pr. Ahmed S. Ettahiri, Pr. Abdal-lah Fili and Pr. Jean-Pierre Van Staëvel, carried out in Îgiliz since 2009. It was supported by the HARGANA interdisciplinary project (Histoire et Archéologie des Ressources biologiques et stratégie de Gestion vivrière de l'Arganeraie médiévale en montagne Anti-atlasique (Maroc)—2014), a PRES Sorbonne-Universités Convergence Program co-directed by Pr. Jean-Pierre Van Staëvel and Dr. Marie-Pierre Ruas. This work was also supported by a postdoctoral fellowship from the HARGANA project granted to Dr. Jérôme Ros. The authors thank the Ministry of Agriculture, the Ministry of Culture and the Kingdom of Morocco, the inhabitants of Tifigit and the local authorities who allowed us to collect Argan wood samples. We thank M. Lemoine (CNRS) for his technical support. We thank both reviewers for their constructive comments. This article is the ISEM contribution no ISEM-SUD 2021-056.

## Author contributions

J.R. and J.F.T. conceived and led the writing of the manuscript; J.R., J.F.T. and S.I. conceived experiments and methodology; J.R., J.F.T., M.P.R., M.A. and S.I. collected the current material; A.S.E., A.F. and J.P.V.S. made available the bioarchaeological material; L.P. managed the geographical data; J.R. and B.L. collected the eco-anatomical data. J.R. and J.F.T. analysed the data. All co-authors provided constructive suggestions.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-90286-4.

**Correspondence** and requests for materials should be addressed to J.-F.T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021