Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Holocene and the Last Glacial Maximum

We measured oxygen isotopes and Mg/Ca ratios in the surface-dwelling planktonic foraminifer Globigerinoides ruber (white s.s.) and the thermocline dweller Pulleniatina obliquiloculata to investigate upper ocean spatial variability in the Indo-Pacific Warm Pool (IPWP). We focused on three critical time intervals: the Last Glacial Maximum (LGM; 18-21.5 ka), the early Holocene (8-9 ka), and the late Holocene (0-2 ka). Our records from 24 stations in the South China Sea, Timor Sea, Indonesian seas, and western Pacific indicate overall dry and cool conditions in the IPWP during the LGM with a low thermal gradient between surface and thermocline waters. During the early Holocene, sea surface temperatures increased by ~3°C over the entire region, indicating intensification of the IPWP. However, in the eastern Indian Ocean (Timor Sea), the thermocline gradually shoaled from the LGM to early Holocene, reflecting intensification of the subsurface Indonesian Throughflow (ITF). Increased surface salinity in the South China Sea during the Holocene appears related to northward displacement of the monsoonal rain belt over the Asian continent together with enhanced influx of saltier Pacific surface water through the Luzon Strait and freshwater export through the Java Sea. Opening of the freshwater portal through the Java Sea in the early Holocene led to a change in the vertical structure of the ITF from surface- to thermocline-dominated flow and to substantial freshening of Timor Sea thermocline waters.

Interstitial-water chemistry and stable isotope record of Prydz Bay sediments

Leg 119 of the Ocean Drilling Program (ODP) provided the first opportunity to study the interstitial-water chemistry of the eastern Antarctic continental margin. Five sites were cored in a northwest-southeast transect of Prydz Bay that extended from the top of the continental slope to within 30 km of the coastline. Geological studies of the cores reveal a continental margin that has evolved through terrestrial, glacial, and glacial-marine environments. Chemical and stable isotopic analyses of the interstitial-waters were performed to determine the types of depositional environments and the diagenetic and hydrologic processes that are operating in this unusual marine environment. Highly compacted glacial sediments provide an effective barrier to the vertical diffusion of interstitial-water solutes. Meteoric water from the Antarctic continent appears to be flowing into Prydz Bay sediments through the sequence of terrestrial sediments that lie underneath the glacial sediments. The large amounts of erosion associated with glacial advances appear to have had the effect of limiting the amount of marine organic matter that is incorporated into the sediments on the continental shelf. Although all of the sites cored in Prydz Bay exhibit depletions in dissolved sulfate with increasing depth, the greatest bacterial activity is associated with a thin layer of diatom ooze that coats the seafloor of the inner bay. Results of alkalinity modeling, thermodynamic calculations, and strontium analyses indicate that (1) ocean bottom waters seaward of Site 740 are undersaturated with respect to both calcite and aragonite, (2) interstitial waters at each site become saturated or supersaturated with respect to calcite and aragonite with increasing depth, (3) precipitation of calcium carbonate reduces the alkalinity of the pore waters with increasing depth, and (4) recrystallization of aragonite to calcite accounts for 24% of the pore-water strontium. Weathering of unstable terrestrial debris and cation exchange between clay minerals and pore fluids are the most probable chemical processes affecting interstitial water cation gradients.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.

Neogene and Quaternary paleoenvironmental history at DSDP Leg 71 Holes

The Neogene and Quaternary sedimentary record of Leg 71 and previously drilled sequences from the Southern Ocean reveal evidence of a major late Miocene change of oceanic and glacial conditions in the southern high latitudes during paleomagnetic Chron 9. The characteristics of late Miocene sedimentation and in particular the study of erosional patterns and ice-rafted debris suggest the following conclusions. 1) In the late Miocene, the Polar Front first migrated to the northern latitudes of the Southern Ocean and surface water temperatures became similar to those of today. 2) Extensive ice shelves or ice tongues were not present along the Antarctic margin until late Chron 9 (about 9.0 Ma). 3) Before Chron 9, West Antarctica was occupied by an archipelago and the West Antarctic Sea. 4) Extensive ice shelves formed in the West Antarctic region, eventually coalescing and thickening to form the grounded West Antarctic ice sheet by Chron 9. 5) The newly formed West Antarctic ice sheet was probably unstable and frequently became an ungrounded ice shelf, thus accounting for the scarcity of late Miocene ice-rafted debris. 6) Extensive erosion or nondeposition of sediment was probably the result of increased Antarctic Bottom Water (AABW) formation in the West Antarctic region during the initial formation of extensive West Antarctic ice shelves and during periods when the West Antarctic ice sheet was ungrounded. 7) In the Southwest Atlantic, AABW velocity waned during the latest Miocene. During the late Gilbert Chron a major and permanent change occurred in the pattern of ice-rafting to the South Atlantic, and after 4.35 Ma the increased IRD accumulation rate and frequency of major episodes of IRD accumulation suggest increased stability of the West Antarctic ice sheet. In addition, radiolarian faunas of Hole 514 record at least eight migrations of the Polar Front to the north of the site during the past 4.07 m.y. An apparent increase in the frequency of Polar Front migrations occurred about 2.7-2.6 Ma, possibly in response to oceanic change induced by fluctuations in glacial conditions of the Northern Hemisphere.

