Pieter Meiert Grootes

East Asian monsoon climate during the Late Pleistocene: high-resolution sediment records from the south China Sea

Abstract Based on the study of 10 sediment cores and 40 core-top samples from the South China Sea (SCS) we obtained proxy records of past changes in East Asian monsoon climate on millennial to bidecadal time scales over the last 220,000 years. Climate proxies such as global sea level, estimates of paleotemperature, salinity, and nutrients in surface water, ventilation of deep water, paleowind strength, freshwater lids, fluvial and/or eolian sediment supply, and sediment winnowing on the sea floor were derived from planktonic and benthic stable-isotope records, the distribution of siliciclastic grain sizes, planktonic foraminifera species, and the UK37 biomarker index. Four cores were AMS-14C-dated. Two different regimes of monsoon circulation dominated the SCS over the last two glacial cycles, being linked to the minima and maxima of Northern Hemisphere solar insolation. (1) Glacial stages led to a stable estuarine circulation and a strong O2-minimum layer via a closure of the Borneo sea strait. Strong northeast monsoon and cool surface water occurred during winter, in part fed by an inflow from the north tip of Luzon. In contrast, summer temperatures were as high as during interglacials, hence the seasonality was strong. Low wetness in subtropical South China was opposed to large river input from the emerged Sunda shelf, serving as glacial refuge for tropical forest. (2) Interglacials were marked by a strong inflow of warm water via the Borneo sea strait, intense upwelling southeast of Vietnam and continental wetness in China during summer, weaker northeast monsoon and high sea-surface temperatures during winter, i.e. low seasonality. On top of the long-term variations we found millennial- to centennial-scale cold and dry, warm and humid spells during the Holocene, glacial Terminations I and II, and Stage 3. The spells were coeval with published variations in the Indian monsoon and probably, with the cold Heinrich and warm Dansgaard–Oeschger events recorded in Greenland ice cores, thus suggesting global climatic teleconnections. Holocene oscillations in the runoff from South China centered around periodicities of 775 years, ascribed to subharmonics of the 1500-year cycle in oceanic thermohaline circulation. 102/84-year cycles are tentatively assigned to the Gleissberg period of solar activity. Phase relationships among various monsoon proxies near the onset of Termination IA suggest that summer-monsoon rains and fluvial runoff from South China had already intensified right after the last glacial maximum (LGM) insolation minimum, coeval with the start of Antarctic ice melt, prior to the δ18O signals of global sea-level rise. Vice versa, the strength of winter-monsoon winds decreased in short centennial steps only 3000–4000 years later, along with the melt of glacial ice sheets in the Northern Hemisphere.

Validity of the temperature reconstruction from water isotopes in ice cores

Well-documented present-day distributions of stable water isotopes (HDO and H218O) show the existence, in middle and high latitudes, of a linear relationship between the mean annual isotope content of precipitation (δD and δ18O) and the mean annual temperature at the precipitation site. Paleoclimatologists have used this relationship, which is particularly well obeyed over Greenland and Antarctica, to infer paleotemperatures from ice core data. There is, however, growing evidence that spatial and temporal isotope/surface temperature slopes differ, thus complicating the use of stable water isotopes as paleothermometers. In this paper we review empirical estimates of temporal slopes in polar regions and relevant information that can be inferred from isotope models: simple, Rayleigh-type distillation models and (particularly over Greenland) general circulation models (GCMs) fitted with isotope tracer diagnostics. Empirical estimates of temporal slopes appear consistently lower than present-day spatial slopes and are dependent on the timescale considered. This difference is most probably due to changes in the evaporative origins of moisture, changes in the seasonality of the precipitation, changes in the strength of the inversion layer, or some combination of these changes. Isotope models have not yet been used to evaluate the relative influences of these different factors. The apparent disagreement in the temporal and spatial slopes clearly makes calibrating the isotope paleothermometer difficult. Nevertheless, the use of a (calibrated) isotope paleothermometer appears justified; empirical estimates and most (though not all) GCM results support the practice of interpreting ice core isotope records in terms of local temperature changes.

The Greenland Ice Sheet Project 2 Depth-age Scale: Methods and Results

