
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39231958
51636
10.1038/s41467-024-51636-8
Article
Photosynthetic light requirement near the theoretical minimum detected in Arctic microalgae
http://orcid.org/0000-0002-2509-0546
Hoppe Clara J. M. Clara.Hoppe@awi.de

1
http://orcid.org/0000-0001-8536-6877
Fuchs Niels 2
http://orcid.org/0000-0003-0365-5654
Notz Dirk 2
http://orcid.org/0000-0001-9196-6289
Anderson Philip 3
Assmy Philipp 4
http://orcid.org/0000-0003-0900-5679
Berge Jørgen 5
http://orcid.org/0000-0001-8388-4945
Bratbak Gunnar 6
Guillou Gaël 7
Kraberg Alexandra 1
http://orcid.org/0000-0001-6927-5537
Larsen Aud 8
Lebreton Benoit 7
Leu Eva 9
http://orcid.org/0000-0003-4276-4781
Lucassen Magnus 1
http://orcid.org/0000-0001-5405-052X
Müller Oliver 6
http://orcid.org/0000-0001-7939-4011
Oziel Laurent 1
http://orcid.org/0000-0001-5452-5505
Rost Björn 110
http://orcid.org/0000-0001-6879-7800
Schartmüller Bernhard 5
http://orcid.org/0000-0002-8283-656X
Torstensson Anders 1112
Wloka Jonas 13
1 https://ror.org/032e6b942 grid.10894.34 0000 0001 1033 7684 Alfred-Wegener-Institut—Helmholtz-Zentrum für Polar- und Meeresforschung, Bremerhaven, Germany
2 https://ror.org/00g30e956 grid.9026.d 0000 0001 2287 2617 Center for Earth System Sustainability, Institute of Oceanography, University of Hamburg, Hamburg, Germany
3 https://ror.org/04ke6ht85 grid.410415.5 0000 0000 9388 4992 Scottish Association for Marine Science, Oban, Scotland
4 grid.417991.3 0000 0004 7704 0318 Norwegian Polar Institute, Fram Centre, Tromsø, Norway
5 https://ror.org/00wge5k78 grid.10919.30 0000 0001 2259 5234 UiT The Arctic University of Norway, Tromsø, Norway
6 https://ror.org/03zga2b32 grid.7914.b 0000 0004 1936 7443 University of Bergen, Bergen, Norway
7 https://ror.org/04mv1z119 grid.11698.37 0000 0001 2169 7335 Joint Research Unit Littoral, Environment and Societies (CNRS—University of La Rochelle), La Rochelle, France
8 https://ror.org/02gagpf75 grid.509009.5 NORCE Norwegian Research Centre, Bergen, Norway
9 https://ror.org/03nrps502 grid.510420.2 0000 0004 7554 3448 Akvaplan-niva, Tromsø, Norway
10 https://ror.org/04ers2y35 grid.7704.4 0000 0001 2297 4381 Faculty of Biology/Chemistry, University Bremen, Bremen, Germany
11 https://ror.org/02yy8x990 grid.6341.0 0000 0000 8578 2742 Department of Aquatic Sciences and Assessment, Swedish University of Agricultural Sciences, Uppsala, Sweden
12 https://ror.org/00hgzve81 grid.6057.4 0000 0001 0289 1343 Swedish Meteorological and Hydrological Institute, Community Planning Services—Oceanography, Västra Frölunda, Sweden
13 grid.5949.1 0000 0001 2172 9288 Independent Researcher, Bremen, Germany
4 9 2024
4 9 2024
2024
15 738521 11 2023
14 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Photosynthesis is one of the most important biological processes on Earth, providing the main source of bioavailable energy, carbon, and oxygen via the use of sunlight. Despite this importance, the minimum light level sustaining photosynthesis and net growth of primary producers in the global ocean is still unknown. Here, we present measurements from the MOSAiC field campaign in the central Arctic Ocean that reveal the resumption of photosynthetic growth and algal biomass buildup under the ice pack at a daily average irradiance of not more than 0.04 ± 0.02 µmol photons m−2 s−1 in late March. This is at least one order of magnitude lower than previous estimates (0.3–5 µmol photons m−2 s−1) and near the theoretical minimum light requirement of photosynthesis (0.01 µmol photons m−2 s−1). Our findings are based on measurements of the temporal development of the under-ice light field and concurrent measurements of both chlorophyll a concentrations and potential net primary production underneath the sea ice at 86 °N. Such low light requirements suggest that euphotic zones where photosynthesis can occur in the world’s oceans may extend further in depth and time, with major implications for global productivity estimates.

This study shows that photosynthesis in the Arctic Ocean can occur at lower light levels than previously thought. This finding is based on unique measurements made on microalgae underneath sea ice in the high Arctic at the end of the polar night.

Subject terms

Microbial biooceanography
Microbial ecology
Marine biology
https://doi.org/10.13039/501100002347 Bundesministerium für Bildung und Forschung (Federal Ministry of Education and Research) 03F0867A Hoppe Clara J. M. https://doi.org/10.13039/501100001659 Deutsche Forschungsgemeinschaft (German Research Foundation) 390683824 Hoppe Clara J. M. https://doi.org/10.13039/501100005416 Norges Forskningsråd (Research Council of Norway) 280292 Hoppe Clara J. M. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The distribution of life on Earth is critically linked to environmental conditions suitable for photosynthesis. In the ocean, photosynthesis can only occur in the uppermost euphotic layer where sufficient light is available. Therefore, the extent of the euphotic zone, where photosynthesis exceeds respiration and other loss terms1, is a key control for oceanic net primary production (NPP) and major biogeochemical cycles2. Despite this importance, the minimum light requirement for photosynthetic growth, and the according maximum depth of the euphotic zone, is still not well constrained. In this study, we address this gap by presenting complementary field measurements of under-ice irradiances and phototrophic biomass from the high Arctic that show NPP to be possible at extremely low light levels.

Traditionally, the lower limit of the euphotic zone is defined as the depth at which light levels are 1% of their surface value1, which, in absolute terms, can be as high as 20 µmol photons m−2 s−1 for midday conditions at the equator with a typical surface value of 2000 µmol photons m−2 s−1. However, it has long been known that this threshold overestimates the minimum light requirement for NPP3, implying that the euphotic zone actually extends to greater depths. For the absolute lower light requirement for photosynthetic growth, theoretical considerations suggest a light level of around 0.01 µmol photons m−2 s−1 4–6. Below this threshold, photosynthetic carbon fixation is considered impossible due to biochemical constraints6. However, it is unclear how far from this threshold photosynthesis indeed occurs in nature. In the ocean, NPP has been observed at light levels as low as 0.3–5 µmol photons m−2 s−1 5,7–10, i.e., 30–500 times higher than the theoretical minimum threshold. For ice algae, thresholds of 0.17 and 0.36 µmol photons m−2 s−1 have been reported4,11, still exceeding the theoretically derived minimum threshold by far. While differences between these measurements may seem small in absolute terms, the light levels of large fractions of the world’s oceans are in these ranges, at least for some part of the year. Therefore, knowledge of the lower light limit at which NPP can occur is a prerequisite to realistically estimate the temporal and spatial distribution of primary production and concurrent ecosystem dynamics.

Here, we address this knowledge gap by analyzing field measurements from the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC), during which an interdisciplinary team of researchers on RV Polarstern drifted across the Arctic Ocean with the sea ice for one year (September 2019–October 202012,13). Such a year-round survey allowed us to study the resumption of photosynthetic activity after the hiatus of the polar night, as sunlight returned to the central Arctic Ocean (at 88–84°N) during the period February to April 2020 (Fig. S1). For our study, data were collected at different locations on the pack ice around the ship, while the ice drifted over the mixed layer of the ocean (Figs. 1, S2). We used three different methods to determine the potential onset of photosynthetic activity during that period: measurements of potential primary production, determination of chlorophyll a (Chl-a) concentrations, and direct cell counts of photosynthetic microalgae. In the following, we discuss the results of these measurements as a function of time before linking them to the related light levels.Fig. 1 Schematic illustration of sampling sites and depths.

