A comparison between multi-proxy and historical data of drift-ice conditions on the East Greenland shelf

The reduction in sea ice along the SE Greenland coast during the last century has severely impacted ice-rafting to this area. In order to reconstruct ice-rafting and oceanographic conditions in the area of Denmark Strait during the last ~150 years, we conducted a multiproxy study on three short (20 cm) sediment cores from outer Kangerdlugssuaq Trough (~300 m water depth). The proxy-based data obtained have been compared with historical and instrumental data to gain a better understanding of the ice sheet-ocean interactions in the area. A robust chronology has been developed based on 210Pb and 137Cs measurements on core PO175GKC#9 (~66.2°N, 32°W) and expanded to the two adjacent cores based on correlations between calcite weight percent records. Our proxy records include sea-ice and phytoplankton biomarkers, and a variety of mineralogical determinations based on the <2 mm sediment fraction, including identification with quantitative x-ray diffraction, ice-rafted debris counts on the 63-150 µm sand fraction, and source identifications based on the composition of Fe oxides in the 45-250 µm fraction. A multivariate statistical analysis indicated significant correlations between our proxy records and historical data, especially with the mean annual temperature data from Stykkishólmur (Iceland) and the storis index (historical observations of sea-ice export via the East Greenland Current). In particular, the biological proxies (calcite weight percent, IP25, and total organic carbon %) showed significant linkage with the storis index. Our records show two distinct intervals in the recent history of the SE Greenland coast. The first of these (ad 1850-1910) shows predominantly perennial sea-ice conditions in the area, while the second (ad 1910-1990) shows more seasonally open water conditions.

Mid-Wisconsin to Holocene permafrost and landscape dynamics based on a drained lake basin core from the northern Seward Peninsula, Northwest Alaska

Permafrost-related processes drive regional landscape dynamics in the Arctic terrestrial system. A better understanding of past periods indicative of permafrost degradation and aggradation is important for predicting the future response of Arctic landscapes to climate change. Here, we used a multi-proxy approach to analyze a ~4 m long sediment core from a drained thermokarst lake basin on the northern Seward Peninsula in western Arctic Alaska (USA). Sedimentological, biogeochemistical, geochronological, micropaleontological (ostracoda, testate amoeba) and tephra analyses were used to determine the long-term environmental Early-Wisconsin to Holocene history preserved in our core for Central Beringia. Yedoma accumulation dominated throughout the Early to Late-Wisconsin but was interrupted by wetland formation from 44.5 to 41.5 ka BP. The latter was terminated by deposition of 1 m of volcanic tephra, most likely originating from the South Killeak Maar eruption at about 42 ka BP. Yedoma deposition continued until 22.5 ka BP and was followed by a depositional hiatus in the sediment core between 22.5 and 0.23 ka BP. We interpret this hiatus as due to intense thermokarst activity in the areas surrounding the site, which served as a sediment source during the Late-Wisconsin to Holocene climate transition. The lake forming the modern basin on the upland initiated around 0.23 ka BP, which drained catastrophically in spring 2005. The present study emphasizes that Arctic lake systems and periglacial landscapes are highly dynamic and permafrost formation as well as degradation in Central Beringia was controlled by regional to global climate patterns and as well as by local disturbances.

Multi-sensor global compilation of mid-Holocene sea surface temperatures, based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst

We present and examine a multi-sensor global compilation of mid-Holocene (MH) sea surface temperatures (SST), based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst census counts. We assess the uncertainties originating from using different methodologies and evaluate the potential of MH SST reconstructions as a benchmark for climate-model simulations. The comparison between different analytical approaches (time frame, baseline climate) shows the choice of time window for the MH has a negligible effect on the reconstructed SST pattern, but the choice of baseline climate affects both the magnitude and spatial pattern of the reconstructed SSTs. Comparison of the SST reconstructions made using different sensors shows significant discrepancies at a regional scale, with uncertainties often exceeding the reconstructed SST anomaly. Apparent patterns in SST may largely be a reflection of the use of different sensors in different regions. Overall, the uncertainties associated with the SST reconstructions are generally larger than the MH anomalies. Thus, the SST data currently available cannot serve as a target for benchmarking model simulations.

Continuous 25-years aerosol records at coastal Antarctica: Part I. Variability and climate implications of ionic compounds

