Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Mid-Pliocene glaciation preceded by a 0.5-million-year North African humid period

Abstract

Past North African humid periods caused expanded vegetation over the Sahara, due to northward tropical African rainbelt displacement, opening migration pathways for hominins. Commonly, these precession-timed humid periods ended within 15,000 years due to rainbelt retreat. During North African humid periods, eastern Mediterranean organic-rich layers called sapropels were deposited at least since 8 Myr. Here we combine climate modelling with palaeoclimate proxy data to show that weakened sapropel preservation during the 5.3–3.3 Myr period resulted from nutrient runoff limitation associated with enhanced North African vegetation cover due to a persistently more northward-located African monsoon front, relative to the mid-Pliocene (3.3–3.0 Myr, when glacial intensity increased). Moreover, sapropel absence within the 3.8–3.3 Myr period coincided with maximum monsoon runoff and extensively humid, vegetated conditions throughout North Africa. Our model results indicate that this 0.5-Myr-long pan-North African humid period ended at ~3.3 Myr because of southward monsoon front displacement with Northern Hemisphere glacial intensification. The 3.8–3.3 Myr humid period coincided with the earliest known evidence for hominin coexistence over eastern and central North Africa. We posit that persistent green corridors during this humid phase facilitated early hominin connectivity and migration, expanding their habitat range over the wider North African territory.

This is a preview of subscription content, access via your institution

Access options

Buy this article

USD 39.95

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Eastern Mediterranean sapropels across the past 5 Myr.
Fig. 2: North Africa drying during the M2 glaciation.
Fig. 3: Evidence for enhanced humid conditions across North Africa during the 3.3–3.8 Myr interval.
Fig. 4: The Pan-NAHP and hominin evolution.

Similar content being viewed by others

Data availability

ODP 967 Nd isotope data from this study have been submitted to Pangaea (www.pangaea.de) under ‘Plio-Pleistocene Neodymium (Nd) isotope data for Ocean Drilling Program (ODP) Site 967, Eastern Mediterranean’ and are also available as online Supplementary Data accompanying this article. The model data presented in this study are available to download via Zenodo at https://doi.org/10.5281/zenodo.11124493 (ref. 69). Source data are provided with this paper.

Code availability

The model code used in this study is available to download via Github at https://github.com/mom-ocean/MOM5 using the CM2M configuration under version 5.1.0. Input files and restart files for each experiment are available via Zenodo at https://doi.org/10.5281/zenodo.11124493 (ref. 69).

