Glacial North Atlantic: Sea surface conditions reconstructed by GLAMAP 2000

The response of the tropical ocean to global climate change and the extent of sea ice in the glacial nordic seas belong to the great controversies in paleoclimatology. Our new reconstruction of peak glacial sea surface temperatures (SSTs) in the Atlantic is based on census counts of planktic foraminifera, using the Maximum Similarity Technique Version 28 (SIMMAX-28) modern analog technique with 947 modern analog samples and 119 well-dated sediment cores. Our study compares two slightly different scenarios of the Last Glacial Maximum (LGM), the Environmental Processes of the Ice Age: Land, Oceans, Glaciers (EPILOG), and Glacial Atlantic Ocean Mapping (GLAMAP 2000) time slices. The comparison shows that the maximum LGM cooling in the Southern Hemisphere slightly preceeded that in the north. In both time slices sea ice was restricted to the north western margin of the nordic seas during glacial northern summer, while the central and eastern parts were ice-free. During northern glacial winter, sea ice advanced to the south of Iceland and Faeroe. In the central northern North Atlantic an anticyclonic gyre formed between 45° and 60°N, with a cool water mass centered west of Ireland, where glacial cooling reached a maximum of >12°C. In the subtropical ocean gyres the new reconstruction supports the glacial-to-interglacial stability of SST as shown by CLIMAP Project Members (CLIMAP) [1981]. The zonal belt of minimum SST seasonality between 2° and 6°N suggests that the LGM caloric equator occupied the same latitude as today. In contrast to the CLIMAP reconstruction, the glacial cooling of the tropical east Atlantic upwelling belt reached up to 6°-8°C during Northern Hemisphere summer. Differences between these SIMMAX-based and published U37[k]- and Mg/Ca-based equatorial SST records are ascribed to strong SST seasonalities and SST signals that were produced by different planktic species groups during different seasons.

(Table 1) Age determination of North Atlantic sediment cores

High-resolution sediment cores from the Vøring Plateau, the North Iceland shelf, and the East Greenland shelf have been studied to investigate the stability of major surface currents in the Nordic Seas during the Holocene. Results from diatom assemblages and reconstructed sea-surface temperatures (SSTs) indicate a division of the Holocene into three periods: the Holocene Climate Optimum (9500-6500 calendar (cal) years BP), the Holocene Transition Period (6500-3000 cal years BP) and the Cool Late Holocene Period (3000-0 cal years BP). The overall climate development is in step with the decreasing insolation on the Northern Hemisphere, but regional differences occur regarding both timing and magnitude of SST changes. Sites under the direct influence of the Norwegian Atlantic Current and the Irminger Current indicate SST cooling of 4-5°C from early Holocene to present, compared to 2°C recorded under the East Greenland Current. Superimposed on the general Holocene cooling trend, there is a high-frequency SST variability, which is in the order of 1-1.5°C for the Vøring Plateau and the East Greenland shelf and 2.5-3°C on the North Iceland shelf.

Investigation of N. pachyderma from the North Atlantic

The shells of the planktonic foraminifer Neogloboquadrina pachyderma have become a classical tool for reconstructing glacial-interglacial climate conditions in the North Atlantic Ocean. Palaeoceanographers utilize its left- and right-coiling variants, which exhibit a distinctive reciprocal temperature and water mass related shift in faunal abundance both at present and in late Quaternary sediments. Recently discovered cryptic genetic diversity in planktonic foraminifers now poses significant questions for these studies. Here we report genetic evidence demonstrating that the apparent 'single species' shell-based records of right-coiling N. pachyderma used in palaeoceanographic reconstructions contain an alternation in species as environmental factors change. This is reflected in a species-dependent incremental shift in right-coiling N. pachyderma shell calcite d18O between the Last Glacial Maximum and full Holocene conditions. Guided by the percentage dextral coiling ratio, our findings enhance the use of d18O records of right-coiling N. pachyderma for future study. They also highlight the need to genetically investigate other important morphospecies to refine their accuracy and reliability as palaeoceanographic proxies.

