Model results of ocean physical field of the reference state and the targets in NetCDF format

In a feasibility study, the potential of proxy data for the temperature and salinity during the Last Glacial Maximum (LGM, about 19 000 to 23 000 years before present) in constraining the strength of the Atlantic meridional overturning circulation (AMOC) with a general ocean circulation model was explored. The proxy data were simulated by drawing data from four different model simulations at the ocean sediment core locations of the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) project, and perturbing these data with realistic noise estimates. The results suggest that our method has the potential to provide estimates of the past strength of the AMOC even from sparse data, but in general, paleo-sea-surface temperature data without additional prior knowledge about the ocean state during the LGM is not adequate to constrain the model. On the one hand, additional data in the deep-ocean and salinity data are shown to be highly important in estimating the LGM circulation. On the other hand, increasing the amount of surface data alone does not appear to be enough for better estimates. Finally, better initial guesses to start the state estimation procedure would greatly improve the performance of the method. Indeed, with a sufficiently good first guess, just the sea-surface temperature data from the MARGO project promise to be sufficient for reliable estimates of the strength of the AMOC.

Morphology of Orbulina universa and Globorotalia scitula across local extinctions during Sapropel S5 (Eastern Mediterranean Sea, c.126-121 ka)

Extinction is a remarkably difficult phenomenon to study under natural conditions. This is because the outcome of stress exposure and associated fitness reduction is not known until the extinction occurs and it remains unclear whether there is any phenotypic reaction of the exposed population that can be used to predict its fate. Here we take advantage of the fossil record, where the ecological outcome of stress exposure is known. Specifically, we analyze shell morphology of planktonic Foraminifera in sediment samples from the Mediterranean, during an interval preceding local extinctions. In two species representing different plankton habitats, we observe shifts in trait state and decrease in variance in association with non-terminal stress, indicating stabilizing selection. At terminal stress levels, immediately before extinction, we observe increased growth asymmetry and trait variance, indicating disruptive selection and bet-hedging. The pre-extinction populations of both species show a combination of trait states and trait variance distinct from all populations exposed to non-terminal levels of stress. This finding indicates that the phenotypic history of a population may allow the detection of threshold levels of stress, likely to lead to extinction. It is thus an alternative to population dynamics in studying and monitoring natural population ecology.

Interactive effects of light, nitrogen source, and carbon dioxide on energy metabolism in the diatom Thalassiosira pseudonana

Due to the ongoing effects of climate change, phytoplankton are likely to experience enhanced irradiance, more reduced nitrogen, and increased water acidity in the future ocean. Here, we used Thalassiosira pseudonana as a model organism to examine how phytoplankton adjust energy production and expenditure to cope with these multiple, interrelated environmental factors. Following acclimation to a matrix of irradiance, nitrogen source, and CO2 levels, the diatom's energy production and expenditures were quantified and incorporated into an energetic budget to predict how photosynthesis was affected by growth conditions. Increased light intensity and a shift from inline image to inline image led to increased energy generation, through higher rates of light capture at high light and greater investment in photosynthetic proteins when grown on inline image. Secondary energetic expenditures were adjusted modestly at different culture conditions, except that inline image utilization was systematically reduced by increasing pCO2. The subsequent changes in element stoichiometry, biochemical composition, and release of dissolved organic compounds may have important implications for marine biogeochemical cycles. The predicted effects of changing environmental conditions on photosynthesis, made using an energetic budget, were in good agreement with observations at low light, when energy is clearly limiting, but the energetic budget over-predicts the response to inline image at high light, which might be due to relief of energetic limitations and/or increased percentage of inactive photosystem II at high light. Taken together, our study demonstrates that energetic budgets offered significant insight into the response of phytoplankton energy metabolism to the changing environment and did a reasonable job predicting them.

Geochemical and physical sediment properties of sediment cores from the Portuguese shelf

