Age determination and stable isotope record of planktonic foraminifera from sediments of the Baffin Shelf

Stable isotopic values on planktonic foraminifera in a suite of cores from basins across the SE Baffin Shelf are used to extract a record of meltwater events during Termination I deglaciation. Resolution and Hatton basins lie on the SE Baffin Shelf at water depths > 500 m, seaward of major conduits for ice drainage from the eastern sector of the Laurentide Ice Sheet (LIS). Accelerator mass spectrometry 14C dates are used to constrain our chronology of events in ten cores. In Resolution Basin, three cores have 14C AMS dates on foraminifera of > 20 ka at their bases; whereas Hatton Basin cores terminate in sediments < 13 kyr. Sedimentation rates varied between 0.1 to 4.5 m/ka. Stable oxygen and carbon isotopic ratios were obtained on 146 samples of the planktonic foraminifera Neogloboquadrina pachyderma (Ehrenberg) sinistral, from seven of the ten cores. No evidence was found to indicate that test morphology or size affected delta18O. Between 7 and 13.5 ka the surface water on the shelf was on average 1 per mil lower than the open ocean signal. Significant temporal variations were found in both delta18O and delta13C. Evidence for significant low delta18O events occurred between 13 and 8 ka. The delta13C record from the planktonic foraminifera suggests a threefold division of events between 13 and 7 ka, with positive values between 10.8 and 13.0 ka, negative values between 9 and 10.8 ka, and positive values from 7 to 9 ka. The delta18O data suggest the presence of meltwater on the shelf some 3,000 years prior to the first late glacial dates on terrestrial deglaciation (at circa 10.4 ka). "Hudson Strait must be the real key to the importance of the calving process during deglaciation, because it is potentially the largest marine outlet for the Laurentide Ice Sheet and because it leads into the very center of the ice sheet.....the rates of calving through Hudson Strait during the period of initial ?18O rise unfortunately are unknown." W. F. Ruddiman (1987, p. 151)

Grain size and diatom valves in sediment cores KS82-30 and KS82-31

Diatoms were studied quantitatively in six latest Quaternary (~70 kyr B.P. to Recent) piston cores from the westernmost Mediterranean, the Alboran Basin, and the Atlantic region immediately to the west of the Straits of Gibraltar. The Atlantic cores completely lack diatoms. In the Alboran Basin, diatoms are common from late Stage 3 (~27.5 kyr B.P.) to Termination lb (9 kyr B.P.) and in Recent core tops, but are absent in the other latest Quaternary intervals. Maximum accumulation of diatoms and highest abundance of species normally in sediments associated with increased productivity occurred during the latest Quaternary deglaciation, in the first phase of Termination I (~14.8 kyr B.P.). In the modern Alboran Basin, a region of high biological productivity occurs immediately east of the Gibraltar Straits. This high productivity results from upwelling associated with the interaction between the Atlantic inflow and the bottom topography near the Spanish coast. The upwelled nutrient-rich waters are then advected to the east and southeast by the surficial anticyclonic gyral circulation. Late Quaternary variations in diatom abundance are considered to reflect changes in this upwelling intensity with highest diatom abundances inferred to result from increased upwelling associated with an intensification of the anticyclonic gyral circulation. Highest inferred upwelling rates occurred during the first phase of latest Quaternary deglaciation. It is possible that an intensification of circulation within the Mediterranean Basin as a whole occurred from late Stage 3 to mid Termination I because widespread hiatus formation has been reported at this time in the Straits of Sicily due to an increase in the formation of intermediate waters. Diatoms were not preserved in other latest Quaternary intervals due to insufficient productivity to counterbalance their dissolution.

(Appendix A) Relative abundance of Emiliania huxleyi in sediment cores of the North Atlantic