References

  1. Crocker, A. J. et al. Astronomically controlled aridity in the Sahara since at least 11 million years ago. Nat. Geosci. 15, 671–676 (2022).

    Article  CAS  Google Scholar 

  2. Skonieczny, C. et al. African humid periods triggered the reactivation of a large river system in western Sahara. Nat. Commun. 6, 8751 (2015).

    Article  CAS  Google Scholar 

  3. Menviel, L. et al. Drivers of the evolution and amplitude of African Humid Periods. Commun. Earth Environ. 2, 237 (2021).

    Article  Google Scholar 

  4. deMenocal, P. et al. Abrupt onset and termination of the African Humid Period: rapid climate responses to gradual insolation forcing. Quat. Sci. Rev. 19, 347–361 (2000).

    Article  Google Scholar 

  5. Larrasoaña, J. C., Roberts, A. P. & Rohling, E. J. Dynamics of Green Sahara Periods and their role in Hominin evolution. PLoS ONE 8, e76514 (2013).

    Article  Google Scholar 

  6. Kuper, R. & Kröpelin, S. Climate-controlled Holocene occupation in the Sahara: motor of Africa’s evolution. Science 313, 803–807 (2006).

    Article  CAS  Google Scholar 

  7. Osborne, A. H. et al. A humid corridor across the Sahara for the migration of early modern humans out of Africa 120,000 years ago. Proc. Natl Acad. Sci. USA 105, 16444–16447 (2008).

    Article  Google Scholar 

  8. Shanahan, T. M. et al. The time-transgressive termination of the African Humid Period. Nat. Geosci. 8, 140–144 (2015).

    Article  CAS  Google Scholar 

  9. Kröpelin, S. et al. Climate-driven ecosystem succession in the Sahara: the past 6000 years. Science 320, 765–768 (2008).

    Article  Google Scholar 

  10. Grant, K. M. et al. Organic carbon burial in Mediterranean sapropels intensified during Green Sahara Periods since 3.2 Myr ago. Commun. Earth Environ. 3, 11 (2022).

    Article  Google Scholar 

  11. Amarathunga, U. et al. Sill-controlled salinity contrasts followed post-Messinian flooding of the Mediterranean. Nat. Geosci. 15, 720–725 (2022).

    Article  CAS  Google Scholar 

  12. Rossignol-Strick, M. African monsoons, an immediate climate response to orbital insolation. Nature 304, 46–49 (1983).

    Article  Google Scholar 

  13. Rohling, E. J., Marino, G. & Grant, K. M. Mediterranean climate and oceanography, and the periodic development of anoxic events (sapropels). Earth Sci. Rev. 143, 62–97 (2015).

    Article  CAS  Google Scholar 

  14. Rohling, E. J. et al. Sea-level and deep-sea-temperature variability over the past 5.3 million years. Nature 508, 477–482 (2014).

    Article  CAS  Google Scholar 

  15. De Lange, G. J. et al. Synchronous basin-wide formation and redox-controlled preservation of a Mediterranean sapropel. Nat. Geosci. 1, 606–610 (2008).

    Article  Google Scholar 

  16. Emeis, K.-C., Sakamoto, T., Wehausen, R. & Brumsack, H.-J. The sapropel record of the eastern Mediterranean Sea—results of Ocean Drilling Program Leg 160. Palaeogeogr. Palaeoclimatol. Palaeoecol. 158, 371–395 (2000).

    Article  Google Scholar 

  17. Garcin, Y. et al. Short-lived increase in erosion during the African Humid Period: evidence from the northern Kenya Rift. Earth Planet. Sci. Lett. 459, 58–69 (2017).

    Article  CAS  Google Scholar 

  18. Rohling, E. J. et al. Sea level and deep-sea temperature reconstructions suggest quasi-stable states and critical transitions over the past 40 million years. Sci. Adv. 7, eabf5326 (2021).

    Article  CAS  Google Scholar 

  19. Routson, C. C. et al. Mid-latitude net precipitation decreased with Arctic warming during the Holocene. Nature 568, 83–87 (2019).

    Article  CAS  Google Scholar 

  20. Chiang, J. C. H. & Friedman, A. R. Extratropical cooling, interhemispheric thermal gradients, and tropical climate change. Annu. Rev. Earth Planet. Sci. 40, 383–412 (2012).

    Article  CAS  Google Scholar 

  21. Schneider, T., Bischoff, T. & Haug, G. H. Migrations and dynamics of the intertropical convergence zone. Nature 513, 45–53 (2014).

    Article  CAS  Google Scholar 

  22. Stepanek, C., Samakinwa, E., Knorr, G. & Lohmann, G. Contribution of the coupled atmosphere–ocean–sea ice–vegetation model COSMOS to the PlioMIP2. Clim 16, 2275–2323 (2020).

    Google Scholar 

  23. Willeit, M., Ganopolski, A. & Feulner, G. On the effect of orbital forcing on mid-Pliocene climate, vegetation and ice sheets. Clim 9, 1749–1759 (2013).

    Google Scholar 

  24. Heslop, D., Amarathunga, U. & Rohling, E. J. Estimating Plio‐Pleistocene North African monsoon runoff into the Mediterranean Sea and temperature impacts. Paleoceanogr. Paleoclimatol. 38, e2023PA004677 (2023).

    Article  Google Scholar 

  25. Laskar, J., Fienga, A., Gastineau, M. & Manche, H. La2010: a new orbital solution for the long-term motion of the Earth. Astron. Astrophys. 532, A89 (2011).

    Article  Google Scholar 

  26. Beuscher, S. et al. End-member modelling as a tool for climate reconstruction—an Eastern Mediterranean case study. PLoS ONE 12, e0185136 (2017).

    Article  Google Scholar 

  27. Garzanti, E., Andò, S., Padoan, M., Vezzoli, G. & El Kammar, A. The modern Nile sediment system: processes and products. Quat. Sci. Rev. 130, 9–56 (2015).

    Article  Google Scholar 

  28. Bastian, L. et al. Co-variations of climate and silicate weathering in the Nile Basin during the Late Pleistocene. Quat. Sci. Rev. 264, 107012 (2021).

    Article  Google Scholar 

  29. Scheuvens, D., Schütz, L., Kandler, K., Ebert, M. & Weinbruch, S. Bulk composition of northern African dust and its source sediments—a compilation. Earth Sci. Rev. 116, 170–194 (2013).

    Article  CAS  Google Scholar 

  30. Jewell, A. M. et al. Three North African dust source areas and their geochemical fingerprint. Earth Planet. Sci. Lett. 554, 116645 (2021).

    Article  CAS  Google Scholar 

  31. Trauth, M. H., Maslin, M. A., Deino, A. & Strecker, M. R. Late Cenozoic moisture history of East Africa. Science 309, 2051–2053 (2005).

    Article  CAS  Google Scholar 

  32. Maslin, M. A. et al. East African climate pulses and early human evolution. Quat. Sci. Rev. 101, 1–17 (2014).

    Article  Google Scholar 

  33. Joordens, J. C. A., Feibel, C. S., Vonhof, H. B., Schulp, A. S. & Kroon, D. Relevance of the eastern African coastal forest for early hominin biogeography. J. Hum. Evol. 131, 176–202 (2019).

    Article  Google Scholar 

  34. Hernández Fernández, M. & Vrba, E. S. Plio-Pleistocene climatic change in the Turkana Basin (East Africa): evidence from large mammal faunas. J. Hum. Evol. 50, 595–626 (2006).

    Article  Google Scholar 

  35. Cerling, T. E. et al. Woody cover and hominin environments in the past 6 million years. Nature 476, 51–56 (2011).

    Article  CAS  Google Scholar 

  36. Levin, N. E. Environment and climate of early human evolution. Annu. Rev. Earth Planet. Sci. 43, 405–429 (2015).

    Article  CAS  Google Scholar 

  37. Lebatard, A.-E. et al. Application of the authigenic 10Be/9Be dating method to continental sediments: reconstruction of the Mio-Pleistocene sedimentary sequence in the early hominid fossiliferous areas of the northern Chad Basin. Earth Planet. Sci. Lett. 297, 57–70 (2010).

    Article  CAS  Google Scholar 

  38. Moussa, A. et al. Lake Chad sedimentation and environments during the late Miocene and Pliocene: new evidence from mineralogy and chemistry of the Bol core sediments. J. Afr. Earth Sci. 118, 192–204 (2016).

    Article  CAS  Google Scholar 

  39. Brierley, C. M. et al. Greatly expanded tropical warm pool and weakened Hadley circulation in the early Pliocene. Science 323, 1714–1718 (2009).

    Article  CAS  Google Scholar 

  40. Herbert, T. D., Peterson, L. C., Lawrence, K. T. & Liu, Z. Tropical ocean temperatures over the past 3.5 million years. Science 328, 1530–1534 (2010).

    Article  CAS  Google Scholar 

  41. Herbert, T. D. et al. Late Miocene global cooling and the rise of modern ecosystems. Nat. Geosci. 9, 843–847 (2016).

    Article  CAS  Google Scholar 

  42. Park, J.-Y., Bader, J. & Matei, D. Northern-hemispheric differential warming is the key to understanding the discrepancies in the projected Sahel rainfall. Nat. Commun. 6, 5985 (2015).

    Article  CAS  Google Scholar 

  43. Haile-Selassie, Y. et al. New species from Ethiopia further expands Middle Pliocene hominin diversity. Nature 521, 483–488 (2015).

    Article  CAS  Google Scholar 

  44. Wood, B. & Boyle, E. K. Hominin taxic diversity: fact or fantasy?: Hominin taxic diversity. Am. J. Phys. Anthropol. 159, 37–78 (2016).

    Article  Google Scholar 

  45. Guy, F. et al. Symphyseal shape variation in extant and fossil hominoids, and the symphysis of Australopithecus bahrelghazali. J. Hum. Evol. 55, 37–47 (2008).

    Article  Google Scholar 

  46. Maxwell, S. J., Hopley, P. J., Upchurch, P. & Soligo, C. Sporadic sampling, not climatic forcing, drives observed early hominin diversity. Proc. Natl Acad. Sci. USA 115, 4891–4896 (2018).

    Article  CAS  Google Scholar 

  47. Brunet, M. et al. The first australopithecine 2,500 kilometres west of the Rift Valley (Chad). Nature 378, 273–275 (1995).

    Article  CAS  Google Scholar 

  48. Harmand, S. et al. 3.3-million-year-old stone tools from Lomekwi 3, West Turkana, Kenya. Nature 521, 310–315 (2015).

    Article  CAS  Google Scholar 

  49. Zeller, E. et al. Human adaptation to diverse biomes over the past 3 million years. Science 380, 604–608 (2023).

    Article  CAS  Google Scholar 

  50. Cuthbert, M. O. et al. Modelling the role of groundwater hydro-refugia in East African hominin evolution and dispersal. Nat. Commun. 8, 15696 (2017).

    Article  CAS  Google Scholar 

  51. Murat, A. & Got, H. Organic carbon variations of the eastern Mediterranean Holocene sapropel: a key for understanding formation processes. Palaeogeogr. Palaeoclimatol. Palaeoecol. 158, 241–257 (2000).

    Article  Google Scholar 

  52. Nijenhuis, I. A. & de Lange, G. J. Geochemical constraints on Pliocene sapropel formation in the eastern Mediterranean. Mar. Geol. 163, 41–63 (2000).

    Article  CAS  Google Scholar 

  53. Castradori, D. Calcareous nannofossils and the origin of eastern Mediterranean sapropels. Paleoceanography 8, 459–471 (1993).

    Article  Google Scholar 

  54. Dymond, J., Suess, E. & Lyle, M. Barium in deep-sea sediment: a geochemical proxy for paleoproductivity. Paleoceanography 7, 163–181 (1992).

    Article  Google Scholar 

  55. Van Os, B. J. H., Lourens, L. J., Hilgen, F. J., De Lange, G. J. & Beaufort, L. The formation of Pliocene sapropels and carbonate cycles in the Mediterranean: diagenesis, dilution, and productivity. Paleoceanography 9, 601–617 (1994).

    Article  Google Scholar 

  56. Möbius, J., Lahajnar, N. & Emeis, K.-C. Diagenetic control of nitrogen isotope ratios in Holocene sapropels and recent sediments from the Eastern Mediterranean Sea. Biogeosciences 7, 3901–3914 (2010).

    Article  Google Scholar 

  57. Jung, M., Ilmberger, J., Mangini, A. & Emeis, K.-C. Why some Mediterranean sapropels survived burn-down (and others did not). Mar. Geol. 141, 51–60 (1997).

    Article  CAS  Google Scholar 

  58. Emeis, K.-C., Robertson, A. H. & Richter, C. Site 966, initial reports. Proc. Ocean Drill. Program 160, 155–213 (1996).

  59. Delworth, T. L. et al. GFDL’s CM2 global coupled climate models. Part I: formulation and simulation characteristics. J. Clim. 19, 643–674 (2006).

    Article  Google Scholar 

  60. Haywood, A. M. et al. The Pliocene Model Intercomparison Project (PlioMIP) Phase 2: scientific objectives and experimental design. Clim. Past 12, 663–675 (2016).

  61. Dowsett, H. et al. The PRISM4 (mid-Piacenzian) paleoenvironmental reconstruction. Clim. Past 12, 1519–1538 (2016).

  62. Haywood, A. M. et al. Pliocene Model Intercomparison Project Phase 3 (PlioMIP3)—science plan and experimental design. Glob. Planet. Change 232, 104316 (2024).

    Article  Google Scholar 

  63. Dolan, A. M. et al. Modelling the enigmatic Late Pliocene Glacial Event—Marine Isotope Stage M2. Glob. Planet. Change 128, 47–60 (2015).

    Article  Google Scholar 

  64. Roy, K. & Peltier, W. R. Relative sea level in the Western Mediterranean basin: a regional test of the ICE-7G_NA (VM7) model and a constraint on late Holocene Antarctic deglaciation. Quat. Sci. Rev. 183, 76–87 (2018).

    Article  Google Scholar 

  65. Golledge, N. R. et al. The multi-millennial Antarctic commitment to future sea-level rise. Nature 526, 421–425 (2015).

    Article  CAS  Google Scholar 

  66. Feng, R. et al. Past terrestrial hydroclimate sensitivity controlled by Earth system feedbacks. Nat. Commun. 13, 1306 (2022).

    Article  CAS  Google Scholar 

  67. Prescott, C. L., Dolan, A. M., Haywood, A. M., Hunter, S. J. & Tindall, J. C. Regional climate and vegetation response to orbital forcing within the mid-Pliocene warm period: a study using HadCM3. Glob. Planet. Change 161, 231–243 (2018).

    Article  Google Scholar 

  68. Salzmann, U., Haywood, A. M., Lunt, D. J., Valdes, P. J. & Hill, D. J. A new global biome reconstruction and data‐model comparison for the Middle Pliocene. Glob. Ecol. Biogeogr. 17, 432–447 (2008).

    Article  Google Scholar 

  69. Hutchinson, D. K. & Amarathunga, U. Data for manuscript ‘A 0.5-million-year North African humid period preceded the mid-Pliocene glaciation’. Zenodo https://doi.org/10.5281/zenodo.11124493 (2024).

  70. Francke, A., Carney, S., Wilcox, P. & Dosseto, A. Sample preparation for determination of comminution ages in lacustrine and marine sediments. Chem. Geol. 479, 123–135 (2018).

    Article  CAS  Google Scholar 

  71. Rohling, E. J. et al. Reconstructing past planktic foraminiferal habitats using stable isotope data: a case history for Mediterranean sapropel S5. Mar. Micropaleontol. 50, 89–123 (2004).

    Article  Google Scholar 

  72. Amies, J. D., Rohling, E. J., Grant, K. M., Rodríguez‐Sanz, L. & Marino, G. Quantification of African monsoon runoff during last interglacial sapropel S5. Paleoceanogr. Paleoclimatol. 34, 1487–1516 (2019).

    Article  Google Scholar 

  73. Athanasiou, M. et al. Sea surface temperatures and environmental conditions during the ‘warm Pliocene’ interval (~4.1–3.2 Ma) in the Eastern Mediterranean (Cyprus). Glob. Planet. Change 150, 46–57 (2017).

    Article  Google Scholar 

  74. Khélifi, N., Sarnthein, M., Frank, M., Andersen, N. & Garbe-Schönberg, D. Late Pliocene variations of the Mediterranean outflow. Mar. Geol. 357, 182–194 (2014).

    Article  Google Scholar 

  75. Mulitza, S. et al. Increase in African dust flux at the onset of commercial agriculture in the Sahel region. Nature 466, 226–228 (2010).

    Article  CAS  Google Scholar 

  76. O’Mara, N. A. et al. Pleistocene drivers of Northwest African hydroclimate and vegetation. Nat. Commun. 13, 3552 (2022).

    Article  Google Scholar 

  77. Larrasoaña, J. C., Roberts, A. P., Rohling, E. J., Winklhofer, M. & Wehausen, R. Three million years of monsoon variability over the northern Sahara. Clim. Dyn. 21, 689–698 (2003).

    Article  Google Scholar 

  78. Wu, J., Böning, P., Pahnke, K., Tachikawa, K. & De Lange, G. J. Unraveling North-African riverine and eolian contributions to central Mediterranean sediments during Holocene sapropel S1 formation. Quat. Sci. Rev. 152, 31–48 (2016).

    Article  Google Scholar 

  79. Wu, J., Filippidi, A., Davies, G. R. & De Lange, G. J. Riverine supply to the eastern Mediterranean during last interglacial sapropel S5 formation: a basin-wide perspective. Chem. Geol. 485, 74–89 (2018).

    Article  CAS  Google Scholar 

  80. Schlitzer, R. Ocean Data View https://odv.awi.de (2023).

  81. Fairhead, J. D. The Mesozoic West and Central Africa Rift System (WCARS) and the older Kandi Shear Zone (KSZ): rifting and tectonics of North Africa and South America and fragmentation of Gondwana based on geophysical investigations. J. Afr. Earth Sci. 199, 104817 (2023).

    Article  Google Scholar 

  82. Makeen, Y. M. et al. Sedimentology, petrography, and reservoir quality of the Zarga and Ghazal formations in the Keyi oilfield, Muglad Basin, Sudan. Sci. Rep. 11, 743 (2021).

    Article  CAS  Google Scholar 

  83. Macgregor, D. History of the development of the East African Rift System: a series of interpreted maps through time. J. Afr. Earth Sci. 101, 232–252 (2015).

    Article  Google Scholar 

