Synthesized grounding line and ice shelf mask for Antarctica

Iceberg calving has been assumed to be the dominant cause of mass loss for the Antarctic ice sheet, with previous estimates of the calving flux exceeding 2,000 gigatonnes per year. More recently, the importance of melting by the ocean has been demonstrated close to the grounding line and near the calving front. So far, however, no study has reliably quantified the calving flux and the basal mass balance (the balance between accretion and ablation at the ice-sheet base) for the whole of Antarctica. The distribution of fresh water in the Southern Ocean and its partitioning between the liquid and solid phases is therefore poorly constrained. Here we estimate the mass balance components for all ice shelves in Antarctica, using satellite measurements of calving flux and grounding-line flux, modelled ice-shelf snow accumulation rates and a regional scaling that accounts for unsurveyed areas. We obtain a total calving flux of 1,321 ± 144 gigatonnes per year and a total basal mass balance of -1,454 ± 174 gigatonnes per year. This means that about half of the ice-sheet surface mass gain is lost through oceanic erosion before reaching the ice front, and the calving flux is about 34 per cent less than previous estimates derived from iceberg tracking. In addition, the fraction of mass loss due to basal processes varies from about 10 to 90 per cent between ice shelves. We find a significant positive correlation between basal mass loss and surface elevation change for ice shelves experiencing surface lowering and enhanced discharge. We suggest that basal mass loss is a valuable metric for predicting future ice-shelf vulnerability to oceanic forcing.

Digital elevation model (DEM), surveyed stream network (SSN), mask and stream sampling points (SSP)

These data are provided to allow users for reproducibility of an open source tool entitled 'automated Accumulation Threshold computation and RIparian Corridor delineation (ATRIC)'

Radiance, irradiance, and remote sensing reflectance during the North Sea Coast Harmful Algal Bloom (NORCOHAB II) HEINCKE cruise HE302

In this study four data quality flags are presented for automated and unmanned above-water hyperspectral optical measurements collected underway in the North Sea, The Minch, Irish Sea and Celtic Sea in April/May 2009. Coincident to these optical measurements a DualDome D12 (Mobotix, Germany) camera system was used to capture sea surface and sky images. The first three flags are based on meteorological conditions, to select erroneous incoming solar irradiance (ES) taken during dusk, dawn, before significant incoming solar radiation could be detected or under rainfall. Furthermore, the relative azimuthal angle of the optical sensors to the sun is used to identify possible sunglint free sea surface zones. A total of 629 spectra remained after applying the meteorological masks (first three flags). Based on this dataset, a fourth flag for sunglint was generated by analysing and evaluating water leaving radiance (LW) and remote sensing reflectance (RRS) spectral behaviour in the presence and absence of sunglint salient in the simultaneously available sea surface images. Spectra conditions satisfying "mean LW (700-950 nm) < 2 mW/m**2/nm/Sr" or alternatively "minimum RRS (700-950 nm) < 0.010/Sr", mask the most measurements affected by sunglint, providing efficient flagging of sunglint in automated quality control. It is confirmed that valid optical measurements can be performed 0° <= theta <= 360° although 90° <= theta <= 135° is recommended.

Sediment geochemistry of ODP Leg 127 sites