Marine sediments from the Portuguese shelf are influenced by environmental changes in the surrounding continental and marine environment. These are largely controlled by the North Atlantic Oscillation, but additional impacts may arise from episodic tsunamis. In order to investigate these influences, a high resolution multi-proxy study has been carried out on a 5.4 m long gravity core and five box cores from the Tagus prodelta on the western Portuguese margin, incorporating geochemical (Corg/Ntotal ratios, d13Corg, d15N, d18O, Corg and CaCO3 content) and physical sediment properties (magnetic susceptibility, grain-size). Subsurface data of the five box cores indicate no major effect of early postdepositional alteration. Surface data show a higher fraction of terrigenous organic material close to the river mouth and in the southern prodelta. Gravity core GeoB 8903 covers the last 3.2 kyrs with a temporal resolution of at least 0.1 cm/yr. Very high sedimentation rates between 69 and 140 cm core depth indicate a possible disturbance of the record by the AD1755 tsunami, although no evidence for a disturbance is observed in the data. Sea surface temperature and salinity on the prodelta, the local budget of marine NO3- as well as the provenance of organic matter remained virtually constant during the past 3.2 kyrs. A positive correlation between magnetic susceptibility (MS) and North Atlantic Oscillation (NAO) is evident for the past 250 years, coinciding with a negative correlation between mean grain-size and NAO. This is assigned to a constant riverine supply of fine material with high MS, which is diluted by the riverine input of a coarser, low-MS component during NAO negative, high-precipitation phases. End-member modelling of the lithic grain-size spectrum supports this, revealing a third, coarse lithic component. The high abundance of this coarse end-member prior to 2 kyr BP is interpreted as the result of stronger bottom currents, concentrating the coarse sediment fraction by winnowing. As continental climate was more arid prior to 2 kyr BP (Subboreal), the coarse end-member may also consist of dust from local sources. A decrease in grain-size and CaCO3 content after 2 kyr BP is interpreted as a result of decreasing wind strength. The onset of a fining trend and a further decrease in CaCO3 around AD900 occurs simultaneous to climatic variations, reconstructed from eastern North Atlantic records. A strong increase in MS between AD1400 and AD1500 indicates higher lithic terrigenous input, caused by deforestation in the hinterland.

Rock magnetic properties and titanomagnetite oxydation state of ODP Hole 109-648B (Table 1)

The titanomagnetite oxidation state of "zero age" ocean floor basalts was investigated. For this purpose the oxidation parameter, z, of Hole 648B basalts was determined by SEM observation of "shrinkage cracks" in individual titanomagnetite grains and by Curie temperature measurements. A mean z-value of 0.1 has been deduced for the Hole 648B basalts. Assuming a linear relationship between titanomagnetite low-temperature oxidation state and age of the oceanic basalt, an age of 0.7 m.y. is deduced for Hole 648B.

Planktonic foraminiferal (Globigerinoides ruber) Mg/Ca ratios and SSTs of sediment cores GeoB9501-4 and GeoB9501-5

Multidecadal variations in Atlantic sea surface temperatures (SST) influence the climate of the Northern Hemisphere. However, prior to the instrumental time period, information on multidecadal climate variability becomes limited, and there is a particular scarcity of sufficiently resolved SST reconstructions. Here, we present an eastern tropical North Atlantic reconstruction of SSTs based on foraminiferal (Globigerinoides ruber pink) Mg/Ca ratios that resolves multidecadal variability over the past 1700 years. Spectral power in the multidecadal band (50 to 70 years period) is significant over several time intervals suggesting that the Atlantic Multidecadal Oscillation (AMO) has been influencing local SST. Since our data exhibit high scatter the absence of multidecadal variability in the remaining record does not exclude the possibility that SST variations on this time scale might have been present without being detected in our data. Cooling by ~0.5 °C takes place between about AD 1250 and AD 1500; while this corresponds to the inception of the Little Ice Age (LIA), the end of the LIA is not reflected in our record and SST remains relatively low. This transition to cooler SSTs parallels the previously reconstructed shift in the North Atlantic Oscillation towards a low pre-20th century mean state and possibly reflects common solar forcing.

Geochemistry and oxidation state of Lau Basin basalts