The coccolithophore species Emiliania huxleyi is characterized by a wide range of sizes, which can be easily distinguished in the light microscope. In this study we have quantified the abundance of large (coccoliths > 4 µm in maximum length) E. huxleyi specimens during the last 25 kyr in sedimentary records from eleven cores and drill sites in the NE Atlantic and W Mediterranean Sea, to prove its usefulness in the reconstruction of water mass dynamics and biostratigraphic potential. During the Last Glacial Maximum this large form, a cold-water indicator, was common in the NE Atlantic and Mediterranean, and its regional variation in abundance indicates a displacement of the climatic zones southwards in agreement with the development of ice sheets and sea ice in the Northern Hemisphere during this period. On the other hand, the gradient between northern and southern surface water masses in the Subtropical Gyre appears to have been more pronounced than at present, while the Portugal and Canary Currents were more intense. In the western Mediterranean basin temperatures were cooler than in the adjacent Atlantic, provoking a quasi-endemism of these specimens until the end of Heinrich Event 1. This may have been due to a restriction in the communication between the Atlantic and Mediterranean through the Strait of Gibraltar, the arrival of cold surface water and the amplification of cooling after the development of ice sheets in the Northern Hemisphere. During the deglaciation, large E. huxleyi specimens decreased in abundance at medium and low latitudes, but were still numerous close to the Subarctic region during the Holocene. In transitional waters this decrease to present day abundances occurred after Termination Ib. The abrupt change in abundance of this large E. huxleyi form is proposed as a new biostratigraphic event to characterize the Holocene in mid- to low-latitude water masses in the North Atlantic, although this horizon seems to be diachronous by 5 kyr from tropical to subarctic regions, in agreement with the gradual onset of warm conditions.

Paleotemperature reconstructions on sediment core GIK161603-3

