Although coccolithophore physiological responses to CO2-induced changes in seawater carbonate chemistry have been widely studied in the past, there is limited knowledge on the variability of physiological responses between populations from different areas. In the present study, we investigated the specific responses of growth, particulate organic (POC) and inorganic carbon (PIC) production rates of three populations of the coccolithophore Emiliania huxleyi from three regions in the North Atlantic Ocean (Azores: six strains, Canary Islands: five strains, and Norwegian coast near Bergen: six strains) to a CO2 partial pressure (pCO2) range from 120 to 2630 µatm. Physiological rates of each population and individual strain increased with rising pCO2 levels, reached a maximum and declined thereafter. Optimal pCO2 for growth, POC production rates, and tolerance to low pH (i.e., high proton concentration) was significantly higher in an E. huxleyi population isolated from the Norwegian coast than in those isolated near the Azores and Canary Islands. This may be due to the large environmental variability including large pCO2 and pH fluctuations in coastal waters off Bergen compared to the rather stable oceanic conditions at the other two sites. Maximum growth and POC production rates of the Azores and Bergen populations were similar and significantly higher than that of the Canary Islands population. This pattern could be driven by temperature–CO2 interactions where the chosen incubation temperature (16 °C) was slightly below what strains isolated near the Canary Islands normally experience. Our results indicate adaptation of E. huxleyi to their local environmental conditions and the existence of distinct E. huxleyi populations. Within each population, different growth, POC, and PIC production rates at different pCO2 levels indicated strain-specific phenotypic plasticity. Accounting for this variability is important to understand how or whether E. huxleyi might adapt to rising CO2 levels.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2022) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2023-07-05.