A distal 140-kyr sediment record of Nile discharge and East African monsoon variability

Abstract

Clay mineral assemblages in a sediment core from the distal Nile discharge plume off Israel have been used to reconstruct the late Quaternary Nile sediment discharge into the eastern Mediterranean Sea (EMS). The record spans the last ca. 140-kyr. Smectite abundances indicate the influence of the Blue Nile and the Atbara River that have their headwaters in the volcanic rocks of the Ethiopian Highlands. Kaolinite abundances indicate the influence of wadis, which contribute periodically to the suspension load of the Nile. <br><br> Due to the geographical position, the climate and the sedimentary framework of the EMS is controlled by two climate systems. The long-term climate regime was governed by the African monsoon that caused major African humid periods (AHPs) with enhanced sediment discharge at 132 to <-126 (AHP 5), 116 to 99 (AHP4), and 89 to 77-ka (AHP3). They lasted much longer than the formation of the related sapropel layers S5 (>-2-kyr), S4 (3.5-kyr), and S3 (5-kyr). During the last glacial period (Marine Isotope Stages (MISs) 4-2), the long-term changes in the monsoonal system were superimposed by millennial-scale changes in an intensified midlatitude glacial system. This climate regime caused short but pronounced drought periods in the Nile catchment, which are linked to Heinrich events and alternate with more humid interstadials. <br><br> The clay mineral record further implies that feedback mechanisms between vegetation cover and sediment discharge of the Nile are detectable but of minor importance for the sedimentary record in the southeastern Mediterranean Sea during the investigated African humid periods.

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/ doi:10.5194/cp-12-713-2016 Author(s) 2016. CC Attribution 3.0 License.

Werner Ehrmann1, Gerhard Schmiedl2, Martin Seidel1, Stefan Krger1, and Hartmut Schulz3

1Universitt Leipzig, Institut fr Geophysik und Geologie, Talstrae 35, 04103 Leipzig, Germany

2Universitt Hamburg, Centrum fr Erdsystemforschung und Nachhaltigkeit, Bundesstrae 55, 20146 Hamburg, Germany

3Universitt Tbingen, Fachbereich Geowissenschaften, Hlderlinstrae 12, 72074 Tbingen, Germany Correspondence to: Werner Ehrmann ([email protected])

Received: 20 July 2015 Published in Clim. Past Discuss.: 7 September 2015

Revised: 22 January 2016 Accepted: 17 February 2016 Published: 21 March 2016

Abstract. Clay mineral assemblages in a sediment core from the distal Nile discharge plume off Israel have been used to reconstruct the late Quaternary Nile sediment discharge into the eastern Mediterranean Sea (EMS). The record spans the last ca. 140 kyr. Smectite abundances indicate the inuence of the Blue Nile and the Atbara River that have their headwaters in the volcanic rocks of the Ethiopian Highlands. Kaolinite abundances indicate the inuence of wadis, which contribute periodically to the suspension load of the Nile.

Due to the geographical position, the climate and the sedimentary framework of the EMS is controlled by two climate systems. The long-term climate regime was governed by theAfrican monsoon that caused major African humid periods (AHPs) with enhanced sediment discharge at 132 to < 126 (AHP 5), 116 to 99 (AHP4), and 89 to 77 ka (AHP3). They lasted much longer than the formation of the related sapropel layers S5 (> 2 kyr), S4 (3.5 kyr), and S3 (5 kyr). During the last glacial period (Marine Isotope Stages (MISs) 42), the long-term changes in the monsoonal system were superimposed by millennial-scale changes in an intensied mid-latitude glacial system. This climate regime caused short but pronounced drought periods in the Nile catchment, which are linked to Heinrich events and alternate with more humid interstadials.

The clay mineral record further implies that feedback mechanisms between vegetation cover and sediment discharge of the Nile are detectable but of minor importance for the sedimentary record in the southeastern MediterraneanSea during the investigated African humid periods.

A distal 140 kyr sediment record of Nile discharge and East African monsoon variability

1 Introduction

The Nile in northern Africa is the longest river of the world and the dominant sediment source for the eastern Mediterranean Sea (EMS; Milliman and Syvitski, 1992). Its drainage basin is about 3 [notdef] 106 km2 (Garzanti et al., 2015) and ex

tends from the equator to ca. 30 N (Fig. 1). The main tributaries of the Nile are the perennial White Nile originating from Lake Victoria in tropical eastern Africa and the highly seasonal Blue Nile and the Atbara River originating from the subtropical Ethiopian Highlands. Downstream of the conuence, the Nile ows for more than 2000 km through the hyperarid Sahara.

Because of this special geographical situation, the runoff of the Nile is sensitive to climatic changes, and its discharge sediments are major recorders of the geological and climatic conditions in the catchment areas. The Nile sediments mainly reect the intensity of rainfall in the headwaters; this intensity has direct control on weathering, erosion, and transport of sediments (Krom et al., 2002; Revel et al., 2010; Box et al., 2011; Garzanti et al., 2015). The present-day summer oods of the Nile are linked to a northward movement of the Intertropical Convergence Zone (ITCZ) and the African rain belt (ARB; Fig. 1), causing especially intense precipitation in July and August in the Ethiopian Highlands.

Prior to the emplacement of the Aswan High Dam in 1964, the Nile carried a suspension load of about 120 160 [notdef] 106 t yr1 to the EMS (Holeman, 1968; Milliman and

Syvitski, 1992; Stanley and Wingerath, 1996a). More than 95 % of the material is derived from the Ethiopian Highlands via the Blue Nile and the Atbara River. In contrast, the sed-

Published by Copernicus Publications on behalf of the European Geosciences Union.

714 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

and Stanley, 1989; Williams et al., 2006, 2015; Garzanti et al., 2015). About 32.5 % of the Nile suspension consists of clay (Quelennec and Kruc, 1976). Most of the silt- and sand-sized sediment accumulates in the Nile Delta, on the Nile cone, and along the Mediterranean shelf, whereas the clay-sized sediment fraction is transported in suspension by the surface currents of the Mediterranean Sea to the east and north (Fig. 1; Venkatarathnam and Ryan, 1971; Foucault and Mlires, 2000; Hamann et al., 2009).

It is well known that past changes in the amount of rainfall in the Ethiopian Highlands were caused by precession-driven shifts in the position of the ITCZ that inuenced the intensity and the spatial extent of the monsoon (e.g.Rossignol-Strick, 1983; Emeis et al., 2003). During African humid periods (AHPs) the spreading of the vegetation cover in northern Africa resulted in a green Sahara (Renssen et al., 2006; Krpelin et al., 2008) and a reduced inux of Saharan dust into the Mediterranean Sea (deMenocal et al., 2000;Ehrmann et al., 2013). The increased runoff of the Nile River and of other northern African river and wadi systems brought huge amounts of freshwater and nutrients to the Mediterranean Sea, leading to enhanced productivity in the surface water, stagnation in the deep marine basin, and formation of sapropel layers (Rossignol-Strick et al., 1982; Rohling, 1994;Emeis et al., 2003; Rohling et al., 2015). Although numerous studies concentrated on the Holocene AHP, very little is known about timing and intensity of former AHPs. Marine proxy data across sapropel S5 suggest that during the Eemian period (Marine Isotope Stage (MIS) 5e) the African monsoon was particularly strong and the African rain belt was shifted even further north than during the early Holocene. This resulted in ooding into the EMS through wadi systems along the wider North African margin (Rohling et al., 2002a; Osborne et al., 2008).

Different opinions exist on whether Nile sediment discharge is primarily controlled by precipitation and runoff, and thus closely linked to changes in the monsoon intensity (Wehausen and Brumsack, 2000; Emeis et al., 2000; Revel et al., 2010), or whether feedback mechanisms play an important role, with vegetation protecting soils and weathering products in the source area from being eroded (Adamson et al., 1980; Krom et al., 2002; Box et al., 2011).

A comprehensive effort to reconstruct Nile discharge during the last 100 kyr was undertaken using the Fe content in sediments recovered from the proximal Nile Delta (Revel et al., 2010; Caley et al., 2011). Most other geochemical investigations concentrated on shorter time spans and used different parameters such as K / Al, Mg / Al, Ti / Al, Ba / Ca to characterize the Nile discharge (Wehausen and Brumsack, 2000;Box et al., 2011; Weldeab et al., 2014). Several studies include the isotopic ngerprint of Nile sediments to reconstruct Nile history (Krom et al., 1999a; Weldeab et al., 2002; Box et al., 2011; Blanchet et al., 2014).

We choose a different approach to reconstructing the Nile River sediment supply throughout the last 140 kyr by mainly

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

30 32 34 36

Levantine Sea

+ MD84-641

GeoT

SL110 +

SL112

Haifa

+ KL83

9509

+ MS27PT

Soreq

+ + Cave

Gaza

R.B.

D.B.

(b)

Nile

33

32

31

30

Mediterranean Sea

GeoT SL110

Nile

White Nile

Eastern

Sahara

desert

The Sudd

Aswan Red

Sea

ARB-S

ARB-W

Atbara

Lake Tana

Blue Nile

Basalt

Sobat

Lake Turkana

Lake Albert

Lake Edward

Lake Victoria

(a)

30

20

0

10

30

40