Download references

Acknowledgements

This work contributes to Australian Research Council projects FL120100050 and DP200101157 (E.J.R.), DE190100042 (K.M.G.), DP190100874 (A.P.R., D.H.), DE220100279 (D.K.H.) and the Australia–New Zealand IODP Consortium (ANZIC) Legacy/Special Analytical Funding grant LE160100067 (K.M.G.). We thank Mawson Analytical Spectrometry Services, University of Adelaide for contributing to Nd isotope measurements.

Author information

Authors and Affiliations

Authors

Contributions

U.A. designed and led the study, prepared samples for Nd isotope analyses, analysed the data and climate model outputs, developed the hypothesis, generated the figures and wrote the paper; A.F. advised on Nd isotope analysis and data interpretation; R.M.K. carried out Nd isotope analyses; D.K.H. generated the climate model ice sheet configurations and performed modelling experiments; D.H. advised on statistical analyses and performed the box model experiments; K.M.G. and J.L. helped with sample selection and preparation; E.J.R. advised on data interpretation; all authors contributed to paper development.

Corresponding author

Correspondence to Udara Amarathunga.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Geoscience thanks Gert De Lange, J. Faith, Christian Stepanek and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Primary Handling Editor: James Super, in collaboration with the Nature Geoscience team.

Additional information

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