MOSAiC sampling sites and depths around RV Polarstern during spring 2020 relevant to this study are illustrated in a schematic way. Biological parameters from the upper mixed layer (circular arrows) were collected from 20 m depth via rosettes deployed through holes in the ice at Ocean City and next to RV Polarstern (see refs. 13,24), as well as the ship’s underway system at 11 m depth. Sea ice cores for sea ice Chlorophyll a and Net Primary Production P were collected at the first-year ice coring site. Light measurements were collected with OptiCALs (Optical Chain And Logger) gg, hh, and ee down to 50 m water depth as well as the lightharp, the latter inside the ice column only. Exact locations of the different sampling sites are illustrated in Fig. S2. Please note that the ice was drifting about six times faster than the underlying water column23 so that specific locations on the pack ice do not represent different sampling locations in the ocean.

Results and Discussion

The most sensitive method we employed to determine the potential onset of photosynthesis underneath the ice is measuring the potential 14C productivity (termed NPP, in the following) of the ambient phytoplankton assemblages for 24 h using 14C-based volume-specific carbon fixation rates. Surface water samples were exposed to reference conditions of a temperature of 1 °C and an irradiance of 10 µmol photons m−2 s−1 in the laboratory. NPP rates significantly above dark controls were measured as early as in January and February, indicating physiologically active overwintering assemblages that can resume photosynthesis as soon as they are exposed to light (Fig. S3). Surprisingly, the onset of an exponential increase of potential NPP occurred in the mixed layer before March 14th (Fig. 2A). The lack of trends in carbon fixation to Chl-a (Fig. S4, Mann–Kendall Trend test) and particulate organic carbon (POC) to Chl-a ratios (Fig. S5, Mann–Kendall Trend test) throughout the study period indicate that this signal does not relate to physiological acclimation. This is in line with field observations7,14,15 and laboratory studies16–18, which found cells to retain their photosynthetic machinery active over the polar night so that it can be instantly used upon the return of light. Even though our measurements represent only the potential capacity of the existing biomass to fix carbon rather than the actual in situ rates, the increase in potential NPP that we identified from March 14th onwards is highly indicative of an actual increase in photosynthetic biomass despite the extremely low light levels this early in the season.Fig. 2 Temporal development in phytoplankton activity and biomass.

Development of (A) potential net 14C productivity (as a proxy for Net Primary Production) under reference conditions, and (B) 11 m water column Chlorophyll a (Chl-a) concentrations (light green) and 50 m depth-integrated Chl-a standing stocks (dark green) as a function of time. Chl-a-specific 14C productivity production is displayed in Fig. S4.

As a more direct measure of in situ biomass production through photosynthesis, we used surface ocean Chl-a concentrations from discrete underway samples. These were obtained from 11 m water depth on a daily basis, which allowed us to derive a highly resolved time series of Chl-a-based biomass in the upper mixed layer of the ocean (Fig. 2B). This conversion of Chl-a concentration into biomass is possible because Chl-a specific carbon fixation and the ratio of POC to Chl-a did not display a trend over the study period (Figs. S4, S5, Mann–Kendall Trend test).

Chl-a concentrations measured during the first half of March were similar to those measured during February (i.e., around 0.01–0.02 µg L−1), but then increased exponentially from mid-March to the end of April, reaching 0.05 µg L−1 by mid-April. We applied a change-point analysis package19 to this data to estimate the onset of photosynthetic growth. In line with a visual inspection of the time series (Fig. 2B), six independent models from this package consistently identified March 27th ± 1 day as the change point in the Chl-a time series. Opting for the more conservative estimate, we identified March 28th as the onset of photosynthetic growth in this dataset. The rate of increase after this date fits well to the rate of change in surface Chl-a concentrations calculated over 5-day intervals in our time series20: While the rate of change varied between slightly negative and slightly positive values during most of March, rates were positive for nine consecutive days from March 25th and remained mainly positive at higher rates for the rest of the study period (Fig. S6). Congruently, accumulation rate constants based on underway Chl-a over 7-day intervals were −0.01 ± 0.04 d−1 before March 27th and +0.06 ± 0.07 d−1 after the change point. Depth-integrated Chl-a from CTD-rosette samples (upper 50 m, Fig. 2B), being available only in lower temporal resolution, followed a similar trend with standing stocks increasing between March 27th and April 5th for the first time, with standing stocks of 0.6 and 0.8 mg m−2, respectively. In summary, as a conservative estimate, March 28th marked the onset of photosynthetic growth in the water column underneath the ice cover as determined from the measurement of Chl-a concentration.

As a third method to determine the increase in primary production and the subsequent onset of photosynthetic growth, we used cell counts via light microscopy and flow cytometry (FCM) of water samples from 20 m water depth (Fig. S7). Given the need to buildup biomass before cell division, we expected to see cell counts to increase sometime after the actual onset of net photosynthesis. Consistent with this expectation, cell counts indeed increased a few weeks after the identified onset of net photosynthesis on March 28th, mainly driven by the pennate diatom Pseudo-nitzschia spp. (Table S1).

Additionally, we also examined the temporal development of ice algal communities21 during the study period. These were sampled in lower temporal resolution owing to the inaccessibility of the ice coring areas under the prevailing harsh weather and ice conditions12. Chl-a concentrations increased ten-fold in the bottom 10 cm of first-year sea ice from mid-February to the end of March (Fig. S8). Similarly, FCM cell counts indicated a 5 to 10-fold increase in photosynthetic cell numbers between these sampling dates, mainly driven by an increase in the nanoplanktonic size class (2–20 µm; Fig. S9). Light microscopic cell counts confirmed this pattern, indicating that different pennate diatoms (mainly Nitzschia spp., and Navicula spp.) dominated this early increase in the ice algal assemblage (Table S1). Using the exponential rate of increase from the end of March to the end of April (r2 = 0.86), we estimate that the first doubling of Chl-a concentrations in first-year bottom ice occurred as early as March 11th. Also, potential ice-algae NPP (Fig. S8) increased by an order of magnitude between February and late March.

Having established that photosynthetic growth started between mid-March and early April, depending on the sensitivity of the detection method used, we will now discuss the light levels under which this early photosynthesis took place. For this analysis, we draw on a number of optical sensors that were deployed during the MOSAiC expedition: above, within, and underneath the ice (Figs. 1, S2). These sensors provided continuous PAR measurements during the transition from polar night to polar day (Fig. 3).Fig. 3 Temporal development in light availability.

Photosynthetically active radiation (PAR) at the sea ice bottom (blue) and in the water column (3–50 m depth-integrated values; green) for the different OptiCALs (Optical Chain And Logger sensors ee, gg, hh) and light harp as a function of time. Shaded areas indicate error bands of 20% and 12.2% uncertainty for OptiCALs and lightharp, respectively. A lead that opened on April 17th near sensor gg strongly increased light availability in the water column at this location for 3–4 days.

To obtain the light intensity in the water column at which NPP set in, we used measurements from three OptiCAL (Optical Chain And Logger) autonomous ice-tethered observatories that sampled downwelling irradiance in the spectrally integrated photosynthetic active radiation (PAR) wavelength range 400–700 nm [µmol photons m−2 s−1]22. These light chains with nominal sensor depths of 1, 2, 3, 4, 5, 6, 7, 9, 14, 21, 32, and 50 m below the ice surface were deployed at three different locations, allowing us to estimate the spatial variability of the under-ice light field. To determine the mean irradiance that primary producers in the water column were exposed to, we averaged the light fields horizontally across these three sites that continuously drifted with the ice over the underlying ocean23,24. Vertically, we averaged our measurements across the upper 50 m of the water column. This vertical averaging is motivated by the fact that the Chl-a profiles varied little over this depth range, indicating that cells were likely mixed homogeneously in the upper 50 m and experienced the average light level across this depth range25,26. In addition, this depth range agrees well with an independently estimated mixed layer depth of 40–65 m23,27 and a Brunt-Väisälä frequency maximum between 50 and 60 m23,27. We consider the mean light field 24 h prior to the biological sampling, as NPP is an integrated signal commonly quantified on a daily basis. To account for potential horizontal inhomogeneities in the translucency of the ice that could have caused an area-averaged PAR above the specific measurement levels at the sensor sites, we have added generous uncertainty margins to all light measurements, thus ensuring that the values reported here provide an absolute upper bound on the true light intensity that the cells were exposed to local shading effects were accounted for by using a correction model for the euphotic zone under sea ice28, while the impact of inhomogeneities in the snow cover was accounted for by adding 33% to all light-intensity measurements, a correction term that has been derived from the statistical distribution of snow thicknesses in the measurement area (see “Methods” for details).

