Although studies of biogeochemical processes in polar sea ice have been increasing, similar research on relatively warm low-salinity sea ice remains sparse. In this study, we investigated biogeochemical properties of the landfast sea ice cover in the brackish Bothnian Bay (Northern Baltic Sea) and the possible role of this sea ice in mediating the exchange of greenhouse gases, including carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) across the water column–sea ice–atmosphere interface. Observations of total alkalinity and dissolved inorganic carbon in both landfast sea ice and the water column suggest that the carbonate system is mainly driven by salinity. While high CH4 and N2O concentrations were observed in both the water column (up to 14.3 and 17.5 nmol L–1, respectively) and the sea ice (up to 143.6 and 22.4 nmol L–1, respectively), these gases appear to be enriched in sea ice compared to the water column. This enrichment may be attributable to the sea ice formation process, which concentrates impurities within brine. As sea ice temperature and brine volume decrease, gas solubility decreases as well, promoting the formation of bubbles. Gas bubbles originating from underlying sediments may also be incorporated within the ice cover and contribute to the enrichment in sea ice. The fate of these greenhouse gases within the ice merits further research, as storage in this low-salinity seasonal sea ice is temporary.
The Baltic Sea is a brackish subpolar semi-enclosed sea, only connected to the world ocean through the shallow Danish straits, with surface salinity ranging from <1 to 9 (Granskog et al., 2006a). Each winter, sea ice forms and lasts from 5 to 7 months, with a maximum sea ice extent usually observed from mid-February to mid-March. The sea ice cover has large interannual variability in both dates of freezing and break up and its extent and thickness (Granskog et al., 2006a). During mild winters, sea ice only forms in coastal areas and the northernmost basin of the Baltic Sea, the Bothnian Bay (BB), whereas during severe winters, most of the Baltic Sea is ice-covered (Vihma and Haapala, 2015). Due to the brackish nature of its parent seawater, sea ice from the Baltic Sea is characterized by lower bulk ice salinity (usually <2; Granskog et al., 2003) and smaller brine channels (reflecting the brine volume), compared to its polar counterparts. However, while sea ice from the Baltic Sea may be structurally similar to sea ice formed in polar waters, its low salinity (and brine volume content) strongly affects the physical and chemical properties of the ice, with further consequences for biology, remote sensing, and transfer of solutes (e.g., salts, nutrients, gases) across the sea ice–seawater interface. The contribution of meteoric ice (i.e., precipitation transformed into ice) to the total thickness of landfast sea ice in the Baltic can be substantial (Kawamura et al., 2001; Granskog et al., 2006a), as meteoric ice formation is a recurring process in the Baltic Sea. The relatively thin sea ice cover in the Baltic Sea can easily be depressed below seawater level by a small amount of snow deposited onto the ice, resulting in surface sea ice flooding, potentially promoting snow-ice formation. In addition, the mild conditions are favorable for superimposed ice formation by successive melting and refreezing cycles of snow (Granskog et al., 2006a). Depending on the season and year, the meteoric ice fraction has been estimated to contribute up to 35% of the total ice thickness (Granskog et al., 2004).
As sea ice forms, brine channels and pockets concentrate impurities such as salts, gases, and inorganic and organic components (Eicken, 2003). Once concentrated within sea ice, gas dynamics and fluxes are largely controlled by sea ice temperature and salinity through brine saturation state and sea ice permeability. The resulting chemical activity of these impurities within sea ice influences the exchange of greenhouse gases (GHGs), including carbon dioxide (CO2), methane (CH4), and N2O (notable GHGs), across the ocean–sea ice–atmosphere interface (Lannuzel et al., 2020). During winter, sea ice acts as a source of CO2 with high brine partial pressure of CO2 due to increased brine concentrations and associated ikaite precipitation (Rysgaard et al., 2011; Geilfus et al., 2013; Fransson et al., 2015; Geilfus et al., 2016), whereas in spring and summer, sea ice shifts to become a sink for atmospheric CO2 due to brine freshening, ikaite dissolution, and the biological carbon pump (Rysgaard et al., 2011; Geilfus et al., 2012; Van der Linden et al., 2020).
Relative to CO2, the impact of sea ice on CH4 exchange across the ocean–sea ice–atmosphere interface is less well understood. Sea ice has been suggested to play a key role on the regulation of CH4 concentrations in the Arctic. In winter, CH4 has been reported to accumulate either within the sea ice matrix (Shakhova et al., 2010; Crabeck et al., 2014) or at its interface with the water column (Kvenvolden et al., 1993; Kitidis et al., 2010). As sea ice temperature increases, CH4 could be released, either directly to the atmosphere (Shakhova et al., 2010) or to the underlying seawater, along with the brine drainage (Damm et al., 2015). Sea ice melt decreases the saturation state of surface water, resulting in a drawdown of atmospheric CH4 (Heeschen et al., 2005). However, sea ice can also enhance CH4 oxidation, potentially reducing CH4 effluxes to the atmosphere (Tison et al., 2017), particularly above continental shelves where sediment represents the main source of CH4 to the ocean (Arctic Monitoring and Assessment Programme, 2015), and significant ebulliative fluxes are reported (Shakhova et al., 2010). The major CH4 sources in marine systems are microbial methanogenesis in anaerobic sediments and abiotic geochemical processes (Reeburgh, 2007). Usually a large fraction (up to 90%) of released CH4 is oxidized in microbially mediated processes before it reaches the upper water layers (Reeburg, 2007; Knittel and Boetius, 2009). In addition to anaerobic methanogenesis, CH4 can also be produced in aerobic conditions via methylotrophic methanogenesis when bacteria use methylphosphonate as their phosphate source, releasing CH4 (Karl et al., 2008; Carini et al., 2014). Highest CH4 emissions occur in coastal areas of high organic matter input from rivers (Humborg et al., 2019). In the Baltic Sea, lateral transport of CH4 via riverine flow can also significantly contribute to high coastal concentrations and CH4 supersaturation (Myllykangas et al., 2021). Eutrophication and global warming are suggested to increase CH4 release to the atmosphere (Borges et al., 2016; Humborg et al., 2019).
N2O production, mainly occurring in sediments, is regulated by the presence of oxygen and nitrogen availability. Microbial N2O production is complex and takes place in multiple oxidative and reductive processes of the nitrogen cycle (Anderson and Levine, 1986; Bakker et al., 2014). Estuaries and coastal areas are known to be areas for high N2O production due to high riverine loads of nitrogen (Bakker et al., 2014). Fewer studies have investigated N2O cycling in the Arctic, and the N2O cycle within sea ice is relatively unknown. Sea ice growth and melt (Kitidis et al., 2010; Randall et al., 2012) as well as denitrification and diazotrophy (Verdugo et al., 2016) have been suggested to affect the N2O cycle and the exchanges of N2O with both the atmosphere and the water column (Kaartokallio, 2001). The sea ice bacterial community harbors potentially denitrifying organisms (Deming and Collins, 2017); however, identification of active denitrifying microbes and their prevalence remain areas of ongoing research.
GHGs are produced and/or consumed by diverse assemblages of auto- and heterotrophic organisms living in or under the ice during the annual cycle of sea ice growth and melt (Arrigo, 2017; Bluhm et al., 2017; Caron et al., 2017; Deming and Collins, 2017; Lannuzel et al., 2020). Bacteria are the main group of organisms responsible for decomposing organic matter and remineralizing the nutrients, and their growth is regulated by multiple simultaneously acting environmental factors, such as salinity, temperature, and nutrients (Pomeroy and Wiebe, 2001; Kuosa and Kaartokallio, 2006), as well as food web interactions, such as protistan grazing (Kaartokallio, 2004; Riedel et al., 2007). Baltic and Arctic sea ice bacterial communities harbor potential methane-oxidizing bacteria, such as Methyloacidophilaceae (E Eronen-Rasimus, personal communication, September 2021; Op den Camp et al., 2009) and the genus Methylomicrobium (Bowman et al., 2012; Puri et al., 2015); however, their abundance is low, and their methane oxidation activity within the ice has not been verified.
While sea ice plays an active role in GHG exchange (Rysgaard et al., 2011; Geilfus et al., 2012), more extensive study is needed to determine the magnitudes and controls of that exchange (Parmentier et al., 2017; Lannuzel et al., 2020). Biogeochemical processes at low salinity and warm temperature potentially differ from those in polar seas but are absent from published data sets. To our knowledge, we present the first observations of GHG dynamics within sea ice in the BB and how sea ice could affect GHG exchanges across the ocean–sea ice–atmosphere interface. Due to the low salinity in the Baltic Sea, its sea ice could potentially be used as a model to investigate future changes in Arctic sea ice due to warming associated with climate change, especially on shallow shelves where high river inputs can significantly reduce sea surface salinity (Granskog et al., 2006a). Therefore, we investigated how landfast sea ice in the BB controls the exchange of GHGs across the ocean–sea ice–atmosphere and relate our observations to the current understanding of microbial production within the ice cover.
2.1. Study area
Sea ice and seawater samples were collected in the northernmost basin of the Baltic Sea, the BB. This shallow basin, with an average depth of 43 m, is separated from the deeper southern basin of the Gulf of Bothnia, the Bothnian Sea (BS), by the Northern Quark, a shallow sill of only 20 m depth (Figure 1). The brackish nature of the seawater in both basins is due to the high riverine discharge from large catchments and limited supply of saline water from lateral intrusions of surface water from the central Baltic Sea. Mean seawater surface salinity is 3.5 in the BB and slightly higher at 6.5 in the BS due to larger inflow of saline waters from the Baltic Proper (Löffler et al., 2012).
Samples were collected from 9 landfast sea ice stations from the Swedish and Finnish coasts (Figure 1; Table 1) from 18 to 25 February 2018. Sampling locations were chosen close to riverine inputs and major cities as part of a research program dedicated to microplastic incorporation within sea ice (Geilfus et al., 2019). In Luleå (Sweden) and Oulu (Finland), sea ice was also collected away from the river geographically (noted as Luleå* and Oulu*). Additional landfast sea ice (stations Ar-1 and Ar-2) and seawater samples (stations RR-7, F-13, and F-18) were collected onboard the R/V Aranda, from 3 to 5 February 2020, as part of the annual winter monitoring cruise COMBINE 1. All sea ice samples were collected in shallow areas, as samples from 2018 were collected by walking from the seashore while the water depths at stations Ar-1 and Ar-2 were 21 and 22 m, respectively.
|Date (dd/mm/yyyy) .||Stations .||Latitude (N) .||Longitude (W) .||Sea Ice Depth (cm) .||Snow Depth (cm) .||Freeboard (cm) .|
|Date (dd/mm/yyyy) .||Stations .||Latitude (N) .||Longitude (W) .||Sea Ice Depth (cm) .||Snow Depth (cm) .||Freeboard (cm) .|
aParameters not measured or not applicable.
2.2. Sampling procedure
Sea ice cores were collected using a MARK II coring system (internal diameter = 9 cm, Kovacs Enterprises®, Roseburg, OR, USA). Directly after extraction of the core, the temperature profile was measured in situ using a calibrated temperature probe (Testo 720®, ±0.1°C precision) inserted into predrilled holes (5-cm intervals), perpendicular to core sides. Cores were then cut into 5–10 cm sections, stored in sealed plastic containers, and vacuum-sealed in a gas-tight plastic bag (nylon/poly, Cabela’s, Sidney, NE, USA; Hu et al., 2018) directly in the field. Back in the laboratory, samples were melted in the dark at 4°C to minimize the possible dissolution of ikaite crystals. Once melted, aliquots of meltwater were transferred into 12-ml gas-tight vials (Exetainers, Labco, High Wycombe, United Kingdom) for total alkalinity (TA) and dissolved inorganic carbon (DIC) measurements, while CH4 and N2O samples were collected into 60-ml serum bottles. Samples were preserved by adding 20 and 60 µL, respectively, of a saturated mercuric chloride (HgCl2, Thermo Fisher [Kandel] GmbH, Germany) solution. Water aliquots (2 mL) were also collected in glass vials without head space for stable isotope composition of oxygen in water, δ18O-H2O, herein referred to as δ18O. Samples were then stored in the dark at room temperature until analysis.
Seawater samples were collected onboard R/V Aranda using a rosette sampling system equipped with 24 12-L Niskin-type bottles (Sea-Bird PVC sample bottle) and probes for conductivity, temperature, and depth (CTD Sea-Bird SBE-911+, calibrated by the manufacturer). Conductivity drift was checked and corrected, based on discrete water samples taken throughout the water column and analyzed onboard using a Guildline Autosal 8400 salinometer calibrated with standard seawater from the International Association for Physical Sciences of the Ocean. Samples for TA, DIC, CH4, N2O, and δ18O were collected, following the same procedure as for melted sea ice, immediately after the rosette was secured onboard.
2.3. Sample analysis
Bulk ice salinity was calculated based on conductivity (Grasshoff et al., 1983) measured using a meter (Orion 3-star, Thermo Scientific, Waltham, MA, USA) coupled with a conductivity cell (Orion 013610MD, Thermo Scientific). Brine volume was estimated from measurements of bulk ice salinity, temperature, and density according to Cox and Weeks (1983) for temperatures below –2°C and according to Leppäranta and Manninen (1988) for ice temperatures within the range of –2°C to 0°C.
TA was determined by Gran titration (Gran, 1952) using a TIM 840 titration system (Radiometer Analytical, ATS Scientific, Burlington, Ontario, Canada), consisting of a Ross sure-flow combination pH glass electrode (Orion 8172BNWP, Thermo Scientific) and a temperature probe (Radiometer Analytical, Lyon, France). A 12-mL sample was titrated with a standard 0.05 M HCl solution (Alfa Aesar, Ward Hill, MA, USA). DIC was measured on a DIC analyzer (Apollo SciTech, Newark, DE, USA) by acidification of a 0.75-mL subsample with 1 mL 10% H3PO4 (Sigma-Aldrich, Saint-Louis, MO, USA), and quantification of the released CO2 with a nondispersive infrared CO2 analyzer (LI-7000, LI-COR, Lincoln, NE). Results were then converted from µmol L–1 to µmol kg–1 based on sample density, which was estimated from salinity and temperature at the time of the analysis. Accuracies of ±3 and ±2 µmol kg–1 were determined for TA and DIC, respectively, from routine analysis of certified reference materials (AG Dickson, Scripps Institution of Oceanography, San Diego, CA, USA). CH4 and N2O concentrations were measured via the headspace equilibrium technique (25 ml N2 headspace in 60-ml serum bottles) and measured with a gas chromatograph (SRI 8610C) with flame ionization detection, and electron capture detection calibrated with CH4: CO2: N2O: N2 mixtures (Air Liquide, Paris, France) of 1, 10, and 30 ppm CH4 and 0.2, 2.0, and 6.0 ppm N2O. The measurement precision was 3.9% for CH4 and 3.1% for N2O. Values for δ18O were measured on an isotope analyzer (L2130-i, Picarro, Santa Clara, CA, USA) with a precision of 0.025‰, calculated as the standard deviation (SD) from 10 repeated measurements of standard materials.
3.1. Sea ice
Sea ice thickness ranged from 30 to 46 cm (Figure 2; Table 1) and exhibited a typical temperature gradient, with colder temperatures observed at the sea ice interface with the atmosphere (minimum of –7.1°C at Ar-2). Bulk ice salinity was typical for the Baltic Sea (Granskog et al., 2006a), with values lower than 1. Due to the low salinity, the brine volume content of the ice was below the permeability threshold (5%; Golden et al., 2007; Figure 2d). Sea ice δ18O ranged from –18.7‰ to –7.8‰. Lowest values were observed in the upper section of the ice cover. The δ18O increased downward throughout the ice cover and remained relatively constant toward the bottom of the cores (mean = –10.1‰, SD = 1.1, n = 22), except at Oulu* where the δ18O profile remained relatively constant throughout the entire ice thickness with averaged δ18O less depleted than at other stations (Oulu*: –8.3‰; Figure 2c).
Concentrations of TA and DIC ranged from 18 to 282 µmol kg–1 and from 12 to 267 µmol kg–1, respectively (Figure 3). These concentrations fall at the low end of previously reported values in Arctic sea ice (Rysgaard et al., 2007; Geilfus et al., 2015; Fransson et al., 2017) but are comparable to values reported from multiyear sea ice (Rysgaard et al., 2009; Geilfus et al., 2021) and expected from the low bulk ice salinity. Oulu* stands out of the general pattern with higher concentrations of both TA and DIC compared to the other stations. No general pattern could be observed in the TA and DIC profiles. The sea ice surface has either maximum concentrations of TA and DIC (Skellefteå, Kemi, and Oulu*) or minimum concentrations (Luleå, Umeå, Oulu, Ar-1, and Ar-2). Finally, concentrations at the Kalix and Vaasa sites appear homogeneous throughout the ice thickness.
Concentrations of CH4 ranged from 2.6 to 143.6 nmol L–1 (Figure 3), up to 10 times the sea ice CH4 concentration reported in a Greenlandic fjord (Crabeck et al., 2014) or in landfast sea ice in Barrow, Alaska (Zhou et al., 2014). Maximum concentrations were observed at the sea ice interface with the atmosphere at Skellefteå, Umeå, Oulu, Vaasa, and Ar-1, which contrasts with Luleå, Kalix, Kemi, and Oulu*, where CH4 profiles remained constant throughout the ice thickness. Bottom CH4 concentrations for all stations were in the same range (mean = 7.9 nmol L–1, SD = 3.6, n = 22), except at Vaasa where the highest concentration was observed (143.6 nmol L–1). N2O concentrations ranged from 3.2 to 22.4 nmol L–1, with the highest concentrations (>10 nmol L–1) observed at Vaasa. Local maxima were observed at the upper layer of the ice surface, except at Luleå* and Oulu*.
3.2. Water column
The water column of the BB was sampled at stations RR-7 and F-13, while the BS was sampled at station F-18 (Figure 1). In the BB, the water column was quite homogeneous, without a thermocline or halocline observed (Figure 4). Seawater temperature was <1°C (mean = 0.64, SD = 0.29, n = 16), and its salinity averaged 3.21 (SD = 0.09, n = 16). However, at F-13, seawater temperature and salinity increased slightly with depth from 0.8°C and 3.21°C, respectively, at the surface to 0.96°C and 3.42°C, respectively, at the bottom. At F-18, both a thermo- and halocline were observed at 40–50 m depth (Figure 4). The upper layer had an average temperature and salinity of 1.73°C (SD = 0.07, n = 7) and 5.29°C (SD = 0.02, n = 7), respectively, while the bottom layer averaged 3.27°C (SD = 0.16, n = 6) and 5.56°C (SD = 0.11, n = 6), respectively. Values for δ18O were stable throughout the water column at each station, with averaged values increasing from –10.0‰ at RR-7 (SD = 0.2, n = 5) to –9.4‰ at F-13 (SD = 0.1, n = 5), and –8.2‰ (SD = 0.1, n = 5) at F-18.
In the BB, TA and DIC concentrations were similar between stations, with TA ranging from 860 to 1002 µmol kg–1, and DIC ranging from 888 to 1031 µmol kg–1. In the BS (F-18), TA and DIC increased with concentrations ranging from 1302 to 1414 µmol kg–1 and from 1241 to 1441 µmol kg–1, respectively. Profiles of CH4 and N2O concentrations were quite similar between stations. CH4 concentrations in the upper 50 m were on average at 6.8 nmol L–1 (SD = 0.6, n = 6). At both F-13 and F-18, an increase in CH4 concentration (up to 14.3 nmol L–1) with depth was observed. Based on seawater temperature, salinity, CH4 concentration, and using constants from Wiesenburg and Guinasso (1979), we estimated that surface waters in both the BB and BS were supersaturated with CH4 (up to 176%) compared to the atmosphere. N2O concentrations were stable between stations with the concentration averaging 14.8 nmol L–1 (SD = 3.7, n = 13), with a low concentration of 3.1 nmol L–1 observed at the bottom of F-13. Based on constants from Weiss and Price (1980), we determined that the surface water was undersaturated in N2O compared to the atmosphere.
4.2. Sea ice properties
The sea ice conditions encountered during this survey were typical for the Baltic Sea (Granskog et al., 2006a), where the low bulk ice salinity is mainly due to the brackish nature of the parent seawater (Figure 4). However, the presence of snow and/or superimposed ice on the sea ice cover also potentially affects the bulk ice salinity. The strongly depleted δ18O values observed in the upper layer of the ice cover (Figure 2) could suggest the presence of meteoric ice in our samples. Meteoric water distributions from riverine sources were not differentiated due to the lack of endmember or ice texture analysis. The upper 10–30 cm of the ice cover exhibited more depleted δ18O values compared to the averaged δ18O of surface seawater in the BB (–9.75‰; dotted line, Figure 2c). This difference could indicate the presence of meteoric ice formation, either due to snow ice formation or surface flooding by local river runoff, except at Oulu* where the δ18O profile within sea ice was more homogeneous and less depleted (mean = –8.3‰) compared to the average surface seawater.
4.2. Biogeochemical processes within sea ice
The low TA and DIC concentrations within sea ice are expected due to low bulk ice salinity, as TA and DIC in both sea ice and seawater exhibited a strong correlation with salinity (black dotted line, Figure 5). Most of the TA: DIC data from both seawater and sea ice align on the slope of 1 (dashed black line, Figure 5c and f), except the upper sea ice layer in Luleå (blue square, Figure 5f). As no ikaite precipitates were observed during the sea ice melt, this alignment confirms a lack of precipitation of calcium carbonate within sea ice, which contrasts with polar sea ice where ikaite is commonly reported (Dieckmann et al., 2008; Geilfus et al., 2013; Rysgaard et al., 2013). Taking into account the low bulk sea ice salinity (<1) and related brine volume fraction below the sea ice permeability threshold (5%; Golden et al., 2007), sea ice in the BB appears to mimic the “closed system” of sea ice formation, unfavorable for ikaite precipitation (Papadimitriou et al., 2013). In such a system, rejection of excess CO2 from sea ice does not occur and thus the low pCO2 needed for ikaite to precipitate is not maintained (Papadimitriou et al., 2013). In addition, precipitation of calcium carbonate within the water column has been suggested to be negligible in the BB and BS due to undersaturation of both calcite and aragonite minerals in winter (Tyrrell et al., 2008).
Within sea ice, both TA and DIC exhibited large variability compared to their relationship with salinity (Figure 5d and e). Based on averaged surface seawater TA (871 µmol kg–1, SD = 13, n = 4) and DIC (905 µmol kg–1, SD = 10, n = 4) in the BB, we calculated the expected surface seawater TA and DIC at salinities observed within sea ice (from 0.1 to 0.9) considering that TA and DIC are conservative with salinity (red dashed line, Figure 5d and e). Most of the sea ice data exhibited higher TA and DIC than expected in seawater at similar salinity values, suggesting an enrichment in the ice cover compared to the underlying seawater. However, the sea ice data covered a large part of the BB, while seawater was only sampled twice in the BB (at RR-7 and F-13). RR-7 was sampled not too far away from Oulu, Oulu*, Ar-1, and Ar-2, and both RR-7 and F-13 exhibited similar properties (Figure 4). This sea ice enrichment was also observed for both CH4 and N2O (Figure 6). CH4 concentrations within sea ice were an order of magnitude higher than those observed within the water column (Figure 6a). To account for CH4 and N2O from parent seawater, we normalized the averaged surface seawater concentration of CH4 (6.8 nmol L–1, SD = 0.1, n = 3) and N2O (16.1 nmol L–1, SD = 1.0, n = 3), observed in the BB to a bulk ice salinity of 0.5 (noted as nCH4_sw and nN2O_sw). CH4 concentrations within sea ice were significantly higher than nCH4_sw (1.0 nmol L–1; red line, Figure 6c), suggesting a CH4 enrichment within sea ice compared to the underlying seawater. N2O concentrations within sea ice are the same order of magnitude as those in the water column (Figure 6b). However, nN2O_sw (2.5 nmol L–1) is lower than observed sea ice concentrations, also suggesting an enrichment of N2O within sea ice compared to the underlying seawater (Figure 6d).
As sea ice forms, it concentrates impurities, including gases, within the ice structure (Eicken, 2003). Therefore, high GHG concentrations within sea ice could originate from the parent seawater. If the carbonate system follows the dilution curve (Figure 5), then high CH4 and N2O concentrations are observed in the water column (Figure 4). CH4 supersaturation has been reported across the Baltic Sea (Bange et al., 1994; Schmale et al., 2010), with CH4 concentrations increasing with seawater depth (Figure 4), suggesting a release of CH4 from the sediment to the water column (Schmale et al., 2010). Another source of GHGs to the marine environment could be the waters discharged from rivers (Frankignoulle et al., 1998; Borges, 2005). In Liminganlahti Bay, south of Oulu, CO2 and CH4 supersaturation has been observed in river water during wintertime and has been associated with the lack of photosynthesis and the presence of river ice impeding gas exchange with the atmosphere (Silvennoinen et al., 2008). Apart from the river water layer, lowest CH4 concentrations were observed during winter, suggesting other controls than solubility (Myllykangas et al., 2021). However, riverine contributions to TA in the Baltic Sea are spatially variable and depend on the bedrock nature of the catchment area (Schneider, 2011). In the Gulf of Bothnia, riverine contributions to TA ranged from 140 to 229 µmol kg–1 (Hjalmarsson et al., 2008, and Beldowski et al., 2010, respectively). This small contribution is mainly due to igneous rock in the catchment areas that does not release weathering products contributing to TA (Löffler et al., 2012). From Figure 5, the contribution of river water (salinity = 0) in TA and DIC can be estimated at 45 and 41 µmol kg–1, respectively, suggesting a minor contribution of riverine water to the carbonate system in the marine environment.
In addition to their enrichment from either riverine or marine waters, GHGs may accumulate in sea ice from bubbles. As sea ice temperature decreases, the brine volume decreases with the potential to affect pressure, bubbles, and the solubility of gases trapped within sea ice (Crabeck et al., 2019). Gas supersaturation leads to bubble nucleation (Zhou et al., 2013). Because gas bubbles move upward due to their buoyancy, while dissolved compounds are subject to gravity drainage (i.e., downward movement), gas bubble formation can drive gas accumulation in sea ice (Zhou et al., 2014). The enrichment of CH4 in sea ice could be due to CH4 bubbles released from the sediment that migrate through the water column due to their own buoyancy and accumulate within the ice structure as sea ice grows.
Vertical profiles of CH4 concentration in sea ice exhibit maximum concentrations in the upper layer (Figure 3), except at Vaasa where an increase in CH4 concentration was observed at the surface but the maximum value was at the sea ice bottom. The increase of gas concentration at the surface is less marked for N2O and DIC than for CH4. Excess of CH4 at the ice surface could be associated with the formation of meteoric ice due to the flooding of CH4-supersaturated river water onto the ice surface. In addition, an enhanced ebullition flux of CH4 early in the ice formation could contribute to the higher CH4 concentrations. In fall, organic matter from the summer phytoplankton bloom and relatively high bottom water temperatures may support significant benthic mineralization (Humborg et al., 2019). Decrease of CH4 ebullition fluxes from fall to winter has been documented for temperate continental aquatic systems (Wilkinson et al., 2015; Tušer et al., 2017) and could explain the decrease in CH4 concentrations with sea ice depth, associated with the slower growth rate of the ice as it thickens, increasing the efficiency of impurity rejection (salts, gases) to the water column as sea ice grows (Granskog et al., 2006b).
Once GHGs are trapped within the ice structure, gas bubbles will not escape to the atmosphere as long as the sea ice permeability is low (e.g., brine volume < 5%; Golden et al., 2007), and GHGs could remain within the ice matrix until the onset of sea ice melt. CH4 temporarily trapped within the ice structure or under-ice water could be subject to microbial oxidation to CO2 within sea ice and meltwater and buffered by the carbonate system of the ocean. However, the prevalence and activity of the potential CH4 oxidizers need to be verified in future studies. Nevertheless, our results suggest that sea ice in the BB could act as a sink for CH4, whereas the increasing temperatures and loss of sea ice cover due to global warming may increase CH4 fluxes in the Baltic Sea (Humborg et al., 2019).
While the role of sea ice in GHG exchanges across the ocean–sea ice–atmosphere interface is poorly constrained, observations over brackish sea ice, as in the BB, are still lacking. Due to ongoing changes in the polar regions, with undefined consequences for sea ice biogeochemical cycles (Lannuzel et al., 2020), investigations of low salinity sea ice in the Baltic Sea may help to constrain future changes in the warming Arctic, especially on shallow shelves where increased freshwater input can significantly reduce sea surface salinity.
The carbonate system of both sea ice and the water column in this study appeared to be driven mainly by salinity. High CH4 and N2O concentrations were observed in both sea ice and the water column, but the ice cover appeared to be enriched in GHGs compared to the water column. This enrichment is due to the concentration effect, associated with sea ice growth, of already enriched parent seawater and/or enriched riverine freshwater. As sea ice temperature and brine volume decreases, the solubility of GHGs trapped within the ice matrix decreases, leading to bubble formation. In addition, gas bubbles originating from the sediment may be incorporated within the ice cover and contribute to the sea ice enrichment in GHGs. Due to the low brine volume of brackish sea ice, bubbles accumulated in the ice matrix are expected to remain until the onset of sea ice melt. This temporary GHG storage in sea ice may have implications for the magnitude and timing of the BB GHG fluxes, which would be further constrained if sea ice supports active microbial CH4 oxidation.
Data accessibility statement
Data are available at HydroShare: https://doi.org/10.4211/hs.13f2735419cf4d2cac8ab353725e3145.
This research was primarily funded by the International Union for Conservation of Nature, with additional financial support from the Canadian Natural Sciences and Engineering Research Council and the Canada Foundation for Innovation. BD is a research associate of the F.R.S-FNR. The authors would like to thank Jody Deming, Mats Granskog, and an anonymous reviewer for their constructive comments that greatly improved the overall quality of the article.
Financial support for this study was provided by the Canada Excellence Research Chairs program, the Canada Research Chairs program, the Natural Sciences and Engineering Research Council, and the Canada Foundation for Innovation. E-ER is funded by the Academy of Finland (PRICE 325140). The study utilized SYKE marine research infrastructure as a part of the national FINMARI RI consortium.
The authors declare no competing interests.
Field sampling was performed by N-XG, KMM, HK, E-ER. Laboratory and data analyses were performed by N-XG, ML, BD. All authors contributed to the article redaction and have reviewed and approved the submitted version for publication.
How to cite this article: Geilfus, N-X, Munson, KM, Eronen-Rasimus, E, Kaartokallio, H, Lemes, M, Wang, F, Rysgaard, S, Delille, B. 2021. Landfast sea ice in the Bothnian Bay (Baltic Sea) as a temporary storage compartment for greenhouse gases. Elementa: Science of the Anthropocene 9(1). DOI: https://doi.org/10.1525/elementa.2021.00028
Domain Editor-in-Chief: Jody W. Deming, University of Washington, Seattle, WA, USA
Associate Editor: Stephen F. Ackley, Department of Geological Sciences, University of Texas at San Antonio, TX, USA
Knowledge Domain: Ocean Science
Part of an Elementa Special Feature: Insights into Biogeochemical Exchange Processes at Sea Ice Interfaces (BEPSII-2)