Journal topic
Biogeosciences, 17, 635–647, 2020
https://doi.org/10.5194/bg-17-635-2020
Biogeosciences, 17, 635–647, 2020
https://doi.org/10.5194/bg-17-635-2020

Research article 07 Feb 2020

Research article | 07 Feb 2020

# The Arctic picoeukaryote Micromonas pusilla benefits from ocean acidification under constant and dynamic light

The Arctic picoeukaryote Micromonas pusilla benefits from ocean acidification under constant and dynamic light
Emily White1, Clara J. M. Hoppe1, and Björn Rost1,2 Emily White et al.
• 1Alfred-Wegener-Institut, Helmholtz-Zentrum für Polar- und Meeresforschung, 27570 Bremerhaven, Germany
• 2Fachbereich Biologie/Chemie, Universität Bremen, Leobener Strasse, 28359 Bremen, Germany

Correspondence: Emily White (ewhite14@msn.com)

Abstract

Compared to the rest of the globe, the Arctic Ocean is affected disproportionately by climate change. Despite these fast environmental changes, we currently know little about the effects of ocean acidification (OA) on marine key species in this area. Moreover, the existing studies typically test the effects of OA under constant, hence artificial, light fields. In this study, the abundant Arctic picoeukaryote Micromonas pusilla was acclimated to current (400 µatm) and future (1000 µatm) pCO2 levels under a constant as well as a dynamic light, simulating more realistic light fields as experienced in the upper mixed layer. To describe and understand the responses to these drivers, growth, particulate organic carbon (POC) production, elemental composition, photophysiology and reactive oxygen species (ROS) production were analysed. M. pusilla was able to benefit from OA on various scales, ranging from an increase in growth rates to enhanced photosynthetic capacity, irrespective of the light regime. These beneficial effects were, however, not reflected in the POC production rates, which can be explained by energy partitioning towards cell division rather than biomass build-up. In the dynamic light regime, M. pusilla was able to optimize its photophysiology for effective light usage during both low- and high-light periods. This photoacclimative response, which was achieved by modifications to photosystem II (PSII), imposed high metabolic costs leading to a reduction in growth and POC production rates when compared to constant light. There were no significant interactions observed between dynamic light and OA, indicating that M. pusilla is able to maintain effective photoacclimation without increased photoinactivation under high pCO2. Based on these findings, M. pusilla is likely to cope well with future conditions in the Arctic Ocean.

1 Introduction

Alterations to the ecosystem caused by climate change are far more pronounced in the Arctic than in the rest of the world (Pörtner et al., 2014). The increase in pCO2 and concomitant decrease in seawater pH, for instance, is particularly fast in the Arctic Ocean, which is mainly due to the higher solubility of CO2 at low water temperatures (Bates and Mathis, 2009). Many studies have investigated the effects of ocean acidification (OA) on phytoplankton and have observed species-specific responses on the physiological level (Lohbeck et al., 2012; Riebesell and Tortell, 2011; Rost et al., 2008). The negative effects of OA are thought to result from disturbed ion homeostasis under decreasing pH, while positive responses seem to be driven by the physiological mechanisms of inorganic carbon uptake (Bach et al., 2013; Rokitta et al., 2012). Photosynthesis requires CO2 as a substrate for the carbon-fixing enzyme RuBisCO, which is however characterized by a poor substrate affinity (Badger et al., 1998). To avoid CO2 limitation arising from this, larger phytoplankton especially typically depend on carbon-concentrating mechanisms (CCMs). These CCMs involve the transport of CO2 and/or ${\mathrm{HCO}}_{\mathrm{3}}^{-}$ into the cell and the prevention of leakage out of the cell, as well as the expression of carbonic anhydrase, an enzyme accelerating the interconversion between CO2 and ${\mathrm{HCO}}_{\mathrm{3}}^{-}$ (Reinfelder, 2011). As CCMs are energetically expensive, a potential downregulation under OA may be beneficial for phytoplankton (Hopkinson et al., 2011).

Also the rates of warming are 2 to 3 times faster than the global average (Trenberth et al., 2007), a phenomenon known as Arctic amplification (Miller et al., 2010). Warming causes changes to mixing regimes in the surface ocean (Houghton et al., 2001), probably leading to a shoaling of the mixed layer due to increased thermal stratification and freshening caused by sea-ice melting (Steinacher et al., 2010). Furthermore, models predict that the Arctic Ocean will become more frequently nearly ice-free during the summer months (Pachauri et al., 2014). With a decrease in sea-ice cover, the primary productivity in the Arctic is expected to increase, due to higher light availability and longer growing seasons for phytoplankton (Arrigo et al., 2008). At the same time, however, annual productivity may become increasingly limited by the low nutrient supply in the Arctic Ocean, which may decrease even further due to reduced upwelling (Tremblay et al., 2015; Wassmann and Reigstad, 2011). These environmental changes in the Arctic have already led to differences in community structure and are expected to cause more dramatic regime shifts in the future (e.g. Nöthig et al., 2015; Li et al., 2009). In addition to accelerated rates of OA and warming (Pörtner et al., 2014; Trenberth et al., 2007), the Arctic is affected by strong changes in wind and weather patterns (Hu and Bates, 2018), indicating that this region is disproportionately affected by overall climate change.

Next to climate-driven changes, phytoplankton growing in the turbulent upper mixed layer must generally acclimate to the dynamics of light availability. Due to the low-light periods in such dynamic light fields, phytoplankton, on the one hand, need to increase the light harvesting efficiency, e.g. by increasing the photosynthetic pigments like chlorophyll a (Chl a; Palmer et al., 2013). In the high-light periods, on the other hand, photoprotective mechanisms need to be activated to prevent photodamage to cells (Ragni et al., 2008). There are species-specific differences in the responses to differing light regimes, which can include changes in the number of photoprotective pigments, different photorepair mechanisms and modifications to the number of reaction centres in PSII (photosystem II; Ragni et al., 2008). Such acclimatory responses may be particularly important in the Arctic shelf seas, where high organic matter loading leads to particularly high light attenuation with depth (Granskog et al., 2012).

In the world's oceans, the picoplankton-size fraction (< 2 µm) is a significant contributor to overall productivity (Worden et al., 2015) and Micromonas-like picoeukaryotes are also highly abundant in the Arctic region (Lovejoy et al., 2007). Micromonas has been described as a shade-adapted species that can persist in Arctic winter darkness (Marquardt et al., 2016) with the help of mixotrophy (McKie-Krisberg and Sanders, 2014). Such low-light-adapted organisms are thought to show a lack of plasticity with regard to Chl a quota and photoacclimation (Talmy et al., 2013). Smaller phytoplankton are furthermore expected to particularly benefit from reduced nutrients under enhanced stratification due to their high surface : volume ratio, which allows them to take up nutrients more efficiently (Brussaard et al., 2013). For the same reasons, picoeukaryotes may benefit from elevated pCO2 levels due to increased CO2 diffusion into the cell. Concurrently, it has been shown in different experiments that smaller phytoplankton will thrive under future OA conditions (Brussaard et al., 2013; Engel et al., 2008; Hoppe et al., 2017, 2018; Maat et al., 2014; Meakin and Wyman, 2011; Schulz et al., 2017) and are regarded as potential “winners” of climate change (e.g. Hoppe et al., 2018; Li et al., 2009; Schulz et al., 2017).

Despite their prevalence in all marine habitats, the effect of natural light variability has been mostly ignored in previous studies on OA effects. Fluctuating light has however been reported to indeed affect phytoplankton photosynthesis and growth (Falkowski, 1980; Huisman, 1999; Köhler et al., 2018; Litchman, 2000; Litchman et al., 2004). Interactive effects between dynamic light regimes and OA have been observed in the coccolithophore Gephyrocapsa oceanica, causing decreased productivity (Jin et al., 2013). Hoppe et al. (2015) reported that the Antarctic diatom Chaetoceros debilis benefitted from OA under static irradiance, while a dynamic light regime reversed this positive response. This was attributed to the fact that OA-dependent downregulation of the CCM can expose cells to oxidative stress during high-light peaks under dynamic light (Hoppe et al., 2015). Oxidative stress occurs when the production of reactive oxygen species (ROS) exceeds the defensive mechanisms for ROS reduction, leading to accumulation in the cells (Apel and Hirt, 2004). In this study, the response of M. pusilla to OA was investigated under constant and dynamic light in order to determine whether there was an interactive effect of the two environmental factors. A particular focus was laid on the physiological mechanisms that determined the observed overall responses.

2 Materials and methods

## 2.1 Experimental setup

Monoclonal cultures of the picoeukaryote M. pusilla (isolated in 2014 by Klara Wolf in Kongsfjorden, Svalbard, 79 N) were grown in l L glass bottles in semi-continuous dilute batch cultures (max 158 000 cells mL−1; diluted every 3–5 d). The temperature remained stable at 2.6±0.2C. The media was composed of Arctic seawater (from Hausgarten observatory, Alfred-Wegener-Institut, 78 N, collected during an RV Maria S. Merian cruise in 2013) filtered through a 0.2 µm membrane filter capsule (Sartorius Stedim Biotech, Sartobran 300) and enriched with vitamins and trace metals in accordance with the F/2 protocol (Guillard and Ryther, 1962), as well as with macronutrients in Redfield proportions (containing 100 µmol L−1 nitrate and silicate and 6.2 µmol L−1 phosphate).