Using this approach, we find that PAR values underneath first-year ice (thickness: 1.14–2.05 m21) were similar to those under second-year ice (thickness 1.71–2.44 m29). This is related to the dominating impact of the snow cover on the under-ice light field30. The snow cover at the MOSAiC site varied little between ice areas of different ages and barely changed between February and the onset of melt in May31. Spatial differences in snow depth and thus light availability were instead primarily related to the dynamic deformation of the ice cover, contradicting the common notion that first-year ice generally supports higher primary production due to the assumed larger light availability32,33. For the purpose of this study, we therefore do not differentiate between first and second-year ice.

The light measurements indicate that our conservative estimate of the first day of positive water column Chl-a accumulation on March 28th was associated with a daily average PAR of not more than 0.04 ± 0.02 µmol photons m−2 s−1. After this, Chl-a concentration increased exponentially with increasing PAR (Fig. 4). Given that alternative energy sources such as respiration of storage compounds or mixotrophic uptake of organic matter have been available throughout winter and are not related to the return of sun light18,34,35, the exponential increase in Chl-a concentrations must be a direct response to the measured increase in PAR.Fig. 4 Relationship between photosynthetic biomass and light availability.

Underway Chlorophyll a (Chl-a) concentrations (11 m depth; log scale) as a function of 3–50 m depth-integrated daily average photosynthetically active radiation (PAR) values with exponential fits of the period before (r2 = 0.03 indicating no relationship between Chl-a and PAR) and after (r2 = 0.66 indicating an exponential increase in Chl-a with increasing PAR) the change point determined after Beaulieu & Killick18. A lead that opened on April 17th at sensor gg strongly increased light availability at this location for 3–4 days (see Fig. 3). The respective data points were excluded from this calculation (see Fig. S9 for full data).

The conservatively estimated PAR value of 0.04 µmol photons m−2 s−1 is two to three orders of magnitude lower than the 0.415 mol photons m−2 d−1 isolume commonly used to constrain the zone of NPP in biogeochemical models36–38, one order of magnitude lower than the lowest previous record for ice algae4, and two orders of magnitude lower than previous estimates from surface waters7. This surprisingly low PAR value is very close to the theoretical minimum photon requirement for oxygenic photosynthesis of 0.01 µmol photons m−2 s−1 4–6, and provides evidence that evolution has optimized the efficiency of photosynthesis in microalgae to an astonishing degree.

Our finding that NPP can occur at extremely low light levels has substantial implications for our understanding and quantification of processes driving ocean ecosystems, not only in the Arctic but potentially in all light-limited marine habitats on Earth. For the Arctic, knowing the lower light limit for microalgae is crucial for understanding the timing of primary production and matching with their zooplankton grazers and higher trophic levels39. Even though the measured early spring Chl-a biomass levels in the MOSAiC study region were very low, they are comparable to typical Arctic summer surface Chl-a concentrations (>80°N, 0.33 ± 0.55 µg L−1; n = 148740–42), and represent a substantial reservoir of organic material that higher trophic levels can rely on. In addition, the early production of biomass that we found in bottom ice led to Chl-a concentrations of more than 5 µg L−1 in early April, substantially more than the average Chl-a concentration of 1.0 ± 1.3 µg L−1 previously observed during summer in sea ice of the Central Arctic Ocean (>80°N, n = 11240,43,44). This implies that the high Arctic Ocean ecosystem may be similarly fueled by such low light productivity in the severely undersampled spring and fall, with the short summer period being less important than previously thought. Our early spring measurements of phototrophic biomass underneath the ice and snow cover illustrate that satellite-derived Chl-a and NPP in ice-covered waters45 might be underestimating early season biomass buildup. Our results also challenge the common notion that summer microalgal growth in the perennially ice-covered Arctic Ocean is light-limited45. Instead, nutrient limitation46 and grazing pressure47 might be the primary controlling factors of biomass buildup during the polar day.

Our data from the Arctic Ocean indicate that net photosynthetic biomass buildup close to the theoretical minimum light requirement is physiologically possible, despite certain losses for example due to grazing. While the specific conditions of the Arctic polar night and the return of sunlight may represent a unique opportunity to quantify this light intensity threshold in situ, the fundamental capacity of photosynthesis to be this efficient likely also applies to other low light environments. Our results thus imply that phytoplankton communities in the world’s oceans during low light seasons and at greater depths throughout the year may be significantly more productive than previously thought. As a consequence, the depth of the euphotic zone in large parts of the world’s oceans may be considerably deeper than previously reported, making our measurements highly relevant for understanding life in the ocean’s twilight zones. To quantify this possible impact, our estimate corresponds to a much larger euphotic zone depth with a threshold of 0.003% of an incoming irradiance of 2000 µmol photons m−2 s−1. In fact, this would translate into euphotic zone depths more than twice as deep compared to the 1% criterion. Using the previous examples and assuming, for example, an attenuation coefficient of 0.2 m−1, our threshold would deepen the bottom of the euphotic zone from 23 to 54 m. This substantially increases the vertical extent and thus the total volume of the euphotic zone in the world’s oceans3 and may change our view on upper twilight zone ecology and biogeochemistry48. As such, our finding that net photosynthesis is possible at extremely low light levels can have consequences for our understanding of the entire world’s oceans and points towards large potentials for hidden and previously overseen marine primary production.

Methods

The measurements reported in this study were conducted in the framework of the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) drift expedition in the Central Arctic Ocean. During the study period February to April 2020, RV Polarstern49 drifted with the sea ice from 88°N to 84°N12,45 (see Figs. S1 and S2).

Water column sampling

Samples from the upper mixed layer were collected approximately once per week at 20 m depth from Niskin bottles either from the ship’s CTD rosette, or a smaller rosette deployed on the ice floe13,24. Additional upper mixed layer samples with higher temporal resolution were taken daily from the ship’s underway system with an intake at 11 m water depth.

Sea ice sampling

Sea ice samples were collected from level first-year ice (FYI) at a site located ~1.5 km away from the ship, where light pollution from the ship was too low to be measurable even during the polar night12,13. Cores were collected using a Kovacs Mark II 9 cm inner diameter corer and were sectioned from both the ice surface and bottom directly in the field. To allow for co-located measurements of key parameters from the same samples, sections from 3–4 ice cores were pooled in Whirl-Pack® bags. Back on the ship, ice samples (n = 1 per sampling event) were further processed in a 1 °C cold container under dim light conditions (NPP samples were only handled under red light). For salinity-buffered melting50, 500 mL of 0.2 µm filtered surface seawater (11 m) was added per 10 cm ice. NPP samples consisted of less consolidated bottom sections only, which were crushed and melted in the dark at room temperature for 1–2 h. Pooled subsamples for all other parameters were melted in the dark at room temperature overnight. The total volume of the melted pools was measured before samples were split for different sample pipelines. Total volumes, filtered seawater, and subsample volumes were used to calculate corrected sample volumes.

Chlorophyll a analysis

