Data and Code archive for the interpolation of surface ocean carbon dioxide

We have developed a statistical gap-filling method adapted to the specific coverage and properties of observed fugacity of surface ocean CO2 (fCO2). We have used this method to interpolate the Surface Ocean CO2 Atlas (SOCAT) v2 database on a 2.5°×2.5° global grid (south of 70°N) for 1985-2011 at monthly resolution. The method combines a spatial interpolation based on a 'radius of influence' to determine nearby similar fCO2 values with temporal harmonic and cubic spline curve-fitting, and also fits long term trends and seasonal cycles. Interannual variability is established using deviations of observations from the fitted trends and seasonal cycles. An uncertainty is computed for all interpolated values based on the spatial and temporal range of the interpolation. Tests of the method using model data show that it performs as well as or better than previous regional interpolation methods, but in addition it provides a near-global and interannual coverage.

High resolution bathymetric compilation for Potter Cove, WAP, Antarctica, with links to data in ArcGIS format

The bathymetry raster with a resolution of 5 m x 5 m was processed from unpublished single beam data from the Argentine Antarctica Institute (IAA, 2010) and multibeam data from the United Kingdom Hydrographic Office (UKHO, 2012) with a cell size of 5 m x 5 m. A coastline digitized from a satellite image (DigitalGlobe, 2014) supplemented the interpolation process. The 'Topo to Raster' tool in ArcMap 10.3 was used to merge the three data sets, while the coastline represented the 0-m-contour to the interpolation process ('contour type option').

Palaeolake Rehwiese, NE-Germany, Younger Dryas varved data set RW09

The first 1400-year floating varve chronology for north-eastern Germany covering the late Allerød to the early Holocene has been established by microscopic varve counts from the Rehwiese palaeolake sediment record. The Laacher See Tephra (LST), at the base of the studied interval, forms the tephrochronological anchor point. The fine laminations were examined using a combination of micro-facies and ?-XRF analyses and are typical of calcite varves, which in this case provide mainly a warm season signal. Two varve types with different sub-layer structures have been distinguished: (I) complex varves consisting of up to four seasonal sub-layers formed during the Allerød and early Holocene periods, and, (II) simple two sub-layer type varves only occurring during the Younger Dryas. The precision of the chronology has been improved by varve-to-varve comparison of two independently analyzed sediment profiles based on well-defined micro-marker layers. This has enabled both (1) the precise location of single missing varves in one of the sediment profiles, and, (2) the verification of varve interpolation in disturbed varve intervals in the parallel core. Inter-annual and decadal-scale variability in sediment deposition processes were traced by multi-proxy data series including seasonal layer thickness, high-resolution element scans and total organic and inorganic carbon data at a five-varve resolution. These data support the idea of a two-phase Younger Dryas, with the first interval (12,675 - 12,275 varve years BP) characterised by a still significant but gradually decreasing warm-season calcite precipitation and a second phase (12,275 - 11,640 varve years BP) with only weak calcite precipitation. Detailed correlation of these two phases with the Meerfelder Maar record based on the LST isochrone and independent varve counts provides clues about regional differences and seasonal aspects of YD climate change along a transect from a location proximal to the North Atlantic in the west to a more continental site in the east.

Age model, raw data and interpolated raw data from proxy records of sediment cores GeoB6007-1 and GeoB6007-2

Here we present a 1200 yr long benthic foraminiferal Mg/Ca based temperature and oxygen isotope record from a ~900 m deep sediment core off northwest Africa to show that atmosphere-ocean interactions in the eastern subpolar gyre are transferred at central water depth into the eastern boundary of the subtropical gyre. Further we link the variability of the NAO (over the past 165 yrs) and solar irradiance (Late Holocene) and their control on subpolar mode water formation to the multidecadal variability observed at mid-depth in the eastern subtropical gyre. Our results show that eastern North Atlantic central waters cooled by up to ~0.8± 0.7 °C and densities decreased by Sigma theta=0.3±0.2 during positive NAO years and during minima in solar irradiance during the Late Holocene. The presented records demonstrate the sensitivity of central water formation to enhanced atmospheric forcing and ice/freshwater fluxes into the eastern subpolar gyre and the importance of central water circulation for cross-gyre climate signal propagation during the Late Holocene.

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.

(Table 3) Carbon and oxygen isotopic data of carbonates from Core SLR 37-1, HYC 25-1 and an aragonite crust from station DR 35-1

