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.

Mean snow accumulation rate, and mean methane sulfonate, chloride and nitrate concentration in snow pits and firn cores of Dronning Maud Land

We quantified postdepositional losses of methane sulfonate (MSA-), nitrate, and chloride at the European Project for Ice Coring in Antarctica (EPICA) drilling site in Dronning Maud Land (DML) (75°S, 0°E). Analyses of four intermediate deep firn cores and 13 snow pits were considered. We found that about 26 ± 13% of the once deposited nitrate and typically 51 ± 20% of MSA- were lost, while for chloride, no significant depletion could be observed in firn older than one year. Assuming a first order exponential decay rate, the characteristic e-folding time for MSA- is 6.4 ± 3 years and 19 ± 6 years for nitrate. It turns out that for nitrate and MSA- the typical mean concentrations representative for the last 100 years were reached after 5.4 and 6.5 years, respectively, indicating that beneath a depth of around 1.2-1.4 m postdepositional losses can be neglected. In the area of investigation, only MSA- concentrations and postdepositional losses showed a distinct dependence on snow accumulation rate. Consequently, MSA- concentrations archived at this site should be significantly dependent on the variability of annual snow accumulation, and we recommend a corresponding correction. With a simple approach, we estimated the partial pressure of the free acids MSA, HNO3, and HCl on the basis of Henry's law assuming that ionic impurities of the bulk ice matrix are localized in a quasi-brine layer (QBL). In contrast to measurements, this approach predicts a nearly complete loss of MSA-, NO3 - , and Cl-.

Biogeochemistry measurements in a mesocosm experiment in the Bay of Hopavagen, Norway, in August 2012

The goal of this work has been to examine the influence of upper ocean food web structure and functioning on both the natural and artificially enhanced sequestration of carbon within the ocean. Data obtained in the mesocosm experiment run in the Bay of Hopavågen in August 2012 are used to assess the extent to which organic matter produced within four different food webs is retained in the upper ocean food web versus remineralized back to carbon dioxide and inorganic nutrients (ammonium, dissolved silicon, phosphate) versus exported from the system in the form of rapidly sinking particles. The experiment was carried out in a set of 12 mesocosms covering, in triplicate, 2 different phytoplankton communities (diatom versus non-diatom) exposed to 2 different zooplankton communities (-copepod and +copepod). These starting conditions were established by first filling the bags, roughly simultaneously, with seawater from the Bay of Hopavågen. Mesozooplankton were then removed to the most complete extent possible immediately removed from half of the mesocosms through repeated vertical hauls of a plankton net (200 µm mesh). Nitrate and phosphate was added to half mesocosms daily to promote the growth of non-siliceous phytoplankton (e.g. dinoflagellates or coccolithophores). To the other half of the mesocosms, nitrate, phosphate, and silicate were added to promote the growth of diatoms. Material was allowed to settle and the two distinct phytoplankton populations were allowed to develop for 4 days, after which copepods collected from the Bay of Hopavågen were added back to the half of the N+P mesocosms and to the half of the N+P+Si mesocosms from which mesozooplankton had not been removed at the beginning. This yielded a set of four initial starting conditions (N+P-copepods, N+P+copepods, N+P+Si-copepods, and N+P+Si+copepods). In the primary mesocosms, samples for a set of core parameters were taken every time the mesocosms were sampled. Samples for particulates (PIC, BSi, POC, PON) were collected on GF/F or 0.4 µm polycarbonate.

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.

Adaptation of a globally important coccolithophore to ocean warming and acidification

