Millennial and orbital variations of El Nino/Southern Oscillation and high-latitude climate in the last glacial period

The El Nino/Southern Oscillation (ENSO) phenomenon is believed to have operated continuously over the last glacial interglacial cycle(1). ENSO variability has been suggested to be linked to millennial-scale oscillations in North Atlantic climate during that time(2,3), but the proposals disagree on whether increased frequency of El Nino events, the warm phase of ENSO, was linked to North Atlantic warm or cold periods. Here we present a high-resolution record of surface moisture, based on the degree of peat humification and the ratio of sedges to grass, from northern Queensland, Australia, covering the past 45,000 yr. We observe millennial-scale dry periods, indicating periods of frequent El Nino events ( summer precipitation declines in El Nino years in northeastern Australia). We find that these dry periods are correlated to the Dansgaard - Oeschger events - millennial-scale warm events in the North Atlantic climate record - although no direct atmospheric connection from the North Atlantic to our site can be invoked. Additionally, we find climatic cycles at a semiprecessional timescale (, 11,900 yr). We suggest that climate variations in the tropical Pacific Ocean on millennial as well as orbital timescales, which determined precipitation in northeastern Australia, also exerted an influence on North Atlantic climate through atmospheric and oceanic teleconnections.

Timing and structure of the 8.2 kyr BP event inferred from delta(18)O records of stalagmites from China, Oman, and Brazil

Oxygen isotope records of stalagmites from China and Oman reveal a weak summer monsoon event, with a double-plunging structure, that started 8.21 +/- 0.02 kyr B. P. An identical but antiphased pattern is also evident in two stalagmite records from eastern Brazil, indicating that the South American Summer Monsoon was intensified during the 8.2 kyr B. P. event. These records demonstrate that the event was of global extent and synchronous within dating errors of <50 years. In comparison with recent model simulations, it is plausible that the 8.2 kyr B. P. event can be tied in changes of the Atlantic Meridional Overturning Circulation triggered by a glacial lake draining event. This, in turn, affected North Atlantic climate and latitudinal position of the Intertropical Convergence Zone, resulting in the observed low-latitude monsoonal precipitation patterns.

Age models of ODP Site 184-1143

Here we present an improved astronomical timescale since 5 Ma as recorded in the ODP Site 1143 in the southern South China Sea, using a recently published Asian summer monsoon record (hematite to goethite content ratio, Hm/Gt) and a parallel benthic d18O record. Correlation of the benthic d18O record to the stack of 57 globally distributed benthic d18O records (LR04 stack) and the Hm/Gt curve to the 65°N summer insolation curve is a particularly useful approach to obtain refined timescales. Hence, it constitutes the basis for our effort. Our proposed modifications result in a more accurate and robust chronology than the existing astronomical timescale for the ODP Site 1143. This updated timescale further enables a detailed study of the orbital variability of low-latitude Asian summer monsoon throughout the Plio-Pleistocene. Comparison of the Hm/Gt record with the d18O record from the same core reveals that the oscillations of low-latitude Asian summer monsoon over orbital scales differed considerably from the glacial-interglacial climate cycles. The popular view that summer monsoon intensifies during interglacial stages and weakens during glacial stages appears to be too simplistic for low-latitude Asia. In low-latitude Asia, some strong summer monsoon intervals appear to have also occurred during glacial stages in addition to their increased occurrence during interglacial stages. Vice versa, some notably weak summer monsoon intervals have also occurred during interglacial stages next to their anticipated occurrence during glacial stages. The well-known mid-Pleistocene transition (MPT) is only identified in the benthic d18O record but not in the Hm/Gt record from the same core. This suggests that the MPT may be a feature of high- and middle-latitude climates, possibly determined by high-latitude ice sheet dynamics. For low-latitude monsoonal climate, its orbital-scale variations respond more directly to insolation and are little influenced by high-latitude processes, thus the MPT is likely not recorded. In addition, the Hm/Gt record suggests that low-latitude Asian summer monsoon intensity has a long-term decreasing trend since 2.8 Ma with increased oscillation amplitude. This long-term variability is presumably linked to the Northern Hemisphere glaciation since then.

(Table 2) Age model of ODP Site 184-1146

