Terrigenous material supply to the Peruvian central continental shelf (Pisco, 14° S) during the last 1000 years: paleoclimatic implications

. In the eastern Paciﬁc, lithogenic input to the ocean responds to variations in the atmospheric and oceanic sys-tem and their teleconnections over different timescales. At-mospheric (e.g., wind ﬁelds), hydrological (e.g., fresh water plumes) and oceanic (e.g., currents) conditions determine the transport mode and the amount of lithogenic material transported from the continent to the continental shelf. Here, we present the grain size distribution of a composite record of two laminated sediment cores retrieved from the Peru-vian continental shelf that record the last ∼ 1000 years at a sub-decadal to centennial time-series resolution. We propose novel grain size indicators of wind intensity and ﬂuvial input that allow reconstructing the oceanic–atmospheric variability modulated by sub-decadal to centennial changes in climatic conditions. Four grain size modes were identiﬁed. Two are linked to aeolian inputs (M3: ∼ 54; M4: ∼ 91 µm on average), the third is interpreted as a marker of sediment discharge (M2: ∼ 10 µm on average), and the last is without an associated origin (M1: ∼ 3 µm). The coarsest components (M3 and M4) dominated during the Medieval Climate Anomaly (MCA) and the Current Warm Period (CWP) periods, suggesting that aeolian transport increased as a consequence of surface wind stress intensiﬁcation. In contrast, M2 displays an opposite behavior, exhibiting an increase in ﬂuvial terrigenous input during the Little Ice Age (LIA) in response to more humid conditions associated with El Niño-like conditions. Comparison with other South American pa-leoclimate records indicates that the observed changes are driven by interactions between meridional displacement of the Intertropical Convergence Zone (ITCZ), the South Paciﬁc Subtropical High (SPSH) and Walker circulation at decadal and centennial timescales.


