Automated Tools for Layer Recognition (BMPix) and Counting (PEAK)

We present tools for rapid and quantitative detection of sediment lamination. The BMPix tool extracts color and gray-scale curves from images at pixel resolution. The PEAK tool uses the gray-scale curve and performs, for the first time, fully automated counting of laminae based on three methods. The maximum count algorithm counts every bright peak of a couplet of two laminae (annual resolution) in a smoothed curve. The zero-crossing algorithm counts every positive and negative halfway-passage of the curve through a wide moving average, separating the record into bright and dark intervals (seasonal resolution). The same is true for the frequency truncation method, which uses Fourier transformation to decompose the curve into its frequency components before counting positive and negative passages. We applied the new methods successfully to tree rings, to well-dated and already manually counted marine varves from Saanich Inlet, and to marine laminae from the Antarctic continental margin. In combination with AMS14C dating, we found convincing evidence that laminations in Weddell Sea sites represent varves, deposited continuously over several millennia during the last glacial maximum. The new tools offer several advantages over previous methods. The counting procedures are based on a moving average generated from gray-scale curves instead of manual counting. Hence, results are highly objective and rely on reproducible mathematical criteria. Also, the PEAK tool measures the thickness of each year or season. Since all information required is displayed graphically, interactive optimization of the counting algorithms can be achieved quickly and conveniently.

Biomarker and XRD results for the last 3.4 Ma from IODP Site 306-U1313

Various types of abrupt/millennial-scale climate variability such as Dansgaard/Oeschger and Heinrich Events characterized the last glacial period. Over the last decade, a number of studies demonstrated that such millennial-scale climate variability was not limited to the last glacial but inherent to Quaternary climate. Here we review the occurrence and origin of millennial ice-rafting events in the North Atlantic during the late Pliocene and Pleistocene (last 3.4 Ma) with a special focus on North Atlantic Hudson Strait (HS) Heinrich(-like) Events. Besides a clear biomarker signature, we show that Heinrich Layers 5, 4, 2, and 1 in marine sediment cores from across the North Atlantic all bear the organic geochemical fingerprint of the Hudson area. Using this framework and combining previously published results, detailed investigations into the organic and inorganic chemistry of ice-rafted debris (IRD) found across the North Atlantic demonstrate that prior to MIS 16 (~ 650 ka) IRD in the North Atlantic did not originate from the Hudson area of northern Canada. The signature of this early IRD is distinctly different compared to that of HS Heinrich Layers. Rather ice-rafting events during the late Pliocene and early Pleistocene predominantly emanated from the calving of the Greenland and Fennoscandian ice sheets and possibly minor contributions from local ice streams from the North American and British ice sheets. Compared to North Atlantic HS Heinrich Events, these early Pleistocene IRD-events had a limited impact on surface water characteristics in the North Atlantic. North Atlantic HS Heinrich(-like) Events first occurred during MIS 16. At the same time, the dominant frequency in silicate-rich IRD accumulation shifted from the obliquity (41-ka) to a 100-ka frequency across the North Atlantic. Iceberg survivability or a change in iceberg trajectory likely did not control this change in IRD-regime. These results lend further support for the existing hypothesis that an increase in size (thickness) of the Laurentide ice sheet controls the occurrence of North Atlantic HS Heinrich Events, favoring an internal dynamic mechanism for their occurrence.

Clay mineral assemblages for the last 10 Ma off northeastern Australia

Clay mineral assemblages for the last 10 m.y. are described for Site 823, at 16°S in the Queensland Trough, to the northeast of Australia. Largely unaffected by diagenetic influences, these mostly express the evolution of northeastern Australian continental environments during the late Neogene: (1) beginning during the late Miocene at about 7.0 Ma is an increase of illite derived from rocky substrates at the expense of smectite from deeply weathered soils; this increase was the result of increasing aridity in the Australian interior and globally cooler temperatures, associated with increases in Antarctic glaciation; (2) concomitant and further increases of kaolinite fluxes to the Queensland Trough during the late Miocene-early Pliocene largely reflect an increase in rainfall in northeastern Australia; (3) increases in both soil- and rock-derived minerals probably intensified as a result of late Neogene uplift of the eastern highlands; (4) clay-mineral associations during the Pliocene and Pleistocene display minor variations only and probably resulted in part from differential settling and sea-level changes; (5) similar trends of clay-mineral variations occur at both ODP Site 823 and DSDP Site 588 (Lord Howe Rise). Less abundant kaolinite relative to illite at Site 588 nevertheless suggests a southward decrease of continental humidity and/or of the eastern highlands uplift; (6) influences of global climate and oceanic and atmospheric circulations on clay-mineral associations dominated during the late Miocene and were progressively replaced by influences of more regional environmental variations during the Pliocene and especially the Pleistocene.