Extended data

Extended Data Fig. 1 Preserved sapropel record of the eastern Mediterranean.

a, Conceptual diagram with modern locations of eastern Mediterranean ODP sites. The present water depth at each site is indicated. b, Preserved sapropels at each site. The ODP Site 964 sediment record only goes back to ~3.8 Myr. At Site 966, sediments beyond 4.5 Myr are tectonically deformed16,58. Prolonged sapropel absence during a period of high-amplitude precession minima (centred at ~3.5 Myr) over the entire Mediterranean is marked by vertical blue shading. The ODP Site 967 Ba record indicates enhanced organic carbon burial after 3.2 Myr, which is observed at all sites. The presence of Ba/Al peaks (with a low magnitude) in the Early Pliocene suggests the presence of organic carbon burial, that were oxidized later. Cross section in a drawn using Ocean Data View version 5.6.3 (ref. 80).

Extended Data Fig. 2 Evidence for enhanced freshwater flux to the eastern Mediterranean during the Pan-NAHP.

a, Mediterranean δ18OCalcite sensitivity to sea surface temperature (SST), monsoon runoff and surface-water temperature concentration (Tc) changes (see Methods). For unchanging Gibraltar sill relative sea-level (RSLGib), 1 °C SST increase can cause a ~ 0.2 ‰ δ18O drop. A freshwater runoff increase equivalent to pre-Aswan Nile runoff can cause a ~ 0.28 ‰ δ18O drop. A 1 °C surface freshwater layer temperature increase can cause a ~ 0.21 ‰ δ18O drop. For example, with unchanging SST and temperature concentration across transition 1 (T1), the observed 0.45 ‰ shift to negative δ18O can be generated with a runoff increase of 1.6 times pre-Aswan Nile runoff equivalents. Shaded intervals indicate the standard error associated with sensitivity estimates. b, SST reconstructions available for the Mediterranean41,73,74. The mean SST variability does not indicate an SST rise that can explain the δ18O drop at T1 (See Methods). c, ODP 967 δ18O and global ice volume change18 across the 2.5-4.5 Myr period. The blue dotted line (upper diagram) comprises all δ18O data points across the interval. The thick grey line (lower diagram) indicates global ice volume evolution. Thick blue and black lines represent the 60-kyr moving average for the δ18O and ice volume curves, respectively. The shading for each curve represents 90% confidence intervals. The shifts to more negative δ18O values at T1 and more positive values at transition 2 (T2) are highlighted in grey vertical bars. While T2 was caused mainly by the global ice volume increase (a minimum of 0.3 ‰, see Methods), T1 was largely caused by a mean freshwater runoff increase to the eastern Mediterranean (no relative global ice volume change across T1, see Methods). Green boxes indicate data gaps.

