Barium measured on water bottle samples during POLARSTERN cruise ARK-XXII/2 (SPACE)

Dissolved barium has been shown to have the potential to distinguish Eurasian from North American (NA) river runoff. As part of the ARK-XXII/2 Polarstern expedition in summer 2007, Ba was analyzed in the Barents, Kara, Laptev seas, and the Eurasian Basins as well as the Makarov Basin up to the Alpha and Mendeleyev Ridges. By combining salinity, d18O and initial phosphate corrected for mineralization with oxygen (PO4*) or N/P ratios we identified the water mass fractions of meteoric water, sea ice meltwater, and marine waters of Atlantic as well as Pacific origin in the upper water column. In all basins inside the lower halocline layer and the Arctic intermediate waters we find Ba concentrations close to those of the Fram Strait branch of the lower halocline (41-45 nM), reflecting the composition of the incoming Atlantic water. A layer of upper halocline water (UHW) with higher Ba concentrations (45-55 nM) is identified in the Makarov Basin. Atop of the UHW, the Surface Mixed Layer (SML), including the summer and winter mixed layers, has high concentrations of Ba (58-67 nM). In the SML of the investigated area of the central Arctic the meteoric fraction can be identified by assuming a conservative behavior of Ba to be primarily of Eurasian river origin. However, in productive coastal regions biological removal compromises the use of Ba to distinguish between Eurasian and NA rivers. As a consequence, the NA river water fraction is underestimated in productive surface waters or waters that have passed a productive region, whereas this fraction is overestimated in subsurface waters containing remineralised Ba, particularly when these waters have passed productive shelf regions. Especially in the Laptev Sea and small regions in the Barents Sea, Ba concentrations are low in surface waters. In the Laptev Sea exceptionally high Ba concentrations in shelf bottom waters indicate that Ba is removed from surface waters to deep waters by biological activity enhanced by increasing ice-free conditions as well as by scavenging by organic matter of terrestrial origin. We interpret high Ba concentrations in the UHW of the Makarov Basin to result from enrichment by remineralisation in bottom waters on the shelf of the Chukchi Sea and therefore the calculated NA runoff is an artefact. We conclude that no NA runoff can be demonstrated unequivocally anywhere during our expedition with the set of tracers considered here. Small contributions of NA runoff may have been masked by Ba depletion and could only be resolved by supportive tracers on the uptake history. We thus suggest that Ba has to be used with care as it can put limits but not yield quantitative water mass distributions. Only if the extra Ba inputs exceed the cumulative biological uptake the signal can be unequivocally attributed to NA runoff.

Dissolved iron measurements from 32 stations in the Central Arctic Ocean and Arctic Shelf Sea during POLARSTERN expedition ARK-XXII/2

Concentrations of dissolved (<0.2 µm) Fe (DFe) in the Arctic shelf seas and in the surface waters of the central Arctic Ocean are presented. In the Barents and Kara seas, near-surface DFe minima indicate depletion of DFe by phytoplankton growth. Below the surface, lower DFe concentrations in the Kara Sea (~0.4-0.6 nM) than in the Barents Sea (~0.6-0.8 nM) likely reflect scavenging removal or biological depletion of DFe. Very high DFe concentrations (>10 nM) in the bottom waters of the Laptev Sea shelf may be attributed to either sediment resuspension, sinking of brine or regeneration of DFe in the lower layers. A significant correlation (R2 = 0.60) between salinity and DFe is observed. Using d18O, salinity, nutrients and total alkalinity data, the main source for the high (>2 nM) DFe concentrations in the Amundsen and Makarov Basins is identified as (Eurasian) river water, transported with the Transpolar Drift (TPD). On the North American side of the TPD, the DFe concentrations are low (<0.8 nM) and variations are determined by the effects of sea-ice meltwater, biological depletion and remineralization and scavenging in halocline waters from the shelf. This distribution pattern of DFe is also supported by the ratio between unfiltered and dissolved Fe (high (>4) above the shelf and low (<4) off the shelf).

Biomarker and x-ray diffraction data from Site U1313 (MIS 16 - 9)