Different proxies for sea surface temperature (SST) often exhibit divergent trends for deglacial warming in tropical regions, hampering our understanding of the phase relationship between tropical SSTs and continental ice volume at glacial terminations. To reconcile divergent SST trends, we report reconstructions of two commonly used paleothermometers (the foraminifera G. ruber Mg/Ca and the alkenone unsaturation index) from a marine sediment core collected in the southwestern tropical Indian Ocean encompassing the last 37,000 years. Our results show that SSTs derived from the alkenone unsaturation index (UK'37) are consistently warmer than those derived from Mg/Ca by ~2-3°C except for the Heinrich Event 1. In addition, the initial timing for the deglacial warming of alkenone SST started at ~15.6 ka, which lags behind that of Mg/Ca temperatures by 2.5 kyr. We argue that the discrepancy between the two SST proxies reflects seasonal differences between summer and winter rather than post-depositional processes or sedimentary biases. The UK'37 SST record clearly mimics the deglacial SST trend recorded in the North Atlantic region for the earlier part of the termination, indicating the early deglacial warming trend attributed to local summer temperatures was likely mediated by changes in the Atlantic Meridional Overturning Circulation at the onset of the deglaciation, In contrast, the glacial to interglacial SST pattern recorded by G. ruber Mg/Ca probably reflects cold season SSTs. This indicates that the cold season SSTs was likely mediated by climate changes in the southern hemisphere, as it closely tracks the Antarctic timing of deglaciation. Therefore our study reveals that the tropical southwestern Indian Ocean seasonal SST was closely linked to climate changes occurring in both hemispheres. The austral summer and winter recorded by each proxy is further supported with seasonal SST trends modeled by AOGCMs for our core site. Our interpretation that the alkenone and Mg/Ca SSTs are seasonally biased may also explain similar proxy mismatches observed in other tropical regions at the onset of the last termination.

Age determination of sediment core MD07-3088

The Southern Ocean plays a prominent role in the Earth's climate and carbon cycle. Changes in the Southern Ocean circulation may have regulated the release of CO2 to the atmosphere from a deep-ocean reservoir during the last deglaciation. However, the path and exact timing of this deglacial CO2 release are still under debate. Here we present measurements of deglacial surface reservoir 14C age changes in the eastern Pacific sector of the Southern Ocean, obtained by 14C dating of tephra deposited over the marine and terrestrial regions. These results, along with records of foraminifera benthic-planktic 14C age and d13C difference, provide evidence for three periods of enhanced upwelling in the Southern Ocean during the last deglaciation, supporting the hypothesis that Southern Ocean upwelling contributed to the deglacial rise in atmospheric CO2. These independently dated marine records suggest synchronous changes in the Southern Ocean circulation and Antarctic climate during the last deglaciation.

Nd isotopes of sediments, foraminifera and fish debris of the western North Atlantic

Understanding changes in ocean circulation during the last deglaciation is crucial to unraveling the dynamics of glacial-interglacial and millennial climate shifts. We used neodymium isotope measurements on postdepositional iron-manganese oxide coatings precipitated on planktonic foraminifera to reconstruct changes in the bottom water source of the deep western North Atlantic at the Bermuda Rise. Comparison of our deep water source record with overturning strength proxies shows that both the deep water mass source and the overturning rate shifted rapidly and synchronously during the last deglacial transition. In contrast, any freshwater perturbation caused by Heinrich event 1 could have only affected shallow overturning. These findings show how changes in upper-ocean overturning associated with millennial-scale events differ from those associated with whole-ocean deglacial climate events.

Mg/Ca values for planktonic foraminifera N. pachyderma and G. bulloides from ocean sediment core MD02-2496 on the Vancouver margin, Northeastern Pacific Ocean

A new, high-resolution planktonic foraminiferal Mg/Ca-based ocean temperature record has been generated for deep sea core MD02-2496, sited offshore of Vancouver Island, Western Canada during the last deglaciation (21-12 ka). The relationship between Cordilleran Ice Sheet (CIS) retreat and changing regional ocean temperatures has been reconstructed through glaciomarine sediments in MD02-2496 that capture tidewater glacier response to surface ocean thermal forcing. At CIS maximum extent, the marine margin of the ice sheet advanced onto the continental shelf. During this interval, ocean temperatures recorded by surface ocean dwelling Globigerina bulloides remained a relatively constant ~7.5°C while subsurface dwelling Neogloboquadrina pachyderma (s.) recorded temperatures of ~5°C. These ocean temperatures were sufficiently warm to induce significant melt along the tidewater ice terminus similar to modern Alaskan tidewater glacial systems. During the deglacial retreat of the CIS, the N. pachyderma temperature record shows two distinct warming steps of ~2 and 2.5°C between 17.2-16 and 15.5-14 ka respectively, coincident with ice rafting events from the CIS, while G. bulloides records an ~3°C warming from 15 to14 ka. We hypothesize that submarine melting resulting from relatively warm ocean temperatures was an important process driving ice removal from CIS tidewater glaciers during the initial stages of deglaciation.

Geochemistry of surface sediment samples from the equatorial Atlantic south of the Azores

In order to recognize the distribution of hydrothermal tracers south of the Azores, a series of cores has been sampled during the GEOFAR cruise. The distribution of the Mn concentrations shows that the hydrothermal influence is maximum within and to the north-west of the Lucky Strike segment. North of the East-Azores Fracture Zone the sediments are enriched in Ba which could be derived from different sources. The chemical composition of the interstitial water shows that water advection is limited. Mn, Cu, Ni fluxes evaluated in one site have increased during isotopic stages 4 and 2 and the deglaciation.

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.

Revised calibration of the MBT-CBT paleotemperature proxy

The MBT-CBT proxy for the reconstruction of paleotemperatures and past soil pH is based on the distribution of branched glycerol dialkyl glycerol tetraether (brGDGT) membrane lipids. The Methylation of Branched Tetraether (MBT) and the Cyclisation of Branched Tetraether (CBT) indices were developed to quantify these distributions, and significant empirical relations between these indices and annual mean air temperature (MAT) and/or soil pH were found in a large data set of soils. In this study, we extended this soil dataset to 278 globally distributed surface soils. Of these soils, 26% contains all nine brGDGTs, while in 63% of the soils the seven most common brGDGTs were detected, and the latter were selected for calibration purposes. This resulted in new transfer functions for the reconstruction of pH based on the CBT index: pH = 7.90-1.97 × CBT (r**2 = 0.70; RMSE = 0.8; n = 176), as well as for MAT based on the CBT index and methylation index based on the seven most abundant GDGTs (defined as MBT'): MAT = 0.81-5.67 × CBT + 31.0 × MBT' (r**2 = 0.59; RMSE = 5.0 °C; n = 176). The new transfer function for MAT has a substantially lower correlation coefficient than the original equation (r**2 = 0.77). To investigate possible improvement of the correlation, we used our extended global surface soil dataset to statistically derive the indices that best describe the relations of brGDGT composition with MAT and soil pH. These new indices, however, resulted in only a relatively minor increase in correlation coefficients, while they cannot be explained straightforwardly by physiological mechanisms. The large scatter in the calibration cannot be fully explained by local factors or by seasonality, but MAT for soils from arid regions are generally substantially (up to 20 °C) underestimated, suggesting that absolute brGDGT-based temperature records for these areas should be interpreted with caution. The applicability of the new MBT'-CBT calibration function was tested using previously published MBT-CBT-derived paleotemperature records covering the last deglaciation in Central Africa and East Asia, the Eocene-Oligocene boundary and the Paleocene-Eocene thermal maximum. The results show that trends remain similar in all records, but that absolute temperature estimates and the amplitude of temperature changes are lower for most records, and generally in better agreement with independent proxy data.

Stable isotope composition and ice rafted debris in early Late Pliocene sediments from Maud Rise, South Atlantic

We have conducted an integrated study of ice-rafted debris (IRD) and oxygen isotopes (measured on Cibicides, Globigerina bulloides, and Neogloboquadrina pachyderma, using identical samples). We used samples from the early Late Pliocene Gauss Chron from ODP Site 114-704 on the Meteor Rise in the subantarctic South Atlantic. During the early Gauss Chron, the oxygen isotopic ratios are generally up to 0.5?-0.6? less than their respective Holocene values. The lowest values in this record can accommodate a warming of about 2.5°C or a sea-level rise of about 50 m, but not both, and probably result from some warming and a small reduction in global ice volume. Starting with isotope stage MG2 [ 3.23 Ma on the Berggren et al. ( 1985) time scale; 3.38 on the Shackleton et al. ( 1995b) time scale] oxygen-isotopic values generally increase (and oscillate about a Holocene mean). The first significant IRD appears at the same time. There is a subsequent increase in IRD amounts upsection. In order to reach the site, this material must have been transported by large, tabular icebergs derived from Antarctic ice shelves or ice tongues, similar to occasional, large modern icebergs. This combined record suggests strongly that the Antarctic ice sheet was essentially intact; some warming at the drill site is indicated, but not a major reduction in ice-volume on Antarctica.

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.

Silt grain-size analysis of ODP holes from the western equatorial Atlantic

Grain-size records of the terrigenous and calcareous silt fraction, preservation of planktic foraminifera, and benthic foraminiferal stable-isotope data (delta13C, delta18O values of C. wuellerstorfi) at ODP Site 927 on the Ceará Rise (5°27.7'N, 44°28.8'W), are used to reconstruct variations in the history of bottom current strength, ventilation, and carbonate corrosiveness of deep waters during the time interval from 0.8 to 0.3 Ma. Glacial periods are characterized by generally smaller mean sizes of the terrigenous sortable silt fraction (mean(SS)), lower delta13C values, and poorer preservation of planktic foraminifera compared to interglacials. This indicates lower bottom current speeds, larger nutrient contents and more corrosive deep water. By contrast, larger mean(SS) sizes, higher delta13C values, and well preserved planktic foraminifera indicate strong circulation and a well ventilated deep-water mass during interglacials. The observed changes are most likely related to the weakening and strengthening of circulation of Lower North Atlantic Deep Water (LNADW). Cross-spectral analysis between the mean(SS) and benthic delta18O records reveals that minima in mean(SS) occur about 7.6 k.y. after the maximum in ice volume. This indicates a considerable lag time between ice-shield induced changes in LNADW production and subsequent changes in the velocity of LNADW flow in the western equatorial Atlantic. Striking changes in bottom current speed occur regularly during glacial to interglacial transitions. Extremely fine mean(SS) minima point to an almost complete shutdown of bottom current vigor in response to a cessation of LNADW production caused by an enhanced melt water release during the initial phases of deglaciation. However, each of the fine minima extremes is followed by a rapid shift to very high mean(SS) values that indicate strong bottom currents, and hence, vigorous LNADW flow during the early interglacials. After the onset of glacial Stage 12, generally poorer carbonate preservation and higher variability is registered. This coincides with a global decrease in carbonate preservation during the mid-Brunhes (mid-Brunhes dissolution event). Detailed grain-size analysis of the calcareous fine fraction (<63 µm) revealed a considerable reduction of particles in the fraction from 7 to 63 µm during periods of enhanced dissolution. This indicates a preferential dissolution of larger planktic foraminiferal fragments which leads to an enrichment of coccoliths in the calcareous fine fraction.

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.