Profiles of diatoms, pollen, XRD, XRF and organic composition of two sediment cores (PG2022 and PG2023) from Lake Kyutyunda in Yakutia, NE Siberia, Russia

Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between ~13,500 and ~8,900 cal. years BP and possibly during the 8,200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8°C during the Holocene Thermal Maximum (HTM) between ~8,900 and ~4,500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3,000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice-sheet vanished 7,000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions.

Granulometry, biogenic silica content, magnetic susceptibility and bulk densities of ODP Site 188-1165

The Prydz Bay area is a key region for studying and understanding the history of the eastern Antarctic Continental Ice Sheet (O'Brien, Cooper, Richter, et al., 2001, doi:10.2973/odp.proc.ir.188.2001). Ocean Drilling Program (ODP) Site 1165 is situated in a water depth of 3357 m on the continental rise offshore from Prydz Bay and lies in front of the outlet for the Lambert Glacier-Amery Ice Shelf system that today drains 22% of East Antarctica. The site was drilled into mixed pelagic and hemipelagic sediments from the southwestern side of the Wild Drift. The drift is an elongate sediment body formed by the interaction of sediment supplied from continental shelf and slope with westward-flowing bottom currents. The sedimentary sequence is characterized by alternations between a generally gray to dark gray facies and a green to greenish gray facies. The greenish facies are structureless diatom-bearing clays with common bioturbation and larger amounts (>15%-20%) of biogenic silica, dispersed clasts, and lonestones than the dark gray facies, which are mostly less bioturbated clay with some silt laminations (Shipboard Scientific Party, 2001, doi:10.2973/odp.proc.ir.188.103.2001). High-quality advanced piston corer and extended core barrel cores containing nearly complete sections of middle Miocene to early Pliocene age allow a detailed characterization of sedimentary cycles and can provide indications for ice advances of the Lambert Glacier system into Prydz Bay, for the extent of sea ice, and for changes in oceanic circulation. The purpose of this work is to provide a data set of coarse-fraction mass percentage (>63, >125, and >250 µm) and biogenic silica content measured on sediments of late Miocene to early Pliocene age drilled at Site 1165. Additionally, high-resolution records of magnetic susceptibility (MS) and gamma ray attenuation (GRA) bulk density are presented. These shipboard data sets were edited postcruise. Furthermore, I provide a high-resolution dry bulk density record that is derived from GRA bulk density and can be used for the calculation of mass accumulation rates. These sedimentological and physical parameters will be used in future work to understand the depositional pattern of alternating biogenic and terrigenous sediments that was observed at Site 1165 (Shipboard Scientific Party, 2001, doi:10.2973/odp.proc.ir.188.103.2001).

Mean air temperatures of Greenland

Melt rate and surface temperature on the Greenland ice sheet are parameterized in terms of snow accumulation, mean annual air temperatur and mean July air temperature. Melt rates are calculated using positive degree-days, and firn warming (i.e. the positive deviation of the temperature at 10-15 m depth from the mean annual air temperature) is estimated from the calculated amount of refrozen melt water in the firn. A comparison between observed and calculated melt rates shows that the parameterization provides a reasonable estimate of the present ablation rates in West Greenland between 61°N and 76°N. The average equilibrium line elevation is estimated to be about 1150 m and 1000 m for West and East Greenland respectively, which is several hundred meter lower than previous estimates. However, the total annual ablation from the ice sheet is found to be about 280 km**3 of water per year which agrees well with most other estimates. The melt-rate model predicts significant melting and consequently significant firn warming even at the highest elevations of the South Greenland ice sheet, whereas a large region of central Greenland north of 70° N experiences little or no summer melting. This agrees with the distribution of the dry snow facics as given by BENSON (1962).

Geochemistry of FeMn oxyhydroxides of IODP Site 303-U1302 and 303-1303

Regional weathering intensity must have changed dramatically at high latitudes during the Quaternary as a consequence of repeated continental glaciation. Investigation of these glacial/interglacial changes at high temporal resolution is possible with the recent development of Pb isotopes in FeMn oxyhydroxide phases as a proxy for region-specific weathering intensity, where increases in the radiogenic component are thought to correspond to increased continental weathering fluxes. Here we present a Pb isotope record sourced from the FeMn oxyhydroxide fraction in marine sediments from IODP Sites U1302/3 on Orphan Knoll (~3500 mbsl, NW Atlantic), spanning the last 37 ka. Located at the eastern edge of the Laurentide Ice Sheet (LIS), Site U1302/3 is well-placed to monitor changes in weathering intensity associated with LIS glacial history. Overall, the data show a close correspondence to local surface water d18O, with least radiogenic values during times of heavy d18O (glacial maximum) and most radiogenic values during times of light d18O (Holocene). This supports the prediction that weathering intensity in glaciated regions of the North Atlantic correlates with the exposure age of glacial debris. Superimposed on these background trends are extreme radiogenic excursions (e.g. variation in 206Pb/204Pb from ~19.2-21.0) contemporaneous with Heinrich events and the Younger Dryas. These data are substantially more radiogenic than existing records from the NW Atlantic, and most likely represent episodes of exceptionally high inputs of pre-formed FeMn oxyhydroxides during drainage of the LIS. Due to its extreme isotope composition, at least in the NW Atlantic region, Pb would appear to be a good proxy for the fluxes of weathered continental material and perhaps, by inference, nutrients to the surface ocean

Elevation Models of Antarctica and Greenland derived from CryoSat-2 in the period 2011 to 2013

This study focuses on the present-day surface elevation of the Greenland and Antarctic ice sheets. Based on 3 years of CryoSat-2 data acquisition we derived new elevation models (DEMs) as well as elevation change maps and volume change estimates for both ice sheets. Here we present the new DEMs and their corresponding error maps. The accuracy of the derived DEMs for Greenland and Antarctica is similar to those of previous DEMs obtained by satellite-based laser and radar altimeters. Comparisons with ICESat data show that 80% of the CryoSat-2 DEMs have an uncertainty of less than 3 m ± 15 m. The surface elevation change rates between January 2011 and January 2014 are presented for both ice sheets. We compared our results to elevation change rates obtained from ICESat data covering the time period from 2003 to 2009. The comparison reveals that in West Antarctica the volume loss has increased by a factor of 3. It also shows an anomalous thickening in Dronning Maud Land, East Antarctica which represents a known large-scale accumulation event. This anomaly partly compensates for the observed increased volume loss of the Antarctic Peninsula and West Antarctica. For Greenland we find a volume loss increased by a factor of 2.5 compared to the ICESat period with large negative elevation changes concentrated at the west and southeast coasts. The combined volume change of Greenland and Antarctica for the observation period is estimated to be -503 ± 107 km**3/yr. Greenland contributes nearly 75% to the total volume change with -375 ± 24 km**3/yr.

Radiocarbon ages and clay minerals in sediment cores obtained during a 2009 Jan Mayen cruise to the Norwegian Sea

The origin of two acoustic sediment units has been studied based on lithological facies, chronology and benthic stable isotope values as well as on foraminifera and clay mineral assemblages in six marine sediment cores from Kveithola, a small trough west of Spitsbergenbanken on the western Barents Sea margin. We have identified four time slices with characteristic sedimentary environments. Before c. 14.2 cal. ka, rhythmically laminated muds indicate extensive sea ice cover in the area. From c. 13.9 to 14.2 cal. ka, muds rich in ice-rafted debris were deposited during the disintegration of grounded ice on Spitsbergenbanken. From c. 10.3 to 13.1 cal. ka, sediments with heterogeneous lithologies suggest a shifting influence of suspension settling and iceberg rafting, probably derived from a decaying Barents Sea Ice Sheet in the inner-fjord and land areas to the north of Kveithola. Holocene deposition was episodic and characterized by the deposition of calcareous sands and shell debris, indicative of strong bottom currents. We speculate that a marked erosional boundary at c. 8.2 cal. ka may have been caused by the Storegga tsunami. Whilst deposition was sparse during the Holocene, Kveithola acted as a sediment trap during the preceding deglaciation. Investigation of the deglacial sediments provides unprecedented details on the dynamics and timing of glacial retreat from Spitsbergenbanken.

Proglacial river dataset from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008 - 2013, version 2.0

The dataset contains measurements of river stage and discharge for one sites along the Akuliarusiarsuup Kuua River's northern tributary, with 30 minute temporal resolution between June 2008 and August 2013 This river is a tributary to the Watson River discharging into Kangerlussuaq Fjord by the town of Kangerlussuaq, Southwest Greenland. Additional data of water temperature, air pressure are also provided. Compared to version 1.0 of the dataset, this dataset used a total of 36 in situ discharge observations collected between 2008 and 2012 to construct the rating curve. Furthermore, data of Station AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 have been removed from version 2.0 because these values were likely caused by backflow when a jokulhlaup from a large glacier dammed lake caused increased water levels in the downstreams lake. Thus, data measured at AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 are not representative for the AK-004 catchment.

Age model and Plio-Pleistocene benthic stable oxygen isotope ratios of ODP Site 184-1143 in the Sourh China Sea

Based on benthic foraminiferal delta18O from ODP Site 1143, a 5-Myr astronomical timescale for the West Pacific Plio-Pleistocene was established using an automatic orbital tuning method. The tuned Brunhes/Matuyama paleomagnetic polarity reversal age agrees well with the previously published age of 0.78 Ma. The tuned ages for several planktonic foraminifer bio-events also agree well with published dates, and new ages for some other bio-events in the South China Sea were also estimated. The benthic delta18O from Site 1143 is highly coherent with the Earth's orbit (ETP) both at the obliquity and precession bands for the last 5 Myr, and at the eccentricity band for the last 2 Myr. In general, the 41-kyr cycle was dominant through the Plio-Pleistocene although the 23-kyr cycle was also very strong. The 100-kyr cycle became dominant only during the last 1 Myr. A comparison of the benthic delta18O between the Atlantic (ODP 659) and the East and West Pacific (846 and 1143) reveals that the Atlantic-Pacific benthic oxygen isotope difference ratio (Delta delta18OAtl-Pac) displays an increasing trend in three time intervals: 3.6-2.7 Ma, 2.7-2.1 Ma and 1.5-0.25 Ma. Each of the intervals begins with a rapid negative shift in Delta delta18OAtl-Pac, followed by a long period with an increasing trend, corresponding to the growth of the Northern Hemisphere ice sheet. This means that all three intervals of ice sheet growth in the Northern Hemisphere were accompanied at the beginning by a rapid relative warming of deep water in the Atlantic as compared to that of the Pacific, followed by its gradual relative cooling. This general trend, superimposed on the frequent fluctuations with glacial cycles, should yield insights into the processes leading to the boreal glaciation. Cross-spectral analyses of the Delta delta18OAtl-Pac with the Earth's orbit suggests that after the initiation of Northern Hemisphere glaciation at about 2.5 Ma, obliquity rather than precession had become the dominant force controlling the vertical structure or thermohaline circulation in the paleo-ocean.

(Table S1) Flux gate discharge and centerline discharge of Helheim, Kangerdlugssuaq and Jakobshavn glacier, Greenland

Acceleration of Greenland's three largest outlet glaciers, Helheim, Kangerdlugssuaq and Jakobshavn Isbræ, accounted for a substantial portion of the ice sheet's mass loss over the past decade. Rapid changes in their discharge, however, make their cumulative mass-change uncertain. We derive monthly mass balance rates and cumulative balance from discharge and surface mass balance (SMB) rates for these glaciers from 2000 through 2010. Despite the dramatic changes observed at Helheim, the glacier gained mass over the period, due primarily to the short duration of acceleration and a likely longer-term positive balance. In contrast, Jakobshavn Isbræ lost an equivalent of over 11 times the average annual SMB and loss continues to accelerate. Kangerdlugssuaq lost over 7 times its annual average SMB, but loss has returned to the 2000 rate. These differences point to contrasts in the long-term evolution of these glaciers and the danger in basing predictions on extrapolations of recent changes.