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.

Calcium carbonate content and anisotropy of magnetic susceptibility of ODP Leg 172 sites

Sediments from Sites 1057 and 1061 of Ocean Drilling Program Leg 172 on the Blake Outer Ridge exhibit nearly isotropic magnetic susceptibility. Resolving the degree of anisotropy of magnetic susceptibility proved difficult in many samples because of the generally weak magnetic susceptibility of the sediments relative to the noise level of the susceptibility meters used. Lineation varies from 1.0 to 1.013 and foliation varies from 1.0 to 1.08 in the samples that pass rejection criteria. In general the foliation is better resolved than the lineation, particularly at Site 1061, where the foliation exhibits long-term trends that mimic the mean susceptibility. The changes in the foliation at this site are likely the result of changes in the magnetic mineralogy of the sediment. The poorly developed or absent magnetic fabric in the sediments overall can be attributed to high carbonate concentrations and to a circulation regime that was diffuse or with currents too weak to effectively align magnetic particles.

Pliocene-Pleistocene stable isotope record of ODP Site 107-653

Planktonic foraminiferal oxygen and carbon isotope analyses from Tyrrhenian Sea Ocean Drilling Program (ODP) Site 653 provide a continuous record of the Pliocene-Pleistocene paleoceanographic history of the Mediterranean. Long-term trends in oxygen isotopes primarily reflect changes in global climatic conditions, with a more local or regional signal superimposed on this record. For example, significant enrichments in 18O due to decreases in surface water temperature and/or increases in continental ice volume occurred at 3.1, 2.7, 2.1, 1.6, and 0.4 Ma. In contrast to most open-ocean results, the early Pliocene 6lsO record of Site 653 exhibits high-amplitude fluctuations indicative of very unstable climatic conditions in this region. Another unique aspect of this Mediterranean d18Orecord is the pronounced cooling at the Pliocene/Pleistocene boundary. The carbon isotope record for Site 653 also exhibits high-amplitude variability throughout the Pliocene-Pleistocene. This variability most probably reflects changes in the carbon isotopic composition of the source of Mediterranean surface waters.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

We present new d13C measurements of atmospheric CO2 covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in d13C(atm) of 0.5 permil occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in d13C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO2]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of d13C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

Chemistry and grain-size distribution of terrigenous sediments deposited at Site GeoB9501

The Sahara Desert is the largest source of mineral dust in the world. Emissions of African dust increased sharply in the early 1970s, a change that has been attributed mainly to drought in the Sahara/Sahel region caused by changes in the global distribution of sea surface temperature. The human contribution to land degradation and dust mobilization in this region remains poorly understood, owing to the paucity of data that would allow the identification of long-term trends in desertification. Direct measurements of airborne African dust concentrations only became available in the mid-1960s from a station on Barbados and subsequently from satellite imagery since the late 1970s: they do not cover the onset of commercial agriculture in the Sahel region ~170 years ago. Here we construct a 3,200-year record of dust deposition off northwest Africa by investigating the chemistry and grain-size distribution of terrigenous sediments deposited at a marine site located directly under the West African dust plume. With the help of our dust record and a proxy record for West African precipitation we find that, on the century scale, dust deposition is related to precipitation in tropical West Africa until the seventeenth century. At the beginning of the nineteenth century, a sharp increase in dust deposition parallels the advent of commercial agriculture in the Sahel region. Our findings suggest that human-induced dust emissions from the Sahel region have contributed to the atmospheric dust load for about 200 years.

(Table 1) Active-layer thickness as measured by the Biocomplexity Experiment (BE) for the flooded, drained and control sections in Barrow, Alaska