Orbital tuning of benthic d18O is a common approach for assigning ages to ocean sediment records. Similar environmental forcing of the northern South China Sea and the southeast Asian cave regions allows for transfer of the speleothem d18O radiometric chronology to the planktonic and benthic d18O records from Ocean Drilling Program Site 1146, yielding a new chronology with 41 radiometrically calibrated datums, spanning the past 350 kyr. This approach also provides for an independent assessment of the accuracy of the orbitally tuned benthic d18O chronology for the last 350 kyr. The largest differences relative to the latest chronology occur in marine isotope stages (MIS) 5.4, 5.5, 6, 7, and 9.3. Prominent suborbital-scale structure believed to be global in nature is identified within MIS 5.4 and MIS 7.2. On the basis of the radiometrically calibrated chronology, the time constant of the ice sheet is found to be 5.4 kyr at the precession band (light d18O lags precession minima by -55.4°) and 10.4 kyr at the obliquity band (light d18O lags obliquity maxima by 57.4°). These values are significantly shorter than the single 17 kyr time constant originally estimated by Imbrie et al. (1984), based primarily on the timing of terminations I and II and the 15 kyr time constant used by Lisiecki and Raymo (2005, doi:10.1029/2004PA001071).

Rockmagnetic investigations of IODP Site U1302-U1303

Integrated Ocean Drilling Program (IODP) Sites U1302-U1303, drilled on the SE flank of Orphan Knoll (Labrador Sea), preserve a record of detrital layers and other proxies of hydrographic change that extend the record of ice-sheet/ocean interactions through most of the Brunhes Chron. The age model is built by tandem matching of relative paleointensity (RPI) and oxygen isotope data (d18O) from Neogloboquadrina pachyderma (sin.) to reference records, indicating a mean Brunhes sedimentation rate of 14 cm/kyr. Sedimentation back to marine isotope stage (MIS) 18 is characterized by detrital layers that are detected by higher than background gamma-ray attenuation (GRA) density, peaks in X-ray fluorescence (XRF) indicators for detrital carbonate (Ca/Sr) and detrital silicate (Si/Sr), and an ice-rafted debris (IRD) proxy (wt.% >106 µm). The age model enables correlation of Site U1302/03 to IODP Site U1308 in the heart of the central Atlantic IRD belt where an age model and a similar set of detrital-layer proxies have already been derived. Ages of Heinrich (H) layers H1, H2, H4, H5 and H6 are within ~2 kyr at the two sites (H0, H3 and H5a are not observed at Site U1308), and agree with previous work at Orphan Knoll within ~3 kyr. At Site U1308, Brunhes detrital layers are restricted to peak glacials and glacial terminations back to marine isotope stage (MIS) 16 and have near-synchronous analogs at Site U1302/03. Detrital layers at Site U1302/03 are distributed throughout the record in both glacial and most interglacial stages. We distinguish Heinrich-like layers associated with IRD from detrital layers marked by multiple detrital-layer proxies (including Ca/Sr) but usually not associated with IRD, that may be attributed to lofted sediment derived from drainage and debris-flow events funneled down the nearby Northwest Atlantic Mid-Ocean Channel (NAMOC). The prominent detrital layers at Sites U1302/03 and U1308 can be correlated to millennial scale features in the Chinese speleothem (monsoon) record over the last 400 kyr, implying a link between monsoon precipitation and Laurentide Ice Sheet (LIS) instability. The detrital-layer stratigraphy at Site U1302/03 provides a long record of LIS dynamics against which other terrestrial and marine records can be compared.

FORAMINIFERAL RECORD OF CHANGES IN SUMMER MONSOON PRECIPITATION AT THE SOUTHEASTERN BRAZILIAN CONTINENTAL MARGIN SINCE THE LAST GLACIAL MAXIMUM

Changes in the oxygen isotopic composition of the planktonic foraminifer Globigerinoides ruber and in the foraminifera faunal composition in a core retrieved from the southeastern Brazilian continental margin were used to infer past changes in the hydrological balance and monsoon precipitation in the western South Atlantic since the Last Glacial Maximum (LGM). The results suggest a first-order orbital (precessional) control on the South American Monsoon precipitation. This agrees with previous studies based on continental proxies except for LGM estimates provided by pollen records. The causes for this disagreement are discussed.

A 300 kyr record of upwelling off Oman during the Late Quaternary: evidence of the Asian southwest monsoon