The backarc glasses recovered during Ocean Drilling Program Leg 135 are unique among submarine tholeiitic glasses with respect to their oxygen fugacity and sulfur concentrations. Unlike mid-ocean-ridge basalt glasses, fO2 in these samples (inferred from ratios Fe3+/Fe2+) is high and variable, and S variations (90-1140 ppm) are not coupled with FeO concentration. Strong correlations occur between the alkali and alkaline-earth elements and both fO2 (positive correlations) and S concentrations (negative correlations). Correlations between fO2 and various trace elements are strongest for those elements with a known affinity for hydrous fluids (perhaps produced during slab dehydration), suggesting the presence of a hydrous fluid with high fO2 and high alkali and alkaline earth element concentrations in the Lau Basin mantle. Concentrations of S and fO2 are strongly correlated; high fO2 samples are characterized by low S in addition to high alkali and alkaline earth element concentrations. The negative correlations between S and these trace elements are not consistent with incompatible behavior of S during crystallization. Mass balance considerations indicate that the S concentrations cannot result simply from mixing between low-S and high-S sources. Furthermore, there is no relationship between S and other trace elements or isotope ratios that might indicate that the S variations reflect mixing processes. The S variations more likely reflect the fact that when silicate coexists with an S-rich vapor phase the solubility of S in the silicate melt is a function of fO2 and is at a minimum at the fO2 conditions recorded by these glasses. The absence of Fe-sulfides and the high and variable vesicle contents are consistent with the idea that S concentrations reflect silicate-vapor equilibria rather than silicate-sulfide equilibria (as in MORB). The low S contents of some samples, therefore, reflect the high fO2 of the supra-subduction zone environment rather than a low-S source component.

Iron isotope and chemical compositions of sediments and basalts drilled seaward of the Mariana Trench

The purpose of this work is to study the mobility and budget of Fe isotopes in the oceanic crust and in particular during low-temperature interaction of seawater with oceanic basalt. We carried out this investigation using samples from Ocean Drilling Program (ODP) Site 801C drilled during Leg 129 and Leg 185 in Jurassic Pacific oceanic crust seaward of the Mariana Trench. The site comprises approximately 450 m of sediment overlying a section of 500 m of basalt, which includes intercalated pelagic and chemical sediments in the upper basaltic units and two low-temperature (10-30°C) ocherous Si-Fe hydrothermal deposits. Fe was chemically separated from 70 selected samples, and 57Fe/54Fe ratios were measured by MC-ICP-MS Isoprobe. The isotopic ratios were measured relative to an internal standard solution and are reported relative to the international Fe-standard IRMM-14. Based on duplicate measurements of natural samples, an external precision of 0.2? (2 sigma) has been obtained. The results indicate that the deep-sea sediment section has a restricted range of d57Fe, which is close to the igneous rock value. In contrast, large variations are observed in the basaltic section with positive d57Fe values (up to 2.05?) for highly altered basalts and negative values (down to ?2.49?) for the associated alteration products and hydrothermal deposits. Secondary Fe-minerals, such as Fe-oxyhydroxides or Fe-bearing clays (celadonite and saponite), have highly variable d57Fe values that have been interpreted as resulting from the partial oxidation of Fe(2+) leached during basalt alteration and precipitated as Fe(3+)-rich minerals. In contrast, altered basalts at Site 801C, which are depleted in Fe (up to 80%), display an increase in d57Fe values relative to fresh values, which suggest a preferential leaching of light iron during alteration. The apparent fractionation factor between dissolved Fe(2+) and Fe remaining in the mineral is from 0.5? to 1.3? and may be consistent with a kinetic isotope fractionation where light Fe is stripped from the minerals. Alternatively, the formation of secondary clays minerals, such as celadonite during basalt alteration may incorporate preferentially the heavy Fe isotopes, resulting in the loss of light Fe isotopes in the fluids. Because microbial processes within the oceanic crust are of potential importance in controlling rates of chemical reactions, Fe redox state and Fe-isotope fractionation, we evaluated the possible effect of this deep biosphere on Fe-isotope signatures. The Fe-isotope systematics presented in this study suggest that, even though iron behavior during seafloor weathering may be mediated by microbes, such as iron-oxidizers, d57Fe variations of more than 4? may also be explained by abiotic processes. Further laboratory experiments are now required to distinguish between various processes of Fe-isotope fractionation during seafloor weathering.

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.

Standing stocks and carbon isotopes of live benthic foraminifera from the South Atlantic