A continuous time series of annual soil thaw records, extending from 1994 to 2009, is available for comparison with the records of thaw obtained from the Biocomplexity Experiment (BE) for the period 2006-2009. Discontinuous records of thaw at Barrow from wet tundra sites date back to the 1960s. Comparisons between the longer records with the BE observations reveal strong similarities. Records of permafrost temperature, reflecting changes in the annual surface energy exchange, are available from the 1950s for comparison with results from measurement programs begun in 2002. The long-term systematic geocryological investigations at Barrow indicate an increase in permafrost temperature, especially during the last several years. The increase in near-surface permafrost temperature is most pronounced in winter. Marked trends are not apparent in the active-layer record, although subsidence measurements on the North Slope indicate that penetration into the ice-rich layer at the top of permafrost has occurred over the past decade. Active-layer thickness values from the 1960s are generally higher than those from the 1990s, and are very similar to those of the 2000s. Analysis of spatial active-layer observations at representative locations demonstrates significant variations in active-layer thickness between different landscape types, reflecting the influence of vegetation, substrate, microtopography, and, especially, soil moisture. Landscape-specific differences exist in the response of active-layer thickness to climatic forcing. These differences are attributable to the existence of localized controls related to combinations of surface and subsurface characteristics. The geocryological records at Barrow illustrate the importance and effectiveness of sustained, well organized monitoring efforts to document long-term trends.

Radiocarbon, stable isotopic compositions and seawater chemical properties of Black Sea core-top sediments and CTD stations

Chloropigments and their derivative pheopigments preserved in sediments can directly be linked to photosynthesis. Their carbon and nitrogen stable isotopic compositions have been shown to be a good recorder of recent and past surface ocean environmental conditions tracing the carbon and nitrogen sources and dominant assimilation processes of the phytoplanktonic community. In this study we report results from combined compound-specific radiocarbon and stable carbon and nitrogen isotope analysis to examine the time-scales of synthesis and fate of chlorophyll-a and its degradation products pheophytin-a, pyropheophytin-a, and 132,173-cyclopheophorbide-a-enol until burial in Black Sea core-top sediments. The pigments are mainly of marine phytoplanktonic origin as implied by their stable isotopic compositions. Pigment ?15N values indicate nitrate as the major uptake substrate but 15N-depletion towards the open marine setting indicates either contribution from N2-fixation or direct uptake of ammonium from deeper waters. Radiocarbon concentrations translate into minimum and maximum pigment ages of approximately 40 to 1200 years. This implies that protective mechanisms against decomposition such as association with minerals, storage in deltaic anoxic environments, or eutrophication-induced hypoxia and light limitation are much more efficient than previously thought. Moreover, seasonal variations of nutrient source, growth period, and habitat and their associated isotopic variability are likely at least as strong as long-term trends. Combined triple isotope analysis of sedimentary chlorophyll and its primary derivatives is a powerful tool to delineate biogeochemical and diagenetic processes in the surface water and sediments, and to assess their precise time-scales.

(Table 1) Annual mean concentrations of atmospheric a- and g-Hexachlorocyclohexane at stations Alert, Canada and Zeppelin, Norway

Twelve year datasets of weekly atmospheric concentrations of alpha- and gamma-HCH were compared between the two Arctic monitoring stations of Alert, Nunavut, Canada, and Zeppelin Mountain, Svalbard, Norway. Time-series analysis was conducted with the use of dynamic harmonic regression (DHR), which provided a very good model fit, to examine both the seasonal behaviour in these isomers and the longer-term, underlying trends. Strong spatial differences were not apparent between the two sites, although subtle differences in seasonal behaviour and composition were identified. For example, the composition of gamma-HCH to total HCH (alpha + gamma) was greater at Zeppelin compared to Alert, probably reflecting this site's proximity to major use regions of lindane. Pronounced seasonality in air concentrations for gamma-HCH was marked by a 'spring maximum event' (SME), confirming earlier studies. For alpha-HCH, the SME was much weaker and only evident at Alert, whereas at Zeppelin, seasonal fluctuations for alpha-HCH were marked by elevated concentrations in summer and lower concentrations during winter, with this pattern most apparent for the years after 2000. We attribute this difference in spatial and temporal patterns to the Arctic oscillation. A similar climatic pattern was not evident at either site in the gamma-HCH data. Seasonally adjusted, long-term trends revealed declining concentrations at both sites for alpha- and gamma-HCH over the entire time-series. Recent legislation affecting lindane use appear to account for this decline in gamma-HCH, with little evidence of a delay or 'lag' between the banning of lindane in Europe (a main source region) or Canada, and a decline in air concentrations observed at both Arctic sites.

