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.

Pollen, DNA and Vegetation records from surface sediments of lacustrine lakes along a forest tundra - tundra transect at the Taymyr peninsula, northern Siberia

Reliable information of past vegetation changes are important to project future changes, especially for areas undergoing rapid transitioning such as the boreal treeline. The application of detailed sedDNA records has the potential to enhance our understanding of vegetation changes gained mainly from pollen studies of lake sediments. This study investigates sedDNA and pollen records from 31 lakes along a gradient of increasing larch forest cover in northern Siberia (Taymyr Peninsula) and compares them with vegetation field surveys within the lake's catchment. With respect to vegetation richness, sedDNA recorded 114 taxa, about half of them to species level, while pollen analyses identified 43 pollen taxa. Both approaches exceed the 31 taxa revealed by vegetation field surveys of 400 m**2 plots. From north to south, Larix percentages increase, as is consistently recorded by all three methods. Furthermore, tundra sites are separated from forested sites in the plots of the principal component analyses. Comparison of ordination results by Procrustes and Protest analyses yields a significant fit among all compared pairs of records. Despite the overall comparability of sedDNA and pollen analyses certain idiosyncrasies in the compositional signal are observed, such as high percentages of Alnus and Betula in all pollen spectra and high percentages of Salix in all sedDNA spectra. In conclusion, our results from the treeline show that sedDNA analyses perform better than pollen in recording site-specific richness (i.e. presence/absence of certain vegetation taxa in the direct vicinity of the lake) and perform as good as pollen in tracing regional vegetation composition.

Total organic carbon in surface sediments of the Ob and Yenisei estuaries and adjacent coastal areas, appendix 1

Two main mechanisms are controlling the accumulation of organic matter in the sediments of the Kara Sea. The large rivers Ob and Yenisei supply significant quantities of freshwater onto the shelf (Lisitsyn and Vinogradov, 1995; Bobrovitskaya et al., 1996; Johnson et al., 1997) and deliver terrigenous organie matter and aquatic algae. Additionally, marine organic matter is produced in the water column. In order to distinguish between the different sources of the organic material maceral analysis, organic-geochemical bulk Parameters and biomarkers (short- and long-chain D-alkanes, fatty acids and pigments) were used to determine the quality (marine vs. terrigenous) and quantity of the organic carbon fraction in the surface sediments taken during the 28th cruise of RV Akademik Boris Petrov (Matthiessen and Stepanets, 1998) (Fig. 1). Previous organic-geochemical investigations (i.e., total organic-carbon content (TOC), hydrogen indices (Hl), CIN-ratios) indicate the importance of terrigenous input of organic matter (Galimov et al., 1996; Stein, 1996). Studies of lipid biomarkers in surface sediments in the Ob estuary show also a predominance of terrestrial constituents and an increase in planktonogenic and bacterial lipids further offshore (Belyaeva and Eglinton, 1997). In complex systems such as the Eurasian continental margin characterized by high input of terrestriallaquatic organic matter and strong seasonal variation in sea-ice Cover and primary productivity, the Interpretation of the organic geochemical data is much more complicated and restricted in comparison to similar data Sets from low-latitude open-ocean environments (Fahl and Stein, 1998). Microscopical studies (maceral analysisl palynology), however, allow a direct visual inspection of the particulate organic matter and allow to differentiate particles of different biological sources. Thus, a combination of both methods as shown in this study, yields a more precise identification of organic-carbon sources.

Late Pliocene sea-surface temperature record of DSDP Hole 80-548 from the Northeast Atlantic Ocean (Table 1)

Application of quantitative micropaleontologic methods of sea-surface temperature (SST) estimation to the planktic foraminiferal faunas at Hole 548 in the northeast Atlantic between 3.5 and 2.0 Ma provide a detailed record of Pliocene surface temperature changes during the onset of major Northern Hemisphere glaciation. SST estimates indicate a relatively warm interval of "enhanced interglacials (warm peaks)" between 3.2 and 2.8 Ma followed by steady cooling with a significant step-like change toward colder temperatures occurring about 2.4 Ma. This record compares well with faunal and SST records from other North Atlantic sites.