Live (Rose Bengal stained) and dead benthic foraminifera of surface and subsurface sediments from 25 stations in the eastern South Atlantic Ocean and the Atlantic sector of the Southern Ocean were analyzed to decipher a potential influence of seasonally and spatially varying high primary productivity on the stable carbon isotopic composition of foraminiferal tests. Therefore, stations were chosen so that productivity strongly varied, whereas conservative water mass properties changed only little. To define the stable carbon isotopic composition of dissolved inorganic carbon (d13CDIC) in ambient water masses, we compiled new and previously published d13CDIC data in a section running from Antarctica through Agulhas, Cape and Angola Basins, via the Guinea Abyssal Plain to the Equator. We found that intraspecific d13C variability of all species at a single site is constantly low throughout their distribution within the sediments, i.e. species specific and site dependent mean values calculated from all subbottom depths on average only varied by +/-0.09 per mil. This is important because it makes the stable carbon isotopic signal of species independent of the particular microhabitat of each single specimen measured and thus more constant and reliable than has been previously assumed. So-called vital and/or microhabitat effects were further quantified: (1) d13C values of endobenthic Globobulimina affinis, Fursenkoina mexicana, and Bulimina mexicana consistently are by between -1.5 and -1.0 per mil VPDB more depleted than d13C values of preferentially epibenthic Fontbotia wuellerstorfi, Cibicidoides pachyderma, and Lobatula lobatula. (2) In contrast to the Antarctic Polar Front region, at all stations except one on the African continental slope Fontbotia wuellerstorfi records bottom water d13CDIC values without significant offset, whereas L. lobatula and C. pachyderma values deviate from bottom water values by about -0.4 per mil and -0.6 per mil, respectively. This adds to the growing amount of data on contrasting cibicid d13C values which on the one hand support the original 1:1-calibration of F. wuellerstorfi and bottom water d13CDIC, and on the other hand document severe depletions of taxonomically close relatives such as L. lobatula and C. pachyderma. At one station close to Bouvet Island at the western rim of Agulhas Basin, we interpret the offset of -1.5 per mil between bottom water d13CDIC and d13C values of infaunal living Bulimina aculeata in contrast to about -0.6 +/- 0.1 per mil measured at eight stations close-by, as a direct reflection of locally increased organic matter fluxes and sedimentation rates. Alternatively, we speculate that methane locally released from gas vents and related to hydrothermal venting at the mid-ocean ridge might have caused this strong depletion of 13C in the benthic foraminiferal carbon isotopic composition. Along the African continental margin, offsets between deep infaunal Globobulimina affinis and epibenthic Fontbotia wuellerstorfi as well as between shallow infaunal Uvigerina peregrina and F. wuellerstorfi, d13C values tend to increase with generally increasing organic matter decomposition rates. Although clearly more data are needed, these offsets between species might be used for quantification of biogeochemical paleogradients within the sediment and thus paleocarbon flux estimates. Furthermore, our data suggest that in high-productivity areas where sedimentary carbonate contents are lower than 15 weight %, epibenthic and endobenthic foraminiferal d13C values are strongly influenced by 13C enrichment probably due to carbonate-ion undersaturation, whereas above this sedimentary carbonate threshold endobenthic d13C values reflect depleted pore water d13CDIC values.

Mg/Ca ratio of Globorotalia tumida and subsurface temperature reconstruction for ODP Site 130-806

The tropical Pacific thermocline strength, depth, and tilt are critical to tropical mean state and variability. During the early Pliocene (~3.5 to 4.5 Ma), the Eastern Equatorial Pacific (EEP) thermocline was deeper and the cold tongue was warmer than today, which resulted in an mean state with a reduced zonal sea surface temperature gradient, or El Padre. However, it is unclear whether the deep thermocline was a local feature of the EEP or a basin-wide condition with global implications. Our measurements of Mg/Ca of Globorotalia tumida in a western equatorial Pacific site indicate Pliocene subsurface temperatures warmer than today; thus, El Padre included a basin-wide thermocline that was relatively warm, deep, and weakly tilted. At ~4 Ma, thermocline steepening was coupled to cooling of the cold tongue. Since ~4 Ma, the basin-wide thermocline cooled/shoaled gradually, with implications for thermocline feedbacks in tropical dynamics and the interpretation of TEX86-derived temperatures.