Figure 1. Location map of the Nile River basin and the eastern Mediterranean Sea. Panel (a): the Nile River and its main tributaries. Blue dashed lines indicate wadis. The green dotted line shows the main outcrop of Cenozoic basalts in Ethiopia. Also shown are the northern summer and winter limits of the African rain belt (red dashed lines; ARB-S, ARB-W, after Blanchet et al., 2014). Black arrows in the eastern Mediterranean Sea indicate the surface currents (after Pickard and Emery, 1990). Panel (b): map of the southeastern Levantine Sea with the location of the investigated sediment core GeoT SL110 and other cores mentioned in the text; not shown are GeoT SL143 and LC21, which are positioned further to the north (cf. Table 1). D.B.: Damietta branch of the Nile in the delta; R.B.: Rosetta branch.

iment discharge of the White Nile can be neglected because its load is stored in Lake Victoria, Lake Albert, and in the Sudd Basin of South Sudan (Adamson et al., 1980; Foucault

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 715

investigating the content of typical Nile-derived clay minerals in sediment core GeoT SL110, which was recovered from the distal Nile discharge plume off Israel (Fig. 1). Some 80 % of the terrigenous sediment fraction of the surface sediments near SL110 is derived from the Nile (Krom et al., 1999b) and is characterized by a high smectite content (e.g.Venkatarathnam and Ryan, 1971; Foucault and Mlires, 2000). By analysing a distal core from the plume, we document the Nile discharge activity as a whole rather than the activity of only a single channel in the delta.

2 Material

The investigated sediment core GeoT SL110 was recovered from the SE Levantine Sea (Fig. 1, Table 1) during cruise M51/3 of the German research vessel Meteor in 2001 (Hemleben et al., 2003). It was retrieved at the Israeli continental slope off Haifa from a water depth of 1437 m. The sediments of the 652 cm long core mainly consist of muds with traces of pteropods and foraminifers. The core top contains basically undisturbed surface sediments as indicated by the presence of an oxidized layer. The muds show frequent diffuse colour changes, mainly between grey and olive grey. Darker sapropel layers (greenish black, dark greenish grey, dark grey) occur at 3561 (S1, with a possible interruption at 4445 cm), 500514 (S3), 558570 (S4), and 598622 cm (S5, including pre-sapropelic layer). The core does not show any indications of slumping, debris ows, or turbidites but indicates rather constant sedimentary conditions in the distal discharge plume of the Nile. No tephra layers were detected.

3 Methods

The lightness (L ) of core SL110 was determined in 1 cm steps with a Minolta colour spectrophotometer immediately after opening the core.

The content of total organic carbon (TOC) was measured on 174 ground bulk sediment samples with an Eltra METALYST-CS-1000-S after removal of carbonate with HCl.

The core was sampled at 1 cm intervals for investigations of the grain size composition of the terrigenous sediment fraction and of the clay mineral assemblages. The samples were oxidized and disaggregated by means of 5 % hydrogen peroxide. Carbonate was dissolved by 10 % acetic acid. We isolated the terrigenous sand fraction by sieving the samples through a 63 m mesh. The silt fraction (263 m) was separated from the clay fraction (< 2 m) in settling tubes.

The analyses of the clay mineral composition followed standard methods (e.g. Ehrmann et al., 2007). We added MoS2 as an internal standard to the clay suspension. Texturally oriented clay mounts were solvated with ethylene glycol vapour at a temperature of 60 C and then X-rayed with a Rigaku MiniFlex system (CoK radiation; 30 kV;

15 mA). We analysed the samples in the range 340 2[Theta1] with a step size of 0.02 2[Theta1] and a measuring time of 2 s step1. In addition, we analysed the range 27.530.6 2[Theta1] with a step size of 0.01 2[Theta1] and a measuring time of 4 s step1 in order to better resolve the (002) peak of kaolinite and the (004)

peak of chlorite. We evaluated the diffractograms by using the MacDiff software (Petschick, 2001). After adjusting the diffractograms to the MoS2 peak at 6.15 , we deconvoluted the peak doublets smectitechlorite (17 and 14 ), palygorskiteillite (10.5 and 10 ), and kaolinitechlorite(3.58 and 3.54 ). We based the semi-quantitative estimations of the clay mineral abundances on the integrated peak areas of the individual peaks and weighting factors (Biscaye, 1964, 1965; Weaver and Beck, 1977).

4 Age model

In order to establish a stable oxygen isotope stratigraphy, we analysed tests from the 250355 m size fraction of the surface-dwelling foraminifer species Globigerinoides ruber (white) of about 130 samples using a Kiel IV online carbonate preparation line connected to a MAT 253 mass spectrometer. The values are reported in per mill relative to VPDB using the delta notation. The reproducibility was better than

0.06 (1). Age points were derived by a graphic correlation of our 18O record with the G. ruber (white) 18O data of sediment core LC21 from the southern Aegean Sea (Grant et al., 2012), aided by the software AnalySeries 2.0 (Paillard et al., 1996) (Table 2; Fig. 2). The LC21 chronology was constructed by using 14C dating and by tuning its 18O record to the U / Th-dated Soreq Cave speleothem 18O record (Grant et al., 2012).

Further age control comes from four 14C accelerator mass spectrometry (AMS) datings performed by the Beta Analytic Radiocarbon Dating Laboratory on well-preserved shells of planktonic foraminifera (G. ruber, G. bulloides, G. sacculifer, O. universa) that represent the age of the surface waters. We applied an eastern Levantine Sea Delta-R of 3

66 years (Marine Reservoir Database). The radiocarbon ages were converted to calendar years using the Marine13 database (Reimer et al., 2013; Table 2).

We identied a hiatus at 598 cm, at the top of the S5 interval, based on several arguments. The lower part of the dark interval between 622 and 604 cm is characterized by only about 0.8 % TOC (Fig. 3a), by the occurrence of a relatively diverse benthic foraminiferal assemblage containingG. orbicularis, G. translucens, and C. laevigata, and by relatively high 18O values (Fig. 3a). Despite the dark colour, this indicates a pre-sapropelic interval at the transition between MISs 6 and 5 rather than a full sapropel (Schmiedl et al., 2003). The TOC concentrations start to rise to ca. 1.4 % at 603 cm. We therefore assume the base of sapropel S5 at 603 cm, dated to 128.0 ka (Grant et al., 2012). The TOC values drop abruptly to 0.5 % at 596 cm. Thus, the TOC con-

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

716 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

Table 1. Metadata for the investigated sediment core GeoT SL110 and other cores mentioned in the text (cf. Fig. 1).Core Latitude Longitude Water depth Reference

( N) ( E) (m)

GeoT SL110 32 38.95[prime] 34 06.22[prime] 1437 this studyGeoT KL83 32 36.87[prime] 34 08.89[prime] 1431 Weldeab et al. (2002)

GeoT SL112 32 44.52[prime] 34 39.02[prime] 892 Hamann et al. (2008, 2009)

GeoT SL143 38 15.71[prime] 25 06.19[prime] 665 Ehrmann et al. (2013)

MD84-641 33 02[prime] 32 38[prime] 1375 Fontugne and Calvert (1992) 9501 34 32[prime] 33 59[prime] 980 Almogi-Labin et al. (2009)

9509 32 01[prime] 34 16[prime] 884 Almogi-Labin et al. (2009)

MS27PT 31 47.90[prime] 29 27.70[prime] 1389 Revel et al. (2010)

LC21 35 40[prime] 26 35[prime] 1522 Grant et al. (2012)

Table 2. Data used for constructing the age model for the investigated sediment core GeoT SL110. Oxygen isotope ages are based on a graphical correlation with the LC21 record (Grant et al., 2012) from the southern Aegean Sea.

Depth Age Data (cm) (date cal. ka BP)

0.00 0.00 Sediment surface34.50 6.43 14C AMS dating

64.50 11.53 14C AMS dating

150.50 24.38 14C AMS dating 189.11 29.29 LC21234.50 33.92 14C AMS dating 254.50 36.29 LC21359.54 53.50 LC21465.30 67.11 LC21509.55 84.99 LC21529.23 89.03 LC21554.51 103.32 LC21569.46 107.79 LC21579.50 110.38 LC21598.00 119.00 Hiatus end 598.00 126.00 Hiatus start 603.00 128.00 LC21624.52 131.94 LC21

Age (ka)

0 20 40 60 80 100 120 140

0

50

150

250

Depth (cm)

100

200

300

400

500

600

650

550

MD84641

350

450

GeoT SL110

MIS

centrations remain below the typical values of 23 % in S1, S3, S4 (Fig. 3a), and in S5 of nearby core KL83 (Weldeab et al., 2003). Also the analysis of the benthic foraminifera showed an extremely thin anoxic phase in SL110, compared to a 27 cm thick S5 in KL83 (Schmiedl et al., 2003). Thus, most of S5 is missing in SL110 due to a hiatus that was probably caused by a slump as indicated by the sharp and curved lithological boundary at 598 cm. By comparing our 18O record with standard curves and the record of MD84-641 (Table 1; Kallel et al., 2000) and considering a duration of S5 spanning the period 128121 ka (Grant et al., 2012), we argue for an extrapolated date of the hiatus from 126 to 119 ka. In nearby core GeoT KL83 (Table 1), Weldeab et

1 2 3 4 5 6

Figure 2. Depthage plot for the investigated sediment core GeoT SL110 (red). The age points are listed in Table 2. For comparison, corresponding data were extracted for core MD84-641 (blue;Fontugne and Calvert, 1992: Table 1). Marine isotope stages (MISs) follow Lisiecki and Raymo (2005).

al. (2002, 2003) identied a hiatus spanning the period 110 70 ka, with the base of the hiatus at 13 cm above S5.

For consistency reasons, we also adapted the age model of core GeoT SL143 (Ehrmann et al., 2013; Fig. 4) by using basal dates of 86.0 and 108.0 ka for sapropels S3 and S4, respectively, according to the LC21 record (Grant et al., 2012)

5 Results

The raw data of our investigations on sediment core GeoT SL110 are presented in Fig. 3, which also shows the positions of the sapropel layers S1, S3, S4, and S5.

The 18O data of G. ruber go back to MIS 6 (Fig. 3a). The Last Glacial Maximum has values of ca. 3 , whereas the minimum values during interglacials MIS 1 and MIS 5 are around 1.5 .

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 717

Figure 3. Basic data of core GeoT SL110. Panel (a): 18O of G. ruber, lightness (L ), % total organic carbon (TOC), and % sand and silt / clay ratio of the terrigenous sediment fraction, and linear sedimentation rates (LSRs; cm ka1). Panel (b): clay mineral percentages of smectite, illite, kaolinite, chlorite, and palygorskite. Grey bars mark the position of sapropel layers S1, S3, S4, S5, and the pre-sapropelic layer beneath S5 (light grey). Marine isotope stages (MISs) are indicated.

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

718 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

Age (ka)

0 10 40 60

30 50

20 70 80 110 130

100 120

90 140

500 400 Insolation (W/m)

2

((a)

(b)

(c)

(d)

(e)

12345

Smectite / (Illite + Chlorite)

0.8 1.2 1.4 1.0

Kaol / Chl SL143

Saharan dust

1.5

1.0 0.5 0.0

1.5 1.0

+umid

$ridHumid

Arid

XRF Fe / total counts

Kaolinite / (Illite + Chlorite)

Hiatus

Hiatus

S 1

S 3

S 4

S 5

H6

H1 ?

AHP1 AHP2 AHP3 AHP4 AHP5

H2

H3

H5 H4

1

2

3

4

5

6

MIS

0 10 40 60

30 50

20 70 80 110 130

100 120

90 140

Age (ka)

Figure 4. Combination of the June insolation at 65 N (a), the Saharan dust record in the central Aegean Sea core GeoT SL143 (b; Ehrmann et al., 2013), the Nile discharge based on Fe data in the Nile Delta (c; Revel et al., 2010; age model by Caley et al., 2011), and the Nile sediment discharge based on Kar (d) and Smr (e) at site GeoT SL110. Grey bars indicate the positions of sapropels S1, S3, S4, S5, and of the pre-sapropelic layer (light grey) beneath S5 in core in SL110. Heinrich events (H) are agged with arrows. Marine isotope stages (MISs)

and African humid periods (AHPs) are indicated at the bottom.

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

and Ryan, 1971; Foucault and Mlires, 2000). Aeolian inux from the Sahara and riverine supply of small Near East rivers and wadis are only of minor importance (Stanley et al., 1997; Sandler and Herut, 2000). The enhanced glacial sedimentation rates observed in core SL110 may be a result of stronger currents and/or a shift of the currents and thus of the high-accumulation areas to a more seaward position due to a lower sea level. These processes may also be responsible for the coarsening of the sediments as indicated by the terrigenous sand content and the silt / clay ratios (Fig. 3a).

The increasing glacial sedimentation rates theoretically also could be explained by increasing aridity resulting in an enhanced inux of aeolian dust. Endmember modelling of the terrigenous silt fraction of the nearby core SL112 (892 m water depth), which is located closer to the coast, indicates increased dust uxes during the Last Glacial Maximum (LGM) and the late glacial (Hamann et al., 2008). However, such an inux is not documented in the clay mineral record of SL110, e.g. by the concentration of windblown palygorskite (Fig. 3b). In addition, glacial sedimentation rates at site SL112 were lower when compared to the Holocene, thus supporting the interpretation of a glacialinterglacial shift of the high-accumulation zone.

6.2 Clay mineral assemblages

Smectite in southeastern Levantine Sea sediments is generally regarded to be of detrital origin and is the typical clay mineral in Nile-derived sediments (Venkatarathnam and Ryan, 1971; Maldonado and Stanley, 1981; Stanley and Wingerath, 1996a; Foucault and Mlires, 2000). According to the compilation by Hamann et al. (2009), the Nile assemblage contains up to 85 % smectite. Kaolinite is the second important clay mineral with concentrations in the order of 20 %. Illite contents are < 10 %, and chlorite occurs in trace amounts.

The riverine suspension load is reected in the modern clay mineral assemblage in seaoor surface sediments of the SE Levantine Sea. Accordingly, the clay mineral assemblage off southern Israel consists of up to > 70 % smectite, 10 25 % kaolinite, 10 % illite, and < 10 % chlorite (Hamann et al., 2009) and thus reects the Nile source. The uppermost samples of core SL110 with ca. 60 % smectite, 22 % kaolinite, 8 % illite, and 9 % chlorite t well into this distribution pattern.

The composition of the aeolian clay mineral assemblage reaching the SE Levantine Sea is variable but generally characterized by high illite and kaolinite concentrations (Chester et al., 1977). Palygorskite is a typical clay mineral in aeolian dust derived from the Sahara (Foucault and Mlires, 2000;Scheuvens et al., 2013).

The Nile incorporates almost its entire suspension load close to its source (Adamson et al., 1980). The smectite comes with the Nile suspension load from the catchment of the Blue Nile and Atbara rivers and originates from weath-

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 719

The lightness data (L , Fig. 3a) reect subtle colour changes within the core. The darkest intervals correlate with the sapropels (L ca. 37). However, large parts of the other muds are also relatively dark (L ca. 4045). The lightest interval occurs during the MIS 6 interval at 132135 ka.

TOC contents generally uctuate between 0.4 and 0.7 %, but increase to 2.53 % in the sapropel layers S1, S3, and S4, and to 1.5 % in the pre-sapropelic layer of S5.

The terrigenous sand content is generally < 4 %, and the silt / clay ratio is < 1.0. The data sets show the same general pattern (Fig. 3a) with ner-grained intervals enclosing the sapropel layers. A distinct coarse interval occurs at 132 135 ka. Between 73 and 18 ka we observe an upward coarsening trend. The high silt / clay ratio at ca. 17 ka reects a thin, light grey, not graded silt layer consisting almost entirely of quartz grains.

The clay mineral assemblage of SL110 (Fig. 3b) is dominated by smectite that uctuates between ca. 30 and 65 %.Maximum values occur around 3, 3137, 46, 5059, 65 69, 7789, 99116, and 127132 ka. Minimum values occur at 132135 ka. Illite shows an opposite distribution pattern with concentrations ranging between 5 and 25 %. Kaolinite concentration uctuates between 20 and 30 % and shows only a few distinct maxima centred at 95, 116, and 133 ka and minima at 05, 35, 87, 111, and 128132 ka. Chlorite is present in low concentrations of 711 %, with a pattern roughly resembling that of kaolinite. Palygorskite occurs in trace amounts throughout the core but jumps to 4 % at 132 135 ka.

6 Discussion

6.1 Sedimentation at SL110

According to our age model (Fig. 2, Table 2), sediment core SL110 has an extrapolated basal date of ca. 137 ka. The linear sedimentation rates are 2.2 to 9.8 cm kyr1 (Fig. 3a). They are low during MIS 5 and distinctly higher in late MIS 6 and

MIS 4 to MIS 1. Sedimentation rates for nearby core MD84-641, extracted from data by Fontugne and Calvert (1992), show less distinct changes (Fig. 2) and are distinctly lower than in SL110 due to its position further offshore.

The non-carbonate sediment fraction of seaoor surface sediments in the region of SL110 is dominantly derived by long-distance sediment transport from the Nile (Venkatarathnam and Ryan, 1971; Stanley et al., 1997, 1998; Krom et al., 1999b; Weldeab et al., 2002; Hamann et al., 2009). Storm waves and coastal currents erode coarse sediments from the delta. Geostrophic and wave-induced nearshore and shelf currents displace this material in a large anticlockwise current along the coast towards the Levantine Sea (Fig. 1) and deposit it in the shallow waters off Israel. The ne fraction, in contrast, is derived from direct input from the Nile into the sea and may be transported also further offshore and further to the north (e.g. Stanley et al., 1997, 1998; Venkatarathnam

720 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

ering of Cenozoic volcanic rocks on the Ethiopian Plateau (e.g. Stanley and Wingerath, 1996a). The White Nile south of the conuence of the Sobat River, also discharging from the Ethiopian Highlands, shows a considerably different clay mineral signature with signicantly lower smectite contents.Alluvial sediments from central Uganda have smectite contents < 20 % (Nyakairu and Koeberl, 2001); no smectite was described from lake-bottom sediments of Lake Victoria (Mothersill, 1976). Tropical weathering in the equatorial region produces only small amounts of smectite (Garzanti et al., 2015). Furthermore, the smectite content in SL110 correlates with the Fe content documented by Revel et al. (2010) in delta sediments of the Nile, especially during humid periods. Because Fe is unequivocally derived from the volcanic rocks of the Ethiopian Highlands, this also suggests that the overwhelming part of the smectite comes from that source.

Sediments of the Egyptian wadis discharging into the Nile have clay mineral assemblages with typically 4055 % kaolinite, 3545 % smectite, and 515 % illite (Stanley and Wingerath, 1996b; Hamann et al., 2009; chlorite not determined; thus, actual kaolinite concentrations may be somewhat lower). The kaolinite is mainly derived from the erosion of kaolinite-rich Eocene, Palaeocene, and Mesozoic sediments and lateritic soils in the wadis (Stanley and Wingerath, 1996b; Bolle et al., 2000). Theoretically, the White Nile could also deliver kaolinite during certain times of enhanced sediment transport during humid periods (Garzanti et al. 2015). This would result in a decrease in smectite at the same time due to dilution. This, however, is not the case (Fig. 3b).Therefore, we argue that the kaolinite accumulating in SL110 predominantly comes from the northern African wadi regions.

According to the standard method, percentages of the clay minerals smectite, illite, chlorite, kaolinite, and palygorskite add up to 100 %. A change in the abundance of one mineral therefore also causes changes in the concentration of the others. Hence, we prefer to discuss clay mineral ratios. We use Smr = smectite / (illite + chlorite) as the main proxy for Nile

discharge provided by the Blue Nile and the Atbara River, and Kar = kaolinite / (illite + chlorite) as an indicator for the

contribution of the wadis.

The good correlation between the Kar and the Smr data (Fig. 4; r2 = 0.73) conrms that the kaolinite in SL110 is

of uvial origin (Venkatarathnam and Ryan, 1971; Stanley and Wingerath, 1996a; Hamann et al., 2009) and not derived by aeolian inux from the Sahara, as in other regions of the EMS (e.g. Chester et al., 1977; Foucault and Mlires, 2000;Ehrmann et al., 2013). The correlation implies that during times of enhanced Nile runoff and smectite transport from Ethiopia, the erosion of kaolinite-rich sediments and soils in the wadis was generally also active and that during dry periods in Ethiopia the wadis were dry, too.

This conrms earlier ndings based on the interpretation of satellite images and sedimentological investigations, which showed that extensive drainage systems were active

in the eastern Sahara during humid phases, especially of Eemian and Holocene time (e.g. Pachur and Hoelzmann, 2000; Rohling et al., 2002b; Osborne et al., 2008). Furthermore, oxygen isotope data from speleothems of Soreq Cave (Bar-Matthews et al., 2000) imply that rainfall extended beyond the Sahara.

The clay mineral data of SL110 document a highly variable smectite and kaolinite input through time (Figs. 3b, 4) and show similar distribution patterns as the Fe record from the Nile Delta (Revel et al., 2010; Caley et al., 2011). Fe is supposed to be derived from weathering of volcanic rocks in the Atbara and Blue Nile headwaters in the Ethiopian Highlands, and thus changing Fe contents may reect changing precipitation regimes. Furthermore, the smectite and kaolinite distribution patterns are inversely correlated to the record of Saharan dust in the central Aegean Sea core GeoT SL143 (Ehrmann et al., 2013).

6.3 Penultimate glacial period

The oldest part of the record displays minimum Smr and Kar values at ca. 135132 ka of MIS 6 (Fig. 4). The lack of sediment derived from the Nile Delta indicates that the site was not within the range of the Nile deep-sea turbidite system postulated for MIS 6 at ca. 33 30[prime] to 33 50[prime] E (Ducassou et al. 2009). The clay mineral data point to extremely dry conditions for the penultimate glacial maximum. The extreme drought possibly correlates with Heinrich event 11, which is known as a very cold period in the western Mediterranean Sea and dated to 135130 ka (Martrat et al., 2014; Marino et al., 2015). Obviously, only minor amounts of Nile suspension load reached the coring site off Israel. Instead, maximum palygorskite concentrations (Fig. 3b) and high proportions of terrigenous silt and sand (Fig. 3a) document an increased aeolian inux from the Sahara. This arid phase seems to have been more severe than that of the Last Glacial Maximum.After the maximum drought of MIS 6 the record of SL110 shows a sharp transition to a sapropelic event embedded in an AHP.

6.4 Sapropel formation

It is well known that sapropel formation was coupled to precession-related maxima in Northern Hemisphere summer insolation causing a northward shift of the ITCZ, an intensication of the African monsoon system, and enhanced precipitation in the Ethiopian Highlands (e.g. Rossignol-Strick, 1983; Rohling, 1994; Emeis et al., 2003; Rohling et al., 2015). In SL110, high Smr, Kar, and low silt / clay ratios occur mainly during interglacial periods with negative 18O

values and indicate phases with enhanced Nile and wadi discharge of ne-grained suspension to the Levantine Sea accompanying formation of sapropels S5, S4, S3, and S1.The enhanced rainfall under a warm climate caused stronger chemical weathering of the volcanic rocks and thus smec-

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

remote site is governed dominantly by northern climatic factors (Tzedakis, 2005; Pross et al., 2015).

6.6 The last glacial period

After a further short dry period at 7770 ka, Smr, Kar, and silt / clay ratios indicate AHP 2 and enhanced Nile runoff in early MIS 4 at 6965 ka. AHP 2 is of similar magnitude as the interglacial AHPs discussed before and also starts during increasing insolation (Fig. 4). In contrast to the previous AHPs, this one is not accompanied by low 18O values. It is also not accompanied by a sapropel, dark sediment colour, or enhanced TOC values (Fig. 3a). However, it coincides with the common occurrence of shallow infaunal Uvigerina peregrina. This highly opportunistic species is adapted to the seasonal deposition of high amounts of fresh phytodetritus (Koho et al., 2008). In analogy to nearby SL112, the occurrence of this species at site SL110 during AHP 2 is likely fuelled by pronounced algal blooms related to seasonally enhanced Nile runoff and nutrient supply (Schmiedl et al., 2010). However, in core SL110 we have strong indications that AHP 2 did not culminate in a sapropel formation.Sediment and nutrient delivery by the Nile was cut-off by a drought linked to Heinrich event 6 occurring just at the time of maximum insolation (Fig. 4). This drought, possibly combined with stronger winds and better ventilation of the deep waters, reduced sediment discharge into the Mediterranean Sea and hampered stagnation of the deep water and sapropel formation.

Lourens et al. (1996) described a thin sapropel layer S2, corresponding to AHP 2, south of Crete, that occurs after the 60 ka insolation maximum. It was argued that this sapropel is not visible any more in most sedimentary sequences of the eastern Mediterranean Sea because of post-depositional oxidation of the organic matter (Emeis et al., 2000; Lwemark et al., 2006).

It is not understood why our clay mineral signal has no counterpart in the Fe signal from the Nile Delta (Revel et al., 2010). Possibly Nile discharge to the eastern Mediterranean Sea at that time was not via the western Rosetta branch, but via the eastern Damietta branch, and thus site MS27PT was starved of sediment. It also cannot be ruled out that a match is obscured because of different age models, insofar as the Fe 73 ka peak corresponds to our well-dated 67 ka clay mineral peak (Fig. 4).

The slightly enhanced Smr and Kar values at 5950 ka possibly document the late phase of AHP 2 after the termination of Heinrich event 6, reecting the re-occurrence of slightly enhanced Nile River runoff. Other indications for somewhat increased temperature, humidity, and Nile inux during this time interval come from the Fe record in the Nile Delta (Revel et al., 2010) and the 18O records of Soreq Cave and of sediment cores 9509 and 9501 in the eastern Levantine Sea (Almogi-Labin et al., 2009). The latter records show a marked drop in planktonic foraminiferal 18O values be-

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 721

tite formation. The enhanced river runoff contributed to the stratication of the water column and facilitated the transport of nutrients and weathering products to the EMS promoting a shoaling of the nutricline and thus development of a deep chlorophyll maximum (Rohling and Gieskes, 1989;Grelaud et al., 2012) particularly in areas close to the Nile Delta. These processes fostered sustained anoxic conditions in the deeper basins and preservation of TOC (Fig. 3a).

Both the clay mineral record of southeastern Levantine Sea core SL110, the Fe record from the Nile Delta (Revel et al., 2010; Caley et al., 2011), and the record of Saharan dust inux to the Aegean Sea (Ehrmann et al., 2013) reveal that each phase of sapropel formation occurred within a longer AHP (Fig. 4; AHPs numbered according to the corresponding sapropels).

6.5 The last interglacial period

The Smr maximum linked to AHP 5 shows a sharp onset after the glacial maximum of MIS 6, simultaneously with the start of the pre-sapropelic layer of S5 (Fig. 4). It indicates a northward-moving rain belt resulting in enhanced monsoonal rainfall in Ethiopia and Blue Nile runoff as early as 132 ka at a time of increasing insolation but under the inuence of penultimate glacial boundary conditions. Also the Kar values clearly increase with the beginning of the pre-sapropelic layer. However, they do not exhibit such distinct maximum values as during the early phases of the younger AHPs. Most of sapropel S5 and the termination of AHP 5, however, is not preserved in SL110 due to the hiatus at ca. 126119 ka. Furthermore, because of the curved top of the sapropel and the technique of sampling the core in parallel slices, we did not receive the typical sapropel signature but a mixed sapropel post-sapropel signal (Fig. 4).

According to the Fe and the dust record (Revel et al., 2010; Ehrmann et al., 2013), AHP 4 started at > 108 ka and lasted until 95 or 98 ka, respectively. The Smr and Kar data of SL110 indicate that AHP 4 started during insolation rise at 116 ka and ended at 99 ka, with sapropel formation occurring at 108104.5 ka (Fig. 4).

AHP 3 started just before the insolation maximum and lasted from ca. 89 to 77 ka based on the Smr and Kar maxima and the minima in the silt / clay ratio and the terrigenous sand content. These ages are in good accordance with those obtained from the dust record in the Aegean Sea (Ehrmann et al., 2013). Sapropel S3 has a shorter duration and spans the time interval 8681 ka.

Despite the hiatus in the SL110 record, by assuming a duration spanning 128121 ka for S5 (Grant et al., 2012), our data indicate that the dry phase between AHP 5 and AHP 4 lasted only ca. 5 kyr, compared with ca. 11 kyr for the one between AHP 4 and AHP 3. It is interesting to note that a similar duration of the dry phases can be deduced from the pollen record of Tenaghi Philippon in Greece, although this

722 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

tween 58 and 49 ka BP. The signal of warm and freshened surface waters was expressed particularly clearly in core 9501 from the northern part of the basin and coincided with DO interstadial 14. The spread of deciduous and other trees during the same time period likely supported the migration of anatomically modern humans from Africa into the Levant (Mller et al., 2011). In the Smr data of core SL110, this time interval documents only a moderate increase in Nile runoff (Fig. 4), suggesting that the observed warming and freshening of the eastern Mediterranean region was predominantly related to a high-latitude climate change.

Thus, the sediments of SL110 document that the long-term climate regime governed by the African monsoon and characterized by pronounced AHPs during the last inter-glacial period was replaced at the occurrence of Heinrich event 6 by a distinctly different glacial climate system governed by short-term, millennial-scale climate changes having their origin in the North Atlantic and Arctic regions. Former studies showed that sudden freshwater inuxes into the northern Atlantic Ocean during the Younger Dryas and the Heinrich events not only weakened the thermohaline circulation and oceanic heat transport to the northern high latitudes (e.g. Bond et al., 1993) but also had major consequences for the Mediterranean Sea environment. Increased aridity and even northern African mega-droughts during such events can be seen in marine sediments of the North Atlantic Ocean (Mulitza et al., 2008), speleothem data, and lake level drops of Lake Lisan (palaeo-Dead Sea; Bar-Matthews et al., 1999; Bartov et al., 2003) and of African lakes (Gasse et al., 2008; Shanahan et al., 2015). Cold water inux from the North Atlantic Ocean to the Mediterranean Sea intensied.Sea surface temperatures and deep-water temperatures decreased and deep-water ventilation strengthened in the northwestern Mediterranean Sea (Cacho et al., 1999, 2000) and probably also in the EMS. Enhanced northward ux of wind-transported sediment from the Sahara is documented for both the western and the eastern Mediterranean Sea (Moreno et al., 2002; Bout-Roumazeilles et al., 2007; Hamann et al., 2008). Strongly decreased Smr and Kar values in core SL110 (Fig. 4) indicate reduced Nile sediment discharge during Heinrich events and thus major droughts in the headwaters of the Atbara River and Blue Nile and diminished erosion and river runoff. During these events the ITCZ and the African rain belt moved southward. Both Lake Tana, the source of the Blue Nile, and Lake Victoria, the source of the White Nile, were desiccated during the Younger Dryas and Heinrich event 1 (Lamb et al., 2007; Marshall et al., 2011). Minimum Smr and Kar values as well as the occurrence of a silt layer of probably aeolian origin (Fig. 3a) imply that Heinrich event 1 was the driest interval since the MIS 6 glacial maximum (cf.Stager et al., 2011).

The 65 N insolation curve shows a further low-amplitude maximum at ca. 37 ka. It is also reected by slightly enhanced smectite and kaolinite abundances in core SL110, ranging from 37 to 31 ka, and by the Fe record in the

Nile Delta (Revel et al., 2010). No sapropel formation was recorded following the weak insolation maximum at 37 ka.Thus, the two potential humid periods linked to the glacial insolation maxima at 60 and 37 ka are only vaguely expressed in the EMS because their effects were much less distinct and because they were suppressed by the glacial climate regime, particularly the impacts of cold and dry stadials.

6.7 The last African humid period and the Holocene

The clay mineral data of SL110 and of nearby core SL112 (Hamann et al., 2009) document a gradual onset of AHP 1 at about 14 ka (Fig. 4), which coincides with the overow of Lake Tana, the source of the Blue Nile (Lamb et al., 2007;Marshall et al., 2011). At the same time, Saharan dust inux to the EMS decreased (Ehrmann et al., 2013). S1 at ca. 116.5 ka and the peak Blue Nile runoff at 109 ka (Weldeab et al., 2014) are only poorly reected in the clay mineral data.Precipitation generally declined gradually after sapropel formation throughout the Holocene, the Nile River runoff responded quasilinearly to the changes in rainfall (Weldeab et al., 2014; Blanchet et al., 2014), and the AHP 1 ended at ca. 5 ka (deMenocal et al., 2000). However, the clay mineral data do not mirror this trend but show continued increasing Smr and Kar values and Nile sediment discharge until ca. 2 ka. Sediment discharge is thus decoupled from water runoff. Possibly, human activity altered the natural signal.Saharan population increased shortly after 11 ka and reached a maximum at 7.5 ka (Manning and Timpson, 2014). Major land-use changes occurred with the advent of Neolithic farming some 8 kyr ago (Williams, 2009).

6.8 Feedback effects

The clay mineral and grain size data of SL110 indicate an enhanced riverine inux of suspension from the Blue Nile and the Atbara River to the Mediterranean Sea during the AHPs and therefore document enhanced weathering of the Ethiopian basalts, erosion of soils, and transport of the suspension load to the sea. This view is supported by several investigations (Wehausen and Brumsack, 1998; Emeis et al., 2000; Revel et al., 2010; Blanchet et al., 2014). However, it is in conict with some other studies (Adamson et al., 1980;Krom et al., 1999a, 2002; Box et al., 2011) that have argued for a lower Blue Nile sediment discharge during the warm and humid Holocene AHP and higher uxes during dry climate phases. They reasoned that at present the limited vegetation cover on the Ethiopian Highlands allows extensive erosion when the summer monsoon rains fall. During increased monsoon activity, in contrast, longer periods of rain, longer growing seasons, and larger vegetation-covered areas reduced soil erosion.

Our data indicate that an increase in precipitation and Nile suspension delivery happened long before sapropel formation, during phases of increasing insolation. However, in

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 723

the case of AHP 4, Smr and Kar maxima occur just before and after sapropel S4 formation, each accompanied by low silt / clay ratios, also indicating enhanced riverine suspension (Fig. 3a). Thus, maximum riverine sediment inux does not coincide with maximum water discharge during sapropel formation. This is possibly due to a feedback mechanism as postulated, for example, by Adamson et al. (1980) and Krom et al. (2002). Vegetation feedbacks were obviously less expressed during glacial conditions, particularly during the last glacial maximum. This is linked to the generally more arid boundary conditions during the last glacial in tropical and northern Africa (Gasse, 2000) and resulting sparse vegetation cover, particularly in the deation areas of the Sahara (Tjallingii et al., 2008). The rapid climate changes linked to DansgaardOeschger cycles and Heinrich events nd an immediate response in the sediment delivery to the EMS.

The different timing between the Saharan dust record and the Nile discharge record of the humid phases (Fig. 4) may reect regional differences in climatic history. The northward-migrating tropical rain belt may rst have reached the Ethiopian Highlands, where the smectite signals has its origin, and with a time lag the central Sahara, there inhibiting effective dust production. The differences may also be due to different processes being at work. The Nile discharge is mainly driven by annual oods, uvial erosion, and weathering. It therefore may have reacted immediately on a northward shift of the ITCZ and the rain belt. The dust record, in contrast, is mainly a function of aridity, wind strength, and, probably most important, the availability of deatable material. The spreading and retreat of vegetation cover plays an important role and may react non-linearly and with a different time lag to a shift of the ITCZ than river discharge (Claussen et al., 2013; Ehrmann et al., 2013). Thus, recently Blanchet et al. (2014) could show for the end of AHP 1 that the vegetation retreat from the Sahara was much faster than the decrease in Nile runoff during the southward migration of the rain belt.

7 Conclusions

Clay minerals are suitable tools to reconstruct the late Quaternary Nile sediment discharge to the EMS. The abundance of smectite reects the runoff of the Blue Nile and the Atbara River, whereas kaolinite abundance reects the discharge from wadis. The correlation of smectite and kaolinite abundances implies that during times of enhanced Nile suspension discharge the wadis were also active sediment sources.

High smectite and kaolinite abundances and low terrigenous silt / clay ratios indicate phases of enhanced precipitation and suspension inux from the Blue Nile, the Atbara River, and the wadis to the Levantine Sea, the so-called African humid periods (AHPs), which are accompanied by the formation of sapropel layers. This largely contradicts the view that sediment delivery during humid periods was re-

duced because an extensive vegetation cover hampered soil

erosion.

Precipitation and suspension discharge started much earlier and ended later than sapropel formation. Maximum sediment discharge does not coincide with maximum water discharge during times of sapropel formation because of minor vegetation feedback.

The time lag between the increase and decrease in Nile suspension discharge and the decrease and increase in the dust record (SL143) can be explained by a successive migration of the rain belt, possibly combined with a vegetation feedback.

The dry periods between the interglacial AHPs were relatively short, each lasting only ca. 510 kyr.

The sedimentary record of core SL110 demonstrates the interplay of two competing climate systems. The long-term climate regime was governed by the African monsoon and characterized by pronounced AHPs during the last inter-glacial period. This system was replaced in the last glacial by short-term, millennial-scale climate changes having their origin in the North Atlantic region. They are characterized by short but pronounced drought periods linked to Heinrich events and more humid interstadials. Precipitation in the Nile catchment, river runoff, and sapropel formation can be hampered or even inhibited under the inuence of the North Atlantic climate system, particularly under glacial boundary conditions.

Data availability

The raw data for this article are stored in the PANGAEA database at the Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany (http://doi.pangaea.de/10.1594/PANGAEA.848291

Web End =http://doi.pangaea.de/10. http://doi.pangaea.de/10.1594/PANGAEA.848291

Web End =1594/PANGAEA.848291 ).

Acknowledgements. We thank the chief scientist C. Hemleben, the master, and the crew of the RV Meteor for their support during cruise 51/3. For technical assistance in the laboratory, we thank S. Haener and M. Haener. C. Kertscher performed TOC analyses. M. Paul is acknowledged for picking foraminifera and establishing a preliminary age model. We acknowledge M. Revel and an anonymous reviewer for their comments, which helped to improve our manuscript. The Deutsche Forschungsgemeinschaft (DFG) nancially supported this study.

Edited by: L. Beaufort

References

Adamson, D. A., Gasse, F., Street, F. A., and Williams, M. A. J.:

Late Quaternary history of the Nile, Nature, 288, 5055, 1980.

Almogi-Labin, A., Bar-Matthews, M., Shriki, D., Kolosovsky, E.,

Paterne, M., Schilman, B., Ayalon, A., Aizenshtat, Z., and Matthews, A.: Climatic variability during the last 90 ka of the

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

724 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

southern and northern Levantine Basin as evident from marine records and speleothems, Quaternary Sci. Rev., 28, 28822896, 2009.

Bar-Matthews, M., Ayalon, A., Kaufman, A., and Wasserburg, G.J.: The Eastern Mediterranean paleoclimate as a reection of regional events: Soreq cave, Israel, Earth Planet. Sci. Lett., 166, 8595, 1999.

Bar-Matthews, M., Ayalon, A., and Kaufman, A.: Timing and hydrological conditions of sapropel events in the Eastern Mediterranean, as evident from speleothems, Soreq cave, Israel, Chem.Geol., 169, 145156, 2000.

Bartov, Y., Goldstein, S. L., Stein, M., and Enzel, Y.: Catastrophic arid episodes in the Eastern Mediterranean linked with the North Atlantic Heinrich events, Geology, 31, 439442, 2003.Biscaye, P. E.: Distinction between kaolinite and chlorite in recent sediments by X-ray diffraction, Am. Mineral., 49, 12811289, 1964.

Biscaye, P. E.: Mineralogy and sedimentation of recent deep-sea clay in the Atlantic Ocean and adjacent seas and oceans, Geol.Soc. Am. Bull., 76, 803832, 1965.

Blanchet, C. L., Frank, M., and Schouten, S.: Asynchronous changes in vegetation, runoff and erosion in the Nile River watershed during the Holocene, PLoS One, 9, e115958, doi:http://dx.doi.org/10.1371/journal.pone.0115958

Web End =10.1371/journal.pone.0115958 http://dx.doi.org/10.1371/journal.pone.0115958

Web End = , 2015.

Bolle, M.-P., Tantawy, A. A., Pardo, A., Adatte, T., Burns, S., and Kassab, A.: Climatic and environmental changes documented in the upper Paleocene and lower Eocene of Egypt, Eclogae Geol.Helv., 93, 3351, 2000.

Bond, G., Broecker, W., Johnsen, S., McManus, J., Labeyrie, L., Jouzel, J., and Bonani, G.: Correlation between climate records from North Atlantic sediments and Greenland ice, Nature, 365, 143147, 1993.

Bout-Roumazeilles, V., Combourieu Nebout, N., Peyron, O., Cortijo, E., Landais, A., and Masson-Delmotte, V.: Connection betweeen South Mediterranean climate and North African atmospheric circulation during the last 50,000 yr BP North Atlantic cold events, Quaternary Sci. Rev., 26, 31973215, 2007.

Box, M. R., Krom, M. D., Cliff, R. A., Bar-Matthews, M., Almogi-Labin, A., Ayalon, A., and Paterne, M.: Response of the Nile and its catchment to millenial-scale climatic change since the LGM from Sr isotopes and major elements of Eastern Mediterranean sediments, Quaternary Sci. Rev., 30, 431442, 2011.

Cacho, I., Grimalt, J. O., Pelejero, C., Canals, M., Sierro, F. J., Flores, J. A., and Shackleton, N.: Dansgaard-Oeschger and Heinrich event imprints in Alboran Sea paleotemperatures, Paleoceanography, 14, 698705, 1999.

Cacho, I., Grimalt, J. O., Sierro, F. J., Shackleton, N., and Canals,M.: Evidence for enhanced Mediterranean thermohaline circulation during rapid climatic coolings, Earth Planet. Sci. Lett., 183, 417429, 2000.

Caley, T., Malaiz, B., Revel, M., Ducassou, E., Wainer, K., Ibrahim, M., Shoeaib, D., Migeon, S., and Marieu, V.: Orbital timing of the Indian, East Asian and African boreal monsoons and the concept of a global monsoon, Quaternary Sci. Rev., 30, 37053715, 2011.

Chester, R., Baxter, G. G., Behairy, A. K. A., Connor, K., Cross,D., Eldereld, H., and Padgham, R. C.: Soil-sized eolian dusts from the lower troposphere of the eastern Mediterranean Sea, Mar. Geol., 24, 201217, 1977.

Claussen, M., Bathiany, S., Brovkin, V., and Kleinen, T.: Simulated climate-vegetation interaction in semi-arid regions affected by plant diversity, Nat. Geosci., 6, 954958, 2013.

deMenocal, P., Ortiz, J., Guilderson, T., Adkins, J., Sarnthein, M.,

Baker, L., and Yarusinsky, M.: Abrupt onset and termination of the African Humid Period: rapid climate responses to gradual insolation forcing, Quaternary Sci. Rev., 19, 347361, 2000.Ducassou, E., Migeon, S., Capotondi, L., and Mascle, J.: Run-out distance and erosion of debris-ows in the Nile deep-sea fan system: Evidence from lithofacies and micropalaeontological analyses, Mar. Petrol. Geol., 39, 102123, 2013.

Ehrmann, W., Schmiedl, G., Hamann, Y., Kuhnt, T., Hemleben, C., and Siebel, W.: Clay minerals in late glacial and Holocene sediments of the northern and southern Aegean Sea, Palaeogeogr.Palaeocl., 249, 3657, 2007.

Ehrmann, W., Seidel, M., and Schmiedl, G.: Dynamics of Late Quaternary North African humid periods documented in the clay mineral record of central Aegean Sea sediments, Glob. Planet.Change, 107, 186195, 2013.

Emeis, K.-C., Sakamoto, T., Wehausen, R., and Brumsack, H.-J.:

The sapropel record of the eastern Mediterranean Sea results of Ocean Drilling Program Leg 160, Palaeogeogr. Palaeocl., 158, 371395, 2000.

Emeis, K.-C., Schulz, H., Struck, U., Rossignol-Strick, M., Erlenkeuser, H., Howell, M. W., Kroon, D., Mackensen, A., Ishizuka, S., Oba, T., Sakamoto, T., and Koizumi, I.: Eastern Mediterranean surface water temperatures and 18O composition during deposition of sapropels in the late Quaternary, Paleoceanography, 18, 1005, doi:http://dx.doi.org/10.1029/2000PA000617

Web End =10.1029/2000PA000617 http://dx.doi.org/10.1029/2000PA000617

Web End = , 2003.Fontugne, M. R. and Calvert, S. E.: Late Pleistocene variability of the carbon isotopic composition of organic matter in the Eastern Mediterranean: Monitor of changes in carbon sources and atmospheric CO2 concentrations, Paleoceanography, 7, 120, 1992.

Foucault, A. and Stanley, D. J.: Late Quaternary paleoclimatic oscillations in East Africa recorded by heavy minerals in the Nile delta, Nature, 3390, 4446, 1989.

Foucault, A. and Mlires, F.: Palaeoclimatic cyclicity in central Mediterranean Pliocene sediments: the mineralogical signal, Palaeogeogr. Palaeocl., 158, 311323, 2000.

Garzanti, E., Ando, S., Padoan, M., Vezzoli, G., and El Kammar,A.: The modern Nile sediment system: Processes and products, Quaternary Sci. Rev., 130, 956, 2015.

Gasse, F.: Hydrological changes in the African tropics since the Last

Glacial Maximum, Quaternary Sci. Rev., 19, 189211, 2000.

Gasse, F., Chalie, F., Vincens, A., Williams, M. A. J., and

Williamson, D.: Climatic patterns in equatorial and southern Africa from 30,000 to 10,000 years ago reconstructed from terrestrial and near-shore proxy data, Quaternary Sci. Rev., 27, 23162340, 2008.

Grant, K. M., Rohling, E. J., Bar-Matthews, M., Ayalon, A.,

Medina-Elizalde, M., Bronk Ramsey, C., Satow, C., and Roberts,A. P.: Rapid coupling between ice volume and polar temperature over the past 150,000 years, Nature, 491, 744747, 2012.Grelaud, M., Marino, G., Ziveri, P., Rohling, E. J.: Abrupt shoaling of the nutricline in response to massive freshwater ooding at the onset of the last interglacial sapropel event, Paleoceanography, 27, PA3208, doi:3210.1029/2012PA002288, 2012.

Hamann, Y., Ehrmann, W., Schmiedl, G., Krger, S., Stuut, J.-B., and Kuhnt, T.: Sedimentation processes in the Eastern Mediter-

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 725

ranean Sea during the Late Glacial and Holocene revealed by end-member modelling of the terrigenous fraction in marine sediments, Mar. Geol., 248, 97114, 2008.

Hamann, Y., Ehrmann, W., Schmiedl, G., and Kuhnt, T.: Modern and late Quaternary clay mineral distribution in the area of the SE Mediterranean Sea, Quaternary Res., 71, 453464, 2009.Hemleben, C., Hoernle, K., Jrgensen, B. B., and Roether, W.

(Eds.): OstatlantikMittelmeerSchwarzes Meer, Cruise No. 51,12 September28 December 2001, Meteor-Berichte, Universitt Hamburg, 03-01, 225 pp., 2003.

Holeman, J. N.: The sediment yield of major rivers of the world,Water Resour. Res., 4, 737747, 1968.

Kallel, N., Duplessy, J.-C., Labeyrie, L., Fontugne, M., Paterne, M., and Montacer, M.: Mediterranean pluvial periods and sapropel formation over the last 200 000 years, Palaeogeogr. Palaeocl., 157, 4558, 2000.

Koho, K. A., Garca, R., de Stigter, H. C., Epping, E., Koning, E., Kouwenhoven, T. J., and Van der Zwaan, G. J.: Sedimentary labile organic carbon and pore water redox control on species distribution of benthic foraminifera: A case study from Lisbon-Setbal Canyon (southern Portugal), Prog. Oceanogr., 79, 5582, 2008.

Krom, M. D., Michard, A., Cliff, R. A., and Strohle, K.: Sources of sediment to the Ionian Sea and western Levantine basin of the Eastern Mediterranean during S-1 sapropel times, Mar. Geol., 160, 4561, 1999a.

Krom, M. D., Cliff, R. A., Eijsink, L. M., Herut, B., and Chester, R.: The characterisation of Saharan dusts and Nile particulate matter in surface sediments from the Levantine basin using Sr isotopes, Mar. Geol., 155, 319330, 1999b.

Krom, M. D., Stanley, J. D., Cliff, R. A., and Woodward, J. C.: Nile

River sediment uctuations over the past 7000 yr an their key role in sapropel development, Geology, 30, 7174, 2002.Krpelin, S., Verschuren, D., Lzine, A.-M., Eggermont, H., Cocquyt, C., Francus, P., Cazet, J.-P., Fagot, M., Rumes, B., Russell,J. M., Darius, F., Conley, D. J., Schuster, M., von Suchodoletz,H., and Engstrom, D. R.: Climate-driven ecosystem succession in the Sahara: The past 6000 years, Science, 320, 765768, 2008.Lamb, H. F., Bates, C. R., Coombes, P. V., Marshall, M. H., Umer,M., Davies, S. J., and Dejen, E.: Late Pleistocene desiccation of Lake Tana, source of the Blue Nile, Quaternary Sci. Rev., 26, 287299, 2007.

Lisiecki, L. E. and Raymo, M. E.: A Pliocene-Pleistocene stack of 57 globally distributed benthic 18O records, Paleoceanography, 20, PA1003, doi:http://dx.doi.org/10.1029/2004PA001071

Web End =10.1029/2004PA001071 http://dx.doi.org/10.1029/2004PA001071

Web End = , 2005.

Lourens, L. J., Antonarakou, A., Hilgen, F. J., Van Hoof, A. A. M., Vergnaud-Grazzini, C., and Zachariasse, W. J.: Evaluation of the Plio-Pleistocene astronomical timescale, Paleoceanography, 11, 391413, 1996.

Lwemark, L., Lin, Y., Chen, H.-F., Yang, T.-N., Beier, C., Werner,F., Lee, C.-Y., Song, S.-R., and Kao, S.-J.: Sapropel burn-down and ichnological response to late Quaternary sapropel formation in two 400 ky records from the eastern Mediterranean Sea,

Palaeogeogr. Palaeocl., 239, 406425, 2006.

Maldonado, A. and Stanley, D. J.: Clay mineral distribution patterns as inuenced by depositional processes in the southeastern Levantine Sea, Sedimentology, 28, 2132, 1981.

Manning, K. and Timpson, A.: The demographic response to

Holocene climate change in the Sahara, Quaternary Sci. Rev., 101, 2835, 2014.

Marino, G., Rohling, E. J., Rodriguez-Sanz, L., Grant, K. M., Heslop, D., Roberts, A. P., Stanford, J. D., and Yu, J.: Bipolar seesaw control on last interglacial sea level, Nature, 522, 197201, 2015.Marshall, M. H., Lamb, H. F., Huws, D., Davies, S. J., Bates, R.,

Bloemendal, J., Boyle, J., Leng, M. J., Umer, M., and Bryant,C.: Late Pleistocene and Holocene drought events at Lake Tana, the source of the Blue Nile, Glob. Planet. Change, 78, 147161, 2011.

Martrat, B., Jimenez-Amat, P., Zahn, R., and Grimalt, J. O.: Similarities and dissimilarities between the last two deglaciations and interglaciations in the North Atlantic region, Quaternary Sci. Rev., 99, 122134, 2014.

Milliman, J. D. and Syvitski, J. P. M.: Geomorphic/tectonic control of sediment discharge to the ocean: the importance of small mountainous rivers, J. Geol., 100, 525544, 1992.

Moreno, A., Cacho, I., Canals, M., Prins, M. A., Sanchez-Goni, M.-F., Grimalt, J. O., and Weltje, G. J.: Saharan dust transport and high-latitude glacial climatic variability: the Alboran Sea record, Quaternary Res., 58, 318328, 2002.

Mothersill, J. S.: The mineralogy and geochemistry of the sediments of northwestern Lake Victoria, Sedimentology, 23, 553 565, 1976.

Mulitza, S., Prange, M., Stuut, J.-B., Zabel, M., von Dobe-neck, T., Itambi, A. C., Nizou, J., Schulz, M., and Wefer,G.: Sahel megadroughts triggered by glacial slowdowns of Atlantic meridional overturning, Paleoceanography, 23, PA4206, doi:http://dx.doi.org/10.1029/2008PA001637

Web End =10.1029/2008PA001637 http://dx.doi.org/10.1029/2008PA001637

Web End = , 2008.

Mller, U. C., Pross, J., Tzedakis, P. C., Gamble, C., Kotthoff, U.,

Schmiedl, G., Wulf, S., and Christanis, K.: The role of climate in the spread of modern humans into Europe, Quaternary Sci. Rev., 30, 273279, 2011.

Nyakairu, G. W. A. and Koeberl, C.: Mineralogical and chemical composition and distribution of rare earth elements in clay-rich sediments from central Uganda, Geochem. J., 35, 1328, 2001.Osborne, A. H., Vance, D., Rohling, E. J., Barton, N., Rogerson,M., and Fello, N.: A humid corridor across the Sahara for the migration of early modern humans out of Africa 120,000 years ago, PNAS, 105, 1644416447, 2008.

Pachur, H.-J. and Hoelzmann, P.: Late Quaternary palaeoecology and palaeoclimates of the eastern Sahara, J. Afr. Earth Sci., 30, 929939, 2000.

Paillard, D., Labeyrie, L., and Yiou, P.: Macintosh program performs time-series analysis, Eos Transactions, American Geophysical Union, 77, p. 379, 1996.

Petschick, R.: MacDiff 4.2.5. http://www.geologie.uni-frankfurt.de/staff/homepages/petschick/petschick.html

Web End =http://www.geologie.uni-frankfurt.de/ http://www.geologie.uni-frankfurt.de/staff/homepages/petschick/petschick.html

Web End =staff/homepages/petschick/petschick.html (last access: 15 December 2015), 2001.

Pickard, G. L. and Emery, W. J.: Descriptive Physical Oceanography, Pergamon Press, 320 pp., 1990.

Pross, J., Koutsodendris, A., Christanis, K., Fischer, T., Fletcher,W. J., Hardiman, M., Kalaitzidis, S., Knipping, M., Kotthoff,U., Milner, A. M., Mller, U. C., Schmiedl, G., Siavalas, G., Tzedakis, P. C., and Wulf, S.: The 1.35-Ma-long terrestrial climate archive of Tenaghi Philippon, northeastern Greece: Evolution, exploration, and perspectives for future research, Newsl.Stratigr., 48, 253276, 2015.

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016

726 W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge

Quelennec, R. E. and Kruc, C. B.: Nile suspended load and its importance for the Nile delta morphology. UNDP/UNESCO Proceedings of the Seminar on Nile Delta Sedimentation. Alexandria, Academy of Scientic Research and Technology, Egypt, 1 257, 1976.

Reimer, P. J., Bard, E., Bayliss, A., Beck, J. W., Blackwell, P. G., Bronk Ramsey, C., Buck, C. E., Cheng, H., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Haidason, H., Hajdas, I., Hatt, C., Heaton, T. J., Hoffmann, D. L., Hogg, A. G., Hughen, K. A., Kaiser, K. F., Kromer, B., Manning, S. W., Niu,M., Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Staff, R. A., Turney, C. S. M., and van der Plicht, J.: IntCal13 and Marine13 radiocarbon age calibration curves 050,000 years cal BP, Radiocarbon, 55, 18691887, 2013.

Renssen, H., Brovkin, V., Fichefet, T., and Goosse, H.: Simulation of the Holocene climate evolution in Northern Africa: The termination of the African Humid Period, Quatern. Int., 150, 95102, 2006.

Revel, M., Ducassou, E., Grousset, F. E., Bernasconi, S. M., Migeon, S., Revillon, S., Mascle, J., Murat, A., Zaragosi, S., and Bosch, D.: 100,000 Years of African monsoon variability recorded in sediments of the Nile margin, Quaternary Sci. Rev., 29, 13421362, 2010.

Rohling, E.: Review and new aspects concerning the formation of eastern Mediterranean sapropels, Mar. Geol., 122, 128, 1994.Rohling, E. J. and Gieskes, W. W. C.: Late Quaternary changes in

Mediterranean intermediate water density and formation rate, Paleoceanography, 4, 531545, 1989.

Rohling, E. J., Mayewski, P. A., Abu-Zied, R. H., Casford, J. S.L., and Hayes, A.: African monsoon variability during the previous interglacial maximum, Earth Planet. Sci. Lett., 202, 6175, 2002a.

Rohling, E. J., Mayewski, P. A., Abu-Zied, R. H., Casford, J. S. L., and Hayes, A.: Holocene atmosphere-ocean interactions: records from Greenland and the Aegean Sea, Clim. Dynam., 18, 587 593, 2002b.

Rohling, E. J., Marino, G., and Grant, K. M.: Mediterranean climate and oceanography, and the periodic development of anoxic events (sapropels), Earth-Sci. Rev., 143, 6297, 2015.Rossignol-Strick, M.: African monsoons, an immediate climate response to orbital insolation, Nature, 304, 4649, 1983.Rossignol-Strick, M., Nesteroff, W., Olive, P., and Vergnaud-

Grazzini, C.: After the deluge: Mediterranean stagnation and sapropel formation, Nature, 295, 105110, 1982.

Sandler, A. and Herut, B.: Composition of clays along the continental shelf off Israel: contribution of the Nile versus local sources, Mar. Geol., 167, 339354, 2000.

Scheuvens, D., Schtz, L., Kandler, K., Ebert, M., and Weinbruch,S.: Bulk composition of northern African dust and its source sediments - A compilation, Earth-Sci. Rev., 116, 170194, 2013.Schmiedl, G., Mitschele, A., Beck, S., Emeis, K.-C., Hemleben, C.,

Schulz, H., Sperling, M., and Weldeab, S.: Benthic foraminiferal record of ecosystem variability in the eastern Mediterranean Sea during times of sapropel S5 and S6 deposition, Palaeogeogr.Palaeocl., 190, 139164, 2003.

Schmiedl, G., Kuhnt, T., Ehrmann, W., Emeis, K.-C., Hamann, Y., Kotthoff, U., Dulski, P., and Pross, J.: Climatic forcing of eastern Mediterranean deep-water formation and benthic ecosystems

during the past 22000 years, Quaternary Sci. Rev., 29, 3006 3020, 2010.

Shanahan, T. M., McKay, N. P., Hughen, K. A., Overpeck, J. T., Otto-Bliesner, B., Heil, C. W., King, J., Scholz, C. A., and Peck,J.: The time-transgressive termination of the African Humid Period, Nat. Geosci., 8, 140144, 2015.

Stager, J. C., Ryves, D. B., Chase, B. M., and Pausata, F. S. R.: Catastrophic drought in the Afro-Asian monsoon region during Heinrich Event 1, Science, 331, 12991302, 2011.

Stanley, D. J. and Wingerath, J. G.: Clay mineral distributions to interpret Nile Cell provenance and dispersal: I. Lower River Nile to delta sector, J. Coast. Res., 12, 911929, 1996a.

Stanley, D. J. and Wingerath, J. G.: Nile sediment dispersal altered by the Aswan High Dam: The kaolinite trace, Mar. Geol., 133, 19, 1996b.

Stanley, D. J., Mart, Y., and Nir, Y.: Clay mineral distributions to interpret Nile Cell provenance and dispersal: II. Coastal plain from Nile Delta to northern Israel, J. Coast. Res., 13, 506533, 1997. Stanley, D. J., Nir, Y., and Galili, E.: Clay mineral distributions to interpret Nile Cell provenance and dispersal: III. Offshore margin between Nile delta and northern Israel, J. Coast. Res., 14, 196 217, 1998.

Tjallingii, R., Claussen, M., Stuut, J.-B. W., Fohlmeister, J., Jahn,A., Bickert, T., Lamy, F., and Rhl, U.: Coherent high- and low-latitude control of the northwest African hydrological balance, Nat. Geosci., 1, 670675, 2008.

Tzedakis, P. C.: Vegetation variability in Greece during the last interglacial, Geol. Mijnbouw, 79, 355367, 2000. Venkatarathnam, K. and Ryan, W. B. F.: Dispersal patterns of clay minerals in the sediments of the eastern Mediterranean Sea, Mar. Geol., 11, 261282, 1971.

Weaver, C. E. and Beck, K. C.: Miocene of the S.E. United States: A model for chemical sedimentation in a peri-marine environment, Sediment. Geol., 17, 1234, 1977Wehausen, R. and Brumsack, H.-J.: The formation of Pliocene

Mediterranean sapropels: constraints from high-resolution major and minor element studies, in: Proceedings of the Ocean Drilling Program, Scientic Results, edited by: Robertson, A. H.F., Emeis, K.-C., Richter, C., and Camerlenghi, A., 160, Ocean Drilling Program, College Station, TX, 207217, 1998. Wehausen, R. and Brumsack, H. J.: Chemical cycles in Pliocene sapropel-bearing and sapropel-barren eastern Mediterranean sediments, Palaeogeogr. Palaeocl., 158, 325352, 2000.

Weldeab, S., Emeis, K.-C., Hemleben, C., and Siebel, W.: Provenance of lithogenic surface sediments and pathways of riverine suspended mattern in the Eastern Mediterranean Sea: evidence from 143Nd / 144Nd and 87Sr / 86Sr ratios, Chem. Geol., 186, 139149, 2002.

Weldeab, S., Emeis, K.-C., Hemleben, C., Schmiedl, G., and

Schulz, H.: Spatial productivity variations during formation of sapropels S5 and S6 in the Mediterranean Sea: evidence from Ba contents, Palaeogeogr. Palaeocl., 191, 169190, 2003. Weldeab, S., Menke, V., and Schmiedl, G.: The pace of East African monsoon evolution during the Holocene, Geophys. Res. Lett., 41, doi:http://dx.doi.org/10.1002/2014GL059361

Web End =10.1002/2014GL059361 http://dx.doi.org/10.1002/2014GL059361

Web End = , 2014.

Williams, M., Talbot, M., Aharon, P., Salaam, Y. A., Williams, F., and Brendeland, K. I.: Abrupt return of the summer monsoon 15,000 years ago: new supporting evidence from the lower White

Clim. Past, 12, 713727, 2016 www.clim-past.net/12/713/2016/

W. Ehrmann et al.: A distal 140 kyr sediment record of Nile discharge 727

Nile valley and Lake Albert, Quaternary Sci. Rev., 25, 2651 2665, 2006.

Williams, M. A. J.: Late Pleistocene and Holocene environments in the Nile basin, Glob. Planet. Change, 69, 115, 2009.

Williams, M. A. J., Usai, D., Salvatori, S., Williams, F. M., Zerboni,A., Maritan, L., and Linseele, V.: Late Quaternary environments and prehistoric occupation in the lower White Nile valley, central Sudan, Quaternary Sci. Rev., 130, 7288, 2015.

www.clim-past.net/12/713/2016/ Clim. Past, 12, 713727, 2016