Age models and stable isotopes of sediment cores from the northern seas

Oxygen and carbon isotope measurements were carried out on tests of planktic foraminifers N. pachyderma (sin.) from eight sediment cores taken from the eastern Arctic Ocean, the Fram Strait, and the lceland Sea, in order to reconstruct Arctic Ocean and Norwegian-Greenland Sea circulation patterns and ice covers during the last 130,000 years. In addition, the influence of ice, temperature and salinity effects on the isotopic signal was quantified. Isotope measurements on foraminifers from sediment surface samples were used to elucidate the ecology of N. pachyderma (sin.). Changes in the oxygen and carbon isotope composition of N. pachyderma (sin.) from sediment surface samples document the horizontal and vertical changes of water mass boundaries controlled by water temperature and salinity, because N. pachyderma (sin.) shows drastic changes in depth habitats, depending on the water mass properties. It was able to be shown that in the investigated areas a regional and spatial apparent increase of the ice effect occurred. This happened especially during the termination I by direct advection of meltwaters from nearby continents or during the termination and in interglacials by supply of isotopically light water from rivers. A northwardly proceeding overprint of the 'global' ice effect, increasing from the Norwegian-Greenland Sea to the Arctic Ocean, was not able to be demonstrated. By means of a model the influence of temperature and salinity on the global ice volume signal during the last 130,000 years was recorded. In combination with the results of this study, the model was the basis for a reconstruction of the paleoceanographic development of the Arctic Ocean and the Norwegian-Greenland Sea during this time interval. The conception of a relatively thick and permanent sea ice cover in the Nordic Seas during glacial times should be replaced by the model of a seasonally and regionally highly variable ice cover. Only during isotope stage 5e may there have been a local deep water formation in the Fram Strait.

Physical oceanography during POLARSTERN cruise ARK-XIX/4b

Exchanges between the North Atlantic and the Arctic Ocean result in the most dramatic water mass conversions in the World Ocean: warm and saline Atlantic waters, flowing through the Nordic Seas into the Arctic Ocean, are modified by cooling, freezing and melting to become shallow fresh waters, ice and saline deep waters. The outflow from the Nordic Seas to the south provides the initial driving of the global thermohaline circulation cell. Knowledge of these fluxes and understanding of the modification processes is a major prerequisite for the quantification of the rate of overturning within the large circulation cells of the Arctic and the Atlantic Oceans, and is also a basic requirement for understanding the role of these ocean areas in climate variability on interannual to decadal time scales. The Fram Strait represents the only deep connection between the Arctic Ocean and the Nordic Seas. Just as the freshwater transport from the Arctic Ocean is of major influence on convection in the Nordic Seas and further south, the transport of warm and saline Atlantic water affects the water mass characteristics in the Arctic Ocean which has consequences for the internal circulation and possibly influences also ice and atmosphere. The West Spitsbergen Current carrying Atlantic Water northward. The East Greenland Current, carrying water from the Arctic Ocean southwards has a concentrated core above the continental slope. It is our aim to measure the oceanic fluxes through Fram Strait and to determine their variability in seasonal to decadal time scales. 53 CTD profiles were taken at 51 stations. Two CTD systems from Sea-Bird Electronics Inc SBE911+ were used. Mainly SN 561 with duplicate T and C sensors (temperature sensors SBE3, SN 2685 and 2678, conductivity sensors SBE4, SN 2325 and 2618 and pressure sensor Digiquartz 410K-105 SN 75659) was in service. For the control of the temperature sensors a SBE35 RT digital reversing thermometer, SN 27 was applied. The CTD was connected to a SBE32 Carousel Water Sampler, SN 273 (24 12-liter bottles). For 3 CTD-Stations (726-3, 727-1, 728-1) the Sea-Bird 911+ probe SN 485 was used with temperature sensor SBE3 SN 2460, conductivity sensor SBE4 SN 2054, pressure sensor Digiquartz 410K SN 68997 and the SBE32 Carousel Water Sampler SN 202. Additionally Benthos Altimeters Model 2110-2, SN 189 and SN 208 and Wetlabs C-Star Transmissiometers SN 403 and SN 267 were mounted on the carousels. During the cruise a total number of 184 water samples were analysed with a Guildline Autosal 8400B salinometer, and IAPSO standard seawater batch number P141, K=0.99993. 20 salinity samples were brought back to AWI for onshore analysis. The CTD sensors were calibrated before and after the cruise by Sea-Bird Electronics.

