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.

Chronology for the EDML ice core from Dronning Maud Land, Antarctica, over the last 150 000 years

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDML1 is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDML1 was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core.

Accumulation rates for western Dronning Maud Land from firn cores and snow pits

As a result of intensive field activities carried out by several nations over the past 15 years, a set of accumulation measurements for western Dronning Maud Land, Antarctica, was collected, based on firn-core drilling and snow-pit sampling. This new information was supplemented by earlier data taken from the literature, resulting in 111 accumulation values. Using Geographical Information Systems software, a first region-wide mean annual snow-accumulation field was derived. In order to define suitable interpolation criteria, the accumulation records were analyzed with respect to their spatial autocorrelation and statistical properties. The resulting accumulation pattern resembles well known characteristics such as a relatively wet coastal area with a sharp transition to the dry interior, but also reveals complex topographic effects. Furthermore, this work identifies new high-return shallow drilling sites by uncovering areas of insufficient sampling density.

Concentrations of PBDE and BDE congeners in adipose tissue of East Greenland polar bears (Ursus maritimus)

Polybrominated diphenyl ethers (PBDEs) were analysed in adipose tissue from 92 East Greenland polar bears (Ursus maritimus) sampled during 1999-2001. Mean SumPBDE concentrations were 70 ng/g lipid weight (lw) (range: 22-192 ng/g lw) and showed no relationship with age or sex. Of the 32 analysed PBDE congeners; BDE47, BDE153, BDE99 and BDE100 dominated, and comprised 99.6% of the SumPBDE concentration. The SumPBDE concentration had a highly significant correlation with SumPCB, SumCHL, dieldrin, HCB and SumHCH concentrations. We found a seasonal pattern of median SumPBDE concentration with 1.2 to 1.8 times higher concentrations in March to July than the rest of the year. The seasonal variation also provides a clue to the seasonal exposure, bio-availability, toxic exposure and degradation. We suggest that future geographical PBDE data comparisons may not need corrections for sex or age, but such data sets should be corrected for seasonal variability, using the presented correctional trigonometric regression.

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').

Contusotruncana contusa abundances and morphotypes in the cretaceous mid-latitude South Atlantic

Spatial and temporal patterns in test size and shape (test conicity and spiral roundness) and absolute abundance (accumulation rate) of the planktonic foraminifer Contusotruncana contusa were studied in the South Atlantic Ocean (DSDP sites 356, 516, 525 and 527) during an interval corresponding to the last 800 kyr of the Cretaceous. The variation in absolute abundance of C. contusa was characterised by alternating periods of high and low abundance; some of these periods were traceable across the entire mid-latitude South Atlantic Ocean. While the mean spiral roundness did not show any interpretable patterns, a sudden increase of the mean test size and mean test conicity occurred between 65.3 and 65.2 Ma (based on linear interpolation within the Cretaceous part of Subchron C29R) at all sites studied, indicating a poleward migration followed by rapid withdrawal of the low-latitude C. contusa morphotypes from the mid-latitude South Atlantic Ocean. We suggest that this event was caused by a short period of surface-water warming in the southern mid-latitudes corresponding to the brief high-latitude warming event and associated faunal migrations in the Boreal and Austral realms.

Geochemistry of DSDP Hole 91-596

Geochemical analyses of sediments from the top 24.5 m of Deep Sea Drilling Project hole 596 (23°51.20'S, 169°39.27'W) show great variability in the composition of pelagic clays accumulated in the South Pacific since the late Cretaceous. Elemental associations indicate that most of this variability can be attributed to variations in abundances of six sediment end-member components: detrital (eolian), andesitic (volcanic), hydrothermal, hydrogenous, phosphate (fish debris), and biogenic silica. We develop a sedimentation model which is used to infer processes that might have influenced the accumulation rates of these components over the last 85 million years. The accumulation of eolian detritus in the South Pacific shows some similarities to that observed in the North Pacific and has been largely controlled by global climate trends in the Cenozoic. Much of the variation in the accumulation of other sediment components likely reflects the paleoceanographic evolution of the South Pacific. The most notable change in the sedimentary environment occurred at about the Paleogene/Neogene boundary. At that time, significant changes in the color, mineralogy, and chemistry of the sediment probably reflect major shifts in climate mode as well as oceanic circulation in the central South Pacific region.

Revised geomagnetic polarity timescale for the late Oligocene to early Miocene of ODP Site 177-1090

At Ocean Drilling Program (ODP) Site 1090 (subantarctic South Atlantic), benthic foraminiferal stable isotope data (from Cibicidoides and Oridorsalis) span the late Oligocene through early Miocene (~24-16 Ma) at a temporal resolution of ~5 ky. Over the same interval, a magnetic polarity stratigraphy can be unequivocally correlated to the geomagnetic polarity time scale (GPTS), thereby providing direct correlation of the isotope record to the GPTS. In an initial age model, we use the newly derived age of the Oligocene/Miocene (O/M) boundary of 23.0 Ma of Shackleton et al. (2000, doi:10.1130/0091-7613(2000)28<447:ACAFTO>2.0.CO;2), revised to the new astronomical calculation (La2003) of Laskar et al (2004, doi:10.1016/j.icarus.2004.04.005) to recalculate the spline ages of Cande and Kent (1995, doi:10.1029/94JB03098). We then tune the Site 1090 dekta18O record to obliquity using La2003. In this manner, we are able to refine the ages of polarity chrons C7n through C5Cn.1n. The new age model is consistent, within one obliquity cycle, with previously tuned ages for polarity chrons C7n through C6Bn from Shackleton et al. (2000) when rescaled to La2003. The results from Site 1090 provide independent evidence for the revised age of the Oligocene/Miocene boundary of 23.0 Ma. For early Miocene polarity chrons C6AAr through C5Cn, our obliquity-scale age model is the first to allow a direct calibration to the GPTS. The new ages are generally within one obliquity cycle of those obtained by rescaling the Cande and Kent (1995) interpolation using the new age of the O/M boundary (23.0 Ma) and the same middle Miocene control point (14.8 Ma) used by Cande and Kent (1995).

Physical oceanography of CTD stations in the Persian Gulf

The following tables show physical and chemical data observed by the "Meteor" in the Persian Gulf and the Strait of Hormus. This study was performed in accordance with the general programme of the International Indian Ocean Expedition (IIOE) during the oeriod from March 25th until April 16th, 1965. The water temperature was measured by reversing thermometers; in most cases two instruments were used simultaneously. The absolute mean temperature difference of this double measurement is 0.0153 °C. The salinity was determined both by salinometer and by titration. In this case of the density, the specific volume anomaly, the dynamic depth anomaly, the sound velocity and the interpolation for standard depths were carried out by the National Oceanographic Data Center (NODC), Washington.