The aerosol climatology at the coastal Antarctic Neumayer Station (NM) was investigated based on continuous, 25-yr long observations of biogenic sulphur components (methanesulfonate and non-sea salt sulphate), sea salt and nitrate. Although significant long-term trends could only be detected for nitrate (-3.6 ± 2.5% per year between 1983 and 1993 and +4.0 ± 3.2% per year from 1993-2007), non-harmonic periodicities between 2 and 5 yr were typical for all species. Dedicated time series analyses revealed that relations to sea ice extent and various circulation indices are weak at best or not significant. In particular, no consistent link between sea ice extent and sea salt loadings was evident suggesting only a rather local relevance of the NM sea salt record. Nevertheless, a higher Southern Annular Mode index tended to entail a lower biogenic sulphur signal. In examining the spatial uniformity of the NM findings we contrasted them to respective 17 yr records from the coastal Dumont d'Urville Station. We found similar long-term trends for nitrate, indicating an Antarctic-wide but not identifiable atmospheric signal, although any significant impact of solar activity or pollution could be ruled out. No inter-site variability on the multiannual scale was evident for the other ionic compounds.

Characterization of late Campanian and Maastrichtian planktonic foraminiferal depth habitats and vital activities based on stable isotopes

Depth habitats of 56 late Cretaceous planktonic foraminiferal species from cool and warm climate modes were determined based on stable isotope analyses of deep-sea samples from the equatorial Pacific DSDP Sites 577A and 463, and South Atlantic DSDP Site 525A. The following conclusions can be reached: Planoglobulina multicamerata (De Klasz) and Heterohelix rajagopalani (Govindan) occupied the deepest plankton habitats, followed by Abathomphalus mayaroensis (Bolli), Globotruncanella havanensis (Voorwijk), Gublerina cuvillieri Kikoine, and Laeviheterohelix glabrans (Cushman) also at subthermocline depth. Most keeled globotruncanids, and possibly Globigerinelliodes and Racemiguembelina species, lived at or within the thermocline layer. Heterohelix globulosa (Ehrenberg) and Rugoglobigerina, Pseudotextularia and Planoglobulina occupied the subsurface depth of the mixed layer, and Pseudoguembelina species inhabited the surface mixed layer. However, depth ranking of some species varied depending on warm or cool climate modes, and late Campanian or Maastrichtian age. For example, most keeled globotruncanids occupied similar shallow subsurface habitats as Rugoglobigerina during the warm late Campanian, but occupied the deeper thermocline layer during cool climatic intervals. Two distinct types of "vital effect" mechanisms reflecting photosymbiosis and respiration effects can be recognized by the exceptional delta13C signals of some species. (1) Photosymbiosis is implied by the repetitive pattern of relatively enriched delta13C values of Racemiguembelina (strongest), Planoglobulina, Rosita and Rugoglobigerina species, Pseudoguembelina excolata (weakest). (2) Enriched respiration 12C products are recognized in A. mayaroensis, Gublerina acuta De Klasz, and Heterohelix planata (Cushman). Isotopic trends between samples suggest that photosymbiotic activities varied between localities or during different climate modes, and may have ceased under certain environmental conditions. The appearance of most photosymbiotic species in the late Maastrichtian suggests oligotrophic conditions associated with increased water-mass stratification.

Gotland Deep (Baltic Sea) independent age model based on Palaeomagnetic Secular Variation (PSV) and lead (Pb) deposition profiles