Extended Data Fig. 3 Topography and bathymetry used for model simulations.

a, c, e, Topography used for the early Pliocene, M2 glacial and mid-Pliocene glacial simulations, respectively. The ice sheets were adjusted for each simulation, following the calculations based on ice volume reconstructions. b, d, f, Bathymetry for the early Pliocene, M2 glacial and mid-Pliocene glacial simulations, respectively. g, h, Difference in topography for 2 scenarios, M2 simulation-early Pliocene simulation and mid-Pliocene glacial simulation-early Pliocene simulation.

Extended Data Fig. 4 Climate model result comparison between the Early Pliocene and the M2 glacial.

a, Boreal summer 600 hPa vertical velocity change (M2 glacial–Early Pliocene), with decreased ascending motion over North Africa at ~20°-25°N during the M2 glacial. Contours represent Early Pliocene climatology; continuous red lines indicate descending motion (positive) and dashed blue lines indicate ascending motion (negative). Decreased northward ascending motion indicates a more southward monsoon front during the M2 glaciation, resulting in the observed precipitation decrease centred at ~20°N. b, The change in 850 hPa zonal wind field (M2 glacial–Early Pliocene), which indicates southward displacement of the easterly (negative/dashed blue contours) -westerly (positive/ solid red contours) convergence zone (solid grey contour) centred at ~18°N over North Africa. Contours indicate Early Pliocene climatology, whereas shading indicates the change from Early Pliocene to the M2 glacial.