In the northwest Arabian Sea upwelling occurs each summer, driven by the strong SW monsoon winds. Upwelling results in high biological productivity and a distinctive assemblage of plankton species in the surface waters off Oman that are preserved in the sediments along the Oman continental margin, creating a geologic record of monsoon-driven upwelling. Sediments recovered from the Oman continental margin during Ocean Drilling Program leg 117 provide an opportunity to examine how upwelling has varied during the late Quaternary, spanning a longer interval than piston cores recovered prior to the ODP cruise. Variations in foraminifer shell accumulation and in the relative abundance of Globigerina bulloides indicate dominant cycles of variation at 1/100 kyr, the dominant frequency of glacial-interglacial variations, and at 1/23 kyr, the frequency of precessionally driven cycles in seasonal insolation. The strongest monsoon winds (indicated by increased upwelling) occurred during interglacial times when perihelion was aligned with the summer solstice, an orbital change that increased the insolation received during summer in the northern hemisphere. During glacial times upwelling was reduced, and although the precessional cycles were still present their amplitude was smaller. At both frequencies the upwelling cycles are in phase with minimum ice volume, evidence that glacial-interglacial climate changes also include changes to the climate system that influence the low-latitude monsoon. We attribute the decrease in the monsoon winds observed during glacial times to changes in bare land albedo over Asia and/or to changes in the areal extent and seasonal cycle in Asian snow cover that decrease the summer land-sea temperature contrast. Other mechanisms may also be involved. These new upwelling time series differ substantially from previous results, however the previous work relied on cores located farther offshore where upwelling is less intense and other physical mechanisms become important. Our results support the observations derived from atmospheric general circulation models of the atmosphere that indicate that both glacial boundary conditions, and the strength of summer insolation are important variables contributing to cycles in the monsoon winds during the late Quaternary.

The effect of doubled CO2 and model basic state biases on the monsoon-ENSO system. II: Changing ENSO regimes

Integrations of a fully-coupled climate model with and without flux adjustments in the equatorial oceans are performed under 2×CO2 conditions to explore in more detail the impact of increased greenhouse gas forcing on the monsoon-ENSO system. When flux adjustments are used to correct some systematic model biases, ENSO behaviour in the modelled future climate features distinct irregular and periodic (biennial) regimes. Comparison with the observed record yields some consistency with ENSO modes primarily based on air-sea interaction and those dependent on basinwide ocean wave dynamics. Simple theory is also used to draw analogies between the regimes and irregular (stochastically forced) and self-excited oscillations respectively. Periodic behaviour is also found in the Asian-Australian monsoon system, part of an overall biennial tendency of the model under these conditions related to strong monsoon forcing and increased coupling between the Indian and Pacific Oceans. The tropospheric biennial oscillation (TBO) thus serves as a useful descriptor for the coupled monsoon-ENSO system in this case. The presence of obvious regime changes in the monsoon-ENSO system on interdecadal timescales, when using flux adjustments, suggests there may be greater uncertainty in projections of future climate, although further modelling studies are required to confirm the realism and cause of such changes.

The effect of doubled CO2 and model basic state biases on the monsoon-ENSO system. I: Mean response and interannual variability

The impact of doubled CO2 concentration on the Asian summer monsoon is studied using a coupled ocean-atmosphere model. Both the mean seasonal precipitation and interannual monsoon variability are found to increase in the future climate scenario presented. Systematic biases in current climate simulations of the coupled system prevent accurate representation of the monsoon-ENSO teleconnection, of prime importance for seasonal prediction and for determining monsoon interannual variability. By applying seasonally varying heat flux adjustments to the tropical Pacific and Indian Ocean surface in the future climate simulation, some assessment can be made of the impact of systematic model biases on future climate predictions. In simulations where the flux adjustments are implemented, the response to climate change is magnified, with the suggestion that systematic biases may be masking the true impact of increased greenhouse gas forcing. The teleconnection between ENSO and the Asian summer monsoon remains robust in the future climate, although the Indo-Pacific takes on more of a biennial character for long periods of the flux-adjusted simulation. Assessing the teleconnection across interdecadal timescales shows wide variations in its amplitude, despite the absence of external forcing. This suggests that recent changes in the observed record cannot be distinguished from internal variations and as such are not necessarily related to climate change.

A record of Holocene climate change from lake geochemical analyses in southeastern Arabia

Lacustrine sediments from southeastern Arabia reveal variations in lake level corresponding to changes in the strength and duration of Indian Ocean Monsoon (IOM) summer rainfall and winter cyclonic rainfall. The late glacial/Holocene transition of the region was characterised by the development of mega-linear dunes. These dunes became stabilised and vegetated during the early Holocene and interdunal lakes formed in response to the incursion of the IOM at approximately 8500 cal yr BP with the development of C3 dominated savanna grasslands. The IOM weakened ca. 6000 cal yr BP with the onset of regional aridity, aeolian sedimentation and dune reactivation and accretion. Despite this reduction in precipitation, the take was maintained by winter dominated rainfall. There was a shift to drier adapted C4 grasslands across the dune field. Lake sediment geochemical analyses record precipitation minima at 8200, 5000 and 4200 cal yr BP that coincide with Bond events in the North Atlantic. A number of these events correspond with changes in cultural periods, suggesting that climate was a key mechanism affecting human occupation and exploitation of this region. (c) 2006 University of Washington. All rights reserved.