Properties of seawater and particulate matter from a WETLabs AC-S spectrophotometer and a WETLabs chlorophyll fluorometer mounted on the continuous surface water sampling system during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) was a global survey of ocean ecosystems aboard the Sailing Vessel Tara. It carried out extensive measurements of evironmental conditions and collected plankton (viruses, bacteria, protists and metazoans) for later analysis using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data publication provides permanent links to original and updated versions of validated data files containing measurements from the Continuous Surface Sampling System [CSSS]. Water was pumped at the front of the vessel from ~2m depth, then de-bubbled and circulated to a WETLabs AC-S spectrophotometer and a WETLabs chlorophyll fluorometer. Systems maintenance (instrument cleaning, flushing) was done approximately once a week and in port between successive legs. All data were stamped with a GPS.

Estimation of sea surface temperature from sediment core PS1768-8

The quantitative diatom analysis of 218 surface sediment samples recovered in the Atlantic and western Indian sector of the Southern Ocean is used to define a base of reference data for paleotemperature estimations from diatom assemblages using the Imbrie and Kipp transfer function method. The criteria which justify the exclusion of samples and species out of the raw data set in order to define a reference database are outlined and discussed. Sensitivity tests with eight data sets were achieved evaluating the effects of overall dominance of single species, different methods of species abundance ranking, and no-analog conditions (e.g., Eucampia Antarctica) on the estimated paleotemperatures. The defined transfer functions were applied on a sediment core from the northern Antarctic zone. Overall dominance of Fragilariopsis kerguelensis in the diatom assemblages resulted in a close affinity between paleotemperature curve and relative abundance pattern of this species downcore. Logarithmic conversion of counting data applied with other ranking methods in order to compensate the dominance of F. kerguelensis revealed the best statistical results. A reliable diatom transfer function for future paleotemperature estimations is presented.

Planktonic foraminiferal faunas from surface-sediment samples from the South China Sea and the western Pacific

We present 30 new planktonic foraminiferal census data of surface sediment samples from the South China Sea, recovered between 630 and 2883 m water depth. These new data, together with the 131 earlier published data sets from the western Pacific, are used for calibrating the SIMMAX-28 transfer function to estimate past sea-surface temperatures. This regional SIMMAX method offers a slightly better understanding of the marginal sea conditions of the South China Sea than the linear transfer function FP-12E, which is based only on open-ocean data. However, both methods are biased toward the tropical temperature regime because of the very limited data from temperate to subpolar regions. The SIMMAX formula was applied to sediment core 17940 from the northeastern South China Sea, with sedimentation rates of 20-80 cm/ka. Results revealed nearly unchanged summer temperatures around 28°C for the last 30 ky, while winter temperatures varied between 19.5°C in the last glacial maximum and 26°C during the Holocene. During Termination 1A, the winter estimates show a Younger Dryas cooling by 3°C subsequent to a temperature optimum of 24°C during the Bölling=Alleröd. Estimates of winter temperature differences between 0 and 100 m water depth document the seasonal variations in the thickness of the mixed layer and provide a new proxy for estimating past changes in the strength of the winter monsoon.

Sea surface and the sea floor regionalization of the Southern Ocean by multivariate cluster analysis, links to ArcGIS project files

This study subdivides the Weddell Sea, Antarctica, into seafloor regions using multivariate statistical methods. These regions are categories used for comparing, contrasting and quantifying biogeochemical processes and biodiversity between ocean regions geographically but also regions under development within the scope of global change. The division obtained is characterized by the dominating components and interpreted in terms of ruling environmental conditions. The analysis uses 28 environmental variables for the sea surface, 25 variables for the seabed and 9 variables for the analysis between surface and bottom variables. The data were taken during the years 1983-2013. Some data were interpolated. The statistical errors of several interpolation methods (e.g. IDW, Indicator, Ordinary and Co-Kriging) with changing settings have been compared for the identification of the most reasonable method. The multivariate mathematical procedures used are regionalized classification via k means cluster analysis, canonical-correlation analysis and multidimensional scaling. Canonical-correlation analysis identifies the influencing factors in the different parts of the cove. Several methods for the identification of the optimum number of clusters have been tested. For the seabed 8 and 12 clusters were identified as reasonable numbers for clustering the Weddell Sea. For the sea surface the numbers 8 and 13 and for the top/bottom analysis 8 and 3 were identified, respectively. Additionally, the results of 20 clusters are presented for the three alternatives offering the first small scale environmental regionalization of the Weddell Sea. Especially the results of 12 clusters identify marine-influenced regions which can be clearly separated from those determined by the geological catchment area and the ones dominated by river discharge.

Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG Last Glacial Maximum