Samples for the analysis of Chlorophyll a (Chl-a) were filtered onto pre-combusted glass fiber filters (GF/F; 0.7 µm nominal pore size) using mild vacuum. Samples (n = 1 per sampling event) were stored in the dark at −80 °C. Water column samples of 4 L were collected at 2, 10, 20, and 50 m depths relative to sea level. Bottom sea ice samples were collected in two 5 cm thick sections. Samples were extracted in 90% acetone overnight at 4 °C and subsequently analyzed on a fluorometer (TD-700; Turner Designs, USA), including an acidification step (1 M HCl) to determine phaeopigments51. One subsample per underway sampling event was measured on board within 3 days after sampling, while a second subsample per underway sampling as well as the weekly samples were analyzed at the Alfred Wegener Institute (AWI) after the campaign had ended (5–9 months after sampling). No systematic difference between the replicates could be detected, indicating that no significant degradation of Chl-a took place during storage and transport. For sampling events with more than one replicate (i.e., the daily underway sampling), daily averages were used for further analyses (change point, rate of change).

Particulate organic carbon measurements

Samples for measurements of POC concentrations were collected by filtering water on pre-combusted (450 °C, 4 h) glass fiber filters (GF/F; 0.7 µm nominal pore size) using mild vacuum. Filters (n = 1 per sampling event) were freeze-dried and carbonates were removed by contact with HCl fumes for 4 h in a vacuum-enclosed system. Filters were packed into tin cups and analyzed with an elemental analyzer (Flash 2000, Thermo Scientific, Milan, Italy) at the University of La Rochelle, France (Littoral, Environment and Societies Joint Research Unit stable isotope facility).

Light microscopic analysis

Samples for light microscopic analysis (n = 1 per sampling event) were fixed with hexamethylenetetramine-buffered formalin solution (1% final concentration). Samples were stored at 4 °C in the dark until analysis. Samples from the water column were analyzed via inverse light microscopy at AWI (Germany), while samples from ice cores were analyzed at IOPAN (Poland).

Flow cytometry

Samples for flow cytometric abundances, collected in triplicate, were preserved by adding 38 µL of glutaraldehyde to 1.5 mL of sample. Samples were kept for ~2 h at 4 °C in the dark, flash-frozen in liquid nitrogen, and then stored at −80 °C until analysis. The abundance of pico- and nano-sized phytoplankton (approximate cell size 0.2–2 μm and 2–20 μm, respectively) were determined using an Attune® NxT, Acoustic Focusing Cytometer (Invitrogen by Thermo Fisher Scientific) with a syringe-based fluidic system and a 20 mW 488 nm (blue) laser. Cells were counted after thawing the sample (n = 1 per sampling event) and the various groups discriminated based on their red fluorescence (BL3) vs. side scatter (SSC), red fluorescence (BL3) vs. orange fluorescence (BL2) and orange fluorescence (BL2) vs. side scatter (SSC). The gating strategy is illustrated in Fig. S12.

Primary production measurements

Potential NPP was determined by duplicate incubation (n = 2 per sampling event) with NaH14CO3 spike (53.1 mCi mmol−1; Perkin Elmer; applied specific activity of 0.5 µCi mL−1) for 24 h under reference conditions (1.0 ± 0.5 °C and 10 ± 3 µmol photons m−2 s−1) together with a dark control. For water column measurements, incubations contained 500 mL per sample, while sea ice incubation had volumes of 250–500 mL. After 24 h, incubated samples were filtered onto GF/F filters, acidified with 200 µL of 1 M HCl, and left to degas overnight. After addition of 10 mL of scintillation cocktail (Ultima Gold AB, PerkinElmer), samples were vortexed and left to stand in the dark for 20–30 h before counting on the liquid scintillation counter, using automatic quench correction and a counting time of 5 min. For blank determination, one replicate was immediately filtered and acidified. Per incubation bottle, 100 µL aliquots were mixed with 100 µL ethanolamine immediately after spiking to determine the total amount of added NaH14CO3. Subtracted blank values were on average 25% (n = 23) of the incubated sample counts. Potential NPP was calculated as1 NPP=(DICDPMsample-DPM0%1.05)/(DPM100%t)

where [DIC] denotes the concentrations of dissolved inorganic carbon in the sample. DPMsample denotes the disintegrations per minute (DPM) in the samples, DPM0% reflects the blank value, DPM100% denotes the DPM of the total amount of NaH14CO3 added to the samples, and t is the duration of the incubation.

For the calculations, DIC data from the MOSAiC expedition52 was used. The value of 1.05 is used to correct for fractionation against 14C relative to 12C53. C fixation in dark controls was subtracted to derive NPP.

Data analysis

To detect the timing of photosynthetic biomass buildup initiation, we calculated the rate of change for daily measurements of underway Chl-a, using a 5-day interval54, and detected the initiation as the first period in which dx/dt was positive for 5 days or more. Changepoint analysis19 for daily underway Chl-a measurements were performed via the EnvCrp package in R55. Out of the twelve models in the package, six returned successful fits (the others failed as they would have required a larger dataset size). The six functional models returned three consecutive values (with corresponding dates) as the change points of the dataset. More details are provided in the supplementary information (Supplementary Note 1). The average of these six models was used with one standard deviation.

Light measurements in the water column

Light measurements in the mixed water layer underneath the ice were obtained from OptiCALs22, which are autonomous instruments that provide vertically resolved downwelling irradiance as a proxy for downward planar irradiance. Three OptiCALs were installed at the MOSAiC site, referred to as “hh”, “gg”, and “ee”. OptiCAL hh was deployed at the SYI dark coring site (LM-site), OptiCAL ee at the so-called Ridge site, and OptiCAL gg at the so-called Monster site (Fig. S2). All three instruments were calibrated to absolute PAR intensity (photosynthetic photon flux density between 400 nm to 700 nm) with a remaining uncertainty of ±20%56.

Because of inhomogeneities of the ice and snow, the raw measurements of the three OptiCALs have to be processed to obtain a robust average light field underneath the ice. In doing so, we took care to rather err towards an estimate of too much light, such that our resulting light intensity provides an upper bound of the true light intensity.

Our processing comprised two steps, with the first accounting for local inhomogeneities at the measuring sites, and the second for inhomogeneities across the ice cover. In the first step, we correct the individual light profiles at the measuring sites. These profiles generally follow an exponential decay with depth, as would be expected for a homogenous light field at the top of the water column and constant light-transmission properties of the water column57. However, since both these conditions are not perfectly fulfilled, slight deviations from an exponential decay curve were apparent in the profiles. For example, the average light intensity across the narrow surface area sampled by the uppermost sensor can differ from the average light intensity across the much wider surface area sampled by the sensors further down. To correct for these local surface inhomogeneities we used the correction model of Laney et al28. The model eliminates near-surface shaded areas in the irradiance profiles. Occasionally, the model was unable to reliably approximate the measured light profile, resulting in an excessively high standard deviation of the fit (i.e., σ > 1 µmol photons m−2 s−1) or an unreasonably high diffuse attenuation coefficient for such low-biomass waters (i.e., κ > 0.06 m−1 58); the affected profiles were left uncorrected (see Fig. S13 for example of a valid and an invalid fit). For the two lowest sensors at 32 m and 50 m depth, we found that a linear decay curve better represented the measured profiles, and therefore used a linear fit for this lowermost part of the profiles.

Having thus obtained profiles that consider the local inhomogeneities at the measurement site, in a second step we corrected for horizontal inhomogeneities in the light field underneath the entire ice cover. These horizontal inhomogeneities are almost exclusively related to differences in the snow thickness, as described in the main paper. The spatio-temporal variability of the snow cover across the MOSAiC site was obtained from snow-thickness transect data59. These transects showed that the snow cover at the MOSAiC site varied little between sea ice of different ages and remained largely unchanged between March and the melt onset in May. Spatial differences in snow thickness were mainly related to the three observed sea-ice surface types level ice, rubble ice, and deformed ice31 of which the ice cover consisted to roughly equal shares. All optical sensors were deployed on level ice, which generally had the thinnest snow cover and the thinnest ice thickness, and therefore the highest light transmissivity. Snow thickness at the optical measurement sites was 0.19–0.21 m throughout March and April as measured by snow buoy 2019S9660, which was deployed at the LM-site along with OptiCAL hh. Snow-thickness measurements at OptiCALs gg and ee gave very similar values. These snow thicknesses of around 20 cm were consistently lower than the mean snow thickness along transect lines S-loop and N-loop59 which varied between 0.25 m and 0.30 m. This in isolation suggests that our light measurements are at the upper end of the true light intensity underneath the ice, as it was on average covered by thicker snow than at the measuring sites. However, this is countered by very high light transmission in regions with particularly thin snow thickness as found along the snow-thickness transect. To account for these inhomogeneities, we determined the ratio between the PAR-transmissivity (Tsnow,PAR) averaged across the snow transects and the PAR-transmissivity averaged across the deployment sites. To obtain Tsnow, PAR, we use a standard radiative transfer model for ice and snow61 given as2 Tsnow,PAR=I0⋅e−κsnow(hsnow−0.03m)