Introduction
The Pisco region (∼ 14-15 • S) hosts one of the most intense coastal upwelling cells off Peru due to the magnitude and persistence of alongshore equatorward winds during the annual cycle (Fig. 1b). Regional winds can also be affected at interannual timescales by El Niño-Southern Oscillation (ENSO) variability (i.e., enhanced or weakened during La Niña and El Niño events, respectively), as well as by the Pacific Decadal Oscillation (PDO) at decadal timescales (Flores-Aqueveque et al., 2015). These factors also affect the inputs of terrigenous material to the Peruvian continental shelf. Saukel et al. (2011) found that wind is the major transport agent of terrigenous material west of the Peru-Chile Trench between 5 and 25 • S. Flores-Aqueveque et al. (2012) showed that, in the arid region of northern Chile, transport of aeolian coarser particles (approximately ∼ 100 µm) is directly related to in-Published by Copernicus Publications on behalf of the European Geosciences Union.
terannual variations in the domain of the strongest winds. The Pisco region is also home to local dust storms called Paracas, which transport dust material to the continental shelf as a response to seasonal erosion and transport events in the Ica Desert (∼ 15 • S). This process reflects atmospheric stability conditions and coastal sea surface temperature connections (Gay, 2005). In contrast, sediment fluvial discharge is more important on the northern coast of Peru, where there are large rivers, and it decreases southward, where arid conditions are dominant (Garreaud and Falvey, 2009;Scheidegger and Krissek, 1982). This discharged material is redistributed southward by coastal currents along the continental shelf (Montes et al., 2010;Smith, 1983). In addition, small rivers exist in our study area, such as the Pisco River, which can increase their flow during strong El Niño events (Bekaddour et al., 2014). It has also been demonstrated that, during El Niño events and coincident positive PDO, there is an increase in precipitation along northern Peru and, consequently, higher river discharge, mainly from the large rivers (e.g., the Santa River), whereas an opposite behavior is observed during La Niña events and the negative phase of PDO (Bekaddour et al., 2014;Böning and Brumsack, 2004;Lavado Casimiro et al., 2012;Ortlieb, 2000;Rein, 2005Rein, , 2007Scheidegger and Krissek, 1982;Sears, 1954).
A significant number of studies have described the climatic, hydrologic and oceanographic changes during the last 1000 years on the Peruvian continental shelf (Ehlert et al., 2015;Gutiérrez et al., 2011;Salvatteci et al., 2014b;Sifeddine et al., 2008). Evidence of changes in the Humboldt Current circulation system and in the precipitation pattern has been reported. Salvatteci et al. (2014b) show that the Medieval Climatic Anomaly (MCA) exhibits two distinct patterns of Peruvian upwelling characterized by weak/intense marine productivity and sub-surface oxygenation, respectively, as a response to the intensity of South Pacific Sub-tropical High (SPSH) linked to the Walker circulation. During the Little Ice Age (LIA), an increased sediment discharge over the Pisco continental shelf was described, as well as a stronger oxygenation and lower productivity (Gutiérrez et al., 2009;Salvatteci et al., 2014b;Sifeddine et al., 2008). In addition, during the Current Warm Period (CWP), the Peruvian Upwelling Ecosystem exhibited (1) an intense oxygen minimum zone (OMZ) and an increase in marine productivity, (2) a significant sea surface temperature cooling (∼0.3-0.4 • C decade −1 ), and (3) an increase in terrigenous material input (Gutiérrez et al., 2011).
Here we present new data regarding the effective mode of transport of mineral fractions to the Pisco shelf during the last millennium, confirming previous work and bringing new knowledge about the climatic mechanism behind Humboldt circulation and atmospheric changes, especially during the MCA. Our results identify wind intensification during the second part of the MCA and CWP, in contrast to a decrease in the wind intensity during the LIA and the first part of the MCA synchronous with fluvial discharge increases. Comparisons with other paleoclimate records indicate that the ITCZ displacement, the SPSH and the Walker circulation were the main drivers for the hydroclimate changes along the coastal Peruvian shelf during the last millennium. Reinhardt et al. (2002), Suess et al. (1987) and Gutiérrez et al. (2006) described the sedimentary facies in the Peruvian shelf and the role of currents in the erosion process as well as the redistribution and favorable hemipelagic sedimentation of material over the continental shelf. These studies showed that high-resolution sediment records are present in specific localities of the Peruvian continental margin. Suess et al. (1987) described the two sedimentary characteristic facies between 6 and 10 • S and between 11 and 16 • S. The first one, 6-10 • S (Salaverry Basin), is characterized by the absence of hemipelagic sediment accumulation, because in this zone the southward poleward undercurrent is strong. The second one, Lima Basin (11-16 • S), is characterized by a lens-shaped depositional center of organic-rich mud facies favored by oceanographic dynamics from the position and low velocity of the southward poleward current on the continental shelf (Reinhardt et al., 2002;Suess et al., 1987). High-resolution sediment echo sounder profiles further characterize the mud lens nature and complement the continental shelf information (Salvatteci et al., 2014a). These upper mud lenses are characterized by fine grain size, a diatomaceous, hemiplegic mud with high organic carbon, and the absence of erosive and bioturbation processes.

Sedimentary settings
The Pisco continental shelf sediments are a composite of laminated structures characterized by an array of more or less dense sections of dark and light millimetric laminae (Brodie and Kemp, 1994;Salvatteci et al., 2014a;Sifeddine et al.,  2008). The laminae structure and composition result from a complex interplay of factors including the terrigenous material input (both aeolian and fluvial), the upwelling productivity, and associated particle export to the seafloor (Brodie and Kemp, 1994;Salvatteci et al., 2014a). The anoxic conditions favored by an intense OMZ (Gutiérrez et al., 2006) and weak current activity in some areas (Reinhardt et al., 2002;Suess et al., 1987) favor the preservation of paleoenvironmental signals and consequently a successful recording of the environmental and climate variability.
Along the Peruvian coast, lithogenic fluvial material is supplied by a series of large rivers that are more significant to the north of the study area (Lavado Casimiro et al., 2012;McPhillips et al., 2013;Morera et al., 2011;Rein, 2005;Scheidegger and Krissek, 1982;Unkel et al., 2007). In fact, Smith (1983) concluded that sedimentary material can be transported for long distances in an opposite direction of prevailing winds and surface currents in upwelling zones. In fact, the coastal circulation off Peru is dominated by the poleward Peru-Chile undercurrent (PCUC), which flows over the outer continental shelf and upper continental slope, whereas the equatorward Peru coastal current is limited to a few dozens of meters in the surface layer (Chaigneau et al., 2013). On the other hand, several works have shown that precipitation, fluvial input discharge (Bekaddour et al., 2014;Bendix et al., 2002;Lavado-Casimiro and Espinoza, 2014), and the PCUC increase during the El Niño events (Hill et al., 1998;Strub et al., 1998;Suess et al., 1987). These observations suggest a potential for the fluvial particles to spread over the continental margin under wet paleoclimatic conditions (e.g., El Niño or El Niño-like). Lithogenic material in the study area might also originate from wind-driven dust storms or vientos Paracas, which are more frequent and intense during austral winters (Escobar Baccaro, 1993;Gay, 2005;Haney and Grolier, 1991) and by the saltation and suspension mechanisms with which this material reaches the continental shelf.