Based on the quantitative study of diatoms and radiolarians, summer sea-surface temperature (SSST) and sea ice distribution were estimated from 122 sediment core localities in the Atlantic, Indian and Pacific sectors of the Southern Ocean to reconstruct the last glacial environment at the EPILOG (19.5-16.0 ka or 23 000-19 000 cal yr. B.P.) time-slice. The statistical methods applied include the Imbrie and Kipp Method, the Modern Analog Technique and the General Additive Model. Summer SSTs reveal greater surface-water cooling than reconstructed by CLIMAP (Geol. Soc. Am. Map Chart. Ser. MC-36 (1981) 1), reaching a maximum (4-5 °C) in the present Subantarctic Zone of the Atlantic and Indian sector. The reconstruction of maximum winter sea ice (WSI) extent is in accordance with CLIMAP, showing an expansion of the WSI field by around 100% compared to the present. Although only limited information is available, the data clearly show that CLIMAP strongly overestimated the glacial summer sea ice extent. As a result of the northward expansion of Antarctic cold waters by 5-10° in latitude and a relatively small displacement of the Subtropical Front, thermal gradients were steepened during the last glacial in the northern zone of the Southern Ocean. Such reconstruction may, however, be inapposite for the Pacific sector. The few data available indicate reduced cooling in the southern Pacific and give suggestion for a non-uniform cooling of the glacial Southern Ocean.

Diatom abundance in surface sediments of the Southern Ocean

Based on the quantitative analysis of diatom assemblages preserved in 274 surface sediment samples recovered in the Pacific, Atlantic and western Indian sectors of the Southern Ocean we have defined a new reference database for quantitative estimation of late-middle Pleistocene Antarctic sea ice fields using the transfer function technique. The Detrended Canonical Analysis (DCA) of the diatom data set points to a unimodal distribution of the diatom assemblages. Canonical Correspondence Analysis (CCA) indicates that winter sea ice (WSI) but also summer sea surface temperature (SSST) represent the most prominent environmental variables that control the spatial species distribution. To test the applicability of transfer functions for sea ice reconstruction in terms of concentration and occurrence probability we applied four different methods, the Imbrie and Kipp Method (IKM), the Modern Analog Technique (MAT), Weighted Averaging (WA), and Weighted Averaging Partial Least Squares (WAPLS), using logarithm-transformed diatom data and satellite-derived (1981-2010) sea ice data as a reference. The best performance for IKM results was obtained using a subset of 172 samples with 28 diatom taxa/taxa groups, quadratic regression and a three-factor model (IKM-D172/28/3q) resulting in root mean square errors of prediction (RMSEP) of 7.27% and 11.4% for WSI and summer sea ice (SSI) concentration, respectively. MAT estimates were calculated with different numbers of analogs (4, 6) using a 274-sample/28-taxa reference data set (MAT-D274/28/4an, -6an) resulting in RMSEP's ranging from 5.52% (4an) to 5.91% (6an) for WSI as well as 8.93% (4an) to 9.05% (6an) for SSI. WA and WAPLS performed less well with the D274 data set, compared to MAT, achieving WSI concentration RMSEP's of 9.91% with WA and 11.29% with WAPLS, recommending the use of IKM and MAT. The application of IKM and MAT to surface sediment data revealed strong relations to the satellite-derived winter and summer sea ice field. Sea ice reconstructions performed on an Atlantic- and a Pacific Southern Ocean sediment core, both documenting sea ice variability over the past 150,000 years (MIS 1 - MIS 6), resulted in similar glacial/interglacial trends of IKM and MAT-based sea-ice estimates. On the average, however, IKM estimates display smaller WSI and slightly higher SSI concentration and probability at lower variability in comparison with MAT. This pattern is a result of different estimation techniques with integration of WSI and SSI signals in one single factor assemblage by applying IKM and selecting specific single samples, thus keeping close to the original diatom database and included variability, by MAT. In contrast to the estimation of WSI, reconstructions of past SSI variability remains weaker. Combined with diatom-based estimates, the abundance and flux pattern of biogenic opal represents an additional indication for the WSI and SSI extent.

Surface Ocean CO2 Atlas (SOCAT) V3

