Benthic foraminiferal fauna and isotope record in Holocene sediments of the northwest Indian Ocean

Historically, the Holocene has been considered an interval of relatively stable climate. However, recent studies from the northern Arabian Sea (Netherlands Indian Ocean Program 905) suggested high-amplitude climate shifts in the early and middle Holocene based on faunal and benthic isotopic proxy records. We examined benthic foraminiferal faunal and stable isotopic data from Ocean Drilling Program (ODP) Site 723 and total organic carbon data from ODP Site 724, Oman Margin (808 and 593 m water depths, respectively). At Site 723 the mid-Holocene shift in d18O values of infaunal benthic species Uvigerina peregrina (1.4 per mil) is 3 times larger than that of epifaunal benthic species Cibicides kullenbergi recorded at Site NIOP 905 off Somalia. However, none of the five other benthic species we measured at Hole 723A exhibits such a shift in d18O. We speculate that the late Holocene d18O decrease in U. peregrina represents species-specific changes in ecological habitat or food preference in response to changes in surface and deep ocean circulation. While the stable isotopic data do not appear to indicate a middle Holocene climatic shift, our total organic carbon and benthic faunal assemblage data do indicate that the early Holocene deep Arabian Sea was influenced by increased ventilation perhaps by North Atlantic Deep Water and/or Circumpolar Deep Water incursions into the Indian Ocean, leading to remineralization of organic matter and a relatively weak early Holocene oxygen minimum zone in the northwest Arabian Sea in spite of strong summer monsoon circulation.

Radiocarbon ages and plateau definitions of sediments from the abyssal ocean

Ice core records demonstrate a glacial-interglacial atmospheric CO2 increase by ~100 ppm, while 14C calibration efforts document a strong decrease in atmospheric 14C concentration during this period. A calculated transfer of ~530 Gt of 14C depleted carbon is required to produce the deglacial coeval rise of carbon in the atmosphere and terrestrial biosphere. This amount is usually ascribed to oceanic carbon release, although the actual mechanisms remained elusive, since an adequately old and carbon-enriched deep-ocean reservoir seemed unlikely. Here we present a new, though still fragmentary, ocean-wide d14C dataset showing that during the Last Glacial Maximum (LGM) and Heinrich Stadial 1 (HS-1) the maximum 14C age difference between ocean deep waters and the atmosphere exceeded the modern values by up to 1500 14C yr, in the extreme reaching 5100 14C yr. Below 2000 m depth the 14C ventilation age of modern ocean waters is directly linked to the concentration of dissolved inorganic carbon (DIC). We propose as working hypothesis that the modern regression of DIC vs d14C also applies for LGM times, which implies that a mean LGM aging by ~600 14C yr corresponded to a global rise of ~85-115 µmol DIC/kg in the deep ocean. Thus, the prolonged residence time of ocean deep waters may indeed have made it possible to absorb an additional ~730-980 Gt DIC, one third of which possibly originated from intermediate waters. We also infer that LGM deep-water O2 dropped to suboxic values of <10µmol/kg in the Atlantic sector of the Southern Ocean, possibly also in the subpolar North Pacific. The outlined deglacial transfer of the extra aged, deep-ocean carbon to the atmosphere via the dynamic ocean-atmosphere carbon exchange would be sufficient to account for two trends observed, (1) for the increase in atmospheric CO2 and (2) for the 190-permil drop in atmospheric d14C during the so-called HS-1 'Mystery Interval', when atmospheric 14C production rates were largely constant.

3He/4He and 4He/20Ne ratios in pore water and gas in Pacific deep-sea sediments

