Multidecadal increase in North Atlantic coccolithophores and the potential role of rising CO2

See allHide authors and affiliations

Science  26 Nov 2015:
DOI: 10.1126/science.aaa8026


As anthropogenic CO2 emissions acidify the oceans, calcifiers generally are expected to be negatively affected. However, using data from the Continuous Plankton Recorder, we show that coccolithophore occurrence in the North Atlantic increased from ~2 to over 20% from 1965 through 2010. We used Random Forest models to examine >20 possible environmental drivers of this change, finding that CO2 and the Atlantic Multidecadal Oscillation were the best predictors, leading us to hypothesize that higher CO2 levels might be encouraging growth. A compilation of 41 independent laboratory studies supports our hypothesis. Our study shows a long-term basin-scale increase in coccolithophores and suggests that increasing CO2 and temperature have accelerated the growth of a phytoplankton group that is important for carbon cycling.

Marine organisms that produce external features made of calcium carbonate are susceptible to harmful consequences from ocean acidification (1). Coccolithophores, the main calcifying phytoplankton, are unicellular algae surrounded by calcite plates called coccoliths, whose photosynthesis is strongly carbon-limited (2). Coccoliths are a major source of oceanic particulate inorganic carbon (PIC) and serve as ballast for sinking aggregates (3), thus accelerating carbon export (4). Given increasing partial pressures of atmospheric CO2 (pCO2), global warming, and ocean acidification, it is expected that coccolithophores will be affected, producing concomitant impacts on ocean carbon fluxes, dimethyl sulfide fluxes (5), carbonate geochemistry (6), and phytoplankton community structure (6). Current evidence regarding how increased pCO2 will affect coccolithophores is contradictory (710). Most laboratory manipulations study how coccolithophores respond to the increased pCO2 levels predicted for the end of the century rather than to the CO2 changes observed in the last five decades.

Here we report changes in the occurrence of coccolithophores in the North Atlantic during the last 45 years, and use Random Forest (RF) statistical models to evaluate the importance of various environmental drivers for these changes.

The in situ Continuous Plankton Recorder (CPR) surveys were developed to sample plankton in the North Atlantic using ships-of-opportunity. The surveys have followed the same methodology since 1946 (11). Sample preservation methods (using Borax-buffered formalin) and analysis have remained unchanged since 1958 (12), producing a unique, consistent, multi-decadal dataset. While the CPR filtering system was designed to sample larger microplankton, coccolithophores are trapped, particularly in the intersection of the silk fibers (12). It is not possible to accurately quantify organisms that are smaller than the mesh-size, but we can use the dataset to estimate the probability of coccolithophore occurrence. While our sampling underestimates natural abundances, this probability is a proxy for changes in coccolithophore abundance (fig. S1).

We calculated the annual probability of coccolithophore occurrence as the fraction of samples per year containing coccolithophores. The CPR data shows an increase in occurrence of coccolithophores across the North Atlantic from ~1% of samples in the 1960s to over 20% of samples with coccolithophores in the 2000s (Fig. 1, A to F, and fig. S2). Regional abundances of coccolithophores in the 2000s are at least ten times higher than those observed at the beginning of the data record. Our observations are supported by a shift in the opal:carbonate ratio in sediment traps in the Atlantic from the 1990s (13), satellite evidence of global poleward expansion of Emiliania huxleyi (14), and recurring blooms in areas where coccolithophores were previously absent or sparse (1417).

Fig. 1 Distribution of CPR samples (red) and coccolithophore observations (green) per decade.

(A) 1960s, (B) 1970s, (C) 1980s, (D) 1990s, (E) 2000s, (F) comparison of annual sampling effort (total number of samples per year in the North Atlantic) versus average probability of coccolithophore occurrence in raw data (= sum of samples with coccolithophore records per year / total number of samples per year × 100). Each sample corresponds to observations found in 10 cm of CPR silk (or ~3 m3 of water).

To evaluate possible top-down and bottom-up drivers for the increase in coccolithophore occurrence in the North Atlantic, we investigated factors that could affect coccolithophore growth rates and biogeography. Temperature, nutrient availability, light levels, competition, and predation are critical on a local scale. In turn, these may be affected by large-scale processes such as climate modes, global warming, and increases in CO2. The CPR sampling survey is irregular in time and space, making classic time series analysis inappropriate for this dataset. Additionally, the effects of different environmental forcings on phytoplankton groups are nonlinear and interdependent. After evaluating a suite of statistical methods (see the supplementary materials), we selected Random Forest (RF) models (18), an increasingly popular method in ecology that characterizes structure in high dimensional data while making no distributional assumptions about the response variable or predictors. RF has the advantage of allowing for nonlinearities, geographically and temporally discontinuous data, and the ability to model complex interactions among predictor variables without overfitting the data.