The Surface Ocean CO2 Atlas (SOCAT) is a synthesis of quality-controlled fCO2 (fugacity of carbon dioxide) values for the global surface oceans and coastal seas with regular updates. Version 3 of SOCAT has 14.5 million fCO2 values from 3646 data sets covering the years 1957 to 2014. This latest version has an additional 4.4 million fCO2 values relative to version 2 and extends the record from 2011 to 2014. Version 3 also significantly increases the data availability for 2005 to 2013. SOCAT has an average of approximately 1.2 million surface water fCO2 values per year for the years 2006 to 2012. Quality and documentation of the data has improved. A new feature is the data set quality control (QC) flag of E for data from alternative sensors and platforms. The accuracy of surface water fCO2 has been defined for all data set QC flags. Automated range checking has been carried out for all data sets during their upload into SOCAT. The upgrade of the interactive Data Set Viewer allows better interrogation of the SOCAT data collection and rapid creation of high-quality figures for scientific presentations. Automated data upload has been launched for version 4 and will enable more frequent SOCAT releases in the future. High-profile scientific applications of SOCAT include quantification of the ocean sink for atmospheric carbon dioxide and its long-term variation, detection of ocean acidification, as well as evaluation of coupled-climate and ocean-only biogeochemical models. Users of SOCAT data products are urged to acknowledge the contribution of data providers, as stated in the SOCAT Fair Data Use Statement. This living data publication documents changes in the methods and data sets used in this new version of the SOCAT data collection compared with previous publications of this data collection (Pfeil et al., 2013; Sabine et al., 2013; Bakker et al., 2014).

(Table 1) Age control points and dating methods for sediment cores of the Sea of Okhotsk

Marine- and terrestrial-derived biomarkers (alkenones, brassicasterol, dinosterol, and long-chain n-alkanes), as well as carbonate, biogenic opal, and ice-rafted debris (IRD), were measured in two sediment cores in the Sea of Okhotsk, which is located in the northwestern Pacific rim and characterized by high primary productivity. Down-core profiles of phytoplankton markers suggest that primary productivity abruptly increased during the global Meltwater Pulse events 1A (about 14 ka) and 1B (about 11 ka) and stayed high in the Holocene. Spatial and temporal distributions of the phytoplankton productivity were found to be consistent with changes in the reconstructed sea ice distribution on the basis of the IRD. This demonstrates that the progress and retreat of sea ice regulated primary productivity in the Sea of Okhotsk with minimum productivity during the glacial period. The mass accumulation rates of alkenones, CaCO3, and biogenic opal indicate that the dominant phytoplankton species during deglaciation was the coccolithophorid, Emiliania huxleyi, which was replaced by diatoms in the late Holocene. Such a phytoplankton succession was probably caused by an increase in silicate supply to the euphotic layer, possibly associated with a change in surface hydrography and/or linked to enhanced upwelling of North Pacific Deep Water.

Dissolved Iron, Aluminium, and Dissolved Inorganic Phosphorus in the (Sub-)Tropical Atlantic Surface Ocean

Inorganic nitrogen depletion restricts productivity in much of the low-latitude oceans, generating a selective advantage for diazotrophic organisms capable of fixing atmospheric dinitrogen (N2). However, the abundance and activity of diazotrophs can in turn be controlled by the availability of other potentially limiting nutrients, including phosphorus (P) and iron (Fe). Here we present high-resolution data (~0.3°) for dissolved iron, aluminum, and inorganic phosphorus that confirm the existence of a sharp north-south biogeochemical boundary in the surface nutrient concentrations of the (sub)tropical Atlantic Ocean. Combining satellite-based precipitation data with results from a previous study, we here demonstrate that wet deposition in the region of the intertropical convergence zone acts as the major dissolved iron source to surface waters. Moreover, corresponding observations of N2 fixation and the distribution of diazotrophic Trichodesmium spp. indicate that movement in the region of elevated dissolved iron as a result of the seasonal migration of the intertropical convergence zone drives a shift in the latitudinal distribution of diazotrophy and corresponding dissolved inorganic phosphorus depletion. These conclusions are consistent with the results of an idealized numerical model of the system. The boundary between the distinct biogeochemical systems of the (sub)tropical Atlantic thus appears to be defined by the diazotrophic response to spatial-temporal variability in external Fe inputs. Consequently, in addition to demonstrating a unique seasonal cycle forced by atmospheric nutrient inputs, we suggest that the underlying biogeochemical mechanisms would likely characterize the response of oligotrophic systems to altered environmental forcing over longer timescales.