Regulating intracellular pH (pHi) is critical for optimising the metabolic activity of corals, yet mechanisms involved in pH regulation and the buffering capacity within coral cells are not well understood. Our study investigated how the presence of symbiotic dinoflagellates affects the response of pHi to pCO2-driven seawater acidification in cells isolated from Pocillopora damicornis. Using the fluorescent dye BCECF-AM, in conjunction with confocal microscopy, we simultaneously characterised the response of pHi in host coral cells and their dinoflagellate symbionts, in symbiotic and non-symbiotic states under saturating light, with and without the photosynthetic inhibitor DCMU. Each treatment was run under control (pH 7.8) and CO2 acidified seawater conditions (decreasing pH from 7.8 - 6.8). After two hours of CO2 addition, by which time the external pH (pHe) had declined to 6.8, the dinoflagellate symbionts had increased their pHi by 0.5 pH units above control levels. In contrast, in both symbiotic and non-symbiotic host coral cells, 15 min of CO2 addition (0.2 pH unit drop in pHe) led to cytoplasmic acidosis equivalent to 0.4 pH units. Despite further seawater acidification over the duration of the experiment, the pHi of non-symbiotic coral cells did not change, though in host cells containing a symbiont cell the pHi recovered to control levels. This recovery was negated when cells were incubated with DCMU. Our results reveal that photosynthetic activity of the endosymbiont is tightly coupled with the ability of the host cell to recover from cellular acidosis after exposure to high CO2 / low pH.

Determination of sulfite and manganese in sediments from the K/T boundary

The solid-state-physics technique of electron spin resonance (ESR) has been employed in an exploratory study of marine limestones and impact-related deposits from Cretaceous-Tertiary (KT) boundary sites including Spain (Sopelana and Caravaca), New Jersey (Bass River), the U.S. Atlantic continental margin (Blake Nose, ODP Leg 171B/1049/A), and several locations in Belize and southern Mexico within -600 km of the Chicxulub crater. The ESR spectra of SO3(1-) (a radiation-induced point defect involving a sulfite ion substitutional for CO3(2-) which has trapped a positive charge) and Mn(2+) in calcite were singled out for analysis because they are unambiguously interpretable and relatively easy to record. ESR signal strengths of calcite-related SO3(1-) and Mn(2+) have been studied as functions of stratigraphic position in whole-rock samples across the KT boundary at Sopelana, Caravaca, and Blake Nose. At all three of these sites, anomalies in SO3(1-) and/or Mn(2+) intensities are noted at the KT boundary relative to the corresponding background levels in the rocks above and below. At Caravaca, the SO3(1-) background itself is found to be lower by a factor of 2.7 in the first 30,000 years of the Tertiary relative to its steady-state value in the last 15,000 years of the Cretaceous, indicating either an abrupt and quasi-permanent change in ocean chemistry (or temperature) or extinction of the marine biota primarily responsible for fixing sulfite in the late Cretaceous limestones. An exponential decrease in the Mn(2+) concentration per unit mass calcite, [Mn(2+)], as the KT boundary at Caravaca is approached from below (1/e characteristic length =1.4 cm) is interpreted as a result of post-impact leaching of the seafloor. Absolute ESR quantitative analyses of proximal impact deposits from Belize and southern Mexico group naturally into three distinct fields in a twodimensional [SO3(1-)]-versus-[Mn(2+)] scatter plot. These fields contain (I) limestone ejecta clasts, (II) accretionary lapilli, and (III) a variety of SO3(1-) -depleted/Mn(2+) enriched impact deposits. Data for the investigated non-impact-related Cretaceous and Tertiary marine limestones (Spain and Blake Nose) fall outside of these three fields. With reference to thes enon-impact deposits, fields I, II, and III can be respectively characterized as Mn(2+) -depleted, SO3(1-) -enhanced, and SO3(1-) -depleted. It is proposed that (1) field I represents calcites from the Yucatin Platform, and that the Mn(2+) -depleted signature can be used as an indicator of primary Chicxulub ejecta in deep marine environments and (2) field II represents calcites that include a component formed in the vapor plume, either from condensation in the presence of CO2/SO3(1-) -rich vapors, or reactions between CaO and CO2/SO3 rich vapors, and that this SO3(1-) -enhanced signature can be used as an indicator of impact vapor plume deposits. Given these two propositions, the ESR data for the Blake Nose deposits are ascribed to the presence of basal coarse calcitic Chicxulub ejecta clasts, while the finer components that are increasingly represented toward the top are interpreted to contain high- SO3(1-) calcite from the vapor plume. The apparently-undisturbed Bass River deposit may contain even higher concentrations of vapor-plume calcite. None of the three components included in field III appear to be represented at distal, deep marine KT-boundary sites; this field may include several types of impact-related deposits of diverse origins and diagenetic histories.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, BATS