We have measured the 3He/3He and 3He/20 Ne ratios of thirty-nine pore water and gas samples in deep-sea sediments collected at twelve sites on the Pacific Ocean bottom during the cruises of Deep Sea Drilling Project Legs 87, 89, 90 and 92. The 3He/4He and 4He/20Ne ratios vary from 0.000000215 to 0.00000165 and from 0.29 to 20, respectively. He in the sample is composed of four components: (1) atmospheric He dissolved in seawater; (2) atmospheric He with mantle-derived He in Pacific bottom water; (3) in situ radiogenic He in the sediment; and (4) crustal He in the basement rock. Assuming that the 20Ne contents are constant with the value of seawater, the depth variations in the 4He/20Ne ratios at five Sites, 583D, 594, 597A, 598A and 504B, may provide useful information on 4He flux at the ocean bottom. The estimated 4He fluxes vary from 2000 to 40000 atoms cm**-2 s**-1 and are one to three orders of magnitude less than those calculated from the excess He in deep ocean water. An overall similarity between the geographical distribution of the 3He/4He ratios and heat flow data is found in the study area, between the East Pacific Rise across the Pacific Ocean and the Japanese Islands. The tendency is well explained by a conventional sea-floor spreading model.

Deep water formation in the North Pacific and deglacial CO2 rise

Deep water formation in the North Atlantic and Southern Ocean is widely thought to influence deglacial CO2 rise and climate change; here we suggest that deep water formation in the North Pacific may also play an important role. We present paired radiocarbon and boron isotope data from foraminifera from sediment core MD02-2489 at 3640 m in the North East Pacific. These show a pronounced excursion during Heinrich Stadial 1, with benthic-planktic radiocarbon offsets dropping to ~350 years, accompanied by a decrease in benthic d11B. We suggest this is driven by the onset of deep convection in the North Pacific, which mixes young shallow waters to depth, old deep waters to the surface, and low-pH water from intermediate depths into the deep ocean. This deep water formation event was likely driven by an increase in surface salinity, due to subdued atmospheric/monsoonal freshwater flux during Heinrich Stadial 1. The ability of North Pacific Deep Water (NPDW) formation to explain the excursions seen in our data is demonstrated in a series of experiments with an intermediate complexity Earth system model. These experiments also show that breakdown of stratification in the North Pacific leads to a rapid ~30 ppm increase in atmospheric CO2, along with decreases in atmospheric d13C and D14C, consistent with observations of the early deglaciation. Our inference of deep water formation is based mainly on results from a single sediment core, and our boron isotope data are unavoidably sparse in the key HS1 interval, so this hypothesis merits further testing. However we note that there is independent support for breakdown of stratification in shallower waters during this period, including a minimum in d15N, younging in intermediate water 14C, and regional warming. We also re-evaluate deglacial changes in North Pacific productivity and carbonate preservation in light of our new data, and suggest that the regional pulse of export production observed during the Bølling-Allerød is promoted by relatively stratified conditions, with increased light availability and a shallow, potent nutricline. Overall, our work highlights the potential of NPDW formation to play a significant and hitherto unrealized role in deglacial climate change and CO2 rise.

Dissolved and dissolvable iron concentrations and ligand characteristics in samples taken during POLARSTERN expedition ARK-XXII/2 to the Arctic Ocean

The speciation of iron was investigated in three shelf seas and three deep basins of the Arctic Ocean in 2007. The dissolved fraction (<0.2 µm) and a fraction < 1000 kDa were considered here. In addition, unfiltered samples were analyzed. Between 74 and 83% of dissolved iron was present in the fraction < 1000 kDa at all stations and depth, except at the chlorophyll maximum (42-64%). Distinct trends in iron concentrations and ligand characteristics were observed from the shelf seas toward the central deep basins, with a decrease of total dissolvable iron ([TDFe] > 3 nM on the shelves and [TDFe] < 2 nM in the Makarov Basin). A relative enrichment of particulate Fe toward the bottom was revealed at all stations, indicating Fe export toward the deep ocean. In deep waters, dissolved ligands became less saturated with Fe (increase of [Excess L]/[Fe]) from the Nansen Basin via the Amundsen Basin toward the Makarov Basin. This trend was explained by the reactivity of the ligands, higher (log alpha > 13.5) in the Nansen and Amundsen basins than in the Makarov Basin (log alpha <13) where the sources of Fe and ligands were limited. The ligands became nearly saturated with depth in the Amundsen and Nansen Basins, favoring Fe removal in the deep ocean, whereas in the deep Makarov Basin, they became unsaturated with depth. Still here scavenging occurred. Although scavenging of Fe was attenuated by the presence of unsaturated organic ligands, their low reactivity in combination with a lack of sources of Fe in the Makarov Basin might be the reason of a net export of Fe to the sediment.

