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.

Sampling intervals, fluxes, percentages of total flux, organic carbon and lithogen for sediment trap CB9

We combined the analysis of sediment trap data and satellite-derived sea surface chlorophyll to quantify the amount of organic carbon export to the deep sea in the upwelling induced high production area off northwest Africa. In contrast to the generally global or basin-wide adoption of export models, we used a regionally fitted empirical model. Furthermore, the application of our model was restricted to a dynamically defined region of high chlorophyll concentration in order to restrict the model application to an environment of more homogeneous export processes. We developed a correlation-based approximation to estimate the surface source area for a sediment trap deployed from 11 June 1998 to 7 November 1999 at 21.25°N latitude and 20.64°W longitude off Cape Blanc. We also developed a regression model of chlorophyll and export of organic carbon to the 1000 m depth level. Carbon export was calculated for an area of high chlorophyll concentration (>1 mg/m**3) adjacent to the coast on a daily basis. The resulting zone of high chlorophyll concentration was 20,000-800,000 km**2 large and yielded a yearly export of 1.123 to 2.620 Tg organic carbon. The average organic carbon export within the area of high chlorophyll concentration was 20.6 mg/m**2d comparable to 13.3 mg/m**2d as found in the sediment trap results if normalized to the 1000 m level. We found strong interannual variability in export. The period autumn 1998 to summer 1999 was exceeding the mean of the other three comparable periods by a factor of 2.25. We believe that this approach of using more regionally fitted models can be successfully transferred even to different oceanographic regions by selecting appropriate definition criteria like chlorophyll concentration for the definition of an area to which it is applicable.

Global database of surface ocean particulate organic carbon export fluxes diagnosed from the 234Th technique

Through the processes of the biological pump, carbon is exported to the deep ocean in the form of dissolved and particulate organic matter. There are several ways by which downward export fluxes can be estimated. The great attraction of the 234Th technique is that its fundamental operation allows a downward flux rate to be determined from a single water column profile of thorium coupled to an estimate of POC/234Th ratio in sinking matter. We present a database of 723 estimates of organic carbon export from the surface ocean derived from the 234Th technique. Data were collected from tables in papers published between 1985 and 2013 only. We also present sampling dates, publication dates and sampling areas. Most of the open ocean Longhurst provinces are represented by several measurements. However, the Western Pacific, the Atlantic Arctic, South Pacific and the South Indian Ocean are not well represented. There is a variety of integration depths ranging from surface to 220m. Globally the fluxes ranged from -22 to 125 mmol of C/m**2/d. We believe that this database is important for providing new global estimate of the magnitude of the biological carbon pump.

Alkenone stable carbon isotopes, planktic foraminfera stable carbon isotopes, alkenone unsaturation indices and calculated sea surface temperatures and atmospheric pCO2 for ODP Site 165-999

Temperature reconstructions indicate that the Pliocene was ~3 °C warmer globally than today, and several recent reconstructions of Pliocene atmospheric CO2 indicate that it was above pre-industrial levels and similar to those likely to be seen this century. However, many of these reconstructions have been of relatively low temporal resolution, meaning that these records may have failed to capture variations associated with the 41 Kyr glacial-interglacial cycles thought to operate in the Pliocene. Here we present a new, high temporal resolution alkenone carbon isotope based record of pCO2 spanning 2.8 to 3.3 million years ago from ODP Site 999. Our record is of high enough resolution (~19 Kyrs) to resolve glacial-interglacial changes beyond the intrinsic uncertainty of the proxy method. The record suggests that Pliocene CO2 levels were relatively stable, exhibiting variation less than 55 ppm. We perform sensitivity studies to investigate the possible effect of changing sea surface temperature, which highlights the importance of accurate and precise SST reconstructions for alkenone palaeobarometry, but demonstrate that these uncertainties do not affect our conclusions of relatively stable pCO2 levels during this interval.

Clay mineralogy and organic carbon content of sediment core GeoB1401-4

Late Quaternary sediments recovered in a core from the area of the Zaire Fan, Central Africa, were analyzed for clay mineral composition in order to reconstruct fluctuations in the sediment input and freshwater discharge of the Zaire River. Clay mineral assemblages are dominated by kaolinite and smectite, which both originate mainly from the Zaire River and contain only minor contributions of eolian dust. Smectite crystallinity and chemical character of illites (Fe, Mg- or Al-rich) are used to track sediment input from the Zaire River and assess fluctuations in the freshwater discharge. Both parameters record a high-latitude forcing of river runoff at 100 ka periodicities reflecting glacial aridity and increased runoff during interglacials 1, 5 and 7. This signal is also observed in kaolinite/smectite ratios which represent the extension and intensity of the freshwater plume of the Zaire River. Clay mineral proxies reveal that river discharge and associated sediment input fluctuated in tune with precessional cycles of African monsoon intensity. Increased eolian input of kaolinite-rich dust with intensified northeast trades during glacials flattens the precessional signal in kaolinite/smectite ratios.

