INTRODUCTION
Archaeological interpretations remain incomplete without an understanding of the people who inhabited the sites (Larsen Reference Larsen2000, Reference Larsen2015), and human remains can inform on aspects of human behavior, demography, and health. Southern Africa (Africa south of the Zambezi River), with its ancient human origins and rich population diversity (Schlebusch et al. Reference Schlebusch, Malmström, Günther, Sjödin, Coutinho, Edlund, Munters, Vicente, Steyn, Soodyall and Lombard2017, Reference Schlebusch, Sjödin, Breton, Günther, Naidoo, Hollfelder, Sjöstrand, Xu, Gattepaille, Vicente and Scofield2020; Fortes-Lima et al. Reference Fortes-Lima, Burgarella, Hammarén, Eriksson, Vicente, Jolly, Semo, Gunnink, Pacchiarotti and Mundeke2023), has much to offer in this regard (e.g., Sealy Reference Sealy, Lee-Thorp and Katzenberg2016a; Pfeiffer et al. Reference Pfeiffer, Sealy, Harrington, Loftus and Maggs2020; Steyn et al. Reference Steyn, Whitelaw, Botha, Vicente, Schlebusch and Lombard2019; Loftus and Pfeiffer Reference Loftus and Pfeiffer2023; Rifkin et al. Reference Rifkin, Vikram, Alcorta, Ramond, Cowan, Jakobsson, Schlebusch and Lombard2023). Human remains represent an integral part of the population record of a region, and South Africa (the current geopolitical republic) has a long history of reporting on human skeletal remains. Initially, this was done to describe the various population groups of the world according to their perceived physical typologies, and later in a more contextual way to help understand past lifeways or aspects of disease and dental health (see Morris Reference Morris2022 for discussion).
Morris (Reference Morris1992) compiled the first inventory of Holocene human remains curated by the various South African museums and universities. He arranged the data according to the South African biomes and provided useful information about where and by whom they were found, some aspects of their archaeological context, their dates where available, any associated grave goods, skeletal parts retrieved and references. Since the publication of his catalogue, radiocarbon (14C) dating methods have been refined to allow more accurate and precise measurements from small samples, and advances have been made in the calibration protocols of radiocarbon dates (Bronk Ramsey Reference Bronk Ramsey2008; Loftus Reference Loftus2023). Consequently, there are ever greater numbers of radiocarbon dates available for regional synthetic analyses and a more widespread awareness of the methods available for interpreting radiocarbon data (e.g., Loftus et al. Reference Loftus, Mitchell and Bronk Ramsey2019; Loftus and Pfeiffer Reference Loftus and Pfeiffer2023).
Revised syntheses of both the Stone Age (Lombard et al. Reference Lombard, Wadley, Deacon, Wurz, Parsons, Mohapi, Swart and Mitchell2012, Reference Lombard, Bradfield, Caruana, Makhubela, Dusseldorp, Kramers and Wurz2022) and Iron Age (Huffman Reference Huffman2007, Reference Huffman2009) archaeological sequences have also been published based on dated assemblages (Table 1). We acknowledge that not everyone agrees fully with these syntheses (e.g., Orton Reference Orton2014; Moffett Reference Moffett2020; Wilkins Reference Wilkins2020), and that no concise overview can capture the richness and variability expressed over the more than two million years of human endeavor represented in the archaeological record of southern Africa. Yet, lacking any other comprehensive overviews, these summaries provide a pragmatic, heuristic shorthand for communicating and learning about complex data and broad spatiotemporal trends. Importantly, the synthesized sequences do not automatically imply culture historical assumptions, nor do they assume a teleological understanding of evolutionary processes (biological or socio-technical). Instead, they provide data-driven catalysts for future discussion and exploration by creating a template to situate large-scale techno-typological patterning within a chronological framework (Lombard et al. Reference Lombard, Bradfield, Caruana, Makhubela, Dusseldorp, Kramers and Wurz2022).
With this contribution we revisit the Holocene human remains of South Africa. Pre-Holocene human remains from the region are less abundant and often highly fragmentary—direct radiocarbon dates from such remains are rare, given patterns of collagen preservation in southern Africa (but see Dusseldorp et al. Reference Dusseldorp, Lombard and Wurz2013 for summary). We provide a revised georeferenced database of all the dated remains we were able to trace (SOM 1). We calibrated or recalibrated all the radiocarbon age estimates and aligned them broadly with the archaeological sequences presented in Table 1. This alignment does not automatically imply that the dated individuals were directly associated with a specific technocomplex or period, but merely that they lived during a chronological phase when most archaeological contexts show broadly similar and recognizable trends in material-culture. We use the resulting data to roughly map the distribution of Stone Age vs Iron Age human remains, record variation in the frequency of human remains and/or burials through time, and broadly associate these variations with the archaeological phases and some climatic indicators. The aim of the paper is not to interpret or explain the complete Holocene human remains record of South Africa, but to highlight possible spatiotemporal trends and other aspects or patterns that require further exploration.
MATERIALS AND METHODS
As a point of departure for our database, we extracted all the radiocarbon dated individuals from the Morris (Reference Morris1992) catalogue and conducted a systematic search for additional and subsequently published radiocarbon dates obtained for individuals in South African institutions. We also included data from unpublished reports and requested information about newly dated skeletons from various curators across the country. Here, we focus on South Africa, but future research could productively include neighbouring southern African countries, given broad similarities in the archaeological frameworks. However, we believe the numbers of dated remains from other countries are likely many fewer, given the generally lower numbers of radiocarbon dates from these countries (Loftus et al. Reference Loftus, Mitchell and Bronk Ramsey2019).
Only the remains of individuals with dates falling within the Holocene (≤11.7 ka) and before contact with European settlers since about CE 1650 were included. Each entry is accompanied by its repository accession number (identification code). Some sets of remains have been reburied in line with the wishes of descendent communities so that our inventory does not reflect all the human remains still under curation. Instead, it is a database of individuals for whom pre-contact Holocene radiocarbon dates are available. We incorporated site names into the unique ID’s created for the purposes of this study to help facilitate future sorting according to site, town, or region, and for calibration purposes indicate whether the Stone Age sites are located in the inland or along the coast (to accommodate the expected marine reservoir). Each entry has one or more associated references that report on excavation, dating, bioanthropological and/or archaeological information (SOM 1).
To gain insight into the broad distribution of the dated human remains, the site where they were found are geographically grouped according to the South African 1:50000 topographic map grid (https://shop.geospatial.com/publication/1TFM8QGY2F2M49Y8BZQ15X7Q75/South-Africa-1-to-50000-Scale-Topographic-Maps, visited September 2023). The first and second sets of two digits reflect latitude and longitude respectively, while the two-letter code ([A-D] [A-D]) reflects the particular map sheet, as divided into quadrants and sub-quadrants, labelled from top [AB] to bottom [CD]. Thus, each sheet spans an area of 15’×15’ (https://ngi.dalrrd.gov.za/index.php/what-we-do/maps-and-geospatial-information/41-sa-mapsheet-referencing, visited September 2023: see also SOM 1). Latitude and longitude are derived from this map ID, reflecting the top left corner of each map on the grid. This coarse-grained level of geographic information is suitable for broad spatial assessment on a regional or subcontinental scale as we aim to do here. More detailed spatial analyses of smaller areas will require more accurate mapping of the precise locations where the human remains were found.
Every entry in our database is associated with a unique radiocarbon ID number assigned by the laboratory that produced the date, along with the original uncalibrated radiocarbon age and measurement error. Most dates were obtained directly from an individual’s sampled bone collagen, but in some instances (∼40) dates were obtained from closely associated charcoals or other bone material—such dates mostly represent the more recent Iron Age human remains. We calibrated or recalibrated all the radiocarbon measurements, using the most recent calibration curve for southern Africa (SHCal20: Hogg et al. Reference Hogg, Heaton, Hua, Palmer, Turney, Southon, Bayliss, Blackwell, Boswijk, Ramsey, Pearson, Petchey, Reimer, Reimer and Wacker2020), and provide the calendar date ranges at the 95.4% range (both cal BP and cal BCE/CE presented, SOM 1). Due to how the calibration procedure works, the true age of the sample can lie anywhere in this range, and the entire range must be considered (see Hare and Loftus Reference Hare and Loftus2018).
Working on a large regional scale, it is necessary to consider the marine radiocarbon reservoir for variable amounts of marine protein in the diet of coastal-dwelling individuals, because dietary protein is preferentially reflected in bone collagen. Radiocarbon from marine environments is generally several hundred years “too old,” meaning that organisms that ingest substantial amounts of marine protein have older radiocarbon ages than equivalently aged organisms that acquire carbon from purely terrestrial or atmospheric sources (e.g., herbivores, plants) (Alves et al. Reference Alves, Macario, Ascough and Bronk Ramsey2018). Individualized estimates of the proportions of dietary marine protein indicate substantial variation amongst Later Stone Age foragers (e.g., Loftus and Pfeiffer Reference Loftus and Pfeiffer2023). We therefore used a simplified procedure to partially compensate for this complication by calibrating each coastal radiocarbon date for foraging people with a mixed calibration curve that incorporates between 10–50% of a marine input to the diet (MarineCal20: Heaton et al. Reference Heaton, Köhler, Butzin, Bard, Reimer, Austin, Ramsey, Grootes, Hughen, Kromer, Reimer, Adkins, Burke, Cook, Olsen and Skinner2020). (Note that no corrections were made for three dates from coastal Iron Age contexts as farming groups are believed to have consumed little, if any, marine food.) This is a conservative and overly generalized method of accounting for marine protein in the diet of coastal dwellers that will result in increased age ranges for each date. A more targeted approach would require a localized delta offset for each subregion (see Alves et al. Reference Alves, Macario, Ascough and Bronk Ramsey2018), and more individualized estimates of marine dietary intake by using, for example, stable carbon isotope ratios (Lewis and Sealy Reference Lewis and Sealy2018; Loftus and Pfeiffer Reference Loftus and Pfeiffer2023). These approaches are, however, unpragmatic and computationally unrealistic for a first attempt at a large-scale inventory of dated human remains as presented here. One of the purposes of this database and our broad-scale analyses thereof is to serve as an updated resource from which to launch more detailed investigations at centennial or sub-regional scales that employ more precise estimates of marine dietary intake and the local marine reservoir.
To evaluate the frequency of dates associated with the Holocene human remains, we used Summed Probability Density (SPD) and Kernel Density Estimation (KDE) methods. The SPD method is a widely used tool for the aggregation of radiocarbon dates, providing a curve with peaks and troughs that reflect the relative frequency of radiocarbon dates through time (see discussions and critiques of SPD approaches in Williams 2012, Carleton and Groucutt Reference Carleton and Groucutt2021, and Crema Reference Crema2022). The approach simply sums the nonnormal distributions of calibrated radiocarbon dates. The resulting distribution curve can be greatly affected by the shape of the calibration curve (e.g. “plateaus” or “jumps” which reflect major shifts in the rate of radiocarbon production through time). It is generally only considered reliable for sample numbers of >500, and care must be taken not to overinterpret minor fluctuations in the shape of the SPD curve, especially given the lack of any measure of uncertainty associated with the SPD (Williams 2012; Crema Reference Crema2022). In response to the limitations of SPD approaches, KDE methods are increasingly used, particularly in the context of open-access calibration and chronological modeling software (e.g. Bronk Ramsey Reference Bronk Ramsey2017). The KDE_Model tool in OxCal (v.4.4) generates multiple KDEs by repeatedly sampling subsets of a set of dates, and aggregating these KDEs into an ensemble, which has the advantage over an SPD of providing a visual estimate of uncertainty, as an envelope around the combined curve (Bronk Ramsey Reference Bronk Ramsey2017; Crema Reference Crema2022). The KDE method used in OxCal derives an appropriate bandwidth from the dataset itself. This tool provides smoother distribution curves than an SPD function and generally offers a more conservative estimate of the frequency of dates through time (Bronk Ramsey Reference Bronk Ramsey2017; Crema Reference Crema2022).
Well-dated, continuous palaeoenvironmental datasets were selected to provide a generalized picture of environmental shifts throughout the Holocene for comparison with the burial record. Considerable variation in precipitation trends is believed to have occurred across the subcontinent over this period, due to the different forcings on subtropical and temperate weather systems that bring moisture to the summer and winter rainfall regions (e.g., Zhao et al. Reference Zhao, Dupont, Schefuß, Meadows, Hahn and Wefer2016). In these regions, more than two-thirds of precipitation occurs in the relevant season, with a year-round rainfall region lying geographically between the two, and receiving rainfall from both systems throughout the year. We thus selected records (e.g., Talma and Vogel Reference Talma and Vogel1992; Scott et al. Reference Scott, Holmgren, Talma, Woodborne and Vogel2003; Fischer et al. Reference Fischer, Fundel, Ruth, Twarloh, Wegner, Udisti, Becagli, Castellano, Morganti and Severi2007; Chevalier and Chase Reference Chevalier and Chase2015; Chase et al. Reference Chase, Boom, Carr, Chevalier, Quick, Verboom and Reimer2019, Reference Chase, Boom, Carr, Quick and Reimer2020), that reflect broad trends in both regions but recognize that the complexity of climatic shifts throughout southern Africa make this synthetic approach highly reductive. Note that the published age models for each record are presented here, with no attempt made to correct or synchronize the various records.
RESULTS
The Morris (Reference Morris1992) catalogue listed 182 dated individuals; our inventory now includes 590 dates for the remains of people who lived during the Holocene in the area now known as South Africa. Of these, 491 dates are associated with the Stone Age, 77 with Iron Age contexts, and for 22 individuals their possible archaeological phase is unknown based on the available dating and distribution information. The earliest 22 dates represent Holocene individuals who lived roughly contemporaneous with the Oakhurst technocomplex, followed by 83 dates roughly contemporaneous with the Wilton (Table 1). Based on their associated contexts, 246 dates fall within the final Later Stone Age. All these individuals represent people who followed hunter-gatherer lifeways. The dated remains of 27 people may be associated with the ceramic final Later Stone Age, representing pastoralist or herding lifeways (Table 1). For 113 dates associated with the Stone Age falling within the last two millennia, we could not distinguish between the final Later Stone Age and the ceramic final Later Stone Age phases, because of the spatiotemporal overlap between these two technocomplexes and inadequate archaeological information—these are indicated as “uncertain final Later Stone Age” in the database.
The Iron Age sample represents African farmers/agriculturalists whose ancestors originated around Nigeria in western Africa and migrated eastwards and southwards over the course of several millennia (Fortes-Lima et al. Reference Fortes-Lima, Burgarella, Hammarén, Eriksson, Vicente, Jolly, Semo, Gunnink, Pacchiarotti and Mundeke2023). The distinction between Early and Late Iron Age archaeological contexts is believed to represent the arrival of new communities in the region, antecedents to Nguni- and Tswana-speaking communities today (Huffman Reference Huffman2007). Dated human remains broadly contemporaneous with the Early Iron Age number 13, with the Middle Iron Age 12, and the Late Iron Age 52. Note that the number of dates grossly underestimates the number of individual sets of remains from some Iron Age sites—for example, more than 120 individuals were recovered from Mapungubwe and K2, two adjacent sites dating to around the turn of the 1st millennium CE. However, destructive sampling of human skeletal material from these sites faces considerable ethical scrutiny over the concerns of descendent communities, and the number of direct dates from these sites is unlikely to increase as the remains have been reburied.
In Figure 1 we show the broad distribution of dated Holocene human remains across South Africa, separated into individuals associated with either their Stone Age or Iron Age contexts. This analysis reveals the large portions of the country from which no recorded, radiocarbon dated human remains have been recovered. The result is probably an artefact of human remains collection practices, archaeological research focus, and preservation conditions on the one hand, and a broad reflection of population distribution and densities on the other (for the distribution of Stone Age vs Iron Age stone walling and occupation see Huffman Reference Huffman2007; Sadr Reference Sadr2012; Lombard et al. Reference Lombard, Caruana, van der Walt and Högberg2020). Stone Age remains cluster around the coastal regions in the Western Cape and Eastern Cape, with some sites also located in the Northern Cape, KwaZulu-Natal, and Free State (Figure 1). One of the northern-most dated Stone Age individual is Bronkhorstspruit TM92-136 (SOM 1), of whom the mummified remains were found in a cave, and her Khoe-San (Stone Age) population affinity was confirmed through mtDNA (Pereira Reference Pereira2007). Most of the dated Iron Age remains were retrieved from the Limpopo Province, with some also from the North West, Free State, KwaZulu-Natal and Mpumalanga. This spatial distribution may reflect the broad southward migration of incoming agriculturalists in the eastern regions of South Africa (Figure 1). Currently, there are no dated Iron Age human remains from any of the Cape provinces, which reflects limited settlement of these regions by agriculturalists due to climatic factors (Russell Reference Russell2020). Table 2 shows the approximate number of dated human remains recovered from each rainfall region. The greater number of dated remains from the winter- and year-round rainfall zones (n = 466 combined) relative to the summer rainfall zone (n = 125) may reflect bone and collagen preservation conditions in these climatic zones and likely also partly reflects patterns of institutional research intensity across the country, and dating practices within hunter-gatherer and agriculturalist archaeology. However, the density of hunter-gatherer burials along the clement southern Cape coast surely speaks to the habitability of this region in the past.
The KDE plot (Figure 2) for the Later Stone Age remains (n = 491) shows at least three major peaks, at ∼5.5 ka cal BP, ∼2 ka cal BP and ∼0.5 ka cal BP. Smaller peaks are apparent in the earlier part of the Holocene, but the small sample size associated with this phase compromises the interpretative power of the fluctuations. The mid-Holocene peak at ∼5.5 ka cal BP seems to be broadly associated with the early phase of the Wilton technocomplex (Table 1). Some authors explained the high degree of standardisation evident in this technocomplex as a response to difficult environmental conditions, and archaeological evidence for markedly lower populations during the Wilton compared to the preceding Oakhurst or the subsequent post-Wilton Later Stone Age phases (Deacon Reference Deacon1984; Wadley Reference Wadley1989; Sealy Reference Sealy2016b). These interpretations require future testing against new fine-grained environmental data (e.g., Chase et al. Reference Chase, Boom, Carr, Chevalier, Quick, Verboom and Reimer2019, Reference Chase, Boom, Carr, Quick and Reimer2020), foraging-fitness assessments (e.g., Lombard and van Aardt Reference Lombard and van Aardt2023), and genetic data for populations sizes (Schlebusch et al. Reference Schlebusch, Malmström, Günther, Sjödin, Coutinho, Edlund, Munters, Vicente, Steyn, Soodyall and Lombard2017; Reference Schlebusch, Sjödin, Breton, Günther, Naidoo, Hollfelder, Sjöstrand, Xu, Gattepaille, Vicente and Scofield2020).
In their work on southwestern Cape burials, Loftus and Pfeiffer (Reference Loftus and Pfeiffer2023) previously also identified the late Holocene peak at ∼2 ka cal BP, and a subsequent decline in human remains. This may have been caused by societal upheavals within foraging groups as a result of the arrival of herders from eastern Africa (e.g., Sampson Reference Sampson2010; Sadr Reference Sadr2013; Breton et al. Reference Breton, Schlebusch, Lombard, Sjödin, Soodyall and Jakobsson2014). The DNA record also indicates a sharp rise in population density amongst the southern African groups at this time (Schlebusch et al. Reference Schlebusch, Malmström, Günther, Sjödin, Coutinho, Edlund, Munters, Vicente, Steyn, Soodyall and Lombard2017, Reference Schlebusch, Sjödin, Breton, Günther, Naidoo, Hollfelder, Sjöstrand, Xu, Gattepaille, Vicente and Scofield2020). The subsequent decrease in human remains may then reflect the establishment of herding groups associated with the ceramic final Later Stone Age (Lombard et al. Reference Lombard, Bradfield, Caruana, Makhubela, Dusseldorp, Kramers and Wurz2022). By ∼500 years ago there is a late peak that may reflect intensified contact between Stone Age and Iron Age groups in the Northern Cape and Eastern Cape and it also coincides with the first European travellers to visit the Cape (Smith Reference Smith1993; Hall and Smith Reference Hall and Smith2000).
The KDE analysis of Iron Age farmer dates shows two marked peaks separated by a trough (Figure 2). Peaks fall roughly at 1.1 ka cal BP during the later phase of the Early Iron Age (Table 1) and at 0.5 ka cal BP during the Late Iron Age, at the same time as a peak in dated Stone Age human remains. An evaluation of the underlying dates (represented as circles underneath the KDE plot) shows that this structured distribution may reflect—at least in part—the small size of the dataset, although the apparent complete lack of dated remains between ∼950–750 cal BP (1150–1250 cal CE) is certainly notable. The decline in dates in the latter half of the 2nd millennium may be due to a tendency amongst archaeologists not to date material attributed to the later part of the Late Iron Age, when dates may be too recent for reliable calibration, and the remains can be closely dated with relative dating methods, for example by ceramic association.
To assess variation between the number of dated human remains and the number of dates obtained from non-human archaeological material, we extracted the “non-human” dates for both the Holocene Stone Age and Iron Age contexts from the Southern African Radiocarbon Database (SARD: Loftus et al. Reference Loftus, Mitchell and Bronk Ramsey2019). For this analysis we used the SPD method, which is computationally more feasible for assessing large datasets than the KDE method (Crema Reference Crema2022). The more complex shape of these curves (Figure 3), compared to the KDE curves (Figure 2), warrants some caution in the interpretation of the results.
Comparison of the human remains with other archaeological Later Stone Age dates confirms the previously observed peaks in excavated human remains at about 5.5 ka cal BP, 2 ka cal BP and 0.5 ka cal BP (Figure 3). For the ∼5.5 ka cal BP and ∼2 ka cal BP peaks in human remains the corresponding peaks in other dated material are, however, relatively lower, whereas the 0.5 ka cal BP peaks in both human remains and other dated material largely overlap with each other. The Iron Age dated human remains curve appears to broadly replicate several of the prominent peaks in the set of non-human archaeological dates (Figure 3), although there is no matching trough ∼950–750 cal BP in the non-human dataset. This correspondence may indicate that the dated human remains reflect broader patterns in Iron Age research and dating endeavour, and/or it reflects the effect of the calibration curve on the SPDs. In general, the comparatively small size of the Iron Age burial dataset makes it difficult to gauge the validity of certain features, such as the near lack of dates between ∼950–750 cal BP (1150–1250 cal CE).
Reconstructions of climatic and environmental conditions during the mid-Holocene indicate variability across the subcontinent (Figure 4; Zhao et al. Reference Zhao, Dupont, Schefuß, Meadows, Hahn and Wefer2016). In general, however, sea-surface temperatures increase during the Holocene compared to the preceding MIS 2 (Caley et al. Reference Caley, Extier, Collins, Schefuß, Dupont, Malaizé, Rossignol, Souron, McClymont, Jimenez-Espejo and García-Comas2018). Several climate indicators may be relevant to the Stone Age groups living in the Cape region. For example, broadly correlating with the ∼5.5 ka cal BP peak observed in dated human remains is a downward trend in the extent of Atlantic sea-ice (Fischer et al. Reference Fischer, Fundel, Ruth, Twarloh, Wegner, Udisti, Becagli, Castellano, Morganti and Severi2007), relatively high levels of humidity in the Eastern Cape (Baviaanskloof record: Chase et al. Reference Chase, Boom, Carr, Quick and Reimer2020), a distinct drop and sharp rise in both the eastern and western regions of the Western Cape (Cango Cave and Pakhuis Pass records: Talma and Vogel Reference Talma and Vogel1992; Chase et al. Reference Chase, Boom, Carr, Chevalier, Quick, Verboom and Reimer2019), and an apparent increase in moisture availability in the Wonderkrater record from the northern Limpopo province (Scott Reference Scott1982; Scott et al. Reference Scott, Holmgren, Talma, Woodborne and Vogel2003) (Figure 4). Broadly correlating with the ∼2 ka cal BP peak is a sharp downward fluctuation in the Atlantic sea-ice record, a distinct rise in the humidity indicators of the Eastern Cape and eastern region of the Western Cape (Baviaanskloof and Cango Cave records), whilst the western region of the Western Cape experienced a drop in humidity (Pakhuis Pass record) (Figure 4). A compilation of multiple pollen records from the summer rainfall zone that covers most of the Northern Cape, and the other non-Cape provinces of South Africa indicates generally moist conditions at ∼2 ka cal BP, declining thereafter (Figure 4; Chevalier and Chase Reference Chevalier and Chase2015). The peak in dated Stone Age human remains at ∼0.5 ka cal BP correlates with upward trends in humidity in several records (Figure 4).
The peaks in dated human remains from Iron Age contexts also show some correlation with the summer rainfall zone precipitation and record and the Wonderkrater moisture index that are relevant to their distribution in northeastern South Africa (Figures 2 and 4). For example, the ∼1.1 ka cal BP peak correlates broadly with a marked downward trend in precipitation in the summer rainfall. The peak at 0.5 ka cal BP—also evident in the Stone Age dated human remains—seem to correspond with peaks in both records. Additionally, shorter, and higher resolution, paleoclimate records (not shown here) attest to regional effects of global climate events over this period, including the Medieval Warming Anomaly (∼1000 to 1300 CE) and the Little Ice Age (∼1500 to 1800 CE) (e.g., Tyson et al. Reference Tyson, Karlen, Holmgren and Heiss2000; Holmgren et al. Reference Holmgren, Tyson, Moberg and Svanered2001).
CONCLUDING DISCUSSION
The main purpose of this contribution is to present an updated database (SOM 1) of human remains excavated from South African Holocene contexts and to align their newly calibrated/recalibrated radiocarbon dates with the published syntheses of the relevant archaeological phases (Table 1). We applied a few analytical methods to detect possible patterns in the spatiotemporal distribution of the dated human remains and possible correlations with some climatic indicators (Figure 4). Subsequent to the Morris (Reference Morris1992) catalogue, this is the first attempt at a large-scale, multi-dimensional representation of dated human remains associated with the South African Holocene. The database, and the coarse-grained trends and correlations highlight aspects that require further detailed and contextual analyses to test their robusticity against the available archaeological records, population histories and climatic indicators through time and across space. Such studies have the potential to deepen our understanding about aspects of human behaviour, interaction, and demography that may not come to the fore through other lines of archaeological or anthropological inquiry.
For example, focusing on the skeletons of late Holocene foragers from the southwestern Cape coast with radiocarbon dates of 3500–900 BP, Loftus and Pfeiffer (Reference Loftus and Pfeiffer2023) hypothesized that socio-economic change coincided with the arrival of pastoralists in the region, supported by a sheep bone age of 2310–1890 cal BP at Spoegrivier (Vogel et al. Reference Vogel, Plug and Webley1997; Coutu et al. Reference Coutu, Taurozzi, Mackie, Jensen, Collins and Sealy2021). Variation in forager behaviour included a shift from marine to more mixed diets, and perimortem trauma from interpersonal violence (Gibbon and Davies Reference Gibbon and Davies2020; but also see Pfeiffer Reference Pfeiffer2016). Some evidence for persistent burial places or “cemeteries” occurs at, or soon after, the introduction of pastoralism to the region (Pfeiffer et al. Reference Pfeiffer, Sealy, Harrington, Loftus and Maggs2020). New burial clusters from this period may reflect attempts to control territory through direct ties with ancestors (Loftus and Pfeiffer Reference Loftus and Pfeiffer2023). This study shows how the arrival of groups with different lifeways and worldviews may disrupt and/or change the lifeways and worldviews of those they encounter on the landscape.
Burial behaviors associated with some of the remains may therefore reflect changes in social practices or societal stresses. However, they are unlikely to be a straightforward reflection of population size, given the complexity of some mortuary behaviors, where intentional internment is only one possible practice. The aggregation of radiocarbon dates to assess changes in archaeological density or past populations through time is a widely deployed approach, which assumes that variations in the frequency of dates through time reflect past human activity (Rick et al. Reference Rick1987; Williams 2012). The approach is, however, subject to several limitations in its interpretative potential (Becerra-Valdivia et al. Reference Becerra-Valdivia, Leal-Cervantes, Wood and Higham2020; Carleton and Groucutt Reference Carleton and Groucutt2021; Crema Reference Crema2022). For example, biases in anthropological and archaeological spatiotemporal research focus, shortcomings in sampling strategies, and varying preservation conditions may all be inherent in the data. What is more, not all human remains have been dated or linked with archaeological phases. For example, the 127 burials from Mapungubwe (n = 28) and K2 (n = 99) (Steyn Reference Steyn1994) are represented by only three directly dated skeletons, and two dates each (maximum and minimum date) for the sites in our inventory. At Matjes River Rock Shelter on the southern Cape coast, poor curation of the human remains has led to extensive commingling of the assemblage. Estimates of the number of individuals range between 86 and 103 (L’Abbé et al. Reference L’abbé, Loots and Keough2008), many more than the 33 unique dates recorded in this dataset, which altogether span ten millennia. Without dates, the full collection of burials from these sites cannot be adequately represented in either behavioural, socio-economical or spaciotemporal models. That said, our inventory of human remains—with their newly calibrated/recalibrated radiocarbon dates and broadly associated archaeological phases—provides a revised, data-driven framework. The next phase of our project will be to use this framework to reintegrate the undated Holocene human remains and explore aspects of demography and behaviour based on the contexts from which they were excavated.
Similar to the syntheses of the archaeological sequences (Table 1), our inventory of dated human remains for the South African Holocene (SOM 1) is not designed to provide detailed interpretations on the observations highlighted here. Instead, it is intended as a heuristic tool for communicating and learning about complex data and broad spatiotemporal trends, highlighting areas that may be of future research interest. Increasingly, the necessity of research that requires the destructive sampling or institutional storage of human remains is being reassessed in the light of widespread and coalescing ethical concerns (see Schroeder Reference Schroeder2020). Yet, many decades of research and curatorial efforts and resources are represented in a dataset such as that presented here and thorough data collations are imperative to maximise the cultural and scientific value of preserved human remains.
ACKNOWLEDGMENTS
We extend a sincere “thank you” to Alan Morris for generously sharing raw data from his original catalogue and some updates that he collected subsequently. Maryna Steyn acknowledges STIAS (Stellenbosch Institute for Advanced Study) for providing the opportunity to work on this project. An Ethics Waiver (Ref: WvN-230913-01; HREC, University of the Witwatersrand) was obtained for this research. Further ethics clearance was not needed as data were taken from published sources and no destructive sampling or new measurements were undertaken specifically for this research.
SUPPLEMENTARY MATERIAL
To view supplementary material for this article, please visit https://doi.org/10.1017/RDC.2024.22