Climate on the Southern Black Sea coast during the Holocene: implications from the Sofular Cave record

We present the updated Holocene section of the Sofular Cave record from the southernBlackSeacoast (northern Turkey); an area with considerably different present-day climate compared to that of the neighboring Eastern Mediterranean region. Stalagmite δ13C, growth rates and initial (234U/238U) ratios provide information about hydrological changes above the cave; and prove to be more useful than δ18O for deciphering Holocene climatic variations. Between ∼9.6 and 5.4 ka BP (despite a pause from ∼8.4 to 7.8 ka BP), the Sofular record indicates a remarkable increase in rainfall amount and intensity, in line with other paleoclimate studies in the Eastern Mediterranean. During that period, enhanced summertime insolation either produced much stronger storms in the following fall and winter through high sea surface temperatures, or it invoked a regional summer monsoon circulation and rainfall. We suggest that one or both of these climatic mechanisms led to a coupling of the BlackSea and the Mediterranean rainfall regimes at that time, which can explain the observed proxy signals. However, there are discrepancies among the Eastern Mediterranean records in terms of the timing of this wet period; implying that changes were probably not always occurring through the same mechanism. Nevertheless, the Sofular Cave record does provide hints and bring about new questions about the connection between regional and large scale climates, highlighting the need for a more extensive network of high quality paleoclimate records to better understand Holoceneclimate.

The speliothem record of climate variability in Southern Arabia.

Uranium-series dated stalagmites from Oman indicate that pluvial conditions prevailed from 6.3 to 10.5, 78 to 82, 120 to 130, 180 to 200 and 300 to 330 kyr B.P.; all of these periods coincide with peak interglacials. Oxygen (δ18O) and hydrogen (δD) isotope ratios of speleothem calcite and fluid inclusions reveal the source of moisture and provide information on the amount of precipitation, respectively. δ18O and δD values of stalagmites deposited during peak interglacials vary between −8 and −4 ‰ (VPDB) and −53 and −20‰ (Vienna Standard Mean Ocean Water [VSMOW]) respectively, whereas modern stalagmites range from −2.6 to −1.1‰ in δ18O (VPDB) and −7.6 and −3.3‰ in δD (VSMOW), respectively. The growth and isotopic records indicate that during peak interglacial periods, the limit of the monsoon rainfall was shifted far north of its present location and each pluvial period was coinciding with an interglacial stage of the marine oxygen isotope record.

A 45 kyr palaeoclimate record from the lowland interior of tropical South America

We present a well-dated, high-resolution, ~ 45 kyr lake sediment record reflecting regional temperature and precipitation change in the continental interior of the Southern Hemisphere (SH) tropics of South America. The study site is Laguna La Gaiba (LLG), a large lake (95 km2) hydrologically-linked to the Pantanal, an immense, seasonally-flooded basin and the world's largest tropical wetland (135,000 km2). Lake-level changes at LLG are therefore reflective of regional precipitation. We infer past fluctuations in precipitation at this site through changes in: i) pollen-inferred extent of flood-tolerant forest; ii) relative abundance of terra firme humid tropical forest versus seasonally-dry tropical forest pollen types; and iii) proportions of deep- versus shallow-water diatoms. A probabilistic model, based on plant family and genus climatic optima, was used to generate quantitative estimates of past temperature from the fossil pollen data. Our temperature reconstruction demonstrates rising temperature (by 4 °C) at 19.5 kyr BP, synchronous with the onset of deglacial warming in the central Andes, strengthening the evidence that climatic warming in the SH tropics preceded deglacial warming in the Northern Hemisphere (NH) by at least 5 kyr. We provide unequivocal evidence that the climate at LLG was markedly drier during the last glacial period (45.0–12.2 kyr BP) than during the Holocene, contrasting with SH tropical Andean and Atlantic records that demonstrate a strengthening of the South American summer monsoon during the global Last Glacial Maximum (~ 21 kyr BP), in tune with the ~ 20 kyr precession orbital cycle. Holocene climate conditions occurred as early as 12.8–12.2 kyr BP, when increased precipitation in the Pantanal catchment caused heightened flooding and rising lake levels in LLG. In contrast to this strong geographic variation in LGM precipitation across the continent, expansion of tropical dry forest between 10 and 3 kyr BP at LLG strengthens the body of evidence for widespread early–mid Holocene drought across tropical South America.