Total organic carbon content and redox-sensitiv trace metal concentrations of selected late Cenoman claystones from DSDP Hole 93-603B

The Cenomanian/Turonian (C/T) intervals at DSDP Sites 105 and 603B from the northern part of the proto-North Atlantic show high amplitude, short-term cyclic variations in total organic carbon (TOC) content. The more pronounced changes in TOC are also reflected by changes in lithology from green claystones (TOC<1%) to black claystones (TOC>1%). Although their depositional history was different, the individual TOC cycles at Sites 105 and 603B can be correlated using stable carbon isotope stratigraphy. Sedimentation rates obtained from the isotope stratigraphy and spectral analyses indicate that these cycles were predominately precession controlled. The coinciding variations in HI, OI, delta13Corg and the abundance of marine relative to terrestrial biomarkers, as well as the low abundance of lignin pyrolysis products generated from the kerogen of the black claystones, indicate that these cyclic variations reflect changes in the contribution of marine organic matter (OM). The cooccurrence of lamination, enrichment of redox-sensitive trace metals and presence of molecular fossils of pigments from green sulfur bacteria indicate that the northern proto-North Atlantic Ocean water column was periodically euxinic from the bottom to at least the base of the photic zone (<150 m) during the deposition of the black claystones. In contrast, the green claystones are bioturbated, are enriched in Mn, do not show enrichments in redox-sensitive trace metals and show biomarker distributions indicative of long oxygen exposure times, indicating more oxic water conditions. At the same time, there is evidence (e.g., abundance of biogenic silica and significant 13C-enrichment for OC of phytoplanktic origin) for enhanced primary productivity during the deposition of the black claystones. We propose that increased primary productivity periodically overwhelmed the oxic OM remineralisation potential of the bottom waters resulting in the deposition of OM-rich black claystones. Because the amount of oxygen used for OM remineralisation exceeded the amount supplied by diffusion and deep-water circulation, the northern proto-North Atlantic became euxinic during these periods. Both Sites 105 and 603B show trends of continually increasing TOC contents and HI values of the black claystones up section, which most likely resulted from both enhanced preservation due to increased anoxia and increased production of marine OM during oceanic anoxic event 2 (OAE2).

Age determination and stable carbon isotope ratios of Neogloboquadrina pachyderma of IODP Hole 302-M0004C

The Integrated Ocean Drilling Program (IODP) Arctic Coring Expedition (ACEX) Hole 4C from the Lomonosov Ridge in the central Arctic Ocean recovered a continuous 18 m record of Quaternary foraminifera yielding evidence for seasonally ice-free interglacials during the Matuyama, progressive development of large glacials during the mid-Pleistocene transition (MPT) ~1.2-0.9 Ma, and the onset of high-amplitude 100-ka orbital cycles ~500 ka. Foraminiferal preservation in sediments from the Arctic is influenced by primary (sea ice, organic input, and other environmental conditions) and secondary factors (syndepositional, long-term pore water dissolution). Taking these into account, the ACEX 4C record shows distinct maxima in agglutinated foraminiferal abundance corresponding to several interglacials and deglacials between marine isotope stages (MIS) 13-37, and although less precise dating is available for older sediments, these trends appear to continue through the Matuyama. The MPT is characterized by nearly barren intervals during major glacials (MIS 12, 16, and 22-24) and faunal turnover (MIS 12-24). Abundant calcareous planktonic (mainly Neogloboquadrina pachyderma sin.) and benthic foraminifers occur mainly in interglacial intervals during the Brunhes and very rarely in the Matuyama. A distinct faunal transition from calcareous to agglutinated foraminifers 200-300 ka in ACEX 4C is comparable to that found in Arctic sediments from the Lomonosov, Alpha, and Northwind ridges and the Morris Jesup Rise. Down-core disappearance of calcareous taxa is probably related to either reduced sea ice cover prior to the last few 100-ka cycles, pore water dissolution, or both.

On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century, links to supplementary material