This model accounts for a surface scattering layer with an initial transmissivity of I0 of the uppermost 3 cm of the snow cover and exponential attenuation below. The surface transmission value I0 is irrelevant for our purposes, as its value cancels when calculating the ratio between different values of Tsnow, PAR. The attenuation coefficient of snow for downwelling PAR irradiance κsnow was obtained from exponential fits to MOSAiC irradiance measurements above the snow62 and below the snow63. These measurements result in a mean attenuation coefficient of κsnow=14.1m−1±1.4m−1 for the snow layer at the MOSAiC site in March/April, including all diurnal solar zenith angles and sky conditions. The average value of Tsnow, PAR across all snow thickness measurements along the transect lines S-loop and N-loop from March and April59, was 33 ± 13% higher than the average value of Tsnow, PAR at the deployment sites of the OptiCALs, primarily because of the high transmissivity in areas of exceptionally thin snow. To account for the related potential underestimation of area-averaged transmitted light, we therefore corrected the raw measurements at the deployment sites by adding 33% to their values. Note that recent more advanced formulations of radiative transfer parametrizations64 yield smaller biases of the order of 10% based on the MOSAiC measurements such that the light intensity after a correction of +33% provides an upper limit of the under-ice PAR.

The mean irradiance in the surface mixed layer was then finally calculated by spatially (3–50 m water depth) and temporally (24 h before ECO sampling) averaging the corrected light profiles. On 28 March, as conservatively identified as the day of onset of photosynthesis, the three measurement sites gave average light intensities across the mixed layer of 0.014, 0.033, and 0.061 µmol photons m−2 s−1, giving an upper bound for the light intensity required for photosynthesis of 0.04 ±0.02 µmol photons m−2 s−1.

Light measurements in the ice

Light measurements in the ice were obtained from a lightharp65, an instrument that measures profiles of upward and downward irradiance in the ice. Calibration of the light harp to absolute PAR intensity was possible with a remaining upper bound for potential biases of 12.2%65. For this study, the calibrated light harp data were used as independent data to qualitatively examine the long-term stability and consistency of the OptiCAL data, to determine PAR directly at the underside of the ice, and to obtain the diffuse attenuation coefficient of snow from measurements directly at the snow—ice interface.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51636-8.

Acknowledgements

Data used in this manuscript was produced as part of the international Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) with the tag MOSAiC20192020 and the Project ID AWI_PS122_00. C.J.M.H., N.F., and D.N. received funding from the German Ministry of Education and Research (BMBF) through the project NiceLABpro (grant no. 03F0867A). D.N. received funding from the Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy (EXC 2037; CLICCS—Climate, Climatic Change, and Society; grant no. 390683824). A.L., E.L., G.B., O.M., and P.A. and were supported by the Research Council of Norway through project HAVOC (grant no. 280292). BL was supported by a Fellowship at the Hanse-Wissenschaftskolleg Institute for Advanced Study, Delmenhorst, Germany. L.O. received funding from the BMBF-funded project nuArctic (grants 03F0918A). We thank the cruise participants, ship’s crew, and logistics support for MOSAiC legs 2 and 3 as well as everyone else who contributed to the realization of MOSAiC (Nixdorf et al. 2021). Particularly, we would like to thank those colleagues contributing to FYI coring and water sampling, i.e., A. Ulfsbo, S. Fons, S. Spahic, C. Marsay, L. Eggers, J. Grosse, G. Castellani, D. V. Divine, L. M. Olsen, S. Torres-Valdez, E. Damm, A. Dumitrascu, and P. Simões Pereira. We thank the MOSAiC leg 2 and 3 CTD teams, particularly J. Schaffer, for their help with water sampling. Sea ice light microscopy data was provided by M. Różańska-Pluta, A. Tatarek, and J. M. Wiktor. We thank L. Heitmann, T. Brenneis, and A. Terbrüggen for help with Chl-a measurements. D.V. Divine, C. Katlein, P. Itkin, and I. Raphael are acknowledged for deploying the OptiCALs and the lightharp. We thank G. Birnbaum for the coordination and recording of aerial images. Y. Novak is acknowledged for producing the schematic figure of the sampling sites. Invaluable support, planning, and coordination were provided by A.A. Fong and R. Gradinger.

Author contributions

C.J.M.H. developed the research question and concept. C.J.M.H., D.N., and N.F. developed the overall storyline. C.J.M.H. and A.T. conducted field measurements. C.J.M.H., A.K., and O.M. measured biological samples. CJMH analyzed biological data. N.F. analyzed the optical measurements. CJMH and DN wrote the paper, with contributions from N.F., P.An., P.As., J.B., G.B., G.G., A.K., A.L., B.L., E.L., M.L., O.M., L.O., B.R., B.S., A.T., and J.W.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Funding

Open Access funding enabled and organized by Projekt DEAL.

Data availability

