INTRODUCTION
Reliable information regarding past climate and vegetation change in the Namib Desert region is extremely difficult to obtain because the arid environment is not conducive to preservation of organic material and associated proxies (see Lancaster, Reference Lancaster2002; Chase and Meadows, Reference Chase and Meadows2007). Indirect evidence from marine sediments has thus been drawn upon extensively to compensate for this lack of terrestrial archives (e.g., Shi et al., Reference Shi, Schneider, Beug and Dupont2001; Stuut et al., Reference Stuut, Prins, Schneider, Weltje, Jansen and Postma2002; Collins et al., Reference Collins, Schefuß, Govin, Mulitza and Tiedemann2014). However, the influences on such marine records of the region's complex circulation and taphonomic dynamics are difficult to assess, and clear discrepancies exist with those terrestrial records that do exist (Scott et al., Reference Scott, Marais and Brook2004, Reference Scott, Gil Romera, Marais and Brook2018; Lim et al., Reference Lim, Chase, Chevalier and Reimer2016; Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019). As a result, a clear understanding of long-term vegetation change in the Namib Desert region remains to be established.
As an iconic feature of southern Africa's diverse geography, the Namib Desert, which extends over 2000 km along the southwestern coast of Africa, has long been a focus of paleoenvironmental study. Arid conditions in the region likely date to at least the late Cretaceous (ca. 80 Ma) (Ward et al., Reference Ward, Seely and Lancaster1983), while intensification of the South Atlantic Anticyclone and the Benguela upwelling system since the Eocene (ca. 34 Ma) is likely to have induced increasingly arid conditions through higher atmospheric pressure and cooler coastal waters (van Zinderen Bakker, Reference van Zinderen Bakker1975; Tankard and Rodgers, Reference Tankard and Rodgers1978; Etourneau et al., Reference Etourneau, Martinez, Blanz and Schneider2009). Within this context and based on these mechanisms, it has been suggested that the last glacial period, particularly Marine Oxygen Isotope Stages (MIS) 2–4 (12–70 ka), which experienced marked increases in anticyclonic intensity and coastal upwelling (Little et al., Reference Little, Schneider, Kroon, Price, Bickert and Wefer1997; Stuut et al., Reference Stuut, Prins, Schneider, Weltje, Jansen and Postma2002; Farmer et al., Reference Farmer, deMenocal and Marchitto2005; Pichevin et al., Reference Pichevin, Cremer, Giraudeau and Bertrand2005), would have been a period of increased aridity in the Namib Desert (van Zinderen Bakker, Reference van Zinderen Bakker1975). However, studies of relict lacustrine deposits from across the Namib, including sites such as Kannikwa (Beaumont, Reference Beaumont1986), Koichab Pan (Lancaster, Reference Lancaster1984), Narabeb (Teller and Lancaster, Reference Teller and Lancaster1986), Khommabes (Teller and Lancaster, Reference Teller and Lancaster1985), and Gobabeb (Vogel and Visser, Reference Vogel and Visser1981) indicate periods of more/more-effective precipitation (Δprecipitation:potential evapotranspiration) during much of MIS 2–3 (37–31 cal ka BP, 30–24 cal ka BP), as well as the during the late glacial and earliest Holocene (14–10 cal ka BP) (Lancaster, Reference Lancaster2002). Much of these early data are derived from paleolacustrine and paleofluvial (Marker and Müller, Reference Marker and Müller1978; Vogel, Reference Vogel1982; Ward, Reference Ward1982; Heine and Heine, Reference Heine and Heine2002) archives associated with hydrologic networks sourced in the relatively mesic uplands of the Great Escarpment, to the east of the Namib Desert. Thus, an often-posed question has been to what extent does evidence from these archives reflect changes in the uplands or in the desert itself (Ward et al., Reference Ward, Seely and Lancaster1983; Heine, Reference Heine1998; Lancaster, Reference Lancaster2002)?
Questions regarding (1) the significance of phases of increased humidity inferred from discrete lacustrine deposits, and (2) the loci of the rainfall that fostered development of the region's lacustrine and fluvial features are being addressed via a wider network of new paleoenvironmental archives. In recent years, it has been shown that rock hyrax middens (Chase et al., Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012) can provide such valuable archives. Several continuous high-resolution records from rock hyrax middens from the Namib Desert have indicated that the last glacial period was indeed relatively humid compared to the Holocene (Lim et al., Reference Lim, Chase, Chevalier and Reimer2016; Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), adding significant detail to our understanding of regional climate change over the last 50,000 years. While some middens, such as those from Pella (Lim et al., Reference Lim, Chase, Chevalier and Reimer2016), are located at the margins of the Namib Desert, and may be considered to reflect conditions in the Great Escarpment, others, such those recovered from Spitzkoppe (Chase et al., Reference Chase, Meadows, Scott, Thomas, Marais, Sealy and Reimer2009, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019) and, particularly, Mirabib (Scott et al., Reference Scott, Gil Romera, Marais and Brook2018), are located to the west and lie within the arid core of the desert. In contrast to the region's fluvio-lacustrine archives, it is generally considered that rock hyrax middens preserve a predominantly local signal that reflects the foraging range of the animals (less than ~60 m from the midden/shelter site) (Sale, Reference Sale1965; Chase et al., Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012) and that, accordingly, the indications of increased glacial-age humidity at Spitzkoppe and Mirabib can be interpreted to reflect conditions within the desert itself.
Strong similarities are evident in the stable nitrogen isotope records (interpreted as a reflection of water availability) from middens recovered from Aba Huab/Austerlitz (Chase et al., Reference Chase, Meadows, Carr and Reimer2010), Spitzkoppe (Chase et al., Reference Chase, Meadows, Scott, Thomas, Marais, Sealy and Reimer2009, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), Zizou (Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), and Pella (Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), indicating a generally common climate response signal across a 900-km transect of the Namib Desert (Fig. 1). Evidence of associated changes in vegetation, however, is often restricted to the Holocene (Scott et al., Reference Scott, Cooremans, de Wet and Vogel1991, Reference Scott, Gil-Romera, Marais and Brook2022; Scott, Reference Scott1996; Gil-Romera et al., Reference Gil-Romera, Scott, Marais and Brook2006, Reference Gil-Romera, Scott, Marais and Brook2007), with the only indications of glacial-age vegetation coming from the continuous 50-kyr record from Pella (Lim et al., Reference Lim, Chase, Chevalier and Reimer2016), and shorter MIS 2–3 time slices from the Brandberg (Scott et al., Reference Scott, Marais and Brook2004) and Mirabib (Scott et al., Reference Scott, Gil Romera, Marais and Brook2018). While each record expresses changes specific to its environmental context, evidence from all of these sites contrasts markedly with pollen records recovered from marine cores off the coast of Namibia (Shi et al., Reference Shi, Dupont, Beug and Schneider2000, Reference Shi, Schneider, Beug and Dupont2001; Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015). A distinctive component in these marine records—both in terms of their content and the associated proxy interpretations—is the prevalence of pollen types associated with the Cape Floristic Region during the Pleistocene. This is evinced in levels of Restionaceae pollen as high as 14–20% during MIS 2–3 (Shi et al., Reference Shi, Dupont, Beug and Schneider2000, Reference Shi, Schneider, Beug and Dupont2001; Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015). In some cases, the addition of other pollen types, such as Anthospermum, Stoebe-type, and Artemisia pollen, that are found in, but are not specific to, fynbos augments the percentage of pollen from “fynbos taxa” to as much as 40% of the pollen sum (e.g., Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015). Such conclusions have played a significant role in theories regarding expansions of the southwest African winter rainfall zone and Fynbos Biome to the north during glacial periods (Chase and Meadows, Reference Chase and Meadows2007; Shi et al., Reference Shi, Dupont, Beug and Schneider2000, Reference Shi, Schneider, Beug and Dupont2001; Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015).
As early as 2004 the reliability of interpretations of marine pollen records in the region was questioned, based on an absence of Restionaceae and Ericaceae pollen in glacial-age samples from the Brandberg hyrax middens (Scott et al., Reference Scott, Marais and Brook2004). As a single site, this could not provide definitive refutation (see Chase and Meadows, Reference Chase and Meadows2007), but the absences of these taxa in subsequent records (Scott et al., Reference Scott, Gil Romera, Marais and Brook2018), including sites on the southern boundary of the Namib (Lim et al., Reference Lim, Chase, Chevalier and Reimer2016), implies that the microfossils from marine cores in the Benguela system may be significantly influenced by long-distance transport. As such, they cannot therefore be considered a direct reflection of environmental change on the adjacent continent (Scott et al., Reference Scott, Marais and Brook2004, Reference Scott, Gil Romera, Marais and Brook2018; Lim et al., Reference Lim, Chase, Chevalier and Reimer2016; Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019).
To characterise the nature of environmental and vegetation change in the Namib Desert, and to further test the hypothesis of substantial expansions of the Fynbos Biome during the last glacial period, we present new pollen, microcharcoal, and stable isotope records from the Zizou rock hyrax midden site. Located on the plains below the Great Escarpment on the eastern margin of the Namib Sand Sea (Figs. 1, 2), the Zizou rock hyrax midden site encompasses much of the last 39,000 years and allows an analysis of changes in vegetation and fire regimes within the context of an independent perspective (via stable isotope data) on climate variability.
Regional setting
The Zizou rock hyrax midden site is located on the eastern margin of the northern Namib Sand Sea, 20 km SSW of the town of Solitaire (24.07°S, 15.97°E, 1120 m asl) (Figs. 1, 2). Situated at the base of the western flank of the Great Escarpment, the site lies within the steep environmental gradient between the hyper-arid regions of the Namib Desert and the semi-arid interior (Fig. 2). The limited precipitation received at the site (130 mm/yr on average, but with a high degree of interannual variability) falls almost exclusively during the austral summer (Fick and Hijmans, Reference Fick and Hijmans2017). Vegetation at the site is classified as xeric Namibian Savanna Woodland (Fig. 2), characterised by sparse C4 grass cover on the plains surrounding the site, with sporadic shrubs of the Capparaceae family (Boscia, Maerua) and isolated small trees (Vachellia and Searsia) along drainage lines (Fig. 1). Shrubs and arboreal taxa are more common on the rocky slopes surrounding the midden site, including the occurrence of Commiphora. To the west, the Namib Desert ecoregion is characterised by sparse cover of xeric grasses, particularly Stipagrostis, if vegetation is present at all. Here, the water table is generally too low to support even deep-rooting arboreal taxa, and only plants capable of exploiting the region's extremely sporadic rainfall occur. To the east, the uplands of the Great Escarpment and interior plateau support the low-shrub vegetation of the Gariep Karoo ecoregion. Shrubs of the Asteraceae and Aizoaceae families are most common, including Eriocephalus, Pentzia, and Ruschia, as are arid-adapted C4 grasses such as Aristida and Stipagrostis.
Material and methods
The middens considered in this study were collected and analyzed according to the techniques described by Chase et al. (Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019). Depending on the architecture of the shelter in which a midden is found, the accumulation may have varying proportions of faecal pellets to urine (known as hyraceum). At Zizou, two midden sections (ZIZ-1-1 and ZIZ-1-3) were found in close proximity (~2 m) as part of the same midden complex. They were selected for analysis because they are composed primarily of hyraceum, and thus have superior stratigraphic integrity compared to more pellet-rich (and thus coarser) middens (Chase et al., Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012). Representative portions of the middens (measuring 227 mm and 145 mm in depth for ZIZ-1-1 and ZIZ-1-3, respectively) were extracted using an angle grinder and transported back to the laboratory, where the sections were cleaned using progressively finer grades of sandpaper prior to sub-sampling for radiocarbon, microfossil, and stable isotope samples.
Radiocarbon dating
Radiocarbon age determinations for the ZIZ-1-1 (n = 14) and ZIZ-1-3b (n = 13) middens were processed at the 14CHRONO Centre, Queen's University Belfast using accelerator mass spectrometry (AMS) (Table 1; Fig. 3). Samples were pre-treated with 2% HCl for one hour at room temperature to remove carbonates and dried at 60°C. They were then weighed into quartz tubes with an excess of CuO, sealed under vacuum, and combusted to CO2. The CO2 was converted to graphite on an iron catalyst using the zinc reduction method (Slota et al., Reference Slota, Jull, Linick and Toolin1987). The radiocarbon ages were corrected for isotope fractionation using the AMS measured δ13C. The ages were calibrated using SHCal20 calibration data (Hogg et al., Reference Hogg, Heaton, Hua, Palmer, Turney, Southon and Bayliss2020), and rbacon v.2.5.8 software package (Blaauw and Christen, Reference Blaauw and Christen2011, Reference Blaauw and Christen2020) was used to generate the age-depth models for each midden, employing in each case calculated mean accumulation rate of 100 yr/mm (Fig. 3).
Stable carbon and nitrogen isotopes
Following the previous works of Chase et al. (Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), the stable carbon and nitrogen isotope composition of 271 hyraceum samples (~2 mg) was determined at the School of Geography, Geology and the Environment, University of Leicester, with contiguous/overlapping samples obtained from two offset tracks using a 1-mm drill. Isotope ratios were measured on a Sercon 20-20 continuous flow isotope ratio mass spectrometer. For the stable isotope analyses, the standard deviation derived from replicate analyses of homogeneous reference material (Casein protein CatNo.B2155 Batch no. 114859 Elemental Microanalysis Ltd UK) was better than 0.2‰ for both carbon and nitrogen. Carbon isotope results are expressed relative to Vienna PDB, and nitrogen isotope results are expressed relative to atmospheric nitrogen (Fig. 5).
Pollen and microcharcoal
Following the protocol of Lim et al. (Reference Lim, Chase, Chevalier and Reimer2016), pollen samples were prepared with standard physical (600-μm sieving and decanting) and chemical (HCl, KOH, HF, and acetolysis) methods (Moore et al., Reference Moore, Webb and Collinson1991). Lycopodium tablets (Batch 1031, containing ~20,848 spores/tablet) were added to the weighed sample to estimate pollen concentrations (Stockmarr, Reference Stockmarr1971). A minimum pollen sum of 400 grains was counted at a magnification of ×400 under a light microscope, and identified using (1) the literature (van Zinderen Bakker, Reference van Zinderen Bakker1953, Reference van Zinderen Bakker1956; van Zinderen Bakker and Coetzee, Reference van Zinderen Bakker and Coetzee1959; Scott, Reference Scott1982); (2) the photographic and slide reference collections of the University of the Free State, University of Cape Town, and University of Montpellier, and (3) the atlas of southern African pollen types and their climatic affinities (Chevalier et al., Reference Chevalier, Chase, Quick, Scott, Runge, Gosling, Lézine and Scott2021).
Microcharcoal particles were identified as black, completely opaque, angular fragments that occurred in the pollen slides (Clark, Reference Clark1988). Only charcoal particles >75 μm2 (or longer than 10 μm) were counted under a light microscope at ×400 magnification (Patterson et al., Reference Patterson, Edwards and Maguire1987; Mooney and Tinner, Reference Mooney and Tinner2011). A minimum count of 200 items (given by the sum of charcoal particles and exotic marker grains) was used. Charcoal particles that exceed the mesh-width size of 600 μm are missing from the microscopic charcoal record and particles of less than ~10 μm were not counted in order to ensure correct identification (Mooney and Tinner, Reference Mooney and Tinner2011). Accordingly, the charcoal signal is related primarily to the regional fire signal, specifically with local fires (large particles) and remote, extra-regional fires (<10 μm particles) being excluded. The TILIA program was used to construct the pollen diagrams, and pollen zones were determined by CONISS cluster analysis (Grimm, Reference Grimm2011).
In comparing the stable-isotope and pollen data from middens, differences in sample sizes should be considered since the stable isotope record is based on 1-mm overlapping samples, while pollen samples from the Zizou middens, limited by the brittleness and density of the hyraceum and the larger sample sizes required for pollen analysis, vary in width from 3.15–12.54 mm.
RESULTS
Chronology
The radiocarbon and age-depth modeling results suggest continuous midden accumulation despite changes in accumulation rates. Midden structure and sampling constraints have provided a stable isotope record extending from 2500–37,100 cal yr BP and pollen and microcharcoal records extending from 6000–39,000 cal yr BP (average age of microfossil samples, which measure 6 ± 2 mm in thickness). Notable in the ZIZ-1-1 midden, which extends into the Holocene, is the substantial decrease in accumulation rate, which renders the Holocene record from the middens of significantly lower resolution (two samples comprising 8100 years of accumulated hyraceum) (Fig. 3).
Stable isotopes
The δ13C and δ15N records from the ZIZ-1-1 and ZIZ-1-3b middens have each been integrated into single records based on sample age. The δ13C values of the Zizou middens vary from −18.6 to −24.3‰ (Fig. 5), with the highest values occurring during the Holocene, while the last glacial period is generally characterised by lower values. The δ15N values range from 7.2–13.6‰ (Fig. 5) with Holocene and glacial periods being characterised by generally higher (mean 12.1‰) and lower (mean 10.9‰ values) respectively. The δ13C and δ15N records exhibit a degree of covariance at centennial–millennial timescales, with higher δ13C values occurring, broadly, when δ15N values are elevated. The lowest δ13C and δ15N values occur around 30,000 cal yr BP and at the end of the Last Glacial Maximum (LGM)/earliest glacial-interglacial transition, ca. 16,000–18,000 cal yr BP, with relatively higher values during the LGM (especially 20,000–23,000 cal yr BP) and most notably during the Holocene after 7000 cal yr BP.
Pollen and microcharcoal
As with the stable isotope records, the pollen and microcharcoal records from the ZIZ-1-1 and ZIZ-1-3b middens were integrated into single records based on sample age. Fifty pollen taxa were identified, which then were clustered into four statistically significant pollen assemblage zones (ZIZ 1–ZIZ 4) spanning much of the last 39,000 years (Fig. 4). Poaceae, Asteraceae, Capparaceae, and Zygophyllum are the most common pollen types, with much greater pollen diversity in the Pleistocene. The Holocene pollen spectra are dominated by grass (Poaceae) pollen. This Pleistocene–Holocene dichotomy is less evident in the microcharcoal record, but charcoal concentrations are highest in Pleistocene zones 1 and 3.
The lowermost pollen Zone (Zone 1; ca. 39,000 to ca. 35,800 cal yr BP) is characterised by early peaks in Searsia-type and Moraceae pollen among the arboreal taxa but is primarily defined by elevated levels of Capparaceae and Zygophyllum pollen (Fig. 4). Asteraceae pollen is relatively rare (average 10%) and Poaceae pollen is at its lowest levels in the sequence (average 31%).
Pollen Zone 2 (ca. 35,800 to ca. 23,900 cal yr BP) is defined by a strong increase in Asteraceae pollen (average 22%) at the expense of arboreal taxa. Poaceae pollen increases slightly (average 36%) but fluctuates around this mean value. Microcharcoal concentrations decrease slightly in Zone 2, from averages of 913 to 704 fragments/gram.
Subsequently, Pollen Zone 3 (ca. 23,900 to ca. 18,000 cal yr BP) sees a strong increase in microcharcoal, while Asteraceae pollen, particularly Felicia, declines. Zygophyllum pollen also declines in this zone, while Capparaceae again becomes a dominant element in the pollen assemblage (average 15%). Poaceae pollen increases again (average 42%) and stabilises around this mean.
Pollen Zone 4 (ca. 18,000 to ca. 5300 cal yr BP) reflects a fundamentally different vegetation that the Pleistocene assemblages. Poaceae pollen is totally dominant, representing 92% of the pollen recovered, with only very low percentages of drought-resistant taxa such as Aizoaceae, Amaranthaceae-type, and Myrothamnus comprising the other apparent elements. Asteraceae and arboreal pollen are found in trace amounts (<2%).
DISCUSSION
Combining stable carbon and nitrogen data with pollen and microcharcoal data from the same archives provides the opportunity to evaluate vegetation and fire dynamics within the context of an independent climate record derived from midden stable isotope content. As described by Chase et al. (Reference Chase, Scott, Meadows, Gil-Romera, Boom, Carr, Reimer, Truc, Valsecchi and Quick2012, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), rock hyrax δ15N records are interpreted as a proxy for water availability, with higher δ15N indicating drier conditions. As described in several previous studies, and supported by strong correlations with a diverse range of proxy climate records (Chase et al., Reference Chase, Meadows, Scott, Thomas, Marais, Sealy and Reimer2009, Reference Chase, Lim, Chevalier, Boom, Carr, Meadows and Reimer2015a, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), drier conditions induce more N flow to inorganic soil nitrogen pools that are subject to gaseous loss of 15N-depleted products (Austin and Vitousek, Reference Austin and Vitousek1998; Handley et al., Reference Handley, Austin, Stewart, Robinson, Scrimgeour, Raven, Heaton and Schmidt1999; Murphy and Bowman, Reference Murphy and Bowman2009). Plants growing in the resulting 15N-enriched soils exhibit higher δ15N values (Craine et al., Reference Craine, Elmore, Aidar, Bustamante, Dawson, Hobbie and Kahmen2009; Hartman and Danin, Reference Hartman and Danin2010), a signal that is transmitted to the tissues (and excreta) of animals that consume them (Murphy and Bowman, Reference Murphy and Bowman2006; Hartman, Reference Hartman2011; Carr et al., Reference Carr, Chase, Boom and Medina-Sanchez2016).
Hyrax midden δ13C reflects animal diet and, in turn, two aspects of the response of vegetation to climate change: (1) plant photosynthetic pathway of C3 (approximately −34 to −24‰), C4 (approximately −16 to −10‰), and CAM (typically within the range −20 to −10‰) plants (Smith and Epstein, Reference Smith and Epstein1971; Smith, Reference Smith1972; Werger and Ellis, Reference Werger and Ellis1981; Rundel et al., Reference Rundel, Esler and Cowling1999; Boom et al., Reference Boom, Carr, Chase, Grimes and Meadows2014); and (2) changes in C3 plant water-use efficiency (Farquhar and Richards, Reference Farquhar and Richards1984; Ehleringer and Cooper, Reference Ehleringer and Cooper1988; Farquhar et al., Reference Farquhar, Hubick, Condon, Richards, Rundel, Ehleringer and Nagy1989a). Pollen data from Zizou indicate that C3 vegetation in the form of trees and shrubs, was at times abundant during the Pleistocene. Today, C4 grasses are the most dominant vegetation, and these data show this has likely been the case throughout much of the Holocene. Although it is possible that C3 grasses occurred at the site during the Pleistocene, midden δ13C values suggest that C4 grasses were still present. CAM plants also have been abundant at the site (Fig. 5), and as generalist feeders (Hoeck, Reference Hoeck1975; Fourie, Reference Fourie1990) these may have formed part of the hyraxes’ diet, although perhaps not in significant quantities (Fourie, Reference Fourie1983; Carr et al., Reference Carr, Boom and Chase2010, Reference Carr, Chase, Boom and Medina-Sanchez2016).
The Zizou δ15N record (Fig. 5) exhibits a pattern of climate variability that is similar to that observed across the Namib Desert region (Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019). Locally and regionally, these data indicate that the last glacial period was generally more humid than in the Late Holocene. Considering data from a 900 km north-south transect of hyrax middens from the Namib Desert, Chase et al. (Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019) ascribed these changes in humidity to: (1) changes in high latitude boreal and low latitude austral summer insolation, which operate in concert to modify land-sea pressure gradients and the advection of cold air from the Benguela Upwelling System; and 2) mediation of the response to these orbital parameters by Heinrich events and related stadials, which warm the SE Atlantic and result in more humid conditions in the Namib region.
The observed changes in vegetation are largely consistent with changes in climate inferred from the midden δ15N record. Considered within the context of regional ecosystems, wetter conditions may be expected to favor the development of vegetation types similar to either the Angolan Mopane Woodlands or the Gariep Karoo (Fig. 2). Notably, Colophospermum (Mopane) pollen was not recorded in the Zizou middens, while the proportions of pollen from Asteraceae and succulent taxa, the primary plant types of the Gariep and Nama-Karoo ecoregions (sensu Dinerstein et al., Reference Dinerstein, Olson, Joshi, Vynne, Burgess, Wikramanayake and Hahn2017), dominated the glacial-age pollen spectra (pollen zones 1–3) (Fig. 4). Pollen zone 1 (ca. 39,000 to ca. 32,900 cal yr BP), pre-dates the phase of peak humidity from 30,700–28,800 cal yr BP and has an assemblage suggesting a dry karroid vegetation dominated by succulents, with arboreal taxa also being well represented. As conditions became wetter, pollen zone 2 (ca. 32,900 to ca. 22,500 cal yr BP) exhibits a marked decline in succulent pollen, concurrent with peak humidity, and a clear increase in Asteraceae pollen. We interpret these changes as indicating a shift to more mesic Karoo vegetation, comprising grasses and asteraceous shrubs, with a reduced succulent component. Based on the climatic affinities of vegetation types in the region today (Fig. 2), humidity may have been 50–100% higher than the modern situation (e.g., aridity index values increasing from ~0.08 to 0.12–0.16).
A similar sequence of vegetation changes is observed in hyrax middens at Mirabib (Scott et al., Reference Scott, Gil Romera, Marais and Brook2018), on the now grass-dominated plains of the Namib Desert to the north of the Kuiseb River (Fig. 1). As at Zizou, MIS-2–MIS-3 pollen assemblages indicate much higher percentages of Asteraceae, Aizoaceae, and arboreal taxa (particularly Searsia, Euclea, and Capparaceae). These changes in vegetation observed in both the eastern and central Namib Desert, along with the stable isotope data presented here and elsewhere (Chase et al., Reference Chase, Meadows, Scott, Thomas, Marais, Sealy and Reimer2009, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), modify the conclusion of Lancaster's (Reference Lancaster2002) review of Namib Desert paleoclimates, in which it was suggested that there was little evidence for significantly increased rainfall in the Namib Desert during the last glacial period. Our findings are not inconsistent with the view that conditions during this time remained “arid,” since even if humidity at Zizou was doubled, it would still leave the region below the threshold for “semi-arid” conditions (aridity index values of 0.2–0.5; UNEP, 1997). A shift towards drier conditions after ca. 29,000 cal yr BP saw a decline in Asteraceae pollen and an increase in succulent pollen, indicating drier karroid vegetation (Fig. 5).
Pollen zone 3 (ca. 22,500 to ca. 17,700 cal yr BP) spans much of the LGM, which is marked by an abrupt decrease in humidity at Zizou at ca. 24,000 cal yr BP. Vegetation responded to this aridification with a decline in Asteraceae pollen and notable increase in the abundance of Poaceae and arboreal pollen. We interpret this as reflecting a lowering water table, and a change in relative competitive advantage associated with plant rooting depth. As conditions became drier, plant-available water in the near-surface sediments would have been reduced, limiting the occurrence of shallower rooting shrubs. However, arboreal taxa, with their deeper root systems, were able to persist through accessing groundwater resources, while grasses were able to thrive, as they do in the present day, by exploiting superficial soil moisture following occasional rainfall events.
Pollen zone 4 (ca. 17,700 to ca. 5300 cal yr BP) represents the most dramatic shift in the vegetation at Zizou. The regional Namib climate record (Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019) indicates that conditions became more humid after ca. 18,000 cal yr BP, culminating in a period of particularly wet conditions from ca. 15,100–14,400 cal yr BP. This is consistent with a slowing of Atlantic Meridional Overturning Circulation and progressive build-up of heat in the South Atlantic associated with Heinrich stadial 1 (Broecker, Reference Broecker1998; Stocker, Reference Stocker1998; Stocker and Johnsen, Reference Stocker and Johnsen2003; Chase et al., Reference Chase, Boom, Carr, Carré, Chevalier, Meadows, Pedro, Stager and Reimer2015b, Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019). The pollen record from Zizou is of insufficient temporal resolution to address this point but does suggest that after the LGM karroid vegetation and arboreal taxa were largely absent, and opportunistic grasses were the only plants able to persist under the warm, dry conditions that became established. This is inferred to reflect lowering of groundwater resources and dominance of growth forms that were capable of exploiting limited and sporadic rainfall. The better-resolved δ15N values from the Zizou middens vary consistently with other regional records (Stute and Talma, Reference Stute and Talma1998; Chase et al., Reference Chase, Meadows, Scott, Thomas, Marais, Sealy and Reimer2009, Reference Chase, Meadows, Carr and Reimer2010; Carr et al., Reference Carr, Boom and Chase2010) that suggest a relatively mesic Early Holocene. The decreased midden accumulation rate and the pollen assemblages suggest a decrease in plant biomass and a decline in the population of the hyrax colony at the site during the Holocene. In the absence of any significantly better-resolved Holocene pollen sequences from the western slope of the Escarpment to facilitate comparison (Scott, Reference Scott1996; Scott et al., Reference Scott, Gil-Romera, Marais and Brook2022), the samples from Zizou presently can provide only an initial indication of vegetation at the site in the post-glacial period.
Throughout the record, changes in vegetation are reflected in the midden δ13C values, which track variations in Poaceae (predominantly C4) and to a lesser degree succulent (predominantly CAM) pollen percentages (Figs 4, 5), supporting the evidence that woody/arboreal (C3) vegetation was dominant during a more humid last glacial period. Together with the observation that grass dominated the Holocene vegetation, our observations indicate a limited role of atmospheric CO2 concentrations as a driver of glacial–interglacial vegetation dynamics in the Namib Desert. This differs from observations and interpretations in other contexts, where atmospheric CO2 concentration is thought to play a central role in vegetation change (Bond and Midgley, Reference Bond and Midgley2000; Bragg et al., Reference Bragg, Prentice, Harrison, Eglinton, Foster, Rommerskirchen and Rullkötter2013; Prentice et al., Reference Prentice, Cleator, Huang, Harrison and Roulstone2017; Dupont et al., Reference Dupont, Caley and Castañeda2019) through the competitive advantage that C3 woody plants have under augmented CO2 conditions (Farquhar et al., Reference Farquhar, Hubick, Condon, Richards, Rundel, Ehleringer and Nagy1989a, Reference Farquhar, Ehleringer and Hubickb; Bond and Midgley, Reference Bond and Midgley2000; Donohue et al., Reference Donohue, Roderick, McVicar and Farquhar2013). The evidence here suggests that Holocene high-CO2 conditions were accompanied by the expansion of C4 grasses and a virtual absence of woody C3 taxa. This pattern of vegetation change has been observed elsewhere, for example in the southern Namib at Pella (Lim et al., Reference Lim, Chase, Chevalier and Reimer2016), and in other dryland regions such as the North American southwest (Holmgren et al., Reference Holmgren, Betancourt, Peñalba, Delgadillo, Zuravnsky, Hunter, Rylander and Weiss2014). In summary, this evidence suggests that, while variation in CO2 may influence C3/C4 vegetation dynamics, it is not necessarily the dominant driver in some contexts over glacial-interglacial timescales (see Nolan et al., Reference Nolan, Overpeck, Allen, Anderson, Betancourt, Binney and Brewer2018).
The microcharcoal content of the Zizou middens facilitates consideration of the region's fire history and ecology. Generally speaking, arid and semi-arid regions are considered to have fire regimes that are biomass limited (Marlon et al., Reference Marlon, Bartlein, Daniau, Harrison, Maezumi, Power, Tinner and Vanniére2013). As a result, charcoal concentrations may be expected to increase in these environments under wetter conditions, with denser vegetation providing more biomass, allowing fire to travel more effectively and thereby producing more charcoal. At Zizou, the relationships between climate, vegetation, and fire indicate that this nexus may be more complex. During the last glacial period, phases of relatively wet conditions are shown here to be associated with reduced microcharcoal concentrations (Fig. 5), an apparent anomaly that may be explained by the relative ability of different vegetation types to carry fire due to spacing of individual plants. For example, the relatively wide spacing between plants in karroid vegetation limits its potential to carry fire (Mucina et al., Reference Mucina, Rutherford, Palmer, Milton, Scott, Lloyd, van der Merwe, Mucina and Rutherford2006; Hebbelmann et al., Reference Hebbelmann, O'Connor and du Toit2022). At Zizou, increased charcoal concentrations correspond to phases with higher percentages of grass pollen, implying either that grass expansion increased charcoal concentrations as a function of increased biomass, or that it may play an important role in carrying fire in such environments (Rahlao et al., Reference Rahlao, Milton, Esler, Van Wilgen and Barnard2009). Elevated percentages of arboreal pollen during these phases likely indicates the presence of an important fuel source, further augmenting charcoal concentrations in what may have been more extensive fires. In contrast, the terminal Pleistocene shift to sparse xeric grassland was accompanied by a decline in microcharcoal concentrations, conforming to predictions based on a marked decreased in biomass.
The Zizou pollen record provides further confirmation of the concerns raised by Scott et al. (Reference Scott, Marais and Brook2004) regarding the interpretation of pollen and other terrestrial material (e.g., microcharcoal and leaf waxes) recovered from marine cores in the Benguela Upwelling System. For example, whereas abundant Restionaceae pollen is evident in regional marine cores during MIS 2–3, including evidence used to infer significant expansions of the Fynbos Biome and winter rainfall zone (Shi et al., Reference Shi, Dupont, Beug and Schneider2000, Reference Shi, Schneider, Beug and Dupont2001; Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015), no similar evidence has been recovered from the Namib Desert (Scott et al., Reference Scott, Marais and Brook2004, Reference Scott, Gil Romera, Marais and Brook2018; Lim et al., Reference Lim, Chase, Chevalier and Reimer2016; this paper).
Modeling studies of wind trajectories indicate the primary source area for terrestrial components in marine cores GeoB 1710 (Dupont and Wyputta, Reference Dupont and Wyputta2003), GeoB 1712, (Shi et al., Reference Shi, Schneider, Beug and Dupont2001), MD96-2098 (Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015), and MD08-3016 (Collins et al., Reference Collins, Schefuß, Govin, Mulitza and Tiedemann2014) lies directly to the east of the core sites (Dupont and Wyputta, Reference Dupont and Wyputta2003). This observation, combined with finding that aeolian transport trajectories deviated little during the LGM (Dupont and Wyputta, Reference Dupont and Wyputta2003), has been used to support the use of marine records as evidence of environmental change on adjacent continental regions (Dupont and Wyputta, Reference Dupont and Wyputta2003; Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015). However, there are no terrestrial records from this putative source region that support this interpretation (Scott et al., Reference Scott, Marais and Brook2004, Reference Scott, Gil Romera, Marais and Brook2018; Lim et al., Reference Lim, Chase, Chevalier and Reimer2016; this paper). Records of terrestrial material preserved in marine sediments from this region, for example pollen (Urrego et al., Reference Urrego, Sánchez Goñi, Daniau, Lechevrel and Hanquiez2015), leaf waxes (Collins et al., Reference Collins, Schefuß, Govin, Mulitza and Tiedemann2014), or charcoal (Daniau et al., Reference Daniau, Sánchez Goñi, Martinez, Urrego, Bout-Roumazeilles, Desprat and Marlon2013), require further confirmation of their possible sources. Without such information, reliable reconstructions of the adjacent terrestrial environments are questionable. Possible alternative hypotheses to explain the anomaly between marine and terrestrial records, such as changes in source area driven by demonstrated changes in wind strength (Chase et al., Reference Chase, Niedermeyer, Boom, Carr, Chevalier, He, Meadows, Ogle and Reimer2019), warrant further consideration.
CONCLUSIONS
• The stable carbon and nitrogen isotope, pollen, and microcharcoal data from the Zizou rock hyrax middens provide evidence for significant differences between MIS 2–3 and Holocene climates and vegetation types in the central Namib Desert.
• While Pleistocene conditions were likely still arid to semi-arid, they were significantly more humid than in the later Holocene.
• During MIS 2–3, shrubby karroid vegetation developed under the most humid conditions, while deeper-rooting arboreal taxa were relatively more abundant under intermediate conditions, in a response to lowering water table. The modern xeric grass-dominant vegetation, an opportunistic vegetation that is capable of exploiting sporadic rains, was established after the last glacial period.
• The Zizou microcharcoal record is most strongly correlated with percentages of grass and arboreal pollen, rather than with humidity. The karroid vegetation that characterises the most humid periods is not conducive to extensive fires because inter-plant spacing does not carry fire well. Under drier conditions, a better-developed grass component carries fire more effectively and the persistence of a significant arboreal component provides fuel, raising charcoal concentrations. This fuel source apparently was absent during the period reflected in our Holocene samples, and charcoal concentrations declined.
• Evolution of the vegetation composition at Zizou does not indicate that changes in atmospheric CO2 played a dominant role in determining plant functional type. The dominance of grasses (primarily C4) during the Holocene and shrubs and arboreal (C3) taxa during the Pleistocene indicates that climate is the primary driver of vegetation composition.
• As at other terrestrial sites from the Namib Desert, no elements of the Cape flora (e.g., Restionaceae) were found in the Zizou pollen assemblages. These results differ markedly from pollen records obtained from marine sediment cores, requiring a fundamental reconsideration of interpretations of environmental change in the region based on marine records.
Acknowledgments
Funding was received from the European Research Council (ERC) under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC Starting Grant “HYRAX”, grant agreement no. 258657, and from the Leverhulme Trust.