We use a suite of eight ocean biogeochemical/ecological general circulation models from the MAREMIP and CMIP5 archives to explore the relative roles of changes in winds (positive trend of Southern Annular Mode, SAM) and in warming- and freshening-driven trends of upper ocean stratification in altering export production and CO2 uptake in the Southern Ocean at the end of the 21st century. The investigated models simulate a broad range of responses to climate change, with no agreement ona dominance of either the SAM or the warming signal south of 44° S. In the southernmost zone, i.e., south of 58° S, they concur on an increase of biological export production, while between 44 and 58° S the models lack consensus on the sign of change in export. Yet, in both regions, the models show an enhanced CO2 uptake during spring and summer. This is due to a larger CO 2 (aq) drawdown by the same amount of summer export production at a higher Revelle factor at the end of the 21st century. This strongly increases the importance of the biological carbon pump in the entire Southern Ocean. In the temperate zone, between 30 and 44° S all models show a predominance of the warming signal and a nutrient-driven reduction of export production. As a consequence, the share of the regions south of 44° S to the total uptake of the Southern Ocean south of 30° S is projected to increase at the end of the 21st century from 47 to 66% with a commensurable decrease to the north. Despite this major reorganization of the meridional distribution of the major regions of uptake, the total uptake increases largely in line with the rising atmospheric CO2. Simulations with the MITgcm-REcoM2 model show that this is mostly driven by the strong increase of atmospheric CO2, with the climate-driven changes of natural CO2 exchange offsetting that trend only to a limited degree (~10%) and with negligible impact of climate effects on anthropogenic CO2 uptake when integrated over a full annual cycle south of 30° S.

Adsorption of Lennard-Jones fluids in carbon slit pores of a finite length. A computer simulation study

The adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length was studied with Canonical Ensemble (NVT) and Gibbs Ensemble Monte Carlo Simulations (GEMC). The Canonical Ensemble was a collection of cubic simulation boxes in which a finite pore resides, while the Gibbs Ensemble was that of the pore space of the finite pore. Argon was used as a model for Lennard-Jones fluids, while the adsorbent was modelled as a finite carbon slit pore whose two walls were composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. The Lennard-Jones (LJ) 12-6 potential model was used to compute the interaction energy between two fluid particles, and also between a fluid particle and a carbon atom. Argon adsorption isotherms were obtained at 87.3 K for pore widths of 1.0, 1.5 and 2.0 nm using both Canonical and Gibbs Ensembles. These results were compared with isotherms obtained with corresponding infinite pores using Grand Canonical Ensembles. The effects of the number of cycles necessary to reach equilibrium, the initial allocation of particles, the displacement step and the simulation box size were particularly investigated in the Monte Carlo simulation with Canonical Ensembles. Of these parameters, the displacement step had the most significant effect on the performance of the Monte Carlo simulation. The simulation box size was also important, especially at low pressures at which the size must be sufficiently large to have a statistically acceptable number of particles in the bulk phase. Finally, it was found that the Canonical Ensemble and the Gibbs Ensemble both yielded the same isotherm (within statistical error); however, the computation time for GEMC was shorter than that for canonical ensemble simulation. However, the latter method described the proper interface between the reservoir and the adsorbed phase (and hence the meniscus).

Flux of organic carbon in a riverine mangrove wetland in the Florida Coastal Everglades

Short-term (daily) and seasonal variations in concentration and flux of dissolved organic carbon (DOC) were examined over 15 tidal cycles in a riverine mangrove wetland along Shark River, Florida in 2003. Due to the influence of seasonal rainfall and wind patterns on Shark River’s hydrology, samplings were made to include wet, dry and transitional (Norte) seasons. We used a flume extending from a tidal creek to a basin forest to measure vertical (vegetated soil/water column) and horizontal (mangrove forest/tidal creek) flux of DOC. We found significant (p < 0.05) variations in surface water temperature, salinity, conductivity, pH and mean concentration of DOC with season. Water temperature and salinity followed seasonal patterns of air temperature and rainfall, while mean DOC concentration was highest during the dry season (May), followed by the wet (October) and ‘Norte’ (December) seasons. This pattern of DOC concentration may be due to a combination of litter production and inundation pattern of the wetland. In contrast to daily (between tides) variation in DOC flux between the mangrove forest and tidal creek, daily variations of mean water quality were not significant. However, within-tide variation of DOC flux, dissolved oxygen content and salinity was observed. This indicated that the length of inundation and water source (freshwater vs. saltwater) variation across tidal cycles influenced water quality and DOC flux in the water column. Net DOC export was measured in October and December, suggesting the mangrove forest was a source of DOC to the adjacent tidal creek during these periods. Net annual export of DOC from the fringe mangrove to both the tidal creek and basin mangrove forest was 56 g C m−2 year−1. The seasonal pattern in our flux results indicates that DOC flux from this mangrove forest may be governed by both freshwater discharge and tidal range.