Chemical composition of volcanic glass of some tephra layers in Northeastern Pacific Ocean

Five widespread upper Cenozoic tephra layers that are found within continental sediments of the western United States have been correlated with tephra layers in marine sediments in the Humboldt and Ventura basins of coastal California by similarities in major-and trace-element abundances; four of these layers have also been identified in deep-ocean sediments at DSDP sites 34, 36, 173, and 470 in the northeastern Pacific Ocean. These layers, erupted from vents in the Yellowstone National Park area of Wyoming and Idaho (Y), the Cascade Range of the Pacific Northwest (C), and the Long Valley area, California (L), are the Huckleberry Ridge ash bed (2.0 Ma, Y), Rio Dell ash bed (ca. 1.5 Ma, C), Bishop ash bed (0.74 Ma, L), Lava Creek B ash bed (0.62 Ma, Y), and Loleta ash bed (ca. 0.4 Ma, C). The isochronous nature of these beds allows direct comparison of chronologic and climatic data in a variety of depositional environments. For example, the widespread Bishop ash bed is correlated from proximal localities near Bishop in east-central California, where it is interbedded with volcanic and glacial deposits, to lacustrine beds near Tecopa, southeastern California, to deformed on-shore marine strata near Ventura, southwestern California, to deep-ocean sediments at site 470 in the eastern Pacific Ocean west of northern Mexico. The correlations allow us to compare isotopic ages determined for the tephra layers with ages of continental and marine biostratigraphic zones determined by magnetostratigraphy and other numerical age control and also provide iterative checks for available age control. Relative age variations of as much as 0.5 m.y. exist between marine biostratigraphic datums [for example, highest occurrence level of Discoaster brouweri and Calcidiscus tropicus (= C. macintyrei)], as determined from sedimentation rate curves derived from other age control available at each of several sites. These discrepancies may be due to several factors, among which are (1) diachronism of the lowest and highest occurrence levels of marine faunal and floral species with latitude because of ecologic thresholds, (2) upward reworking of older forms in hemipelagic sections adjacent to the tectonically active coast of the western United States and other similar analytical problems in identification of biostratigraphic and magnetostratigraphic datums, (3) dissolution of microfossils or selective diagenesis of some taxa, (4) lack of precision in isotopic age calibration of these datums, (5) errors in isotopic ages of tephra beds, and (6) large variations in sedimentation rates or hiatuses in stratigraphic sections that result in age errors of interpolated datums. Correlation of tephra layers between on-land marine and deep-ocean deposits indicates that some biostratigraphic datums (diatom and calcareous nannofossil) may be truly time transgressive because at some sites, they are found above and, at other sites, below the same tephra layers.

A deep-sea coral record of the North Atlantic

Our record of Younger Dryas intermediate-depth seawater D14C from North Atlantic deep-sea corals supports a link between abrupt climate change and intermediate ocean variability. Our data show that northern source intermediate water (~1700 m) was partially replaced by 14C-depleted southern source water at the onset of the event, consistent with a reduction in the rate of North Atlantic Deep Water formation. This transition requires the existence of large, mobile gradients of D14C in the ocean during the Younger Dryas. The D14C water column profile from Keigwin (2004) provides direct evidence for the presence of one such gradient at the beginning of the Younger Dryas (~12.9 ka), with a 100 per mil offset between shallow (<~2400 m) and deep water. Our early Younger Dryas data are consistent with this profile and also show a D14C inversion, with 35 per mil more enriched water at ~2400 m than at ~1700 m. This feature is probably the result of mixing between relatively well 14C ventilated northern source water and more poorly 14C ventilated southern source intermediate water, which is slightly shallower. Over the rest of the Younger Dryas our intermediate water/deepwater coral D14C data gradually increase, while the atmosphere D14C drops. For a very brief interval at ~12.0 ka and at the end of the Younger Dryas (11.5 ka), intermediate water D14C (~1200 m) approached atmospheric D14C. These enriched D14C results suggest an enhanced initial D14C content of the water and demonstrate the presence of large lateral D14C gradients in the intermediate/deep ocean in addition to the sharp vertical shift at ~2500 m. The transient D14C enrichment at ~12.0 ka occurred in the middle of the Younger Dryas and demonstrates that there is at least one time when the intermediate/deep ocean underwent dramatic change but with much smaller effects in other paleoclimatic records.