Extended Data Fig. 5 Climate model result comparison between the Early Pliocene and mid-Pliocene glacials.

a, Boreal summer 600 hPa vertical velocity change (Mid Pliocene glacial mean-Early Pliocene mean), with decreased ascending motion over North Africa centred at ~20°N compared to the Early Pliocene. Contours represent Early Pliocene climatology; continuous red lines indicate descending motion (positive) and dashed blue lines indicate ascending motion (negative). Decreased northward ascending motion indicates a more southward monsoon front during mid-Pliocene glacials, resulting in the observed precipitation decrease centred at ~20°N. b, The change in 850 hPa zonal wind field (Mid Pliocene glacial mean-Early Pliocene mean), indicating a southward displacement of the easterly -westerly convergence zone centred at ~18°N over North Africa. Contour arrangement (type and colour) is similar to Extended Data Fig. 4. Vertical and horizontal wind vigour is smaller during the Mid Pliocene glacials compared to the M2 glaciation (see Extended Data Fig. 4).

Extended Data Fig. 6 Data for ODP sites 967 and 659.

a, Northern Hemispheric ice volume evolution through the past 10 Myr18. b, ODP 967 Ti/Al data, plotted separately for periods of sapropel (blue) and marl (orange) deposition10. Thick blue and orange lines represent the 60-kyr moving average for sapropel and marl intervals. Both sapropelic and marl Ti/Al values indicate an abrupt slope change to higher relative fluvial sediment inputs at 3.8 Myr. At 3.4 Myr, only the marl Ti/Al values indicate the return of increased aeolian input to the eastern Mediterranean. c, ODP 967 δ18O record back to the Miocene–Pliocene boundary. The thick blue line represents the 60-kyr moving average. d, ODP 659 [Al+Fe]/[Si+K+Ti] record back to 10 Myrref. 1. The thick brown line represents the 60-kyr moving average. MSC, Messinian salinity crisis11.