An area of massive barite precipitations was studied at a tectonic horst in 1500 m water depth in the Derugin Basin, Sea of Okhotsk. Seafloor observations and dredge samples showed irregular, block- to column-shaped barite build-ups up to 10 m high which were scattered over the seafloor along an observation track 3.5 km long. High methane concentrations in the water column show that methane expulsion and probably carbonate precipitation is a recently active process. Small fields of chemoautotrophic clams (Calyptogena sp., Acharax sp.) at the seafloor provide additional evidence for active fluid venting. The white to yellow barites show a very porous and often layered internal fabric, and are typically covered by dark-brown Mn-rich sediment; electron microprobe spectroscopy measurements of barite sub-samples show a Ba substitution of up to 10.5 mol% of Sr. Rare idiomorphic pyrite crystals (1%) in the barite fabric imply the presence of H2S. This was confirmed by clusters of living chemoautotrophic tube worms (1 mm in diameter) found in pores and channels within the barite. Microscopic examination showed that micritic aragonite and Mg-calcite aggregates or crusts are common authigenic precipitations within the barite fabric. Equivalent micritic carbonates and barite carbonate cemented worm tubes were recovered from sediment cores taken in the vicinity of the barite build-up area. Negative ?13C values of these carbonates (>?43.5? PDB) indicate methane as major carbon source; ?18O values between 4.04 and 5.88? PDB correspond to formation temperatures, which are certainly below 5°C. One core also contained shells of Calyptogena sp. at different core depths with 14C-ages ranging from 20 680 to >49 080 yr. Pore water analyses revealed that fluids also contain high amounts of Ba; they also show decreasing SO42- concentrations and a parallel increase of H2S with depth. Additionally, S and O isotope data of barite sulfate (?34S: 21.0-38.6? CDT; ?18O: 9.0-17.6? SMOW) strongly point to biological sulfate reduction processes. The isotope ranges of both S and O can be exclusively explained as the result of a mixture of residual sulfate after a biological sulfate reduction and isotopic fractionation with 'normal' seawater sulfate. While massive barite deposits are commonly assumed to be of hydrothermal origin, the assemblage of cheomautotrophic clams, methane-derived carbonates, and non-thermally equilibrated barite sulfate strongly implies that these barites have formed at ambient bottom water temperatures and form the features of a Giant Cold Seep setting that has been active for at least 49 000 yr.

Spatial distribution of key zooplankton species using continuous plankton recorder (CPR) data from the North Atlantic (2000-2009)

The continuous plankton recorder (CPR) survey is an upper layer plankton monitoring program that has regularly collected samples, at monthly intervals, in the North Atlantic and adjacent seas since 1946. Water from approximately 6 m depth enters the CPR through a small aperture at the front of the sampler and travels down a tunnel where it passes through a silk filtering mesh of 270 µm before exiting at the back of the CPR. The plankton filtered on the silk is analyzed in sections corresponding to 10 nautical miles (approx. 3 m**3 of seawater filtered) and the plankton microscopically identified (Richardson et al., 2006 and reference therein). In the present study we used the CPR data to investigate the current basin scale distribution of C. finmarchicus (C5-C6), C. helgolandicus (C5-C6), C. hyperboreus (C5-C6), Pseudocalanus spp. (C6), Oithona spp. (C1-C6), total Euphausiida, total Thecosomata and the presence/absence of Cnidaria and the Phytoplankton Colour Index (PCI). The PCI, which is a visual assessment of the greenness of the silk, is used as an indicator of the distribution of total phytoplankton biomass across the Atlantic basin (Batten et al., 2003). Monthly data collected between 2000 and 2009 were gridded using the inverse-distance interpolation method, in which the interpolated values were the nodes of a 2 degree by 2 degree grid. The resulting twelve monthly matrices were then averaged within the year and in the case of the zooplankton the data were log-transformed (i.e. log10 (x+1).

Global data compilation of benthic data sets II

In this study we present a global distribution pattern and budget of the minimum flux of particulate organic carbon to the sea floor (J POC alpha). The estimations are based on regionally specific correlations between the diffusive oxygen flux across the sediment-water interface, the total organic carbon content in surface sediments, and the oxygen concentration in bottom waters. For this, we modified the principal equation of Cai and Reimers [1995] as a basic monod reaction rate, applied within 11 regions where in situ measurements of diffusive oxygen uptake exist. By application of the resulting transfer functions to other regions with similar sedimentary conditions and areal interpolation, we calculated a minimum global budget of particulate organic carbon that actually reaches the sea floor of ~0.5 GtC yr**-1 (>1000 m water depth (wd)), whereas approximately 0.002-0.12 GtC yr**-1 is buried in the sediments (0.01-0.4% of surface primary production). Despite the fact that our global budget is in good agreement with previous studies, we found conspicuous differences among the distribution patterns of primary production, calculations based on particle trap collections of the POC flux, and J POC alpha of this study. These deviations, especially located at the southeastern and southwestern Atlantic Ocean, the Greenland and Norwegian Sea and the entire equatorial Pacific Ocean, strongly indicate a considerable influence of lateral particle transport on the vertical link between surface waters and underlying sediments. This observation is supported by sediment trap data. Furthermore, local differences in the availability and quality of the organic matter as well as different transport mechanisms through the water column are discussed.