Clay mineralogy and multi-element chemistry of surface sediments on the Aiberian-Arctic Shelf

Clay mineral and bulk chemical (Si, Al, K, Mg, Sr, La, Ce, Nd) analyses of terrigenous surface sediments on the Siberian-Arctic shelf indicate that there are five regions with distinct, or endmember, sedimentary compositions. The formation of these geochemical endmembers is controlled by sediment provenance and grain size sorting. (1) The shale endmember (Al, K and REE rich sediment) is eroded from fine-grained marine sedimentary rocks of the Verkhoyansk Mountains and Kolyma-Omolon superterrain, and discharged to the shelf by the Lena, Yana, Indigirka and Kolyma Rivers. (2) The basalt endmember (Mg rich) originates from NE Siberia's Okhotsk-Chukotsk volcanic belt and Bering Strait inflow, and is prevalent in Chukchi Sea Sediments. Concentrations of the volcanically derived clay mineral smectite are elevated in Chukchi fine-fraction sediments, corroborating the conclusion that Chukchi sediments are volcanic in origin. (3) The mature sandstone endmember (Si rich) is found proximal to Wrangel Island and sections of the Chukchi Sea's Siberian coast and is derived from the sedimentary Chukotka terrain that comprises these landmasses. (4) The immature sandstone endmember (Sr rich) is abundant in the New Siberian Island region and reflects inputs from sedimentary rocks that comprise the islands. (5) The immature sandstone endmember is also prevalent in the western Laptev Sea, where it is eroded from sedimentary deposits blanketing the Siberian platform that are compositionally similar to those on the New Siberian Islands. Western Laptev can be distinguished from New Siberian Island region sediments by their comparatively elevated smectite concentrations and the presence of the basalt endmember, which indicate Siberian platform flood basalts are also a source of western Laptev sediments. In certain locations grain size sorting noticeably affects shelf sediment chemistry. (1) Erosion of fines by currents and sediment ice rafting contributes to the formation of the coarse-grained sandstone endmembers. (2) Bathymetrically controlled grain size sorting, in which fines preferentially accumulate offshore in deeper, less energetic water, helps distribute the fine-grained shale and basalt endmembers. An important implication of these results is that the observed sedimentary geochemical endmembers provide new markers of sediment provenance, which can be used to track sediment transport, ice-rafted debris dispersal or the movement of particle-reactive contaminants.

Deuterium analysis of ice core GISP2

The Greenland Ice Sheet Project 2 (GISP2) core can enhance our understanding of the relationship between parameters measured in the ice in central Greenland and variability in the ocean, atmosphere, and cryosphere of the North Atlantic Ocean and adjacent land masses. Seasonal (summer, winter) to annual responses of dD and deuterium excess isotopic signals in the GISP2 core to the seesaw in winter temperatures between West Greenland and northern Europe from A.D. 1840 to 1970 are investigated. This seesaw represents extreme modes of the North Atlantic Oscillation, which also influences sea surface temperatures (SSTs), atmospheric pressures, geostrophic wind strength, and sea ice extents beyond the winter season. Temperature excursions inferred from the dD record during seesaw/extreme NAO mode years move in the same direction as the West Greenland side of the seesaw. Symmetry with the West Greenland side of the seesaw suggests a possible mechanism for damping in the ice core record of the lowest decadal temperatures experienced in Europe from A.D. 1500 to 1700. Seasonal and annual deuterium excess excursions during seesaw years show negative correlation with dD. This suggests an isotopic response to a SST/ land temperature seesaw. The isotopic record from GISP2 may therefore give information on both ice sheet and sea surface temperature variability. Cross-plots of dD and d show a tendency for data to be grouped according to the prevailing mode of the seesaw, but do not provide unambiguous identification of individual seesaw years. A combination of ice core and tree ring data sets may allow more confident identification of GA and GB (extreme NAO mode) years prior to 1840.

(Appendix A) Geochemical composition of snow and eolian mineral dust from Berkner Island ice sheet, Antarctica

