An ~15,000-Year Record of El Niño-Driven Alluviation in Southwestern Ecuador

See allHide authors and affiliations

Science  22 Jan 1999:
Vol. 283, Issue 5401, pp. 516-520
DOI: 10.1126/science.283.5401.516


Debris flows have deposited inorganic laminae in an alpine lake that is 75 kilometers east of the Pacific Ocean, in Ecuador. These storm-induced events were dated by radiocarbon, and the age of laminae that are less than 200 years old matches the historic record of El Niño events. From about 15,000 to about 7000 calendar years before the present, the periodicity of clastic deposition is greater than or equal to 15 years; thereafter, there is a progressive increase in frequency to periodicities of 2 to 8.5 years. This is the modern El Niño periodicity, which was established about 5000 calendar years before the present. This may reflect the onset of a steeper zonal sea surface temperature gradient, which was driven by enhanced trade winds.

The dramatic effects of the 1997–98 El Niño event have highlighted several shortcomings in our understanding of the El Niño–Southern Oscillation (ENSO) phenomenon (1). These shortcomings include the age of onset of modern ENSO variability, long-term (>103 years) changes in the frequency of past extreme El Niños and their relation with varying oceanic and atmospheric states, and the frequency and magnitude of the ENSO in Earth's greenhouse future. High-resolution records of prehistoric El Niños are needed to address these questions.

Proxy records of prehistoric El Niños have been obtained from a variety of archives, including corals, ice cores, tree rings, flood deposits, beach ridges, archeological middens, and soils (2,3). However, high-resolution coral and ice core records (4) have been limited to the past two millennia, and longer proxy records of El Niños from the tropical Pacific region are inherently discontinuous (3, 5). Here, we present a high-resolution record of storm-derived clastic sedimentation that spans the past 15,000 years and appears to record El Niño events.

Sea surface temperatures (SSTs) near Guayaquil, Ecuador, are some of the first to warm in the region of upwelling along the coasts of Perú and Ecuador during the onset of an El Niño event (1), and typically, rainfall in this region greatly increases over background levels during the onset of strong to very strong (hereafter, severe) (6) El Niños. The mature phase of El Niño, in which the zonal SST gradient is at a minimum, occurs just as the normal austral summer rainy season in western Ecuador begins. During stronger El Niños, the Walker Circulation is altered so that rising motions occur in parts of the eastern Pacific troposphere that are normally characterized by subsidence and inversion. Under these circumstances, the coastal regions of Ecuador and Perú experience intraseasonal bursts of deep convection and torrential rainfall. For example, during the 1982–83 El Niño, the thermocline began to deepen off the south coast of Ecuador in October 1982, and ocean temperatures abruptly increased 6.2°C at the surface and 9°C at a depth of 50 m (7). Convection-driven precipitation was substantially increased between 2.5° and 6°S at sites that are as far as several hundred kilometers inland and was increased at elevations from sea level to >3000 m above sea level (masl) (8, 9). Beginning in January 1983, repetitive bursts of extreme rainfall were responsible for much of the devastating flooding that occurred in the region (9). During normal and La Niña years, the Walker Circulation returns to normal, and conditions required for deep convection seldom exist (10).

We analyzed a 9.2-m-long core (obtained in June 1993) that spans the past ∼15,000 calendar years from Laguna Pallcacocha (4060 masl and 2°46′S, 79°14′W) (Fig. 1). The Holocene section contains hundreds of light-colored inorganic, clastic laminae (<0.1 to 1.0 cm thick) that are interbedded with massive organic-rich laminae (Fig. 2). This stratigraphy is reflected in variable bulk density, magnetic susceptibility, C content, and color (11) (Fig. 3, A through C). Sediment records from most similar deglacial lakes in the tropical Andes reveal an abrupt rise in organic C, which correlates with deglaciation, and a sedimentologically uneventful Holocene interval (12). The 0.05-km2 lake is located ∼500 m east of the continental divide in a cirque basin that has been ice free for the past ∼14,000 calendar years. The volcanic bedrock of the region erodes easily, and debris flows and talus are abundant along the divide and are located within several hundred meters of the lake.

Figure 1

Location of Laguna Pallcacocha is 500 m east of the continental divide in southwestern Ecuador.

Figure 2

Photographs of selected sections of the Pallcacocha core. Depth scale is meters below the top of the recovered core. The transition from light-colored inorganic sediment, which comprises the Late Glacial interval, to dark organic-rich gyttja occurs ∼12,300 cal yr B.P. in the base of drive 10. A progressive increase in the frequency of clastic laminae that are light colored and inorganic occurs throughout the Holocene (for example, drives 10 to 1).

Figure 3