Both the constant and the dynamic light regime consisted of a 20 : 4 h light : dark cycle with an average light intensity of 83±5µmol photons m2 s−1. The dynamic light regime varied between 0 and 590 µmol photons m2 s−1 (Fig. 1). These light levels were calculated based on conditions typically observed in the Arctic Kongsfjorden (Svalbard, 79 N) in late spring, using maximum surface irradiance of 905 µmol photons m2 s−1, a mixed-layer depth of 20 m, an extinction coefficient of 0.35 m−1 (Clara J. M. Hoppe, unpublished results) and a vertical mixing rate of 0.011 m s−1 (Denman and Gargett, 1983). The dynamic light field acts to approximate natural light conditions, assuming homogenous mixing and diurnal changes in incoming irradiance but no weather-associated variability. Light was supplied through LED lamps (ECONLUX, Solar Stinger Sunstrip), and the dynamic light regime was regulated using a daylight controller (LED scaping light control, LiWeBe). In both setups, the light levels were monitored using a ULM-500 universal light meter with a 4π sensor (Effeltrich), and light intensity was adjusted with neutral density screens.

Figure 1Light regimes plotted as a function of time over a 24 h period. Indicated are the dynamic light cycle (solid line) and the constant light cycle (dashed line). Time points 1 and 2 are displayed at midday and in the evening at the start of the dark period (×).

The CO2 partial pressures (pCO2) were achieved through aeration of the incubation bottles with two different pCO2 levels (400 and 1000 µatm) for at least 12 h prior to inoculation. The gas mixtures were created using a gas flow controller (CGM 2000, Umwelttechnik MCZ), which mixed pure CO2 with CO2-free air to the desired pCO2 level. The pCO2 levels were monitored using a non-dispersive infrared analyser (LI-6252, LI-COR Biosciences). The humidified gas mixtures were bubbled through a glass frit and supplied via a sterile 0.2 µm PTFE filter (Midistart 2000, Sartorius Stedim). Cultures were acclimated to the respective pCO2 levels for at least five generations prior to the experiment. To minimize shifts in carbonate chemistry due to biomass production, cell densities were kept low at between 5000 and 158 000 cells mL−1.

## 2.2 Carbonate chemistry

Samples for seawater pH were collected at the start and middle of the experiment (during mid-light phase) and at the end of the experiment (during the first hour of the dark phase). Seawater pH was determined potentiometrically, using a 2-point-calibrated glass reference electrode (IoLine, SCHOTT Instruments) and pH meter (826 pH mobile, Metrohm), and was reported on the NBS scale for incubation temperatures. pH values were temperature-corrected using the CO2SYS program (Pierrot et al., 2006) to 2 C.

The samples for dissolved inorganic carbon (CT) measurements were gently filtered through a sterile 0.2 µm Nalgene syringe filter (Nalgene, Thermo Scientific) and stored in the dark at 2 C in 5 mL borosilicate bottles. The sample was subsequently analysed colorimetrically in duplicate using an autoanalyser (Seal Analytical; Stoll et al., 2001) with a reproducibility of ±8µmol kg−1 (Table 1). A certified reference standard material (CRM) was used to correct for measurement errors (Dickson et al., 2007). The final average CT values were 2141±23µmol kg−1 at ambient pCO2 levels and 2209±11µmol kg−1 under high pCO2 (Table 1).

The total alkalinity (AT) samples were gently filtered through pre-combusted 25 mm GF/F filters (glass microfiber filter, Whatman, GE Healthcare Life Sciences) and stored in 125 mL dark borosilicate bottles at 2 C. Standards and samples were equilibrated to room temperature prior to potentiometric titrations (Brewer et al., 1986) of two 25 mL subsamples with an autoanalyser (TitroLine alpha plus, SCHOTT Instruments). An internal standard was applied to correct for systematic errors based on measurements of CRMs, and the data were processed using TitriSoft 2.71 software. The corrected final AT values ranged between 2194±8 and 2215±5µmol kg−1 (Table 1).

The full carbonate system was calculated with a salinity of 32.2 and a temperature of 2 C using the pH and AT data with the CO2SYS program (Pierrot et al., 2006), following suggestions by Hoppe et al. (2012; Table 1). The calculations used constants of Mehrbach et al. (1973) with a refit by Dickson and Millero (1987) and a (B)T value according to Uppström (1974). The carbonate system remained stable for the duration of the experiment; i.e. the average daily pH values, calculated using the start, middle and end measurements, were 8.12±0.06 for ambient conditions and 7.82±0.03 for high pCO2 levels.

Table 1Carbonate chemistry measurements for each light and pCO2 treatment (n=4; mean ± 1 SD). The measured values are CT (dissolved inorganic carbon), AT (total alkalinity) and pH (NBS scale). pCO2 was calculated using the CO2SYS program, with pH and AT as input values. The values were calculated for 2 C, with a salinity of 32.2. The nutrient levels were 6.5 and 100 µmol kg−1 for PO4 and Si(OH)4, respectively.

## 2.3 Growth and cellular composition

Cell densities of M. pusilla were quantified using a flow cytometer (FCM; Accuri C6, BD Biosciences). Samples were analysed using live cells, where 490 µL of sample was added to 10 µL of 1 µm microspheres fluorescent-beads solution (Fluoresbrite YG, Polysciences Inc), which acted as an internal standard. Cells were identified and counted using the FL3 and FL4 channels as well as forward scatter for 2 min on slow speed with a maximum of 50 000 events. Specific growth rate constants (μ) were calculated from exponential fits of cell numbers over time for each replicate bottle. Samples were measured daily within a 1 h time frame for consistency.

Samples for particulate organic carbon (POC) and particulate organic nitrogen (PON) were collected at the end of the batch-culture experiment during the dark phase; samples were gently filtered onto pre-combusted 25 mm GF/F filters. Before analysis, 200 µL of hydrochloric acid (HCl, 0.2 M) was added to each filter and the samples were dried at 60 C for at least 12 h to remove any inorganic carbon. The samples were analysed using an elemental analyser (Euro EA 3000, HEKAtech). The POC and PON data were corrected by subtracting blank measurements, and values were normalized using the specific cell density and volume filtered to yield cell quotas. Subsequently, production rates were calculated by multiplying the quota with the specific growth rate constant (μ) of the respective incubation.

Samples for the Chl a quota were obtained at the end of the batch-culture experiment during the dark phase, by gentle filtration onto 25 mm GF/F filters, and were immediately stored at −20C until analysis. For chlorophyll extraction, 8 mL of 90 % acetone was added to the filters and subsequently stored at 4 C for at least 4 h in darkness. After centrifugation (4500 rpm for 5 min, Sigma 4K10), samples were measured on a fluorometer (TD-700 Fluorometer, Turner Designs) before and after acidification with HCl (1 M). Chl a concentrations (µg L−1) were calculated as in Knap et al. (1996).

## 2.4 Physiological responses

Photophysiological parameters were measured using a fast-repetition-rate fluorometer (FRRf; FastOcean sensor, Chelsea Technologies) in combination with a FastAct system (Chelsea Technologies). The fluorometer's light-emitting diodes (LEDs) were set to an emission wavelength of 450 nm. A saturation phase of 100 flashlets on a pitch of 2 µs was used, with a relaxation phase comprising 40 flashlets and a pitch of 50 µs. Prior to measurements, samples were dark acclimated for 15 min, and measurements were conducted in a temperature-controlled chamber at 3 C. The maximum (Fm, ${F}_{\mathrm{m}}^{\prime }$) and minimum (F0, F) chlorophyll fluorescence in the dark and light were estimated according to iterative algorithms for induction (Kolber et al., 1998) and relaxation phases (Oxborough et al., 2012). The PSII quantum yield efficiency was estimated as FvFm using the following equation:

$\begin{array}{}\text{(1)}& {F}_{\mathrm{v}}/{F}_{\mathrm{m}}=\left({F}_{\mathrm{m}}-{F}_{\mathrm{0}}\right)/{F}_{\mathrm{m}}.\end{array}$

Additional parameters were measured after dark acclimation, including the absorption cross-section size of PSII (σPSII; ${\mathit{\text{Å}}}^{\mathrm{2}}\cdot {q}^{-\mathrm{1}}$), the connectivity of PSII (ρ) and the PSII reopening rate (τ; ms), according to Kolber et al. (1998). ${F}_{\mathrm{0}}^{\prime }$ was estimated after Oxborough and Baker (1997) as