Swath sonar bathymetry during POLARSTERN cruise ARK-VII/3a (PS17) with links to multibeam raw data files

In multibeam echosounder and subbottom profiler data acquired during R/V Polarstern cruise ARK-VII/3a from the Hovgaard Ridge (Fram Strait), we found evidence for very deep (>1200 m) iceberg scouring. Five elongated seafloor features have been detected that are interpreted to be iceberg scours. The scours are oriented in north-south/south-north direction and are about 15 m deep, 300 m wide, and 4 km long crossing the entire width of the ridge. They are attributed to multiple giant paleo-icebergs that most probably left the Arctic Ocean southward through Fram Strait. The huge keel depths are indicative of ice sheets extending into the Arctic Ocean being at least 1200 m thick at the calving front during glacial maxima. The deep St. Anna Trough or grounded ice observed at the East Siberian Continental Margin are likely source regions of these icebergs that delivered freshwater to the Nordic Seas.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

(Table S1) Age determination of Atlantic Ocean sediment cores

Well-dated benthic foraminifer oxygen isotopic records (d18O) from different water depths and locations within the Atlantic Ocean exhibit distinct patterns and significant differences in timing over the last deglaciation. This has two implications: on the one hand, it confirms that benthic d18O cannot be used as a global correlation tool with millennial-scale precision, but on the other hand, the combination of benthic isotopic records with independent dating provides a wealth of information on past circulation changes. Comparing new South Atlantic benthic isotopic data with published benthic isotopic records, we show that (1) circulation changes first affected benthic d18O in the 1000-2200 m range, with marked decreases in benthic d18O taking place at ~17.5 cal. kyr B.P. (ka) due to the southward propagation of brine waters generated in the Nordic Seas during Heinrich Stadial 1 (HS1) cold period; (2) the arrival of d18O-depleted deglacial meltwater took place later at deeper North Atlantic sites; (3) hydrographic changes recorded in North Atlantic cores below 3000 m during HS1 do not correspond to simple alternations between northern- and southern-sourced water but likely reflect instead the incursion of brine-generated deep water of northern as well as southern origin; and (4) South Atlantic waters at ~44°S and ~3800 m depth remained isolated from better-ventilated northern-sourced water masses until after the resumption of North Atlantic Deep Water (NADW) formation at the onset of the Bølling-Allerod, which led to the propagation of NADW into the South Atlantic.

Epibenthic foraminiferal d13C in the Recent deep Arctic Ocean

Low planktic and benthic d18O and d13C values in sediments from the Nordic seas of cold stadials of the last glaciation have been attributed to brines, formed similar to modern ones in the Arctic Ocean. To expand on the carbon isotopes of this hypothesis I investigated benthic d13C from the modern Arctic Ocean. I show that mean d13C values of live epibenthic foraminifera from the deep Arctic basins are higher than mean d13C values of upper slope epibenthic foraminifera. This agrees with mean high d13C values of dissolved inorganic carbon (DIC) in Arctic Bottom Water (ABW), which are higher than mean d13CDIC values from shallower water masses of mainly Atlantic origin. However, adjustments for oceanic 13C-Suess depletion raise subsurface and intermediate water d13CDIC values over ABW d13CDIC ones. Accordingly, during preindustrial Holocene times, the d13CDIC of ABW was as high or higher than today, but lower than the d13CDIC of younger subsurface and intermediate water. If brine-enriched water significantly ventilated ABW, brines should have had high d13CDIC values. Analogously, high-d13CDIC brines may have been formed in the Nordic seas during warm interstadials. During cold stadials, when most of the Arctic Ocean was perennially sea-ice covered, a cessation of high-d13CDIC brine rejection may have lowered d13CDIC values of ABW, and ultimately the d13CDIC in Nordic seas intermediate and deep water. So, in contrast to the idea of enhanced brine formation during cold stadials, the results of this investigation imply that a cessation of brine rejection would be more likely.