The Greenland Ice Sheet Project 2 (GISP2) depth-age scale is presented based on a multiparameter continuous count approach, to a depth of 2800 m, using a systematic combination of parameters that have never been used to this extent before. The ice at 2800 m is dated at 110,000 years B.P. with an estimated error ranging from 1 to 10% in the top 2500 m of the core and averaging 20% between 2500 and 2800 m. Parameters used to date the core include visual stratigraphy, oxygen isotopic ratios of the ice, electrical conductivity measurements, laser-light scattering from dust, volcanic signals, and major ion chemistry. GISP2 ages for major climatic events agree with independent ages based on varve chronologies, calibrated radiocarbon dates, and other techniques within the combined uncertainties. Good agreement also is obtained with Greenland Ice Core Project ice core dates and with the SPECMAP marine timescale after correlation through the δ 18 O of O 2 . Although the core is deformed below 2800 m and the continuity of the record is unclear, we attempted to date this section of the core on the basis of the laser-light scattering of dust in the ice.

Changes in Atmospheric Circulation and Ocean Ice Cover over the North Atlantic During the Last 41,000 Years

High-resolution, continuous multivariate chemical records from a central Greenland ice core provide a sensitive measure of climate change and chemical composition of the atmosphere over the last 41,000 years. These chemical series reveal a record of change in the relative size and intensity of the circulation system that transported air masses to Greenland [defined here as the polar circulation index (PCI)] and in the extent of ocean ice cover. Massive iceberg discharge events previously defined from the marine record are correlated with notable expansions of ocean ice cover and increases in PCI. During stadials without discharge events, ocean ice cover appears to reach some common maximum level. The massive aerosol loadings and dramatic variations in ocean ice cover documented in ice cores should be included in climate modeling.

The Little Ice Age as Recorded in the Stratigraphy of the Tropical Quelccaya Ice Cap

The analyses of two ice cores from a southern tropical ice cap provide a record of climatic conditions over 1000 years for a region where other proxy records are nearly absent. Annual variations in visible dust layers, oxygen isotopes, microparticle concentrations, conductivity, and identification of the historical (A.D. 1600) Huaynaputina ash permit accurate dating and time-scale verification. The fact that the Little Ice Age (about A.D. 1500 to 1900) stands out as a significant climatic event in the oxygen isotope and electrical conductivity records confirms the worldwide character of this event.

Correlation of marine 14C ages from the nordic seas with the GISP2 isotope record : Implications for 14C calibration beyond 25 ka BP

We present two new high-resolution sediment records from the southwestern Iceland and Norwegian Seas that were dated by numerous (super 14) C ages up to 54 (super 14) C ka BP. Based on various lines of evidence, the local (super 14) C reservoir effect was restricted to 400-1600 yr. The planktic stable isotope records reveal several meltwater spikes that were sampled with an average time resolution of 50 yr in PS2644 and 130 yr in core 23071 during isotope stage 3. Most of the delta (super 18) O spikes correlate peak-by-peak to the stadials and cold rebounds of the Dansgaard-Oeschger cycles in the annual-layer counted GISP2 ice core, with the major spikes reflecting the Heinrich events 1-6. This correlation indicates large fluctuations in the calibration of (super 14) C ages between 20 and 54 (super 14) C ka BP. Generally the results confirm the (super 14) C age shifts as predicted by the geomagnetic model of Laj, Mazaud and Duplessy (1996). However, the amplitude and speed of the abrupt decrease and subsequent major increase of our (super 14) C shifts after 45 (super 14) C ka BP clearly exceed the geomagnetic prediction near 40-43 and 32-34 calendar (cal) ka BP. At these times, the geomagnetic field intensity minima linked to the Laschamp and the Mono Lake excursions and confirmed by a local geomagnetic record, probably led to a sudden increase in cosmogenic (super 14) C and (super 10) Be production, giving rise to excess (super 14) C in the atmosphere of up to 1200 per mil.

Thermohaline instability in the North Atlantic during meltwater events: Stable isotope and ice‐rafted detritus records from Core SO75‐26KL, Portuguese Margin