Stacked record
The B040506 (hereafter "B06"; 14 • 07.90 S, 76 • 30.10 W; 299 m water depth) and the G10-GC-01 (hereafter "G10"; 14 during the Galathea-3 cruise, respectively (Fig. 1a). We compared the age models and performed a laminae crosscorrelation between the two cores in order to develop a continuous record for the last millennium (Salvatteci et al., 2014a) (Fig. S1 in the Supplement). The choice of these two cores was based on previous detailed stratigraphic investigations and available complementary multi-proxy reconstructions (Gutiérrez et al., 2006(Gutiérrez et al., , 2009Salvatteci et al., 2012Salvatteci et al., , 2014aSifeddine et al., 2008). The box core B06 (0.75 m length) is a laminated core with a visible slump at ∼ 52 cm and three thick homogeneous deposits (1.5 to 5.0 cm thick) identified in the SCOPIX images. These intervals were not considered in our study (Fig. S1). The presence of filaments of the giant sulfur bacteria Thioploca spp. in the top of core B06 confirms the successful recovery of the sediment water interface.
According to the biogeochemical analysis in Gutiérrez et al. (2009) (i.e., palynofacies, oxygen index (Rock-Eval), total organic carbon and δ 13 C), B06 is characterized by a distinctive shift at ∼ 30 cm, more details are provided by Sifeddine et al. (2008) and Salvatteci et al. (2014a). The age model of B06 was inferred from five 14 C-calibrated accelerator mass spectrometry (AMS) age distributions (Fig. S1), showing that this core covers the last ∼ 700 years. For the last century, which is recorded only by B06, the age model was based on downcore natural excess 210 Pb and 137 Cs distributions and supported by bomb-derived 241 Am distributions ( Fig. S2 and Gutiérrez et al., 2009). The mass accumulation rate after ca. AD 1950 was 0.036 ± 0.001 g cm −2 yr −1 and before ca. AD 1820 was 0.022 ± 0.001 g cm −2 yr −1 . On the other hand, G10 is a gravity-laminated sediment core of 5.22 m presenting six units and exhibiting some minor slumping. The G10 age model was based on 31 samples of 14 C-calibrated AMS age distributions, showing that the core covers the Holocene period (Salvatteci et al., 2014b. Here we used only a laminated section between ∼ 18 and 45 cm that chronologically covered part of the MCA period (from ∼ AD 1050 to 1500) and presented no slumps (Fig. S1).
The spatial regularity of the initial core sampling combined with the naturally variable sedimentation rate implied variable time rates between samples (150 samples in total). Each sample is 0.5 cm thick in B06 and usually includes 1-2 laminae. On the other hand, in core G10, each sample is 1 cm thick, including 3-4 laminae. The results considering the sedimentation rates showed that the intervals during MCA, LIA and CWP span 18, 7, and 3 years, respectively. Because of differences in the subsampling thickness between cores and variable sedimentation rates, results are binned by 20-year intervals (the lowest time resolution among samples) after linear interpolation and 20-year running mean of the original data set.

Grain size analyses
To isolate the mineral terrigenous fraction, organic matter, calcium carbonate and biogenic silica were successively removed from approximately 100 mg of bulk sediment sample using H 2 O 2 (30 % at 50 • C for 3 to 4 days), HCl (10 % for 12 h) and Na 2 CO 3 (1 M at 90 • C for 3 h), respectively. Between each chemical treatment, samples were repeatedly rinsed with deionized water and centrifuged at 4000 rpm until neutral pH. After pre-treatment, the grain size distribution was determined with an automated image analysis system (model FPIA3000, Malvern Instruments). This system is based on a CCD (charge-coupled device) camera that captures images of all of the particles homogeneously suspended in a dispersal solution by rotation (600 rpm) in a measurement cell. After magnification (×10), particle images are digitally processed and the equivalent spherical diameter (defined as the diameter of the spherical particle having the same surface as the measured particle) is determined. The optical magnification used (×10) allows the counting of particles with equivalent diameters between 0.5 and 200 µm. Prior to the FPIA analysis, all samples were sieved with a 200 µm mesh in order to recover coarser particles. Since particles > 200 µm were never found in any samples, the grain size distribution obtained by the FPIA method reliably represents the full particle size range in the sediment. A statistically significant number of particles (in the hundreds of thousands, up to 300 000) are automatically analyzed by FPIA, providing particle size information comparable to that obtained with a laser granulometer along with images of the individual particles. Using the images to check the efficiency of the pretreatments, we ensured that both organic matter and biogenic silica had been completely removed from all the samples. Finally, particle counts were binned into 45 different size bins between 0.5 and 200 µm instead of the 225 set by the FPIA manufacturer in order to reduce errors related to the presence of very few particles in some of the preselected narrow bins. Grain size distributions are expressed as (%) volume distributions.

Determining sedimentary components and the de-convolution fitting model
As different particle transport/deposition processes are known to influence the grain size distribution of the lithic fraction of sediment (e.g., Holz et al., 2007;Pichevin et al., 2005;Prins et al., 2007;Stuut et al., 2005Stuut et al., , 2002Sun et al., 2002;Weltje andPrins, 2003, 2007;Weltje, 1997), identifying the individual components of the polymodal grain size distribution is decisive for paleoenvironmental reconstructions. The numerical characteristics (i.e., amplitude (A), geometric mean diameter (Gmd), and geometric standard deviation (Gsd) of the individual grain size populations whose combination forms the overall grain size distribution) were determined for all samples using the iterative least-squares  (1990). This fitting method aims to minimize the squared difference between the measured volume grain size distribution and the one computed from a mathematical expression based on lognormal function. The number of individual grain size populations to be used is determined by the operator, and all statistical parameters (A, Gmd and Gsd) are allowed to change from one sample to another. This process presents a strong advantage compared to end-member modeling (e.g., Weltje, 1997), in which the individual grain size distributions are maintained constant over the whole time series, the only fitting parameter being the relative amplitude, A. Indeed, it is unlikely that the parameters that govern both transport and deposition of lithogenic material, and therefore grain size of particles, remain constant over time. In turn, variations in these parameters are expected to induce changes in the grain size distribution parameters such as Gmd and Gsd.

Basis for interpretation
Both sediment cores (B06 and G10) exhibit a roughly bimodal grain size distribution presenting significant variation in amplitude and width. These modes correspond to finegrain-size classes from ∼ 3 to 15 µm and coarser grain size classes between ∼ 50 and 120 µm (Fig. S3). A principal component analysis (PCA) based on the Wentworth (1922) grain size classification identifies four modes that could explain the total variance of the data set (Fig. S4). The measured and computed grain size distributions show high correlations ranging from R 2 = 0.75 to 0.90, demonstrating that the use of four grain size modes is well adapted to our sediment samples and that the computed ones may be reliable for further interpretation (Fig. 2). Lower correlations only occurred for six samples that are characterized by small proportions of terrigenous material compared to biogenic silica, organic matter and carbonates. In these cases, the number of lithic particles remaining after chemical treatments was small, which increased the associated relative error. However, these samples have been included in the data set since they all presented a high contribution of coarser particles.
Grain size parameters are presented in Table 1. The first mode (M1), with a Gmd of approximately 3 ± 1 µm, and the second one (M2), with a Gmd of 10 ± 2 µm, are characterized by large Gsd (∼ 2σ ), indicating a low degree of sorting. Such a low degree of sorting suggests a slow and continuous depositional process as occurs in other environments (Sun et al., 2002). The coarsest modes, M3 and M4, display mean Gmd values of 54 ± 12 and 91 ± 13 µm, respectively. These modes present Gsd values close to 1σ . The Gmd values of the two coarsest modes are consistent with the optimal grain size transported under conditions favorable to soil erosion (lack of vegetation, low threshold friction velocity, surface roughness and low soil moisture) and low wind friction velocity Figure 2. Comparison between a measured grain size distribution and the fitted curve using lognormal function and its partitioning into four individual grain size modes. The measured data are a mean grain size distribution from all samples of B6 and G10 cores. (Iversen and White, 1982;Kok et al., 2012;Marticorena and Bergametti, 1995;Marticorena, 2014;Shao and Lu, 2000). Such conditions prevail in the studied area because central coastal Peru consists of a sand desert area characterized by the absence of rain, a lack of vegetation and persistent wind (Gay, 2005;Haney and Grolier, 1991).
In the vicinity of desert areas, where wind-blown transport prevails, particles with grain size as high as ∼ 100 µm can accumulate in marine sediments (e.g., Flores-Aqueveque et al., 2015;Stuut et al., 2007) or even in lacustrine sediments . Indeed, Stuut et al. (2007) reported the presence of distributions typical of wind-blown particles with ∼ 80 µm grain size (∼ 29 • S North Chile), which is consistent with our results. In the studied area, the emission and transport of mineral particles are related to the strong wind events called Paracas. Paracas dust emission is a local seasonal phenomenon that preferentially occurs in winter (July-September) and is due to an intensification of the local surface winds (Escobar Baccaro, 1993;Haney and Grolier, 1991;Schweigger, 1984). The pressure gradient of sea level between 15 and 20 • S, 75 • W is the controlling factor of Paracas winds (Quijano, 2013), along with local topography (Gay, 2005). Coarse particles found in continental sediments off Pisco cannot have a fluvial origin because substantial hydrodynamic energy is necessary to mobilize particles of this size (50-100 µm), and this region is devoid of large rivers (Reinhardt et al., 2002;Scheidegger and Krissek, 1982;Suess et al., 1987). Therefore, the coarsest modes (M3 and M4) can be interpreted as markers of aeolian transport resulting from surface winds and emission processes (Flores-Aqueveque et al., 2015;Hesse and McTainsh, 1999;Marticorena and Bergametti, 1995;McTainsh et al., 1997;Sun et al., 2002) and www.clim-past.net/12/787/2016/ Clim. Past, 12, 787-798, 2016 Table 1. Averaged parameters (geometric mean diameter (Gmd), amplitude (A) and geometric standard deviation (Gsd)) of the four lognormal modes (components) identified from measured size distributions of sediment samples (B6 and G10 cores).

M1 M2 M3 M4
Gmd (µm) A (%) Gsd Gmd (µm) A (%) Gsd Gmd (µm) A (%) Gsd Gmd (µm) A (%) Gsd 3 ± 1 16 ± 7 1.9 ± 0.2 10 ± 2 43 ± 15 1.9 ± 0.2 54 ± 12 20 ± 10 1.4 ± 0.2 90 ± 13 20 ± 13 1.2 ± 0.2 indicate a local and proximal aeolian source (i.e., Paracas winds). This interpretation is in contrast to the Atacama Desert source suggested by Ehlert et al. (2015) and Molina-Cruz (1977). Ehlert et al. (2015), who used the same sediment core (B06), and also indicated difficulties in the interpretation of the detritical Sr isotopic signatures as an indicator of the terrigenous sources. These difficulties can be associated with the variability in the 87 Sr / 86 Sr due to grain size (Meyer et al., 2011). The finest M1 component (∼ 3 µm) may be linked to both aeolian and fluvial transport mechanisms. Thus, because its origin is difficult to determine, and because its trend appears to be relatively independent of the other components, we do not use it further. The M2 component (∼ 10 µm) is interpreted as an indicator of fluvial transport (Koopmann, 1981;McCave et al., 1995;Stuut and Lamy, 2004;Stuut et al., 2002Stuut et al., , 2007. Indeed, this is consistent with the report by Stuut et al. (2007) for the fluvial mud (∼ 8 µm) in the south of Chile (> 37 • S), where the terrigenous input is dominated by fluvial origins. A fluvial origin of this M2 component is also supported by showing the same trend in the geochemical proxies, such as radiogenic isotope compositions of detrital components (Ehlert et al., 2015), mineral fluxes (Sifeddine et al., 2008) or %Ti (Salvatteci et al., 2014b), indicating more terrigenous transport during the LIA, when humid conditions were dominant. The M2 component is interpreted as being linked to river material discharge, mostly from the north Peruvian coast, and redistribution by the PCUC and bottom currents (Montes et al., 2010;Rein et al., 2004;Scheidegger and Krissek, 1982;Unkel et al., 2007).

Aeolian and fluvial input variability during the past ∼ 1000 years
Grain size component (M2, M3 and M4; Table 1) variations in the composite records (B06 and G10) express changes in wind stress and fluvial runoff at multidecadal to centennial scales during the last millennium. The sediments deposited during the MCA exhibit two contrasting patterns of grain size distributions. A first sequence dated from AD 1050 to 1170 has low values of D50 (i.e., median grain size) that vary around 16 ± 6 µm and are explained by 50 ± 14 M2, 16 ± 8 M3, 21 ± 5 M4 and 13 ± 5 % M1 contributions. A second sequence, dated from AD 1170 to 1450, was marked by high values of D50 in the range of 34 ± 18 µm, with average contributions of 36 ± 8 for M2, 21 ± 10 for M3, 29 ± 15 for M4 and 14 ± 6 % for M1 (Table 2). These results indicate high variability in transport of particles during the MCA, with more fluvial sediment discharge from 1050 to 1170, followed by an aeolian transport increase between AD 1170 and 1450 (Fig. 3).
During the LIA (AD 1450-1800), the deposited particles were dominated by fine grain sizes with a D50 varying around an average of 15 µm, explained by 53 ± 15 % M2 contribution. In contrast, the contribution of M3 averaged 19 ± 9 % and ranged from 4 to 45 %, whereas M4 showed an average contribution of 14 ± 11 % and varied from 0 to 44 % during the same period. The dominant contribution of the finest-sized particles of M2 suggests a high fluvial terrigenous input to the Peruvian continental shelf. It is important to note that M2 contributions increased from the beginning to the end of the LIA at ∼ AD 1800, suggesting a gradual increase in fluvial sediment discharge input related to the enhancement of the continental precipitation (Fig. 3c). Indeed, during the LIA, our results confirm previous interpretations of wet conditions along the Peruvian coast (Gutiérrez et al., 2009;Salvatteci et al., 2014b;Sifeddine et al., 2008). These results also imply that this period was characterized by weak surface winds and hence a weaker coastal upwelling.
Subsequently, D50 variations show multidecadal variability during the last ∼ 200 years that is divided into three distinctive periods. The first one, from ∼ AD 1800 to 1850, shows dominance of coarse particles around 50 µm, explained by the high contribution of M3 and M4 (up to 45 and 50 %, respectively) during this period. These results suggest a period of drier climate and very strong wind conditions. The second one, from AD 1850 to 1900, displays values around ∼ 20 µm explained by ∼ 40 of M2, ∼ 20 of M3 and ∼ 20 % of M4, suggesting that fluvial sediment discharge was the dominant transport mechanism, although not as significant as during the LIA. The third period spans from AD 1900 to the final part of record and covers the CWP. Our results reveal a dominance of coarse particles during the most of this period (D50 up to of 80 µm) that arise from high contributions of M3 and M4 (∼ 40 and ∼ 50 %, respectively). However, a clear decrease in the D50 is displayed at the end of this period that is explained by a decrease in contributions of the aeolian component M4 (∼ 20 %), although the contribution of M3 and M2 remain relatively stable (∼ 25 and 30 %, respectively). These conditions display no clear dominance of a given transport mode during this time. In addition, markedly coarser particles in the M4 component were very Clim. Past, 12, 787-798, 2016 www.clim-past.net/12/787/2016/  16 ± 8 6-28 21 ± 10 0-39 19 ± 9 4-45 23 ± 10 6-44 M4 21 ± 5 12-30 29 ± 15 10-55 14 ± 11 0-44 25 ± 13 0-56 common during this time (the last 200 years), indicating a strong probability of extreme wind stress events (Fig. 3f).

Climatic interpretations
Our findings suggest a combination of regional and local atmospheric circulation mechanism changes that controlled the pattern of sedimentation in the study region. Our record is located under the contemporary seasonal Paracas dust storm path, but it also records discharged fluvial muds that are supplied by the rivers along the Peruvian coast. Hence, this record is particularly well suited for a reconstruction of continental runoff/wind intensity in the central Peruvian continental shelf during the last millennium. The interpretation of the changes in the single records of the components (M2, M3 and M4) and their associations (e.g., ratios) can reflect pale-oclimatic variations in response to changes in atmospheric conditions. Here, we used the ratio between the aeolian components, defined as the contribution of the stronger winds over total wind variability: M4 / (M3 + M4). We consider this ratio to be a proxy of the local wind surface intensity and thus of the SPSH atmospheric circulation (Fig. 4a). Previous studies have similarly and successfully used grain size fraction ratios as paleoclimate proxies of atmospheric conditions and circulation to explain other sediment records (Holz et al., 2007;Huang et al., 2011;Prins, 1999;Shao et al., 2011;Stuut et al., 2002;Sun et al., 2002;Weltje and Prins, 2003). As explained above, the MCA was characterized by a sinelike peak structure that depicts two different climate stages. During the first stage, spanning from ∼ AD 1050 to 1170, the fluvial input show a peak centered at AD 1120 that was linked to a precipitation increase accompanied by a decrease www.clim-past.net/12/787/2016/ Clim. Past, 12, 787-798, 2016  (M2) anomaly reconstruction on the continental shelf, and records of (c) terrigenous flux (total minerals) in Pisco continental shelf by Sifeddine et al. (2008), (d) OMZ activity (Re / Mo anomalies) negative values indicate more anoxic conditions (the axis was reversed) (Salvatteci et al., 2014b), (e) ITCZ migration (%Ti) (Peterson and Haug, 2006), (f) SAMS activity reconstruction (δ 18 O Palestina Cave) (Apaéstegui et al., 2014), (g) eastern temperatures reconstruction (Rustic et al., 2015) and (h) Indo-Pacific temperatures reconstruction (Oppo et al., 2009). in wind intensity. Those results suggest a southward ITCZ displacement (Fig. 4e) as a response to more El Niño-like conditions as suggested by Rustic et al. (2015) (Fig. 4g and h). In contrast, during the second stage the surface winds had their greatest intensity with a peak centered at ∼ AD 1200 as a consequence of displacement of the ITCZ-SPSH system. The displacement of the SPSH core towards the eastern South American coast intensified alongshore winds as a regional response to stronger Walker circulation. These features are in agreement with the ocean thermostat mechanism proposed by Clement et al. (1996). This mechanism produces a shallow thermocline in the eastern Pacific ( Fig. 4g and h) and consequently more intense upwelling conditions and a stronger OMZ offshore of Pisco recorded in low values of the Re / Mo ratio (Fig. 4d). These two patterns (i.e., enhanced fluvial transport/enhanced wind intensity) might have been triggered by the expression of Pacific variability at multidecadal timescales with the combined action of the Atlantic Multidecadal Oscillation (AMO). Indeed, other works provide evidence during the MCA for low South American monsoon system (SAMS) activity at multidecadal timescales driven by the AMO (Fig. 4f) (Apaéstegui et al., 2014;Bird et al., 2011;Reuter et al., 2009). Thus, besides the displacement of the ITCZ, the AMO could have modulated Walker circulation at a multidecadal variability during the MCA through mechanisms such as those described by McGregor et al. (2014) and Timmermann et al. (2007).
Our results, combined with other paleo-reconstructions, suggest that the LIA was accompanied by a weakening of the regional atmospheric circulation and of the upwelling favorable winds. During the LIA, the mean climate state was controlled by a gradual intensification of the fluvial input of sediments to the continental shelf, thus indicating more El Niño-like conditions (Fig. 4b). These features are confirmed by an increase in the terrigenous sediment flux, as described by Sifeddine et al. (2008) (Fig. 4c) and Gutiérrez et al. (2009) and by changes of the radiogenic isotopic composition of the terrigenous fraction (Ehlert et al., 2015). These wet conditions are also marked by an intensification of the SAMS and the southern meridional displacement of the ITCZ, as evidenced by paleo-precipitation records in the Andes and in the Cariaco Trench (Apaéstegui et al., 2014;Haug et al., 2001;Peterson and Haug, 2006) (Fig. 4e). At the same time, a prevalence of weak surface winds (Fig. 4a) and an increase in subsurface oxygenation driving sub-oxic conditions in the surface sediment are recorded (Fig. 4d). These characteristics also support the hypothesis of the ITCZ-SPSH southern meridional displacement and are consistent with a weakening of the Walker circulation (Fig. 4g).
The transition period between the LIA and CWP appears as an abrupt event showing a progressive positive anomaly in the wind intensity synchronous with a rapid decrease in fluvial input to the continental shelf ( Fig. 4a and b). This transition suggests a rapid change of meridional (ITCZ-SPSH) and zonal (Walker) circulation interconnection, which controls the input of terrigenous material (fluvial/aeolian). Gutiérrez et al. (2009) found evidence of a large reorganization in the tropical Pacific climate with immediate effects on ocean biogeochemical cycling and ecosystem structure at the transition between the LIA and CWP. The increase in the regional wind circulation that favors aeolian erosive processes simultaneously leads to an increase in the OMZ intensity related to upwelling intensification.
Finally, during the CWP (∼ AD 1900 to present), a trend to steadying of low fluvial input (Fig. 4b) was combined with an increase in wind intensity (Fig. 4a) that was coupled to a strong OMZ. This setting suggests the northernmost ITCZ-SPSH system position. This hypothesis is supported by other studies on the continental shelf of Peru (Salvatteci et al., 2014b) and also in the eastern Andes, where a decrease in rainfall of between ∼ 10 and 20 % relative to the LIA was reported for the last century (Reuter et al., 2009). Enhancement of wind intensity is also consistent with the multidecadal coastal cooling and increase in upwelling productivity since the late nineteenth century (Gutiérrez et al., 2011;Salvatteci et al., 2014b;Sifeddine et al., 2008) and confirms the relations between the intensification of the upwelling activity induced by the variability in the regional wind intensity from SPSH displacement.
The increase in the wind intensity over the past two centuries likely represents a result of the modern positioning of the ITCZ-SPSH system and the associated intensification of the local and regional winds (Fig. 4a). The contributions of aeolian deposition material ( Fig. 3e and f) and, as a consequence, the wind intensity and its variability during the last 100 years are stronger than during the second sequence of the MCA (Fig. 4a) under similar conditions (i.e., position of the ITCZ-SPSH system). This variability implies a forcing mechanism in addition to the enhancement of the wind intensity, one that may be related to the current climate change conditions (Bakun, 1990;England et al., 2014;Sydeman et al., 2014). Moreover, during the CWP, the wind intensity showed a direct relation with OMZ strength (Fig. 4a and d) that suggests an increase in the zonal gradient and thus in the Walker circulation on a multidecadal scale.
Our record shows that on a centennial scale, the fluvial input changes are driven by the meridional ITCZ position and a weak gradient of the Walker circulation, consistent with El Niño-like conditions. In contrast, variations in the surface wind intensity are linked to the position of the SPSH modulated by both the meridional variation in the ITCZ and the intensification of the zonal gradient temperature related to the Walker circulation and expressing La Niña-like conditions. A clear relation between the zonal circulation and wind intensity at a centennial timescale is displayed. All these features modulate the biogeochemical behavior of the Peruvian upwelling system.

Conclusions
Study of the grain size distribution in laminated sediments from the Pisco Peruvian shelf has allowed the reconstruction of changes in wind intensity and terrigenous fluvial input at centennial and multidecadal timescales during the last millennium. The long-term variation in the M2 (∼ 10 µm) mode is an indicator of hemipelagic fluvial input related to the regional precipitation variability. At the same time, the M3 (54 ± 11 µm) and M4 (91 ± 11 µm) components are related to aeolian transport and thus with both local and regional wind intensity. The temporal variations in these fractions indicate that the MCA and CWP periods were characterized by an increment in the coarse-particle transport (M3 and M4) and thus an enhancement of the surface wind intensity, whereas the LIA was characterized by stronger fluvial input as evidenced from an increase in fine (M2) particles. Comparison between records reveals a coherent match between the meridional displacement of the ITCZ-SPSH system and the regional fluvial and aeolian terrigenous input variability. The ITCZ-SPSH system northern displacement during the second period of the MCA and the CWP was associated with the intensification of the Walker cell and La Niña Like conditions, resulting in stronger winds, upwelling-favorable conditions, enhanced marine productivity and greater oxygen depletion in the water column. In contrast, the southward migrations of the ITCZ-SPSH system during the LIA correspond to an enhancement to the South American monsoon circulation and El Niño-like conditions, driving the increase in the precipitation and the terrigenous fluvial input to the Pisco continental shelf, lower productivity and increased oxygenation. Two patterns observed during the MCA, respectively marked by fluvial intensification and wind intensification, could have been forced by Pacific Ocean variability at multidecadal timescales. Further studies of the paleowind reconstruction at high time resolution, combined with model simulation, are needed to better understand the interplay between the Pacific and Atlantic Ocean connection on climate variability as evidenced by McGregor et al. (2014) in the modern Pacific climate pattern.
The Supplement related to this article is available online at doi:10.5194/cp-12-787-2016-supplement.