Seawater carbonate chemistry and calcification rate of eastern oyster Crassostrea virginica, 2011

Anthropogenic carbon dioxide (CO2) emissions reduce pH of marine waters due to the absorption of atmospheric CO2 and formation of carbonic acid. Estuarine waters are more susceptible to acidification because they are subject to multiple acid sources and are less buffered than marine waters. Consequently, estuarine shell forming species may experience acidification sooner than marine species although the tolerance of estuarine calcifiers to pH changes is poorly understood. We analyzed 23 years of Chesapeake Bay water quality monitoring data and found that daytime average pH significantly decreased across polyhaline waters although pH has not significantly changed across mesohaline waters. In some tributaries that once supported large oyster populations, pH is increasing. Current average conditions within some tributaries however correspond to values that we found in laboratory studies to reduce oyster biocalcification rates or resulted in net shell dissolution. Calcification rates of juvenile eastern oysters, Crassostrea virginica, were measured in laboratory studies in a three-way factorial design with 3 pH levels, two salinities, and two temperatures. Biocalcification declined significantly with a reduction of ~0.5 pH units and higher temperature and salinity mitigated the decrease in biocalcification.

(Table 2) Lipid fluxes derived from sediment trap samples

Phytoplankton is a sentinel of marine ecosystem change. Composed by many species with different life-history strategies, it rapidly responds to environment changes. An analysis of the abundance of 54 phytoplankton species in Galicia (NW Spain) between 1989 and 2008 to determine the main components of temporal variability in relation to climate and upwelling showed that most of this variability was stochastic, as seasonality and long term trends contributed to relatively small fractions of the series. In general, trends appeared as non linear, and species clustered in 4 groups according to the trend pattern but there was no defined pattern for diatoms, dinoflagellates or other groups. While, in general, total abundance increased, no clear trend was found for 23 species, 14 species decreased, 4 species increased during the early 1990s, and only 13 species showed a general increase through the series. In contrast, series of local environmental conditions (temperature, stratification, nutrients) and climate-related variables (atmospheric pressure indices, upwelling winds) showed a high fraction of their variability in deterministic seasonality and trends. As a result, each species responded independently to environmental and climate variability, measured by generalized additive models. Most species showed a positive relationship with nutrient concentrations but only a few showed a direct relationship with stratification and upwelling. Climate variables had only measurable effects on some species but no common response emerged. Because its adaptation to frequent disturbances, phytoplankton communities in upwelling ecosystems appear less sensitive to changes in regional climate than other communities characterized by short and well defined productive periods.

Clay mineralogy of ODP Leg 178 sites

The clay mineral compositions of upper Miocene to Quaternary sediments recovered at Ocean Drilling Program (ODP) Leg 178, Sites 1095 and 1096, from the continental rise west of the Antarctic Peninsula were analyzed in order to reconstruct the Neogene and Quaternary Antarctic paleoclimate and ice dynamics. The clay mineral assemblages are dominated by smectite, illite, and chlorite. Kaolinite occurs only in trace amounts. Analysis of a surface-sample data set facilitates the assignment of these clay minerals to particular source areas on the Antarctic Peninsula and, thus, the reconstruction of transport pathways. In the ODP cores, clay mineral composition cyclically alternates between two end-member assemblages. One assemblage is characterized by <20% smectite and >40% chlorite. The other assemblage has >20% smectite and <40% chlorite. Illite fluctuates between 30% and 50% without a significant affinity to one end-member assemblage. By comparison with a Quaternary sediment sequence from gravity core PS1565, the clay mineral fluctuations can be ascribed to glacial and interglacial periods, respectively. The cyclic changes in the clay mineral composition suggest that glacial-interglacial cycles, repeated ice advances and retreats, and changes in the Antarctic ice volume were already a main control of the global climate in late Miocene time. Throughout the late Neogene and Quaternary, the clay mineral records in the drift sediments exhibit only slight long-term changes predominantly attributed to local changes in glacial erosion and supply of source rocks. The absence of clear long-term trends associated with major climatic or glaciological changes points to an onset of vast glaciation in the Antarctic Peninsula region before ~9 Ma and to relative stability of the Antarctic ice sheet since then.