The relative effects of paleoceanographic and paleogeographic variations, sediment lithology, and diagenetic processes on the final preserved chemistry of Japan Sea sediments are evaluated by investigating the rare earth element (REE), major element, and trace element concentrations in 59 squeeze-cake whole-round and 27 physical-property sample residues from Sites 794, 795, and 797, cored during ODP Leg 127. The most important variation in sedimentary chemical composition is the increase in SiO2 concentration through the Pliocene diatomaceous sequences, which dilutes most other major and trace element components by various degrees. This biogenic input is largest at Site 794 (Yamato Basin), moderately developed at Site 797 (Yamato Basin), and of only minor importance at Site 795 (Japan Basin), potentially reflecting basinal contrasts in productivity with the Yamato Basin recording greater biogenic input than the Japan Basin and with the easternmost sequence of Site 794 lying beneath the most productive waters. There are few systematic changes in solid-phase chemistry resulting from the opal-A/opal-CT or opal-CT/quartz silica phase transformations. Most major and trace element concentrations are controlled by the aluminosilicate fraction of the sediment, although the effects of diagenetic silica phases and manganese carbonates are of localized importance. REE total abundances (Sum REE) in the Japan Sea are strongly dependent upon the paleoceanographic position of a given site with respect to terrigenous and biogenic sources. REE concentrations at Site 794 overall correspond well to aluminosilicate chemical indices and are strongly diluted by SiO2 within the upper Miocene-Pliocene diatomaceous sequence. Eu/Eu* values at Site 794 reach a maximum through the diatomaceous interval as well, most likely suggesting an association of Eu/Eu* with the siliceous component, or reflecting slight incorporation of a detrital feldspar phase. SumREE at Site 795 also is affiliated strongly with aluminosilicate phases and yet is diluted only slightly by siliceous input. At Site 797, SumREE is not as clearly associated with the aluminosilicate fraction, is correlated moderately to siliceous input, and may be sporadically influenced by detrital heavy minerals originating from the nearby rifted continental fragment composing the Yamato Rise. Ce/Ce* profiles at all three sites increase essentially monotonically with depth and record progressive diagenetic LREE fractionation. The observed Ce/Ce* increases are not responding to changes in the paleoceanographic oxygenation state of the overlying water, as there is no independent evidence to suggest the proper oceanographic conditions. Ce/Ce* correlates slightly better with depth than with age at the two Yamato Basin sites. The downhole increase in Ce/Ce* at Sites 794 and 797 is a passive response to the diagenetic transfer of LREE (except Ce) from sediment to interstitial water. At Site 795, the overall lack of correlation between Ce/Ce* and Lan/Ybn suggests that other processes mask the diagenetic behavior of all LREEs. First-order calculations of the Ce budget in Japan Sea waters and sediment indicate that ~20% of the excess Ce adsorbed by settling particles is recycled within the water column and that an additional ~38% is recycled at or near the seafloor. Thus, because the remaining excess Ce is only ~10% of the total Ce, there is not a large source of Ce to the deeply buried sediment, further suggesting that the downhole increase in Ce/Ce* is a passive response to diagenetic behavior of the other LREEs. The REE chemistry of Japan Sea sediment therefore predicts successive downhole addition of LREEs to deeply buried interstitial waters.

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.

Paleocene and Eocene calcareous nannofossils at DSDP Legs 25 and 40 Holes

Paleocene and Eocene nannofossil flora from Deep Sea Drilling Project Legs 25 and 40 were analyzed in order to provide a basis of comparison with DSDP Legs 36 and 71 and with other South Atlantic assemblages. A mid-latitude biostratigraphic zonation, using previously described zonal markers, was adopted for the southwest Indian Ocean. Various diagenetic effects were noted in the sedimentary sequences. Some of these mask to some extent paleoecologic signals, particularly those generated by the Discoaster/Chiasmolithus ratio.

Rare earth element data for sediments from ODP Leg 127