Dating of sediment cores from the Baltic Sea has proven to be difficult due to uncertainties surrounding the 14C reservoir age and a scarcity of macrofossils suitable for dating. Here we present the results of multiple dating methods carried out on cores in the Gotland Deep area of the Baltic Sea. Particular emphasis is placed on the Littorina stage (8 ka ago to the present) of the Baltic Sea and possible changes in the 14C reservoir age of our dated samples. Three geochronological methods are used. Firstly, palaeomagnetic secular variations (PSV) are reconstructed, whereby ages are transferred to PSV features through comparison with varved lake sediment based PSV records. Secondly, lead (Pb) content and stable isotope analysis are used to identify past peaks in anthropogenic atmospheric Pb pollution. Lastly, 14C determinations were carried out on benthic foraminifera (Elphidium spec.) samples from the brackish Littorina stage of the Baltic Sea. Determinations carried out on smaller samples (as low as 4 µg C) employed an experimental, state-of-the-art method involving the direct measurement of CO2 from samples by a gas ion source without the need for a graphitisation step - the first time this method has been performed on foraminifera in an applied study. The PSV chronology, based on the uppermost Littorina stage sediments, produced ten age constraints between 6.29 and 1.29 cal ka BP, and the Pb depositional analysis produced two age constraints associated with the Medieval pollution peak. Analysis of PSV data shows that adequate directional data can be derived from both the present Littorina saline phase muds and Baltic Ice Lake stage varved glacial sediments. Ferrimagnetic iron sulphides, most likely authigenic greigite (Fe3S4), present in the intermediate Ancylus Lake freshwater stage sediments acquire a gyroremanent magnetisation during static alternating field (AF) demagnetisation, preventing the identification of a primary natural remanent magnetisation for these sediments. An inferred marine reservoir age offset (deltaR) is calculated by comparing the foraminifera 14C determinations to a PSV & Pb age model. This deltaR is found to trend towards younger values upwards in the core, possibly due to a gradual change in hydrographic conditions brought about by a reduction in marine water exchange from the open sea due to continued isostatic rebound.

Abundance of zooplankton based on a long-term study at the Galata transect and covers the spring-summer periods from 1967 till 2005

The dataset is based on a long-term study (38 years) at the Galata transect and covers the spring-summer periods from 1967 till 2005. The whole dataset is composed of 360 data of total zooplankton biomass and abundance . Samples were collected in discrete layers 0-10m, 10-20m, 10-25m, 25-50m, 50-70m, 50-100m, 100-150. Mesozooplankton abundance: the collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber for taxomomic identification and count. Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Fishery Resource by Prof. Asen Konsulov and Institute of Oceanology by Prof. Asen Konsulov, Lyudmila Kamburska and Kremena Stefanova using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972). Taxon-specific mesozooplankton abundance: The collected material was analysed using the method of Domov (1959). Samples were brought to volume of 25-30 ml depending upon zooplankton density and mixed intensively until all organisms were distributed randomly in the sample volume. After that 5 ml of sample was taken and poured in the counting chamber for taxomomic identification and count. Copepods and Cladoceras were identified and enumerated; the other mesozooplankters were identified and enumerated at higher taxonomic level (commonly named as mesozooplankton groups). Large (> 1 mm body length) and not abundant species were calculated in whole sample. Counting and measuring of organisms were made in the Dimov chamber under the stereomicroscope to the lowest taxon possible. Taxonomic identification was done at the Institute of Fishery Resource by prof. Asen Konsulov and Institute of Oceanology by Prof. Asen Konsulov, Lyudmila Kamburska and Kremena Stefanova using the relevant taxonomic literature (Mordukhay-Boltovskoy, F.D. (Ed.). 1968, 1969,1972).

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

(Figure 3d) Export production reconstruction based on planktonic foraminiferal sampled in surface sediments

To reconstruct Export Productivity (Pexp), 27 taxonomic categories of the planktonic foraminifera census data were used with the modern analog technique SIMMAX 28 (Pflaumann et al., 1996, doi:10.1029/95PA01743; 2003, doi:10.1029/2002PA000774). To the 26 taxonomic groups widely used and listed in Kucera et al. (2005, doi:10.1016/j.quascirev.2004.07.014), Turborotalita humilis was added in our calibration as it is associated with the PCC source region (Meggers et al., 2002, doi:10.1016/S0967-0645(02)00103-0). The modern analog file is based on the Iberian margin database (Salgueiro et al., 2008, doi:10.1016/j.marmicro.2007.09.003) combined with the North Atlantic surface samples used by the MARGO project (Kucera et al., 2005). This results in a total of 999 analogs for Pexp. Modern oceanic primary productivity (PP) is obtained for each site by averaging 12 monthly primary productivity values for a 8-year period (1978-1986) that were estimated from satellite color data (CZCS) and gridded at 0.5° latitude - longitude fields (Antoine et al., 1996, doi:10.1029/95GB02832). Export Productivity (Pexp) was calculated from the PP values following the empirical relationship Pexp = PP**2/400 for primary production below 200 gC/m**2/yr, and Pexp = PP/2 for primary production above 200 gC/m2/yr (Eppley and Peterson, 1979, doi:10.1038/282677a0; Sarnthein et al., 1988, doi:10.1029/PA003i003p00361). The residuals gives the differences between satellite based Pexp and foraminiferal Pexp.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