Global Ocean Particulate Organic Carbon flux merged with satellite parameters

The efficiency of the biological pump of carbon to the deep ocean depends largely on the biologically mediated export of carbon from the surface ocean and its remineralization with depth. Global satellite studies have primarily focused on chlorophyll concentration and net primary production (NPP) to understand the role of phytoplankton in these processes. Recent satellite retrievals of phytoplankton composition now allow for the size of phytoplankton cells to be considered. Here, we improve understanding of phytoplankton size structure impacts on particle export, remineralization and transfer. Particulate organic carbon (POC) flux observations from sediment traps and 234Th are compiled across the global ocean. Annual climatologies of NPP, percent microplankton, and POC flux at four time series locations and within biogeochemical provinces are constructed, and sinking velocities are calculated to align surface variables with POC flux at depth. Parameters that characterize POC flux vs. depth (export flux ratio, labile fraction, remineralization length scale) are then fit to the aligned dataset. Times of the year dominated by different size compositions are identified and fit separately in regions of the ocean where phytoplankton cell size showed enough dynamic range over the annual cycle. Considering all data together, our findings support the paradigm of high export flux but low transfer efficiency in more productive regions and vice versa for oligotrophic regions. However, when parsing by dominant size class, we find periods dominated by small cells to have both greater export flux and lower transfer efficiency than periods when large cells comprise a greater proportion of the phytoplankton community.

Objective Mapping of Argo data in the Weddell Gyre: a gridded dataset of upper ocean water properties, link to data files in NetCDF format

The Weddell Gyre plays a crucial role in the regulation of climate by transferring heat into the deep ocean through deep and bottom water mass formation. However, our understanding of Weddell Gyre water mass properties is limited to regions of data availability, primarily along the Prime Meridian. The aim is to provide a dataset of the upper water column properties of the entire Weddell Gyre. Objective mapping was applied to Argo float data in order to produce spatially gridded, time composite maps of temperature and salinity for fixed pressure levels ranging from 50 to 2000 dbar, as well as temperature, salinity and pressure at the level of the sub-surface temperature maximum. While the data are currently too limited to incorporate time into the gridded structure, the data are extensive enough to produce maps of the entire region across three time composite periods (2002-2005, 2006-2009 and 2010-2013), which can be used to determine how representative conclusions drawn from data collected along general RV transect lines are on a gyre scale perspective. The time composite data sets are provided as netCDF files; one for each time period. Mapped fields of conservative temperature, absolute salinity and potential density are provided for 41 vertical pressure levels. The above variables as well as pressure are provided at the level of the sub-surface temperature maximum. Corresponding mapping errors are also included in the netCDF files. Further details are provided in the global attributes, such as the unit variables and structure of the corresponding data array (i.e. latitude x longitude x vertical pressure level). In addition, all files ending in "_potTpSal" provide mapped fields of potential temperature and practical salinity.

Mid-Pliocene deep-sea bottom water temperatures based on ostracode Mg/Ca ratios