The relative effects of paleoceanographic and paleogeographic variations, sediment lithology, and diagenetic processes on the recorded rare earth element (REE) chemistry of Japan Sea sediments are evaluated by investigating REE total abundances and relative fractionations in 59 samples from Ocean Drilling Program Leg 127. REE total abundances (Sum REE) in the Japan Sea are strongly dependent upon the paleoceanographic position of a given site with respect to terrigenous and biogenic sources. REE concentrations at Site 794 (Yamato Basin) overall correspond well to aluminosilicate chemical indices and are strongly diluted by SiO2 within the late Miocene-Pliocene diatomaceous sequence. Eu/Eu* values at Site 794 reach a maximum through the diatomaceous interval as well, most likely suggesting an association of Eu/Eu* with the siliceous component, or reflecting slight incorporation of a detrital feldspar phase. Sum REE at Site 795 (Japan Basin) also is affiliated strongly with aluminosilicate phases, yet is diluted only slightly by siliceous input. At Site 797 (Yamato Basin), REE is not as clearly associated with the aluminosilicate fraction, is correlated moderately to siliceous input, and may be sporadically influenced by detrital heavy minerals originating from the nearby rifted continental fragment composing the Yamato Rise. The biogenic influence is largest at Site 794, moderately developed at Site 797, and of only minor importance at Site 795, reflecting basinal contrasts in productivity such that the Yamato Basin records greater biogenic input than the Japan Basin, while the most productive waters overlie the easternmost sequence of Site 794. Ce/Ce* profiles at all three sites increase monotonically with depth, and record progressive diagenetic LREE fractionation. The observed Ce/Ce* record does not respond to changes in oxygenation state of the overlying water, and Ce/Ce* correlates slightly better with depth than with age. The downhole increase in Ce/Ce* at Site 794 and Site 797 is a passive response to diagenetic transfer of LREE (except Ce) from sediment to interstitial water. At Site 795, the overall lack of correlation between Ce/Ce* and La_n/Yb_n suggests that other processes are occurring which mask the diagenetic behavior of all LREEs. First-order calculations of the Ce budget in Japan Sea waters and sediment indicate that ~20% of the excess Ce adsorbed by settling particles is recycled within the water column, and that an additional ~38% is recycled at or near the seafloor (data from Masuzawa and Koyama, 1989). Thus, because the remaining excess Ce is only ~10% of the total Ce, there is not a large source of Ce to the deeply buried sediment, further suggesting that the downhole increase in Ce/Ce* is a passive response to diagenetic behavior of the other LREEs. The REE chemistry of Japan Sea sediment therefore predicts successive downhole addition of LREEs to deeply-buried interstitial waters.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Iron isotope and trace metal record of Cenomanian-Turonian sediments

he global carbon cycle during the mid-Cretaceous (~125-88 million years ago, Ma) experienced numerous major perturbations linked to increased organic carbon burial under widespread, possibly basin-scale oxygen deficiency and episodes of euxinia (anoxic and H2S-containing). The largest of these episodes, the Cenomanian-Turonian boundary event (ca. 93.5 Ma), or oceanic anoxic event (OAE) 2, was marked by pervasive deposition of organic-rich, laminated black shales in deep waters and in some cases across continental shelves. This deposition is recorded in a pronounced positive carbon isotope excursion seen ubiquitously in carbonates and organic matter. Enrichments of redox-sensitive, often bioessential trace metals, including Fe and Mo, indicate major shifts in their biogeochemical cycles under reducing conditions that may be linked to changes in primary production. Iron enrichments and bulk Fe isotope compositions track the sources and sinks of Fe in the proto-North Atlantic at seven localities marked by diverse depositional conditions. Included are an ancestral mid-ocean ridge and euxinic, intermittently euxinic, and oxic settings across varying paleodepths throughout the basin. These data yield evidence for a reactive Fe shuttle that likely delivered Fe from the shallow shelf to the deep ocean basin, as well as (1) hydrothermal sources enhanced by accelerated seafloor spreading or emplacement of large igneous province(s) and (2) local-scale Fe remobilization within the sediment column. This study, the first to explore Fe cycling and enrichment patterns on an ocean scale using iron isotope data, demonstrates the complex processes operating on this scale that can mask simple source-sink relationships. The data imply that the proto-North Atlantic received elevated Fe inputs from several sources (e.g., hydrothermal, shuttle and detrital inputs) and that the redox state of the basin was not exclusively euxinic, suggesting previously unknown heterogeneity in depositional conditions and biogeochemical cycling within those settings during OAE-2.

Occurrence of nannoplankton observed in Early Pliocene samples from Zanclean stratotype of Capo Rossello, Sicily (Table I)