Ground-based electromagnetic (EM) and drill-hole ice and snow thickness and melt pond depth measurements during Polarstern cruise ARK-IX/4, ARK-XI/1, and ARK-XII

Drillhole-determined sea-ice thickness was compared with values derived remotely using a portable small-offset loop-loop steady state electromagnetic (EM) induction device during expeditions to Fram Strait and the Siberian Arctic, under typical winter and summer conditions. Simple empirical transformation equations are derived to convert measured apparent conductivity into ice thickness. Despite the extreme seasonal differences in sea-ice properties as revealed by ice core analysis, the transformation equations vary little for winter and summer. Thus, the EM induction technique operated on the ice surface in the horizontal dipole mode yields accurate results within 5 to 10% of the drillhole determined thickness over level ice in both seasons. The robustness of the induction method with respect to seasonal extremes is attributed to the low salinity of brine or meltwater filling the extensive pore space in summer. Thus, the average bulk ice conductivity for summer multiyear sea ice derived according to Archie's law amounts to 23 mS/m compared to 3 mS/m for winter conditions. These mean conductivities cause only minor differences in the EM response, as is shown by means of 1-D modeling. However, under summer conditions the range of ice conductivities is wider. Along with the widespread occurrence of surface melt ponds and freshwater lenses underneath the ice, this causes greater scatter in the apparent conductivity/ice thickness relation. This can result in higher deviations between EM-derived and drillhole determined thicknesses in summer than in winter.

The Lena River Delta - Land cover classification of tundra environments based on Landsat 7 ETM+ data and its application for upscaling of methane emissions

The Lena River Delta, situated in Northern Siberia (72.0 - 73.8° N, 122.0 - 129.5° E), is the largest Arctic delta and covers 29,000 km**2. Since natural deltas are characterised by complex geomorphological patterns and various types of ecosystems, high spatial resolution information on the distribution and extent of the delta environments is necessary for a spatial assessment and accurate quantification of biogeochemical processes as drivers for the emission of greenhouse gases from tundra soils. In this study, the first land cover classification for the entire Lena Delta based on Landsat 7 Enhanced Thematic Mapper (ETM+) images was conducted and used for the quantification of methane emissions from the delta ecosystems on the regional scale. The applied supervised minimum distance classification was very effective with the few ancillary data that were available for training site selection. Nine land cover classes of aquatic and terrestrial ecosystems in the wetland dominated (72%) Lena Delta could be defined by this classification approach. The mean daily methane emission of the entire Lena Delta was calculated with 10.35 mg CH4/m**2/d. Taking our multi-scale approach into account we find that the methane source strength of certain tundra wetland types is lower than calculated previously on coarser scales.

Table 1. Tree growth index and reconstructed precipitation series from 515 BC to AD 1993 based on the archaeological and living tree ring datasets from Dulan and Shenge, China

Annual precipitation for the last 2,500 years was reconstructed for northeastern Qinghai from living and archaeological juniper trees. A dominant feature of the precipitation of this area is a high degree of variability in mean rainfall at annual, decadal, and centennial scales, with many wet and dry periods that are corroborated by other paleoclimatic indicators. Reconstructed values of annual precipitation vary mostly from 100 to 300 mm and thus are no different from the modern instrumental record in Dulan. However, relatively dry years with below-average precipitation occurred more frequently in the past than in the present. Periods of relatively dry years occurred during 74-25 BC, AD 51-375, 426-500, 526-575, 626-700, 1100-1225, 1251-1325, 1451-1525, 1651-1750 and 1801-1825. Periods with a relatively wet climate occurred during AD 376-425, 576-625, 951-1050, 1351-1375, 1551-1600 and the present. This variability is probably related to latitudinal positions of winter frontal storms. Another key feature of precipitation in this area is an apparently direct relationship between interannual variability in rainfall with temperature, whereby increased warming in the future might lead to increased flooding and droughts. Such increased climatic variability might then impact human societies of the area, much as the climate has done for the past 2,500 years.