The data used in this study have been deposited in the PANGAEA repository (https://www.pangaea.de/). Cruise track data used in Fig. S1 is available as Kanzow (202066, 10.1594/PANGAEA.924681). Continuous light data from OptiCALs are available as Anderson et al., (2023a67, 10.1594/PANGAEA.928495), Anderson et al., (2023b57, 10.1594/PANGAEA.955045), and Anderson et al., (2023c68, 10.1594/PANGAEA.954849). Continuous light data from the lightharp are available as Fuchs et al., (2024a,b63,65, 10.1594/PANGAEA.963743). Chl-a concentration data from the underway system are available as Hoppe et al., (2023a69, 10.1594/PANGAEA.963277) and the CTD rosette casts are available as Hoppe et al., (2023b70, 10.1594/PANGAEA.962597). Flow cytometric cell count data from the water column are available as Müller et al., (2023a71, 10.1594/PANGAEA.963430), and data from sea ice are available as Müller et al., (2023b72, 10.1594/PANGAEA.963560). Data on microscopic cell counts are available as Kraberg (202473, 10.1594/PANGAEA.965913); the summarized data as used in this study are provided in the Supplementary Information (Table S1). Source data are provided with this paper.

Code availability

The Python code to process PAR data in the described way is available under Fuchs (202474, 10.5281/zenodo.12772364). The R package EnvCPT (v1.1.3) used for change point analyses is available on GitHub (https://github.com/rkillick/EnvCpt/).

Competing interests

The Authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Falkowski PG The role of phytoplankton photosynthesis in global biogeochemical cycles Photosynth. Res. 1994 39 235 258 10.1007/BF00014586 24311124
Falkowski, P. G. The role of phytoplankton photosynthesis in global biogeochemical cycles. Photosynth. Res. 39, 235–258 (1994).24311124 10.1007/BF00014586
2. Field CB Behrenfeld MJ Randerson JT Falkowski P Primary production of the biosphere: integrating terrestrial and oceanic components Science 1998 281 237 240 10.1126/science.281.5374.237 9657713
Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237–240 (1998).9657713 10.1126/science.281.5374.237
3. Marra JF Lance VP Vaillancourt RD Hargreaves BR Resolving the ocean’s euphotic zone Deep Sea Res. Part Oceanogr. Res. Pap. 2014 83 45 50 10.1016/j.dsr.2013.09.005
Marra, J. F., Lance, V. P., Vaillancourt, R. D. & Hargreaves, B. R. Resolving the ocean’s euphotic zone. Deep Sea Res. Part Oceanogr. Res. Pap. 83, 45–50 (2014).10.1016/j.dsr.2013.09.005
4. Hancke K Extreme low light requirement for algae growth underneath sea ice: a case study from station Nord, NE Greenland J. Geophys. Res. Oceans 2018 123 985 1000 10.1002/2017JC013263
Hancke, K. et al. Extreme low light requirement for algae growth underneath sea ice: a case study from station Nord, NE Greenland. J. Geophys. Res. Oceans 123, 985–1000 (2018).10.1002/2017JC013263
5. Randelhoff A Arctic mid-winter phytoplankton growth revealed by autonomous profilers Sci. Adv. 2020 6 eabc2678 10.1126/sciadv.abc2678 32978152
Randelhoff, A. et al. Arctic mid-winter phytoplankton growth revealed by autonomous profilers. Sci. Adv. 6, eabc2678 (2020).32978152 10.1126/sciadv.abc2678
6. Raven JA Kübler JE Beardall J Put out the light, and then put out the light J. Mar. Biol. Assoc. UK. 2000 80 1 25 10.1017/S0025315499001526
Raven, J. A., Kübler, J. E. & Beardall, J. Put out the light, and then put out the light. J. Mar. Biol. Assoc. UK. 80, 1–25 (2000).10.1017/S0025315499001526
7. Hoppe CJM Always ready? Primary production of Arctic phytoplankton at the end of the polar night Limnol. Oceanogr. Lett. 2022 7 167 174 10.1002/lol2.10222
Hoppe, C. J. M. Always ready? Primary production of Arctic phytoplankton at the end of the polar night. Limnol. Oceanogr. Lett. 7, 167–174 (2022).10.1002/lol2.10222
8. Behrenfeld MJ Revaluating ocean warming impacts on global phytoplankton Nat. Clim. Change 2016 6 323 330 10.1038/nclimate2838
Behrenfeld, M. J. et al. Revaluating ocean warming impacts on global phytoplankton. Nat. Clim. Change 6, 323–330 (2016).10.1038/nclimate2838
9. Smetacek V Passow U Spring bloom initiation and Sverdrup’s critical‐depth model Limnol. Oceanogr. 1990 35 228 234 10.4319/lo.1990.35.1.0228
Smetacek, V. & Passow, U. Spring bloom initiation and Sverdrup’s critical‐depth model. Limnol. Oceanogr. 35, 228–234 (1990).10.4319/lo.1990.35.1.0228
10. Robinson DH Arrigo KR Iturriaga R Sullivan CW Microalgal light-harvesting in extreme low-light environments in McMurdo Sound, Antarctica J. Phycol. 1995 31 508 520 10.1111/j.1529-8817.1995.tb02544.x
Robinson, D. H., Arrigo, K. R., Iturriaga, R. & Sullivan, C. W. Microalgal light-harvesting in extreme low-light environments in McMurdo Sound, Antarctica. J. Phycol. 31, 508–520 (1995).10.1111/j.1529-8817.1995.tb02544.x
11. Mock, R. & Gradinger, T. Changes in photosynthetic carbon allocation in algal assemblages of Arctic sea ice with decreasing nutrient concentrations and irradiance. Mar. Ecol. Prog. Ser. 202, 1–11 (2000).
12. Nicolaus M Overview of the MOSAiC expedition: Snow and sea ice Elem. Sci. Anthr. 2022 10 000046 10.1525/elementa.2021.000046
Nicolaus, M. et al. Overview of the MOSAiC expedition: Snow and sea ice. Elem. Sci. Anthr. 10, 000046 (2022).10.1525/elementa.2021.000046
13. Fong, A. A. et al. Overview of the MOSAiC expedition: Ecosystem. Elementa: Science of the Anthropocene, 12, 10.1525/elementa.2023.00135 (2024).
14. Berge J Unexpected levels of biological activity during the polar night offer new perspectives on a warming Arctic Curr. Biol. 2015 25 2555 2561 10.1016/j.cub.2015.08.024 26412132
Berge, J. et al. Unexpected levels of biological activity during the polar night offer new perspectives on a warming Arctic. Curr. Biol. 25, 2555–2561 (2015).26412132 10.1016/j.cub.2015.08.024
15. Kvernvik AC Fast reactivation of photosynthesis in Arctic phytoplankton during the polar night J. Phycol. 2018 54 461 470 10.1111/jpy.12750 29723414
Kvernvik, A. C. et al. Fast reactivation of photosynthesis in Arctic phytoplankton during the polar night. J. Phycol. 54, 461–470 (2018).29723414 10.1111/jpy.12750
16. Morin, P.-I. et al. Response of the sea-ice diatom Fragilariopsis cylindrus to simulated polar night darkness and return to light. Limnol. Oceanogr. 65, 1041–1060 (2019).
17. Lacour T Decoupling light harvesting, electron transport and carbon fixation during prolonged darkness supports rapid recovery upon re-illumination in the Arctic diatom Chaetoceros neogracilis Polar Biol. 2019 42 1787 1799 10.1007/s00300-019-02507-2
Lacour, T. et al. Decoupling light harvesting, electron transport and carbon fixation during prolonged darkness supports rapid recovery upon re-illumination in the Arctic diatom Chaetoceros neogracilis. Polar Biol. 42, 1787–1799 (2019).10.1007/s00300-019-02507-2
18. Kennedy F Martin A Bowman JP Wilson R McMinn A Dark metabolism: a molecular insight into how the Antarctic sea-ice diatom Fragilariopsis cylindrus survives long-term darkness N. Phytol. 2019 223 675 691 10.1111/nph.15843
Kennedy, F., Martin, A., Bowman, J. P., Wilson, R. & McMinn, A. Dark metabolism: a molecular insight into how the Antarctic sea-ice diatom Fragilariopsis cylindrus survives long-term darkness. N. Phytol. 223, 675–691 (2019).10.1111/nph.15843
19. Beaulieu C Killick R Distinguishing trends and shifts from memory in climate data J. Clim. 2018 31 9519 9543 10.1175/JCLI-D-17-0863.1
Beaulieu, C. & Killick, R. Distinguishing trends and shifts from memory in climate data. J. Clim. 31, 9519–9543 (2018).10.1175/JCLI-D-17-0863.1
20. Behrenfeld MJ Boss ES Student’s tutorial on bloom hypotheses in the context of phytoplankton annual cycles Glob. Change Biol. 2018 24 55 77 10.1111/gcb.13858
Behrenfeld, M. J. & Boss, E. S. Student’s tutorial on bloom hypotheses in the context of phytoplankton annual cycles. Glob. Change Biol. 24, 55–77 (2018).10.1111/gcb.13858
21. Oggier, M. et al. First-year sea-ice salinity, temperature, density, oxygen and hydrogen isotope composition from the main coring site (MCS-FYI) during MOSAiC legs 1 to 4 in 2019/2020. PANGAEA [data] 10.1594/PANGAEA.956732 (2023).
22. Schartmüller B Development and calibration of a high dynamic range and autonomous ocean-light instrument to measure sub-surface profiles in ice-covered waters Appl. Opt. 2023 62 8308 8315 10.1364/AO.502437 38037934
Schartmüller, B. et al. Development and calibration of a high dynamic range and autonomous ocean-light instrument to measure sub-surface profiles in ice-covered waters. Appl. Opt. 62, 8308–8315 (2023).38037934 10.1364/AO.502437
23. Schulz K The Eurasian Arctic Ocean along the MOSAiC drift in 2019–2020: an interdisciplinary perspective on physical properties and processes Elem. Sci. Anthr. 2024 12 00114 10.1525/elementa.2023.00114
Schulz, K. et al. The Eurasian Arctic Ocean along the MOSAiC drift in 2019–2020: an interdisciplinary perspective on physical properties and processes. Elem. Sci. Anthr. 12, 00114 (2024).10.1525/elementa.2023.00114
24. Rabe B Overview of the MOSAiC expedition: physical oceanography Elem. Sci. Anthr. 2022 10 00062 10.1525/elementa.2021.00062
Rabe, B. et al. Overview of the MOSAiC expedition: physical oceanography. Elem. Sci. Anthr. 10, 00062 (2022).10.1525/elementa.2021.00062
25. MacIntyre HL Kana TM Geider RJ The effect of water motion on short-term rates of photosynthesis by marine phytoplankton Trends Plant Sci. 2000 5 12 17 10.1016/S1360-1385(99)01504-6 10637656
MacIntyre, H. L., Kana, T. M. & Geider, R. J. The effect of water motion on short-term rates of photosynthesis by marine phytoplankton. Trends Plant Sci. 5, 12–17 (2000).10637656 10.1016/S1360-1385(99)01504-6
26. Kvernvik AC Arctic sea ice algae differ markedly from phytoplankton in their ecophysiological characteristics Mar. Ecol. Prog. Ser. 2021 666 31 55 10.3354/meps13675
Kvernvik, A. C. et al. Arctic sea ice algae differ markedly from phytoplankton in their ecophysiological characteristics. Mar. Ecol. Prog. Ser. 666, 31–55 (2021).10.3354/meps13675
27. Schulz, K., Koenig, Z. & Muilwijk, M. The Eurasian Arctic Ocean along the MOSAiC drift (2019–2020): Core hydrographic parameters. 10.18739/A21J9790B (2023).
28. Laney SR Krishfield RA Toole JM The euphotic zone under Arctic Ocean sea ice: vertical extents and seasonal trends Limnol. Oceanogr. 2017 62 1910 1934 10.1002/lno.10543
Laney, S. R., Krishfield, R. A., & Toole, J. M. The euphotic zone under Arctic Ocean sea ice: vertical extents and seasonal trends. Limnol. Oceanogr. 62, 1910–1934 (2017).10.1002/lno.10543
29. Oggier, M. et al. Second-year sea-ice salinity, temperature, density, oxygen and hydrogen isotope composition from the main coring site (MCS-SYI) during MOSAiC legs 1 to 4 in 2019/2020. PANGAEA [data] 10.1594/PANGAEA.959830 (2023).
30. Petrich C Eicken H Growth, structure and properties of sea ice Sea Ice 2010 2 23 77
Petrich, C. & Eicken, H. Growth, structure and properties of sea ice. Sea Ice 2, 23–77 (2010).
31. Itkin P Sea ice and snow characteristics from year-long transects at the MOSAiC Central Observatory Elem. Sci. Anthr. 2023 11 00048 10.1525/elementa.2022.00048
Itkin, P. et al. Sea ice and snow characteristics from year-long transects at the MOSAiC Central Observatory. Elem. Sci. Anthr. 11, 00048 (2023).10.1525/elementa.2022.00048
32. Nicolaus M Katlein C Maslanik J Hendricks S Changes in Arctic sea ice result in increasing light transmittance and absorption Geophys. Res. Lett. 2012 39 L24501 10.1029/2012GL053738
Nicolaus, M., Katlein, C., Maslanik, J. & Hendricks, S. Changes in Arctic sea ice result in increasing light transmittance and absorption. Geophys. Res. Lett. 39, L24501 (2012).10.1029/2012GL053738
33. Arrigo KR van Dijken G Pabi S Impact of a shrinking Arctic ice cover on marine primary production Geophys. Res. Lett. 2008 35 L19603 10.1029/2008GL035028
Arrigo, K. R., van Dijken, G. & Pabi, S. Impact of a shrinking Arctic ice cover on marine primary production. Geophys. Res. Lett. 35, L19603 (2008).10.1029/2008GL035028
34. Joli, N. et al. Hypometabolism to survive the long polar night and subsequent successful return to light in the diatom Fragilariopsis cylindrus. New Phytol. 241, 2193–2208 (2023).
35. White AW Growth of two facultatively heterotrophic marine centric diatoms J. Phycol. 1974 10 292 300 10.1111/j.0022-3646.1974.00292.x
White, A. W. Growth of two facultatively heterotrophic marine centric diatoms. J. Phycol. 10, 292–300 (1974).10.1111/j.0022-3646.1974.00292.x
36. Oziel L Faster Atlantic currents drive poleward expansion of temperate phytoplankton in the Arctic Ocean Nat. Commun. 2020 11 1705 10.1038/s41467-020-15485-5 32249780
Oziel, L. et al. Faster Atlantic currents drive poleward expansion of temperate phytoplankton in the Arctic Ocean. Nat. Commun. 11, 1705 (2020).32249780 10.1038/s41467-020-15485-5
37. Boss, E. & Behrenfeld, M. In situ evaluation of the initiation of the North Atlantic phytoplankton bloom. Geophys. Res. Lett. 37, L18603 (2010).
38. Castellani G Losch M Lange BA Flores H Modeling Arctic sea-ice algae: physical drivers of spatial distribution and algae phenology J. Geophys. Res. Oceans 2017 122 7466 7487 10.1002/2017JC012828
Castellani, G., Losch, M., Lange, B. A. & Flores, H. Modeling Arctic sea-ice algae: physical drivers of spatial distribution and algae phenology. J. Geophys. Res. Oceans 122, 7466–7487 (2017).10.1002/2017JC012828
39. Søreide JE Leu E Berge J Graeve M Falk-Petersen S Timing of blooms, algal food quality and Calanus glacialis reproduction and growth in a changing Arctic Glob. Change Biol. 2010 16 3154 3163 10.1111/j.1365-2486.2010.02175.x
Søreide, J. E., Leu, E., Berge, J., Graeve, M. & Falk-Petersen, S. Timing of blooms, algal food quality and Calanus glacialis reproduction and growth in a changing Arctic. Glob. Change Biol. 16, 3154–3163 (2010).10.1111/j.1365-2486.2010.02175.x
40. Assmy P Leads in Arctic pack ice enable early phytoplankton blooms below snow-covered sea ice Sci. Rep. 2017 7 40850 10.1038/srep40850 28102329
Assmy, P. et al. Leads in Arctic pack ice enable early phytoplankton blooms below snow-covered sea ice. Sci. Rep. 7, 40850 (2017).28102329 10.1038/srep40850
41. Nöthig E-M Summertime chlorophyll a and particulate organic carbon standing stocks in surface waters of the Fram Strait and the Arctic Ocean (1991–2015) Front. Mar. Sci. 2020 7 350 10.3389/fmars.2020.00350
Nöthig, E.-M. et al. Summertime chlorophyll a and particulate organic carbon standing stocks in surface waters of the Fram Strait and the Arctic Ocean (1991–2015). Front. Mar. Sci. 7, 350 (2020).10.3389/fmars.2020.00350
42. Zhang F He J Lin L Jin H Dominance of picophytoplankton in the newly open surface water of the central Arctic Ocean Polar Biol. 2015 38 1081 1089 10.1007/s00300-015-1662-7
Zhang, F., He, J., Lin, L. & Jin, H. Dominance of picophytoplankton in the newly open surface water of the central Arctic Ocean. Polar Biol. 38, 1081–1089 (2015).10.1007/s00300-015-1662-7
43. Ehrlich J Sea-ice associated carbon flux in Arctic spring Elem. Sci. Anthr. 2021 9 00169 10.1525/elementa.2020.00169
Ehrlich, J. et al. Sea-ice associated carbon flux in Arctic spring. Elem. Sci. Anthr. 9, 00169 (2021).10.1525/elementa.2020.00169
44. Lund-Hansen LC Effects of sea-ice light attenuation and CDOM absorption in the water below the Eurasian sector of central Arctic Ocean (>88°N) Polar Res 2015 34 23978 10.3402/polar.v34.23978
Lund-Hansen, L. C. et al. Effects of sea-ice light attenuation and CDOM absorption in the water below the Eurasian sector of central Arctic Ocean (>88°N). Polar Res 34, 23978 (2015).10.3402/polar.v34.23978
45. Ardyna, M. & Arrigo, K. R. Phytoplankton dynamics in a changing Arctic Ocean. Nat. Clim. Change10.1038/s41558-020-0905-y (2020).
46. Tremblay J-É Global and regional drivers of nutrient supply, primary production and CO2 drawdown in the changing Arctic Ocean Prog. Oceanogr. 2015 139 171 196 10.1016/j.pocean.2015.08.009
Tremblay, J.-É. et al. Global and regional drivers of nutrient supply, primary production and CO2 drawdown in the changing Arctic Ocean. Prog. Oceanogr. 139, 171–196 (2015).10.1016/j.pocean.2015.08.009
47. Olli K The fate of production in the central Arctic Ocean—top–down regulation by zooplankton expatriates? Prog. Oceanogr. 2007 72 84 113 10.1016/j.pocean.2006.08.002
Olli, K. et al. The fate of production in the central Arctic Ocean—top–down regulation by zooplankton expatriates? Prog. Oceanogr. 72, 84–113 (2007).10.1016/j.pocean.2006.08.002
48. Buesseler KO Revisiting carbon flux through the ocean’s twilight zone Science 2007 316 567 570 10.1126/science.1137959 17463282
Buesseler, K. O. et al. Revisiting carbon flux through the ocean’s twilight zone. Science 316, 567–570 (2007).17463282 10.1126/science.1137959
49. Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung. Polar research and supply vessel POLARSTERN operated by the Alfred-Wegener-Institute. J. Large Scale Res. Facil. 3, A119 (2017).
50. Garrison DL Buck KR Organism losses during ice melting: a serious bias in sea ice community studies Polar Biol. 1986 6 237 239 10.1007/BF00443401
Garrison, D. L. & Buck, K. R. Organism losses during ice melting: a serious bias in sea ice community studies. Polar Biol. 6, 237–239 (1986).10.1007/BF00443401
51. Knap, A., Michaels, A., Close, A., Ducklow, H. & Dickson, A. (eds). Protocols for the Joint Global Ocean Flux Study (JGOFS) Core Measurements. 170 (UNESCO-IOC, 1996).
52. Ulfsbo, A. et al. Dissolved inorganic carbon and total alkalinity of seawater samples during RV POLARSTERN expedition PS122 – MOSAiC. PANGAEA [data] 10.1594/PANGAEA.954969 (2023).
53. Nielsen ES The interaction of photosynthesis and respiration and its importance for the determination of 14C-discrimination in photosynthesis Physiol. Plant. 1955 8 945 953 10.1111/j.1399-3054.1955.tb07790.x
Nielsen, E. S. The interaction of photosynthesis and respiration and its importance for the determination of 14C-discrimination in photosynthesis. Physiol. Plant. 8, 945–953 (1955).10.1111/j.1399-3054.1955.tb07790.x
54. Behrenfeld MJ Boss ES Resurrecting the ecological underpinnings of ocean plankton blooms Annu. Rev. Mar. Sci. 2014 6 167 194 10.1146/annurev-marine-052913-021325
Behrenfeld, M. J. & Boss, E. S. Resurrecting the ecological underpinnings of ocean plankton blooms. Annu. Rev. Mar. Sci. 6, 167–194 (2014).10.1146/annurev-marine-052913-021325
55. R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria https://ropensci.org/blog/2021/11/16/how-to-cite-r-and-r-packages/ (2021).
56. Anderson, P. et al Upwelling and downwelling visible radiation measurements of the autonomous ice-tethered OptiCAL ‘hh’ buoy deployed during MOSAiC in the Dark site Second Year Ice. PANGAEA [data] 10.1594/PANGAEA.955045 (2023).
57. Duntley SQ Light in the sea J. Opt. Soc. Am. 1963 53 214 233 10.1364/JOSA.53.000214
Duntley, S. Q. Light in the sea. J. Opt. Soc. Am. 53, 214–233 (1963).10.1364/JOSA.53.000214
58. Saulquin, B. et al. Estimation of the diffuse attenuation coefficient KdPAR using MERIS and application to seabed habitat mapping. Remote Sens. Environ. 128, 224–233 (2013).
59. Itkin, P. et al. Magnaprobe snow and melt pond depth measurements from the 2019-2020 MOSAiC expedition. PANGAEA [data] 10.1594/PANGAEA.937781 (2021).
60. Nicolaus, M. et al. Snow height on sea ice, meteorological conditions and drift of sea ice from autonomous measurements from buoy 2019S96, deployed during MOSAiC 2019/20. Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven PANGAEA [data] 10.1594/PANGAEA.925326 (2020).
61. Maykut GA Untersteiner N Some results from a time-dependent thermodynamic model of sea ice J. Geophys. Res. 1896-1977 1971 76 1550 1575
Maykut, G. A. & Untersteiner, N. Some results from a time-dependent thermodynamic model of sea ice. J. Geophys. Res. 1896-1977 76, 1550–1575 (1971).
62. Nicolaus, M., Anhaus, P., Hoppmann, M., Tao, R. & Katlein, C. Spectral radiation fluxes, albedo and transmittance from autonomous measurement from Radiation Station 2020R11, deployed during MOSAiC 2019/20. Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, PANGAEA [data] 10.1594/PANGAEA.948828 (2023).
63. Fuchs, N. et al. Measurements of in-ice light profiles and optical properties of Arctic sea ice on MOSAiC 2020: Unified NetCDF data from the lightharp and lightchain instruments. PANGAEA [data] 10.1594/PANGAEA.963743 (2024).
64. Roche RBH King MD Calculations of Arctic ice-ocean interface photosynthetically active radiation (PAR) transmittance values Earth Space Sci. 2024 11 e2023EA002948 10.1029/2023EA002948
Roche, R. B. H. & King, M. D. Calculations of Arctic ice-ocean interface photosynthetically active radiation (PAR) transmittance values. Earth Space Sci. 11, e2023EA002948 (2024).10.1029/2023EA002948
65. Fuchs N In-ice light measurements during the MOSAiC expedition Sci. Data 2024 11 702 10.1038/s41597-024-03472-0 38937491
Fuchs, N. et al. In-ice light measurements during the MOSAiC expedition. Sci. Data 11, 702 (2024).38937491 10.1038/s41597-024-03472-0
66. Kanzow, T. Links to master tracks in different resolutions of POLARSTERN cruise PS122/3, Arctic Ocean - Longyearbyen, 2020-02-24–2020-06-04 (Version 2). Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven. PANGAEA [data] 10.1594/PANGAEA.924681 2020.
67. Anderson, P. et al. Upwelling and downwelling visible radiation measurements of the autonomous ice-tethered OptiCAL ‘gg’ buoy deployed during MOSAiC in the Monster Bay area. PANGAEA [data] 10.1594/PANGAEA.954849 (2023).
68. Anderson, P. et al. Upwelling and downwelling visible radiation measurements of the autonomous ice-tethered OptiCAL ‘ee’ buoy deployed during MOSAiC in the Sea Ice Ridge Observatory area. PANGAEA [data] 10.1594/PANGAEA.928495 (2023).
69. Hoppe, C. J. M. et al. Year-round discrete underway water column Chlorophyll a concentrations from the central Arctic. PANGAEA [data] 10.1594/PANGAEA.962597 (2023).
70. Hoppe, C. J. M. et al. Water column Chlorophyll a concentrations during the MOSAiC expedition (PS122) in the Central Arctic Ocean 2019-2020. PANGAEA [data] 10.1594/PANGAEA.963277 (2023).
71. Müller, O. et al. Flow cytometry dataset from CTD casts showing the abundance of microorganisms (smaller than 20 µm) during the Arctic MOSAiC expedition. PANGAEA [data] 10.1594/PANGAEA.963430 (2023).
72. Müller, O. et al. Flow cytometry dataset from first year sea ice (FYI) core bottom 5 cm sections showing the abundance of microorganisms (<20 µm) during leg 2, 3 (February, March and April 2020) of the Arctic MOSAiC expedition. PANGAEA [data] 10.1594/PANGAEA.963560 (2023).
73. Kraberg, A. C. Abundance of the main protist groups in a subset of microscopy samples from leg 3 of the MOSAiC campaign. PANGAEA [data] 10.1594/PANGAEA.965913 (2024).
74. Fuchs, N. Derivation of under-ice, mixed layer PAR values from profile measurements during MOSAiC. Zenodo [Code]. 10.5281/zenodo.12772364 (2024).