The Sr and Nd isotopic composition of dust extracted from recent snow layers at the top of Berkner Island ice sheet (located within the Filchner-Ronne Ice Shelf at the southern end of the Weddell Sea) enables us, for the first time, to document dust provenance in Antarctica outside the East Antarctic Plateau (EAP) where all previous studies based on isotopic fingerprinting were carried out. Berkner dust displays an overall crust-like isotopic signature, characterized by more radiogenic 87Sr/86Sr and much less radiogenic 143Nd/144Nd compared to dust deposited on the EAP during glacial periods. Differences with EAP interglacial dust are not as marked but still significant, indicating that present-day Berkner dust provenance is distinct, at least to some extent, from that of the dust reaching the EAP. The fourteen snow-pit sub-seasonal samples that were obtained span a two-year period (2002-2003) and their dust Sr and Nd isotopic composition reveals that multiple sources are at play over a yearly time period. Southern South America, Patagonia in particular, likely accounts for part of the observed spring/summer dust deposition maxima, when isotopic composition is shifted towards 'younger' isotopic signatures. In the spring, possible additional inputs from Australian sources would also be supported by the data. Most of the year, however, the measured isotopic signatures would be best explained by a sustained background supply from putative local sources in East Antarctica, which carry old-crust-like isotopic fingerprints. Whether the restricted East Antarctic ice-free areas produce sufficient eolian material has yet to be substantiated however. The fact that large (> 5 µm) particles represent a significant fraction of the samples throughout the entire time-series supports scenarios that involve contributions from proximal sources, either in Patagonia and/or Antarctica (possibly including snow-free areas in the Antarctic Peninsula and other areas as well). This also indicates that additional dust transport, which does not reach the EAP, must occur at low-tropospheric levels to this coastal sector of Antarctica.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Under ice shelf photographs from Drescher Inlet (Riiser Larsen Ice Shelf) taken by seal mounted cameras during expedition DRE2003

While modern sampling techniques, such as autonomous underwater vehicles, are increasing our knowledge of the fauna beneath Antarctic sea ice of only a few meters in depth, greater sampling difficulties mean that little is known about the marine life underneath Antarctic ice shelves over 100 m thick. In this study, we present underwater images showing the underside of an Antarctic ice shelf covered by aggregated invertebrate communities, most likely cnidarians and isopods. These images, taken at an average depth of 145 m, were obtained with a digital still camera system attached to Weddell seals Leptonychotes weddellii foraging just beneath the ice shelf. Our observations indicate that, similar to the sea floor, ice shelves serve as an important habitat for a remarkable amount of marine invertebrate fauna in Antarctica.

Iron (Fe) concentrations and fluxes in TALDICE ice core 0.6 to 314 ky BP

Atmospheric fluxes of iron (Fe) over the past 200 kyr are reported for the coastal Antarctic Talos Dome ice core, based on acid leachable Fe concentrations. Fluxes of Fe to Talos Dome were consistently greater than those at Dome C, with the greatest difference observed during interglacial climates. We observe different Fe flux trends at Dome C and Talos Dome during the deglaciation and early Holocene, attributed to a combination of deglacial activation of dust sources local to Talos Dome and the reorganisation of atmospheric transport pathways with the retreat of the Ross Sea ice shelf. This supports similar findings based on dust particle sizes and fluxes and Rare Earth Element fluxes. We show that Ca and Fe should not be used as quantitative proxies for mineral dust, as they all demonstrate different deglacial trends at Talos Dome and Dome C. Considering that a 20 ppmv decrease in atmospheric CO2 at the coldest part of the last glacial maximum occurs contemporaneously with the period of greatest Fe and dust flux to Antarctica, we confirm that the maximum contribution of aeolian dust deposition to Southern Ocean sequestration of atmospheric CO2 is approximately 20 ppmv.

Glaciological characteristics and carbon dioxide record of ice cores from Antarctica

We report a decadally resolved record of atmospheric CO2 concentration for the last 1000 years, obtained from the West Antarctic Ice Sheet (WAIS) Divide shallow ice core. The most prominent feature of the pre-industrial period is a rapid ~7 ppm decrease of CO2 in a span of ~20-50 years at ~1600 A.D. This observation confirms the timing of an abrupt atmospheric CO2 decrease of ~10 ppm observed for that time period in the Law Dome ice core CO2 records, but the true magnitude of the decrease remains unclear. Atmospheric CO2 variations over the time period 1000-1800 A.D. are statistically correlated with northern hemispheric climate and tropical Indo-Pacific sea surface temperature. However, the exact relationship between CO2 and climate remains elusive due to regional climate variations and/or uneven geographical data density of paleoclimate records. We observe small differences of 0 ~2% (0 ~ 6 ppm) among the high-precision CO2 records from the Law Dome, EPICA Dronning Maud Land and WAIS Divide Antarctic ice cores. However, those records share common trends of CO2 change on centennial to multicentennial time scales, and clearly show that atmospheric CO2 has been increasing above preindustrial levels since ~1850 A.D.

(Table 1) Deuterium excess record during different time periods in Lomonosovfonna ice core