Enhanced 3-dimensional carbon nanotube based anodes for lithium-ion battery applications

A prototype 3-dimensional (3D) anode, based on multiwall carbon nanotubes (MWCNTs), for Li-ion batteries (LIBs), with potential use in Electric Vehicles (EVs) was investigated. The unique 3D design of the anode allowed much higher areal mass density of MWCNTs as active materials, resulting in more amount of Li+ ion intake, compared to that of a conventional 2D counterpart. Furthermore, 3D amorphous Si/MWCNTs hybrid structure offered enhancement in electrochemical response (specific capacity 549 mAhg–1 ). Also, an anode stack was fabricated to further increase the areal or volumetric mass density of MWCNTs. An areal mass density of the anode stack 34.9 mg/cm2 was attained, which is 1,342% higher than the value for a single layer 2.6 mg/cm2. Furthermore, the binder-assisted and hot-pressed anode stack yielded the average reversible, stable gravimetric and volumetric specific capacities of 213 mAhg–1 and 265 mAh/cm3, respectively (at 0.5C). Moreover, a large-scale patterned novel flexible 3D MWCNTs-graphene-polyethylene terephthalate (PET) anode structure was prepared. It generated a reversible specific capacity of 153 mAhg–1 at 0.17C and cycling stability of 130 mAhg –1 up to 50 cycles at 1.7C.

Radiogenic neodymium, strontium and lead isotopes of authigenic and detrital sedimentary phases and stable oxygen and carbon isotopes of benthic foraminifera from gravity cores SO213-59-2 and SO213-60-1 (central South Pacific).

The South Pacific is a sensitive location for the variability of the global oceanic thermohaline circulation given that deep waters from the Atlantic Ocean, the Southern Ocean, and the Pacific Basin are exchanged. Here we reconstruct the deep water circulation of the central South Pacific for the last two glacial cycles (from 240,000 years ago to the Holocene) based on radiogenic neodymium (Nd) and lead (Pb) isotope records complemented by benthic stable carbon data obtained from two sediment cores located on the flanks of the East Pacific Rise. The records show small but consistent glacial/interglacial changes in all three isotopic systems with interglacial average values of -5.8 and 18.757 for epsilon Nd and 206Pb/204Pb, respectively, whereas glacial averages are -5.3 and 18.744. Comparison of this variability of Circumpolar Deep Water (CDW) to previously published records along the pathway of the global thermohaline circulation is consistent with reduced admixture of North Atlantic Deep Water to CDW during cold stages. The absolute values and amplitudes of the benthic delta13C variations are essentially indistinguishable from other records of the Southern Hemisphere and confirm that the low central South Pacific sedimentation rates did not result in a significant reduction of the amplitude of any of the measured proxies. In addition, the combined detrital Nd and strontium (87Sr/86Sr) isotope signatures imply that Australian and New Zealand dust has remained the principal contributor of lithogenic material to the central South Pacific.

Stable carbon and oxygen isotope ratios of planktonic foraminifera from the equatorial Atlantic

Carbon isotopic records of nutrient-depleted surface water place constraints on the past fertility of the oceans and on past atmospheric pCO2 levels. The best records of nutrient-depleted delta13C are obtained from planktonic foraminifera living in the thick mixed layers of the western equatorial and tropical Atlantic Ocean. We have produced a composite, stacked Globigerinoides sacculifer delta13C record from the equatorial Atlantic, which exhibits significant spectral power at the 100,000- and 41,000-year Milankovitch periods, but no power at the 23,000-year period. Similar to the record presented by Shackleton and Pisias [1985], surface-deep ocean Delta delta13C produced with the G. sacculifer record leads the delta18O ice volume record. However, the glacial-interglacial amplitudes of Delta delta13C differ between our record and Shackleton and Pisias [1985] record. Although large changes in Delta delta13C occur in the equatorial Atlantic during early stages of the last three glacial cycles, surface-deep Delta delta13C at glacial maxima (18O stage 2, late stage 6, and late stage 8) was only about 0.2? greater than during the subsequent interglacial. Our results imply that nutrient-driven pCO2 changes account for about one third of the pCO2 decrease observed in ice cores, and consequently, Delta delta13C should not be used as a proxy pCO2 index. Enough variance in the ice core pCO2 records remains to be explained that conclusions about pCO2 and ice volume phase relationships should also be reexamined. As much as 40 ppm pCO2 change still has not been accounted for by models of past physics and chemistry of the ocean.