Gueldaman GLD1 Cave is one of a series of karstic caves formed within the southeastward slope of the Adrar Gueldaman ridge, western part of Babor mountains in Northern Algeria (Kherbouche et al., 2014). It is located close to the large Soummam River, 5–6 km from the Akbou town, and approximately 65 km southern inland from the Western-Mediterranean Sea (36∘26′ N, 4∘34′ E, 507 m a.s.l.; Fig. 1). Gueldaman GLD1 is a relatively short cave (total extension of ∼ 80 m) that developed in Jurassic limestone. The entrance, facing to the SE, is a semi-circular ∼ 6 m large arch, leading to a dome-shaped ∼ 10 m high and 6 m wide corridor which ends with the main chamber “Grande Salle” at a depth of 30–40 m. Previous archaeological excavation and field investigation suggest that Gueldaman GLD1 Cave has probably been open since at least ca. 7 ka ago (Kherbouche et al., 2014). The area is covered by a thin layer (< 10 cm) of soil derived from the limestone bedrock, wind-blown silicate dust, and organic matter from local vegetation such as Pistacia lentiscus, Quercus ilex, Buxus sempervirens, typical Mediterranean Garrigue type plant assemblage (C3 dominated).

Local climate is Mediterranean semi-arid type, characterized by hot-dry summers and mild-wetter winters. From the ERA-interim reanalysis data between 1979 and 2013 (http://apps.ecmwf.int/datasets/) the annual total rainfall is 516 mm, and the annual mean temperature is 17.2 ∘C. Rainfall occurs rarely in the summer (37 mm) but relatively evenly through the autumn (155 mm), winter (178 mm) and spring (147 mm). Gueldaman GLD1 Cave is well ventilated with the outside atmosphere due to its larger opening and shorter extension. Hobo logger data at 10 min resolution from November 2013 to April 2015 shows significant variations in cave temperature ranging from 13.7 to 19.5 ∘C. The relative humidity varies from 56 to 94 %. Carbon dioxide has not been measured, but it is likely to be close to the atmospheric value.

Two stalagmites and three modern calcites samples were collected in 2012 and 2013 from the main chamber of Gueldaman GLD1 Cave. Both stalagmites were laid on the cave floor during collection. Stalagmite GLD1-stm2 is 350 mm long and 100–200 mm wide; stalagmite GLD1-stm4 is 203 mm long with a diameter of 50–120 mm (Fig. 2). They were halved and polished along the longitudinal axis. It was noticed from the sectioned sample that stalagmite GLD1-stm2 was broken at the depth of 5 mm from the top and covered by the calcite deposited at a certain time after that point. The top 5 mm was not analyzed in this study. Both stalagmites show well-marked laminae with several shifts in the drip apex of the lower parts (Fig. 2). Black bandings, with visible incorporations of charcoal particles, are found throughout both stalagmite profiles.

U/Th dating – Seventeen powder samples were drilled from the two stalagmites and dated by a multi-collector inductively coupled plasma mass spectrometer. The procedure to separate uranium and thorium was referred to in Edwards et al. (1987) and Cheng et al. (2013). The dating work was carried out at the University of Minnesota (USA) and the Xi'an Jiaotong University (China). One dating from the base part of stalagmite GLD1-stm2, for the exploration of preliminary age frame, was done at the Laboratoire des Sciences du Climat et de l'Environnement (LSCE, France). The U/Th dates were reported in years before 2000 AD. (Fig. 2, Table 1). The age model for both stalagmites was developed using the StalAge program (Scholz and Hoffmann, 2011) where a linear interpolation between depth and age is made through each progressive triplet of adjacent U/Th dates (Fig. 3). This procedure provides a quantitative estimate of age uncertainty continuously along the record despite having analytical constraints only at locations where the U/Th dates exist. Stalagmite growth rates were calculated based on the StalAge age model (Fig. 4).

Stable isotopes – Four hundred and thirty samples were drilled every 1 to 2 mm along the stalagmite central growth axis (Fig. 2). Stable carbon and oxygen isotopes compositions of both stalagmites and modern calcites were measured using a VG-OPTIMA mass spectrometer at the LSCE. For each analysis, 60 to 80 µg calcite powder is reacted with phosphoric acid at 90 ∘C, and the resultant CO2 is measured relative to a reference gas that has been calibrated against a series of isotopic standards. Duplicates were run every 10 to 20 samples to check replicability. All values are reported in ‰ relative to the V-PDB (Fig. 4). The error is 0.08 ‰ for δ18O and 0.05 ‰ for δ13C.

Archaeological excavations were carried out at two sectors S2 and S3 inside the Gueldaman GLD1 Cave during the 2010–2012 campaign (Fig. 1). This work consisted mainly in collecting, identifying, and referencing the archaeological materials found in stratigraphic layers (refer to Kherbouche et al., 2014 for details). More than 7000 anthropogenic remains were collected, consisting mainly of faunal remains, ceramic, and lithic and bone tools. Besides, all sediments were water screened through 1.5 and 4 mm mesh and subjected systematically to flotation with collection in a 250 µm mesh yielding a huge amount of charcoals. Initial radiocarbon dating of upper stratigraphic sequences from S2 and S3 gave the median ages ranging from ca. 6800 to 1500 cal yr BP (Kherbouche et al., 2014). In order to refine the chronology of these deposits, in this study, six new charcoal samples were collected from the key archaeological layers in excavation area MN 47/48 of S2. These samples were dated using the AMS radiocarbon method at the CEA Saclay (France). Detailed procedures of the chemical preparation and the dating in the lab were referred to Cottereau et al. (2007). The dates were calibrated using the IntCal13 data set (Reimer et al., 2013) and reported in years before 2000 AD (Table 2).

The uranium contents of measured stalagmites samples are relatively high ranging from 95 to 225 ppb (Table 1). The 2 sigma U/Th errors vary from 20 to 210 years with an average of 77 years (1.6 %). The U/Th date (5070 ± 194 yr BP) of sample GLD1-stm4-47 was detected as a major outlier by the StalAge program, thus, it was not used to calculate the final age model. Calculated StalAge age model for stalagmite GLD1-stm4 shows large errors up to 500 years during ca. 4900–4200 yr BP (Fig. 3). This may partly be due to relatively large errors of adjacent U-series dates and/or impropriate hypotheses applied in the algorithm. Based on individual StalAge age model, stalagmite GLD1-stm2 grew from ca. 6200 to 4000 yr BP (4100 yr BP if excluding the top 5 mm), whereas stalagmite GLD1-stm4 grew from ca. 5800 to 3200 yr BP (Fig. 3).

Stalagmite GLD1-stm2 shows high and variable growth rates (mean = 180 µm yr-1) with higher values ∼ 400 µm yr-1 at ca. 4800–4600 yr BP; whereas stalagmite GLD1-stm4 shows relatively lower and less variable growth rates (mean = 120 µm yr-1) with higher values ∼ 200 µm yr-1 at ca. 3800–3200 yr BP (Fig. 4).

The isotopic compositions of modern calcite vary from -5.40 to -5.56 ‰ for the δ18O and from -8.43 to -10.34 ‰ for the δ13C. The δ18O values from stalagmites GLD1-stm2 and GLD1-stm4 range from -7.8 to -2.8 ‰ and from -7.3 to -0.6 ‰, respectively; the δ13C values range from -10.6 to -3.3 ‰ and from -11.9 to -0.6 ‰, respectively. The δ18O and δ13C significantly correlate in both stalagmites: R= 0.87, P < 0.01 for GLD1-stm2 and R= 0.92, P< 0.01 for GLD1-stm4. Albeit the different amplitudes, the isotopic profiles of the two stalagmites show similarities during their common development of ca. 5800–4000 yr BP: relatively elevated isotope values are found at ca. 5700–5400, ca. 5200, and ca. 4500 yr BP (Fig. 4). Two other isotopically enriched periods in stalagmite GLD-stm2 are found at ca. 6200 and ca. 4900 yr BP (Fig. 4). There is a common isotopic enrichment trend since ca. 4800–4600 yr BP (depending on individual age model; abrupt in stalagmite GLD1-stm4 whereas more gradual in GLD1-stm2). Toward the end of this trend, the most prominent anomaly occurs in stalagmite GLD1-stm4 at ca. 4400–3800 yr BP during which the δ18O values are enriched by approximately 3.5 ‰ relative to the background values of that time as well as the modern calcite values for a period of ∼ 500 years (Fig. 4). Specifically within this anomalous period, there is a mild isotopic depletion at ca. 4200–4000 yr BP, blanketed by two significant enrichments at ca. 4400–4200 yr BP and ca. 4000–3800 yr BP (Fig. 4). The last part, ca. 3800–3200 yr BP, of GLD1-stm4, is characterized by a δ18O recovery of about -3 ‰, synchronous with increased growth rates (Fig. 4).

Excavations inside Gueldaman GLD1 Cave revealed a large variety of archeological remains and, among them, are numerous precious macro charcoals that have been used for establishing the chronology of the deposits. In the ∼ 7 m2 total excavated area of S2, more than 7000 archaeological objects were identified and consisted of faunal remains, lithic artifacts and grinding equipment, potteries, bone tools, ornaments, and ochre. In addition, a fragment of a human mandible and two isolated teeth were found during the excavation 2010–2012 in Gueldaman GLD1 Cave (Kherbouche et al., 2014). These deposits belong mainly to the Neolithic; only the top level of the sequence contains potsherds of the historic period. In the lower Neolithic levels, identified domestic species (i.e. sheep and goats) represented ∼ 25 % of total faunal assemblages (N = 2378) suggesting a partly pastoral based economy. The potteries (N = 825) are mostly related to cooking vessels of 25-40 cm rim diameter. Hundreds of black charcoals (> 1 cm) were found and always associated with ceramic concentrations suggesting evidence of cooking activities.

Determined radiocarbon dates give the median ages of the sequence between 7002 and 1482 cal yr BP, with their 2σ-error intervals varying from 132 to 374 years (Table 2). These dates provide a first chronology for the archaeological deposits excavated from sector S2 (Fig. 6; Kherbouche et al., 2014): anthropogenic remains (i.e. charcoals, bones, teeth and potteries) are numerous during ca. 7002–6003 cal yr BP (depths of ∼ 150-120 cm), decreased at ca. 6003–4918 cal yr BP (depths of ∼ 120–105 cm), most abundant in the period of ca. 4918–4403 cal yr BP (depths of 105–75 cm), significantly diminished during the long interval of ca. 4403–1484 cal yr BP (depths of 75–60 cm), and finally, numerous again from ca. 1484 cal yr BP (depths of ∼ 60–50 cm). With an overall decrease in archeological materials, there are two levels clearly marked by their poverty in charcoal and pottery during the periods of ca. 6003–4918 and ca. 4403–1484 cal yr BP (Fig. 6).

Under isotopic equilibrium precipitation, stalagmite calcite δ18O depends mainly on the temperature of calcite-water fractionation and on the δ18O of drip water that is controlled by local rainfall δ18O (Genty et al., 2014; Lachniet, 2009). Observations from the IAEA network show that the rainfall δ18O at many Mediterranean stations (including one in Algiers, Algeria) are partly controlled by the amount of rainfall (IAEA, 2005), which is coherent with previous studies that stalagmite δ18O records from the Mediterranean regions were interpreted to primarily reflect changes in rainfall amount (inversely correlated; e.g., Bar-Matthews and Ayalon, 2011; Bar-Matthews et al., 2003, 1997; Drysdalea et al., 2004, 2006; Zanchetta et al., 2014). Therefore, higher stalagmite δ18O values are expected during periods of decrease in rainfall (drier). The temperature effect on calcite-water fractionation, on the other hand, is partly counteracted by the condensation temperature effect on rainfall δ18O (Drysdale et al., 2006). We notice that this interpretation may particularly hold true for the present study because the regional temperature seems to have been relatively constant since the Mid-Holocene (Martrat et al., 2004).

Stalagmite δ13C variations have several potential causes, the most likely, considering the studied location and time interval, being variations in soil CO2 input and water flow rate (Genty et al., 2001a; McDermott, 2004). Despite the fact that soil biogenic CO2 production varies according to both temperature and moisture level, moisture is likely to be a major controlling factor due to low temperature variability of the considered time interval and limited water availability under semi-arid climates. Moisture also influences water flow rate and thus the CO2 loss during the prior calcite precipitation (Fairchild et al., 2000). Longer residence time due to lower flow rate, under diminished moisture condition, enhances CO2 loss and preferential 12C removal from solution causing enrichments in stalagmite δ13C (Johnson et al., 2006; Mickler et al., 2004). Eventually, atmospheric rainfall largely determines the moisture level and controls the δ13C variations. Therefore, the significant correlation of stalagmite δ13C and δ18O may suggest a common control of rainfall.

The rainfall signal imprinted in the Gueldaman stalagmite isotopes (inversely correlated) is probably amplified by two other processes – evaporation and disequilibrium isotopic fractionation (Mickler et al., 2006) that are very likely to have occurred at the Gueldaman GLD1 Cave due to its large entrance. Longer residence time during drier (lower rainfall) periods would allow extended evaporative and non-equilibrium fractionations, which drives stalagmite isotopes further higher (Mickler et al., 2004). It has recently been observed that evaporation in semi-arid caves could cause 4–5 ‰ δ18O enrichments of a wide range of drip waters (Cuthbert et al., 2014). The Hendy test (i.e. studying the isotopic variation in contemporaneous laminae; Hendy, 1971) carried out at three different depths in stalagmite GLD1-stm2 show that the δ18O and δ13C significantly correlate and simultaneously increase by up to ∼ 1 ‰ from the apex to the edge probably due to the presence of progressive kinetic fractionations when feeding waters flow outwards from the apex, suggesting that the stalagmite formed out of isotopic equilibrium (see the Supplement). These two processes may partly account for the significant correlation of the δ18O and δ13C profiles in both stalagmites (R= 0.87 for GLD1-stm2; R= 0.92 for GLD1-stm4). Comparisons of two stalagmites show different amplitudes in their isotopic profiles. This might be due to the following: (1) having been fed by reservoirs that smooth rainfall signal differently, and (2) having been suffered from variable evaporative and kinetic enrichments associated with different recharge features. Despite this discrepancy the isotopic profiles of two stalagmites broadly show similar patterns. Consequently, synchronous δ18O and δ13C variations in the two stalagmites can be interpreted in terms of humidity change. A prolonged severe drought is inferred from the most elevated isotopes values during ca. 4400–3800 yr BP together with drier periods from higher values at ca. 6200, ca. 5700–5400, ca. 5200, and ca. 4500 yr BP (Fig. 4), while wetter periods are suggested from depleted isotopes at ca. 4600–5100 and ca. 5300 yr BP.

Moreover, stalagmite growth phase has long been regarded as an environmental indicator (e.g., Genty et al., 2001b; Stoll et al., 2013) and could be used to complement and/or test the climatic interpretation from isotope records. Water availability is an essential factor for stalagmite growth in arid-to-semi-arid areas, as shown by Vaks et al. (2013) who correlated growth periods well with periods of effective rainfall regimes. The growth cessation of stalagmite GLD1-stm2 at ca. 4000 yr BP may suggest a phase of increased aridity, which is consistent with the dryness inferred from elevated isotopes (Fig. 4). However, the continuous growth of stalagmite GLD1-stm4 at that time suggests that the growth cessation might also have been caused by shifts in dripping position under deteriorating climates (Fairchild and Baker, 2012). Moreover, fast stalagmite growths together with wide diameters usually are associated with high drip rates suggesting humid conditions. A wet period ca. 4800–4600 yr BP is indicated by the high growth rates of stalagmite GLD1-stm2, which is broadly in line with the humid period inferred from depleted isotopes (Fig. 4). Another wetter period ca. 3800–3200 yr BP is suggested by faster growths and relatively depleted isotopes of stalagmite GLD1-stm4 (Fig. 4). The discrepancy in the growth rate profiles of the two stalagmites is possibly due in part to (1) the lack of substantial U/Th dating especially for stalagmite GLD1-stm4 between 5023 and 4197 yr BP, and/or (2) site-specific processes due to different reservoirs play a key role in controlling stalagmite growth.

Finally, the modern calcite isotopic values fall in the range of two stalagmite records (Fig. 4), which suggests that the current humidity condition in Northern Algeria seems to be within the range of its Mid-Holocene variability. The average δ18O value during ca. 4400–3800 yr BP in stalagmite GLD1-stm4 is more enriched by about 3.5 ‰ relative to the modern values (Fig. 4), therefore, the 4400–3800 yr BP climate anomaly may be considered analogous to end numbers of the most recent and ongoing drying.

High-resolution absolute-dated Mid-Holocene climate records are rare in the Western-Mediterranean basin, however, there are number of paleoenvironment studies using sediment cores that document large oscillations in vegetation ecology and provide clues of climate change during the Mid-Holocene. A drying trend from ca. 4600 cal yr BP onwards was inferred, based on the decreasing pollen ratio of deciduous broad-leaf vs. evergreen sclerophyllous taxa at Capestang in the Mediterranean, Southern France (Jalut et al., 2000). At a nearby site in Northeastern Spain, the most arid Mid-Holocene condition at ca. 4800–4000 cal yr BP was interpreted using maximum salinity values, more positive organic carbon isotope values, and decreased algal productivity in Estanya Lake (Morrellon et al., 2009). In the Mediterranean, Southern Spain, desertification phases at ca. 5200 and ca. 4100 cal yr BP were inferred using multiple paleoecological indicators including pollen, microcharcoal, spores of terrestrial plants, fungi, non-siliceous algae, and other microfossils in Siles Lake (Carrión, 2002). Similar environment changes have also been observed in the Central Mediterranean basin, as shown by a synthesis study of lacustrine palynological data which suggested a dryness peaking at ca. 4000 cal yr BP (Sadori et al., 2011). Evaluated carbonate oxygen isotopes from Lake Shkodra were argued to be an indicator of dryness at ca. 4100–4000 cal yr BP (Zanchetta et al., 2012b). Increasing aridity at ca. 5000–4000 cal yr BP was suggested to explain the increases in non-tree pollen percentage and micro charcoal content in the Lago di Pergusa Lake, Sicily (Roberts et al., 2011; Sadori and Giardini, 2007; Sadori and Narcisi, 2001). Closer to our site, a study at Preola Lake, Eastern Sicily, documented a significant low stand lake level at ca. 4500–4000 cal yr BP suggesting extreme aridity (Magny et al., 2011). Moreover, six tephra layers were carefully studied to correlate climate anomalies at ca. 4500–3800 cal yr BP from different archives in the central Mediterranean (Zanchetta et al., 2012a). In the east, Finné et al. (2011) reviewed the climate history of the Eastern Mediterranean over the last 6000 years and concluded with much evidence of drying conditions at ca. 4600–3800 yr BP. Although the sampling and dating resolutions in most of the above studies are low, they are in good agreement with the present study regarding the 5200 yr BP dry event, the drying trend from 4800–4600 yr BP onward, and the 4400–3800 yr BP drought.

Recently, high frequency Mid-Holocene climate change has been well documented by multiple detailed speleothem studies from the Central and Eastern-Mediterranean basin (Fig. 5). Drysdale et al. (2006) demonstrated a severe drought peaked at 4400–3800 yr BP through multiproxy analyses on a U/Th-dated flowstone from Renella Cave, Central Italy. This finding is supported by following work at nearby Corchia Cave. Zanchetta et al. (2007) established a high-resolution speleothem isotope record for the entire Holocene and interpreted it as an indicator of rainfall change, in which a drier period persisting from ca. 4800 to 3800 yr BP could be inferred from relatively elevated δ18O values. Elemental analysis on the same stalagmite by Regattieri et al. (2014) presented more clear signals of dryness peaked at ca. 4700–4200 yr BP. In the Eastern-Mediterranean basin, Bar-Matthews and Ayalon (2011) explicitly discussed the Mid-Holocene climate change by high-resolution dating and isotopic analyses on speleothems from Soreq Cave in Israel; they identified multiple dry events at 6250–6180, 5700–5600 and 5250–5170 yr BP as well as a long drying trend since ca. 4700 yr BP peaked at 4200–4050 yr BP. Zanchetta et al. (2014) carefully compared the speleothem isotope records from Corchia Cave and Soreq Cave in the Central and Eastern-Mediterranean basin, and found two coeval dry events at ca. 5600 and ca. 5200 yr BP from the comparable δ18O enrichments in two speleothems. Detailed comparisons of these speleothem records with our new records from Gueldaman GLD1 Cave reveal many consistencies. In particular, periods with elevated δ18O values at ca. 5700–5400, 5200 and 4400–3800 yr BP observed in Gueldaman stalagmites are all identifiable in the speleothems from Corchia, Renella and Soreq Cave (Fig. 5), suggesting that anomalous dryness at these periods synchronously developed across the Mediterranean basin. These observations indicate that climates in the Mediterranean basin might have been under an identical regional scale atmospheric regime during the Mid-Holocene.

It has been suggested that climate change in mid-latitude Europe and the Mediterranean might arise from a perturbation of the Westerlies from a high-latitude trigger (i.e. the North Atlantic; Bond et al., 2001; Drysdale et al., 2006; Zanchetta et al., 2014) or from dynamics within the tropics (Booth et al., 2005; Hoerling and Kumar, 2003). The three dry periods in the Mediterranean are broadly in phase with the ice rafting events in the subpolar North Atlantic (Bond et al., 2001), which suggests some links with the North Atlantic circulation. Based on the coincidence with the elevated wind strength in Iceland (Jackson et al., 2005), Zanchetta et al. (2014) argued that the dry events at ca. 5700–5400 and ca. 5200 yr BP might be caused by reduction of vapor advection into the Mediterranean, due to the intensification and northward displacement of the North Atlantic Westerlies. However, lacking evidence of strengthened wind in the fourth millennium BP argues for a different forcing of the 4400–3800 yr BP drought. The considerably lower amplitude of the Bond ice rafting event at ca. 4200 yr BP than at the fifth millennium BP also indicates a varied ocean-atmosphere circulation state. The modern mid-latitude droughts (1998–2002) have been linked to the increased warmth in equatorial oceans (Booth et al., 2005). During this event, SST changes lead to persistent high pressure over the Northern Hemisphere's mid-latitudes, causing widespread synchronous drought (Hoerling and Kumar, 2003). However, the challenge in applying the dynamics under the 1988–2002 drought towards an understanding of the 4400–3800 yr BP climate anomaly is that while the mechanism operates effectively on short timescales, it has never been tested as to whether they could produce an anomalous climate mode for several centuries (Berkelhammer et al., 2013). General circulation model simulations that begin with realistic boundary conditions and are perturbed with a variety of forcings have been successfully undertaken to understand potential mechanisms that lead to the 8200 yr BP event (Tindall and Valdes, 2011). Similar efforts would be a useful starting point to produce hypotheses for the dynamical underpinnings of the 4200 yr BP event.

A regional drought ca. 4200 yr BP has been widely linked to ancient cultural changes in the Eastern Mediterranean and Asia (Staubwasser and Weiss, 2006), though, in many cases, climatic inferences have been derived from sites that are distant to these human settlements. For instance, evidence of reduced precipitation from elevated δ18O of Soreq stalagmites (Israel) and increased dust input into the Gulf of Oman sediment core has been suggested to contribute to the collapse of the Akkadian empire in Mesopotamia (Bar-Matthews and Ayalon, 2011; Cullen et al., 2000; Weiss et al., 1993). Similarly, a dry period inferred from reduced discharge of the Indus river and elevated δ18O of a Northeast-Indian stalagmite has been linked with the Indus Valley de-urbanization (Berkelhammer et al., 2013; Staubwasser et al., 2003; Staubwasser and Weiss, 2006). A recent study in the ancient Near East, however, revealed that the regional impact of the drought on ancient civilizations, being influenced by geographic factors and human technology, were highly diverse even within spatially limited cultural units (Riehla et al., 2014). This highlights the need for caution when linking evidence of human activities from a site to evidence of climate oscillations from another one.

The present study in the Gueldaman GLD1 Cave provides an opportunity to test climate–culture relations by comparing in situ archeological sequences and high-resolution paleoclimate records, thereby avoiding the uncertainty of inter-site correlation arising from complex spatial heterogeneity in climate and demography.

To facilitate the comparison, stalagmite-inferred climate changes at the cave site during ca. 6200–3200 yr BP are separated into four periods 1–4 (Table 3):

Period 1 (∼ 6200–5100 yr BP): wet, superimposed by several centennial-scale drier events;

Moreover, the sole piece of the bone from large ungulate (supposed to be elephant or rhinoceros) found at the depth of ∼ 110 cm (Kherbouche et al., 2014) was anchored by two calibrated 14C dates from the present study between 6003 and 4913 cal yr BP, which is in line with the latest survival of large mammal species (e.g., S. antiquus) at the proximate sites of Northern Algeria during the Mid-Holocene (Faith, 2014 and references wherein). The extinction of the Mid-Holocene large mammal in Northern Algeria was attributed to the competition with pastoralists and livestock for increasingly scarce water, corresponding with an abrupt climatic shift toward extreme aridity in the Sahara region ca. 5500 cal yr BP (i.e. the end of the Humid Africa Period (deMenocal et al., 2000; Faith, 2014). Recently, the timing of this climatic transition was refined to ca. 4900 yr BP ± 200 yrs (McGee et al., 2013). In addition, a paleoenvironmental study in the Sahara revealed that the Mid-Holocene deteriorations of terrestrial ecosystem and climate culminated at ca. 4200–3900 cal yr BP (e.g., Kröpelin et al., 2006). Therefore, it is more likely based on evidence from the Gueldaman GLD1 Cave and the proximate sites (Faith, 2014) that extinction of large mammal/ungulate in Northern Algeria occurred during the prolonged drought ca. 4400–3800 yr BP.