(Table 1) Age determination of composite sediment record LO09-14

A sediment core from Reykjanes Ridge has been studied at 10- to 50-year time resolution to document variability of Holocene surface water conditions in the western North Atlantic and to evaluate effects of Holocene ice-rafting episodes. Diatom assemblages are converted to quantitative sea surface temperatures (SST) using three different transfer functions. Spectral and scale-space methods are also applied on the records to explore variability at different timescales. Diatom assemblage and SST records clearly show that decaying remnants of the Laurentide ice sheet strongly influenced early Holocene climate in the western North Atlantic. This overrode the predominance of Milankovitch forcing, which played a key role in the development of Holocene climate in the eastern North Atlantic and Nordic Seas. Superimposed on general Holocene climate change is high-frequency SST variability on the order of 1°-3°C. The record also documents climatic oscillations with 600- to 1000-, ~1500-, and 2500-year periodicities, with a time-dependent dominance of different periodicities through the Holocene; a clear change in variability occurred about 5 ka BP. The SST record also provides evidence for Holocene cooling events (HCE) that, in some cases, correlate to documented southward intrusions of ice into the North Atlantic.

(Table 1) Age determination of Barents Sea sediment cores

he separate roles of oceanic heat advection and orbital forcing on influencing early Holocene temperature variability in the eastern Nordic Seas is investigated. The effect of changing orbital forcing on the ocean temperatures is tested using the 1DICE model, and the 1DICE results are compared with new and previously published temperature reconstructions from a transect of five cores located underneath the pathway of Atlantic water, from the Faroe-Shetland Channel in the south to the Barents Sea in the north. The stronger early Holocene summer insolation at high northern latitudes increased the summer mixed layer temperatures, however, ocean temperatures underneath the summer mixed layer did not increase significantly. The absolute maximum in summer mixed layer temperatures occurred between 9 and 6 ka BP, representing the Holocene Thermal Maximum in the eastern Nordic Seas. In contrast, maximum in northward oceanic heat transport through the Norwegian Atlantic Current occurred approximately 10 ka BP. The maximum in oceanic heat transport at 10 ka BP occurred due to a major reorganization of the Atlantic Ocean circulation, entailing strong and deep rejuvenation of the Atlantic Meridional Overturning Circulation, combined with changes in the North Atlantic gyre dynamic causing enhanced transport of heat and salt into the Nordic Seas.

(Table 2) Age determination of Irminger Basin sediment cores

High-resolution records of coarse lithic content and oxygen isotope have been obtained in a piston core from the Irminger Basin. The last glacial period is characterized by numerous periods of increased iceberg discharges originating partly from Iceland and corresponding to millennial-scale instabilities of the coastal ice sheets and ice shelves in the Nordic area. A comparison with midlatitude sediment cores shows that ice-rafted material corresponding to the Heinrich events was deposited synchronously from 40° to 60°N. There are thus two oscillating systems: every 5-10 kyr massive iceberg armadas are released from large continental ice caps, whereas more frequent instabilities of the coastal ice sheets in the high latitude regions occur every 1.2-3.8 kyr. At the time of the Heinrich events the synchroneity of the response from all the northern hemisphere ice sheets attests the existence of strong interactions between the two systems.