Phytoliths as an indicator of early modern humans plant gathering strategies, fire fuel and site occupation intensity during the Middle Stone Age at Pinnacle Point 5-6 (south coast, South Africa)

The study of plant remains in archaeological sites, along with a better understanding of the use of plants by prehistoric populations, can help us shed light on changes in survival strategies of hunter-gatherers and consequent impacts on modern human cognition, social organization, and technology. The archaeological locality of Pinnacle Point (Mossel Bay, South Africa) includes a series of coastal caves, rock-shelters, and open-air sites with human occupations spanning the Acheulian through Middle Stone Age (MSA) and Later Stone Age (LSA). These sites have provided some of the earliest evidence for complex human behaviour and technology during the MSA. We used phytoliths—amorphous silica particles that are deposited in cells of plants—as a proxy for the reconstruction of past human plant foraging strategies on the south coast of South Africa during the Middle and Late Pleistocene, emphasizing the use and control of fire as well as other possible plant uses. We analysed sediment samples from the different occupation periods at the rock shelter Pinnacle Point 5–6 North (PP5-6N). We also present an overview of the taphonomic processes affecting phytolith preservation in this site that will be critical to conduct a more reliable interpretation of the original plant use in the rock shelter. Our study reports the first evidence of the intentional gathering and introduction into living areas of plants from the Restionaceae family by MSA hunter-gatherers inhabiting the south coast of South Africa. We suggest that humans inhabiting Pinnacle Point during short-term occupation events during Marine Isotope Stage (MIS) 5 built fast fires using mainly grasses with some wood from trees and/or shrubs for specific purposes, perhaps for shellfish cooking. With the onset of MIS 4 we observed a change in the plant gathering strategies towards the intentional and intensive exploitation of dry wood to improve, we hypothesise, combustion for heating silcrete. This human behaviour is associated with changes in stone tool technology, site occupation intensity and climate change.