Extended Data Fig. 7 Nd isotope data for ODP Sites 967 and 659, compared with North African sediment source regions.

a, Map of the main aeolian and fluvial sediment source areas for the eastern Mediterranean (ODP Site 967) and west Africa (ODP Site 659). Land areas highlighted in yellow, purple and grey indicate preferential dust source areas (PSAs) for western, central and eastern North African regions (PSAWest, PSACen. and PSAEast, respectively) (data from ref. 30). Main modern fluvial sediment contributors to the eastern Mediterranean (Main, Blue and White Nile) are shown. b, Comparison of data for sites 967 and 659, with modern εNd values from each sediment source region shown in a. Shaded horizontal bars indicate the range of values observed from each source region. Dots within the bars represent εNd values for individual samples, retrieved from published databases. ODP Site 659 values mostly correspond to εNd values from PSAWest and PSACen, as interpreted in ref. 1. ODP Site 967 εNd data from the present study indicate a largely constant sediment origin through the past 5 Myr. Samples representing periods of enhanced organic carbon burial (black diamonds) prior to 3.2 Myr were chosen based on geochemical data (Ba/Al)10.

Extended Data Fig. 8 North Atlantic meridional temperature gradient evolution.

a, Northern Hemisphere ice volume evolution through the past 5 Myr18. Thick red line represents the 60-kyr moving average. Upper and lower shadings represent 95th percentile confidence intervals. b, Global sea level relative to present18. Red circles with error bars represent global mean sea-level (GMSL) benchmarks from Mallorca18. Error bars indicate age uncertainties (horizontal) and 16th and 84th percentile sea-level ranges (vertical). c, Meridional temperature gradient (ΔT) evolution between North Atlantic ODP Site 982 and available tropical Atlantic temperature records40,41. Thick lines represent the 60-kyr moving average for each record. Relatively stable Early Pliocene ΔT drops to the lowest values at 3.6-3.4 Myr. TS, Tropical stack. d, ODP 967 δ18O record back to the Miocene–Pliocene boundary. Thick black line represents the 60-kyr moving average. Upper and lower shadings represent 95th percentile confidence intervals. Blue arrows in c and d demarcate the stepwise ΔT drop and simultaneous δ18O decrease. The eastern Mediterranean sapropel stack is shown above the time axis. e, Map of locations of ODP sites used for the North Atlantic meridional temperature gradient evaluation. Orange labels indicate sites considered for the tropical stack (ODP sites 662, 722, 846, 1146 = alkenone-based temperature reconstructions; ODP Site 806 = Mg/Ca-based temperature reconstruction)40,41.