$\begin{array}{}\text{(2)}& {F}_{\mathrm{0}}^{\prime }={F}_{\mathrm{0}}/\frac{{F}_{\mathrm{v}}}{{F}_{\mathrm{m}}}+\frac{{F}_{\mathrm{0}}}{{F}_{\mathrm{m}}^{\prime }}.\end{array}$

Thereafter, the coefficient of photochemical quenching qL was calculated after Kramer et al. (2004) as

$\begin{array}{}\text{(3)}& \mathrm{qL}=\left(\left({F}_{\mathrm{m}}^{\prime }-{F}^{\prime }\right)\phantom{\rule{0.125em}{0ex}}/\left({F}_{\mathrm{m}}^{\prime }-{F}_{\mathrm{0}}^{\prime }\right)\right)\cdot \left({F}_{\mathrm{0}}^{\prime }/{F}^{\prime }\right).\end{array}$

The electron transport rates through PSII (ETR; mol e (mol RCII)−1 s−1) were calculated after Xu et al. (2017) using the following equation:

$\begin{array}{}\text{(4)}& \mathrm{ETR}={\mathit{\sigma }}_{\mathrm{PSII}}×\mathrm{qL}×\mathrm{PAR},\end{array}$

where σPSII is the absorption cross-section size of PSII, and PAR is the photosynthetically active radiation. Photosynthesis–irradiance (PI) curves were estimated at eight irradiances between 0 and 589 µmol photons m−2 s−1. According to the model by Webb et al. (1974), the light harvesting efficiency (α; mol e m2 (mol RCII)−1 (mol photons)−1) and the maximum relative electron transport rate (ETRmax; mol e (mol RCII)−1 s−1) were estimated using the following equation:

$\begin{array}{}\text{(5)}& \mathrm{ETR}={\mathrm{ETR}}_{max}\left[\mathrm{1}-e\left(\frac{-\mathit{\alpha }I}{{\mathrm{ETR}}_{max}}\right)\right].\end{array}$

The light saturation index (Ik; µmol photons m−2 s−1) was calculated as ${\mathrm{ETR}}_{max}/\mathit{\alpha }$.

At the light level of 506 µmol photons m−2 s−1, non-photochemical quenching (NPQ) was calculated as Y(NPQ) using calculations as described in Klughammer and Schreiber (2008):

$\begin{array}{}\text{(6)}& Y\left(\mathrm{NPQ}\right)=\frac{F}{{F}_{\mathrm{m}}^{\prime }}-\frac{F}{{F}_{\mathrm{m}}}.\end{array}$

Measurements of ROS for both $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ free radicals and H2O2 were assessed using the FCM with the fluorochromes dihydroethidium (HE; D7008, Sigma-Aldrich) and dihydrorhodamine 123 (DHR123; D1054, Sigma-Aldrich), respectively. Methods were adapted from Prado et al. (2012), with final dye concentrations adjusted to 158.5 µM for the fluorochrome HE and 28.87 mM for DHR123 and an optimized incubation time of 30 min in the dark at 2 C. After entering the cell, the fluorochrome HE is oxidized by superoxides and subsequently binds with DNA, whereas DHR123 diffuses into the cell and accumulates in the mitochondria (Benov et al., 1998; Prado et al., 2012). To determine cell-specific concentrations, gated FL1 (505–550 nm) and FL3 (600–645 nm) detection channels were used to analyse the relative concentration of $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ free radicals and H2O2, respectively. The ROS measurements were corrected using blank measurements and normalized to cell size using the forward scatter.

The ROS measurements were taken at two specific time points on the last day of incubation, whereas the photophysiological measurements were taken solely at time point 2 (Fig. 1). The midday measurements, referred to as time point 1, were conducted at the highest light intensity (590 µmol photons m−2 s−1) in the dynamic-light cycle and at the same time under the average light intensity (83 µmol photons m−2 s−1) in the constant-light cycle. The evening measurements, referred to as time point 2, were conducted at the start of the dark period (0 µmol photons m−2 s−1) in both the constant- and dynamic-light cycles.

## 2.5 Statistical analysis

The results are presented as means of the n=4 replicates with a single standard deviation. To identify significant differences between the experimental runs, two-way analysis of variance (ANOVA) tests were performed with a significance level of P≤0.05. The tests were completed using the Minitab Express statistical software (Minitab).

3 Results

## 3.1 Growth and cellular composition

In this study, the growth rates of M. pusilla were affected by both light regime and pCO2 level (Fig. 2). Growth was reduced by at least 50 % in dynamic versus constant light, irrespective of pCO2 level (ANOVA, F(13)=1840.4, p<0.0001). In addition, growth rates significantly increased (>4 %) under elevated pCO2 levels in both light regimes (ANOVA, F(13)=21.9, p=0.0004). The POC and PON quotas were not altered by changes in light regime or pCO2 levels (Table 2). The POC production rates were significantly higher in constant versus dynamic light (ANOVA, F(12)=31.2, p=0.0001), irrespective of the pCO2 level applied. The C:N ratio was not significantly affected by the applied treatments. While Chl a quotas decreased significantly under elevated pCO2 levels (ANOVA, F(13)=26.4, p=0.0002), there was no significant response to the light treatments applied. The applied treatments did not have a significant effect on the C:Chl a ratio. For all these parameters, no significant interactive effects between the applied light and pCO2 conditions could be detected (Fig. 2; Table 2).

Figure 2(a) Growth rate (d−1), (b) POC production (fmol cell d−1), (c) C : N ratio (mol mol−1), (d) chlorophyll a quota (fg cell−1) of Micromonas pusilla under constant light and dynamic light and pCO2 levels of 400 µatm (white) and 1000 µatm (black; n=4; mean ± 1 SD). The letters indicate significant differences between treatments (p<0.05), represented as (a) light and (b) pCO2.

Table 2Growth and cellular composition of M. pusilla (n=4; mean ± 1 SD), including the growth rate, POC production, POC quota, PON quota, chlorophyll a quota, C : N ratio and POC : chlorophyll a ratio. Treatments include constant light and dynamic light and the two pCO2 levels of 400 and 1000 µatm. The letters indicate significant differences between treatments (p<0.05) represented as (a) light and (b) pCO2.

## 3.2 Photophysiological measurements

The FRRf measurements yielded a number of physiological parameters, most of which were significantly affected by the different light and/or pCO2 treatments applied (Table 3). The PSII quantum yield efficiency (FvFm) under dynamic light was significantly higher compared to the constant light treatment (ANOVA, F(13)=88.5, p<0.0001). Even though to a lesser extent, high pCO2 levels also significantly increased FvFm (ANOVA, F(13)=4.8, p=0.0480; Table 3). The connectivity of PSIIs (ρ) was higher under dynamic versus constant light (ANOVA, F(13)=17.6, p=0.0011), while there was no significant effect of pCO2. Similarly, the absorption cross section of PSII photochemistry (σPSII) was significantly higher in dynamic compared to constant light (ANOVA, F(12)=7.0, p<0.0001), irrespective of the applied pCO2 level. In addition, there was significantly less NPQ under dynamic compared to constant light (ANOVA, F(13)=110.3, p<0.0001), but there was no significant pCO2 response in NPQ (p>0.05). Under dynamic light, the PSII reopening rate (τ) was significantly reduced and >5 % lower when compared to constant light (ANOVA, F(13)=18.6, p=0.0008), while τ did not display a significant response to pCO2 (Table 3).

Table 3FRRf-based photophysiological parameters for M. pusilla (n=4; mean ± 1 SD). Displayed is the FvFm (dimensionless), the connectivity between PSIIs (ρ; dimensionless), the absorption cross section of PSII photochemistry (σPSII; ${\mathit{\text{Å}}}^{\mathrm{2}}\cdot {q}^{-\mathrm{1}}$), the non-photochemical quenching (Y(NPQ); dimensionless), the PSII reopening rate (τ; ms), the maximum photosynthetic rate (ETRmax; mol e (mol RCII)−1 s−1), the light harvesting efficiency (α; mol e m2 (mol RCII)−1 (mol photons)−1) and the light saturation constant (Ik; µmol photons m−2 s−1) for both light regimes and pCO2 levels (µatm). The letters indicate significant differences between treatments (p<0.05) represented as (a) light and (b) pCO2.

The model by Webb et al. (1974) was used to estimate PI parameters from the FRRf data. The light saturation index (Ik) and maximum photosynthetic rate (ETRmax) both increased significantly by >10 % under elevated pCO2 levels (ANOVA, F(13)=11.8, p=0.0047 for Ik and F(13)=6.8, p=0.0214 for ETRmax; Table 3). While there was no significant response of Ik to the two light treatments, ETRmax was significantly higher under dynamic light compared to constant light (ANOVA, F(13)=41.2, p<0.0001). The light harvesting efficiency (α) was significantly reduced by high pCO2 versus ambient pCO2 levels (ANOVA, F(13)=9.6, p=0.0084) and significantly higher under dynamic light versus constant light (ANOVA, F(13)=36.0, p<0.0001; Table 3).

## 3.3 ROS levels