The contributions of total organic carbon and nitrogen to elemental cycling in the surface layer of the Sargasso Sea are evaluated using a 5-yr time-series data set (1994-1998). Surface-layer total organic carbon (TOC) and total organic nitrogen (TON) concentrations ranged from 60 to 70 µM C and 4 to 5.5 µM N seasonally, resulting in a mean C : N molar ratio of 14.4±2.2. The highest surface concentrations varied little during individual summer periods, indicating that net TOC production ceased during the highly oligotrophic summer season. Winter overturn and mixing of the water column were both the cause of concentration reductions and the trigger for net TOC production each year following nutrient entrainment and subsequent new production. The net production of TOC varied with the maximum in the winter mixed-layer depth (MLD), with greater mixing supporting the greatest net production of TOC. In winter 1995, the TOC stock increased by 1.4 mol C/m**2 in response to maximum mixing depths of 260 m. In subsequent years experiencing shallower maxima in MLD (<220 m), TOC stocks increased <0.7 mol C/m**2. Overturn of the water column served to export TOC to depth (>100 m), with the amount exported dependent on the depth of mixing (total export ranged from 0.4 to 1.4 mol C/m**2/yr). The exported TOC was comprised both of material resident in the surface layer during late summer (resident TOC) and material newly produced during the spring bloom period (fresh TOC). Export of resident TOC ranged from 0.5 to 0.8 mol C/m**2/yr, covarying with the maximum winter MLD. Export of fresh TOC varied from nil to 0.8 mol C/m**2/yr. Fresh TOC was exported only after a threshold maximum winter MLD of ~200 m was reached. In years with shallower mixing, fresh TOC export and net TOC production in the surface layer were greatly reduced. The decay rates of the exported TOC also covaried with maximum MLD. The year with deepest mixing resulted in the highest export and the highest decay rate (0.003 1/d) while shallow and low export resulted in low decay rates (0.0002 1/d), likely a consequence of the quality of material exported. The exported TOC supported oxygen utilization at dC : dO2 molar ratios ranging from 0.17 when TOC export was low to 0.47 when it was high. We estimate that exported TOC drove 15-41% of the annual oxygen utilization rates in the 100-400 m depth range. Finally, there was a lack of variability in the surface-layer TON signal during summer. The lack of a summer signal for net TON production suggests a small role for N2 fixation at the site. We hypothesize that if N2 fixation is responsible for elevated N : P ratios in the main thermocline of the Sargasso Sea, then the process must take place south of Bermuda and the signal transported north with the Gulf Stream system.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.

Nitrate limitation and ocean acidification interact with UV-B to reduce photosynthetic performance in the diatom

It has been proposed that ocean acidification (OA) will interact with other environmental factors to influence the overall impact of global change on biological systems. Accordingly we investigated the influence of nitrogen limitation and OA on the physiology of diatoms by growing the diatom Phaeodactylum tricornutum Bohlin under elevated (1000 µatm; high CO2- HC) or ambient (390 µatm; low CO2-LC) levels of CO2 with replete (110 µmol/L; high nitrate-HN) or reduced (10 ?mol/L; low nitrate-LN) levels of NO3- and subjecting the cells to solar radiation with or without UV irradiance to determine their susceptibility to UV radiation (UVR, 280-400 nm). Our results indicate that OA and UVB induced significantly higher inhibition of both the photosynthetic rate and quantum yield under LN than under HN conditions. UVA or/and UVB increased the cells' non-photochemical quenching (NPQ) regardless of the CO2 levels. Under LN and OA conditions, activity of superoxide dismutase and catalase activities were enhanced, along with the highest sensitivity to UVB and the lowest ratio of repair to damage of PSII. HC-grown cells showed a faster recovery rate of yield under HN but not under LN conditions. We conclude therefore that nutrient limitation makes cells more prone to the deleterious effects of UV radiation and that HC conditions (ocean acidification) exacerbate this effect. The finding that nitrate limitation and ocean acidification interact with UV-B to reduce photosynthetic performance of the diatom P. tricornutum implies that ocean primary production and the marine biological C pump will be affected by OA under multiple stressors.