Shell size variation of Neogloboquadrina pachyderma sin. in the Norwegian-Greenland Sea during the last 1.3 Myr

We present measurements of the maximum diameter of the planktonic foraminifer Neogloboquadrina pachyderma sin. from six sediment cores (Ocean Drilling Program sites 643, 644, 907, 909, 985 and 987) from the Norwegian-Greenland Sea. Our data show a distinct net increase in mean shell size of N. pachyderma sin. at all sites during the last 1.3 Ma, with largest shell sizes reached after 0.4 Ma. External factors such as glacial-interglacial variability and carbonate dissolution alone cannot account for the observed variation in mean shell size of N. pachyderma sin. We consider the observed shell size increase to mirror an evolutionary trend towards better adaptation of N. pachyderma sin. to the cold water environment after 1.1-1.0 Ma. Probably, the Mid Pleistocene climate shift and the associated change of amplitude and frequency of glacial-interglacial fluctuations have triggered the evolution of this planktonic foraminifer. Oxygen and carbon stable isotope analyses of different shell size classes indicate that the observed shell size increase could not be explained by the functional concept that larger shells promote increasing sinking velocities during gametogenesis. For paleoceanographic reconstructions, the evolutionary adaptation of Neogloboquadrina pachyderma sin. to the cold water habitat has significant implications. Carbonate sedimentation in highest latitudes is highly dependent on the presence of this species. In the Norwegian-Greenland Sea, carbonate-poor intervals before 1.1 Ma are, therefore, not necessarily related to severe glacial conditions. They are probably attributed to the absence of this not yet polar-adapted species. Further, transfer function and modern analog techniques used for the reconstruction of surface water conditions in high latitudes could, therefore, contain a large range of errors if they were applied to samples older than 1.1-1.0 Myrs.

DATED-1: compilation of dates and time-slice reconstruction of the build-up and retreat of the last Eurasian (British-Irish, Scandinavian, Svalbard-Barents-Kara Seas) Ice Sheets 40-10 ka

DATED-1 comprises a compilation of dates related to the build-up and retreat of the Eurasian (British-Irish, Scandinavian, Svalbard-Barents-Kara Seas) Ice Sheets, and time-slice maps of the Eurasian Ice sheet margins. Dates are sourced from the published literature. Ice margins are based on published geological and chronological data and include uncertainty bounds (maximum, minimum) as well as what we consider to be the most-credible (mc) based on the available evidence. DATED-1 has a census date of 1 January 2013. Full description and caveats for use are given in: Hughes, A.L.C., Gyllencreutz, R., Lohne, Ø.S., Mangerud, J., Svendsen, J.I. (2015) The last Eurasian Ice Sheets - a chronological database and time-slice reconstruction, DATED-1.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

We use the fully coupled atmosphere-ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 yr). By using a model that is able to explicitly simulate the sensor (d18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean d18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice-sheet. We observe a general pattern toward more enriched values for continental calcite d18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data-model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data-model comparison for calcite d18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a trong mean annual cooling between the LGM and present (>6°C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data-model comparison also reveals that large positive calcite d18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in d18Osw between the LGM and present are not spatially homogenous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Benthic foraminifera stable isotope data from the South Pacific sediment cores PS75/059-2 and PS75/056-1

Investigating the inter-basin deep water exchange between the Pacific and Atlantic Oceans over glacial-interglacial climate cycles is important for understanding circum-Antarctic Southern Ocean circulation changes and their impact on the global Meridional Overturning Circulation. We use benthic foraminiferal d13C records from the southern East Pacific Rise to characterize the d13C composition of Circumpolar Deep Water in the South Pacific, prior to its transit through the Drake Passage into the South Atlantic. A comparison with published South Atlantic deep water records from the northern Cape Basin suggests a continuous water mass exchange throughout the past 500 ka. Almost identical glacial-interglacial d13C variations imply a common deep water evolution in both basins suggesting persistent Circumpolar Deep Water exchange and homogenization. By contrast, deeper abyssal waters occupying the more southern Cape Basin and the southernmost South Atlantic have lower d13C values during most, but not all, stadial periods. We conclude that these values represent the influence of a more southern water mass, perhaps AABW. During many interglacials and some glacial periods, the gradient between Circumpolar Deep Water and the deeper southern Cape Basin bottom water disappears suggesting either no presence of AABW or indistinguishable d13C values of both water masses.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.