The relative concentrations of $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ free radicals and H2O2 were used as an indication of oxidative stress under the applied treatments (Fig. 3). At time point 1 (midday), the production of $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ free radicals was not significantly changed in response to pCO2 levels or light regimes (p>0.05). However, there was a significant increase in H2O2 production in high pCO2 versus ambient pCO2 conditions (ANOVA, F(12)=4.8, p=0.0488), irrespective of the light treatment applied. The applied treatments did not have a significant effect on the ROS levels at time point 2 (Fig. S1, in the Supplement).

Figure 3The relative production of (a) oxygen free radicals (O2) and (b) hydrogen peroxide (H2O2) in Micromonas pusilla under constant light and dynamic light and pCO2 levels of 400 µatm (white) and 1000 µatm (black; n=4; mean ± 1 SD) at time point 1. The letter b indicates a significant difference (p<0.05) between pCO2 treatments.

4 Discussion

## 4.1 Effective acclimation towards dynamic light imposes high metabolic costs

In their natural environment, phytoplankton need to cope with varying light in the upper mixed layer (MacIntyre et al., 2000). Next to variation in insolation, the light fields are critically dependent on the mixed-layer depth, the light attenuation and the vertical mixing rate. In laboratory experiments, however, they are often exposed to an artificially constant light (Köhler et al., 2018). Simulating a dynamic light field to be more representative of an Arctic fjord, we could show that M. pusilla can photoacclimate to these more realistic variations in light availability without showing signs of high-light stress. This is supported by significantly higher PSII quantum yield efficiency (FvFm) under dynamic light (Table 3), which is commonly used as a health indicator of photosynthetic organisms, indicating successful photoacclimation to varying light intensities (Van Leeuwe and Stefels, 2007).

To achieve this photoacclimation, M. pusilla can apparently adjust its PSII physiology to balance photoprotection during high-light periods with sufficient absorption during low-light periods of the dynamic light field. More specifically, there were a number of changes to PSII, including a significant increase in the cross-section size of the antenna in PSII (σPSII), an increase in the connectivity between PSIIs (ρ) and quicker PSII reopening rates (τ; Table 3). An increase in σPSII acts to increase the absorption of light (Suggett et al., 2007), which would have been beneficial within the low-light periods of the dynamic light cycle (Schuback et al., 2017) and is supported by a significant increase in the light harvesting efficiency under low light (α; Table 3). The observed increase in ρ under dynamic light allows for higher flexibility in capturing electrons during low-light phases while at the same time allowing excess excitation energy to be redistributed among PSII centres during high-light phases. This increases energy capture efficiency while protecting the PSII centres from damage through migration of excitation energy between different PSIIs, also termed the lake model (Blankenship, 2014; Trimborn et al., 2014), highlighting that M. pusilla has high potential for photoprotection. Additionally, the higher τ under dynamic versus constant light indicates more efficient drainage of electrons downstream of PSII (Kolber et al., 1998). A faster PSII reopening rate can also compensate for deactivation of functional PSII reaction centres during the high-light periods of the dynamic light field (Behrenfeld et al., 1998). The significantly lower NPQ in combination with higher ETRmax under dynamic versus constant light reflects photoacclimation to a higher light intensity under dynamic light, allowing effective utilization of high excitation energy without initiating high-light stress (Ragni et al., 2008). Consequently, M. pusilla exhibits the physiological plasticity needed to prevent photodamage, which otherwise can disturb the balance between production and scavenging of reactive oxygen species (ROS), causing oxidative stress and accumulation of ROS (Apel and Hirt, 2004). Indeed, dynamic light did not cause increased ROS accumulation in response to the dynamic light field (Fig. 3). Overall, M. pusilla had the capacity to sufficiently acclimate its PSII physiology to deal with dynamic light, displaying photoprotection strategies during high-light phases and upregulated light harvesting during low-light phases.

The described photoacclimation strategies appear to come at a cost, namely lowered energy transfer efficiency to biomass build-up, which is supported by significantly lower growth rates and POC production, despite an increase in FvFm and ETRmax under dynamic compared to constant light (Fig. 2; Table 3). Our findings agree with previous studies, which also found lowered growth under a dynamic light regime (Hoppe et al., 2015; Jin et al., 2013; Mills et al., 2010; Shatwell et al., 2012; Su et al., 2012; Wagner et al., 2006). In previous studies, acclimation to a dynamic light regime has reduced growth rates from 17 % (Hoppe et al., 2015) to 58 % (Boelen et al., 2011), which is comparable to the 47 % reduction in growth rate reported in this study (Table 2; Fig. 2). It thus seems likely that such metabolic costs generally occur and that they are not particularly high in the current study. Changes in light regime strongly influence relationships between photochemistry, carbon fixation and downstream metabolic processes, optimizing light harvesting to sustain growth (Behrenfeld et al., 2008). In view of this, the significant changes to PSII physiology (Table 3) suggest that resources were channelled towards light harvesting rather than protein synthesis and biomass build-up (Talmy et al., 2013). Therefore, it can be concluded that the lower growth rates in dynamic light were caused by the high metabolic costs associated with photoacclimation to the varying light intensities and were not due to photoinhibition. Thus, our results stand in contrast to previous evidence based on which Micromonas was considered as a shade-adapted genus (Lovejoy et al., 2007) as such low-light-adapted species are generally expected to possess limited plasticity in photoacclimative capabilities (Talmy et al., 2013).

## 4.2 Picoeukaryotes benefit from ocean acidification irrespective of the light regime

The low seawater temperatures in the Arctic enhance CO2 solubility and therefore increase OA, from which photosynthetic organisms may benefit due to increased CO2 availability for photosynthesis (AMAP, 2018). This seems true for picoeukaryotes, as in this study M. pusilla showed increased growth rates and photophysiological efficiency under elevated pCO2 (Table 2; Fig. 2). These results are in line with various studies that have reported picoeukaryotes to benefit from increasing pCO2 (Brussaard et al., 2013; Hoppe et al., 2018; Meakin and Wyman, 2011; Newbold et al., 2012; Schaum et al., 2012; Schulz et al., 2017). In the current study, however, there was no observed increase in POC production (Fig. 2) under higher pCO2 levels, which could have been expected assuming lowered costs due to CCM downregulation (Iglesias-Rodriguez et al., 1998; Rost et al., 2008). The observed increase in growth rates nonetheless indicates beneficial OA effects potentially due to reallocation of energy liberated by eased carbon acquisition. Alternatively, the large surface : volume ratio of M. pusilla (cell size of 2–3 µm) may generally lower the need for an active CCM, allowing cells to depend more strongly on diffusive CO2 uptake (Falkowski and Raven, 2013). As the latter is directly linked to the pCO2 level, it could likewise explain the higher growth rates observed under elevated pCO2. In any case, the growth strategy of the investigated strain involves energy allocation into cell division rather than biomass build-up. Whether picoeukaryotes such as M. pusilla benefit from OA due to increased diffusive CO2 uptake, lowered CCM costs or both remains to be tested.

To further explain the increase in growth rate under elevated pCO2, it is essential to look into the upstream physiological parameters. There was a significant increase in ETRmax under OA (Table 3), which indicates an increase in photosynthetic capacity. Previous studies on picoeukaryotes have reported variable results, displaying either no change or an increase in ETRmax in response to OA (Brading et al., 2011; Fu et al., 2007; Kim et al., 2013). In our current study, Ik increased in concert with ETRmax, with increasing pCO2 (Table 3). This Ik-dependent behaviour is known as acclimation to higher light levels in order to optimize balanced growth (Behrenfeld et al., 2008). In the current case, the increase in Ik under OA could indicate that eased carbon acquisition shifted the balance of energy acquisition and its sinks towards saturation at higher irradiances, which fits with the reduced Chl a quota under these conditions (Fig. 2). At the same time, the light harvesting efficiency at low light (α; please note this unit is per photosystem) also decreased in response to OA (Table 3). Such Ik-independent behaviour is influenced by changes in the relative contribution of different sinks of photosynthetic energy, namely carbon fixation, direct use, and ATP generation via cyclic electron transport and other mechanisms (Behrenfeld et al., 2008). Both photoacclimative strategies mimicked acclimation to high light in response to increasing pCO2, which may be a general OA response of phytoplankton (e.g. Hoppe et al., 2015; Rokitta et al., 2012). At the same time, reduced Chl a quotas indicate that such efficient photosystems decrease the need to invest into the total number of them. This potentially balances the reductive pressure on the entire cell, as we did not observe any high-light stress, even during peaks (Fig. 3). Although there was a significant increase in H2O2 concentration under OA relative to ambient pCO2 levels, no change in $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ concentration was observed (Fig. 3). Thus, even if ROS production was enhanced under OA, efficient detoxification mechanisms (e.g. reduction of $\mathrm{•}{\mathrm{O}}_{\mathrm{2}}^{-}$ to H2O2; Asada, 1999) seem to be in place. Additionally, changes to H2O2 concentration have been linked to changes in growth metabolism under non-stressful conditions (Kim et al., 2004), which would fit with the Ik-independent behaviour observed here and suggests that sufficient sinks for the enhanced flow of photosynthetic energy were present. Thus, there is ample evidence that, despite no effect on biomass build-up, elevated pCO2 facilitated carbon acquisition and led to faster and eased photosynthetic energy generation and higher rates of cell division.