Downcore variation in (A) bulk density, magnetic susceptibility, and organic C content; (B)14C ages; and (C) gray scale. The 14 AMS14C ages are from terrestrial plant macrofossils and have been converted to the calendar year time scale with the calibration program of Stuiver and Reimer (27). The ages with an asterisk are from a nearby lake and are correlated to this core by their stratigraphic position in relation to geochemically distinct tephra, which are denoted as T in the left-hand depth scale. Short, stippled horizontal lines in the depth scale denote core section boundaries. The gray scale ranges from 0 (white) to 256 (black) and is primarily controlled by the organic C content of the sediment (14). The age (14) of clastic laminae that are light-colored, inorganic, and <200 years old agrees with the timing of instrumentally or historically documented occurrences (or both) of severe or moderate El Niños (or both) (solid and open boxes, respectively, in the right-hand depth scale) (6), thus establishing ENSO as a major control on sedimentation in Laguna Pallcacocha.

The inorganic laminae likely represent deposition from density-driven undercurrents that carried terrestrially derived clastic sediment and some organic matter (13) into the 15-m-deep basin in response to brief storm events that mobilized the abundant loose sediment in the headwaters of the drainage basin as debris flows. The inorganic laminae have <2% organic C, whereas the organic-rich laminae have >10% organic C; the inorganic laminae tend to have a coarser modal grain size and contain less biogenic silica than organic-rich laminae (1 to 3% versus 5 to 10%, respectively). Most of the inorganic layers have abrupt basal contacts and fine upward.

The core was dated on the basis of 14 accelerator mass spectrometer (AMS) 14C ages of terrestrial macrofossils and the identification of seven tephra layers (Fig. 3B). We developed an age model (14) by assuming that the sediment-water interface is modern (1993 A.D.) (Fig. 4) and that, between AMS 14C-dated intervals (Fig. 3B), the rate of background C deposition was nearly constant and was interrupted by abrupt clastic depositional events, which diluted the C concentration of the sediment. Linear interpolation between dated intervals would result in a substantial underestimation of the time separating the rapidly deposited clastic events. Our constant C flux model (14) essentially removes the clastic laminae from the record, assigns time to the core, and then reinserts the clastic laminae in their correct temporal position.

Figure 4

(below).Age-depth plot for Laguna Pallcacocha. Ages are based on AMS14C dating of terrestrial macrofossils; symbol size includes a ±1σ uncertainty in the calendar year calibration (27). The sediment-water interface is assumed to be modern (1993 A.D.). The four solid dots that are plotted along the curve are the 14C ages (denoted with an asterisk in Fig. 3B) that are correlated to this core from a nearby lake on the basis of their stratigraphic position in relation to geochemically distinct tephra (Fig. 3A). Because of a ±10-cm uncertainty in the position that these14C ages correspond with in the Pallcacocha core, we use them only as supporting ages for the chronology of the core, which is based on 9 of the 10 14C ages that were obtained from Laguna Pallcacocha; because of nonunique calendar year solutions to the youngest 14C date (100 ± 60 14C yr B.P.) (Fig. 3B), we do not use this date in the age model (14). The progressive increase in sedimentation rates through the early Holocene is due to the increase in input of storm-induced clastic sediment, which reflects the progressive increase in frequency of El Niños through the middle and late Holocene.

Analysis of the gray-scale record of the past 660 calendar years (Fig. 2) confirms that the clastic layers reflect El Niño events (15) and depositional events with longer periodicities. The Blackman-Tukey spectrum of the gray-scale record (16) (Fig. 5) reveals concentrations of variance at 1/50 to 1/11, at 1/6.8, and at 1/5 cycles/year; singular spectrum analysis identifies oscillatory pairs with frequencies of 1/25 and 1/10 cycles/year; multitaper and maximum entropy methods (17) also indicate ENSO and decadal band variance. We used the multitaper method (17, 18) to test the null hypothesis that the peaks are red noise; for the ENSO band, we rejected the hypothesis at the 99% level, but we did not reject the red noise hypothesis for the decadal frequency variance.

Figure 5

Results of time series analysis (16) of the gray-scale record (Fig. 3C); age ranges are in calendar years before the present. Because of breaks in the core (Fig. 3A), we performed the time series analysis on individual sections of the core and treated these as discrete time series. There is a clear spectral evolution from the Late Glacial through the Holocene with the ENSO band (19), which progressively achieves its modern strength by ∼5000 cal yr B.P.