We studied magnesium:calcium (Mg/Ca) ratios in shells of the deep-sea ostracode genus Krithe from a short interval in the middle Pliocene between 3.29 and 2.97 Ma using deep-sea drilling sites in the North and South Atlantic in order to estimate bottom water temperatures (BWT) during a period of climatic warmth. Results from DSDP and ODP Sites 552A, 610A, 607, 658A, 659A, 661A and 704 for the period Ma reveal both depth and latitudinal gradients of mean Mg/Ca values. Shallower sites (552A, 610A and 607) have higher mean Mg/Ca ratios (10.3, 9.7, 10.1 mmol/mol) than deeper sites (661A, 6.3 mmol/mol), and high latitude North Atlantic sites (552A, 610 and 607) have higher Mg/Ca ratios than low latitude (658A: 9.8 mmol/mol, 659A: 7.7 mmol/mol, 661A: 6.3 mmol/mol) and Southern Ocean (704: 8.0 mmol/mol) sites. Converting Mg/Ca ratios into estimated temperatures using the calibration of Dwyer et al. (1995) [Dwyer, G.S., Cronin, T.M., Baker, P.A., Raymo, M.E., Buzas, J.S., Corrège, T., 1995. North Atlantic deepwater temperature change during late Pliocene and late Quaternary climatic cycles. Science 270, 1347-1351] suggests that mean middle Pliocene bottom water temperatures at the study sites in the deep Atlantic were about the same as modern temperatures. However, brief pulses of elevated BWT occurred several times between 3.29 and 2.97 Ma in both the North and South Atlantic Ocean suggesting short-term changes in deep ocean circulation.

Global distributions of epipelagic macrozooplankton abundance and biomass - Gridded data product (NetCDF) - Contribution to the MAREDAT World Ocean Atlas of Plankton Functional Types

Macrozooplankton are an important link between higher and lower trophic levels in the oceans. They serve as the primary food for fish, reptiles, birds and mammals in some regions, and play a role in the export of carbon from the surface to the intermediate and deep ocean. Little, however, is known of their global distribution and biomass. Here we compiled a dataset of macrozooplankton abundance and biomass observations for the global ocean from a collection of four datasets. We harmonise the data to common units, calculate additional carbon biomass where possible, and bin the dataset in a global 1 x 1 degree grid. This dataset is part of a wider effort to provide a global picture of carbon biomass data for key plankton functional types, in particular to support the development of marine ecosystem models. Over 387 700 abundance data and 1330 carbon biomass data have been collected from pre-existing datasets. A further 34 938 abundance data were converted to carbon biomass data using species-specific length frequencies or using species-specific abundance to carbon biomass data. Depth-integrated values are used to calculate known epipelagic macrozooplankton biomass concentrations and global biomass. Global macrozooplankton biomass has a mean of 8.4 µg C l-1, median of 0.15 µg C l-1 and a standard deviation of 63.46 µg C l-1. The global annual average estimate of epipelagic macrozooplankton, based on the median value, is 0.02 Pg C. Biomass is highest in the tropics, decreasing in the sub-tropics and increasing slightly towards the poles. There are, however, limitations on the dataset; abundance observations have good coverage except in the South Pacific mid latitudes, but biomass observation coverage is only good at high latitudes. Biomass is restricted to data that is originally given in carbon or to data that can be converted from abundance to carbon. Carbon conversions from abundance are restricted in the most part by the lack of information on the size of the organism and/or the absence of taxonomic information. Distribution patterns of global macrozooplankton biomass and statistical information about biomass concentrations may be used to validate biogeochemical models and Plankton Functional Type models.

Gene expression changes in the coccolithophore Emiliania huxleyi after 500 generations of selection to ocean acidification

Coccolithophores are unicellular marine algae that produce biogenic calcite scales and substantially contribute to marine primary production and carbon export to the deep ocean. Ongoing ocean acidification particularly impairs calcifying organisms, mostly resulting in decreased growth and calcification. Recent studies revealed that the immediate physiological response in the coccolithophore Emiliania huxleyi to ocean acidification may be partially compensated by evolutionary adaptation, yet the underlying molecular mechanisms are currently unknown. Here, we report on the expression levels of 10 candidate genes putatively relevant to pH regulation, carbon transport, calcification and photosynthesis in E. huxleyi populations short-term exposed to ocean acidification conditions after acclimation (physiological response) and after 500 generations of high CO2 adaptation (adaptive response). The physiological response revealed downregulation of candidate genes, well reflecting the concomitant decrease of growth and calcification. In the adaptive response, putative pH regulation and carbon transport genes were up-regulated, matching partial restoration of growth and calcification in high CO2-adapted populations. Adaptation to ocean acidification in E. huxleyi likely involved improved cellular pH regulation, presumably indirectly affecting calcification. Adaptive evolution may thus have the potential to partially restore cellular pH regulatory capacity and thereby mitigate adverse effects of ocean acidification.