The described changes in photoacclimation were not partnered with a significant increase in POC production, despite an increase in growth rate (Fig. 2). These findings contrast with Hoppe et al. (2018), who reported that POC production rates of M. pusilla were generally increased under OA. In this earlier study, however, the Chl a quota of M. pusilla remained relatively constant over a large range of pCO2 levels at two temperatures, so OA effects on the ratio between energy allocated into photosynthesis (i.e. Chl a) and biomass build-up (i.e. POC) in both studies actually agree. Furthermore, if only the pCO2 levels investigated in the current study are considered from Hoppe et al. (2018), varying OA responses (i.e. decreasing vs. increasing for POC production, and constant vs. increasing for the Chl a quota) were observed depending on the applied temperatures. This hints at the well-known fact that even small changes in the environmental conditions can greatly modulate OA responses of phytoplankton (Riebesell and Gattuso, 2015; Rost et al., 2008). In fact, differences between the two studies could also be caused by differences in the average irradiances (approx. 80 vs. 150 µmol photons m−2 s−1). Despite these differences, one should note, however, that high growth rates were obtained under the various OA treatments. As growth rate is the best available fitness indicator for single-strain studies (Schaum and Collins, 2014), our findings are indicative of improved fitness of M. pusilla under OA.

## 4.3M. pusilla's response does not indicate interactions between light regime and pCO2

The interaction between light field and OA has been investigated for the coccolithophore Gephyrocapsa oceanica (Jin et al., 2013) and the Antarctic diatom Chaetoceros debilis (Hoppe et al., 2015). In both studies, the species increased their photochemical performance in response to elevated pCO2 under constant light. Dynamic light fields reversed the positive effect of high pCO2, which was explained by increased high-light stress under OA and a reduction in the energy transfer efficiency from photochemistry to biomass build-up (Hoppe et al., 2015). In the current study, there was no significant interaction between light regime and pCO2 (p>0.05; Figs. 2, 3; Tables 2, 3). These opposing responses could be caused by group- or species-specific differences in carbon acquisition. Diatoms, for example, have highly effective CCMs (Burkhardt et al., 2001), which are energetically expensive (Hopkinson et al., 2011). As CCMs allow cells to efficiently sink energy under sudden high light (Rost et al., 2006), their downregulation in response to high pCO2 can reduce the ability of cells to deal with high-light stress under OA (Hoppe et al., 2015). In contrast to other groups or taxa, which were often found to lose their ability to cope with excess energy under OA and dynamic light (e.g. Gao et al., 2012), M. pusilla maintained effective acclimation without photoinactivation under these conditions. This could be attributed to its size, making it less reliant on CCMs as a mechanism to reduce reductive pressure under high light, as well as to the observed high plasticity in its photophysiological characteristics under dynamic light (Table 3).

In conclusion, the photoacclimation strategies of M. pusilla were optimized for the dynamic light field, and, as this species seems less dependent on CCMs, the previously described interaction between pCO2 and dynamic light (Gao et al., 2012; Hoppe et al., 2015; Jin et al., 2013) was not observed here. This highlights that, depending on their various physiological traits, phytoplankton groups may display different types of interactive responses. It is therefore crucial to understand the underlying physiological mechanisms of observed multi-driver responses in order to judge whether generalizations based on individual studies are feasible or not.

## 4.4 Implications for the future Arctic Ocean

The findings of this study highlight the importance of considering a dynamic light field in laboratory studies. While the interaction between OA and other factors, such as higher temperature, can easily be tested in the lab (Hoppe et al., 2018), light treatments are generally less representative of in situ conditions. The difficulty of measuring and simulating more realistic variations in light has led to the common use of constant light fields, which may substantially alter numerous parameters including growth rates and underestimate the energetic costs of photoacclimation under in situ conditions (Köhler et al., 2018). Therefore, dynamic light fields need to be included when predicting future ecosystem functioning. If the responses of the strain used in this study are representative for this species, M. pusilla can be expected to cope well with a dynamic light field typical of the surface mixed layer (Tables 2, 3). While phytoplankton were often found to suffer from OA under dynamic light (Gao et al., 2012; Hoppe et al., 2015; Jin et al., 2013), M. pusilla benefitted slightly from OA irrespective of the light treatment applied. As beneficial effects by OA were also evident under different temperatures (Hoppe et al., 2018), we can conclude that M. pusilla has a high plasticity towards OA, warming and difference in light regimes, making it well adapted for conditions expected for the future Arctic Ocean. The observed high physiological plasticity, i.e. the ability to adjust physiologically to maintain high growth and/or biomass build-up under all tested scenarios, may thus also explain why picoeukaryotes are often found to dominate mesocosm assemblages under OA (Brussaard et al., 2013; Engel et al., 2008; Schulz et al., 2017).

Global warming is, due to the phenomenon of Arctic amplification (Screen and Simmonds, 2010), a particularly important driver for Arctic phytoplankton. M. pusilla has been shown to synergistically benefit from OA and warming (Hoppe et al., 2018), but the results of this study suggest that future phytoplankton studies should also investigate whether responses differ under dynamic light and determine the mechanisms, metabolic costs and trade-offs associated with interacting physiological processes. Furthermore, warming causes ocean freshening (Peterson et al., 2002) and enhanced stratification that further reduce nutrient availability (Steinacher et al., 2010). Picoeukaryotes may also benefit from these anticipated changes in nutrient supply due to their high surface : volume ratio, allowing for effective nutrient uptake (Li et al., 2009). Additionally, nutrient uptake may be facilitated by lower pH under elevated pCO2 (Bach et al., 2017). Nutrient deficiency was not addressed in this study as the experimental design was aiming to mimic non-limiting nutrient conditions before the spring bloom. Nonetheless, the often limiting nutrient supply in the Arctic sets the trophic status of each region and limits annual productivity (Tremblay et al., 2015) and thus is an important factor to consider in future studies. Changes in the community size structure are biogeochemically important as picoplankton-dominated systems tend to be less efficient with respect to carbon export to depth (Worden et al., 2015). If smaller phytoplankton become more dominant in the Arctic pelagic food web, this may benefit smaller grazers. With these additional steps in the food web, energy transfer efficiency to top predators as well as into the deep ocean will likely decrease (Brussaard et al., 2013). Based on the current study, an increased abundance of M. pusilla under future pCO2 levels can be expected not only for the more stable low-light environments but also for the productive mixed layer in springtime with its dynamic light fields.

Data availability
Data availability.

The dataset for this study is available from the PANGAEA data library, with the identifier https://doi.org/10.1594/PANGAEA.908691 (White et al., 2019).

Supplement
Supplement.

Author contributions
Author contributions.

CJMH designed and supervised the study. EW conducted the research and wrote the paper with contributions from CJMH and BR.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

We are grateful for field support from the 2014/2015 station team of the AWIPEV Base in Ny-Ålesund (Svalbard) as well as for Klara K. E. Wolf's help with strain isolation and maintenance of M. pusilla cultures. Laura Wischnewski, Christine Schallenberg, Tina Brenneis and Marcel Machnik are acknowledged for assistance in the laboratory.

Financial support
Financial support.

The article processing charges for this open-access publication were covered by a research centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by Carol Robinson and reviewed by Douglas Campbell, Lennart Bach and one anonymous referee.

References

AMAP: AMAP Assessment 2018: Arctic Ocean Acidifcation, Arctic Monitoring and Assessment Programme (AMAP), Tromsø, Norway, 187 pp., 2018.

Apel, K. and Hirt, H.: Reactive oxygen species: metabolism, oxidative stress, and signal transduction, Annu. Rev. Plant Biol., 55, 373–399, https://doi.org/10.1146/annurev.arplant.55.031903.141701, 2004.

Arrigo, K. R., van Dijken, G., and Pabi, S.: Impact of a shrinking Arctic ice cover on marine primary production, Geophys. Res. Lett., 35, L19603, https://doi.org/10.1029/2008GL035028, 2008.

Asada, K.: The water-water cycle in chloroplasts: Scavenging of Active Oxygens and Dissipation of Excess Photons, Annu. Rev. Plant Phys., 50, 601–639, https://doi.org/10.1146/annurev.arplant.50.1.601, 1999.

Bach, L. T., Mackinder, L. C., Schulz, K. G., Wheeler, G., Schroeder, D. C., Brownlee, C., and Riebesell, U.: Dissecting the impact of CO2 and pH on the mechanisms of photosynthesis and calcification in the coccolithophore Emiliania huxleyi, New. Phytol., 199, 121–134, https://doi.org/10.1111/nph.12225, 2013.

Bach, L. T., Alvarez-Fernandez, S., Hornick, T., Stuhr, A., and Riebesell, U.: Simulated ocean acidification reveals winners and losers in coastal phytoplankton, PloS One, 12, e0188198, https://doi.org/10.1371/journal.pone.0188198, 2017.