Nd bulk isotopes and 40Ar/39Ar ages of individual minerals from Aurora and Wilkes sub-glacial basin sediments

The Wilkes and Aurora basins are large, low-lying sub-glacial basins that may cause areas of weakness in the overlying East Antarctic ice sheet. Previous work based on ice-rafted debris (IRD) provenance analyses found evidence for massive iceberg discharges from these areas during the late Miocene and Pliocene. Here we characterize the sediments shed from the inferred areas of weakness along this margin (94°E to 165°E) by measuring40Ar/39Ar ages of 292 individual detrital hornblende grains from eight marine sediment core locations off East Antarctica and Nd isotopic compositions of the bulk fine fraction from the same sediments. We further expand the toolbox for Antarctic IRD provenance analyses by exploring the application of 40Ar/39Ar ages of detrital biotites; biotite as an IRD tracer eliminates lithological biases imposed by only analyzing hornblendes and allows for characterization of samples with low IRD concentrations. Our data quadruples the number of detrital 40Ar/39Ar ages from this margin of East Antarctica and leads to the following conclusions: (1) Four main sectors between the Ross Sea and Prydz Bay, separated by ice drainage divides, are distinguishable based upon the combination of 40Ar/39Ar ages of detrital hornblende and biotite grains and the e-Nd of the bulk fine fraction; (2) 40Ar/39Ar biotite ages can be used as a robust provenance tracer for this part of East Antarctica; and (3) sediments shed from the coastal areas of the Aurora and Wilkes sub-glacial basins can be clearly distinguished from one another based upon their isotopic fingerprints.

Changes in bottom water corrosiveness and carbonate ion concentrations in the sub-Antarctic Atlantic during the last glacial period

The glacial climate system transitioned rapidly between cold (stadial) and warm (interstadial) conditions in the Northern Hemisphere. This variability, referred to as Dansgaard-Oeschger variability, is widely believed to arise from perturbations of the Atlantic Meridional Overturning Circulation. Evidence for such changes during the longer Heinrich stadials has been identified, but direct evidence for overturning circulation changes during Dansgaard-Oeschger events has proven elusive. Here we reconstruct bottom water [CO3]2- variability from B/Ca ratios of benthic foraminifera and indicators of sedimentary dissolution, and use these reconstructions to infer the flow of northern-sourced deep water to the deep central sub-Antarctic Atlantic Ocean. We find that nearly every Dansgaard-Oeschger interstadial is accompanied by a rapid incursion of North Atlantic Deep Water into the deep South Atlantic. Based on these results and transient climate model simulations, we conclude that North Atlantic stadial-interstadial climate variability was associated with significant Atlantic overturning circulation changes that were rapidly transmitted across the Atlantic. However, by demonstrating the persistent role of Atlantic overturning circulation changes in past abrupt climate variability, our reconstructions of carbonate chemistry further indicate that the carbon cycle response to abrupt climate change was not a simple function of North Atlantic overturning.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Analysis of organic matter from DSDP Leg 81 Holes

Organic geochemical and visual kerogen analyses were carried out on approximately 50 samples from Leg 81 (Rockall Plateau, North Atlantic). The sediments are from four sites (Sites 552-555), Pleistocene to Paleocene in age, and represent significantly different depositional environments and sources of organic matter. The Pleistocene glacial-interglacial cycles show differences in sedimentary organic matter based on Rock-Eval pyrolysis, organic phosphorus, and pyrolysis/mass-spectrometry analyses. Glacial samples contain more organic carbon, with a larger proportion of reworked organic matter. This probably reflects increased erosion of continental and shelf areas as a result of low sea level stands. Inter glacial samples contain a larger proportion of marine organic matter as determined by organic phosphorus and pyrolysis analyses. This immature, highly oxidized marine organic matter may be associated with the skeletal organic matrix of calcareous organisms. In addition, Rock-Eval data indicate no significant inorganic-carbonate contribution to the S3 pyrolysis peak. The Pliocene-Miocene sediments consist of pelagic, biogenic carbonates. The organic matter is similar to that of the Pleistocene interglacial periods; a mixture of oxidized marine organic matter and reworked, terrestrial detritus. The Paleocene-Oligocene organic matter reflects variations in source and depositional factors associated with the isolation of Rockall from Greenland. Paleocene sediments contain primarily terrestrial organic matter with evidence of in situ thermal stress resulting from interbedded lava flows. Late Paleocene and early Eocene organic matter suggests a highly oxidized marine environment, with major periods of deposition of terrestrially derived organic matter. These fluctuations in organic-matter type are probably the result of episodic shallowing and deepening of Rockall Basins. The final stage of Eocene/Oligocene sedimentation records the accelerating subsidence of Rockall and its isolation from terrestrial sources (Rockall and Greenland). This is shown by the increasingly marine character of the organic matter. The petroleum potential of sediments containing more than 0.5% organic carbon is poor because of their thermal immaturity and their highly oxidized and terrestrial organic-matter composition.