The south coast of South Africa is located within the Greater Cape Floristic Region (GCFR) [22], a floral biodiversity centre with high levels of endemic plant species. It has been argued that the GCFR provided uniquely diverse resources that may have supported human populations throughout the Pleistocene, and especially during glacial periods when most terrestrial environments would have lowered productivity [23,24]. These include a triumvirate of a high variety of edible plants such as fruits and geophytes (plants with edible underground storage organs) [25,26], a rich marine ecosystem [27,28], and ungulate grazer populations occupying the now submerged Palaeo-Agulhas plain [29,30].
Since 2000, the South African Coast Paleoclimate, Paleoenvironment, Paleoecology, Paleoanthropology (SACP4) team has conducted multidisciplinary studies focused on the Mossel Bay region (Western Cape Province, South Africa) and specifically the Pinnacle Point caves and rock shelters on the central south coast of South Africa. Pinnacle Point (-34.207767, 22.085489) is a rocky headland on a cliffed coast on the Indian Ocean, approximately 10 km west of the Mossel Bay point (Fig 1). Although Pinnacle Point is a formal geographic location, where currently there is a large golf and holiday resort settled above the cliffs, we consider it as an area around which are concentrated a wide variety of archaeological, paleontological, and geological localities of scientific interest (Fig 1).
The archaeological work conducted at Pinnacle Point has focused on Pinnacle Point 13B (PP13B), with a discontinuous sequence from~170-90 ka [31], that connects with the high resolution and continuous sequence of Pinnacle Point 5-6 (PP5-6) that spans~92-49 ka [32,33]. PP5-6 includes a northern area (PP5-6N) and a southern area (PP5-6S) that we think originally connected stratigraphically, but we have yet to conduct the excavations required to determine this. PP5-6N is partially under a rock shelter and PP5-6S includes a cave. We originally separated PP5-6N into three disconnected sedimentary packets-the Northwest Remnant, the Long Section, and the Southwest Remnant. In excavations in 2017 successfully connected these three sedimentary packets and we now refer to them as PP5-6N. Together these sites make Pinnacle Point a long composite Middle Stone Age (MSA) record of human occupation spanning the early parts of Marine Isotope Stage (MIS) 6 through MIS 3. This long excavated record is rare on the South African coast.
Plants are one of the natural resources most commonly used in the daily activities of hunter-gatherer societies; as a source of food and water, as a fuel for fire, to create tools for studied, hearths associated with ochre powder production were recorded during the post-Howiesons Poort occupations [40]. At PP13B, Albert and Marean [39] interpreted the abundance of dicotyledonous leaf phytoliths during MIS 6 occupations as indicative of fires with specific properties, short-term fire activities or an unknown way of using fire in activities such it could be cooking.
We know little about the use of plants as food in the South African MSA. Charred seeds have been recovered from a few MSA South African sites, and these are Boomplaas Cave [41,42], Die Kelders Cave 1 [43], Sibudu Cave [40,[44][45][46] and Wonderkrater [47]. Geophyte remains have been found in a variety of Holocene archaeological sites from South Africa, and some examples are Melkhoutboom [48], Scott's Cave [49], Sehonghong [50] and Boomplaas [51,52], and De Hangen, Andriesgrond Cave and Diepkloof [53][54][55][56][57] in the West Coast, where evidence of consumption of Watsonia sp., Moraea sp., Babiana sp. and Hesperantha sp. (Iridaceae), and leaf sheaths of Hypoxis sp. (Hypoxidaceae), among other less common taxa, have been recorded. Deacon [48,58,59] proposed that the presence of geophyte remains in archaeological sites dating from the Holocene located in the Cape Floral Region indicated that people were aware of the high production of underground organs by geophytic plants and managed their productivity with fire. In Strathalan Cave B plant remains of Watsonia sp. were found, dating to the Terminal Pleistocene [60,61]. Despite the fact that there is one site with evidence of geophyte consumption in the MSA (Bushman Rock Shelter, see [62] various researchers have argued that geophytes would have been an important component of MSA hunter-gatherer diet [24,26,29,48,57,[63][64][65][66]).
Plant remains in archaeological sites are most commonly identified as charcoal, seeds, pollen, starches, phytoliths and non-pollen palynomorphs. Here we use phytolith analysis-microscopic particles of amorphous silica formed in plant tissues-to study the use of plants during the MSA at Pinnacle Point. Phytolith analyses applied to MSA and Middle and Upper Palaeolithic sites have proved to be a powerful tool to study fuel use in fires [39,[67][68][69][70], plant consumption [71] and the construction of bedding and living floors [72][73][74]. Due to their inorganic nature, phytoliths are resistant to most biostratinomic and post-depositional processes, including burning. Their characteristic morphologies allow for the identification of the original plant to taxonomic and anatomical level [75][76][77][78][79][80][81]. Notwithstanding the rich archaeological record for the South African MSA, only three sites have been studied through phytolith analysis to date and these are Sibudu Cave (KwaZulu-Natal) [69,82], Wonderkrater (Limpopo) [47] and PP13B at Pinnacle Point (southern Cape) [39].
This paper aims to investigate the evolution of plant exploitation strategies by early modern humans in the south coast of South Africa in the long composite Pinnacle Point sequence, spanning~120,000 years. For that purpose, we use phytolith analysis at a high temporal-spatial resolution to identify changes in human behaviour and adaptation across this long time span by adding new data from PP5-6N to the previously published data from PP13B [39].

Archaeological background
The caves and rock shelters at Pinnacle Point occur in quartzite of the Table Mountain Sandstone formation, which is overlain by calcrete formations across much of the area. All of the caves and rock shelters were eroded by ancient high sea levels. Uranium-lead dating on speleothem and thermal transfer optically stimulated luminescence (TT-OSL) on cemented marine sands provide concordant evidence that some of the caves were forming at least bỹ 1.1 Ma [83,84]. The peak of the MIS 11 high sea stand at~400 ka in this area was at~+13 msl [85] which would have washed out most of the >400 ka sediments from the caves, though some of the sites have portions that are sufficiently high to preserve older sediments.
Pinnacle Point has a long record of MSA occupation that has been excavated with high-resolution techniques including total station 3D plotting of all archaeological materials, stratigraphy, and samples, which include samples for phytolith analysis. So far, the earliest numerically dated human occupations began~170-160 ka at PP13B and this is the only coastal occupation dated to MIS 6 on the south coast of South Africa. The record at PP13B extends to~92 ka, at which time it is closed from further occupation by a dune [86,87]. The Pinnacle Point record continues then at PP5-6N, in a sequence that forms on top of a dune, dated~92 ka and continues until~49 ka [32,33].
The excavations at PP5-6N reveal a~15 vertical meter excavated section spanning~40 meters long horizontally. Our age estimates provided here are from a Bayesian age model calculated from 90 OSL age estimates, the identification and age of the Younger Toba Tuff via cryptotephra analysis (87), and U-Th ages of speleothem intercalated with sediments from a range of sites at PP (31,134). The ranges mostly reflect the 1 sigma range around the oldest and youngest OSL ages rounded to a thousand years, though some estimates are more precise where other techniques (such as cryptotephra correlation to Toba) offers more precision. The sequence documents a turnover from aeolian-dominated sedimentation (YBS,~92 ka) to roofspall-dominated sedimentation (YBSR and LBSR,~92-75 ka), to aeolian-dominated sedimentation (ALBS to OBS2,~74-60 ka), back to roofspall-dominated sedimentation (BAS and BBCSR,~63-51 ka), with a final aeolian sedimentary event at~52-49 ka [32,33]. Occupation intensities rise with the glacial MIS 4 (SADBS) [32], while at the same time there was a shift to intensive use of heat-treated silcrete [4] and the world's earliest known microlithic technology [21]. Evidence for mollusc exploitation fluctuates through the sequence in sync with sea level and coastline changes [88].

Materials
A total of 206 sediment samples were collected from PP5-6N belonging mainly to combustion features and other general anthropogenic layers. In addition, 23 control samples were collected from different geogenic layers in order to establish whether the plant phytolith remains found in archaeologically-rich layers represent anthropogenic input. Samples are stored at the Laboratory of Archaeology of the University of Barcelona (Spain), except for samples from the BBCSR and BAS Stratigraphic Aggregates (hereafter-StratAggs), which are stored at the Laboratory of Palaeobotany of the Evolutionary Studies Institute of the University of Witwatersrand (South Africa). Fig 2 shows the StratAggs of PP5-6N -representing a generally homogeneous set of sedimentary processes (modified after 32)-and the location of the 229 phytolith samples and concentration of phytoliths in samples.
Samples from each combustion feature were collected from three different layers within the combustion feature, which are from top to bottom the white, black and red layer. Samples from above and outside the identified hearths were also collected whenever possible as control samples. The red layer is understood as the sample from below hearths. This layer is produced by the reddening of the underlying sediments when the combustion feature is active [68,89]. Fig 3 shows the schematic drawing of an ideal combustion feature and layers. This sampling procedure was mainly possible in the YBSR and LBSR StratAggs where practically intact combustion features were preserved as the result of high and relatively constant rate of geogenic input, particularly roofspall, together with low rates of anthropogenic activity, such as trampling, that prevented their destruction [32]. Conversely, in the SADBS, complete single combustion features were not easily discernible since this StratAgg consists of a thick sequence of overlapping trampled combustion microfacies (microscopic facies) with a few isolated lenses of intact combustion features identified mostly microscopically [32].

Sampling methods
PP5-6N has been excavated since 2006. Excavation permits were issued by Heritage Western Cape and landowner permission was granted by Pinnacle Point Estate. Excavation methods are described in Dibble et al. [90], Marean et al. [31], Bernatchez and Marean [91] Oestmo and Marean [92], and Fisher et al [93]. All sediment samples recovered for phytolith analysis were plotted directly with a total station in three dimensions. A bar-code number is given to every sample plotted. This is the sample number used for sample identification.

Phytolith extraction
The phytolith extraction from the archaeological sediments was carried out at the Laboratory of Archaeology of the University of Barcelona (Spain) and partially at the Laboratory of Palaeobotany of the Evolutionary Studies Institute of the University of Witwatersrand (South Africa) following the fast extraction procedure of Katz et al. [94]. An initial sediment weight of between 30 and 50 mg was required. Carbonate minerals were dissolved adding 50 μl of hydrochloric acid (6 N HCl). After the bubbling ceased, 450 μl of 2.4 g/ml sodium polytungstate solution Na 6 (H 2 W 12 O 40 )ÁH 2 O] was added. The tube was vortexed, sonicated and centrifuged for 5 min at 5000 rpm (MiniSpin plus, Eppendorf). The supernatant was subsequently removed to a new 0.5 ml centrifuge tube and vortexed. For examination under the optical microscope, an aliquot of 50 μl of the supernatant was placed on a microscope slide and covered with a 24 x 24 mm cover-slip. Quantification of the total phytoliths was based on 20 fields at 200x magnification whereas morphological identification of phytoliths took place at 400x magnification using an optical microscope (Olympus BX41). A minimum of 200 phytoliths were counted for the morphological analysis and when this was not possible only those samples with a minimum number of 50 phytoliths were analysed in order to obtain as much information as possible [95]. However, the error margin in the interpretation of the phytolith representation when using a minimum of 50 phytoliths is very high (40%, Albert and Weiner, 2001).
To determine the degree of preservation of the phytolith assemblages at PP5-6N we used non-parametric Spearman's correlations to measure the strength of the association and the direction of the relationship, in association with a p-value computation, between the phytolith concentration (per gram of sediment) and three taphonomic indicators, i) the total number of morphotypes identified [112,113], ii) the percentage of weathered morphologies-non-recognizable phytoliths with signs of chemical dissolution and/or mechanical damaged-and iii) the percentage of fragile morphologies-which are the first ones to disappear under certain postdepositional processes-such as hair cells, stomata, papillae [114,115], tracheids, epidermal cells (articulated phytoliths) and parenchyma strands from dicotyledonous plants (hereafterdicots) (e.g, [76,77]) for each StratAgg.
All analyses were computed using JMP-SAS 13.2.1 software.

Mineralogical analysis
Fourier Transform Infrared Spectroscopy (FT-IR) was used to identify the bulk mineral components of the archaeological sediments in order to, 1) understand the conditions that may have affected the state of preservation of phytoliths and, 2) further shed light on fire use and site occupation patterns. Infrared spectra were obtained using KBr pellets at 4 cm-1 resolution with a Nicolet iS5 spectrometer. In order to assess the origin of the calcite, we have applied the infrared grinding curve method developed by Regev et al. [116] based on the measurement of the ratio of v 2 /v 4 heights (1420 cm -1 and 713 cm -1 , respectively) normalized to a v 3 height (874 cm -1 ). Clays exposed to high temperatures were identified using changes in the specific absorptions bands of the clay spectrum [117]. Table 1 lists the sixty-three samples with a minimum number of recognizable phytolith morphotypes (>50), together with their stratigraphic location and description, and the main phytolith and mineralogical results. The description of samples with insufficient number of identifiable phytoliths is given in Supplementary Material (S1 Table). The phytolith morphotypes identified, taxonomic association and its frequencies in samples from the different Stra-tAggs are listed in S2 Table. Sixty-two phytolith morphotypes were identified [see phytolith morphotype descriptions in Esteban [118], which were later grouped by plant types and plant parts into twelve general categories: grasses (Poaceae), restios (Restionaceae), sedges (Cyperaceae), palms (Arecaceae), leaves, wood/bark and fruits of dicots, spheroids (non-decorated margins), stomata, elongates with and without decorated margins, and irregular and indeterminate morphologies. Note that the term dicot was used for all non-monocotyledonous (hereafter-monocots) angiosperms because based on phytoliths it is difficult to distinguish between early-diverging angiosperms and eudicots. Despite non-decorated spheroid morphologies have been typically associated to the wood/bark of dicot and other non-flowering plants (i.e. Gymnosperms) (e.g., [95,109,119,120]), they were grouped separately as they constitute an important component also in restios [96,97].

Mineralogy
FT-IR analyses indicated that samples were composed mainly of clay, quartz, calcite, aragonite and dahllite (carbonated hydroxyapatite) in different proportions (Table 1 and S1 Table). Nitrate salts were detected in some of the samples (Table 1 and S1 Table). The presence of nitrate minerals, identified in spectra through the sharp peak at 1,384 −1 , might come from the precipitation due to the evaporation of groundwater at the sediment surface [121]. The main differences in the mineralogical composition were related to the stratigraphic location of samples as well as to the sample types (sediment colours and hearth facies). In the lower StratAggs (YBSR and LBSR), clay, quartz and calcite were the main minerals identified, whereas dahllite and aragonite were barely present, as opposed to ALBS and SADBS. In the ALBS, quartz and clay predominated and calcite was mostly present in the white layers. Traces of dahllite and aragonite were also identified in some samples but absorption peaks were low. Samples from the SADBS (both at the and inside the dripline) differed from other StratAggs by showing the highest absorption peaks of calcite and aragonite. Dahllite was also detected in all the samples List of the sixty-three samples with sufficient number of recognizable phytoliths to be interpreted, together with their stratigraphic location and description, and the main phytolith and mineralogical results: total number of phytoliths morphologically identified, relative number of phytoliths per gram of sediment (/g sed), percentage of weathered morphologies, diatoms and sponge spicules and FTIR results. WM = weathered morphologies. Arg, aragonite. Cal, calcite. Cl, clay (b = burned), (nb = not from this StratAgg (S1 Table). The origin of aragonite in archaeological sites usually derives from land snails or marine molluscs [121], such as mussels. Nonetheless, aragonite can also be formed from aqueous solution at high-temperatures [122]. Fragments of seashells, recovered from those sediment samples where aragonite was detected, were also analysed through FT-IR in order to assess the origin of the aragonite. The majority of the seashells analysed preserved aragonite, suggesting that the most probable source of aragonite in the sediment samples is marine molluscs.  (Table 1 and S1 Table). Aragonite and aragonite transforming into calcite was also identified in a few samples from inside the dripline area (S1 Table).

Phytolith concentration
Phytoliths were detected along all the samples of the PP5-6N sequence but with different concentrations (Table 1 and Fig 2). The 23 control samples collected from geogenic/sand dune layers had few or no phytoliths independently of their provenance in the sequence, and none of them reached the minimum number of phytoliths necessary to conduct a reliable interpretation of the data (S1 Table). All the samples containing enough phytoliths for a reliable morphological interpretation belonged to combustion features. Phytolith concentrations varied from the lowermost StratAggs (YBSR, LBSR, ALBS and SADBS), with LBSR and YBSR showing the highest concentration, to the uppermost StratAggs at the dripline area (OBS2 and OBS1, SGS, DBCS, BBCSR, RBSR), which were sterile with the only exception of sample 418155 from BBCSR (Table 1 and S1 Table). A new area excavated in 2017 inside the dripline at the uppermost StratAggs (BAS and BBCSR) showed the presence of phytoliths in high concentrations. Thus, the PP5-6N phytolith sequence crosses the inter-glacial to glacial transition of MIS 5 to MIS 4 where there are important documented changes in stone tools and site occupation [4,21,32], and the MIS4 to MIS3 transition that is at this time less well-known in southern Africa, though very well preserved at PP5-6N and the subject of future studies and publication.

Phytolith preservation
Phytoliths with signs of weathering (weathered morphologies) were identified in all the samples in moderate frequencies with some exceptions as described below (Table 1). Because phytoliths mostly come from combustion features the presence of these weathered morphologies might be partly related to the effect of fire [114]. Moreover, fragile morphologies were also preserved in most of the samples (Table 2). Strikingly, samples 418155 (55%) from the BBCSR, and 162778 (22.5%) and 356463 (20.8%) from the LBSR presented high frequencies of weathered morphologies, what might be indicative of some sort of post-depositional processes affecting the phytolith assemblage, along with fragile morphologies, what might be indicative of good preservation conditions. The exceptions were samples 356487, 356475 and 162750 from the LBSR, which did not have fragile morphologies. The latter come from an area that was affected by water processes and some of the combustion features show signs of decalcification [32]. Samples from the ALBS and SADBS had low frequencies of weathered morphologies, with the exception of a white layer (356487) from SADBS, with 72.88% of weathered morphologies. In this sample the only phytoliths morphologically identified corresponded to non-fragile morphologies such as elongates, parallelepiped blockys, irregular morphologies and grass silica short cell (GSSC) rondels which are resistant morphotypes [115]. Restio phytoliths were also present, which are also suspected to be resistant morphologies [96,97]. We used Spearman's correlation coefficients analysis to measure the strength of the association between the phytolith concentration and three taphonomic indicators (% weathered morphotypes, % delicate morphologies and number of morphotypes identified), so as to investigate the degree of preservation of the phytolith assemblages at PP5-6N. No significant p-values were recorded for most of the measurements implying that there is inconclusive evidence about the significance of the association between the variables ( Table 2). In other words, we cannot conclude that the correlation is different from 0. Correlation coefficients and p-values for ALBS and ALBS inside the dripline are suspected of error due to the small sample size.

Combustion features: Lateral and vertical variation
Black layers had the highest phytolith concentration, followed by white layers (Table 1). Conversely, red layers had the lowest phytolith concentration (Table 1). We expect this since red layers represent heated zones below hearths. Samples from above and outside hearths showed the lowest phytolith concentration.
Calcite was the dominant mineral component in the white layers, clay and calcite in the black layers, and samples from the red layers and above hearths contained little or no calcite, being mostly quartz and clay (Table 1 and S1 Table). FT-IR also showed traces of burned clay in most of the samples from white and black layers (Fig 4A and Table 1), and some of them at high temperatures ( Fig 4B and Table 1) (following Berna et al., 2007). Samples from white and black layers also showed high absorption peaks of calcite (Fig 4C) which derives from woodash (following Regev et al., 2010). Wood-ash calcite was an important mineral component only in two out of the nine samples analysed from outside hearths. Phytoliths were found in low concentrations in these two mentioned samples (S1 Table). Similarly, some samples from the red layers contained calcite as the main mineral component. The presence of high calcite concentration in some of the samples from outside hearths and from the red layers might be explained by the reworking of the sediments with the black and white layers. Plant gathering strategies during the South African MSA Samples from white layers had the highest frequencies of irregular morphotypes (phytoliths that cannot be geometrically or taxonomically/anatomically described) and the lowest frequencies of grass phytoliths (S2 Table). Samples from black layers had high frequencies of restio phytoliths in comparison with other hearth facies.

The stratigraphic and temporal variation in phytolith morphological distribution along PP5-6N
Here we describe the distribution of phytolith morphologies in respect to their stratigraphic location along the PP5-6N sequence.  Table). Spheroid echinates are known to come from different plant families but only two occur in the GCFR and these are Strelitziaceae and Orchidaceae. Although Arecaceae plants produce the highest numbers of spheroid echinate phytoliths in comparison to the other subfamilies, this plant family does not occur in our study area. However, there is no reference collection on these plant families from South Africa.
The phytolith assemblages from PP5-6N, indistinct of the StratAgg or the sample type, were dominated by grass characteristic morphotypes and elongates without decorated margins, followed by dicot leaf phytoliths and irregular and indeterminate morphologies (Figs 5 and 6A-6C). Among grasses, GSSCs were identified in high frequencies, mostly from the rondel type (Fig 6D-6F) (S2 Table). Dicot leaves, spheroids, blocky morphotypes, epidermal appendages, elongate with decorated margins (mostly sinuate) and restio morphotypes (Fig 6G-6L) were identified in moderate frequencies. Articulated phytoliths from epidermal cells of dicot leaves showed a high variety of shape outlines (Fig 7). Spheroid echinates (palms), sedge phytolith morphotypes, fruits and wood/bark phytoliths from dicot plants, and stomata had the lowest frequencies among all the StratAggs at PP5-6N. Statistically, sedges, dicot leaves and wood/ bark, spheroids, stomata, elongates without decorated margins, blockys and irregular/indeterminate morphologies showed significant differences among StratAggs (Fig 4 and S3 Table). It is important to note that the distribution of grass and restio phytoliths is not statistically different from random among StratAggs.
Specific details on the distribution of the phytolith assemblages (plant types and plant parts) among StratAggs from PP5-6N, based on the Kruskal-Wallis and the pairwise comparison (Dunn's) tests, are given below: • The highest frequencies of stomata and blocky morphotypes are found in the earliest occupation in the YBSR (~92-89 ka) (Dunn's pairwise comparison test-S3 and S4 Tables).
• The LBSR (~89-75 ka) differs from the other StratAggs by having high frequencies of elongates with no decorated margins and by low frequencies of sedge, dicot leaf phytoliths and blocky morphologies (Dunn's pairwise comparison test-S3 and S4 Tables). Although not statistically significant, spheroid echinate phytoliths were identified in three samples although in very low frequencies, and these were not observed in other StratAggs at PP5-6N (S3 Table).
• During the transition from MIS 5 to MIS 4, throughout the occupation of ALBS (~74ka), we did not observe statistically significant differences between the phytolith assemblages from ALBS and other StratAggs (S4 Table).
• Two sets of samples (at and inside the dripline) were collected from SADBS, dating to the beginning of MIS 4 (~73-68 ka). The phytolith assemblages varied depending on sample provenance, and this was statistically significant. At the dripline area irregular and indeterminate morphologies dominated and this constitutes the defining feature of this StratAgg (Dunn's pairwise comparison test-S3 and S4 Tables). Inside the dripline area, dicot leaves dominated the phytolith assemblages constituting the defining feature (Dunn's pairwise comparison test-S3 and S4 Tables).
• None of the samples analysed from OBS2 (~65-60 ka) contained enough phytoliths for a reliable interpretation of the data. • There were no significant differences between the phytolith assemblages of the BBCSR (MIS3,~63-51 ka), and other (Dunn's pairwise comparison test-S3 and S4 Tables).

Phytolith provenience and preservation
In order to better interpret the phytolith record, and thus the plant presence in archaeological sites, there are two factors that need to be taken into account. The first is whether phytoliths represent plants that were introduced by humans, this being even more significant in the case of rock shelters than caves because of the higher chances of receiving water and daylight necessary for the growth of plants. The second is the phytolith state of preservation. Understanding phytolith dissolution is critical for avoiding bias when interpreting the data. In relation to the former, at PP5-6N, the identification of phytoliths in anthropogenic layers (dominated by combustion features) compared to the scarcity or absence of phytoliths in geogenic layers implies that plants were intentionally introduced to the site by human action and the phytoliths underwent little or no post-depositional movement. The identification through FT-IR of wood-ash derived calcite attests to the past presence of firewood even in those samples with an absence or a low phytolith concentration, or the absence of wood/bark phytolith morphotypes. But this leads us to question whether this absence of phytoliths in some combustion features is due to post-depositional processes affecting phytolith preservation or to the use of plants that do not produce phytoliths. Cabanes and Shahack-Gross [115] showed that partial dissolution is responsible for the decrease of phytoliths in absolute quantities, and for the increase of morphotypes showing weathering and rugulate textures. This is because the stability of certain morphotypes differs even though mineral composition and solubility is the same. When sediments are chemically unstable, it should also be expected that the most fragile morphologies would be the first to disappear. Here we considered the association of the mineralogical composition of the sediments together with correlation measurements of the phytolith concentration against three taphonomic indicators, these being the total number of morphotypes identified, the percentage of weathered morphologies and the frequencies of fragile morphologies, in order to understand the state of preservation of phytoliths at PP5-6N (Table 2). No significant p-values were recorded for most of the measurements implying that there is inconclusive evidence about the significance of the association between the variables. We interpret these results as indicative of the low effect of taphonomic processes on the phytolith assemblages as a whole. These results lead us to consider the phytolith assemblages at PP5-6N as representative of the original plant phytolith input due to human intervention.
Samples from the ALBS and SADBS StratAggs contained very low phytolith concentrations and in most of the samples we did not identify the minimum number of phytoliths (>50) needed for a reliable interpretation of the plant component [95]. Aragonite was detected through FT-IR in most of the samples and this is one of the more unstable minerals that calcium carbonate forms. Its preservation in these StratAggs attests to the chemical stability of those sediments. The micromorphological study observed little chemical alteration of the sediments in the SADBS [32] and thus phytoliths should be relatively stable. Then, if dissolution processes did not take place here, the most plausible reason that explains the low phytolith concentration might relate to the type of fuel used rather than preservation and this will be a focus of interpretation below. Because aragonite was found in sediments and related to the presence of molluscs in the context of combustion features, there are two potential explanations as to why aragonite was preserved. The first explanation is the presence of low-intensity fires with a burning temperature below 400˚C, since this is the minimum temperature at which aragonite starts transforming into calcite [123]. The second explanation accounts for a behavioural trait when, once the molluscs were consumed, the shells were discarded by disposing them in the now-extinguished hearths, favouring the preservation of the aragonite. But this remains an open question for future studies.
Finally, the absence of phytoliths in samples from the black layers at the uppermost Stra-tAggs (OBS2, OBS1, SGS, DBCS and RBSR) together with the lack of calcite in sediments might be indicative of some postdepositional processes affecting mineral stability in the uppermost levels of PP5-6N. This is related to the location of these decalcified layers that are near to and partially below and outside the dripline [32]. Nonetheless, this association between absence of calcite and lack of phytoliths does not account for most of the samples from the BAS and BBCSR StratAggs, where most of the samples showed high phytolith concentration and low presence of calcite. Therefore, we suggest that phytolith preservation should be interpreted in combination with FT-IR spectra and micromorphology.

Fire uses, foraging strategies and site occupation at PP5-6N
The presence of a large morphological variety of phytoliths is indicative of the diversity of plants used by the inhabitants of the shelter. This is particularly true with dicots, where the high variety of articulated phytoliths in hearths suggests a wide-ranging selection of trees and shrubs foraged by the inhabitants of the site and of a wide variety of plant fuel used during the oldest occupation at PP5-6N (YBSR and LBSR) during the last phases of MIS 5 (~92-75 ka).
We report for the first time the identification of restios in the South African archaeological record and this occurs~90 ka at PP5-6N during the occupation of the YBSR StratAgg. The presence of restios is constant through the occupation of the site (~92 to 49 ka). The identification of restios in anthropogenic layers indicates the gathering and intentional introduction into the site of these plants by MSA hunter-gatherers inhabiting the south coast of South Africa. Restio plants are a diagnostic family of the Fynbos biome. Restios (together with other graminoids) have been, and still are, used traditionally in southern Africa for building, thatching and for the construction of brooms and sleeping mats ( [124], and references therein). After the arrival of European settlers to the Cape, restios were mainly used for building and thatching, with Thamnochortus insignis Mast. being the most widely used species because of its long culms [124,125]. Evidence for the use of restios by San people during the contact period with European settlements is not well known. In contrast, Khoe-Khoe people, who were herders of both cattle and sheep, used mats made with different graminoids for housing and safekeeping the cattle, as these can be easily dismantled and loaded onto cattle for transport when wandering in search of richer pasture areas and water sources [126]. Such ethnographic evidence for the use of restios and other graminoids by past and current southern African inhabitants raises the question as to their use by MSA populations inhabiting the south coast of South Africa. At PP5-6N, restio phytoliths were found in contexts of combustion features so their use might be related to fire purposes and maybe with the aim of controlling the temperature. As for grasses, phytoliths from restios can also be found attached to the bark of certain trees [96]. Hence, their presence in hearth contexts can also be related to some extent as the result of wood contamination, although marginal as noted by the study of the modern plants from the area [96]. At Sibudu Cave, the presence of remains of burnt bedding was reported at 77ka [127,128], and these were constructed mainly from sedges, other monocots and topped with aromatic, insecticidal and larvicidal dicot leaves (e.g. Cryptocarya woodii) [128]. Restios have similar characteristics to sedges in terms of shape, thickness and resistance, as well as in the traditional uses of the plants, which overlap considerably [124]. Taking into account the context, the ethnographic knowledge and other palaeoarchives, an additional explanation for the presence of restios in fire contexts might be that people inhabiting the south coast of South Africa during the Late Pleistocene used them as sleeping mats and these could have been placed by the fires so they could have been burned accidentally at some point, and mixed with the rest of the plant fuel, or intentionally as a way of disposing of old or infested bedding. Overall, the presence of restio phytoliths as far back as 92 ka at Pinnacle Point is indicative that South African MSA coastal foragers were aware of the distinctiveness of restio plants, which they might have been for fire purposes and possibly as mats.
The combustion features in the lowermost StratAgg deposits (mainly from the YBSR and LBSR) are generally intact, single hearths showing relatively small disturbance and trampling except in their periphery [32]. This is corroborated through phytoliths and FT-IR since woodash derived calcite was the main mineral component in only four out of the nine samples analysed from outside hearths, and only one contained high phytolith concentration, which indicates slight dispersion of ashes and burned remains. Few samples from below hearths showed high phytolith concentration and the presence of calcite, what must be explained by the reworking of the black layer and the soil occupation surface beneath the fire. This might be the reflection of a site maintenance behaviour by sweeping and/or dumping hearth deposits. However, some downward translocation of phytoliths from the black layer might also account for this. Irregular morphotypes were mainly detected in white layers from the LBSR and these have been traditionally associated with the wood and bark of trees and shrubs [95,129]. Collura and Neumann [108] noticed the presence of these irregular morphologies mainly in bark. Our plant reference collection from the GCFR showed a rather higher presence of irregular phytoliths in the wood/bark of trees and shrubs than in other plant types [96]. Because they were mainly observed in samples from white layers, which also contained the highest concentration of wood-ash derived calcite, it is plausible that the majority of them came from wood used as fuel. Conversely, the calcite concentration in most of the samples collected from the black layers is minimal and irregular morphologies also detected but in much lesser amounts. Black layers from the LBSR were characterized by a high phytolith concentration with grasses being the dominant vegetal component, while woody phytoliths (both from the leaves and wood/bark) and restio phytoliths were also detected but in lesser frequencies. Previous results on modern reference material from other Mediterranean environments (Israel and Greece) suggested that the presence of grass phytoliths in archaeological hearths may be related, at least partially, to contamination since between 30 and 50% of the phytolith assemblage detected in the wood and bark of trees belonged to grasses as these are commonly attached to the bark of trees [95,129]. The results of our modern plant reference material also showed that the wood/bark and leaves of dicot plants contained grass phytoliths as a result of contamination, but in much lesser amounts (mean: 16,5%) [96]. This study also showed that, quantitatively, the grass phytolith presence as contamination in other plants but grasses are on average 90 times lower [see Table 2 in Esteban et al. [96]]. This suggests the grass phytolith contribution in the studied hearths has minimal contamination from wood, implying the intentional introduction of grasses into the shelter by past inhabitants.
The fire experiments conducted by Albert and Cabanes [130] showed that the micro-charcoal fraction concentrated the major presence of dicot leaf and monocot (mainly from the grass family) phytoliths and this was noteworthy when the fresh fuel was used. Furthermore, the authors also showed that when fresh fuel was used, less calcite was produced. Samples from black layers, mostly from the lowermost StratAggs (YBSR and LBSR) concur with these findings. Karkanas et al. [32] interpreted the internal microstratigraphy of combustion features from the LBSR as suggestive of the existence of several combustion events but for a relatively short period of time. In this scenario, we suggest that humans inhabiting Pinnacle Point during these short-term occupation events during MIS 5 may have built fast fires and used mainly grasses with some fresh wood from trees and/or shrubs. Because these combustion features are always associated with shellfish, it is plausible that we are faced with fires built for specific purposes, perhaps for shellfish cooking. This interpretation of the data regarding cooking shellfish on fire is further reinforced by an ethnographical study conducted on the Anbarra Australian aborigines; as it resembles one behaviour used to cook shellfish underneath hearths by building fast fires with small sticks and mostly grasses which only lasted for a few minutes [131].
Changes in the mineralogical component, and to a lesser extent the vegetal component, of the sediments through the occupation of PP5-6N were noteworthy between the lower Stra-tAggs (YBSR, LBSR, and ALBS dated from~92 to~74 ka) and the SADBS (dating to~73-68 ka). The micromorphological study also showed that the occupation character of PP5-6N changed from low-intensity activities, suggesting relatively small groups and/or relatively short occupations in the lowermost StratAggs, to larger groups and/or longer occupations in the upper SADBS StratAgg [32]. The phytolith assemblages at the dripline area of the SADBS differed from the rest of the StratAggs in the higher presence of irregular morphologies probably derived, as discussed above, from wood/bark [e.g. [95,96]]. The SADBS also differed from the rest of the StratAggs in the mineralogical composition of the sediments, with the SADBS presenting the highest concentration of wood-ash derived calcite among StratAggs. Albert and Cabanes [130] showed through experimental fires that the percentage of carbonates (ashes) was higher when fires were made with dry wood. This is because dry wood produces a more complete combustion of the plant material [130]. At Sibudu Cave in KwaZulu-Natal, higher presence of calcite in hearths was also related to high amounts of wood fuel [69,132]. Phytolith production in wood from South African trees and shrubs is generally low and thus wood may contribute little to the phytolith record in archaeological deposits located at the GCFR [96]. Thus, although phytoliths were not identified in high concentrations in the SADBS, the high concentration of wood-ash derived calcite indicates the extensive presence of fire events, and hearths being continuously fed, most probably with dry wood, and this is not seen at the older occupations of the site (YBSR and LBSR). This is indicative of a different exploitation strategy of plant resources and a different way of making fires compared to older YBSR and LBSR. These changes in the vegetal and mineralogical component are associated with the changes in human behaviour detected through the cultural material beginning around 74 ka ago, during the transition from the ALBS to SADBS, when a major change in stone tool raw material from predominantly quartzite to heat-treated silcrete occurs [4,21]. Despite the coastline being further away from Pinnacle Point (14 km on average) [133], the cultural material recovered showed that people were still exploiting the coast and transporting shellfish back to PP5-6N [32]. Micromorphology showed that the SADBS deposits consisted of overlapping combustion features with some in situ fine combustion features but in general trampling and raked out hearth remains dominated, indicating intense combustion activities and also intense human occupation [32]. This frequency of combustion features was also seen as indicative of frequent visits to the site, perhaps even on a seasonal basis [32]. In this scenario, it is likely that past inhabitants of Pinnacle Point during early MIS 4 (~73-70 ka) were collecting large amounts of dry wood intentionally, as it is easier to gather and transport, and intensively, to invest in fire production in order to create the proper conditions for the continuous practice of heating silcrete and the subsequent knapping process [4]. The phytolith analysis shows that people during MIS 4 at Pinnacle Point targeted dry wood at the same time as a major change in stone tool technology [4,21] and intensifying human occupation [32] on the heels of significant climate change [134]. This may reflect advanced planning and organization of multiple activities (e.g. collection of appropriate wood, collection of silcrete, heat treatment of silcrete) to support consistent heat treatment of silcrete to produce the bladelet and microlithic technology that appears at Pinnacle Point at this time.

Comparison with Pinnacle Point 13B
There are notable differences in the vegetal component between PP13B [39] and PP5-6N, and this is suggestive of important changes in the patterns of plant foraging strategies by past hunter-gatherers inhabiting the south coast of South Africa through the Middle and Late Pleistocene. Our modern reference material from the study area [96,97] affords the opportunity to shed light on some puzzling results described by Albert and Marean [39] at PP13B. This site is a true cave and very well protected from the wind and rain with hearths located inside and within the dripline. The cave has a restricted floor space and therefore only one human band of people could fit comfortably inside. At PP13B the front and back of the cave have different preservation conditions [87], with better phytolith preservation at the back [39]. This results in a bias in the plant representation in the site because of the possible different plants used in different site areas [39]. Conversely, PP5-6N is a rock shelter and during the occupation of the site it would have been relatively exposed. It presents a much wider open area available for occupation so more people could conceivably inhabit it and/or spread themselves out more widely. Such differences might also be influenced by the mode of occupation of the two sites, shaping the strategies of the patterns of exploitation of the vegetal resources. At PP13B the vegetal component identified in some of the hearths from the back of the cave were mainly leaves from dicot plants. Noteworthy was that the identification of these dicot leaves took place during the oldest periods of occupation (middle MIS 6, DB Sand 4c StratAgg) [39]. Wood/bark phytoliths dominated the record in certain samples from the MIS 5c occupations of the Upper Roof Spall and DB Sand 3 [39] and these corresponded to periods of high-density occupation. Keeping in mind that the wood and bark of trees and shrubs from the central south coast of the GCFR produce very low numbers of phytoliths per gram of plant material [96], these high frequencies of wood/bark phytolith morphotypes is suggestive of the deliberate introduction of large quantities of wood in the cave and this is not observed at PP5-6N (Fig 8). Albert and Marean [39] interpreted the high abundance of dicot leaf phytoliths during DB Sand 4c (MIS 6) as the result of the production of fires with specific properties, including short-term activities or to a different use of activities such as cooking. Unfortunately, we have not been able to identify the taxonomic provenance of these dicot phytoliths in our plant reference material [96]. However, the high variability observed in the articulated phytoliths from dicot leaves found at both PP13B and PP5-6N might imply different plant provenance and thus a wideranging selection of trees and shrubs available and foraged by the inhabitants of Pinnacle Point through time.
With the exception of the samples from the Shelly Brown Sand (dating to MIS 5c), one common trait among the samples from PP13B deposits was the general low abundance of characteristic grass phytoliths [39]. Despite the bad preservation conditions of phytoliths in some of the areas of PP13B, Albert and Marean [39] interpreted the low grass phytolith abundance noted in well preserved samples as this family not being common in the surrounding areas of the cave. Based on the speleothem [135], macromammals [136] and strontium isotopes [30] records, we could expect the presence of extensive grasslands somewhere on the Palaeo-Agulhas Plain that was extensively exposed during MIS 6. If grasslands were present in the surrounding areas of PP13B during MIS 6, it is not detected in the phytolith record. Conversely, the high presence of dicot plants detected in the phytolith assemblage is representative of the presence of some sort of shrubby vegetation occurring in the vicinity areas of PP13B, and that this might have been the area preferred for plant foraging practices by the inhabitants of Pinnacle Point during MIS 6. Conversely, PP5-6N now shows that at Pinnacle Point people used grasses intensively in later periods. Thus, the results observed in both caves can indicate that the strategies of plant exploitation by people inhabiting Pinnacle Point changed overtime and possibly related to plant availability.
The high abundance of elongate morphologies with non-decorated margins-morphotypes that are considered to be originating from monocots-in most of the StratAggs of PP13B was ascribed to monocot plants, and related to grasses in those samples where characteristic grass phytoliths were identified [39]. At PP5-6N, although elongate morphologies with non-decorated margins were also identified in high frequencies, they are not as abundant as at PP13B (Fig 8). This is true even though phytoliths characteristic of restios and grasses are an important component of the phytolith assemblages at PP5-6N (see Fig 8) and these plants are also high producers of elongate morphologies [96]. Previous studies from Mediterranean plants [95,129] showed that elongate morphologies are common in monocots, and mostly in grasses. In the GCFR, Esteban et al. [96] analysed grasses, restios and geophytes all of them belonging to the monocot class. Geophytes lack diagnostic phytoliths although the leaves produced the highest numbers of elongate morphologies, and in particular those from the genera Moraea from the Iridaceae family, and this was statistically significant [96]. These results might point towards a monocot plant source of these elongate morphologies other than that of grasses and restios, and of its exploitation during the occupation of PP13B and mostly during the Upper Roof Spall and the Shelly Brown Sand StratAggs, dating to~120-95 ka [39] (Fig 8). Could it be possible that this higher presence of elongate morphologies was related to the presence of geophytes? The GCFR is a rich flora of nutritious and forager-accessible USOs, most of which are cormous species belonging to the Iridaceae family [137]. It has been suggested that USOs, especially those belonging to the Iridaceae family [138], of the southern Cape could have supported the carbohydrate needs of a hunter-gatherer community over most of the year [138] and could also be used as a fallback food [26]. Despite the poor recording of the ethnobotany of the San and Khoe-Khoe people, it is known that geophytes, and specially Moraea fugax, were and still are, an important resource for food and water for San people from southern Africa [124,139,140]. Although records of plant species consumed is relatively well known [25,124,[141][142][143], inferring plant-diet activities and harvesting strategies are very difficult as much of the huntergatherer tradition has been lost. Plant foods, and specifically the edible parts of tubers and corms, are usually fully consumed reducing the chances to find their remains in archaeological sites. We suggest that different species of geophytes were harvested in two different ways, by harvesting their entirety allowing the leaves to be passively transported and brought back into the cave, or by discarding the aerial parts during the processing in the field. We hypothesize that the first method of harvesting could have taken place at PP13B during MIS 5e, evinced by the high frequencies of elongate morphologies. Nonetheless, this hypothetical scenario is tentative and identifies the need for further comparative studies on highly diverse and endemic floras to increase knowledge and chances of plant identification in the fossil record, as well as the study of other archaeological sites from the southern coast of South Africa.

Conclusions
This study explored the strategies of exploitation of vegetal resources by past hunter-gatherers inhabiting Pinnacle Point during the time span of the origins of modern humans. The main results of this study shed light on the fire fuel used, the mode of occupation of PP5-6N and changes in the strategies of plant exploitation by modern humans inhabiting the south coast of South Africa from~160-49 ka.
We observed that plant phytoliths at PP5-6N were introduced into the site via anthropogenic input and that the phytolith assemblages can be considered representative of the original plant input. Even though partial dissolution might have affected the phytolith assemblages in some of the samples, this was not extensive enough to bias the interpretation of the whole data set.
This study reports the first evidence of the intentional gathering and introduction into living areas of restios by MSA hunter-gatherers inhabiting the south coast of South Africa during the Late Pleistocene. The restio phytoliths were identified in fire related contexts, and thus their presence may be related to fire purposes, or as sleeping mats that could have been accidentally burned and mixed with the rest of the plant fuel, or intentionally burned as a way of disposing of old or infested bedding. Nevertheless, the reason for the gathering and introduction of these plants into the shelter needs further research, as well as the study of phytoliths in other MSA archaeological sites for comparison.
Our study supports evidence for changes in plant foraging strategies that correlate with changes in lithic tool activities, site occupation intensity and climate change. We observed that during the last phases of MIS 5 the inhabitants of Pinnacle Point during these short-term occupation events used mainly grasses and woody plants to a lesser extent. Such a pattern is consistent with fast fires made mainly from grasses and some fresh branches from trees and/or shrubs collected from the vicinity of the site. Based on the phytolith and archaeological record and reinforced by ethnographical references we hypothesize that these combustion features were built to cook shellfish in hearths.
During MIS 4, together with the switch to technologies including regular use of heat treated silcrete and microlithic technology, we observe that people also changed plant foraging strategies by gathering large amounts of dry wood intentionally and intensively, and we infer this was done to improve combustion and to create the proper conditions for the regular practice of heating silcrete.
Further studies on contemporary archaeological South African sites are necessary to obtain a more complete picture on the uses of plants by modern human populations and their relation to the environment.  Mark Difford for general guidance with the statistics used in this work, and A. Thoms, J. Harris, E. Singels, J. De Vynck and S. Botha for helpful comments given on hunter-gatherers geophyte harvesting strategies in South Africa. We thank the two reviewers, Marco Madella and Lyn Wadley, for their constructive comments and suggested edits on the manuscript. The opinions expressed in this publication are those of the author(s) and do not necessarily reflect the views of any of the acknowledged colleagues and reviewers or funding organizations.