Badger, M. R., Andrews, T. J., Whitney, S. M., Ludwig, M., Yellowlees, D. C., Leggat, W., and Price, G. D.: The diversity and coevolution of Rubisco, plastids, pyrenoids, and chloroplast-based CO2-concentrating mechanisms in algae, Can. J. Bot., 76, 1052–1071, https://doi.org/10.1139/b98-074, 1998.

Bates, N. R. and Mathis, J. T.: The Arctic Ocean marine carbon cycle: evaluation of air-sea CO2 exchanges, ocean acidification impacts and potential feedbacks, Biogeosciences, 6, 2433–2459, https://doi.org/10.5194/bg-6-2433-2009, 2009.

Behrenfeld, M. J., Prasil, O., Kolber, Z. S., Babin, M., and Falkowski, P. G.: Compensatory changes in photosystem II electron turnover rates protect photosynthesis from photoinhibition, Photosynth. Res., 58, 259–268, https://doi.org/10.1023/A:1006138630573, 1998.

Behrenfeld, M. J., Halsey, K. H., and Milligan, A. J.: Evolved physiological responses of phytoplankton to their integrated growth environment, Philos. T. R. Soc. Lon. B, 363, 2687–2703, https://doi.org/10.1098/rstb.2008.0019, 2008.

Benov, L., Sztejnberg, L., and Fridovich, I.: Critical evaluation of the use of hydroethidine as a measure of superoxide anion radical, Free Radical Bio. Med., 25, 826–831, https://doi.org/10.1016/S0891-5849(98)00163-4, 1998.

Blankenship, R. E.: Molecular mechanisms of photosynthesis, 2nd Edn., Wiley Blackwell, USA, 2014.

Boelen, P., van de Poll, W. H., van der Strate, H. J., Neven, I. A., Beardall, J., and Buma, A. G.: Neither elevated nor reduced CO2 affects the photophysiological performance of the marine Antarctic diatom Chaetoceros brevis, J. Exp. Mar. Biol. Ecol., 406, 38–45, https://doi.org/10.1016/j.jembe.2011.06.012, 2011.

Brading, P., Warner, M. E., Davey, P., Smith, D. J., Achterberg, E. P., and Suggett, D. J.: Differential effects of ocean acidification on growth and photosynthesis among phylotypes of Symbiodinium (Dinophyceae), Limnol. Oceanogr., 56, 927–938, https://doi.org/10.4319/lo.2011.56.3.0927, 2011.

Brewer, P. G., Bradshaw, A., and Williams, R.: Measurements of total carbon dioxide and alkalinity in the North Atlantic Ocean in 1981, in: The changing carbon cycle, Springer, New York, NY, 348–370, https://doi.org/10.1007/978-1-4757-1915-4_18, 1986.

Brussaard, C. P. D., Noordeloos, A. A. M., Witte, H., Collenteur, M. C. J., Schulz, K., Ludwig, A., and Riebesell, U.: Arctic microbial community dynamics influenced by elevated CO2 levels, Biogeosciences, 10, 719–731, https://doi.org/10.5194/bg-10-719-2013, 2013.

Burkhardt, S., Amoroso, G., Riebesell, U., and Sültemeyer, D.: CO2 and ${\mathrm{HCO}}_{\mathrm{3}}^{-}$ uptake in marine diatoms acclimated to different CO2 concentrations, Limnol. Oceanogr., 46, 1378–1391, https://doi.org/10.4319/lo.2001.46.6.1378, 2001.

Denman, K. L. and Gargett, A. E.: Time and space scales of vertical mixing and advection of phytoplankton in the upper ocean, Limnol. Oceanogr., 28, 801–815, https://doi.org/10.4319/lo.1983.28.5.0801, 1983.

Dickson, A. G. and Millero, F. J.: A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media, Deep-Sea Res., 34, 1733–1743, https://doi.org/10.1016/0198-0149(87)90021-5, 1987.

Dickson, A. G., Sabine, C. L., and Christian, J. R.: Guide to best practices for ocean CO2 measurements, North Pacific Marine Science Organization, Sidney, British Columbia, 2007.

Engel, A., Schulz, K. G., Riebesell, U., Bellerby, R., Delille, B., and Schartau, M.: Effects of CO2 on particle size distribution and phytoplankton abundance during a mesocosm bloom experiment (PeECE II), Biogeosciences, 5, 509–521, https://doi.org/10.5194/bg-5-509-2008, 2008.

Falkowski, P. G.: Light-shade adaptation in marine phytoplankton, in: Primary productivity in the sea, Vol. 19, Springer, Boston, USA, 1980.

Falkowski, P. G. and Raven, J. A.: Aquatic photosynthesis, 2nd Edn., Princeton University Press, USA, 2013.

Fu, F.-X., Warner, M. E., Zhang, Y., Feng, Y., and Hutchins, D. A.: Effects of Increased Temperature and CO2 on Photosynthesis, Growth, and Elemental Ratios in Marine Synechococcus and Prochlorococcus (Cyanobacteria), J. Phycol., 43, 485–496, https://doi.org/10.1111/j.1529-8817.2007.00355.x, 2007.

Gao, K., Helbling, E. W., Häder, D. P., and Hutchins, D. A.: Responses of marine primary producers to interactions between ocean acidification, solar radiation, and warming, Mar. Ecol.-Prog. Ser., 470, 167–189, https://doi.org/10.3354/meps10043, 2012.

Granskog, M. A., Stedmon, C. A., Dodd, P. A., Amon, R. M. W., Pavlov, A. K., de Steur, L., and Hansen, E.: Characteristics of colored dissolved organic matter (CDOM) in the Arctic outflow in the Fram Strait: Assessing the changes and fate of terrigenous CDOM in the Arctic Ocean, J. Geophys. Res.-Oceans, 117, C12021, https://doi.org/10.1029/2012jc008075, 2012.

Guillard, R. R. and Ryther, J. H.: Studies of marine planktonic diatoms. I. Cyclotella nana Hustedt, and Detonula confervacea (cleve) Gran, Can. J. Microbiol., 8, 229–239, https://doi.org/10.1139/m62-029, 1962.

Hopkinson, B. M., Dupont, C. L., Allen, A. E., and Morel, F. M.: Efficiency of the CO2-concentrating mechanism of diatoms, P. Natl. Acad. Sci. USA, 108, 3830–3837, https://doi.org/10.1073/pnas.1018062108, 2011.

Hoppe, C. J., Holtz, L. M., Trimborn, S., and Rost, B.: Ocean acidification decreases the light-use efficiency in an Antarctic diatom under dynamic but not constant light, New Phytol., 207, 159–171, https://doi.org/10.1111/nph.13334, 2015.

Hoppe, C. J. M., Langer, G., Rokitta, S. D., Wolf-Gladrow, D. A., and Rost, B.: Implications of observed inconsistencies in carbonate chemistry measurements for ocean acidification studies, Biogeosciences, 9, 2401–2405, https://doi.org/10.5194/bg-9-2401-2012, 2012.

Hoppe, C. J., Schuback, N., Semeniuk, D. M., Maldonado, M. T., and Rost, B.: Functional Redundancy Facilitates Resilience of Subarctic Phytoplankton Assemblages toward Ocean Acidification and High Irradiance, Front. Mar. Sci., 4, 229, https://doi.org/10.3389/fmars.2017.00229, 2017.

Hoppe, C. J. M., Flintrop, C. M., and Rost, B.: The Arctic picoeukaryote Micromonas pusilla benefits synergistically from warming and ocean acidification, Biogeosciences, 15, 4353–4365, https://doi.org/10.5194/bg-15-4353-2018, 2018.

Houghton, J., Ding, Y., and Griggs, D.: Climate Change: Scientific Basis, IPCC TAR Working Group 1, The Press Syndicate of Cambridge University, UK, 2001.

Hu, A. and Bates, S. C.: Internal climate variability and projected future regional steric and dynamic sea level rise, Nat. Commun., 9, 1068, https://doi.org/10.1038/s41467-018-03474-8, 2018.

Huisman, J.: Population dynamics of light-limited phytoplankton: microcosm experiments, Ecology, 80, 202–210, https://doi.org/10.1890/0012-9658(1999)080[0202:PDOLLP]2.0.CO;2, 1999.

Iglesias-Rodriguez, M., Nimer, N., and Merrett, M.: Carbon dioxide-concentrating mechanism and the development of extracellular carbonic anhydrase in the marine picoeukaryote Micromonas pusilla, New Phytol., 140, 685–690, https://doi.org/10.1046/j.1469-8137.1998.00309.x, 1998.

Jin, P., Gao, K., Villafañe, V. E., Campbell, D. A., and Helbling, W.: Ocean acidification alters the photosynthetic responses of a coccolithophorid to fluctuating UV and visible radiation, Plant Physiol., 162, 2084–2094, https://doi.org/10.1104/pp.113.219543, 2013.