The sediment record from 1800 to 1976 A.D. reveals a close match (Fig. 3C) between the timing of clastic laminae (low gray-scale value) and moderate to severe El Niño events (6). We estimated that the age uncertainty for this part of the record is ≥5% of the age of laminae. Of the 17 severe El Niños that occurred in this time period, 11 correlate within 2 years of major clastic laminae, and 1 is within 3 years of a lamina; we define major clastic laminae as having gray-scale values that are below the mean for the period from 1800 to 1976. The other five severe El Niños of this interval occurred within 2 years of relatively minor clastic laminae. The interval from 1976 to 1993 was disturbed during coring, and thus, we do not have a signal of the 1982–83 event. The eight severe El Niños of the past 100 years (6) correlate precisely with clastic laminae in the core (Fig. 3C). Twenty-six of the 33 moderate El Niños that occurred from 1800 to 1976 A.D. (6) correlate precisely with at least minor clastic layers, as indicated by low or falling gray scale. The moderate El Niños of 1806–07, 1821–24, 1860, 1897, 1904–05, and 1907 correlate with major clastic laminae. Five of these six moderate events follow relatively minor sedimentologic responses to severe events, and thus, it is possible that these moderate El Niños had an anomalously large impact on sediment delivery to Laguna Pallcacocha because of an unusual abundance of fine-grained material left in debris flow channels or along the margin of the lake basin during the preceding severe events.

The full sediment record shows that the frequency of clastic depositional events, which is apparent in the visible stratigraphy (Fig. 2) and time series analysis (Fig. 5), has increased progressively. The periodicity of clastic sedimentation has evolved from ≥15 years during the Late Glacial and early Holocene [∼15,000 to 7000 calendar years before the present (cal yr B.P.)] (Fig. 5). Beginning at ∼7000 cal yr B.P., clastic events were spaced 10 to 20 and 2 to 8.5 years apart; the 2- to 8.5-year periodicity is most apparent after ∼5000 cal yr B.P. and is consistent with the periodicity of modern El Niños (19). This increase in clastic input is reflected in average bulk sedimentation rates (Fig. 4). During the Late Glacial (∼10,000 to 15,000 cal yr B.P.), sedimentation rates averaged 5.2 cm/(100 years). From ∼10,000 to 7000 cal yr B.P., average bulk sedimentation rates dropped to their lowest values of the past 15,000 years [2.7 cm/(100 years)] because the input of glacially derived sediment had ceased and storm–derived sediment input was low. By ∼7000 cal yr B.P., sedimentation rates began to progressively increase and reached 5.5 cm/(100 years) from ∼2400 to 1200 cal yr B.P.

The lack of variance in the ENSO band from ∼15,000 to 7000 cal yr B.P. implies that the ENSO was weak, perhaps because the zonal SST gradient was subdued. Because the tropical Pacific behaves as a coupled ocean-atmosphere system (1), a reduction in trade wind circulation is simultaneously the cause of and the consequence of a reduced zonal SST gradient. A weakened zonal SST gradient may have been caused by a smaller and less intense western Pacific warm pool, by warmer SSTs in the equatorial and coastal upwelling zones of the eastern Pacific, or by both.

Recent Sr/Ca ratios and ∂18O records from the Great Barrier Reef, Australia, indicate that, during the early Holocene, the western Pacific SST was ∼1.2°C higher than the present SST (20). Thus, if the apparent reduction in ENSO during the early Holocene was caused by a reduced zonal SST gradient, the reduced gradient must have been caused by elevated SSTs in the eastern Pacific Ocean. Holocene molluscan assemblages in natural deposits and shell middens that were 14C-dated at >5000 years B.P. (∼5700 cal yr B.P.) along the Peruvian coast north of 10°S and along the southern coast of Ecuador are anomalous in that they are dominated by tropical species rather than the temperate species that have inhabited these coastal waters during non–El Niño years for the past 500014C years (5). The occurrence of these thermally anomalous molluscan assemblages (TAMAs) has been used to infer that the zonal Walker Circulation that is characteristic of normal or La Niña years was reduced during the early Holocene and that the oscillation between El Niño and La Niña states was muted (5, 21). This result is consistent with pollen records from Australasia, which record a middle Holocene onset of the highly variable climate that is produced by ENSO in that region today (22). A warmer eastern Pacific during the early Holocene is also consistent with the ∂18O record of tropical precipitation obtained from an ice core at Nevado Huascarán in the western cordillera of the Peruvian Andes (9°S and 6048 masl) (23).

Recent modeling results of the coupled ocean-atmosphere system over the tropical Pacific corroborate the occurrence of two oceanic modes in the tropical Pacific (24): one is warm and steady, and the other is cold and oscillating. The model results indicate that development of the cold oscillating mode is dependent on the achievement of a strong temperature gradient between surface water and the main thermocline of the tropical east-central Pacific. Because the temperature of the main thermocline in the east-central Pacific is driven by higher latitude SSTs (25), we hypothesize that the temperature gradient may have increased since 10,000 cal yr B.P. as a consequence of the steadily increasing contrast in seasonal insolation between the equator and high latitudes (26). This model (24) implies that if high-latitude SSTs warm more than tropical SSTs in an enhanced-greenhouse future, ENSO would be reduced with the development of a warm and steady tropical Pacific. Insofar as the early Holocene can be used as an analog for Earth's greenhouse future, our results and those of others (5, 22) lend credence to this scenario.

  • * To whom correspondence should be addressed. E-mail: rodbelld{at}


View Abstract

Stay Connected to Science

Navigate This Article