A reconstruction of Milankovitch to millennial-scale variability of sea-surface temperature (SST) and sea-surface productivity in the Pleistocene mid-latitude North Atlantic Ocean (MIS 16-9) and its relationship to ice sheet instability was carried out on sediments from IODP Site U1313. This reconstruction is based on alkenone and n-alkane concentrations, Uk37' index, total organic carbon (TOC) and carbonate contents, X-Ray diffraction (XRD) data, magnetic susceptibility, and accumulation rates. Increased input of ice-rafted debris (IRD) occurred during MIS 16, 12, and 10, characterized by high concentrations of dolomite, quartz, and feldspars and elevated accumulation rates of terrigenous matter. Minimum input values of terrigenous matter, on the other hand, were determined for MIS 13 and 11. Peak values of dolomite, coinciding with quartz, plagioclase, and kalifeldspar peaks and maxima in long-chain n-alkanes indicative for land plants, are interpreted as Heinrich-like Events related to sudden instability of the Laurentide Ice Sheet during early and late (deglacial) phases of the glacials. The coincidence of increased TOC values with elevated absolute concentrations of alkenones suggest increased glacial productivity, probably due to a more southern position of the Polar Front. Alkenone-based SST reached absolute maxima of about 19°C during MIS 11.3 and absolute minima of <10°C during MIS 12 and 10. Within MIS 11, prominent cooling events (MIS 11.22 and 11.24) occurred. The absolute SST minima recorded directly before and after the glacial maxima MIS 10.2 and 12.2, are related to Heinrich-like Event meltwater pulses, as supported by the coincidence of SST minima and maxima in C37:4 alkenones and dolomite. These sudden meltwater pulses - especially during Terminations IV and V - probably caused a collapse of phytoplankton productivity as indicated by the distinct drop in alkenone concentrations. Ice-sheet disintegration and subsequent surges and outbursts of icebergs and meltwater discharge may have been triggered by increased insolation in the Northern High Latitudes.

Physical properties and alpha-/gamma-HCH concentrations in sea-ice and snow samples obtained during a CCGS Amundsen cruise, eastern Beaufort Sea

The alpha- and gamma-hexachlorocyclohexanes (HCHs) are being scavenged from the atmosphere by falling snow, with the average total scavenging ratios (WT) of 3.8 x 10**4 and 9.6 x 10**3, respectively. After deposition, HCH snow concentrations can decrease by 40% because of snowpack ventilation and increase by 50% because of upward migration of brine from the ice. HCH vertical distribution in sufficiently cold winter sea ice, which maintains brine volume fractions <5%, reflects the ice growth history. Initially, the entrapment of brine (and HCHs) in ice depends on the rates of ice growth and desalination. However, after approximately the first week of ice formation, ice growth rate becomes dominant. Deviations of HCH concentrations from the values predicted by the ice bulk salinity (rate of brine entrapment) can be explained by spatial variability of HCHs in surface water. HCH burden in the majority of the ice column remains locked throughout most of the season until the early spring when snow meltwater percolates into the ice, delivering HCHs to the upper ocean via desalination by flushing. Percolation can lead to an increase in alpha- and gamma-HCH in the sea ice by up to 2%-18% and 4%-32%, respectively.

Stable carbon and oxygen isotope record of IODP Site 306-U1313, MIS 23-19

Since the seminal work by Hays et al. (1976), a plethora of studies has demonstrated a correlation between orbital variations and climatic change. However, information on how changes in orbital boundary conditions affected the frequency and amplitude of millennial-scale climate variability is still fragmentary. The Marine Isotope Stage (MIS) 19, an interglacial centred at around 785 ka, provides an opportunity to pursue this question and test the hypothesis that the long-term processes set up the boundary conditions within which the short-term processes operate. Similarly to the current interglacial, MIS 19 is characterised by a minimum of the 400-kyr eccentricity cycle, subdued amplitude of precessional changes, and small amplitude variations in insolation. Here we examine the record of climatic conditions during MIS 19 using high-resolution stable isotope records from benthic and planktonic foraminifera from a sedimentary sequence in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313) in order to assess the stability and duration of this interglacial, and evaluate the climate system's response in the millennial band to known orbitally induced insolation changes. Benthic and planktonic foraminiferal d18O values indicate relatively stable conditions during the peak warmth of MIS 19, but sea-surface and deep-water reconstructions start diverging during the transition towards the glacial MIS 18, when large, cold excursions disrupt the surface waters whereas low amplitude millennial scale fluctuations persist in the deep waters as recorded by the oxygen isotope signal. The glacial inception occurred at ~779 ka, in agreement with an increased abundance of tetra-unsaturated alkenones, reflecting the influence of icebergs and associated meltwater pulses and high-latitude waters at the study site. After having combined the new results with previous data from the same site, and using a variety of time series analysis techniques, we evaluate the evolution of millennial climate variability in response to changing orbital boundary conditions during the Early-Middle Pleistocene. Suborbital variability in both surface- and deep-water records is mainly concentrated at a period of ~11 kyr and, additionally, at ~5.8 and ~3.9 kyr in the deep ocean; these periods are equal to harmonics of precession band oscillations. The fact that the response at the 11 kyr period increased over the same interval during which the amplitude of the response to the precessional cycle increased supports the notion that most of the variance in the 11 kyr band in the sedimentary record is nonlinearly transferred from precession band oscillations. Considering that these periodicities are important features in the equatorial and intertropical insolation, these observations are in line with the view that the low-latitude regions play an important role in the response of the climate system to the astronomical forcing. We conclude that the effect of the orbitally induced insolation is of fundamental importance in regulating the timing and amplitude of millennial scale climate variability.

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.