In this paper we present a deuterium excess (d) record from an ice core drilled on a small ice cap in Svalbard in 1997. The core site is located at Lomonosovfonna at 1255 m asl, and the analyzed time series spans the period 1400-1990 A.D. The record shows pronounced multidecadal to centennial-scale variations coherent with sea surface temperature changes registered in the subtropical to southern middle-latitude North Atlantic during the instrumental period. We interpret the negative trend in the deuterium excess during the 1400s and 1500s as an indication of cooling in the North Atlantic associated with the onset of the Little Ice Age. Consistently positive anomalies of d after 1900, peaking at about 1950, correspond with well-documented contemporary warming. Yet the maximum values of deuterium excess during 1900-1990 are not as high as in the early part of the record (pre-1550). This suggests that the sea surface temperatures during this earlier period of time in the North Atlantic to the south of approximately 45°N were at least comparable with those registered in the 20th century before the end of the 1980s. We examine the potential for a cold bias to exist in the deuterium excess record due to increased evaporation from the local colder sources of moisture having isotopically cold signature. It is argued that despite a recent oceanic warming, the contribution from this local moisture to the Lomonosovfonna precipitation budget is still insufficient to interfere with the isotopic signal from the primary moisture region in the midlatitude North Atlantic.

Sedimentology on 19 cores from the South China Sea

Based on the study of 10 sediment cores and 40 core-top samples from the South China Sea (SCS) we obtained proxy records of past changes in East Asian monsoon climate on millennial to bidecadal time scales over the last 220,000 years. Climate proxies such as global sea level, estimates of paleotemperature, salinity, and nutrients in surface water, ventilation of deep water, paleowind strength, freshwater lids, fluvial and/or eolian sediment supply, and sediment winnowing on the sea floor were derived from planktonic and benthic stable-isotope records, the distribution of siliciclastic grain sizes, planktonic foraminifera species, and the UK37 biomarker index. Four cores were AMS-14C-dated. Two different regimes of monsoon circulation dominated the SCS over the last two glacial cycles, being linked to the minima and maxima of Northern Hemisphere solar insolation. (1) Glacial stages led to a stable estuarine circulation and a strong O2-minimum layer via a closure of the Borneo sea strait. Strong northeast monsoon and cool surface water occurred during winter, in part fed by an inflow from the north tip of Luzon. In contrast, summer temperatures were as high as during interglacials, hence the seasonality was strong. Low wetness in subtropical South China was opposed to large river input from the emerged Sunda shelf, serving as glacial refuge for tropical forest. (2) Interglacials were marked by a strong inflow of warm water via the Borneo sea strait, intense upwelling southeast of Vietnam and continental wetness in China during summer, weaker northeast monsoon and high sea-surface temperatures during winter, i.e. low seasonality. On top of the long-term variations we found millennial- to centennial-scale cold and dry, warm and humid spells during the Holocene, glacial Terminations I and II, and Stage 3. The spells were coeval with published variations in the Indian monsoon and probably, with the cold Heinrich and warm Dansgaard-Oeschger events recorded in Greenland ice cores, thus suggesting global climatic teleconnections. Holocene oscillations in the runoff from South China centered around periodicities of 775 years, ascribed to subharmonics of the 1500-year cycle in oceanic thermohaline circulation. 102/84-year cycles are tentatively assigned to the Gleissberg period of solar activity. Phase relationships among various monsoon proxies near the onset of Termination IA suggest that summer-monsoon rains and fluvial runoff from South China had already intensified right after the last glacial maximum (LGM) insolation minimum, coeval with the start of Antarctic ice melt, prior to the d18O signals of global sea-level rise. Vice versa, the strength of winter-monsoon winds decreased in short centennial steps only 3000-4000 years later, along with the melt of glacial ice sheets in the Northern Hemisphere.

Chlorophyll a concentration and zooplankton abundance and biomass below melting sea ice in the Amundsen Gulf

Early summer in the Arctic with extensive ice melt and break-up represents a dramatic change for sympagic-pelagic fauna below seasonal sea ice. As part of the International Polar Year-Circumpolar Flaw Lead system study (IPY-CFL), this investigation quantified zooplankton in the meltwater layer below landfast ice and remaining ice fauna below melting ice during June (2008) in Franklin Bay and Darnley Bay, Amundsen Gulf, Canada. The ice was in a state of advanced melt, with fully developed melt ponds. Intense melting resulted in a 0.3- to 0.5-m-thick meltwater layer below the ice, with a strong halocline to the Arctic water below. Zooplankton under the ice, in and below the meltwater layer, was sampled by SCUBA divers. Dense concentrations (max. 1,400 ind./m**3) of Calanus glacialis were associated with the meltwater layer, with dominant copepodid stages CIV and CV and high abundance of nauplii. Less abundant species included Pseudocalanus spp., Oithona similis and C. hyperboreus. The copepods were likely feeding on phytoplankton (0.5-2.3 mg Chl-a/m**3) in the meltwater layer. Ice amphipods were present at low abundance (<10 ind./m**2) and wet biomass (<0.2 g/m**2). Onisimus glacialis and Apherusa glacialis made up 64 and 51% of the total ice faunal abundance in Darnley Bay and Franklin Bay, respectively. During early summer, the autochthonous ice fauna becomes gradually replaced by allochthonous zooplankton, with an abundance boom near the meltwater layer. The ice amphipod bust occurs during late stages of melting and break-up, when their sympagic habitat is diminished then lost.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.