Our RF model predicted the probability of coccolithophore occurrence, defined as the percentage of samples containing coccolithophores in a 1-by-1-degree area each month, as a function of over 20 biological and physical predictors. Since the CPR dataset is already complex and discontinuous, we only used in situ measurements of biological and physical parameters without interpolating data. The complete dataset included 81,340 observations from 1965-2010. The importance of each variable in predicting coccolithophore occurrence is ranked in a variable importance plot. Importance is measured as the percent increase of Mean Squared Error (MSE) in the prediction due to that variable. Partial dependence plots graphically represent the marginal effect of each variable on the response variable. Here we discuss the top predictors and how these could have driven the observed increase in coccolithophore occurrence.

The evaluation of CO2 as a driver for observed increases in coccolithophore occurrence required special attention. While pCO2 and CPR data were not collected simultaneously, there is information about spatial and long-term CO2 variability. We created three random forests, each with a different predictor for pCO2: 1) global pCO2 estimated from Mauna Loa that varies only in time (“Global_pCO2” variable in the RF_GLOBAL model), 2) climatological ΔpCO2 from Takahashi et al. (19) that has no interannual variability (“delta_pCO2” variable in the RF_CLIM model), and 3) local pCO2 estimated as the sum of the Mauna Loa and Takahashi dataset varying in time and space (“Local_pCO2” variable in the RF_LOCAL model). Local pCO2 in the North Atlantic can be lower than the global pCO2 by 100 ppm in certain months (supplementary materials, fig. S3).

Local and global pCO2 were the best predictors of coccolithophore occurrence in our RF models (Fig. 2, A and B). In the RF_GLOBAL analysis, the partial dependence plot shows low coccolithophore occurrence (3-5%) at pCO2 ranges between 320-360 ppm. This corresponds to atmospheric CO2 observed at Mauna Loa from 1960 to ~1996. The increase in coccolithophore probability accelerated at pCO2 > 370 ppm (~1997), reaching 22% at 400 ppm. In the RF_LOCAL model (Fig. 2, B and E, and fig. S4), the partial dependence plot shows a three-fold range, with the largest coccolithophore probabilities corresponding to the highest CO2 values. Maximum probabilities of finding coccolithophores coincided at global and local CO2 levels above 360 ppm (Fig. 2, D and E). As discussed below, such dependence falls within the envelope of laboratory responses of coccolithophore growth rates to increased CO2.

Fig. 2 Random Forest analyses.

(A to C) Variable importance plots for RF_GLOBAL, RF_LOCAL, and RF_CLIM models, respectively. A large increase in percent mean squared error (%inc MSE) means a variable is a better predictor of coccolithophore occurrence). Boxes highlight the ranking of CO2 parameters within each model. (D to F) Partial dependence plots for CO2 parameters in RF_GLOBAL, RF_LOCAL, and RF_CLIM models, respectively. (G) Partial dependence plot for Atlantic Multidecadal Oscillation (AMO) based on the probability of coccolithophore occurrence in each of the random forest models. This is the marginal effect of the predictor on coccolithophore probability. Global_pCO2, atmospheric pCO2 based on Mauna Loa records; delta_pCO2, climatological ΔpCO2 (19); local_pCO2, estimate of local pCO2 based on the sum of spatial and long-term trends. MEI, Multivariate ENSO Index; U and Vcomp, zonal and meridional wind components; SST, sea-surface-temperature; tintinnid_pres, tintinnid occurrence.

Within RF_CLIM, the climatological “delta_pCO2” variable is not a strong predictor (Fig. 2, C and F). Instead, the Atlantic Multidecadal Oscillation (AMO) becomes the top predictor in this model and other climate modes (AO, Arctic Oscillation, MEI, Multivariate ENSO Index) rise in their ranking relative to RF_LOCAL and RF_GLOBAL, reflecting the importance of interannual variability. Recent studies have linked the AMO with phytoplankton (20) and coccolithophore variability (15). The AMO index tracks temperature anomalies in the North Atlantic and its positive trend in recent decades could mask global warming or enhance CO2 effects (Fig. 3). AMO ranked 13th in the RF_GLOBAL model, second in RF_LOCAL and first in the RF_CLIM analyses (Fig. 2, A to C, respectively). We propose two explanations: either AMO has a true effect on coccolithophore abundance, or in the absence of “global CO2” or “local CO2”, AMO is the only other variable with a similar long-term trend in the last two decades (Fig. 3). The highest coccolithophore probabilities are found during the recent positive AMO phases (Fig. 2G). If the coccolithophore increase were due to positive AMO, their occurrence should have been high during previous positive phases i.e., 1960s (Fig. 3). Unfortunately, the scarce information prior to 1965 limits our ability to draw strong conclusions, but based on the variability explained and the ranking in partial dependence plots (Fig. 2, A to G), we propose that AMO could be an important secondary driver. AMO has been related to changes in the Meridional Overturning Circulation that would alter nutrient supply (21). Positive AMO periods are associated with greater upward transport of nutrients in convective regions, but lower upward transport elsewhere. If AMO were the only mechanism responsible for the coccolithophore increase, it would be expected to produce opposite effects in the northwestern and eastern regions. Instead, all regions show increasing trends.