Kim, D., Watanabe, M., Nakayasu, Y., and Kohata, K.: Production of superoxide anion and hydrogen peroxide associated with cell growth of Chattonella antiqua, Aquat. Microb. Ecol., 35, 57–64, https://doi.org/10.3354/ame035057, 2004.

Kim, J.-H., Kim, K. Y., Kang, E. J., Lee, K., Kim, J.-M., Park, K.-T., Shin, K., Hyun, B., and Jeong, H. J.: Enhancement of photosynthetic carbon assimilation efficiency by phytoplankton in the future coastal ocean, Biogeosciences, 10, 7525–7535, https://doi.org/10.5194/bg-10-7525-2013, 2013.

Klughammer, C. and Schreiber, U.: Complementary PS II quantum yields calculated from simple fluorescence parameters measured by PAM fluorometry and the Saturation Pulse method, PAM Appl. Notes, 1, 27–35, 2008.

Knap, A., Michaels, A., Close, A., Ducklow, H., and Dickson, A.: Protocols for the Joint Global Ocean Flux Study (JGOFS) Core Measurements, UNESCO, Paris, France, 170 pp., 1996.

Köhler, J., Wang, L., Guislain, A., and Shatwell, T.: Influence of vertical mixing on light-dependency of phytoplankton growth, Limnol. Oceanogr., 63, 1156–1167, https://doi.org/10.1002/lno.10761, 2018.

Kolber, Z. S., Prášil, O., and Falkowski, P. G.: Measurements of variable chlorophyll fluorescence using fast repetition rate techniques: defining methodology and experimental protocols, BBA-Bioenergetics, 1367, 88–106, https://doi.org/10.1016/S0005-2728(98)00135-2, 1998.

Kramer, D. M., Johnson, G., Kiirats, O., and Edwards, G. E.: New fluorescence parameters for the determination of QA redox state and excitation energy fluxes, Photosynth. Res., 79, 209–218, https://doi.org/10.1023/B:PRES.0000015391.99477.0d, 2004.

Li, W. K., McLaughlin, F. A., Lovejoy, C., and Carmack, E. C.: Smallest algae thrive as the Arctic Ocean freshens, Science, 326, 539–539, https://doi.org/10.1126/science.1179798, 2009.

Litchman, E.: Growth rates of phytoplankton under fluctuating light, Freshwater Biol., 44, 223–235, https://doi.org/10.1046/j.1365-2427.2000.00559.x, 2000.

Litchman, E., Klausmeier, C. A., and Bossard, P.: Phytoplankton nutrient competition under dynamic light regimes, Limnol. Oceanogr., 49, 1457–1462, https://doi.org/10.4319/lo.2004.49.4_part_2.1457, 2004.

Lohbeck, K. T., Riebesell, U., and Reusch, T. B.: Adaptive evolution of a key phytoplankton species to ocean acidification, Nat. Geosci., 5, 346, https://doi.org/10.1038/ngeo1441, 2012.

Lovejoy, C., Vincent, W. F., Bonilla, S., Roy, S., Martineau, M. J., Terrado, R., Potvin, M., Massana, R., and Pedrós-Alió, C.: Distribution, phylogeny, and growth of cold-adapted picoprasinophytes in Arctic Seas, J. Phycol., 43, 78–89, https://doi.org/10.1111/j.1529-8817.2006.00310.x, 2007.

Maat, D. S., Crawfurd, K. J., Timmermans, K. R., and Brussaard, C. P.: Elevated pCO2 and phosphate limitation favor Micromonas pusilla through stimulated growth and reduced viral impact, Appl. Environ. Microb., 80, 3119–3127, https://doi.org/10.1128/AEM.03639-13., 2014.

MacIntyre, H. L., Kana, T. M., and Geider, R. J.: The effect of water motion on short-term rates of photosynthesis by marine phytoplankton, Trends Plant Sci., 5, 12–17, https://doi.org/10.1016/S1360-1385(99)01504-6, 2000.

Marquardt, M., Vader, A., Stübner, E. I., Reigstad, M., and Gabrielsen, T. M.: Strong seasonality of marine microbial eukaryotes in a high-arctic fjord (Isfjorden, in West Spitsbergen, Norway), Appl. Environ. Microb., 82, 1868–1880, https://doi.org/10.1128/AEM.03208-15, 2016.

McKie-Krisberg, Z. M. and Sanders, R. W.: Phagotrophy by the picoeukaryotic green alga Micromonas: implications for Arctic Oceans, ISME J., 8, 1953–1961, https://doi.org/10.1038/ismej.2014.16, 2014.

Meakin, N. G. and Wyman, M.: Rapid shifts in picoeukaryote community structure in response to ocean acidification, ISME J., 5, 1397–1405, https://doi.org/10.1038/ismej.2011.18, 2011.

Mehrbach, C., Culberson, C., Hawley, J., and Pytkowicx, R.: Measurement of the apparent dissociation constants of carbonic acid in seawater at atmospheric pressure, Limnol. Oceanogr., 18, 897–907, https://doi.org/10.4319/lo.1973.18.6.0897, 1973.

Miller, G. H., Alley, R. B., Brigham-Grette, J., Fitzpatrick, J. J., Polyak, L., Serreze, M. C., and White, J. W.: Arctic amplification: can the past constrain the future?, Quaternary Sci. Rev., 29, 1779–1790, https://doi.org/10.1016/j.quascirev.2010.02.008, 2010.

Mills, M. M., Kropuenske, L. R., Van Dijken, G. L., Alderkamp, A. C., Berg, G. M., Robinson, D. H., Welschmeyer, N. A., and Arrigo, K. R.: Photophysiology in two southern ocean phytoplankton taxa: Photosynthesis of Phaeocystis Antarctica (prymnesiophyceae) and Fragilariopsis Cylindrus (bacillariophyceae) under simulated mixed-layer irradiance, J. Phycol., 46, 1114–1127, https://doi.org/10.1111/j.1529-8817.2010.00923.x, 2010.

Newbold, L. K., Oliver, A. E., Booth, T., Tiwari, B., Desantis, T., Maguire, M., Andersen, G., van der Gast, C. J., and Whiteley, A. S.: The response of marine picoplankton to ocean acidification, Environ. Microbiol., 14, 2293–2307, https://doi.org/10.1111/j.1462-2920.2012.02762.x, 2012.

Nöthig, E.-M., Bracher, A., Engel, A., Metfies, K., Niehoff, B., Peeken, I., Bauerfeind, E., Cherkasheva, A., Gäbler-Schwarz, S., and Hardge, K.: Summertime plankton ecology in Fram Strait – a compilation of long-and short-term observations, Polar Res., 34, 23349, https://doi.org/10.3402/polar.v34.23349, 2015.

Oxborough, K. and Baker, N. R.: Resolving chlorophyll a fluorescence images of photosynthetic efficiency into photochemical and non-photochemical components – calculation of qP and ${F}_{\mathrm{v}}^{\prime }/{F}_{\mathrm{m}}^{\prime }$; without measuring ${F}_{\mathrm{o}}^{\prime }$, Photosynth. Res., 54, 135–142, https://doi.org/10.1023/A:1005936823310, 1997.

Oxborough, K., Moore, C. M., Suggett, D. J., Lawson, T., Chan, H. G., and Geider, R. J.: Direct estimation of functional PSII reaction center concentration and PSII electron flux on a volume basis: a new approach to the analysis of Fast Repetition Rate fluorometry (FRRf) data, Limnol. Oceanogr.-Meth., 10, 142–154, https://doi.org/10.4319/lom.2012.10.142, 2012.

Pachauri, R. K., Allen, M. R., Barros, V. R., Broome, J., Cramer, W., Christ, R., Church, J. A., Clarke, L., Dahe, Q., and Dasgupta, P.: Climate change 2014: synthesis report. Contribution of Working Groups I, II and III to the fifth assessment report of the Intergovernmental Panel on Climate Change, IPCC, Geneva, Switzerland, 2014.

Palmer, M. A., Van Dijken, G. L., Mitchell, B. G., Seegers, B. J., Lowry, K. E., Mills, M. M., and Arrigo, K. R.: Light and nutrient control of photosynthesis in natural phytoplankton populations from the Chukchi and Beaufort seas, Arctic Ocean, Limnol. Oceanogr., 58, 2185–2205, https://doi.org/10.4319/lo.2013.58.6.2185, 2013.

Peterson, B. J., Holmes, R. M., McClelland, J. W., Vorosmarty, C. J., Lammers, R. B., Shiklomanov, A. I., Shiklomanov, I. A., and Rahmstorf, S.: Increasing river discharge to the Arctic Ocean, Science, 298, 2171–2173, https://doi.org/10.1126/science.1077445, 2002.

Pierrot, D., Lewis, E., and Wallace, D.: MS Excel program developed for CO2 system calculations, ORNL/CDIAC-105a. Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy, Oak Ridge, Tennessee, 2006.