Comparative responses of two dominant Antarctic phytoplankton taxa to interactions between ocean acidification, warming, irradiance, and iron availability

We investigated the responses of the ecologically dominant Antarctic phytoplankton species Phaeocystis antarctica (a prymnesiophyte) and Fragilariopsis cylindrus (a diatom) to a clustered matrix of three global change variables (CO2, mixed-layer depth, and temperature) under both iron (Fe)-replete and Fe-limited conditions based roughly on the Intergovernmental Panel on Climate Change (IPCC) A2 scenario: (1) Current conditions, 39 Pa (380 ppmv) CO2, 50 µmol photons/m**2/s light, and 2°C; (2) Year 2060, 61 Pa (600 ppmv) CO2, 100 µmol photons/m**2/s light, and 4°C; (3) Year 2100, 81 Pa (800 ppmv) CO2, 150 µmol photons/m**2/s light, and 6°C. The combined interactive effects of these global change variables and changing Fe availability on growth, primary production, and cell morphology are species specific. A competition experiment suggested that future conditions could lead to a shift away from P. antarctica and toward diatoms such as F. cylindrus. Along with decreases in diatom cell size and shifts from prymnesiophyte colonies to single cells under the future scenario, this could potentially lead to decreased carbon export to the deep ocean. Fe : C uptake ratios of both species increased under future conditions, suggesting phytoplankton of the Southern Ocean will increase their Fe requirements relative to carbon fixation. The interactive effects of Fe, light, CO2, and temperature on Antarctic phytoplankton need to be considered when predicting the future responses of biology and biogeochemistry in this region.

Sinking rates of particles measured from deployments in the Atlantic Ocean based on seasonal data

The flux of materials to the deep sea is dominated by larger, organic-rich particles with sinking rates varying between a few meters and several hundred meters per day. Mineral ballast may regulate the transfer of organic matter and other components by determining the sinking rates, e.g. via particle density. We calculated particle sinking rates from mass flux patterns and alkenone measurements applying the results of sediment trap experiments from the Atlantic Ocean. We have indication for higher particle sinking rates in carbonate-dominated production systems when considering both regional and seasonal data. During a summer coccolithophorid bloom in the Cape Blanc coastal upwelling off Mauritania, particle sinking rates reached almost 570 m per day, most probably due the fast sedimentation of densely packed zooplankton fecal pellets, which transport high amounts of organic carbon associated with coccoliths to the deep ocean despite rather low production. During the recurring winter-spring blooms off NW Africa and in opal-rich production systems of the Southern Ocean, sinking rates of larger particles, most probably diatom aggregates, showed a tendency to lower values. However, there is no straightforward relationship between carbonate content and particle sinking rates. This could be due to the unknown composition of carbonate and/or the influence of particle size and shape on sinking rates. It also remains noticeable that the highest sinking rates occurred in dust-rich ocean regions off NW Africa, but this issue deserves further detailed field and laboratory investigations. We obtained increasing sinking rates with depth. By using a seven-compartment biogeochemical model, it was shown that the deep ocean organic carbon flux at a mesotrophic sediment trap site off Cape Blanc can be captured fairly well using seasonal variable particle sinking rates. Our model provides a total organic carbon flux of 0.29 Tg per year down to 3000 m off the NW African upwelling region between 5 and 35° N. Simple parameterisations of remineralisation and sinking rates in such models, however, limit their capability in reproducing the flux variation in the water column.