A benthic isotope record has been measured for core SO75-26KL from the upper Portuguese margin (1099 m water depth) to monitor the response of thermohaline overturn in the North Atlantic during Heinrich events. Evaluating benthic δ18O in TS diagrams in conjunction with equilibrium δc fractionation implies that advection of Mediterranean outflow water (MOW) to the upper Portuguese margin was significantly reduced during the last glacial (< 15% compared to 30% today). The benthic isotope record along core SO75-26KL therefore primarily monitors variability of glacial North Atlantic conveyor circulation. The 14C-accelerator mass spectrometry ages of 13.54±.07 and 20.46±.12 ka for two ice-rafted detritus (IRD) layers in the upper core section and an interpolated age of 36.1 ka for a third IRD layer deeper in the core are in the range of published 14C ages for Heinrich events H1, H2, and H4. Marked depletion of benthic δ13C by 0.7–1.1‰ during the Heinrich events suggests reduced thermohaline overturn in the North Atlantic during these events. Close similarity between meltwater patterns (inferred from planktonic δ18O) at Site 609 and ventilation patterns (inferred from benthic δ13C) in core SO75-26KL implies coupling between thermohaline overturn and surface forcing, as is also suggested by ocean circulation models. Benthic δ13C starts to decrease 1.5–2.5 kyr before Heinrich events Hl and H4, fully increased values are reached 1.5–3 kyr after the events, indicating a successive slowdown of thermohaline circulation well before the events and resumption of the conveyors full strength well after the events. Benthic δ13C changes in the course of the Heinrich events show subtle maxima and minima suggesting oscillatory behavior of thermohaline circulation, a distinct feature of thermohaline instability in numerical models. Inferrred gradual spin-up of thermohaline circulation after Hl and H4 is in contrast to abrupt wanning in the North Atlantic region that is indicated by sudden increases in Greenland ice core δ18O and in marine faunal records from the northern North Atlantic. From this we infer that thermohaline circulation can explain only in part the rapid climatic oscillations seen in glacial sections of the Greenland ice core record.

The Leibniz-Labor AMS facility at the Christian-Albrechts University, Kiel, Germany

The AMS facility of the Leibniz-Labor fur Altersbestimmung und Isotopenforschung of the Christian-Albrechts Universitat is based on a 3 MV Tandetron from High Voltage Engineering Europa (HVEE) with a single cesium sputter ion source and a separator/recombinator for simultaneous injecton of the three isotopic carbon beams. The AMS system is similar to those at the Woods Hole Oceanographic Institution, USA, and the University of Groningen, The Netherland, but it has some new features based on experience at these two facilities. These include improved vacuum seals, beam diagnostics, X-ray and background suppression as well as a more reliable system control through a PLC-unit with a serial line to the main system computer. The open system design of the beam optics allows significant horizontal and vertical movement of the ion beams without loss to the walls of the system. This leads to plateaus in the response of the isotope beams and ratios to changing values of various ion optical elements. Combined with highly stable power supplies, this gives reproducible measurements. The acceptance tests, e.g., showed that Poisson counting statistics at 0.15% and 0.22% respectively, determined the statistical uncertainty in the 14C12C ratios measured for the individual samples of two test series. Strong discrimination of unwanted ions results in low background count rates in the detector, equivalent to an apparent age of 75000 years at present, in spite of the open architecture. Routine measurements since late January 1996 (to late May 1996) have dated 127 unknown samples, mostly foraminifera. The prototype of the carbonate to CO2 conversion system and the graphite system used for the measurements are also described.

Fundamental Modes and Abrupt Changes in North Atlantic Circulation and Climate over the last 60 ky — Concepts, Reconstruction and Numerical Modeling