Strong shift from HCO3- to CO2 uptake in Emiliania huxleyi with acidification: new approach unravels acclimation versus short-term pH effects

Effects of ocean acidification on Emiliania huxleyi strain RCC 1216 (calcifying, diploid life-cycle stage) and RCC 1217 (non-calcifying, haploid life-cycle stage) were investigated by measuring growth, elemental composition, and production rates under different pCO2 levels (380 and 950 µatm). In these differently acclimated cells, the photosynthetic carbon source was assessed by a (14)C disequilibrium assay, conducted over a range of ecologically relevant pH values (7.9-8.7). In agreement with previous studies, we observed decreased calcification and stimulated biomass production in diploid cells under high pCO2, but no CO2-dependent changes in biomass production for haploid cells. In both life-cycle stages, the relative contributions of CO2 and HCO3 (-) uptake depended strongly on the assay pH. At pH values =< 8.1, cells preferentially used CO2 (>= 90 % CO2), whereas at pH values >= 8.3, cells progressively increased the fraction of HCO3 (-) uptake (~45 % CO2 at pH 8.7 in diploid cells; ~55 % CO2 at pH 8.5 in haploid cells). In contrast to the short-term effect of the assay pH, the pCO2 acclimation history had no significant effect on the carbon uptake behavior. A numerical sensitivity study confirmed that the pH-modification in the (14)C disequilibrium method yields reliable results, provided that model parameters (e.g., pH, temperature) are kept within typical measurement uncertainties. Our results demonstrate a high plasticity of E. huxleyi to rapidly adjust carbon acquisition to the external carbon supply and/or pH, and provide an explanation for the paradoxical observation of high CO2 sensitivity despite the apparently high HCO3 (-) usage seen in previous studies.

Late Quaternary siliciclastic/carbonate sedimentation in the Cabricorn Channel

A model is presented for hemipelagic siliciclastic and carbonate sedimentation during the last glacial-interglacial cycle in the Capricorn Channel, southern Great Barrier Reef (GBR). Stable isotope ratios, grainsize, carbonate content and mineralogy were analysed for seven cores in a depth transect from 166 to 2892 m below sea level (mbsl). Results show variations in the flux of terrigenous, neritic and pelagic sediments to the continental slope over the last sea level cycle. During the glacial lowstand terrigenous sediment influenced all the cores down to 2000 mbsl. The percentages of quartz and feldspar in the cores decreased with water depth, while the percentage of clay increased. X-ray diffraction analysis of the glacial lowstand clay mineralogy suggests that the siliciclastic sediment was primarily sourced from the Fitzroy River, which debouched directly into the northwest sector of the Capricorn Channel at this time. The cores also show a decrease in pelagic calcite and an increase in aragonite and high magnesium calcite (HMC) during the glacial. The influx of HMC and aragonite is most likely from reworking of coral reefs exposed on the continental shelf during the glacial, and also from HMC ooids precipitated at the head of the Capricorn Channel at this time. Mass accumulation rates (MARs) are high (13.5 g/cm**/kyr) during the glacial and peak at ~20 g/cm** 3/kyr in the early transgression (16-14 ka BP). MARs then decline with further sea level rise as the Fitzroy River mouth retreats from the edge of the continental shelf after 13.5 ka BP. MARs remain low (4 g/cm**3/kyr) throughout the Holocene highstand. Data for the Holocene highstand indicate there is a reduction in siliciclastic influx to the Capricorn Channel with little quartz and feldspar below 350 mbsl. However, fine-grained fluvial sediments, presumably from the Fitzroy River, were still accumulating on the mid slope down to 2000 mbsl. The proportion of pelagic calcite in the core tops increases with water depth, while HMC decreases, and is present only in trace amounts in cores below 1500 mbsl. The difference in the percentage of HMC in the deeper cores between the glacial and Holocene may reflect differences in supply or deepening of the HMC lysocline during the glacial. Sediment accumulation rates also vary between cores in the Capricorn Channel and do not show the expected exponential decrease with depth. This may be due to intermediate or deep water currents reworking the sediments. It is also possible that present bathymetry data are too sparse to detect the potential role that submarine channels may play in the distribution and accumulation of sediments. Comparison of the Capricorn Channel MARs with those for other mixed carbonate/siliciclastic provinces from the northeast margin of Australia indicates that peak MARs in the early transgression in the Capricorn Channel precede those from the central GBR and south of Fraser Island. The difference in the timing of the carbonate and siliciclastic MAR peaks along the northeast margin is primarily related to differences in the physiography and climate of the provinces. The only common trend in the MARs from the northeast margin of Australia is the near synchronicity of the carbonate and siliciclastic MAR peaks in individual sediment cores, which supports a coeval sedimentation model.