Age determination and model of sediment core RAPiD-15-4P

Greenland ice core records indicate that the last deglaciation (~7-21 ka) was punctuated by numerous abrupt climate reversals involving temperature changes of up to 5°C-10°C within decades. However, the cause behind many of these events is uncertain. A likely candidate may have been the input of deglacial meltwater, from the Laurentide ice sheet (LIS), to the high-latitude North Atlantic, which disrupted ocean circulation and triggered cooling. Yet the direct evidence of meltwater input for many of these events has so far remained undetected. In this study, we use the geochemistry (paired Mg/Ca-d18O) of planktonic foraminifera from a sediment core south of Iceland to reconstruct the input of freshwater to the northern North Atlantic during abrupt deglacial climate change. Our record can be placed on the same timescale as ice cores and therefore provides a direct comparison between the timing of freshwater input and climate variability. Meltwater events coincide with the onset of numerous cold intervals, including the Older Dryas (14.0 ka), two events during the Allerød (at ~13.1 and 13.6 ka), the Younger Dryas (12.9 ka), and the 8.2 ka event, supporting a causal link between these abrupt climate changes and meltwater input. During the Bølling-Allerød warm interval, we find that periods of warming are associated with an increased meltwater flux to the northern North Atlantic, which in turn induces abrupt cooling, a cessation in meltwater input, and eventual climate recovery. This implies that feedback between climate and meltwater input produced a highly variable climate. A comparison to published data sets suggests that this feedback likely included fluctuations in the southern margin of the LIS causing rerouting of LIS meltwater between southern and eastern drainage outlets, as proposed by Clark et al. (2001, doi:10.1126/science.1062517).

Stable isotope record of benthic foraminifera from the Gulf of California

A high-resolution, accelerator radiocarbon dated climate record of the interval 8,000-18,000 years B.P. from Deep Sea Drilling Project site 480 (Guaymas Basin, Gulf of California) shows geochemical and lithological oscillations of oceanographic and climatic significance during deglaciation. Nonlaminated sediments are associated with cooler climatic conditions during the late glacial (up to 13,000 years B.P.), and from 10,300 to 10,800 years B.P., equivalent to the Younger Dryas event of the North Atlantic region. We propose that the changes from laminated (varved) to nonlaminated sediments resulted from increased oxygen content in Pacific intermediate waters during the glacial and the Younger Dryas episodes, and that the forcing for the latter event was global in scope. Prominent events of low delta18O are recorded in benthic foraminifera from 8,000 to 10,000 and at 12,000 years B.P.; evidence for an earlier event between 13,500 and 15,000 years B.P. is weaker. Maximum delta18O is found to have occurred 10,500, 13,500, and 15,000 years ago (and beyond). Oxygen isotopic variability most likely reflects changing temperature and salinity characteristics of Pacific waters of intermediate depth during deglaciation or environmental changes within the Gulf of California region. Several lines of evidence suggest that during deglaciation the climate of the American southwest was marked by increased precipitation that could have lowered salinity in the Gulf of California. Recent modelling studies show that cooling of the Gulf of Mexico due to glacial meltwater injection, which is believed to have occurred at least twice during deglaciation, would have resulted in increased precipitation with respect to evaporation in the American southwest during summertime. The timing of deglacial events in the Gulf of Mexico and the Gulf of California supports such an atmospheric teleconnection.

Investigation of the seasonality of N. pachyderma in the North Atlantic Ocean during Heinrich Stadial 1 using an ecosystem modeling approach

Fossil shells of planktonic foraminifera serve as the prime source of information on past changes in surface ocean conditions. Because the population size of planktonic foraminifera species changes throughout the year, the signal preserved in fossil shells is biased towards the conditions when species production was at its maximum. The amplitude of the potential seasonal bias is a function of the magnitude of the seasonal cycle in production. Here we use a planktonic foraminifera model coupled to an ecosystem model to investigate to what degree seasonal variations in production of the species Neogloboquadrina pachyderma may affect paleoceanographic reconstructions during Heinrich Stadial 1 (~18-15 cal. ka B.P.) in the North Atlantic Ocean. The model implies that during Heinrich Stadial 1 the maximum seasonal production occurred later in the year compared to the Last Glacial Maximum (~21-19 cal. ka B.P.) and the pre-industrial era north of 30 ºN. A diagnosis of the model output indicates that this change reflects the sensitivity of the species to the seasonal cycle of sea-ice cover and food supply, which collectively lead to shifts in the modeled maximum production from the Last Glacial Maximum to Heinrich Stadial 1 by up to six months. Assuming equilibrium oxygen isotopic incorporation in the shells of N. pachyderma, the modeled changes in seasonality would result in an underestimation of the actual magnitude of the meltwater isotopic signal recorded by fossil assemblages of N. pachyderma wherever calcification is likely to take place.