Mean sortable silt values from Holocene and glacial Scotia Sea sediments

The Antarctic Circumpolar Current is key to the mixing and ventilation of the world's oceans. This current flows from west to east between about 45° and 70° S connecting the Atlantic, Pacific and Indian oceans, and is driven by westerly winds and buoyancy forcing. High levels of productivity in the current regulate atmospheric CO2 concentrations. Reconstructions of the current during the last glacial period suggest that flow speeds were faster or similar to present, and it is uncertain whether the strength and position of the westerly winds changed. Here we reconstruct Antarctic Circumpolar Current bottom speeds through the constricting Drake Passage and Scotia Sea during the Last Glacial Maximum and Holocene based on the mean grain size of sortable silt from a suite of sediment cores. We find essentially no change in bottom flow speeds through the region, and, given that the momentum imparted by winds, and modulated by sea-ice cover, is balanced by the interaction of these flows with the seabed, this argues against substantial changes in wind stress. However, glacial flow speeds in the sea-ice zone south of 56° S were significantly slower than present, whereas flow in the north was faster, but not significantly so. We suggest that slower flow over the rough topography south of 56° S may have reduced diapycnal mixing in this region during the last glacial period, possibly reducing the diapycnal contribution to the Southern Ocean overturning circulation.

Glacial sedimentation in ODP Leg 119 samples

Although scientific evidence prior to that from ODP Leg 119 indicates the presence of an ice sheet on East Antarctica by at least the earliest Oligocene, the question as to the size and stability of that initial ice sheet is still contested. Current hypotheses include (1) the presence of a small ice sheet in the earliest Oligocene with stepwise growth during the Neogene, (2) the presence of a continental-sized ice sheet in the late middle Eocene with no major evidence of subsequent deglaciation, and (3) the presence of glacial ice in the earliest Oligocene with a major ice sheet during the mid-Oligocene, followed by growth and decay of several ice sheets with characteristics similar to the temperate ice sheets of the Pleistocene of North America but with changes over a longer time scale (millions of years vs. 100,000 yr). Principal results from Leg 119 suggest the presence of significant late middle and late Eocene glaciation in East Antarctica and the presence of a continental-size ice sheet in East Antarctica during the earliest Oligocene. Although the Leg 119 results provide only glimpses of the Neogene glacial history of East Antarctica, they do provide evidence of fluctuations in the extent of the ice sheet and the waxing and waning of glaciers across the Prydz Bay shelf during the later part of the late Miocene and Pliocene.

Stable carbon and oxygen isotope ratios of Planulina wuellerstorfi from sediments of the Sierra Leone Rise

Causes of change in deep water delta13C can be either global or local in extent. Global causes include (1) climatically-induced changes in the amount of terrestrial biomass which alter the average carbon isotopic composition of the oceanic reservoir (Shackleton, 1977), and (2) erosion and deposition of organic-rich, continental shelf sediments during sea level fluctuations which change the mean oceanic carbon: phosphorus ratio (Broecker, 1982 doi:10.1016/0079-6611(82)90007-6). Regional gradients of delta13C are created by remineralization of organic detritus within the deep ocean itself thus reflecting the distribution of water masses and modern thermohaline flow. Changes in a single geological record of benthic foraminiferal delta13C can result from any combination of these global and abyssal circulation effects. By sampling a large number of cores collected over a wide bathymetric range yet confined to a small geographical region we have minimized the ambiguity. We can assume that each delta13C record was equally affected by global causes of delta13C variation. The differences seen between the delta13C records must, therefore, reflect changes in the distribution of delta13C in the deep ocean. We interpret these differences in distribution in terms of changes in the ocean's abyssal circulation. Benthic foraminiferal carbon isotopic evidence from a suite of Sierra Leone Rise cores indicates that the deeper parts of the eastern Atlantic basins underwent a reduction in [O2] during the maximum of the last glaciation. Reduced advection of O2-rich deep water through low-latitude fracture zones, associated with increased delivery of organic matter to the deep ocean, lowered the delta13C of deep water SumCO2 at all depths below the sill separating the eastern and western Atlantic basins (Metcalf et al., 1964 doi:10.1016/0011-7471(64)91078-2). This decreased advection into the eastern Atlantic Ocean coincides with the overall decrease in deep water production in the North Atlantic during the last glacial maximum (Curry and Lohmann, 1982 doi:10.1016/0033-5894(82)90071-0; Boyle and Keigwin, 1982 doi:10.1126/science.218.4574.784; Schnitker, 1979 doi:10.1016/0377-8398(79)90020-3; Streeter and Shackleton, 1979 doi:10.1126/science.203.4376.168).