Photophysiological responses of Southern Ocean phytoplankton to changes in CO2 concentrations: Short-term versus acclimation effects

The present study examines how different pCO2 acclimations affect the CO2- and light-dependence of photophysiological processes and O2 fluxes in four Southern Ocean (SO) key phytoplankton species. We grew Chaetoceros debilis (Cleve), Pseudo-nitzschia subcurvata (Hasle), Fragilariopsis kerguelensis (O'Meara) and Phaeocystis antarctica (Karsten) under low (160 µatm) and high (1000 ?atm) pCO2. The CO2- and light-dependence of fluorescence parameters of photosystem II (PSII) were determined by means of a fluorescence induction relaxation system (FIRe). In all tested species, nonphotochemical quenching (NPQ) is the primary photoprotection strategy in response to short-term exposure to high light or low CO2 concentrations. In C. debilis and P. subcurvata, PSII connectivity (p) and functional absorption cross-sections of PSII in ambient light (sigma PSII') also contributed to photoprotection while changes in re-oxidation times of Qa acceptor (tQa) were more significant in F. kerguelensis. The latter was also the only species being responsive to high acclimation pCO2, as these cells had enhanced relative electron transport rates (rETRs) and sigma PSII' while tQa and p were reduced under short-term exposure to high irradiance. Low CO2-acclimated cells of F. kerguelensis and all pCO2 acclimations of C. debilis and P. subcurvata showed dynamic photoinhibition with increasing irradiance. To test for the role and presence of the Mehler reaction in C. debilis and P. subcurvata, the light-dependence of O2 fluxes was estimated using membrane inlet mass spectrometry (MIMS). Our results show that the Mehler reaction is absent in both species under the tested conditions. We also observed that dark respiration was strongly reduced under high pCO2 in C. debilis while it remained unaltered in P. subcurvata. Our study revealed species-specific differences in the photophysiological responses to pCO2, both on the acclimation as well as the short-term level.

Experiment: Adaptation of a globally important coccolithophore to ocean warming and acidification

Although oceanwarming and acidification are recognized as two major anthropogenic perturbations of today's oceanswe know very little about how marine phytoplankton may respond via evolutionary change.We tested for adaptation to ocean warming in combination with ocean acidification in the globally important phytoplankton species Emiliania huxleyi. Temperature adaptation occurred independently of ocean acidifcation levels. Exponential growth rates were were up to 16% higher in populations adapted for one year to warming when assayed at their upper thermal tolerance limit. Particulate inorganic (PIC) and organic (POC) carbon production was restored to values under present-day ocean conditions, owing to adaptive evolution, and were 101% and 55% higher under combined warming and acidification, respectively, than in non-adapted controls. Cells also evolved to a smaller size while they recovered their initial PIC:POC ratio even under elevated CO2. The observed changes in coccolithophore growth, calcite and biomass production, cell size and elemental composition demonstrate the importance of evolutionary processes for phytoplankton performance in a future ocean. At the end of a 1-yr temperature selection phase, we conducted a reciprocal assay experiment in which temperature-adapted asexual populations were compared to the respective non-adapted control populations under high temperature, and vice versa (1. Assay Data, Dataset #835336). Mean exponential growth rates ? in treatments subjected to high temperature increased rapidly under all high temperature-CO2 treatment combinations during the temperature selection phase (2. time series, Dataset #835339).