Planktonic foraminifera of sediment core MD01-2390

Sea-surface temperature (SST) estimates in the sediment core MD01-2390 based on planktonic foraminiferal species abundances using five different transfer function techniques suggest nearly unchanged or unusually higher temperatures in the tropical southern South China Sea (SCS) during the Last Glacial Maximum (LGM) relative to modern temperatures. These results are in contrast to substantial cooling of 2-5 °C inferred by geochemical (Uk'37, Mg/Ca ratios) and terrestrial proxies from the western tropical Pacific region. Using multivariate statistics we show that the glacial southern SCS harboured unique planktonic foraminiferal assemblages that have no modern analogs. Analyses of faunal variation through the core reveal that planktonic foraminiferal assemblages responded to temperature changes inferred from Mg/Ca data but that this signal is subdued by superimposed variations in the relative abundance of Pulleniatina obliquiloculata and Neogloboquadrina pachyderma (dextral). These species occur in glacial samples at proportions that are not observed in the calibration data set. The glacial high abundance of N. pachyderma (dextral) are interpreted to reflect a seasonal (winter) inflow of cold surface water from the northeast via the Bashi Strait due to the combined effects of an intensified winter monsoon, a southward shift of the polar front and the eastward migration of the Kuroshio Current. In contrast, processes controlling the high relative abundances of P. obliquiloculata during the LGM may be unique to the southern SCS. We propose a scenario involving a stronger (winter) mixing or enhanced upwelling due to an intensified winter monsoon that prevented shallow-dwelling, warm indicators to establish larger populations during the LGM. Our results indicate that a no-analog behaviour of planktonic foraminifera faunas is responsible for the warm glacial conditions in this part of the western Pacific warm pool as implied by foraminiferal transfer functions and that a more significant surface cooling in the region as implied by terrestrial and geochemical (Mg/Ca ratios; alkenone unsaturation index) marine proxies is a more likely scenario.

Long sedimentary records of sediment cores from the Chilean margin

The Antarctic Intermediate Water (AAIW) is a key player in global-scale oceanic overturning processes and an important conduit for heat, fresh water, and carbon transport. The AAIW past variability is poorly understood mainly due to the lack of sedimentary archives at intermediate water depths. We present records of benthic stable isotopes from sediments retrieved with the seafloor drill rig MARUM-MeBo at 956 m water depth off northern Chile (GeoB15016, 27°29.48'S, 71°07.58'W) that extend back to 970 ka. The sediments at this site are presently deposited at the boundary between AAIW and Pacific Deep Water (PDW). For previous peak interglacials, our results reveal similar benthic d13C values at site GeoB15016 and of a newly generated stack of benthic d13C from various deep Pacific cores representing the "average PDW." This suggests, unlike today, the absence of AAIW at the site and the presence of nearly pure PDW. In contrast, more positive d13C values at site GeoB15016 compared to the stack imply a considerable AAIW contribution during cold phases of interglacials and especially during glacials. Besides, we used three short sediment cores to reconstruct benthic d13C values from the AAIW core during the last glacial and found a d13C signature similar to today's. Assuming that this was the case also for the past 970 kyr, we demonstrate that sea level changes and latitudinal migrations of the AAIW formation site can only account for about 50% of the full range of past d13C increases at site GeoB15016 during cold periods. Other processes that could explain the remaining of the positive d13C anomalies are increases in glacial AAIW production and/or deeper convection of the AAIW with respect to preceding interglacials.

Stable isotope record, and abundances of ice rafted debris and Neogloboquadrina pachyderma in sediments of the Iberian margin

This paper documents the migration of the Polar Front (PF) over the Iberian margin during some of the cold climatic extremes of the last 45 ka. It is based on a compilation of robust and coherent paleohydrological proxies obtained from eleven cores distributed between 36 and 42°N. Planktonic delta18O (Globigerina bulloides), ice-rafted detritus concentrations, and the relative abundance of the polar foraminifera Neogloboquadrina pachyderma sinistral were used to track the PF position. These three data sets, compared from core to core, show a consistent evolution of the sea surface paleohydrology along the Iberian margin over the last 45 ka. We focused on five time slices representative of cold periods under distinct paleoenvironmental forcings: the 8.2 ka event and the Younger Dryas (two recent cold events occurring within high values of summer insolation), Heinrich events 1 and 4 (reflecting major episodes of massive iceberg discharges into the North Atlantic), and the Last Glacial Maximum (typifying the highest ice volume accumulated in the Northern Hemisphere). For each event, we generated schematic maps mirroring past sea surface hydrological conditions. The maps revealed that the Polar Front presence along the Iberian margin was restricted to Heinrich events. The sea surface conditions during the Last Glacial Maximum were close to those at present day, except for the northern sites which briefly experienced subarctic conditions.

Late Pleistocene stable isotope record of the Northeastern Pacific Margin

Late Pleistocene intermediate water ventilation history in the northeastern Pacific has been inferred from benthic foraminiferal carbon isotopic data from seven California margin basins. Secular variations in oceanic d13C recorded at North Pacific ODP Site 849 were subtracted from each basin record leaving a residual d13C history that reflects variations in ventilation. During the previous interglacial intermediate waters above 2000m contained less oxygen than they do today or Pacific deep water at Site 849 was better ventilated. Intermediate water ventilation began to improve during oxygen isotope stage 4 and continued to improve throughout stages 3 and 2. This enhanced ventilation was not contemporaneous at all depths and appears to have progressed upwards through the water column. The diachronous nature of these changes suggest that there was not simply an "on"/"off" mechanism which resulted in higher or lower ventilation in the North Pacific during the last glacial cycle.

Averages and ranges of paleoenvironmental proxy records of ODP Hole 108-658C (Table 1)

Several high-resolution proxy environmental records have been obtained for the last 35 kyr from ODP Hole 658C, a well-studied site ca. 200 km off Cap Blanc, NW Africa. The collective assessment based on the marine proxies (UK'37 SST, contents of TOC and chlorins, Upwelling Radiolarian Index and the percentage of Florisphaera profunda), surprisingly indicates that the last glacial maximum (LGM) was characterized by warmer sea surface temperature (SST), weaker upwelling, and lower marine productivity, compared with the preceding older glacial and subsequent deglaciation periods. Of the terrigenous proxies, the mean grain size of the non-carbonate fraction and the terrigenous alkane content indicate that wind strength and aridity were high. The weaker upwelling at the 658 site during the LGM may have resulted from changes in the strength and direction of the wind systems and/or shifts in the position and geometry of the upwelling cell.

Coccolithophore taxa and PCA factor analysis of ODP Site 1233

We provide high-resolution sea surface temperature (SST) and paleoproductivity data focusing on Termination 1. We describe a new method for estimating SSTs based on multivariate statistical analyses performed on modern coccolithophore census data, and we present the first downcore reconstructions derived from coccolithophore assemblages at Ocean Drilling Project (ODP) Site 1233 located offshore Chile. We compare our coccolithophore SST record to alkenone-based SSTs as well as SST reconstructions based on dinoflagellates and radiolaria. All reconstructions generally show a remarkable concordance. As in the alkenone SST record, the Last Glacial Maximum (LGM, 19-23 kyr B.P.) is not clearly defined in our SST reconstruction. After the onset of deglaciation, three major warming steps are recorded: from 18.6 to 18 kyr B.P. (~2.6°C), from 15.7 to 15.3 kyr B.P. (~2.5°C), and from 13 to 11.4 kyr B.P. (~3.4°C). Consistent with the other records from Site 1233 and Antarctic ice core records, we observed a clear Holocene Climatic Optimum (HCO) from ~8-12 kyr B.P. Combining the SST reconstruction with coccolith absolute abundances and accumulation rates, we show that colder temperatures during the LGM are linked to higher coccolithophore productivity offshore Chile and warmer SSTs during the HCO to lower coccolithophore productivity, with indications of weak coastal upwelling. We interpret our data in terms of latitudinal displacements of the Southern Westerlies and the northern margin of the Antarctic Circumpolar Current system over the deglaciation and the Holocene.

Bulk properties and isotopic data of sediment cores from the Bounty Trough and Tasman Sea off New Zealand

Glacial/interglacial changes in Southern Ocean's air-sea gas exchange have been considered as important mechanisms contributing to the glacial/interglacial variability in atmospheric CO2. Hence, understanding past variability in Southern Ocean intermediate- to deep-water chemistry and circulation is fundamental to constrain the role of these processes on modulating glacial/interglacial changes in the global carbon cycle. Our study focused on the glacial/interglacial variability in the vertical extent of southwest Pacific Antarctic Intermediate Water (AAIW). We compared carbon and oxygen isotope records from epibenthic foraminifera of sediment cores bathed in modern AAIW and Upper Circumpolar Deep Water (UCDW; 943 - 2066 m water depth) to monitor changes in water mass circulation spanning the past 350,000 years. We propose that pronounced freshwater input by melting sea ice into the glacial AAIW significantly hampered the downward expansion of southwest Pacific AAIW, consistent with climate model results for the Last Glacial Maximum. This process led to a pronounced upward displacement of the AAIW-UCDW interface during colder climate conditions and therefore to an expansion of the glacial carbon pool.

Geochemistry and foraminiferal stable isotoe record of sediment core SCS90-36

Changes in the Southeast Asia monsoon winds and surface circulation patterns since the last glaciation are inferred using multiple paleoceanographic indicators including planktic foraminifer faunal abundances, fauna and alkenones sea-surface temperature (SST) estimates, oxygen and carbon isotopes of planktic and benthic foraminifers, and sedimentary fluxes of carbonates and organic carbon obtained from deep-sea core SCS90-36 from the South China Sea (SCS) (17°59.70'N, 111°29.64'E at water depth 2050 m). All these paleoceanographic evidences indicate marked changes in the SCS ocean system over the last glacial toward the Holocene. Planktic foraminiferal faunal SST estimates show stable warm-season SST of 28.6°C, close to the modern value, throughout the glacial-interglacial cycle. In contrast, cold-season SST increases gradually from 23.6°C in the last glacial to a mean value of 26.4°C in the Holocene with a fluctuation of about 3°C during 13-16 ka. SST estimates by UK'37 method reveal less variability and are in average 1-3°C lower than the fauna-derived winter-season SST. These patterns reveal that the seasonality of the SST is not only higher by about 3-4°C in the glacial, but also a function of the winter season SST. Sedimentation rates decrease from the last glacial-deglacial stage to the Holocene due to a reduction in supply of terrigenous components, which led to an increase of carbonate contents. Total organic carbon (TOC) contents of primarily marine sources decrease from the last glacial-deglacial to the Holocene. The last deglaciation is also characterized by high surface productivity as indicated by increased ketones abundances and high mass accumulation rates (MAR) of the TOC and carbonates. The gradient of planktic foraminifer ocygen and carbon isotopes of between surface dwellers and deep dwellers increases significantly toward Termination I and Holocene, and is indiscernibly small in the carbon isotope gradient of between 14 and 24 ka, revealing a deep-mixing condition in surface layers prior to 10 ka. The glacial-interglacial fluctuation of the carbon isotope value of a benthic foraminifer is 0.61%. which is significantly larger than a global mean value. The large carbon isotope fluctuation indicates an amplification of marginal-sea effects which is most likely resulted from an increase in surface productivity in the northern SCS during the last glacial-deglacial stage. The multiple proxies consistently indicate that the last glacial-deglacial stage winter monsoon in the Southeast Asia was probably strengthened in the northern SCS, leading to a development of deep-mixing surface layer conditions and a more efficient nutrient cycling which supports more marine organic carbon production.

A Late Pleistocene sea level stack

Late Pleistocene sea level has been reconstructed from ocean sediment core data using a wide variety of proxies and models. However, the accuracy of individual reconstructions is limited by measurement error, local variations in salinity and temperature, and assumptions particular to each technique. Here we present a sea level stack (average) which increases the signal-to-noise ratio of individual reconstructions. Specifically, we perform principal component analysis (PCA) on seven records from 0-430 ka and five records from 0-798 ka. The first principal component, which we use as the stack, describes ~80 % of the variance in the data and is similar using either five or seven records. After scaling the stack based on Holocene and Last Glacial Maximum (LGM) sea level estimates, the stack agrees to within 5 m with isostatically adjusted coral sea level estimates for Marine Isotope Stages 5e and 11 (125 and 400 ka, respectively). When we compare the sea level stack with the d18O of benthic foraminifera, we find that sea level change accounts for about ~40 % of the total orbital-band variance in benthic d18O, compared to a 65 % contribution during the LGM-to-Holocene transition. Additionally, the second and third principal components of our analyses reflect differences between proxy records associated with spatial variations in the d18O of seawater.

Late Weichselian fluvial evolution of the southern Kara Sea

Glaciations had a profound impact on the global sea-level and particularly on the Arctic environments. One of the key questions related to this topic is, how did the discharge of the Siberian Ob and Yenisei rivers interact with a proximal ice sheet? In order to answer this question high-resolution (1-12 kHz), shallow-penetration seismic profiles were collected on the passive continental margin of the Kara Sea Shelf to study the paleo-drainage pattern of the Ob and Yenisei rivers. Both rivers incised into the recent shelf, leaving filled and unfilled river channels and river canyons/valleys connecting to a complex paleo-drainage network. These channels have been subaerially formed during a regressive phase of the global sea-level during the Last Glacial Maximum. Beyond recent shelf depths of 120 m particle transport is manifested in submarine channel-levee complexes acting as conveyor for fluvial-derived fines. In the NE area, uniform draping sediments are observed. Major morphology determining factors are (1) sea-level fluctuations and (2) LGM ice sheet influence. Most individual channels show geometries typical for meandering rivers and appear to be an order of magnitude larger than recent channel profiles of gauge stations on land. The Yenisei paleo-channels have larger dimensions than the Ob examples and could be originated by additional water release during the melt of LGM Putoran ice masses. Asymmetrical submarine channel-levee complexes with channel depths of 60 m and more developed, in some places bordered by glacially dominated morphology, implying deflection by the LGM ice masses. A total of more than 12,000 km of acoustic profiles reveal no evidence for an ice-dammed lake of greater areal extent postulated by several workers. Furthermore, the existence of the channel-levee complexes is indicative of unhindered sediment flow to the north. Channels situated on the shelf above 120-m water depth exhibit no phases of ponding and or infill during sea-level lowstand. These findings denote the non-existence of an ice sheet on large areas of the Kara Sea shelf.

Age determination and geochemical composition of sediment core DA04-31P

There is limited knowledge pertaining to the history of the Greenland Ice Sheet (GIS) during the last glacial-interglacial transition as it retreated from the continental margins to an inland position. Here we use multiproxy data, including ice-rafted debris (IRD); planktonic isotopes; alkenone temperatures; and tephra geochemistry from the northern Labrador Sea, off southwest Greenland, to investigate the deglacial response of the GIS and evaluate its implications for the North Atlantic deglacial development. The results imply that the southern GIS retreated in three successive stages: (1) early deglaciation of the East Greenland margins, by tephra-rich IRD that embrace Heinrich Event 1; (2) progressive retreat during Allerød culminating in major meltwater releases (d18O depletion of 1.2 per mil) at the Allerød-Younger Dryas transition (12.8-13.0 kyr B.P.); and (3) a final stage of glacial recession during the early Holocene (~9-11 kyr B.P.). Rather than indicating local temperatures of ambient surface water, the alkenones likely were transported to the core site by the Irminger Current. We attribute the timing of GIS retreat to the incursion of warm intermediate waters along the base of grounded glaciers and below floating ice shelves on the continental margin. The results lend support to the view that GIS meltwater presented a forcing factor for the Younger Dryas cooling.

Regional ocean model for the Mediterranean, links to model files

A regional ocean general circulation model of the Mediterranean is used to study the climate of the Last Glacial Maximum. The atmospheric forcing for these simulations has been derived from simulations with an atmospheric general circulation model, which in turn was forced with surface conditions from a coarse resolution earth system model. The model is successful in reproducing the general patterns of reconstructed sea surface temperature anomalies with the strongest cooling in summer in the northwestern Mediterranean and weak cooling in the Levantine, although the model underestimates the extent of the summer cooling in the western Mediterranean. However, there is a strong vertical gradient associated with this pattern of summer cooling, which makes the comparison with reconstructions complicated. The exchange with the Atlantic is decreased to roughly one half of its present value, which can be explained by the shallower Strait of Gibraltar as a consequence of lower global sea level. This reduced exchange causes a strong increase of salinity in the Mediterranean in spite of reduced net evaporation.

Late Quaternary siliciclastic/carbonate sedimentation in the Cabricorn Channel

A model is presented for hemipelagic siliciclastic and carbonate sedimentation during the last glacial-interglacial cycle in the Capricorn Channel, southern Great Barrier Reef (GBR). Stable isotope ratios, grainsize, carbonate content and mineralogy were analysed for seven cores in a depth transect from 166 to 2892 m below sea level (mbsl). Results show variations in the flux of terrigenous, neritic and pelagic sediments to the continental slope over the last sea level cycle. During the glacial lowstand terrigenous sediment influenced all the cores down to 2000 mbsl. The percentages of quartz and feldspar in the cores decreased with water depth, while the percentage of clay increased. X-ray diffraction analysis of the glacial lowstand clay mineralogy suggests that the siliciclastic sediment was primarily sourced from the Fitzroy River, which debouched directly into the northwest sector of the Capricorn Channel at this time. The cores also show a decrease in pelagic calcite and an increase in aragonite and high magnesium calcite (HMC) during the glacial. The influx of HMC and aragonite is most likely from reworking of coral reefs exposed on the continental shelf during the glacial, and also from HMC ooids precipitated at the head of the Capricorn Channel at this time. Mass accumulation rates (MARs) are high (13.5 g/cm**/kyr) during the glacial and peak at ~20 g/cm** 3/kyr in the early transgression (16-14 ka BP). MARs then decline with further sea level rise as the Fitzroy River mouth retreats from the edge of the continental shelf after 13.5 ka BP. MARs remain low (4 g/cm**3/kyr) throughout the Holocene highstand. Data for the Holocene highstand indicate there is a reduction in siliciclastic influx to the Capricorn Channel with little quartz and feldspar below 350 mbsl. However, fine-grained fluvial sediments, presumably from the Fitzroy River, were still accumulating on the mid slope down to 2000 mbsl. The proportion of pelagic calcite in the core tops increases with water depth, while HMC decreases, and is present only in trace amounts in cores below 1500 mbsl. The difference in the percentage of HMC in the deeper cores between the glacial and Holocene may reflect differences in supply or deepening of the HMC lysocline during the glacial. Sediment accumulation rates also vary between cores in the Capricorn Channel and do not show the expected exponential decrease with depth. This may be due to intermediate or deep water currents reworking the sediments. It is also possible that present bathymetry data are too sparse to detect the potential role that submarine channels may play in the distribution and accumulation of sediments. Comparison of the Capricorn Channel MARs with those for other mixed carbonate/siliciclastic provinces from the northeast margin of Australia indicates that peak MARs in the early transgression in the Capricorn Channel precede those from the central GBR and south of Fraser Island. The difference in the timing of the carbonate and siliciclastic MAR peaks along the northeast margin is primarily related to differences in the physiography and climate of the provinces. The only common trend in the MARs from the northeast margin of Australia is the near synchronicity of the carbonate and siliciclastic MAR peaks in individual sediment cores, which supports a coeval sedimentation model.

Radiocarbon age determination of 16 sites

The thermal diffusion enrichment apparatus in use in Amsterdam before 1967, has been rebuilt in the Groningen Radiocarbon Dating Laboratory. It has been shown to operate reliably and reproducibly. A reasonable agreement exists between the theoretical calculations and the experimental results. The 14C enrichment of a CO sample is deduced from the simultaneous mass 30 enrichment, which is measured with a mass spectrometer. The relation between both enrichments follows from a series of calibration measurements. The over-all accuracy in the enrichment is a few percent, equivalent to a few hundred years in age. The main problem in dating very old samples is their possible contamination with recent carbon. Generally, careful sample selection and rigorous pretreatment reduce sample contamination to an acceptable value. Also, it has been established that laboratory contamination, due to a memory effect in the combustion system and to impurities in the oxygen and nitrogen gas used for combustion, can be eliminated. A detailed analysis shows that the counter background in our set-up is almost exclusively caused by cosmic ray muons. The measurement of 28 early glacial samples, mostly from North-west Europe, has yielded a consistent set of ages. These indicate the existence of three early glacial interstadials; using the Weichselian definitions: Amersfoort starting at 68 200 ± 1100, Brørup at 64 400 ± 800 and Odderade at 60 500 ± 600 years BP. This 14C chronology shows good agreement with the Camp Century chronology and the dated palaeo sea levels. The discrepancy in the age of the early part of the Last Glacial on the 14C time scale and on that adopted for the deep-sea d18 record, must probably be attributed to the use of a generalized d18 curve and a wrong interpretation of this curve in terms of three Barbados terraces.