Extended Data Fig. 9 Pan-NAHP hominin sites and North African landscape.

a, Map of major hominin fossil sites for early hominin coexisting during the Pliocene Pan-NAHP. White dots indicate discovery sites43,44,45,46,47,48 of A. bahrelgazali (Koro Toro, Chad Basin), A. deyiremeda (Woranso–Mille, central Afar) and K. platyops (Lake Turkana, Kenya). Numbered black dots indicate A. afarensis discovery sites (1, Woranso–Mille; 2, Hadar; 3, Omo basin; 4, Koobi Fora; 5, Lomekwi; 6, Kantis; 7, Laetoli). East and Central African rift systems81,82,83 are highlighted in red and purple, respectively. Data from refs. 81,83. b, Cross-sections drawn to scale, along A-B-C and D-E lines on the map in a. The two sections have contrasting geographies along the East African rift valley and the Central African rift system. Lake names are in blue. Vertical grey shading indicates the intersection of two rift systems, where the Turkana depression is located32. The rift intersection was proposed as a bottleneck, where a possible hominin migration path toward the Chad Basin was proposed34 across the CASZ. CASZ, Central African shear zone.

Extended Data Table 1 Ice sheet boundary conditions used for the early Pliocene, M2 glacial and Mid Pliocene glacial simulations

Source data

Source Data Fig. 3 (download XLSX )

ODP 967 Nd isotope data.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Amarathunga, U., Rohling, E.J., Grant, K.M. et al. Mid-Pliocene glaciation preceded by a 0.5-million-year North African humid period. Nat. Geosci. 17, 660–666 (2024). https://doi.org/10.1038/s41561-024-01472-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • Version of record:

  • Issue date:

  • DOI: https://doi.org/10.1038/s41561-024-01472-8

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing