Abstract
Laminated sediments of the continental slope off the Makran coast in the northern Arabian Sea are well-known climate archives and record productivity, as well as supply of material from land. Here, we studied sediment core 275KL off Pakistan in concert with sediment trap, dust and river samples in order to characterize and quantify land-derived material deposited in varves and event layers. We analysed grain sizes, mineral assemblages, bulk components and stable isotopes (δ13C, δ18O) of carbonates. In winter, enhanced river discharge is the main source of lithogenic matter contributing the major amounts to the total annual sedimentation of the northern Arabian Sea. During the late summer season, lithogenic matter accumulation is slightly enhanced, probably carried along with the south-eastward blowing Levar winds from the Balochistan and the Sistan Basins and the summer monsoon discharge maximum of perennial streams. C/N ratios and stable carbon and oxygen isotopes could not be used to distinguish between organic matter produced on land and in the ocean, whereas stable carbon and oxygen isotope ratios of carbonates suggest that sedimentation of event layers is dominated by direct inputs from land. Catastrophic denudation and storm events occur on average once every 50 years and lead to sedimentation rates that exceed the mean annual sedimentations of 983 g m−2 yr−1 by 6 to 10 times. Nevertheless, due to their rare occurrence, they contributed only 7% to the total sedimentation during the last ca. 5000 years. End-member modelling of grain sizes in accordance with lithogenic matter accumulation rates and event layer frequencies showed that arid conditions prevailed between 4000 and 5000 a BP while more humid conditions commenced around 2000 ka BP in accordance with the Pacific ENSO record.
Keywords
Introduction
The Indian summer monsoon is the main driver of climate variability feeding the perennial streams of southern Asia and despite strong efforts still difficult to predict (Gadgil, 2003; Jayaraman, 2002). A source of precipitation in central southwest Asia are the so-called western disturbances (WD) entailing episodic winter rains (Syed et al., 2010) and strongly affecting winter discharge of the Makran rivers (Lückge et al., 2001, 2002). Utilizing the excellent highly resolved marine sedimentary archives of the northern Arabian Sea helps to better understand the monsoon system and WD intensities. Northern Arabian Sea sediments are well known for their continuous seasonal varve-like lamination that is characterized by an alternation of dark and bright laminae with the potential of seasonally resolved paleo-reconstructions (Lückge et al., 2002; Von Rad et al., 2002b). Until now, these archives have been used to reconstruct sea surface temperatures and salinities (Böll et al., 2014; Doose-Rolinski et al., 2001; Munz et al., 2015), marine productivity (Böll et al., 2014; Von Rad et al., 1999a, 2002b), as well as inputs of pollen (Ivory and Lézine, 2009) and lithogenic matter from land (Lückge et al., 2001, 2002, 2012; Von Rad et al., 1999a, 2002b). These studies revealed strong teleconnections with North Atlantic climate variations reflecting, for example, Dansgaard–Oeschger and Heinrich events during the last glacial period (Deplazes et al., 2014; Reichart et al., 1997; Schulz et al., 1998).
Contrary to the western and central Arabian Sea, where at present regional dust storms are the main source of lithogenic matter (e.g. Kolla and Biscaye, 1977; Sirocko and Sarnthein, 1989; Stewart et al., 1965), local river discharges are assumed to be the main source of lithogenic matter to the northern and northeastern Arabian Sea (Lückge et al., 2001, 2002; Schulz et al., 2002; Staubwasser and Sirocko, 2001; Von Rad et al., 2002b). However, satellite imagery (Kaskaoutis et al., 2014, 2015; Rashki et al., 2013, 2015) in combination with modelled dust deposition rates (Mahowald et al., 2005) suggest that eolian dust might also contribute to the annual lithogenic matter sedimentation off Pakistan.
Within the laminated sediments, thicker whitish (referred to as C-Event layers) and reddish layers (F-Event layers) frequently interrupt the normal sedimentation of laminae described in detail by Schulz et al. (1996). C-Event layers were identified by Lückge et al. (2002) to be triggered by flash flood events following heavy rainstorms in the hinterland and the coast of Makran, as observed in March 1998 in the Turbat district of southern Balochistan (IFRC, 1999; Lückge et al., 2002). F-Event layers are coarser than C-layers, sometimes graded and are clearly distinct from C-layers by their reddish colour (Schulz et al. 1996). They were suggested to be deposited by turbidity flows (Schulz et al., 1996; Von Rad et al., 1999a) with their main depot centres on the continental slope and abyssal plain throughout the Holocene but with changing frequencies, thickness and grain sizes induced by climate change in the hinterland (Bourget et al., 2010, 2011, 2013). According to Bourget et al. (2010), turbidites are deposited along the Makran coast even at mid- and late-Holocene sea-level highstands but with reduced frequency and thicker turbidites compared with the glacial and early Holocene. During the rising sea levels of the early Holocene, much of the terrestrial material contributed mainly by rivers was stored on the shelf and was only occasionally transferred to the deep sea. This trend to thicker and coarser turbidites became enhanced after the mid-Holocene humid period when the decrease in vegetation led to enhanced erosion (Bourget et al., 2010).
In this study, we investigated sediment samples obtained from sediment core SO130-275KL taken at a water depth of 782 m offshore Pakistan (Figure 1) close to the well-studied core SO90-56KA (Von Rad et al., 1999a). The aim of this study is to characterize the variability of lithogenic matter sources of the last 5 ka at the Makran coast by analysing grain-size distributions (GSDs) in varves and event layers and determine sediment sources by end-member modelling and comparison of results with local eolian dust, Hingol river bed sediments as well as sinking particles from a near-by sediment trap investigation (Schulz et al., 2002; Schulz and Von Rad, 2014). In combination with biogeochemical and mineralogical methods, we identify sources of organic matter, carbonate and lithogenic components and try to better understand the source and depositional processes of land-derived material. We use these results to quantify the significance of sedimentation events compared with annual varve deposition of the past 5000 years as recorded in sediment core 275KL.

Study area in the context of the Arabian Sea realm. Arrows indicate the major wind systems during the southwest monsoon season. The track of the mineral dust-laden Levar wind, which distributes huge amounts of Sistan dust over the Sistan Basin, is redrawn from Rashki et al. (2015). Grey areas signify bedrocks, stippled yellow areas potential dust sources characterized by hydrologic generated sediments (i.e. alluvial and salt soils). Major dust sources are located in the desert basins of Sistan, Kharan, Lut and Jaz-Murian. Shallow hamuns, rivers and wadi streams are illustrated in dark blue, river basins in bright blue. The white circle and black triangle constitute the studied sites of core 275KL and sediment trap EPT-II, respectively. Lineament of red outlined circles and triangles mark coordinates for point data sampling sites illustrated in Figure 2.
Study area
The Arabian Sea is strongly influenced by the Indian monsoon, which is driven by the differential summer-heating and winter-cooling of the Asian land mass and the Indian Ocean (Ramage, 1987). In summer (SW monsoon season), strong heating forms an atmospheric low over central Asia that attracts the SE trade winds to form a low-level cross-equatorial jet stream, blowing from the southwest over the Arabian Sea (Figure 1, Findlater, 1977). This so-called Findlater Jet causes upwelling off Oman, which is the major process in the Arabian Sea controlling biological productivity and the associated export of organic matter from the sunlit surface ocean into the deep sea during the summer monsoon (Brock et al., 1991, 1992; Rixen et al., 2000, 2009). The monsoon-driven export of organic matter and its respiration in combination with a restricted circulation within the permanent thermocline form and maintain an oxygen minimum zone, which is most pronounced at water-depths between 100 and 1000 m in the northern part of the Arabian Sea (Bogus et al., 2012; Fischer et al., 2013; Gaye et al., 2013; Rixen et al., 2014; Rixen and Ittekkot, 2005; Sen Gupta et al., 1976; Wilhelms-Dick et al., 2012). Off Pakistan, where the oxygen-deficient water flushes the continental slope, dark-olive and bright alternating laminated sediments are formed due to a reduced bioturbation in the absence of higher organisms avoiding such oxygen-deficient waters (e.g. Von Rad et al., 1999b).
Our studied core is located about 60 km south of the eastern Makran coast, at the easternmost part of the active Makran accretionary wedge that is formed by the subduction of the Arabian under the Eurasian Plate. The Makran subduction zone is an area of frequent earthquake activity with a recurrence time of high-magnitude earthquakes of 125–250 and a presumed part of the two main tsunamigenic zones in the Indian Ocean (e.g. Fischer et al., 2013; Heidarzadeh et al., 2008, 2009; Page et al., 1979), which left traces in the sedimentary records along the shores of the northern Arabian Sea (e.g. Donato et al., 2009; Koster et al., 2014). However, sediment supply to the submarine environment of the active Makran margins is preferentially controlled by seasonal short river systems that are very sensitive to climate changes (Bourget et al., 2010, 2011, 2013).
The bathymetry off Makran is composed of a dense dendritic network of prominent river-erosive channels or canyons that cut a high tectonic relief of accreted ridges and intra-slope basins. These channels were most probably incised during sea-level lowstands and are located directly off the many ephemeral streams and the few perennial rivers such as the Hingol and Hab (Bourget et al., 2011). Even at the present sea-level highstand, these channels funnel gravity flows of material intermittently stored on the shelf and upper slope (Bourget et al., 2010, 2011). Multi-beam bathymetry reveals that our core study site is located on the interfluve between two small canyons that were possibly incised by the river systems of the Phor and Porali rivers (compare core site 56KA of Figure 2 in Bourget et al., 2010). Besides the Phor and Porali rivers, the perennial Hingol river is suggested as the main sediment supplier to the coring site (Lückge et al., 2001), whereas the Hab and Indus rivers do not contribute significantly as the prevailing clockwise surface circulation diverts their river plumes to the southeast during the SW-monsoon induced discharge maximum (Lückge et al., 2001; Von Rad and Tahir, 1997).

Comparison of averaged monthly (a) TOMS AI of northern and central Makran, (b) total lithogenic matter (LM) flux of EPT-II between 1993 and 1996, (c) sum of clay and fine silt (red dots) and medium and coarse silt (black dots) as proxies for distal dust and fluvial mud. Blue bars = NE monsoon season, red bar = SW monsoon season.
The land area north of the Arabian Sea (Balochistan) covers large parts of Pakistan and Iran (Figure 1). Balochistan and the northerly adjacent area of Sistan are important sources of dust and river runoff because of their dry climate, tectonic activities, low or absent soil protection by vegetation, weakly consolidated geological strata and the steep topographic gradient between the coastal mountain ranges and the coast (e.g. Field et al., 2009; Lückge et al., 2002; Milliman and Farnsworth, 2011). Contrary to the monsoon-dominated areas to the south of Balochistan, which experience most rainfall during the summer monsoon, the rainy season in the northern and western parts of Balochistan occurs in winter between December and February/March. This is caused by the southward shift of the climate zones in winter, which brings Balochistan into the influence of the upper tropospheric subtropical westerly jet (SWJ). The SWJ carries humid cyclonic storms from the Mediterranean and the Atlantic Ocean eastward and feeds the winter rainfall over Balochistan often referred to as WD (e.g. Dimri et al., 2015; Hastenrath, 1996; Prasad and Enzel, 2006; Roberts and Wright, 1993; Weyhenmeyer et al., 2000). At present, the Hingol river marks the frontier between areas influenced by summer and winter rains and is accordingly fed by precipitations during both seasons (Snead, 1968). Nevertheless, winter precipitation predominates and occasionally reaches alarming quantities causing highly erosive flash floods that dump their loadings into the Arabian Sea. Processes triggering these events are still poorly understood, but there are attempts to link their occurrence with anomalies in the coupled ocean atmospheric circulation such as the North Atlantic Oscillation (NAO) and the El Niño Southern Oscillation (ENSO) (e.g. Dimri, 2013; Dimri et al., 2013; Syed et al., 2006, 2010). Positive phases of NAO are assumed to cause drought in the Eastern Mediterranean region as well as in central and southern Iran, but lead to excess precipitation further to the north-east over the Hindu Kush Mountains in Afghanistan and the north-western parts of Pakistan (Syed et al., 2010). Excess precipitation over this region and the western Himalayas are also associated with positive ENSO phases, that is, El Niño events (Dimri et al., 2015; Syed et al., 2006, 2010; Yadav et al., 2013). During a strong El Niño-associated precipitation event in March 1998, a devastating flood that formed a thick C-Event layer off Ormara could be observed (IFRC, 1999; Lückge et al., 2002).
During the summer season, a northerly low-level jet, called the Levar or 120-day winds (up to 200 km/h; Alizadeh-Choobari et al., 2013), develops due to the strong summer heating that enforces a convective low-pressure system (Figure 1; Rashki et al., 2015). The Levar winds mobilize easily erodible silts and clays from dried-up rivers and lakes, which often leads to heavy dust storms (e.g. Prospero et al., 2002; Rashki et al., 2013; Rezazadeh et al., 2013; Whitney, 2006). In the Sistan Basin, which is considered to be one of the main dust sources in the region, the number of dust storms often exceeds 80 during a single summer monsoon season (Goudie and Middleton, 2000; Middleton, 1986). The dust storm intensities in the Sistan Basin are strongly linked to the aerosol loading over the northern Arabian Sea (Kaskaoutis et al., 2014). At the confluence of the Levar winds carrying the dust southward towards the Makran coast and the Findlater Jet, decreasing wind speeds favour dust deposition and restrict dust dispersion within the monsoon system (Figure 1; Sirocko and Sarnthein, 1989).
Materials and methods
Material and accumulation rates
We investigated the upper 428 cm of the piston core 275KL (24°49.31′N, 65°54.60′E; 782 m water depth), which was retrieved in 1998 from the oxygen minimum zone at the Pakistani continental slope during cruise SO130 with the German research vessel RV Sonne (MAKRAN II Cruise Report, 1998). The sediments of the piston core show no indications of erosion by turbidity or bottom currents and thus characterize 275KL as an almost undisturbed record of Holocene sedimentation off Pakistan. The PAKOMIN sediment trap EPT-II (24°45.6′N, 65°48.7′E; 537 m water depth) was moored close to the core site of 275KL (Figure 1, about 12 km distance) from May 1995 to October 1996. Due to technical problems (clogging), the obtained sediment trap record covered only the period between May 1995 and February 1996 (Figure 2; Andruleit et al., 2000). Annual averages used in this work were calculated based on the available data set of about 10 months. As we will discuss below, the effect of undertrapping by shallow sediment traps is more significant than the effect of missing 2 months of sampling. As an example of Makran rivers, we analysed lower Hingol river flood terrace sediments provided by the BGR Hannover, Germany (for further references, see Lückge et al., 2001). Geochemically and mineralogically, these flood terrace sediments strongly resemble lithogenic matter of central Makran obtained from sediment core GeoB12309-5 off the Dasht river (24°25.3’N, 62°59.9’E; 956 m water depth; Wilhelms-Dick et al., 2012).
The sediment core 275KL was sampled in approximately 3-cm intervals and the white and reddish events were sampled separately. Chronostratigraphy is based on eye-to-eye correlations with 99 event layers of core SO90-56KA from the same location (distance to 275KL < 1.5 km), which was dated by ‘varve’ counting in combination with the AMS 14C method and conventional radiocarbon ages (detailed methodology described by Munz et al., 2015; Von Rad et al., 1999a). The standard deviations of the AMS 14C dates and conventional radiocarbon ages were between 31–51 years and 140–205 years (Von Rad et al., 1999a), respectively, but the counting of annual varves and two additional exact dates of the Krakatoa and Tambora ash layers resulted in virtually exact ages for SO90-56KA. The obtained ages are displayed in years before present (yrs BP), meaning before AD 1950 (Supplements S1 and S2, available online). The ages of C- and F-layers of SO90-56 KA were transferred to the event layers found in core SO130-275 KL and the ages between these events were linearly interpolated. The age differences between individual event layers are between 2 and 441 years. The depths and ages of C- and F-Event layers are given in Supplement S1 (available online) and the core photograph with ages of C- and F-layers indicated is provided in Supplement S2 (available online).
Total mass accumulation rates (MAR) of the sediments were calculated by multiplying the linear sedimentation rate (LSR) with the dry bulk density (DBD). LSR (cm/yr) was determined using age model depth and age interval values, whereas DBD (g/cm3) was calculated as described by Avnimelech et al. (2001):
where
For the calculation of sediment particle density, we assumed 2.65 g/cm3 as conventionally used for inorganic sediment particles (Blake and Hartge, 1986; Boyd, 1995) and a density of 1 g/cm3 was taken for water. Although organic matter weight percentages are below 5%, we also performed a density correction for organic matter contents by assuming a density of 1.25 g/cm3 for organic particles (Boyd, 1995):
Prior to analysis of bulk components (organic matter, carbonate, biogenic opal and lithogenic matter) and mineral assemblage of lithogenic matter, the sediments were wet weighed, freeze-dried, dry weighed and finally an aliquot was ground and homogenized for analyses (see Böll et al. 2014 for details). Accumulation rates of the bulk components were determined by multiplying the total MAR with the analysed weight percentages of individual bulk components.
Mineral assemblages
X-ray diffraction (XRD) analysis on six C-Event samples of the piston core was done at the University of Bremen (Central Laboratory for Crystallography and Applied Material Sciences, ZEKAM, Department of Geosciences). An amount of 1 g of dried and ground bulk sample material was prepared with the Philips backloading system for measurements of x-ray diffractograms. Measurements were conducted on a Philips X’Pert Pro multipurpose diffractometer with a Cu-tube (1.541 ka, 45 kV, 40 mA), a fixed divergence slit of ¼°, a secondary Ni filter and the X’Celerator detector system. Measurements were done as continuous scans from 3 to 85° 2 ɵ (calculated step size of 0.016° 2ɵ with 50 s calculation time per step). For mineral identification, Philips software X’Pert HighScoreTM was used; for identification of sheet silicates, MacDiff 4.25 was used (see Petschick et al., 1996 for details). Quantification of mineral assemblages was done with the full-pattern QUAX© software package (Vogt et al., 2002). Precision for the bulk mineral determination is 1% for quartz, 2–5% for feldspar and 5–10% for clay minerals.
GSD parameters
For grain-size analysis of sediment core material and Hingol river bed sediment, 250 mg of the bulk sediment was freed from organic matter and carbonates by treatment with H2O2 (boiled up) and HCl, which act as main cements of most sediments (e.g. Sun et al., 2002 and references therein). Afterwards, the sediments were boiled in water and Calgon® (Na6P6O18) dispersant to avoid particle coagulations. Ultrasonification was conducted before measurements by application of the particle size analyser. Grain size was measured on a Horiba LA-950 V2 laser-scattering particle size distribution analyser. GSDs were calculated for 93 grain-size classes within a measuring range of 10 nm to 3 mm. Measurements were conducted three times per sample and then averaged. Mean grain-size spectra were calculated for varves, C- and F-Events layers and the standard deviation (1σ) from the mean grain-size spectra is ±19% for varves, ±25% for C- and ±42% for F-Event layers, respectively. The GSD parameters were calculated using the GRADISTAT programme (Blott and Pye, 2001). End-member modelling was conducted using an end-member modelling analysis (EMMA) script for MATLAB (Dietze et al., 2012).
To test a potential impact of sediment alteration by the treatment method used on GSDs, 12 samples were measured without former chemical treatment (Supplementary material S3–S5, available online). We observed no systematic difference between the GSD in HCl-treated and untreated varve samples and surmise that this is due to the presence of carbonaceous shells of different species and sizes in addition to the more sorted lithogenic or pedogenic carbonates. There is a systematic difference between HCl-treated and untreated C- and F-layers. C- and F-layers lose material in the clay and fine-silt fraction at 0.3–0.4 µm and 2–3 µm after HCl treatment. F-layers and one C-layer loose material at 15–30 µm and F-layers lose material at 70 and 500 µm after acid treatment. We surmise that the loss in the sand fraction is related to destruction of large aggregates, whereas carbonate dissolution occurs in the clay and fine-silt fraction. We did not remove biogenic opal from the sediments as the percentage was below 1% in all analysed samples with an average of 0.6% (see method and results in Böll et al., 2014) and, therefore, the effect of biogenic opal on the grain-size analyses is negligible.
Grain-size analysis of EPT settling sediments was conducted at the University of Hamburg, Germany, with a GALAI CIS 100 laser-scattering particle size analyser on particle sizes between 0.5 and 100 µm on a resolution of 0.5 µm (Tiemann, 2001). Measured samples were collected on grain-size ranges only (<2, 2–6, 6–20 and 20–63 µm) hindering direct comparison of GSDs with the other sediment types due to lacking resolution and the different instruments used.
The aerosol GSD data for Karachi were obtained from the Aerosol Robotic Network (AERONET; https://aeronet.gsfc.nasa.gov/new_web/index.html).
Note that due to the different methodologies, we can only discuss qualitative differences and trends in GSDs when comparing sediment trap and dust samples with Hingol river bed and sediment core samples.
Stable oxygen and carbon isotope ratios
Carbonate stable oxygen (δ18O) and carbon (δ13C) isotope measurements were conducted on aliquots of homogenized bulk C-Event layer material with a Finnigan Kiel IV Carbonate Device coupled to a Finnigan MAT251 mass spectrometer. For δ13C measurements on organic material of bulk sediment samples, we used a Finnigan Delta Plus gas isotope ratio mass spectrometer after removal of carbonates by 1-N HCL and subsequent drying at 40°C. δ13C and δ18O of carbonates and organic material are given as ‰ deviation from the composition of the VPDB standard. The precisions measured by internal standards were −0.01‰ for δ18Ocarb, 0.11‰ for δ13Ccarb and 0.14‰ for δ13Corg.
Nimbus-7 aerosol index
We used the mean monthly Total Ozone Mapping Spectrometer (TOMS) aerosol index data provided by the World Climate Research Programme CLIVAR (Climate Variability and Predictability; www.ouce.ox.ac.uk/~clivar) over central Makran to investigate the mean monthly distribution of aerosols over Balochistan. The original data for index calculations are based on multispectral measurements in the ultraviolet (UV) wavelengths on board of the Nimbus-7 satellite between 1978 and 1993, received from the TOMS data page of the NASA (http://disc.sci.gsfc.nasa.gov/aerosols/data-access). The calculations for the aerosol index are described in detail by Herman et al. (1997).
Results and discussion
Varve accumulation and composition
The Levar winds develop parallel to the summer monsoon and control the dust concentration in the atmosphere over Balochistan and the coast of Makran (e.g. Rashki et al., 2015; Rezazadeh et al., 2013; Snead, 1968) with dust maxima between May and August (Figure 2a). Back-trajectory analysis identified the Sistan Basin as a major source region for aerosols over Karachi (Alam et al., 2012; Rashki et al., 2015). Grain-size analyses of aerosols collected in summer over Karachi show a unimodal distribution with a mean mode at 2.6 µm (Figure 3a), which is generally considered to be characteristic for distal eolian dust (e.g. Mahowald et al., 2011, 2013; Pye, 1987; Seinfeld and Pandis, 2006). Contrary to distal Sistan dust, the Hingol river flood terrace sediment is coarser with a bimodal GSD (Figure 3a). The mean modes at 9 and 42.1 µm probably reflect coarse suspended matter (Prins and Weltje, 1999) and bed load material (e.g. Billi, 2011; Reid and Laronne, 1995), respectively. Flood terrace sediments from other Makran rivers (not shown) are even coarser than the Hingol sample. However, the fine fraction is possibly not deposited on terraces and the river bed and also passes by the proximal cores off the river mouths. Therefore, it is underestimated in the river bed samples and proximal cores.

Comparison of grain-size distributions of (a) distal eolian dust measured by the ground-based remote-sensing method over Karachi (source: www.aeronet.gsfc.nasa.gov) and of the lower Hingol river flood terraces, (b) varves, (c) C-Event layers and (d) F-Event layers. Bold curves represent averages and grey curves represent individual samples.
Mineral particles that are carried as eolian dust and as river suspended matter into the surface layer form the lithogenic matter fraction of sinking particles in the ocean (Honjo et al., 1982). The sediment trap EPT-II collected particles at 590 m water depth, which is approximately 1400 m above the seafloor, so that we can assume that its major source is particulate matter export from the surface layer (Schulz and von Rad, 2014). The lithogenic matter fluxes show a bimodal pattern with the maximum in winter (January–February) and a secondary maximum in late summer (August–September; Figure 2b). During the observed summer monsoon of 1995, the lithogenic matter fluxes peaked in September 1 to 2 months after the atmospheric dust concentrations reached their maximum over marine regions close to the Makran Coast (Figure 2a). On the one hand, the delay could be explained by a weakening of the Levar winds during the declining summer monsoon that allows the dust to fall out (Rashki et al., 2015). On the other hand, the perennial Makran rivers as well as the Indus river (Ittekkot and Arain, 1986) also have SW-monsoon fed discharge maxima during August–September so that a riverine source of the fine material deposited in late summer cannot be ruled out. During winter, grain-size maxima suggest that the high lithogenic fluxes were associated with precipitation events and resulting discharges with flash floods of the Makran rivers (Lückge et al., 2002) (Figure 2c).
However, the seasonality of sedimentation of the fine lithogenic fraction is not solely controlled by the supply from eolian or fluvial sources but is driven by the seasonality of primary productivity as the export of fine-grained lithogenic material depends on its incorporation into fast sinking aggregates (Degens and Ittekkot, 1984; Ittekkot, 1993). These aggregates, referred to as marine snow and faecal pellets, are formed during plankton blooms (Alldredge and Silver, 1988; Osterberg et al., 1963; Silver and Alldredge, 1981; Suzuki and Kato, 1953), which are triggered by mixed-layer deepening caused by winter-cooling, and advection of upwelled water from the Oman coast in summer (Rixen et al., 1996). The basic pattern of lithogenic sedimentation is thus a combined effect of the lithogenic matter supply and the monsoon-driven plankton blooms producing aggregates, which are the vehicles for the export of fine lithogenic matter into the deep sea.
Sediment traps are so far the only tool to measure the export of particles from the surface into the deep ocean and quantitative comparisons between trap results and sediments are biased by a variety of hydrodynamic, biological and chemical constraints (Buesseler et al., 2007 and references therein). Short-lived radionuclides and inverse modelling techniques used to estimate the accuracy of sediment traps indicate a generally strong undertrapping at continental slopes at depths between 200 and 1000 m, with trapping efficiencies of 30–60% (Buesseler et al., 2007; Schlitzer, 2000; Scholten et al., 2005; Usbeck et al., 2003; Yu et al., 2001). By taking a potential undertrapping of EPT-II between 30% and 60% into account, the corrected measured annual mean flux of lithogenic matter would rise from 357 g m−2 yr−1 to values between 510 and 893 g m−2 yr−1, respectively. This is well within the range of calculated varve MARs during the last 3700 years, which have a mean annual lithogenic matter MAR of 781 g m−2 yr−1 (Table 1). The composition of EPT-II samples is similar to varve compositions with, however, significant reductions in organic carbon and biogenic opal due to degradation and, respectively, dissolution in sediment (Table 1). Therefore, we assume that varves depict the normal sedimentation dominated by vertical flux of particles exported from the surface layer.
Averaged lithogenic matter, carbonate, biogenic opal, organic matter and organic carbon (Corg) components (wt%), mass accumulation and flux rates (g m−2 yr−1) for 275KL (varves, C-layers and F-layers), EPT-II trap sediments, and Hingol river samples (HR), respectively, and stable isotopic ratios of carbonate and organic carbon.
MAR: mass accumulation rates; NA: not analysed.
275KL.
Averaged global values of Cerling (1984), Schlesinger et al. (1989) and Cerling and Quade (1993).
Data of sediment core 53KG; for further information, see Von Rad et al. (1996).
Data from Böll et al. (2014) covering the upper 188 cm of core 275KL (=2370 a BP).
Biogenic matter and carbonates in varves and C-Event layers
C-Event layers differ from normal varve sedimentation by higher lithogenic matter and lower organic matter contents (Table 1). The visual inspection of smear slides shows that remains of marine plankton are abundant in varves, but absent in C-layers, suggesting that C-layers are of terrigenous and varves are of a mixed marine and terrigenous origin (e.g. Lückge et al., 2002; Schulz et al., 1996). This view was supported by microscopic studies (Lückge, 1997), which showed that reworked material of higher land plants are the predominant contributors to the organic matter fraction in C-Event layers. However, our mean C/N values of 4.9 in combination with a mean organic matter δ13C value of −22.1‰ traditionally point to a marine algal source (e.g. Meyers, 1997) (Table 1). Since C/N ratios can be affected by nitrogen adsorbed to clay minerals (Meyers, 1997; Müller, 1977), a nitrogen correction was performed by plotting the organic carbon/clay size (particles < 2µm) ratios versus C/N ratios following the method of Hartmann et al. (1971) (Figure 4). The Corg-poor C- and F-Event layers, as well as the Hingol river material, clearly show a significant decrease in C/N below a Corg/clay ratio of about 0.14 probably due to an enrichment of nitrogen sorbed within clay minerals. C/N ratios of samples with a Corg/clay ratio >0.14 are relatively stable and can thus be assumed to represent C/N values of unprotected organic matter (Müller, 1977). However, varves and EPT-II sediments show no relationship between C/N and Corg/clay ratios and thus indicate no or only minor impact of organic matter associated with clay minerals on nitrogen contents (Figure 4). The nitrogen correction results in C/N ratios of 6 to 7.5 in C- and F-Event layers, which are close to those of phytoplankton and of the organic matter in varves and EPT-II (see Table 1). However, these values are also within the range of C/N ratios measured in top soils of the north-eastern Tibetan Plateau (ratios of 5–25; Zhao et al., 2015). The δ13C values of degraded Tibetan soil organic matter varied between −20‰ and −26‰ so that the δ13C values of −22‰ and −21‰ measured in C-Event layers and varves are within this range. C/N ratios and δ13C values can, therefore, not be used to distinguish between organic matter produced by phytoplankton and eroded from terrestrial soils.

Relationships between C/N wt% ratios and organic carbon/clay ratios (×103) as graphical correction for ammonium ions adsorbed on particles of the carbonate free clay size fraction <2 µm, for (a) Hingol river sediments, C- and F-Event layers and (b) varves and EPT sediments. The red arrow indicates the organic carbon/clay threshold between values affected by nitrogen fixed within clay minerals to the left and independent values reflecting C/N ratios of unprotected organic matter to the right. Method adopted from Hartmann et al. (1971).
Carbonates, which contribute approximately between 12% and 16% to the varves, the C-Event layers and the sediment trap material offer a possibility to distinguish between marine and terrestrial sources as their stable isotope ratios differ (Table 1). The main carbonate producers in the Arabian Sea are coccolithophorids and foraminifera (Andruleit et al., 2000; Ramaswamy and Gaye, 2006; Schiebel et al., 2004). The stable isotope ratios (δ18O, δ13C) of foraminifera picked from varves (unpublished data by P Munz) and other cores from the Arabian Sea (e.g. Lückge et al., 2001; Staubwasser et al., 2003) are ~ −1.5‰ and −0.8‰, respectively (Table 1). These values correspond with other stable isotope measurements of planktic foraminifera and coccolithophorids of the Holocene which basically preserve the stable isotopic values of water and the CO2 source with some species-specific differences due to vital effects (e.g. Fink et al., 2010; Rostek et al., 1997). Fluctuations in foraminiferal δ18O can be affected by ice volume, temperature, salinity and circulation changes (Rajeev et al., 2012; Waelbroeck et al., 2006). The carbonates of the C-Event layers are, however, more depleted in 18O and 13C than observable in recent marine carbonates with δ18O and δ13C values of −6.4 ‰ and −3.1‰, respectively (Table 1). Stable isotope measurements of bulk carbonates in varved sediments from box core SO90-53KG (24°48.6′N, 65°54.9′E; water depth 850 m; Von Rad et al., 1996), taken about 1.5 km from site 275KL, show δ18O ratios that are identical to those measured in C-Event layers and δ13C ratio of −1.8‰ between measured δ13C in foraminifera (−0.8‰) and C-Event layers (−3.1‰, Table 1). The differences between isotopic signatures seen in foraminifera and the bulk carbonates and the similarity between those seen in bulk carbonates and C-Event layers suggest that the carbonate sedimentation is substantially driven by isotopically light carbonates from land. Potential terrigenous sources for carbonates from Balochistan are marine outcrops along the Makran Ranges, evaporitic precipitates in hamuns or playas, and pedogenic soil carbonates (also known as caliche). Carbonates that formed under high evaporative conditions are generally enriched in the heavier isotopes 18O and 13O and thus reveal higher δ18O and δ13C values than those observed in C-Events (e.g. Quade et al., 2007; Stiller et al., 1985). Accordingly, they can be excluded as potential source material. In contrast, caliche of modern arid to sub-humid environments show a good accordance with the carbonates of varves and C-Event layers with a global mean δ18O value of −7.6‰, and a global mean δ13 C value of −4.7‰ (average values of Cerling, 1984; Cerling and Quade, 1993; Schlesinger et al., 1989). Soils from modern fluvial systems on the Indian subcontinent (Basu et al., 2018) have similar δ18O and δ13C values.
Mineral assemblage and grain size of varves and event layers
Silicates of C-Event layers are predominated by clay minerals with illite and chlorite as main contributors resembling the Hingol river sediments and other sediments of the Gulf of Oman (Table 2; Prins et al., 2000; Sirocko and Lange, 1991; Stewart et al., 1965). Chlorite is the characteristic of the Makran region and is restricted to areas close to its metamorphic or sedimentary source region (Kolla et al., 1981). The absence of palygorskite and only low contributions of smectites indicate in line with other results that dust inputs from central Arabia via the Findlater Jet into the northeastern Arabian Sea are of minor importance (Kolla et al., 1976, 1981; Sirocko et al., 1991; Sirocko and Lange, 1991).
Average mineralogical composition of potential sediment sources, mean annual EPT-II samples, surface samples (0–2 cm), varves of sediment core GeoB12309-5 and C-Event samples of 275KL. Values are separated into principal and clay components and normalized to 100%.
Average of three sample locations of airborne dust in the vicinity of Gwadar.
Sample of lower Hingol river flood terrace sediments.
Average of two sample locations of the lower of the Indus river plains.
Average of 101 varve samples of sediment core GeoB12309-5.
Average of six C-Event samples of core 275KL.
Average of four mixed-surface samples taken close to the EPT and 275KL sites during METEOR cruise 32/2.
Average of 27 EPT samples representing SW and NE monsoon seasons 1995/96.
Mixed organic and terrigenous calcite.
Chlorite decreases and illite increases over time and converge to values of the recent Hingol river sediments.
GSDs of silicates show uni- and bimodal patterns in varves and C-Event layers (Figure 3b and c). The fractions of silt and clay predominate in varves and C-Event layers by 97% (SD = 0.05%) and 99% (SD = 0.02%), respectively. Varves have the first or only mode (if unimodal) at about 4.8 µm and the second mode at about 0.2 µm. Some varve samples are coarser with modes between 6 and 10 µm. C-Event layers have slightly coarser modes than varves at 0.3 µm and around 7.2 µm. F-Events are mostly unimodal with modes between 8 and 11 µm. The only exception is a sample from the base of the youngest F-layer which has a second mode at 37 µm in the coarse silt fraction corroborating that the F-layer is sorted and may originate from a turbidity flow.
For the end-member modelling with EMMA (see details of the method in Dietze et al., 2012) carried out on varve samples, three samples were excluded and are also not plotted in Figure 3b as they are much coarser and have modes diverging from all samples and would thus have biased the end-member reconstruction. EMMA produces three main end-members which explain 87% of the total variance. End-member (EM) 1 has modes at 0.2 and 5 µm, EM2 at 2.3 and 8 µm and EM3 at 0.3 and 2.3 µm, respectively (Figure 5a). As Sistan dust has its mode at 2.6 µm (Figure 5c), we assume that EM3 represents the distal dust source. The mode of EM1 is close to the mode of about 6 µm found in the EPT trap during the summer when monsoon impacted rivers have their discharge maxima (Tiemann, 2001). Therefore, we believe that EM1 represents river suspension. The mode of EM2 is close to the mode of F-layers (Figure 5b) and to the first mode of Hingol river bed sediments (Figure 5c; 9.5 µm; second mode: 42 µm) supposing that EM2 represents material originating from high energy transport mobilizing coarser material from the river beds or from the continental shelf and transporting it to the depot centres on the slope.

(a) Factor loadings of three end-members modelled with EMMA (Dietze et al., 2012) using grain sizes of varves compared with (b) mean grain-size distribution of varves, C-Event layers and F-Event layers and (c) grain-size distribution of Sistan Dust and Hingol river flood terrace samples.
Grain-size spectra are not only affected by lithogenic matter sources but also by subaquatic transport energy (Bass et al., 2002; McCave, 2017; McCave and Hall, 2006). The absence of marine organisms and abundance of terrestrial organic matter imply an almost exclusive riverine origin for C-layers in line with results of other authors using different methods (Lückge et al., 2001; Von Rad et al., 1999a, 2002a). Therefore, they must be transported directly from land by episodic river floods and probably reach the slope by under- or interflows (Von Rad et al., 2002a). This assumption was also made by Bourget et al. (2010) based on investigations of sediment cores west of 275KL. Observations of hyperpycnal flows triggered by high-density river discharges into the northern Red Sea (Katz et al., 2015) suggest a likely mechanism for fast accumulation of terrigenous material that probably can explain the C-Event formation.
In line with results derived from bulk composition and stable isotopes, our results further imply that varves represent the ‘normal’ annual sedimentation with mixed marine and terrestrial inputs from both transport mechanisms fluvial and eolian into the surface water, subsequently subjected to packing and scavenging into biogenic settling aggregates (see above).
F-Event layers
Contrary to Schulz et al. (1996), who assumed F-Event layers to be exclusively of terrigenous origin, we also observed residues of marine plankton such as foraminifera and diatoms in the reddish sediments pointing to mass redistribution of marine and coastal deposits. Staubwasser and Sirocko (2001) stated that the geochemistry of the F-Event layers closely resembles that of sediments from the lower Indus river plains. Von Rad et al. (2002b), on the other hand, argued that the Indus river is a rather unlikely source for F-Event layer formation due to, for example, the long distance between the Indus mouth and the core site, and the prevailing current direction that diverts the Indus load southward. Due to their restricted occurrence in the northeastern-most part of the Arabian Sea (Schulz et al., 1996), F-Event sediments must be sourced by very local sediment deposits. Mean GSD parameters of F-Event layers point to a coarser composition than varves and C-Event layers and show a close relation to the lower Hingol river terrace material (Figures 3 and 5). The Windar and Porali river located between the Hingol and Indus river (Figure 1) carry reddish sediments from the Kirthar range to the Miani Lagoon to form massive subaquatic sand banks. Snead (1966) reported the complete disappearance of such massive sand banks sometimes between AD 1842 and 1890, which would agree to the most recent F-Event formation in the core 275 KL at around AD 1887. Local fishermen reported that heavy surf during the SW monsoon season and also large floods causing strong currents in the lagoon regularly remove and redistribute the bar sediments. A direct transport of that huge amount of sediment to the slope during one such flood event seems unlikely due to (1) the distance of approximately 90 km between the lagoon outlet and site 275KL, (2) the F-Event layer thicknesses of up to 9 cm and (3) the low mean event layer recurrence interval of ca. 491 (between 76 and 1041 years) during the past ca. 5000 years. Thus, we propose primary deposition of the reddish material on the shelf facilitated by alongshore currents that converge in front of the lagoon outlet from the north and the south creating an offshore suction (compare Figure 2 in Snead, 1966). Submarine landslides triggered by large-magnitude earthquakes (Fischer et al., 2013; Hampton et al., 1996) have mean calculated recurrence intervals between 175 and >300 years along the Makran subduction zone (Byrne et al., 1992) and are therefore a possible trigger for material re-deposition resulting in F-Event layer formation on the slope. Bourget et al. (2010) identified late-Holocene turbidite recurrence intervals of 200–235 years in regions receiving material directly via canyons, of 224 years in the deep abyssal plain, and of 308–500 years in two western isolated piggy-back basins, which resemble those of F-Event layers. Tsunami and tropical cyclone events may also trigger F-Event layer formation but there is no corresponding historical record connected with the youngest F-Event layer formation (Knapp et al., 2010; Rastogi and Jaiswal, 2006).
Sediment mass budget
In total, about 4.6 t/m2 sediments were accumulated at site 275KL during the past 4882 years producing a mean MAR of 983 g/m2 yr. Lithogenic matter dominates the total sedimentation with 78 wt%, followed in decreasing order by carbonates (14 wt%), biogenic opal (6 wt%) and organic matter (2 wt%). Varve sedimentation is the major contributor and accounted for 89.3 wt% of the total sediment accumulation with a mean rate of 839 g/m2 yr. Although a single F-Event layer comprises mean lithogenic matter and carbonate accumulation of 10 years of varve sedimentation, it contributes only 3 wt% to the total sediment amount of the past 5000 years. A C-Event layer deposits on average adequate to 6 years of normal varve sedimentation, but accounts for 7 wt% of the total sediments due to a higher recurrence interval of every 50 years compared to 491 years for F-Events.
Variations in mid- and late-Holocene sediment sources and implications for climate change
Lithogenic matter accumulation rates were extremely high with maxima between 10,000 and 18,000 g m−2 a−1 and C-Event layers were thick and frequently deposited prior to about 4000 a BP (Figure 6a and b). The results of the end-member mixing model showed that the varves were dominated by distal dust during this period (Figure 6c). This coincides with much drier conditions in the westerly wind belt as well as in the Indian monsoon zone proposed for the period from about 5000 to 4000 yr BP (Lauterbach et al., 2014; Prasad et al., 2014). This prolonged dry spell following the Holocene climatic optimum was explained by a reduced import of summer moisture by mid-latitude westerlies (Lauterbach et al., 2014) and with a weakening of the Indian monsoon precipitation due to warming of the Indo-Pacific warm pool (Prasad et al., 2014). More arid conditions in both climatic realms, the westerly wind and the Indian monsoon region, could explain an enhanced supply of mineral matter from land due to enhanced erosion (e.g. Milliman and Farnsworth, 2011). After a period of strongly reduced C-Event layer deposition and a reduction in total accumulation rates to values below 1000 g m−2 a−1, river suspension became the dominant constituent of varves after 2800 a BP indicating the return to more humid conditions (Figure 6c). As C-Event layers started to increase later at about 2300 a BP, we surmise that SW-monsoon precipitation increased earlier than the episodic winter rains producing C-Events layers. The second phase of elevated lithogenic matter accumulation rates between 1000 and 2000 a BP consisted mainly of fluvial material and coincides with a phase of enhanced occurrence of El Niño events (Figure 6d; Conroy et al., 2008) and with probably at the same time occurring positive NAO (Fletcher et al., 2012). Warm ENSO as well as positive NAO phases seem to significantly increase winter precipitation in the Makran region as it was suggested by models and recent observations (Dimri et al., 2015; Syed et al., 2006, 2010; Yadav et al., 2013). Enhanced precipitation during positive phases of the NAO may be triggered by a strengthened SWJ enhancing moisture transport from the Mediterranean and Atlantic (Syed et al., 2010). Warming of the Eastern Pacific during El Niño drags warm water from the western Pacific warm pool into the western Indian Ocean. This leads to a strengthening of the Hadley cell and at the same time to a stronger SWJ at the poleward end of the Hadley cell so that the WD induced winter rains over central southwest Asia are enhanced (Dimri, 2013; Dimri et al., 2016). These rains are additionally fed by enhanced moisture transport from a strong convective cell evolving from the abnormally warm western Indian Ocean (Dimri et al., 2015; Yadav et al., 2013). The establishment of a link with ENSO was also observed in the core monsoon zone at 2000 a BP where the positive ENSO phase leads to severe droughts (Prasad et al., 2014). It is very likely that the recent fluctuations in accumulation rates may be relatable to the intensities of winter precipitation and the West Pacific ENSO records (Figure 6a and d), whereas the accumulation rate maximum between 4000 and 5000 a BP was caused by draught conditions.

(a) Total lithogenic mass accumulation rates (LM MAR) of varves (g m−2 a−1); (b) occurrence and thickness of C-Event layers (cm); (c) end-member scores of grain-size samples for the three end-members modelled with EMMA with EM1 scores in blue, EM2 scores in red and EM3 scores in black; (d) ENSO record from a Galapagos lake with higher percentages of sand indicating wetter (positive ENSO) conditions (Conroy et al., 2008).
Summary
Our results show that in summer, the south-eastward blowing Levar winds carry lithogenic matter, predominantly originating from Balochistan and the Sistan Basins, southward as dust to the northern Arabian Sea realm. In winter, river discharge is the predominant source of lithogenic matter, also on annual time scales. River discharge is fed by humid air masses that are transported eastward from the Mediterranean and the Atlantic Ocean and vary in concert with the NAO and ENSO due to associated changes in the configuration of the SWJ stream over southwest Asia in winter. End-member modelling of grain sizes and variations in lithogenic matter accumulation rates suggest that arid conditions prevailed between 4000 and 5000 a BP, while humid conditions commenced about 2000 a BP and winter precipitation feeding Makran rivers fluctuated in accordance with the ENSO record.
The whitish C-Event layers represent catastrophic denudation events related to flash floods that occurred on the coastal plains of Makran and on averaged recurrence intervals of 50 years. The reddish F-Event layers seem to be associated with massive relocation events of near shore sediments due to heavy earthquakes that resulted in turbidity flows and occurred only about every 500 years. Contrary to the F-Layers that contain material of marine and terrestrial origin, the C-Event layers are composed of mobilized river deposits and soil material. The soil material consists of pedogenic carbonates (referred to as caliche) and degraded soil organic matter. Stable carbon isotope and C/N ratios measured in the degraded soil organic matter resemble those measured in marine plankton showing that these biogeochemical indicators cannot be used to quantify the contribution of terrestrial organic matter in marine sediments. Stable carbon and oxygen isotope ratios measured in the bulk carbonate of varves resemble those of carbonates in the C-Event layers and differ from those measured in shells of marine calcifying organisms, which suggest that lithogenic carbonate dominates within all sediment types along the Makran coast. Although individual C- and F-Event layers supply high amounts of material, their total contribution to sediment accumulation is only about 10% due to their low frequencies.
Supplemental Material
Supplementary_Material_S1-S5 – Supplemental material for Sources of laminated sediments in the northeastern Arabian Sea off Pakistan and implications for sediment transport mechanisms during the late Holocene
Supplemental material, Supplementary_Material_S1-S5 for Sources of laminated sediments in the northeastern Arabian Sea off Pakistan and implications for sediment transport mechanisms during the late Holocene by Sven Forke, Tim Rixen, Nicole Burdanowitz, Andreas Lückge, Venkitasubramani Ramaswamy, Philipp Munz, Dorothee Wilhelms-Dick, Christoph Vogt, Sabine Kasten and Birgit Gaye in The Holocene
Supplemental Material
Supplementary_Material_S6 – Supplemental material for Sources of laminated sediments in the northeastern Arabian Sea off Pakistan and implications for sediment transport mechanisms during the late Holocene
Supplemental material, Supplementary_Material_S6 for Sources of laminated sediments in the northeastern Arabian Sea off Pakistan and implications for sediment transport mechanisms during the late Holocene by Sven Forke, Tim Rixen, Nicole Burdanowitz, Andreas Lückge, Venkitasubramani Ramaswamy, Philipp Munz, Dorothee Wilhelms-Dick, Christoph Vogt, Sabine Kasten and Birgit Gaye in The Holocene
Footnotes
Acknowledgements
The authors thank the BGR Hannover for provision of the untreated SO130-275KL sediment core material and Dr Jörg Tiemann for grain-size data of sediment trap EPT-II. Katharina Pilgram, Dorothee Dasbach, Matthias Birkicht and Dr André Klicpera are thankfully acknowledged for their support in the laboratories. The authors also thank all members of the CARIMA project for their general assistance and openness to discussion. They are grateful to the two anonymous reviewers for their constructive comments and suggestions to improve the quality of the manuscript.
Funding
This work was funded by the German Ministry of Education and Research (BMBF) grant 03G0806C (CARIMA) and 03G0864A (CAHOL). The data are available at PANGAEA (Gaye, Birgit, 2018: Composition and grain size distribution of laminated sediment core SO130_275KL from northeastern Arabian Sea off Pakistan.
).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