Calcareous nannoplankton biostratigraphy has been worked out in the eastern Mediterranean utilizing deep-sea sediments recovered from DSDP Leg 42A Sites 375 and 376. These two drill sites were located approximately 55 km west of Cyprus on the Florence Rise. Sediments, ranging in age from early Miocene (Helicosphaera ampliaperta Zone) through Holocene, contain sufficient age-diagnostic species to recognize essentially all of the lowlatitude nannoplankton zones described by Bukry, although regional, secondary marker species are needed to define some zonal boundaries. Reworked Cretaceous and Paleogene nannoplankton occur throughout the stratigraphic interval studied, but not in quantities large enough to mask indigenous species. Sedimentation rates at Sites 375 and 376 were highest in the late Miocene and late Pleistocene. Open-marine, warm-water species of discoasters are present in significant numbers throughout the Miocene and Pliocene. Earliest Pliocene assemblages contain numerous specimens of ceratoliths. Nannoplankton in post-Messinian sediments at the drill sites and the Zanclean stratotype at Capo Rossello, Sicily, indicate that the base of the Amaurolithus tricorniculatus Zone (base of Triquetrorhabdulus rugosus Subzone) corresponds with the Miocene-Pliocene boundary.

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.

Global 2009 1km MODIS (MOD35/MOD09) cloud frequency and MOD35 processing path

Identifying cloud interference in satellite-derived data is a critical step toward developing useful remotely sensed products. Most MODIS land products use a combination of the MODIS (MOD35) cloud mask and the 'internal' cloud mask of the surface reflectance product (MOD09) to mask clouds, but there has been little discussion of how these masks differ globally. We calculated global mean cloud frequency for both products, for 2009, and found that inflated proportions of observations were flagged as cloudy in the Collection 5 MOD35 product. These erroneously categorized areas were spatially and environmentally non-random and usually occurred over high-albedo land-cover types (such as grassland and savanna) in several regions around the world. Additionally, we found that spatial variability in the processing path applied in the Collection 5 MOD35 algorithm affects the likelihood of a cloudy observation by up to 20% in some areas. These factors result in abrupt transitions in recorded cloud frequency across landcover and processing-path boundaries impeding their use for fine-scale spatially contiguous modeling applications. We show that together, these artifacts have resulted in significantly decreased and spatially biased data availability for Collection 5 MOD35-derived composite MODIS land products such as land surface temperature (MOD11) and net primary productivity (MOD17). Finally, we compare our results to mean cloud frequency in the new Collection 6 MOD35 product, and find that landcover artifacts have been reduced but not eliminated. Collection 6 thus increases data availability for some regions and land cover types in MOD35-derived products but practitioners need to consider how the remaining artifacts might affect their analysis.

Foraminifera in sediment cores of the Lomonosov Ridge and Morris Jesup Rise

X-ray fluorescence (XRF) scanning of sediment cores from the Lomonosov Ridge and the Morris Jesup Rise reveals a distinct pattern of Ca intensity peaks through Marine Isotope Stages (MIS) 1 to 7. Downcore of MIS 7, the Ca signal is more irregular and near the detection limit. Virtually all major peaks in Ca coincide with a high abundance of calcareous microfossils; this is particularly conspicuous in the cores from the central Arctic Ocean. However, the recorded Ca signal is generally caused by a combination of biogenic and detrital carbonate, and in areas influenced by input from the Canadian Arctic, detrital carbonates may effectively mask the foraminiferal carbonates. Despite this, there is a strong correlation between XRF-detected Ca content and foraminiferal abundance. We propose that in the Arctic Ocean north of Greenland a common palaeoceanographic mechanism is controlling Ca-rich ice-rafted debris (IRD) and foraminiferal abundance. Previous studies have shown that glacial periods are characterized by foraminfer-barren sediments. This implies that the Ca-rich IRD intervals with abundant foraminifera were most likely deposited during interglacial periods when glaciers left in the Canadian Arctic Archipelago were still active and delivered a large amount of icebergs. At the same time, conditions were favourable for planktic foraminifera, resulting in a strong covariance between these proxies. Therefore, we suggest that the XRF scanner's capability to efficiently map Ca concentrations in sediment cores makes it possible to systematically examine large numbers of cores from different regions to investigate the palaeoceanographic reasons for the calcareous microfossils' spatial and temporal variability.