Observations of late-winter marine microbial activity in an ice-covered fjord, west Greenland

Direct observations of marine microbial metabolism are sparse in the Arctic, particularly under sea ice during winter. This paper presents the first observations of Arctic winter microbial activity under sea ice in a west Greenland fjord (Lillefjord, ∼ 70 N). Here, measured changes in dissolved oxygen (DO) content in light and dark in situ incubations were used to calculate net community productivity, respiration and photosynthesis rates. Data were collected at two fully ice covered sites during February 2013, shortly after the end of the polar night. Averaged over the full study period, dark incubations showed statistically significant decreases in DO of −0.36± 0.24 (near shore) and −0.09± 0.07 g O2 m−3 d−1 (fjord centre), indicating respiration rates that were 2–20 times greater than rates previously reported under sea ice in the Arctic. Meanwhile, a lack of significant evidence for photosynthesis suggests that the rate of photosynthesis – if it was occurring – was much lower than that of respiration. The data also show no significant evidence of a temporal trend in metabolism rates over the study period; however, ambient seawater DO increased significantly at the fjord centre (0.023± 0.013 g O2 m−3 d−1), possibly attributable to processes not occurring in the incubations (such as sea ice algal photosynthesis). These data may improve our understanding of microbial activity in the fjord during winter, and its contribution to Arctic ecosystems under present and future conditions. The data are archived at PANGAEA (https://doi.org/10.1594/PANGAEA.906332, Chandler and Mackie, 2019; https://doi.org/10.1594/PANGAEA.912677, Chandler and Mackie, 2020).


Introduction
There is increasing evidence for rapid climate change in the Arctic, with wide-reaching impacts in both terrestrial and marine environments (Wassmann et al., 2011;McMeans et al., 2013;Post et al., 2013;Comiso and Hall, 2014). The observed reduction in sea ice cover (duration, extent and/or thickness), and the corresponding increase in solar illumination in the upper layers of the Arctic Ocean, is of particular interest. While estimates of marine net primary productivity (NPP) based on satellite retrievals of chlorophyll a have shown a link between reductions in sea ice cover and increases in NPP across much of the Arctic during 1998-2009, details of the processes associated with this change and its effects on higher levels of the food chain remain uncertain (Hansen et al., 2003;Arrigo et al., 2008;Brown and Arrigo, 2012;Vancoppenolle et al., 2013).
The logistical challenges associated with making direct observations of Arctic marine microbial metabolism mean that very few field data are available with which to assess metabolism magnitudes and controlling factors (Matrai et al., 2013;Vaquer-Sunyer et al., 2013). Satellite retrievals of chlorophyll a can provide excellent temporal and spatial coverage for monitoring NPP but have significant limitations. The data processing algorithms depend on multiple assumptions that may not be justified or appropriate in all cases (Arrigo et al., 2008); for example there may not be a direct relationship between retrieved chlorophyll a concentration and NPP (Flynn et al., 2013); data are unavailable for ocean wa-Published by Copernicus Publications. 898 D. Chandler and S. Mackie: Observations of late-winter marine microbial activity ter under sea ice and for sea ice itself, where productivity can be significant (Gosselin et al., 1997). Furthermore, the spatial resolution is generally too coarse to resolve smaller-scale features such as fjords, where the combination of nutrient inputs and buoyant mixing driven by subglacial meltwater discharge from marine-terminating glaciers can stimulate particularly high levels of productivity (Meire et al., 2017). Field observations of biological processes are therefore extremely valuable, both for improving and validating the parameterisations used in satellite retrieval algorithms and for providing information that cannot be measured remotely (e.g. in regions too small to be resolved by current remote sensing methods; or where observations of the individual components of microbial metabolism, and/or the variability of these with depth, are required).
Two approaches are generally followed for quantifying microbial metabolism: first, measuring the dissolved oxygen (DO) content of seawater in situ (Pomeroy, 1997;Rysgaard et al., 2001;; or second, measuring changes in the concentration of chemical tracers in closed incubation experiments. The former method enables observations at high spatial and/or temporal resolution, but their interpretation is often challenging because the system is open: changes in oxygen concentration due to biological activity must be separated from those of physical processes such as mixing and air-water gas exchange. Monitoring ambient DO in this way only quantifies net community productivity (NCP). In the latter technique, changes in DO or radioisotope concentrations can be used to infer rates of biological processes (Smith, 1994(Smith, , 1995Gosselin et al., 1997;Rysgaard et al., 1999Rysgaard et al., , 2001Hill and Cota, 2005;Regaudie-de-Gioux and Duarte, 2010;Vaquer-Sunyer et al., 2013). This requires samples to be collected and incubated, potentially involving complex analytical procedures, and while in situ incubations are unlikely to fully replicate natural conditions, they allow for more controlled conditions. Comparison of simultaneous incubations of samples exposed to light and samples kept in the dark yields estimates of community respiration (CR, measured in the dark samples) and gross primary productivity (GPP, interpreted as the difference between the light and dark samples) in addition to NCP (light samples). Ideally, the two approaches are used in tandem Cottrell et al., 2006). A recent approach combined in situ chlorophyll a and irradiance observations with a numerical model to estimate NPP under pack ice (Assmy et al., 2017); while less direct than ambient DO measurements or incubations, this method can yield good spatial coverage, is not affected by oxygen exchange with the atmosphere, and provides detailed in situ observations that will help address the limitations of remotely sensed chlorophyll a observations noted above.
Observations based on the above field methods have shown that several physical factors, notably dissolved nutrients and irradiance, as well as biological factors such as species composition and abundance, are each likely to play important roles in different environments and seasons, leading to a diverse range of measured metabolism rates (see Table 1 in this paper and Table 3 in Vaquer-Sunyer et al., 2013). The contribution of sea ice algae to the Arctic Ocean's annual primary production has also been observed to vary widely, for example ranging from 2 %-57 % (mean 17 %) in summer 1994 (Gosselin et al., 1997) to less than 1 % in Young Sound (NE Greenland) in summer 2000 (Rysgaard et al., 2001). The different methods and sampling strategies that are implemented by different teams make like-for-like comparison of observational data difficult. Despite both this and the high variability of the measured processes, some patterns have emerged. Notably, so-called blooms have been observed at, or shortly after, the break-up of the sea ice in summer, when microbial populations, chlorophyll a concentrations and microbial metabolism (both GPP and CR) in the surface layers are seen to increase rapidly Belzille et al., 2008;Mikkelsen et al., 2008;Terrado et al., 2008;Vaquer-Sunyer et al., 2013;Assmy et al., 2017).
While most measurements of metabolism have been carried out in ice-free summer conditions, there is evidence that microbial populations persist in both the sea ice and surface waters throughout the polar night (Berge et al., 2015;Vader et al., 2015) and can respond within a few days to increases in illumination (Zhang et al., 1998). Metabolism measurements in the water column under continuous sea ice, of which there are very few (Table 1), have detected community respiration during the polar night  but have yielded mixed results (positive, negative and insignificant NCP) during spring and summer (Gosselin et al., 1997;Cottrell et al., 2006, Seuthe et al., 2011Vacquer-Sunyer et al., 2013). In Franklin Bay (70 • N), chlorophyll a concentrations in sea ice algae and in the upper 11 m of the water column started to increase in mid-February, despite the persistence of continuous sea ice cover up to 2 m thick (Belzille et al., 2008), demonstrating how increasing activity by primary producers sometimes begins even under thick ice as daylight returns, well before ice break-up.
Models have predicted a strong ecological response to changing sea ice conditions along Greenland's west coast (Hansen et al., 2003), yet there are very few direct observations from the fjords that dominate Greenland's coastline (Rysgaard et al., 1999(Rysgaard et al., , 2001Mikkelsen et al., 2008;Matrai et al., 2013). Although only accounting for a small fraction of the total sea surface area in the Arctic, fjord waters have the potential to make a disproportionately strong contribution to Arctic marine productivity. This is partly due to the extensive area of shallow water along the long fjord coastlines, where benthic production can be important (Glud et al., 2002;Attard et al., 2014), and partly to the large nutrient fluxes transported to the fjords in meltwater runoff from the Greenland Ice Sheet (Hawkings et al., 2014;Lawson et al., 2014;Meire et al., 2017).
There is a need for more observations of microbial metabolism in Greenland's fjords and under sea ice. Such Table 1. Measurements of microbial metabolism under Arctic sea ice. Net community production (NCP) and gross primary production (GPP) are positive if oxygen is being released. Community respiration (CR) is positive if oxygen is being consumed.   Table 2 in this paper.
measurements will allow us to better understand marine productivity and quantify its contribution to Arctic marine ecosystems. Here we present in situ observations of microbial metabolism made under continuous sea ice cover at ∼ 70 • N in a west Greenland fjord (Lillefjord), derived from changes in DO measured in incubation experiments and in ambient seawater during February-March 2013, shortly after the transition from polar night to spring conditions on 21 January.

Field site and methods
Measurements were made in Lillefjord, west Greenland (70 • 30 N, 50 • 40 W). Lillefjord is 16 km long branch of the Uummannaq Fjord system, which opens to Baffin Bay approximately 70 km from the field site. The fjord system (including Lillefjord itself) receives meltwater runoff and calving icebergs from several outlet glaciers that drain the Greenland Ice Sheet, in common with many similar fjords in Greenland. In the winter of 2012/2013, continuous sea ice in Lillefjord had not formed until late January, which, although similar to several immediately preceding winters, was con- and fjord centre (FC) study sites. Coastline data are from https://www.soest.hawaii.edu/pwessel/gshhg/ (last access: 17 April 2020, Wessel and Smith, 1996). sidered locally to be unusually late (fishermen in Uummannaq Fjord, personal communication, 2013). Data were collected at two sites approximately mid-way between the calving front at the head of Lillefjord and the confluence of Lillefjord with Uummannaq Fjord (Fig. 1). Hole 1 (fjord edge) was approximately 50 m from the shore, in water 5-10 m deep; hole 2 (fjord centre) was located centrally in the fjord in water ∼ 300 m deep (Jakobsson et al., 2012). The sea ice thickness was initially measured as approximately 27 cm at both sites and increased slightly (by less than 10 cm) during the study period. Snow was absent from the sea ice until 13 February and then present in variable amounts thereafter (changes in these conditions are reported in Table 2).
Rates of photosynthesis and respiration were quantified using in situ incubation experiments in the uppermost ∼ 30 cm of the water column under the sea ice, based on measured changes in the DO content of seawater samples. A total of 13 experiments were carried out between 6 February and 6 March 2013. In each experiment, up to 10 samples of seawater were collected and incubated in situ under the sea ice in 250 mL biological oxygen demand (BOD) bottles. Half the bottles were wrapped in tin foil to make them opaque to light (dark bottles), and the remainder were left unwrapped and transparent to light (light bottles). It was assumed that no photosynthesis took place in the dark bottles, so any changes in DO between the start and end of the experiment ( DO dark ) are attributed solely to community respiration (CR). Both respiration and photosynthesis can occur in light bottles, so the change in DO ( DO light ) is assumed to indicate net community production (NCP). Rates of gross primary productivity (GPP), inferred to be photosynthesis, are estimated using the difference in DO between the light and dark bottles, i.e. NCP-CR. This is a standard and well-established method for measuring rates of microbial metabolism in freshwater and marine ecosystems Cottrell et al., 2006;Vaquer-Sunyer et al., 2013).
To begin each experiment, a hole of approximately 30 cm diameter was cut in the sea ice, using hand tools to avoid oil contamination. Water salinity and temperature were measured using a WTW handheld electrical conductivity (EC) meter (manufacturer's stated accuracy: temperature ±0.1 • C; EC ± 0.5 %), and the approximate ice thickness and overlying undisturbed snow depth were measured using a ruler. Both the ice thickness and snow depth were disturbed by the opening and reopening of the hole, so the measurements made at the start of each experiment should only be interpreted as indicative of the general ambient conditions. A metal sieve was used throughout sample collection to remove ice debris from the water surface in the hole, to prevent ice fragments from entering the sample bottles. Due to the typically cold air temperatures (−25 to −5 • C), the bottles were kept warm before use by adding ∼ 20 mL of boiling seawater to each bottle prior to transport to the field site; the bottles were then kept in an insulated box until needed. This was important to avoid the seawater freezing directly onto the cold Table 2. Summary of results obtained in the incubation experiments at the two holes FE (fjord edge) and FC (fjord centre). Quoted values are means with 95 % confidence intervals, except when N = 2 (as indicated by * ), where the error bounds are simply the range of the two observations. T is the duration of the incubation, N is the number of samples, and L − D is the difference between the change in dissolved oxygen in the light bottles, and in the dark bottles ( DO light − DO dark ).

Start time
Hole glass, which could have caused formation of ice inside the bottle or compromised the seal around the stopper. Immediately prior to sampling, each bottle was rinsed three times with seawater taken from the hole. The bottles were then refilled with water from the hole and suspended just under the water surface (to prevent ice from forming on the inside of the bottle). The water temperature and DO content in the bottle were measured using a PreSens Fibox3 fibre-optic oxygen meter (manufacturer's stated accuracy: ±1 %), which outputs data every 1 s. To measure the DO in each bottle, the sensor was allowed to stabilise (normally within 60 s), and readings were then taken for a further 20 s. The mean of these readings was recorded as the initial DO for the sample. After making the measurements, the bottle was immediately sealed with a glass stopper. The stoppers are buoyant, so a small piece of tin foil was wrapped over the stopper to keep it in place. The Fibox3 sensor control unit often stops functioning at cold temperatures, so it was kept warm in the insulated box with the preheated bottles. Bottles were checked carefully once filled and sealed to ensure that no air bubbles were present. They were then left suspended on nylon ropes approximately 50 cm below the ice surface. Bottles were left in place for periods of 24-96 h. After the allocated time, the hole was carefully reopened and the bottles retrieved and transferred to an insulated box for transport back to the field base. It was not practical to make the final DO measurements at the incubation site because, after being removed from the sea, the water in the bottles would have started to freeze in the time taken to record the measurements. Therefore, the sealed bottles were transported in an insulated carrier to the field base, where the final DO and temperature were measured using the same Fi-box3 sensor and probe. As for the initial DO measurements, the mean of readings made for 20 s after the sensor had stabilised were recorded as the final DO for the sample. The time between extraction from the hole and DO measurement was approximately 1 h, and temperature data showed the water temperature in the bottles to have increased by less than 2 • C between removal from the incubation site and completion of the last measurement at the field base. The main problems encountered during sample installation and recovery were associated with the cold air temperatures, which sometimes caused equipment failure (Fibox sys-D. Chandler and S. Mackie: Observations of late-winter marine microbial activity tem and/or netbook) or caused ice crystals to form in bottles. On 8 February we were interrupted when filling the bottles, by a calving event which threatened to cause break-up of the sea ice. Therefore, we were not always able to obtain results from the full set of 10 bottles (as indicated by N L and N D in Table 2).
Incubation times were initially chosen as 24 h, corresponding to one diurnal cycle. Since we were finding high variability in the early incubations, we also carried out some longer incubations (2 or 4 d) later in the study period. While these longer incubations allowed for potentially greater changes in DO (i.e. a lower signal-to-noise ratio), we also note that longer incubations increase the effects of methodological artefacts associated with the incubations being a closed rather than open system. We do not have enough data to confidently assess what the optimum incubation time would be.
When converting between oxygen demand and carbon storage in Table 1, a 1 : 1 stoichiometric ratio between CO 2 and O 2 was used, although this ratio is noted to be subject to some uncertainty (Telling et al., 2010). We note that interpretation of the measurements would benefit from simultaneous measurements of microbial biomass; however, facilities for measuring biomass were not available at the field site.

Uncertainty calculations
After each experiment, the difference between initial and final DO was calculated for each bottle, and the mean (µ) and standard deviation (σ ) of the differences were used to infer the change in DO ( DO). The number of bottles (n) in any one experiment was small, so it is appropriate to use a t distribution when calculating the 95 % confidence interval for DO. Treating light and dark bottles separately, DO was divided by the incubation time (T ) to give the rate of change in DO, DO/T as shown in Eq. (1), where t is the critical value of the t distribution at the 95 % confidence level.
The confidence intervals for DO light and DO dark were propagated through the calculations for rates of NCP, CR and photosynthesis. Each of these rates is therefore reported with an uncertainty corresponding to the limits of the 95 % confidence interval and is considered significant if zero lies outside the interval.
The mean and standard deviation of the initial DO measured in all the bottles (light and dark) were used to quantify the ambient DO and associated 95 % confidence interval for the seawater at each experiment start time, again using the t distribution as in Eq. (1). A linear fit was then applied to the time series of ambient DO at each study site, using linear least-squares regression. The gradient of the fit represents the mean rate of change in ambient DO over the study period and is reported with the 95 % confidence interval calculated using the t distribution and regression parameters.

Results
Throughout the study period, the seawater temperature and salinity varied between −1.5 and −1.7 • C and between 32.6 and 32.8 psu, respectively. The mean ±1σ DO of ambient seawater was 12.24 ± 0.23 g O 2 m −3 (fjord edge) and 12.38 ± 0.20 g O 2 m −3 (fjord centre). Linear regression analysis yielded no significant change in DO with time during the study period at the fjord edge (−0.001 ± 0.031 g O 2 m −3 d −1 ), while at the fjord centre there was a statistically significant increase in DO at a rate of 0.023 ± 0.013 g O 2 m −3 d −1 , equivalent to 720 ± 410 nM O 2 d −1 (Fig. 2a).
For each incubation experiment, the changes in DO ( DO) measured for each of the individual light and dark bottles were averaged to give a mean DO light and DO dark for the experiment. For three out of the five incubation experiments at the fjord edge, and three out of the eight experiments at the fjord centre, DO light showed a significant decrease (Table 2, Fig. 2b). DO dark showed a significant decrease for two out of four experiments at the fjord edge and for two out of seven experiments at the fjord centre (Fig. 2c).
Results for all the bottles in all the incubations were grouped together (averaging the rates of change in DO for all light bottles and all dark bottles separately) to reflect mean conditions over the whole study period. This was done separately for the two study sites. We found a significant decrease in DO for the dark bottles at the fjord edge and for both the light and dark bottles at the fjord centre (−0.36 ± 0.24 g O 2 m −3 d −1 for dark bottles at the fjord edge; −0.10 ± 0.07 and −0.09 ± 0.07 g O 2 m −3 d −1 for light and dark bottles at the fjord centre). There was no significant change in DO for light bottles at the fjord edge.
For each incubation experiment, the difference between DO calculated for the light and for the dark bottles is interpreted as GPP. None were found to be significantly different from zero ( Fig. 2d) except for the first experiment at the fjord edge, where the 95 % confidence interval for GPP was 0.31 ± 0.28 g O 2 m −3 d −1 .

Discussion
Significant decreases in DO in the dark incubations at both sites are attributed to microbial respiration ( Fig. 2c and Table 2). This is consistent with the few previous observations of microbial metabolism under sea ice cover, which have found significant CR (Table 1), and is not unexpected given the persistence of microbial communities through the polar night (Berge et al., 2015); however, in Lillefjord the measured rates (particularly at the fjord edge) are considerably higher than those at other ice-covered sites (  The rate of change in DO in the light bottles (considered to represent NCP) was either weakly negative or insignificant. These results for NCP at Lillefjord are consistent with observations from the one other west Greenland fjord studied during February-March, to our knowledge (NCP < 0.001 g O 2 m −3 d −1 at Kangerluarsunnguaq in sub-Arctic SW Greenland; Mikkelsen et al., 2008). They are also consistent with some studies of ice-covered open-ocean sites Hill and Cota, 2005) but contrast with others (weakly positive NCP was reported by Cottrell et al., 2006, in the Arctic Ocean andby Vaquer-Sunyer et al., 2013, in the Fram Strait). With the exception of , measurements at these ice-covered ocean sites were collected later in the year (mid-April to June) than those at Lillefjord.
The lack of any significant difference between DO in the light and dark bottles means that there was no significant evidence for GPP (or photosynthesis). This should not be interpreted as significant evidence for no photosynthesis, particularly given the high variance in the data indicated by wide 95 % confidence intervals; however, it does show that the rate of photosynthesis -if it was occurring -must have been much smaller than that of respiration. For comparison, Rysgaard et al. (1999) and Mikkelsen et al. (2008) both found evidence for very low rates of photosynthesis under sea ice in fjords shortly before ice breakup (GPP < 0.003 g O 2 m −3 d −1 in Young Sound, Rysgaard et al., 1999; NCP = +0.001 g O 2 m −3 d −1 in Kangerluarsunnguaq, Mikkelsen et al., 2008). Other studies have found evidence for significant GPP in the largely (> 80 %) icecovered open ocean in the Fram Strait during April-May, which contributed to an overall positive NCP (Seuthe et al., 2011;Vaquer-Sunyer et al., 2013). Under continuous firstyear sea ice in Franklin Bay (a coastal site at a similar latitude to Lillefjord), Terrado et al. (2008) observed an increase in the abundance of photosynthetic organisms as early as February in response to increasing surface irradiance. Similarly, chlorophyll a concentrations were observed to reach a minimum in January and to begin increasing in February, within first-year sea ice in the open Arctic Ocean in the Canada Basin (Melnikov et al., 2002). In Lillefjord, it is not clear whether the photosynthetic activity in the surface waters had not commenced during the study period or was masked by the stronger and highly variable respiration signal.
In contrast to the incubation results, there was a significant increase in ambient DO at the fjord centre of 0.023 ± 0.013 g O 2 m −3 d −1 . This differs from the findings of , where a decrease was observed in ambient DO under sea ice during winter in the western Arctic Ocean. Lillefjord was completely ice covered during the study period, preventing any air-water gas exchange, and no decrease in DO was observed in the incubation experiments. Therefore, the increase in ambient DO may have been due to processes not occurring in the incubations (such as sea algal photosynthesis on the underside of the sea ice). These contrasting results from simultaneous incubation and in situ experiments demonstrate the advantage of using both closed and open techniques when there is continuous ice cover.
Despite the increasing surface irradiance (longer daylight hours and less shading by surrounding topography at higher solar elevations), the incubation experiments provide no evidence for temporal changes in metabolism rates. It is possible that some of the increase in incident radiation at the snow/ice surface did not reach the water below the ice because of increases in snow cover and ice thickness over this same period. Without under-ice irradiance measurements this is necessarily uncertain; however, the radiation intensity S reaching the water column (as a fraction of surface incident radiation intensity S 0 ) can be estimated using S/S 0 = (1 − α) exp(−k s z s − k i z i ), where α is the surface albedo, and z s,i and k s,i are the thicknesses and extinction coefficients for snow and ice, respectively. Assuming extinction coefficients of 4.8 and 0.9 m −1 for snow and sea ice, and albedos of 0.90 and 0.65 for fresh snow and sea ice (following Mikkelsen et al., 2008), the under-ice irradiance is estimated as 31 % of the surface irradiance before snowfall on 14 February and 11 % afterwards. Therefore, increases in surface irradiance in early February as experienced under the ice could have been considerably reduced following snowfall in mid February.

Conclusions
These data provide a first indication of winter microbial metabolism beneath sea ice in an Arctic fjord in west Greenland. Thirteen in situ incubation experiments provide strong evidence for microbial respiration at rates 2-20 times higher than those reported under sea ice elsewhere in the Arctic (see Table 1). The high variance in the NCP and CR results (both between individual bottles in one experiment and between incubations) is a common characteristic of marine microbial metabolism measurements under sea ice (Table 1) and presents a challenge to accurate calculation of GPP or temporal trends. This variance should be carefully accounted for when considering uncertainties associated with estimates of the regional-scale contributions of microbial activity, which are necessarily based on the limited data that are currently available. In future studies this could be addressed by increasing the number of bottles and/or conducting more frequent experiments and by extending the study period to obtain a longer time series. Finally, the contrast between the increasing trend in ambient seawater DO and the net oxygen decrease in the incubation experiments highlights potential differences between controlled and open experiments. In this study, the difference is most likely attributable to net production by sea ice algae (which would increase DO in the ambient seawater), contrasting with net respiration in the underlying water (which would decrease the DO in the closed incubations). This suggests an earlier onset of photosynthesis at the underside of the sea ice than in the underlying water.
Author contributions. Both authors contributed equally to data collection, data analysis and preparation of the paper.
Competing interests. The authors declare that they have no conflict of interest.