Pörtner, H.-O., Karl, D. M., Boyd, P. W., Cheung, W., Lluch-Cota, S. E., Nojiri, Y., Schmidt, D. N., Zavialov, P. O., Alheit, J., and Aristegui, J.: Ocean systems, in: Climate change 2014: impacts, adaptation, and vulnerability. Part A: global and sectoral aspects, Contribution of working group II to the fifth assessment report of the intergovernmental panel on climate change, Cambridge University Press, 411–484, Cambridge, UK and NY, USA, 2014.

Prado, R., Rioboo, C., Herrero, C., Suarez-Bregua, P., and Cid, A.: Flow cytometric analysis to evaluate physiological alterations in herbicide-exposed Chlamydomonas moewusii cells, Ecotoxicology, 21, 409–420, https://doi.org/10.1007/s10646-011-0801-3, 2012.

Ragni, M., Airs, R. L., Leonardos, N., and Geider, R. J.: Photoinhibition of PSII in Emiliania Huxleyi (Haptophyta) under High Light Stress: The Roles of Photoacclimation, Photoprotection, and Photorepair, J. Phycol., 44, 670–683, https://doi.org/10.1111/j.1529-8817.2008.00524.x, 2008.

Reinfelder, J. R.: Carbon concentrating mechanisms in eukaryotic marine phytoplankton, Annu. Rev. Mar. Sci., 3, 291–315, https://doi.org/10.1146/annurev-marine-120709-142720, 2011.

Riebesell, U. and Gattuso, J.-P.: Lessons learned from ocean acidification research, Nat. Clim. Change, 5, 12–14, 2015.

Riebesell, U. and Tortell, P. D.: Effects of ocean acidification on pelagic organisms and ecosystems, Ocean acidification, Oxford University Press, Oxford, UK, 2011.

Rokitta, S. D., John, U., and Rost, B.: Ocean acidification affects redox-balance and ion-homeostasis in the life-cycle stages of Emiliania huxleyi, PLoS One, 7, 1–10, e52212, https://doi.org/10.1371/journal.pone.0052212, 2012.

Rost, B., Riebesell, U., and Sültemeyer, D.: Carbon acquisition of marine phytoplankton: effect of photoperiod length, Limnol. Oceanogr., 51, 12–20, https://doi.org/10.4319/lo.2006.51.1.0012, 2006.

Rost, B., Zondervan, I., and Wolf-Gladrow, D.: Sensitivity of phytoplankton to future changes in ocean carbonate chemistry: current knowledge, contradictions and research directions, Mar. Ecol.-Prog. Ser., 373, 227–237, https://doi.org/10.3354/meps07776, 2008.

Schaum, C. E. and Collins, S.: Plasticity predicts evolution in a marine alga, Proc. Biol. Sci., 281, 20141486, https://doi.org/10.1098/rspb.2014.1486, 2014.

Schaum, E., Rost, B., Millar, A. J., and Collins, S.: Variation in plastic responses of a globally distributed picoplankton species to ocean acidification, Nat. Clim. Change, 3, 298–302, https://doi.org/10.1038/nclimate1774, 2012.

Schuback, N., Hoppe, C. J., Tremblay, J. É., Maldonado, M. T., and Tortell, P. D.: Primary productivity and the coupling of photosynthetic electron transport and carbon fixation in the Arctic Ocean, Limnol. Oceanogr., 62, 898–921, https://doi.org/10.1002/lno.10475, 2017.

Schulz, K. G., Bach, L. T., Bellerby, R. G., Bermúdez, R., Büdenbender, J., Boxhammer, T., Czerny, J., Engel, A., Ludwig, A., and Meyerhöfer, M.: Phytoplankton blooms at increasing levels of atmospheric carbon dioxide: Experimental evidence for negative effects on prymnesiophytes and positive on small picoeukaryotes, Front. Mar. Sci., 4, 64, https://doi.org/10.3389/fmars.2017.00064, 2017.

Screen, J. A. and Simmonds, I.: The central role of diminishing sea ice in recent Arctic temperature amplification, Nature, 464, 1334, https://doi.org/10.1038/nature09051, 2010.

Shatwell, T., Nicklisch, A., and Köhler, J.: Temperature and photoperiod effects on phytoplankton growing under simulated mixed layer light fluctuations, Limnol. Oceanogr., 57, 541–553, https://doi.org/10.4319/lo.2012.57.2.0541, 2012.

Steinacher, M., Joos, F., Frölicher, T. L., Bopp, L., Cadule, P., Cocco, V., Doney, S. C., Gehlen, M., Lindsay, K., Moore, J. K., Schneider, B., and Segschneider, J.: Projected 21st century decrease in marine productivity: a multi-model analysis, Biogeosciences, 7, 979–1005, https://doi.org/10.5194/bg-7-979-2010, 2010.

Stoll, M., Bakker, K., Nobbe, G., and Haese, R.: Continuous-flow analysis of dissolved inorganic carbon content in seawater, Anal. Chem., 73, 4111–4116, https://doi.org/10.1021/ac010303r, 2001.

Su, W., Jakob, T., and Wilhelm, C.: The Impact of Nonphotochemical Quenching of Fluorescence on the Photon Balance in Diatoms under Dynamic Light Conditions, J. Phycol., 48, 336–346, https://doi.org/10.1111/j.1529-8817.2012.01128.x, 2012.

Suggett, D. J., Le Floc'H, E., Harris, G. N., Leonardos, N., and Geider, R. J.: Different strategies of photoacclimation by two strains of Emiliania huxleyi (Haptophyta), J. Phycol., 43, 1209–1222, https://doi.org/10.1111/j.1529-8817.2007.00406.x, 2007.

Talmy, D., Blackford, J., Hardman-Mountford, N. J., Dumbrell, A. J., and Geider, R. J.: An optimality model of photoadaptation in contrasting aquatic light regimes, Limnol. Oceanogr., 58, 1802–1818, https://doi.org/10.4319/lo.2013.58.5.1802, 2013.

Tremblay, J.-É., Anderson, L. G., Matrai, P., Coupel, P., Bélanger, S., Michel, C., and Reigstad, M.: Global and regional drivers of nutrient supply, primary production and CO2 drawdown in the changing Arctic Ocean, Prog. Oceanogr., 139, 171–196, https://doi.org/10.1016/j.pocean.2015.08.009, 2015.

Trenberth, K., Jones, P., Ambenje, P., Bojariu, R., Easterling, D., Klein Tank, A., Parker, D., Rahimzadeh, F., Renwick, J., and Rusticucci, M.: Observations: surface and atmospheric climate change, in: Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, 235–336, Cambridge, UK and NY, USA, 2007.

Trimborn, S., Thoms, S., Petrou, K., Kranz, S. A., and Rost, B.: Photophysiological responses of Southern Ocean phytoplankton to changes in CO2 concentrations: short-term versus acclimation effects, J. Exp. Mar. Biol. Ecol., 451, 44–54, https://doi.org/10.1016/j.jembe.2013.11.001, 2014.

Uppström, L. R.: The boron/chlorinity ratio of deep-sea water from the Pacific Ocean, Deep-Sea Res. Oceanogr. Abstr., 21, 161–162, https://doi.org/10.1016/0011-7471(74)90074-6, 1974.

Van Leeuwe, M. and Stefels, J.: Photosynthetic responses in Phaeocystis antarctica towards varying light and iron conditions, in: Phaeocystis, major link in the biogeochemical cycling of climate-relevant elements, Springer, 61–70, https://doi.org/10.1007/978-1-4020-6214-8_6, 2007.

Wagner, H., Jakob, T., and Wilhelm, C.: Balancing the energy flow from captured light to biomass under fluctuating light conditions, New Phytol., 169, 95–108, https://doi.org/10.1111/j.1469-8137.2005.01550.x, 2006.

Wassmann, P. and Reigstad, M.: Future Arctic Ocean seasonal ice zones and implications for pelagic-benthic coupling, Oceanography, 24, 220–231, https://doi.org/10.5670/oceanog.2011.74, 2011.

Webb, W. L., Newton, M., and Starr, D.: Carbon dioxide exchange of Alnus rubra: A mathematical model, Oecologia, 17, 281–291, https://doi.org/10.1007/BF00345747, 1974.

White, E., Hoppe, C. J., and Rost, B.: The Arctic picoeukaryote Micromonas pusilla benefits from Ocean Acidification under constant and dynamic light, PANGAEA, https://doi.org/10.1594/PANGAEA.908691, 2019.

Worden, A. Z., Follows, M. J., Giovannoni, S. J., Wilken, S., Zimmerman, A. E., and Keeling, P. J.: Rethinking the marine carbon cycle: factoring in the multifarious lifestyles of microbes, Science, 347, 1257594, https://doi.org/10.1126/science.1257594, 2015.

Xu, K., Grant-Burt, J. L., Donaher, N., and Campbell, D. A.: Connectivity among Photosystem II centers in phytoplankters: Patterns and responses, BBA-Bioenergetics, 1858, 459–474, https://doi.org/10.1016/j.bbabio.2017.03.003, 2017.