Centennial- to millennial-scale changes in global climate over the last 60 ky were first documented in ice cores from Greenland, with ice sheets around the North Atlantic and its thermohaline circulation (THC) as prime candidates for a potential trigger mechanism. To reach a new quality in understanding the origin and causal links behind these changes, two strategies were intimately tied together in this synthesis, high-resolution 3-D ocean modeling and paleoceanographic reconstructions. Here, five time series with a time resolution of several decades and various time slices of surface and deep-water paleoceanography were established from hundreds of deep-sea cores for the purpose of monitoring rapid changes across the North Atlantic and testing or initiating model results. Three fundamental modes were found to operate Atlantic THC. Today, mode I shows intensive formation of North Atlantic Deep Water (NADW) and strong heat and moisture fluxes to the continents adjacent to the North Atlantic. Peak glacial mode II leads to a reduction in NADW formation by 30-50%, in line with a clear drop in heat flux to Europe. The glacial Nordic Seas, however, remain ice-free during summer and little influenced by meltwater, in contrast to the sea west ofIreland, where iceberg meltwater blocks an eastbound flow into the Norwegian Sea and induces a cold longshore current from Faeroe to the Pyrenees. The subsequent Heinrich 1 (HI) meltwater mode III leads to an entire stop in NADW and intermediate-water production as well as a reversed pattern of THC, stopping any heat advection from the central and South Atlantic to the north. In contrast to earlier views, the Younger Dryas, possibly induced by Siberian meltwater, began with mode I and ended with mode III, continuing into the Preboreal. Modeling the impact of modes I to III on the global carbon budget, we find that the atmosphere has lost 34-54 ppmv CO2 from interglacial to glacial times, but has gained 23-62 ppmv CO2 at the end of HI within a few decades, equivalent to 33-90% of modem, man-made CO2 release. The robust 1500-y Dansgaard- Oeschger (D-O) cycles and their multiples of as much as 7200 years, the Heinrich event cycles, are tied to periodical changes between THC modes I/II and II/III. In the Irminger Sea rapid D-O coolings are in phase with initial meltwater injections from glaciers on East Greenland, here suggesting an internal trigger process in accordance with binge-purge models. Ice rafting from East Greenland and Iceland occurs only 240-280 y later, probably inducing a slight sea-level rise and, in tum, Heinrich ice rafting from the Laurentian ice sheet during H1, H2, H4, H5. At H1 a major surge from the Barents shelf has lagged initial cooling by 1500 y and entails the most prominent and extended reversal in Atlantic THC over the last 60 ky (probably also at the end of glacial stage 4, at H6). Meltwater stratification in the Inninger Sea reaches its maximum only 640 y after initial meltwater injection and induces, via seasonal sea-ice formation, brine-water injections down to 4 km water depth, signals leading the classic D-O jump to maximum warmth by only 125 y. It may be inferred from this short-phase lag that brine water-controlled deep-water formation probably entrains warm water from further south, thereby forming the key trigger mechanism for the final tum-on of the Atlantic THC mode II roughly within a decade (or mode I, in case of favorable Milankovitch forcing).

Wisconsinan and Holocene Climate History from an Ice Core at Taylor Dome, Western Ross Embayment, Antarctica

Geochemical data and geophysical measurements from a 554‐m ice‐core from Taylor Dome, East Antarctica, provide the basis for climate reconstruction in the western Ross Embayment through the entire Wisconsinan and Holocene. In comparison with ice cores from central East and West Antarctica, Taylor Dome shows greater variance of temperature, snow accumulation, and aerosol concentrations, reflecting significant variability in atmospheric circulation and air mass moisture content. Extreme aridity during the last glacial maximum at Taylor Dome reflects both colder temperatures and a shift in atmospheric circulation patterns associated with the advance of the Ross Sea ice sheet and accounts for regional alpine glacier retreats and high lake levels in the Dry Valleys. Inferred relationships between spatial accumulation gradients and ice sheet configuration indicate that advance of the Ross Sea ice sheet began in late marine isotope stage 5 or early stage 4. Precise dating of the Taylor Dome core achieved by trace‐gas correlation with central Greenland ice cores shows that abrupt deglacial warming at Taylor Dome was near‐synchronous with the ∼14.6 ka warming in central Greenland and lags the general warming trend in other Antarctic ice cores by at least 3000 years. Deglacial warming was following by a warm interval and transient cooling between 14.6 and 11.7 ka, synchronous with the Bølling/Allerød warming and Younger Dryas cooling events in central Greenland, and out of phase with the Antarctic Cold Reversal recorded in the Byrd (West Antarctica) ice core. Rapid climate changes during marine isotope stages 4 and 3 at Taylor Dome are similar in character to, and may be in phase with, the Northern Hemisphere stadial–interstadial (Dansgaard–Oeschger) events. Results from Taylor Dome illustrate the importance of obtaining ice cores from multiple Antarctic sites, to provide wide spatial coverage of past climate and ice dynamics.