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.

(Supplementary Material 2) Landmark coordinates of the analyzed specimens with associated metadata

The use of planktonic foraminifera in paleoceanographic studies relies on the assumption that morphospecies represent biological species with ecological preferences that are stable through time and space. However, genetic surveys unveiled a considerable level of diversity in most morphospecies of planktonic foraminifera. This diversity is significant for paleoceanographic applications because cryptic species were shown to display distinct ecological preferences that could potentially help refine paleoceanographic proxies. Subtle morphological differences between cryptic species of planktonic foraminifera have been reported, but so far their applicability within paleoceanographic studies remains largely unexplored. Here we show how information on genetic diversity can be transferred to paleoceanography using Globorotalia inflata as a case study. The two cryptic species of G. inflata are separated by the Brazil-Malvinas Confluence (BMC), a major oceanographic feature in the South Atlantic. Based on this observation, we developed a morphological model of cryptic species detection in core top material. The application of the cryptic species detection model to Holocene samples implies latitudinal oscillations in the position of the confluence that are largely consistent with reconstructions obtained from stable isotope data. We show that the occurrence of cryptic species in G. inflata, can be detected in the fossil record and used to trace the migration of the BMC. Since a similar degree of morphological separation as in G. inflata has been reported from other species of planktonic foraminifera, the approach presented in this study can potentially yield a wealth of new paleoceanographical proxies.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Simultaneous high-resolution pH and spectrophotometric recordings of oxygen binding in native blood microvolumes, link to supplementary material

Oxygen equilibrium curves have been widely used to understand oxygen transport in numerous organisms. A major challenge has been to monitor oxygen binding characteristics and concomitant pH changes as they occur in vivo, in limited sample volumes. Here we report a technique allowing highly resolved and simultaneous monitoring of pH and blood pigment saturation in minute blood volumes. We equipped a gas diffusion chamber with a broad range fibre optic spectrophotometer and a micro-pH optode and recorded changes of pigment oxygenation along PO2 and pH gradients to test the setup. Oxygen binding parameters derived from measurements in only 15 µl of haemolymph from the cephalopod Octopus vulgaris showed low instrumental error (0.93%) and good agreement with published data. Broad range spectra, each resolving 2048 data points, provided detailed insight into the complex absorbance characteristics of diverse blood types. After consideration of photobleaching and intrinsic fluorescence, pH optodes yielded accurate recordings and resolved a sigmoidal shift of 0.03 pH units in response to changing PO2 from 0-21 kPa. Highly resolved continuous recordings along pH gradients conformed to stepwise measurements at low rates of pH changes. In this study we showed that a diffusion chamber upgraded with a broad range spectrophotometer and an optical pH sensor accurately characterizes oxygen binding with minimal sample consumption and manipulation. We conclude that the modified diffusion chamber is highly suitable for experimental biologists who demand high flexibility, detailed insight into oxygen binding as well as experimental and biological accuracy combined in a single set up.

Data S1. Results of the analyses considering developmental time, supplementary figures and tables & Video S1

Our ability to project the impact of global change on marine ecosystem is limited by our poor understanding on how to predict species sensitivity. For example, the impact of ocean acidification is highly species-specific, even in closely related taxa. The aim of this study was to test the hypothesis that the tolerance range of a given species to decreased pH corresponds to their natural range of exposure. Larvae of the green sea urchin Strongylocentrotus droebachiensis were cultured from fertilization to metamorphic competence (29 days) under a wide range of pH (from pHT = 8.0/pCO2 ~ 480 ?atm to pHT = 6.5/pCO2 ~ 20 000 ?atm) covering present (from pHT 8.7 to 7.6), projected near-future variability (from pHT 8.3 to 7.2) and beyond. Decreasing pH impacted all tested parameters (mortality, symmetry, growth, morphometry and respiration). Development of normal, although showing morphological plasticity, swimming larvae was possible as low as pHT >= 7.0. Within that range, decreasing pH increased mortality and asymmetry and decreased body length (BL) growth rate. Larvae raised at lowered pH and with similar BL had shorter arms and a wider body. Relative to a given BL, respiration rates and stomach volume both increased with decreasing pH suggesting changes in energy budget. At the lowest pHs (pHT <= 6.5), all the tested parameters were strongly negatively affected and no larva survived past 13 days post fertilization. In conclusion, sea urchin larvae appeared to be highly plastic when exposed to decreased pH until a physiological tipping point at pHT = 7.0. However, this plasticity was associated with direct (increased mortality) and indirect (decreased growth) consequences for fitness.

Impacts of seawater acidification on mantle gene expression patterns of the Baltic Sea blue mussel: implications for shell formation and energy metabolism, link to supplementary material

Marine organisms have to cope with increasing CO2 partial pressures and decreasing pH in the oceans. We elucidated the impacts of an 8-week acclimation period to four seawater pCO2 treatments (39, 113, 243 and 405 Pa/385, 1,120, 2,400 and 4,000 µatm) on mantle gene expression patterns in the blue mussel Mytilus edulis from the Baltic Sea. Based on the M. edulis mantle tissue transcriptome, the expression of several genes involved in metabolism, calcification and stress responses was assessed in the outer (marginal and pallial zone) and the inner mantle tissues (central zone) using quantitative real-time PCR. The expression of genes involved in energy and protein metabolism (F-ATPase, hexokinase and elongation factor alpha) was strongly affected by acclimation to moderately elevated CO2 partial pressures. Expression of a chitinase, potentially important for the calcification process, was strongly depressed (maximum ninefold), correlating with a linear decrease in shell growth observed in the experimental animals. Interestingly, shell matrix protein candidate genes were less affected by CO2 in both tissues. A compensatory process toward enhanced shell protection is indicated by a massive increase in the expression of tyrosinase, a gene involved in periostracum formation (maximum 220-fold). Using correlation matrices and a force-directed layout network graph, we were able to uncover possible underlying regulatory networks and the connections between different pathways, thereby providing a molecular basis of observed changes in animal physiology in response to ocean acidification.

Stable isotopic composition of Poritres lutea from New Caledonia

A 335 year stable isotope record from a New Caledonia coral (22°S, 166°E) helps fill a large gap in historical climate reconstructions. Although the long-term coral s18O-based sea surface temperature (SST) trend is one of warming, there are notable decadal fluctuations, especially in the early 18th and early 19th centuries. Mean annual SSTs between 1658 and 1900 are estimated to be ~0.3°C lower than the 20th century average, with interdecadal excursions of 0.5°-0.8°C. Time series analyses of the coral isotope record reveals significant concentrations of variance in the El Niño band; an inderdecadal spectral peak is present, but its robustness requires additional statistical evaluation. A secular but irregular decrease in coral d13C values begins in the mid-1800s and may reflect the anthropogenic perturbation of the carbon reservoir. These and other results indicate that the New Caledonia coral isotope record is a valuable source of information on southwest Pacific climate history.