Fig. 3 Long-term trends.

(A) Annual basin-averaged coccolithophore probability in CPR samples (= sum of samples with coccolithophore records per year / total number of samples per year). (B) Global atmospheric CO2 measured from Mauna Loa. (C) Atlantic Multidecadal Oscillation (AMO). (D and E) Annual mean basin-averaged diatom and dinoflagellate counts in CPR samples per year. Vertical lines marking years 1965, 1997 are included for reference.

We hypothesize that synergistic effects due to CO2, AMO, and global warming differentially accelerated coccolithophore growth rates, driving recent increases in their occurrence. Compared to other phytoplankton groups, coccolithophore photosynthesis is severely carbon-limited (2) and sedimentary records show a predominance of coccolithophores during interglacial (22) and high CO2 periods (23). Many studies agree that coccolithophores respond to an increase in CO2 by decreasing PIC and increasing POC but there is disagreement with respect to the effects on growth rates (supplementary materials). We assembled a compendium of published growth rates as a function of CO2 (41 laboratory experiments from 16 independent publications, Fig. 4, fig. S5, supplementary materials). Results show a quasi-hyperbolic increase in coccolithophore growth rates with pCO2, with scatter partly produced by differences in experimental treatments: temperature, species, strain, nutrients, and irradiance. Our local pCO2 estimates between 1965-2010 (blue box in Fig. 4) correspond to the ranges in pCO2 where, based on Fig. 4, we would expect changes in coccolithophore growth rates. This compilation reconciles previous contradictory conclusions on the effects of CO2 on coccolithophore growth rates (supplementary materials), buttressing our hypothesis that CO2 enhances coccolithophore growth. Additional RF analyses of other top-down and bottom-up processes (grazing, nutrients, temperature, etc.) are included in the supplementary materials.

Fig. 4 Relationship between CO2 and coccolithophore growth rates based on 41 experiments from 16 publications and 4 species (details in supplementary materials).

Results are binned in 50 ppm intervals, with minimum, quartiles, median, and maximum per interval used to construct box-and-whisker plots. Discontinuous lines represent schematic boundaries of this relationship depending on irradiance, nutrient, and temperature levels. Blue box, range of local oceanic pCO2 values observed across the North Atlantic during this time period; red arrow, global atmospheric equivalent.

In order to project future coccolithophore abundances under elevated CO2 levels, we need to reassess the baseline. Our results show that today’s numbers are an order of magnitude greater than those in the 1960s and will likely continue to increase before growth rates stabilize at ~500 ppm. This is critical for understanding changes in the export ratio, biological pump, and alkalinity pump. Our compilation suggests that the changes seen in the North Atlantic may represent a global trend. Contrary to the generalized assumption of negative effects of ocean acidification on calcifiers, coccolithophores may be capable of adapting to a high-CO2 world (24), especially given evidence of highly calcified coccolithophores in areas with seasonally high pCO2 or low pH (7, 8). Coccolithophores show outstanding competitive abilities under the stratified, warm, nutrient-depleted conditions projected for the future ocean (6). Nevertheless, with increasing pCO2, we might expect changes in community composition and a decrease in calcification, leading to changes in rain ratio, export efficiency, and trophic effects higher in the marine food web.

Supplementary Materials

Materials and Methods

Supplementary Text

Figs. S1 to S7

Table S1

References (2547)

References and Notes

  1. Acknowledgments: We would like to thank D. Johns, the SAHFOS Associated External Researcher fund, and everyone involved in the CPR survey. We thank APL for financial support, and B. Zaitchik, D. Waugh, J. Dunne, and three anonymous reviewers for helpful comments. Support for W.M.B. came from NSF and NASA. CPR data used in this study can be obtained from SAHFOS upon request.
View Abstract


Navigate This Article