Estimating sea–air CO2 fluxes in coastal seas remains a source of uncertainty in global carbon budgets because processes like primary production, upwelling, water mixing, and freshwater inputs produce high spatial and temporal variability of CO2 partial pressure (pCO2). As a result, improving our pCO2 baseline observations in these regions is important, especially in sub-Arctic and Arctic seas that are experiencing strong impacts of climate change. Here, we show the patterns and main processes controlling seawater pCO2 and sea–air CO2 fluxes in Hudson Bay during the 2018 spring and early summer seasons. We observed spatially limited pCO2 supersaturation (relative to the atmosphere) near river mouths and beneath sea ice and widespread undersaturated pCO2 in offshore and ice-melt-influenced waters. pCO2 was highly correlated with salinity and temperature, with a limited but statistically significant relationship with chlorophyll a and fluorescent dissolved organic matter. Hudson Bay on average was undersaturated with respect to atmospheric CO2, which we attribute mainly to the dominance of sea-ice meltwater. We calculated an average net CO2 flux of about –5mmol CO2 m–2 day–1 (–3.3 Tg C) during the spring and early summer seasons (92 days). Combining this result with extrapolated estimates for late summer and fall seasons, we estimate the annual CO2 flux of Hudson Bay during the open water season (184 days) to be –7.2 Tg C. Our findings indicate that the bay on average is a weaker CO2 sink than most other Arctic seas, emphasizing the importance of properly accounting for seasonal variability in the Arctic coastal shelves to obtain reliable sea–air CO2 exchange budgets.

Global oceans slow the effects of climate change by absorbing anthropogenic carbon dioxide (CO2) emitted by the burning of fossil fuels, land-use change, and cement production. However, the invasion of this anthropogenic CO2 into the ocean has caused a decline in surface ocean pH leading to ocean acidification (Orr et al., 2005; Gattuso and Hansson, 2011). Friedlingstein et al. (2019) quantified the major components of the global carbon budget (atmosphere, land, and ocean reservoirs) over the previous decade (2009–2018) and estimated that the global ocean absorbed about 23% of the total anthropogenic CO2 emissions. However, this global budget neglected the carbon cycle in coastal ocean waters largely because the spatial and temporal heterogeneity of biogeochemical processes in these waters make it difficult to accurately scale up results from observational studies. Despite great efforts over the last decade to assess the coastal contribution to the global carbon budget (e.g., Chen and Borges, 2009; Bauer et al., 2013; Regnier et al., 2013; Laruelle et al., 2014; Bourgeois et al., 2016; Roobaert et al., 2019), the CO2 budgeting of these waters is still subject to uncertainty and lacks consensus. The coastal ocean has been impacted significantly by human activities (hydroelectric damming, eutrophication, overfishing, etc.) and is now considered one of the most sensitive parts of the marine environment to climate change (e.g., Arctic Climate Impact Assessment, 2005; Liu et al., 2010; Liao et al., 2015; Breitburg et al., 2018). Furthermore, the coastal shelves play a disproportionately significant role in the global carbon cycle in comparison with its surface area (about 7%–10% of the oceanic surface area) by linking the oceanic, terrestrial, and atmospheric carbon reservoirs (Wollast, 1998; Liu et al., 2010). Hence, studying the carbon cycle in these waters is essential to understand future shifts in response to climate change.

The coastal ocean consists of diverse but closely connected ecosystems that include estuaries, rivers, tidal wetlands, and the continental shelf. The continental (or coastal) shelves are shallow and have been the focus of many studies due to their key role in the global biogeochemical cycling of carbon and nutrients (e.g., Chen et al., 2013; Bourgeois et al., 2016) as well as their importance to local livelihoods. For instance, the ocean’s biological pump operates more efficiently in coastal shelves relative to the open ocean as a result of high terrestrial inputs and efficient use of nutrients (Gattuso et al., 1998; Chavez et al., 1999). The high biological production in these waters usually leads to atmospheric CO2 drawdown and subsequent CO2 export to the subsurface waters via the biological pump (Ducklow et al., 2001). The eventual outflow of CO2-enriched deep water from coastal shelf regions to the subsurface waters of the adjacent deep oceans constitutes “the continental shelf pump,” which is thought to contribute significantly to the global ocean’s atmospheric CO2 uptake (Tsunogai et al., 1999). Furthermore, the production, export, and degradation of organic matter in these waters are several times higher than in the open ocean (Wollast, 1998; Borges et al., 2005). The coastal shelves contribute 10%–30% of global marine primary production, 80% of organic carbon, and 30%–50% of inorganic carbon burial in sediments, and about 50% of the organic carbon supplied to the deep open ocean (Wollast, 1998; Liu et al., 2010). As a result, the sea–air CO2 exchanges are expected to be more intense in these waters than in the open ocean.

Coastal shelves show a latitudinal pattern in CO2 fluxes, with the mid-to-high-latitude (30°–90°) shelves acting as oceanic sinks of CO2, whereas the low latitude shelves (equator to 30°) act as sources of CO2 to the atmosphere or are nearly neutral (Borges et al., 2005; Cai et al., 2006). This latitudinal pattern can be attributed to a complex interaction among temperature (which alters the solubility of CO2), the amount of biological activity, the magnitude of upwelling and mixing from deeper waters, and the export of terrestrial organic and inorganic carbon from rivers (Takahashi et al., 2009; Dai et al., 2012). Adding further complexity, the processes controlling CO2 dynamics vary with proximal distance from shore (Wollast, 1998; Bauer et al., 2013). Inner shelf waters close to land tend to be sources of CO2 due to respiration of terrestrial organic carbon and lateral transport of CO2-enriched waters from the adjacent inshore system, whereas the mid-to-outer shelf waters generally act as a sink of CO2 due to increased primary production as a result of nutrient availability from upwelling and mixing across the shelf break (Walsh, 1991; Jiang et al., 2008). Moreover, coastal shelves show a temporal variability in sea–air CO2 fluxes as a result of variability in riverine input, phytoplankton blooms, ice formation and melting, and sea surface temperature (e.g., Takahashi et al., 2009; Rysgaard et al., 2012; Ahmed et al., 2019). For example, Borges et al. (2005) documented seasonal changes of surface seawater CO2 partial pressure (pCO2) as high as 400 μatm in the coastal shelves as a result of riverine input and other coastal processes.

In high-latitude coastal seas, changes in sea-ice coverage (Notz and Stroeve, 2016), water temperature (Shepherd, 2016), riverine input (Carmack et al., 2016), glacial melt input (Pierre et al., 2019), and primary production (Arrigo and van Dijken, 2015) will have a significant impact on the distribution and magnitude of pCO2 and thus the strength of atmospheric CO2 uptake. Given these potential changes, as well as the spatial and temporal variability in pCO2, studying the spatiotemporal patterns of pCO2 independently in the Arctic and sub-Arctic coastal shelves is important to improve our knowledge of the main drivers affecting pCO2 and to increase the reliability of sea–air CO2 flux budgets.

The Hudson Bay system (including James Bay, Foxe Strait, Ungava Bay, and Hudson Strait) is a sub-Arctic/Arctic continental shelf sea that receives nearly one-third of Canada’s river discharge from more than 42 rivers (Déry et al., 2011). Hudson Bay and its watershed have been recognized as highly sensitive to climate variability relative to other Arctic shelves (Arctic Climate Impact Assessment, 2005; Markus et al., 2009). For instance, several studies have documented changes in timing and volume of river discharge (Déry et al., 2011, 2018), the timing of sea-ice formation and melting (Saucier et al., 2004; Hochheim et al., 2011), and biological production (Sakshaug, 2004; Sibert et al., 2011). Permafrost melting due to climate change may also profoundly affect the quality, quantity, and composition of organic matter in Hudson Bay (Gough and Leung, 2002; Mallory et al., 2010; Godin et al., 2017). These potential changes could have a drastic impact on sea–air CO2 fluxes in Hudson Bay.

Although included in several global coastal ocean CO2 budgets (e.g., Chen and Borges, 2009; Chen et al., 2013), Hudson Bay has not been included in other global and Arctic-specific budgets (Cai et al., 2006; Bates and Mathis, 2009; Laruelle et al., 2014). This limited inclusion is due to a lack of observations on spatial and temporal scales. Only two studies investigating sea–air CO2 fluxes in Hudson Bay have been published (Else et al., 2008a, 2008b). Using pCO2 measurements at discrete stations collected in September and October, Else et al. (2008a) found that the bay acts as a minor sink in October, while Else et al. (2008b) applied a remote sensing approach to determine that Hudson Bay acts as a weak CO2 source during August and September. The limited temporal and spatial extent of those observations means that we have no reliable estimate of sea–air CO2 exchange in Hudson Bay over the complete open water period. Therefore, our baseline estimates in this large inland sea need to be extended, especially into the spring and early summer seasons, when the ecosystem is likely subject to rapid and pronounced physical and biological changes (e.g., Else et al., 2012).

Here, we report the first spring and early summer measurements of pCO2 and estimates of sea–air CO2 fluxes in Hudson Bay and Hudson Strait. The main objectives of this study were to (1) examine the spatial and temporal patterns of pCO2; (2) identify the key processes driving the variability in the observed pCO2; and (3) estimate the sea–air CO2 fluxes during the 2018 spring and early summer seasons in Hudson Bay and Hudson Strait.

Hudson Bay is a large shallow inland sea located in central Canada, with an average depth of 125 m and a maximum depth of 250 m, while adjacent areas of Foxe and Hudson Straits can reach depths of 400 m (Prinsenberg, 1986; Saucier et al., 2004). Hudson Bay transitions from complete ice cover in the winter to open water in summer, which adds about 140–160 cm of freshwater when sea-ice melts from May to July (Prinsenberg, 1988). Freshwater from river runoff contributes about 85 cm, distributed around the perimeter of the bay (Prinsenberg, 1986), and largely constrained near the coast (Granskog et al., 2007). Peak river runoff occurs in May, whereas the minimum occurs in March, with a notable change in recent decades toward higher winter discharge for hydroelectric power generation (Déry et al., 2011). This massive seasonal freshwater input to the bay leads to a strong stratification in the upper surface layer (Granskog et al., 2011; Burt et al., 2016; Ahmed et al., 2020), which impacts the mixing of heat and nutrients between the surface and deep waters (Drinkwater and Jones, 1987; Prinsenberg, 1988). Water flowing out of the Canadian Arctic Archipelago enters Hudson Bay at its northwestern end through Foxe Strait, while the northern portion of Hudson Strait allows water inflow of both Atlantic and Arctic origin from the Labrador Sea. Water exits Hudson Bay via the southern portion of Hudson Strait (Ingram and Prinsenberg, 1998). The general circulation in the bay is cyclonic (counterclockwise), which can be attributed to wind-driven transport and estuarine circulation buoyancy resulting from river runoff and ice melt/formation (Prinsenberg, 1986; Saucier et al., 2004). In this study, we focus on Hudson Bay (approximately 0.8 × 106 km2) and Hudson Strait (approximately 0.2 × 106 km2), where we conducted the majority of our field measurements.

We collected our field measurements during a 6-week survey from May 25 to July 13, 2018, onboard the research icebreaker CCGS Amundsen (Figure 1). This study was part of the BaySys project (see other papers in this special feature of Elementa), a comprehensive study examining the relative influences of climate change and hydroelectric regulation on Hudson Bay’s physical, chemical, and biological systems. The field program was the first-ever research cruise in Hudson Bay in the spring and early summer, at a time when sea-ice cover was still partially intact, and sea-ice melt and river runoff are near their peaks.

Figure 1.

Map of the study area. Ship track visualized with the Julian day (i.e., May 25 to July 13) during the 2018 cruise in Hudson Bay. The hatched area shows the region with sea-ice concentration greater than 90% on May 28, 2018, based on the ice charts provided by the Canadian Ice service. pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f1

Figure 1.

Map of the study area. Ship track visualized with the Julian day (i.e., May 25 to July 13) during the 2018 cruise in Hudson Bay. The hatched area shows the region with sea-ice concentration greater than 90% on May 28, 2018, based on the ice charts provided by the Canadian Ice service. pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f1

### 3.1. Surface water measurements

Surface water pCO2 was measured during the cruise using an underway pCO2 system (General Oceanics model 8050; Pierrot et al., 2009). This automated pCO2 system continuously sampled water through a shower-type equilibrator system at a nominal flow rate of 2.5 L min–1 from a through-hull inlet located approximately at 7 m beneath the water surface. The precision of the underway pCO2 system is 0.1 μatm, with an overall accuracy estimated at 2 μatm (Pierrot et al., 2009). The system was calibrated every 8 h using four certified standard gases (ultra-high purity N2 as a zero gas and three CO2/air mixtures at concentrations of 295.8, 448.6, and 551.7 ppm) traceable to World Meteorological Organization standards. The system was very stable during our sampling time, with drift typically less than 0.2 μatm between calibrations, contributing to uncertainty in our measurements less than 0.5%. We calculated the seawater pCO2 using the mole fraction of CO2 (xCO2) based on Dickson et al. (2007) and following the recommendation of Pierrot et al. (2009) for pCO2 calculations and xCO2 corrections using the measured standards.

Using discrete water samples, Ahmed et al. (2020) found that the sea-ice melt created a shallow, stratified freshwater layer that played a major role in diluting the surface pCO2 across Hudson Bay. To correct for vertical pCO2 gradients resulting from freshwater stratification during our study, we applied a linear regression that predicts the pCO2 gradient as a function of observed pCO2 (fully described in Ahmed et al., 2020) to about 45% of the total underway measurements (those measured 1–5 weeks after ice breakup with wind speed less than 6 m s–1). For pCO2 in areas measured less than 1 week after ice breakup, stratification was considered to be too strong for confident corrections, and those records (about 10% of the total measurements) were discarded. For pCO2 in areas measured more than 5 weeks after ice breakup or areas 1–5 weeks after ice breakup with wind speeds greater than 6 m s–1, no correction was applied (45% of the measurements) as the near-surface vertical pCO2 and salinity gradients were not likely significant under those conditions (Ahmed et al., 2020).

The underway pCO2 system has flow-through temperature, salinity, and oxygen sensors (Idronaut model Ocean Seven 315), which were factory-calibrated before deployment. However, we opted to use the temperature and salinity recorded by a separate thermosalinograph (TSG) that draws water from the same intake line but is installed closer to the ship inlet compared to the underway pCO2 system (distance difference between two instruments in the engine room was about 3 m). A chlorophyll a (Chl a) fluorescence sensor (WestStar fluorometer with an accuracy of 0.03 µg L–1) and a fluorescing dissolved organic matter (FDOM) sensor (WETLabs ECOFLD; 370 nm excitation, 460 nm emission; Guéguen et al., 2015) were also installed on the TSG line. Both the FDOM and Chl a sensors were also factory-calibrated before deployment. Measurements of surface nitrate (NO3) concentrations and surface photosynthetic active radiation (SPAR) across Hudson Bay were collected using the Seabird Satlantic MBARI-ISUS V3 nitrate and the Biospherical Instruments SPAR (QCR 2200) sensors, respectively, which were mounted on the ship’s rosette. Additional details about sensors, calibrations, and data processing of the TSG and rosette data are reported in Amundsen Science Data Collection (2019).

Despite the proximity of the pCO2 system to the inlet (<5 m), we observed a slight increase in the recorded water temperature in the equilibrator relative to the on-station measurements made by the ship’s conductivity-temperature-depth (CTD) sensor on the rosette. Therefore, we corrected the pCO2 measurements for this temperature increase following the procedure described by Takahashi et al. (1993). Using the maximum uncertainty in the in situ seawater temperature measurements through the pCO2 temperature correction, we estimate that the total uncertainty in the final pCO2 measurements is less than 2%. We processed the underway data to remove any measurements collected when the system was undergoing cleaning, calibration, or restricted water flows from the intake line (i.e., water flow <1.5 L m–1) during icebreaking operations. Further details about the data processing of the underway pCO2 system are described in Ahmed et al. (2019).

### 3.2. Atmospheric measurements

Atmospheric CO2 mixing ratio was measured every minute using a closed-path infrared gas analyzer (Li-COR model LI-7000). The analyzer was installed in a container on the foredeck and connected to a gas sampling system that drew air from an inlet located on the ship’s foredeck meteorological tower at a height of 17 m above the sea surface. The LI-7000 was calibrated twice a day with two certified gas standards traceable to World Meteorological Organization standards, and the mixing ratio measurements were converted into pCO2atm following Dickson et al. (2007). Wind speed and air temperature were measured using a conventional propeller anemometer (RM Young Co. model 15106MA) and HMP45C212 temperature sensor, respectively. These sensors were mounted on the meteorological tower at a height of approximately 16 m above sea level. Wind speed was corrected to a height of 10 m assuming a log-linear wind profile and a neutral surface layer (Stull, 1988). We discarded the data when instruments on the tower were cleaned and/or calibrated or when the wind was blowing from the stern of the vessel.

For logistical reasons, the foredeck tower and associated sensors were removed 4 weeks into the cruise (on July 2, 2018). As a result, we opted to utilize wind velocities at 10 m height and surface air temperature provided by the National Centers for Environmental Prediction North American Regional Reanalysis (NARR) to have complete measurements along the ship track. The NARR wind and air temperature data are available eight times daily (i.e., every 3 h) with a horizontal grid spacing of about 32 km × 32 km (Mesinger et al., 2006). We clipped and resampled the NARR data using a bilinear algorithm to 1-km spatial resolution and then matched it with the ship track based on spatial location and closest estimate to our sampling time. The NARR and in situ measurements show good agreement over the period when the in situ sensors were deployed, especially for the wind speed data (Table S1). The atmospheric pCO2 concentrations from May 25 to July 2, 2018, ranged from 400 to 419 μatm, with a mean of 408 ± 2.8 (1σ) μatm (n = 17,308).

### 3.3. Sea-ice data

Daily sea-ice observations along the ship track were obtained from the Advanced Microwave Scanning Radiometer 2 (AMSR2) passive microwave ice product (approximately 3-km spatial resolution) provided by the University of Bremen. The AMSR2 characteristics and the ARTIST Sea Ice algorithm used for calculating the daily ice concentrations are described in detail by Spreen et al. (2008). We used the spatial locations of our ship track to index the daily ice maps during our study time and then retrieve the sea-ice concentration from the closest temporally located observation. We compared the AMSR2 ice data with weekly ice charts from the Canadian Ice Service Digital Archive (CISDA) to ensure that the AMSR2 product did not underestimate ice concentrations along the ship track due to melt ponds on ice and/or weather conditions. The CISDA data are based on integration and analysis of data from satellite imagery (i.e., synthetic aperture radar and optical), weather and oceanographic information, and visual observations from ship and aircraft. Past studies have demonstrated the reliability of CIS data to monitor sea-ice cover change in the Canadian Arctic (Tivy et al., 2011) and its application to sea–air CO2 exchange studies (Ahmed and Else, 2019). Figure S1 shows the sea-ice observations obtained from AMSR2, CISDA, and ship observations recorded by the coast guard officers on the ship. Although the AMSR2 product showed relatively better agreement with the ship observations compared to the CISDA (Figure S1), we noticed that the AMSR2 underestimated ice concentrations in a few locations in Hudson Bay, especially along the southwest and east coasts of Hudson Bay (June 27–30, 2018), where we observed more melt ponds on the ice. For these reasons and to make sure we accounted for any biases in sea-ice concentrations, we decided to use both AMSR2 and CISDA products as upper and lower estimates, respectively, when calculating the sea–air CO2 fluxes in Hudson Bay.

### 4.1. Underway observations

Figure 2 shows the spatial variations of sea-ice concentration, salinity, temperature, and pCO2 throughout the study area. At the start of our cruise, we observed a large polynya on the northwest side of Hudson Bay, with higher ice concentrations (>50%) encountered mainly in the south of Hudson Strait and the central and eastern side of the bay (Figures 2a and S2). The northwest polynya (hereafter NW polynya) in Hudson Bay has been observed occasionally in winter and spring (e.g., Gough et al., 2004; Landy et al., 2017), resulting from strong offshore westerly winds that advect the existing ice cover eastward and enhance local ice production (Saucier et al., 2004). However, sea-ice melt was actively underway throughout our study period, especially close to the ice edge and coastal areas. The rapid change in the sea-ice cover is evident from maps of sea-ice concentrations throughout the bay at the start, middle, and end of our cruise (Figure S2). At the end of our cruise, sea ice was mainly limited to the eastern side of the bay (Figure S2).

Figure 2.

Spatial patterns of sea-ice concentration and underway measurements of salinity, temperature, and pCO2. Spatial variability of (a) sea-ice concentration using the AMSR2 ice data and underway measurements of (b) sea surface salinity, (c) sea surface temperature, and (d) surface seawater pCO2 along the ship track in Hudson Bay from May 25 to July 13, 2018. pCO2 = CO2 partial pressure; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f2

Figure 2.

Spatial patterns of sea-ice concentration and underway measurements of salinity, temperature, and pCO2. Spatial variability of (a) sea-ice concentration using the AMSR2 ice data and underway measurements of (b) sea surface salinity, (c) sea surface temperature, and (d) surface seawater pCO2 along the ship track in Hudson Bay from May 25 to July 13, 2018. pCO2 = CO2 partial pressure; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f2

The spatial patterns of salinity and temperature were associated with the spatial patterns of sea-ice concentrations (Figure 2). Salinity ranged from minima of <20 close to the sea-ice edges and the river-dominated coastal region to a high value of 33 in offshore ice-free regions (Figure 2b and Table 1). We observed higher surface seawater temperatures (>8 °C) close to the Churchill and Nelson rivers, and a low value of surface temperature (about –1.7 °C) close to the sea-ice edge (Figure 2c and Table 1). Typically, colder (<1 °C) and less saline waters (salinity between 10 and 25) were encountered in the southern and eastern side of the bay near the ice edge, while relatively warmer (>4 °C) more saline (salinity from 15 to 29) waters were observed in the river-dominated coastal zone (Figure 2a and c). We also observed high salinity (>30) associated with cold water (<0 °C) in Hudson Strait and the southwest side of Baffin Island. The observed variability in salinity matches the general cyclonic (counterclockwise) circulation in Hudson Bay (Prinsenberg, 1986; Ingram and Prinsenberg, 1998), as the marine water entering Hudson Bay from Hudson Strait and Foxe Basin is usually diluted by sea-ice melt and river runoff before exiting the bay to the Labrador Sea via Hudson Strait.

Table 1.

Summary statistics of the main measured and calculated surface seawater variables in Hudson Bay.a DOI: https://doi.org/10.1525/elementa.2020.00130.t1

StatisticsMeasured VariablesCalculated Variables
pCO2TemperatureSalinityChl aFDOMO2 SaturationWind SpeedpCO2nonthermalCO2 Fluxb
(μatm)(°C)(μg L–1)(Count)(%)(m s–1)(μatm)(mmol m–2 day–1)
Mean 316.8 0.8 29.9 1.1 685.5 111.1 4.4 319.1 –5.1
Standard deviation 60.9 1.7 2.8 1.2 404.6 9.1 2.9 66.5 9.3
Minimum 125.3 –1.8 11.9 0.0 51.0 90.0 0.01 124.5 –83.9
Maximum 649.7 11.1 33.1 10.4 9,039.0 199.4 15.8 548.8 8.9
Median 338.4 0.2 30.7 0.7 694.0 109.1 3.3 333.8 –1.6
StatisticsMeasured VariablesCalculated Variables
pCO2TemperatureSalinityChl aFDOMO2 SaturationWind SpeedpCO2nonthermalCO2 Fluxb
(μatm)(°C)(μg L–1)(Count)(%)(m s–1)(μatm)(mmol m–2 day–1)
Mean 316.8 0.8 29.9 1.1 685.5 111.1 4.4 319.1 –5.1
Standard deviation 60.9 1.7 2.8 1.2 404.6 9.1 2.9 66.5 9.3
Minimum 125.3 –1.8 11.9 0.0 51.0 90.0 0.01 124.5 –83.9
Maximum 649.7 11.1 33.1 10.4 9,039.0 199.4 15.8 548.8 8.9
Median 338.4 0.2 30.7 0.7 694.0 109.1 3.3 333.8 –1.6

pCO2 = CO2 partial pressure; Chl a = chlorophyll a; FDOM = fluorescing dissolved organic matter; AMSR2 = Advanced Microwave Scanning Radiometer 2; NARR = National Centers for Environmental Prediction North American Regional Reanalysis.

aThe values presented are area-weighted based on a 1 × 1 km grid across the study area (n = 7,817 for all the variables, except for wind speed where n = 5,580).

bCalculated sea–air CO2 fluxes using the AMSR2 sea-ice data and NARR wind data in Hudson Bay.

We calculated the area-weighted average and standard deviation of these observations based on a 1 km × 1 km grid to minimize the effect of sampling biases on our calculation of the bay-wide flux resulting from long occupations of some station locations. This binning resulting in decreasing the number of ship observations on the waterside from about 19,500 to 7,817 and on the atmospheric side from about 17,500 to 5,580. Hereafter, all the statistics reported in this study are area-weighted and based on n = 7,817 for the sea and n = 5,580 for the atmosphere. A summary of the range, average, and variability of the observed underway measurements is available in Table 1. The pCO2 ranged between 125 and 650 μatm (Figure 2d, Table 1), with the lowest measurements observed close to the ice edge on the east side of the bay, the NW side of the bay, and in Hudson Strait (125–280 μatm), while the highest pCO2 values were mainly observed in ice-covered waters (380–550 μatm) and close to the Churchill and Nelson Estuaries (390–465 μatm). We calculated the average of pCO2 measurements across Hudson Bay to be 317 ± 61 μatm, which is significantly lower (p = 0) than the average of atmospheric pCO2 (408 ± 3 μatm), indicating that overall Hudson Bay was undersaturated in pCO2 relative to the atmosphere.

The measured pCO2 showed no clear evidence of temporal changes during the study period, although it was clearly correlated in various locations with variations in salinity and temperature (Figure 3a, b). This result suggests that the observed variability was linked to the ship’s spatial proximity to sea-ice or riverine input rather than the timing of sampling. However, we observed a relatively narrow range of salinity and temperature values (cooler, more saline water) during the first 2 weeks (May 30 to June 15), with higher variability and frequency of warmer, less saline water thereafter. This pattern is consistent with the heavy sampling in the NW polynya at the beginning of the cruise and increasing freshwater abundance as the ship entered the southern coastal domain where most of the runoff enters the bay. We also observed an increase in pCO2 along the SW coast, associated with a decrease in salinity, an increase in temperature, and low sea-ice concentrations. On the east coast, we observed a decrease in pCO2 values associated with low temperature and salinity. Low pCO2 values were also observed in Hudson Strait associated with relatively higher temperatures and salinity. Additionally, the pCO2 values showed some correlation with variability in sea-ice concentrations (Figure 3c). For instance, we observed that the high spikes of pCO2 were usually associated with a higher sea-ice concentration (>80%), whereas the lowest pCO2 observations were mainly associated with lower sea-ice concentration (<40%; Figure 3c).

Figure 3.

Temporal patterns of underway surface seawater measurements of pCO2 visualized with salinity, temperature, and sea-ice concentration. Temporal variability of underway surface seawater pCO2 visualized with sea surface (a) salinity, (b) temperature, and (c) AMSR2 sea-ice concentration during our study in Hudson Bay. The black and red dashed lines represent the area-weighted average of pCO2 in seawater and atmosphere, respectively. pCO2 = CO2 partial pressure; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f3

Figure 3.

Temporal patterns of underway surface seawater measurements of pCO2 visualized with salinity, temperature, and sea-ice concentration. Temporal variability of underway surface seawater pCO2 visualized with sea surface (a) salinity, (b) temperature, and (c) AMSR2 sea-ice concentration during our study in Hudson Bay. The black and red dashed lines represent the area-weighted average of pCO2 in seawater and atmosphere, respectively. pCO2 = CO2 partial pressure; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f3

### 4.2. Processes affecting surface pCO2 variability

#### 4.2.1. Seawater temperature

The influence of seawater temperature on pCO2 is governed by the thermodynamic equilibrium of the inorganic carbon system (Takahashi et al., 1993). For example, if pCO2 is only controlled by temperature change with no other physical or biogeochemical processes playing a role, the pCO2 in water will double for every 16 °C increase in temperature (Takahashi et al., 1993). The weighted-area average of seawater temperature observed during our sampling period was 0.8 °C ± 1.7 °C, with values ranging from –1.8 °C to 11.1 °C (Table 1). In order to assess the relative importance of the competing effects of thermal (i.e., temperature) and nonthermal (i.e., biological activity and/or physical water mixing) effects on changes in pCO2 values, we calculated the temperature-normalized pCO2 (pCO2nonthermal) and the effects of temperature on the pCO2 (pCO2thermal) using the decomposition used by Takahashi et al. (2002). More details about the thermal and nonthermal pCO2 calculations are described in Text S1.

In Figure 4, we show the temporal variability of pCO2nonthermal and pCO2thermal. The high temporal variability of pCO2nonthermal compared to pCO2thermal in the first 3 weeks of our sampling reflects the relatively stable temperatures and significant contribution of nonthermal processes (i.e., biological activity, water mixing, and sea–air CO2 exchange) to the observed variations in pCO2. For example, the high pCO2nonthermal (>450 μatm) in the NW region could be attributed to the inflow of pCO2-rich water from rivers, mixing with subsurface waters, and/or respiratory CO2 production.

Figure 4.

Temporal variability of calculated pCO2nonthermal and pCO2thermal. Temporal variability of calculated pCO2nonthermal (orange), pCO2thermal, (blue), and ΔpCO2ratio (pCO2thermal / pCO2nonthermal, gray) in Hudson Bay. The horizontal dashed line indicates a ΔpCO2ratio ratio of 1, with the arrows indicating ratios greater than 1, due to the dominance of thermal processes (thermal), and less than 1, due to the dominance of biological, mixing, and sea–air exchange processes (nonthermal). pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f4

Figure 4.

Temporal variability of calculated pCO2nonthermal and pCO2thermal. Temporal variability of calculated pCO2nonthermal (orange), pCO2thermal, (blue), and ΔpCO2ratio (pCO2thermal / pCO2nonthermal, gray) in Hudson Bay. The horizontal dashed line indicates a ΔpCO2ratio ratio of 1, with the arrows indicating ratios greater than 1, due to the dominance of thermal processes (thermal), and less than 1, due to the dominance of biological, mixing, and sea–air exchange processes (nonthermal). pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f4

Starting at about the 4th week of our study period (June 27), during the shipping transit along the SW coast of the bay, pCO2thermal became more variable (Figure 4). This increased variability is not surprising as we measured higher water temperatures during this time (Figure 3b), which were mostly linked to an increase in the air temperature with time (Figure S3a) and/or the presence of high-temperature riverine waters from Nelson and Churchill Rivers (Figure 2c). The ratio of thermal to nonthermal effects on pCO2pCO2ratio) shown in Figure 4 illustrates that the surface water pCO2 in Hudson Bay and Hudson Strait was controlled mainly by nonthermal processes at the start of the sampling period, with the thermal effects more dominant at the end of our study.

As pCO2nonthermal varies substantially through the entire study, we show the spatial variability of pCO2nonthermal in Figure 5a. The average of pCO2nonthermal is 319 ± 66 μatm, with a somewhat smaller range (124–548 μatm) than the observed pCO2 values (Table 1). The largest differences between pCO2nonthermal and observed pCO2 were close to river mouths. For example, the differences between pCO2nonthermal values and the observed pCO2 measurements in the Nelson and Churchill Estuaries on average were about 55 μatm (Figures 5a and 2d), indicating that mixing with relatively warm water introduced by Nelson and Churchill Rivers was contributing to the observed high pCO2.

Figure 5.

Spatial patterns of calculated pCO2nonthermal and underway measurements of biological variables. Spatial variability of (a) temperature-normalized surface water pCO2 (pCO2nonthermal), (b) FDOM, (c) oxygen saturation, and (d) Chl a concentrations in Hudson Bay. pCO2 = CO2 partial pressure; FDOM = fluorescing dissolved organic matter; Chl a = chlorophyll a. DOI: https://doi.org/10.1525/elementa.2020.00130.f5

Figure 5.

Spatial patterns of calculated pCO2nonthermal and underway measurements of biological variables. Spatial variability of (a) temperature-normalized surface water pCO2 (pCO2nonthermal), (b) FDOM, (c) oxygen saturation, and (d) Chl a concentrations in Hudson Bay. pCO2 = CO2 partial pressure; FDOM = fluorescing dissolved organic matter; Chl a = chlorophyll a. DOI: https://doi.org/10.1525/elementa.2020.00130.f5

#### 4.2.2. Biological production and FDOM

The distributions of oxygen saturation, Chl a, and FDOM (Figure 5b-d) indicated that for the most part, biological processes and terrestrial organic matter were not a major driver of pCO2 distribution in Hudson Bay at the time of sampling. This interpretation is supported by the low Pearson’s correlation coefficients (r < .2; p < .05) among pCO2, Chl a, and FDOM across the study area (see Section 4.3). However, there were some notable, localized exceptions. For example, the FDOM signal was relatively high close to the coast, especially in the SW region of Hudson Bay (Figures 5b and S5), associated with high terrestrial input in this region (Guéguen et al., 2011), while low FDOM signals (<700 counts) were observed in the offshore ice-free water and the NW polynya, likely due to photobleaching in the clearer waters (Dainard et al., 2015). The highest FDOM values (up to 9,000 counts) in our study area were measured east of the Nelson Estuary close to sediment-laden sea ice. Similar to our observations, Granskog et al. (2007) and Guéguen et al. (2011) found high chromophoric dissolved organic matter and FDOM values near the coast in Hudson Bay due to terrestrial input by rivers. Overall, we did not observe any consistent relationship between FDOM and pCO2 in our study area, indicating that FDOM played a minor role in pCO2 variability within Hudson Bay during the study period.

The average of Chl a during our sampling time was 1.1 ± 1.2 µg L–1, with minimum and maximum observations of 0 µg L–1 and approximately 10 µg L–1, respectively (Table 1). O2 saturation averaged 111 ± 9% with a range varying from 90% to 199% across the study area (Table 1). Typically, the response of primary production to ice melt during the spring varied both spatially and temporally as a result of local differences in nutrient supply (e.g., Tremblay et al., 2015). We measured relatively high surface Chl a concentration (>4 µg L–1; Figure 5d) in Hudson Strait (especially in our transit at the end of the cruise) and in many locations in Hudson Bay (particularly close to the Nelson and Churchill Rivers and in the NW Polynya), associated with oxygen supersaturation (>120%; Figure 5c) and depleted nitrate (NO3) concentrations (<1 mmol m–3; Figure S4), but variable pCO2 values (Figure 2d). Because dissolved O2 equilibrates much more quickly with the atmosphere than dissolved CO2, we cannot use O2 saturation to rule out primary production. However, we can use the observed O2 supersaturation as evidence of primary production. For example, the strong correlation among low pCO2 (< 350 μatm), high Chl a, and O2 supersaturation measured in the NW polynya and Hudson Strait, where the measurements were made before significant gas exchange (i.e., the sea ice was actively melting; Figure 2a), suggest that the O2 supersaturation is related to ice-associated blooms and that primary production played a significant role in lowering the pCO2 values in these locations.

By contrast, in the Nelson Estuary, we observed a high Chl a concentration associated with high pCO2 values (>440 μatm), suggesting that the relatively high productivity in this area was not large enough to overcome other processes acting to increase pCO2. This is not surprising as the Nelson River has the highest annual discharge to Hudson Bay (Déry et al., 2011), with large loads of organic compounds and nutrients (e.g., Granskog et al., 2007), which potentially can sustain both the high biological production and the respiration that generates high pCO2 values. Indeed, Demarty et al. (2009) observed the annual average of pCO2 in the Kettle generating station reservoir located on the Nelson River (about 150 km upstream of Hudson Bay) to vary between 550 and 650 μatm, showing the highly supersaturated nature of this river that might not be completely offset with high biological productivity. Furthermore, we observed variable Chl a concentrations close to the ice edge and under the ice, with values varying from less than 1–4 µg L–1, and associated with O2 saturation varying from saturated to supersaturated values.

Similar to our high Chl a observations near the coast, Roff and Legendre (1986) observed higher primary production in these waters compared to offshore waters in Hudson Bay and attributed differences to nutrient replenishment of the upper productive layer through coastal resuspension processes (Griffiths et al., 1981), internal waves (Ingram et al., 1989), or freshwater plume dynamics (Monti et al., 1996). As well, Kuzyk et al. (2010) found that river runoff supports new primary production around the margin of the bay, where river flow is constrained. In contrast, the low surface Chl a (<1 µg L–1) in many locations in Hudson Bay (Figure 5d), suggests either low productivity or we sampled either before or after the spring bloom. Recently, Matthes et al. (2020) observed that 30% of the annual primary production occurs during the ice melt season (May to June) with stratification likely playing a significant role in delaying the bloom onset into the summer months. That said, our findings of low biological production in many areas of Hudson Bay are consistent with several previous studies that described Hudson Bay as an oligotrophic basin characterized by relatively weak biological production due to a short ice-free season and strong stratification (e.g., Legendre et al., 1996; Ferland et al., 2011; Sibert et al., 2011; Lapoussiere et al., 2013).

Although we cannot quantify whether deep primary production affects our surface pCO2 measurements, it seems unlikely to have an impact given the depth of the deep Chl a maximum (which varied from 20 to 40 m based on similarities in CTD profiles across the study area) and its position beneath the strong and persistent mixed layer that exists in the region during spring (Ahmed et al., 2020). Similarly, Ferland et al. (2011) observed deep Chl a maxima at 40–60 m during summer and noted a negative correlation between primary production and stratification strength in the upper water column across Hudson Bay. We therefore conclude that biological production had a limited effect on the observed pCO2 values during our sampling period, except in the NW polynya and Hudson Strait where primary production may have contributed to the observed pCO2 undersaturation.

#### 4.2.3. Ice melt and riverine input

To study the impact of freshwater from sea-ice melt and river runoff on pCO2 variability in Hudson Bay, we plotted temperature versus salinity visualized with pCO2 values throughout the study area (Figure 6). Past studies have shown that sea-ice meltwater is initially undersaturated in pCO2 (e.g., Geilfus et al., 2015) with low seawater temperature and salinity, whereas Arctic river waters tend to be relatively warmer and often supersaturated in pCO2 as a result of degrading terrestrial organic carbon (e.g., Semiletov et al., 2011; Denfeld et al., 2013).

Figure 6.

Temperature–salinity diagram visualized with surface water pCO2. Temperature–salinity scatterplot of all of our data from Hudson Bay, visualized with surface seawater pCO2. The solid lines show the weighted-area average of salinity and temperature, the dashed lines are mixing lines of freshwater sources, and the rectangle highlights the source of high-pCO2 marine waters. pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f6

Figure 6.

Temperature–salinity diagram visualized with surface water pCO2. Temperature–salinity scatterplot of all of our data from Hudson Bay, visualized with surface seawater pCO2. The solid lines show the weighted-area average of salinity and temperature, the dashed lines are mixing lines of freshwater sources, and the rectangle highlights the source of high-pCO2 marine waters. pCO2 = CO2 partial pressure. DOI: https://doi.org/10.1525/elementa.2020.00130.f6

We measured low surface salinity seawater (i.e., <28) associated with colder water (<1 °C) and low pCO2 (<350 μatm) in many locations close to the ice edge in Hudson Bay (Figures 6 and 2b and d), indicating that sea-ice melt contributed to pCO2 undersaturation in these regions. We also measured high temperature (>5 °C) associated with relatively low salinity observations (<29) and high pCO2 (>420 μatm) in some of the coastal regions in Hudson Bay (especially near the Nelson and Churchill Estuaries; Figures 6 and 2b and d), suggesting that river water contributed to pCO2 supersaturation in these regions. Similarly, Capelle et al. (2020) noted that the inputs of terrestrial organic matter into Hudson Bay were a significant driver of CO2 evasion in the coastal surface waters. Interestingly, the observations from the Churchill Estuary were warmer (0.5–11.2 °C) and had relatively lower pCO2 (averaged 350 μatm with a range of 230–480 μatm) compared to the temperature (0.9–8.3 °C) and pCO2 (averaged 416 μatm with a range of 370–470 μatm) observations recorded in the Nelson Estuary (Figure 6). These differences between both estuaries could be attributed to differences in watershed chemistry, associated with the different physiographic regions, which likely influences the biology and chemistry of the two rivers (Rosenberg et al., 2005); river discharge (Déry et al., 2016), which is expected to influence carbon loading and marine biogeochemistry; and/or atmospheric conditions (e.g., air temperature and solar radiation) prior to sampling. For example, water and air temperatures in the Churchill Estuary were about 4 °C and 9 °C warmer, respectively, than in the Nelson Estuary during sampling (Figures 6 and S3a). Also, the Churchill River has higher concentrations of dissolved organic carbon than the Nelson River (Mundy et al., 2010), which is consistent with the Churchill watershed having a lower carbonate rock content and a higher permafrost abundance than the Nelson watershed.

Results from past studies in Hudson Bay are consistent with our observations, with low dissolved inorganic carbon associated with sea-ice melt (Burt et al., 2016), which typically results in low pCO2 values (e.g., Ahmed et al., 2020). Additionally, our observations of high pCO2 values close to river-dominated nearshore waters match the measurements of Else et al. (2008a). Overall, our results suggest that Hudson Bay was undersaturated in pCO2 relative to the atmosphere during our study time as a result of the dominance of low-pCO2 ice-melt waters, with few observations of supersaturated pCO2.

To more closely examine the effect of ice melt and riverine input on pCO2 in Hudson Bay, we visualized the ship transects from ice-covered waters to coastal areas in the NW polynya and the Nelson and Churchill Estuaries (Figure 7). This plot shows typically high pCO2 (>410 μatm) in the coastal domain and low pCO2 (<330 μatm) observed mostly in the offshore waters. Figure 7a confirms that the Nelson and Churchill Rivers were responsible for the supersaturated pCO2 during our sampling time, but their impact diminishes rapidly offshore toward the ice-melt-influenced waters. Similar to the transects in Nelson and Churchill Estuaries, off Chesterfield Inlet, we observed relatively high pCO2 values that decreased rapidly as we moved offshore across the NW polynya (Figure 7b). The undersaturated pCO2 observations in offshore waters of the NW polynya are probably a result of ice melt and/or biological productivity promoted by the delivery of upwelled and/or riverine nutrients. Many of these processes have been observed in past studies of Arctic polynyas (e.g., Yager et al., 1995; Else et al., 2012).

Figure 7.

pCO2 in the Nelson Estuary and NW polynya. pCO2 values in (a) the Nelson and Churchill Estuaries with the sea-ice cover as of July 2, 2018; and (b) the NW polynya with the sea-ice cover as of June 16, 2018. The average atmospheric pCO2 was 408 μatm during our study time. NW polynya = northwest polynya. DOI: https://doi.org/10.1525/elementa.2020.00130.f7

Figure 7.

pCO2 in the Nelson Estuary and NW polynya. pCO2 values in (a) the Nelson and Churchill Estuaries with the sea-ice cover as of July 2, 2018; and (b) the NW polynya with the sea-ice cover as of June 16, 2018. The average atmospheric pCO2 was 408 μatm during our study time. NW polynya = northwest polynya. DOI: https://doi.org/10.1525/elementa.2020.00130.f7

#### 4.2.4. Under-ice processes

We often observed high pCO2 (>415 μatm) associated with relatively high salinity (>31), cold (< –0.5 °C) waters when the ship was navigating through the heavy ice cover (Figure 7). The observed high pCO2 measurements were also associated with high nitrate concentration (Figure S4), suggesting that the under-ice phytoplankton bloom had not started yet when we were sampling at these locations. These high under-ice pCO2 observations are likely remnants of biogeochemical and physical processes during the winter season. For example, community respiration has been observed to exceed photosynthesis during the winter as a result of low incoming solar radiation attenuated by sea ice and the overlying snowpack (Shadwick et al., 2011; Else et al., 2012). Another potential source of the observed high pCO2 might be the expulsion of CO2-enriched brines into the water column during sea-ice formation and the retention of alkalinity by calcium carbonate crystals within the sea ice (Rysgaard et al., 2007; Dieckmann et al., 2008). Overall, Figure 7 shows both CO2 supersaturation and undersaturation in the marginal ice zone, suggesting very high spatial variability in the distributions of sea-ice melt at these locations.

#### 4.2.5. Upwelling

The observed relatively high salinity (>32) and low temperature (<0 °C) on the east side of Southampton Island in Hudson Bay (Figure 7b) and the southwest side of Baffin Island indicate that upwelling contributed to high pCO2 values in these regions (>420 μatm). Upwelling is a wind-driven process that brings the deep, cold, and usually high pCO2 water to the surface, often resulting in oceanic CO2 outgassing. Several past studies (e.g., Prinsenberg, 1986; Kuzyk et al., 2010) have observed upwelling events in the NW region of Hudson Bay during the summer and fall seasons, suggesting that upwelling in this area may be a long-lasting event rather than an episodic phenomenon. Therefore, although the pCO2 supersaturation we observed in association with upwelling during our sampling was limited in spatial extent, it may persist through much of the open water season and possibly occur in other locations within Hudson Bay.

### 4.3. pCO2 regression analysis

To explain statistically the key variables affecting pCO2 in Hudson Bay, we performed a linear regression analysis between pCO2 as the dependent variable and salinity, temperature, Chl a, and FDOM as independent or explanatory variables. We limited our spatial analysis to the ice-free regions (i.e., ice concentrations less than 10%) in Hudson Bay where the majority of our field measurements were collected. As well, we excluded the pCO2 values impacted by the upwelling event in the NW region to ensure that our regression relationship is not biased by this spatially limited phenomenon.

We found that salinity described the major variability in pCO2 (Pearson’s product-moment: r = –.71; p < .05), followed by temperature (r = .65; p < .05; Figure 8 and Table 2). Although Chl a and FDOM measurements showed a significant linear relationship with pCO2, the Pearson’s correlation coefficients were low (Table 2), showing that their impact on the pCO2 variability is minor but statistically significant.

Figure 8.

pCO2 versus salinity and temperature. Scatterplots of the measured pCO2 in the open water regions in Hudson Bay versus sea surface (a) salinity and (b) temperature. pCO2 = CO2 partial pressure; RMSE = root mean square error. DOI: https://doi.org/10.1525/elementa.2020.00130.f8

Figure 8.

pCO2 versus salinity and temperature. Scatterplots of the measured pCO2 in the open water regions in Hudson Bay versus sea surface (a) salinity and (b) temperature. pCO2 = CO2 partial pressure; RMSE = root mean square error. DOI: https://doi.org/10.1525/elementa.2020.00130.f8

Table 2.

Correlation matrix of Pearson’s r correlation coefficientsa for measured surface water variables throughout the open water regions of Hudson Bay.a DOI: https://doi.org/10.1525/elementa.2020.00130.t2

VariablespCO2TemperatureSalinityChl aFDOM
pCO2 — — — —
Temperature 0.65 — — —
Salinity –0.71 –0.87 — —
Chl a 0.18 0.05 –0.14 —
FDOM 0.11 –0.06 0.04 0.10
VariablespCO2TemperatureSalinityChl aFDOM
pCO2 — — — —
Temperature 0.65 — — —
Salinity –0.71 –0.87 — —
Chl a 0.18 0.05 –0.14 —
FDOM 0.11 –0.06 0.04 0.10

pCO2 = CO2 partial pressure; Chl a = chlorophyll a; FDOM = fluorescing dissolved organic matter.

aAll coefficients are based on 10,000 observations and are significant at a 95% confidence interval.

The regression relationship between salinity and pCO2, statistically significant at a 95% confidence level, was

$pCO2=−14.16(±0.17)×Salinity+753.58(±5.22).$
1

A root mean square error (RMSE) of 36 μatm and relatively strong coefficient of determination (R2 = 0.51) indicate that this algorithm is explaining more than half of the variability in pCO2, which is a good result, given that this algorithm operates over a range of about 270 μatm (i.e., pCO2 ranged from 198 to 468 μatm in the ice-free regions). To examine the significance of other variables on pCO2, we added temperature, FDOM, and Chl a, sequentially, and selected the best model based on the lowest RMSE, Akaike information criterion, and Bayesian information criterion. The resulting multivariate relationship combined salinity, temperature, FDOM, and Chl a, with RMSE equal to 35 μatm and R2 of .55:

2

### 4.4. Sea–air CO2 fluxes

To estimate CO2 fluxes, we used a bulk flux equation:

3

where FCO2 is the CO2 flux (mmol CO2 m–2 day–1), α is the solubility of CO2 in seawater (mol m–3 atm–1) as a function of temperature and salinity (Weiss, 1974), and ΔpCO2 is the difference between the measured surface seawater pCO2 and the average atmospheric pCO2 (in μatm). We calculated sea–air CO2 flux scaled by sea-ice concentration (Ci, reported as a percentage) following the studies of Butterworth and Miller (2016) and Prytherch et al. (2017). The gas transfer velocity k (cm h–1) was calculated according to Wanninkhof (2014):

$k=0.251U10n2(Sc/660)−0.5,$
4

where $U10n2$ is the square of wind speed (m s–1) corrected to a height of 10 m, and Sc is the Schmidt number (unitless). Negative sea–air CO2 fluxes represent atmospheric CO2 uptake (oceanic sink), whereas positive fluxes represent oceanic CO2 outgassing to the atmosphere (oceanic source).

As the shipboard wind speed data were not available after July 2, we estimated the CO2 fluxes in Hudson Bay using both shipboard wind speed data and NARR wind data (Figure 9). The NARR wind speeds were usually underestimated compared to the in situ measurements (Figure 9a), with the highest wind speeds (>11 m s–1) during our study time observed in the NW polynya and the Nelson Estuary (Figure 10a). The shipboard wind speeds ranged up to 15 m s–1 with an average of 5.8 ± 2.9 m s–1 (Table S1 and Figure 9a), whereas the NARR wind data showed a similar range with a lower spatial average of 4.4 ± 2.9 m s–1 (Table S1 and Figure 10a). Consequently, CO2 fluxes calculated using the shipboard wind data had a higher magnitude than those calculated using the NARR wind data (Figure 9b). The average CO2 fluxes using the NARR and the ship wind data were –5.8 and –7.8 mmol CO2 m–2 day–1, respectively (Table S1). Therefore, we consider the CO2 fluxes calculated based on the NARR wind data as lower-bound estimates of the oceanic CO2 sink in Hudson Bay.

Figure 9.

Temporal variability of sea–air CO2 fluxes in Hudson Bay. Temporal variability of (a) shipboard 10-m wind speed (blue) and the coincident NARR wind data (black); and (b) sea–air CO2 fluxes based on both wind speed data sets and the AMSR2 daily sea-ice concentrations (gray). NARR = National Centers for Environmental Prediction North American Regional Reanalysis; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f9

Figure 9.

Temporal variability of sea–air CO2 fluxes in Hudson Bay. Temporal variability of (a) shipboard 10-m wind speed (blue) and the coincident NARR wind data (black); and (b) sea–air CO2 fluxes based on both wind speed data sets and the AMSR2 daily sea-ice concentrations (gray). NARR = National Centers for Environmental Prediction North American Regional Reanalysis; AMSR2 = Advanced Microwave Scanning Radiometer 2. DOI: https://doi.org/10.1525/elementa.2020.00130.f9

Figure 10.

Spatial variability of sea–air CO2 fluxes in Hudson Bay. Spatial variability of (a) NARR wind speed and (b) sea–air CO2 fluxes based on NARR wind data. NARR = National Centers for Environmental Prediction North American Regional Reanalysis. DOI: https://doi.org/10.1525/elementa.2020.00130.f10

Figure 10.

Spatial variability of sea–air CO2 fluxes in Hudson Bay. Spatial variability of (a) NARR wind speed and (b) sea–air CO2 fluxes based on NARR wind data. NARR = National Centers for Environmental Prediction North American Regional Reanalysis. DOI: https://doi.org/10.1525/elementa.2020.00130.f10

The calculated ΔpCO2 varied from about –283 μatm to 242 μatm, with an average of –91.5 ± 60.9 μatm (i.e., the surface waters were generally undersaturated, with respect to the atmosphere). During our study period in Hudson Bay, we observed supersaturated values of ΔpCO2 (>10 μatm) mainly close to the Nelson and Churchill Estuaries and along the southeast of Southampton Island (Figure S3b). Whereas, we observed strongly undersaturated values of ΔpCO2 (< –100 μatm) mostly in Hudson Strait and along the east side of the bay and in the NW polynya.

The high variation in the CO2 fluxes shown in Figure 10b is attributed mostly to the variability in seawater pCO2 resulting from the regional and dynamic variability of ice melt and riverine input (as discussed in Section 4.2). However, the high CO2 uptake observed in the NW polynya was associated with high wind speeds (Figure 10a) over ice-free waters (Figure 2a). To examine this spatial variability of sea–air CO2 fluxes more closely, we visualized the ship transects in the Nelson and Churchill Estuaries and the NW polynya (Figure S6). Overall, the strong CO2 outgassing observed in these locations (Figure S6) was correlated with the supersaturated pCO2 (Figure 2d), resulting in positive ΔpCO2 values (Figure S3b) calculated for these ice-free areas.

Using all of the observations collected during our study, we found that Hudson Bay overall acted as a moderate to weak net sink for atmospheric CO2 with an average CO2 uptake of 5.1 (±9.3) mmol CO2 m–2 day–1 (Table 1). We estimated the uncertainty in the CO2 flux to be 38% (±1.9) mmol CO2 m–2 day–1 based on the cumulative uncertainties derived from the field observations, the gas exchange parameterization, and sea-ice concentration (see Text S2 for additional details on the uncertainty calculations). The moderate average sink of CO2 that we calculated is a result of the dominance of uptake versus outgassing over the study area. Using the 1-km spatial bin analysis, 90% of the sampled area acted as a sink for CO2, while only 10% of the area acted as a source of CO2. Varying the source of sea-ice data did not seem to have a significant effect on these results; using the CISDA ice data in Equation 3, we calculated the average CO2 uptake to be 4.8 (±9.0) mmol CO2 m–2 day–1.

### 4.5. CO2 budget and comparisons to other Arctic coastal seas

While reporting the instantaneous CO2 fluxes based on ship observations, understanding the hourly to daily trends and patterns in spring and early summer seasons is important (Figures 9b and 10b); otherwise, estimates might not be representative of the overall source-sink status of the bay. The problem with extrapolating ship observations to calculate CO2 fluxes over space and time is that the processes influencing seawater pCO2 are variable, hence biasing estimates to the sampling times and locations. For example, the biological activity and the magnitude of river inputs almost certainly change over time and space beyond our sampling survey. Although this spatial–temporal variability should be kept in mind, we expect that our resulting errors will be relatively low, given that Hudson Bay is known as an oligotrophic basin with low biological activity (e.g., Ferland et al., 2011) and that river discharge to Hudson Bay has relatively consistent seasonal hydrographs (Déry et al., 2016; 2018). Assuming relatively consistent or slow change of freshwater input and biological production, we can extrapolate the average CO2 fluxes over space and time to estimate the sea–air CO2 fluxes in the open water areas in Hudson Bay. This approach of extrapolating ship CO2 fluxes has been adopted by many studies in the Arctic regions to calculate the regional and annual CO2 fluxes (e.g., Bates and Mathis, 2009, and references therein).

Considering the sea-ice coverage during our sampling time, we estimated the area of open water in each month separately (i.e., May, June, and July in 2018) using the CISDA ice data. Using the average sink of –5.1 mmol CO2 m–2 day–1 and the calculated open water areas, we estimated the total CO2 fluxes for Hudson Bay in May, June, and July to be –0.3, –1.1, and –1.9 Tg C (1012g C), respectively, for a total net sea–air CO2 flux of –3.3 (± 1.2) Tg C.

Compared to other Arctic shelf regions, Hudson Bay during the spring and early summer seasons acts as a relatively moderate to weak CO2 sink (Table 3). Our calculated average uptake rates are similar to rates estimated for the Laptev and East Siberian Seas (Pipko et al., 2017), which are shelf seas similarly affected by strong river runoff. The calculated rates are also close to estimates for the Canadian Arctic Archipelago, which like Hudson Bay overall experiences low primary production (Ahmed and Else, 2019). In contrast, the Chukchi Sea (which has high rates of primary and net community production during the ice-free period) experiences much stronger CO2 uptake (Bates, 2006). The high CO2 uptake observed in the Kara sea (a shelf sea that receives a lot of riverine input) compared to Hudson Bay was more likely due to sampling close to sea-ice meltwater under high wind speeds (Pipko et al., 2011).

Table 3.

Comparison of CO2 fluxes in this study with fluxes in other Arctic shelf seas. Areas, depths, river inflow, and sea–air CO2 flux of the Arctic continental shelves. DOI: https://doi.org/10.1525/elementa.2020.00130.t3

Shelf SeaArea (103 km2)Mean Depth (m)River Inflow (Km3 yr–1)Sea–Air CO2 Flux (mmol m–2 day–1)aSeasonReference
Barents Sea 1,512 200 463 –11.1 Annual Lauvset et al. (2013)
Kara Sea 926 131 1,133 –18.3 to –32.8 Summer–Fall Pipko et al. (2017)
Laptev Sea 498 48 767 –0.8 to –15.7 Summer–Fall Pipko et al. (2017)
E. Siberian Sea 987 58 213 0.8 to 11.5 Summer Pipko et al. (2011)
Chukchi Sea 620 80 78 –14.8 Annual Bates (2006)
Beaufort Sea 178 124 330 –10.0 Summer Murata and Takizawa (2003)
Canadian Archipelago 1,490 290 270 –3.0 Annual Ahmed and Else (2019)
Hudson Bay 841 150 900 –0.73 Fall Else et al. (2008a)
Hudson Bay 841 150 900 1.98 Summer–Fall Else et al. (2008b)
Hudson Bay and Hudson Strait 1,041 150 900 –5.1 Spring–Early summer Present study
Hudson Bay and Hudson Strait 1,041 150 900 –4.3 Open water Present study
Shelf SeaArea (103 km2)Mean Depth (m)River Inflow (Km3 yr–1)Sea–Air CO2 Flux (mmol m–2 day–1)aSeasonReference
Barents Sea 1,512 200 463 –11.1 Annual Lauvset et al. (2013)
Kara Sea 926 131 1,133 –18.3 to –32.8 Summer–Fall Pipko et al. (2017)
Laptev Sea 498 48 767 –0.8 to –15.7 Summer–Fall Pipko et al. (2017)
E. Siberian Sea 987 58 213 0.8 to 11.5 Summer Pipko et al. (2011)
Chukchi Sea 620 80 78 –14.8 Annual Bates (2006)
Beaufort Sea 178 124 330 –10.0 Summer Murata and Takizawa (2003)
Canadian Archipelago 1,490 290 270 –3.0 Annual Ahmed and Else (2019)
Hudson Bay 841 150 900 –0.73 Fall Else et al. (2008a)
Hudson Bay 841 150 900 1.98 Summer–Fall Else et al. (2008b)
Hudson Bay and Hudson Strait 1,041 150 900 –5.1 Spring–Early summer Present study
Hudson Bay and Hudson Strait 1,041 150 900 –4.3 Open water Present study

aNegative CO2 flux indicates oceanic sink.

Our estimates of CO2 flux through May, June, and July complement past research in Hudson Bay that calculated CO2 fluxes in late summer and autumn (Else et al., 2008a, 2008b). Else et al. (2008a) found that warm sea surface temperature in late summer (August) causes weak outgassing, which eventually reverts to a weak CO2 sink as the surface waters cool through the fall. Combined with the moderate uptake we observed during spring and early summer, a fairly complete picture emerges that is consistent with the typical seasonal CO2 cycle for high-latitude shelf seas described by Bates and Mathis (2009, and references therein). However, the temporal separation of the two studies (13 years) limits our confidence in any calculation of the total sea–air CO2 flux throughout Hudson Bay over the entire open water season solely from the observational results.

Instead, Else et al. (2008b) found that the variability in seawater temperature is a good predictor for the pCO2 variability during the summer and fall seasons in Hudson Bay. Therefore, to compute the total CO2 fluxes during the entire 2018 open water season in Hudson Bay, we can use satellite-derived sea surface temperature and the average pCO2 during our study (approximately 317 μatm), corrected for temperature following the approach of Takahashi et al. (2002) to estimate pCO2 during the 2018 summer and fall seasons. To estimate the average seawater temperature from August to October in 2018, we used the monthly Level-3 MODIS Aqua (https://oceancolor.gsfc.nasa.gov/), clipped the images to the study area, and resampled them to 1-km spatial resolution. The calculated average water temperatures for August, September, and October months were 5.6 °C, 4.3 °C, and 1.9 °C, respectively. The resulting estimated pCO2 values for August, September, and October were about 402, 380, and 344 μatm, respectively. We calculated k (gas transfer velocity) using the NARR wind data and Equation 4 and estimated the sea–air CO2 fluxes using Equation 3. The estimated CO2 fluxes during the 2018 August, September, and October were –0.6, –1.2, and –2.2 Tg C, respectively. Adding these values to those we estimated for the spring and early summer seasons, we obtain a rough estimate of the total open water CO2 uptake in Hudson Bay of –7.2 Tg C. While a complete annual budget should account for the sea–air CO2 exchange during winter, through-ice fluxes are orders of magnitude smaller than open water fluxes and are not expected to have a large impact on the annual budget (Butterworth and Else, 2018). At an areal extent of about 1 × 106 km2, Hudson Bay and Hudson Strait make up approximately 8% of the Arctic Ocean and its associated shelf seas (approximately 13 × 106 km2). Here, we estimated a net sea–air CO2 flux of –7.2 Tg C during the open water season, suggesting that Hudson Bay and Hudson Strait are likely to account for 4–6% of the Arctic Ocean carbon sink, which previous estimates have shown is on the order of 118–180 Tg C year–1 (e.g., Manizza et al., 2019, and references therein). As a result, Hudson Bay is considered a slightly weaker sink than might be expected based on its surface area likely due to the influence of high-pCO2 river water along the coastal regions. Yet, this contribution remains significant at the pan-Arctic scale and underscores the importance of including Hudson Bay in carbon budgets for the Arctic Ocean.

This study presents the first spring and early summer measurements of pCO2 and estimates of sea–air CO2 fluxes in Hudson Bay and Hudson Strait and thus fills up a major gap in our understanding of sea–air CO2 fluxes in the Arctic/Sub-Arctic shelf seas. During our study, sea surface pCO2 exhibited large spatial variability, with localized pCO2 supersaturation in areas influenced by rivers and upwelling, as well as under the sea ice, and more generalized pCO2 undersaturation in ice-melt-influenced waters. Overall, Hudson Bay was mostly undersaturated in pCO2. Temperature variations, and thereby their effect on pCO2, were small at the beginning of our study, when water mixing and biology activity were the primary drivers of pCO2 variability. Along the southern coast, where most of the relatively warm riverine waters enter the bay, the thermal effect on pCO2 became more important. The contribution of biological production to the variability in pCO2 was minor but nonnegligible, especially in the NW polynya and Hudson Strait. In the open water areas, pCO2 showed the highest correlations with salinity and temperature, with minor but statistically significant correlations with Chl a and FDOM. Our observations in the Nelson and Churchill Estuaries highlight the need to better understand variability in the biogeochemical makeup of rivers entering Hudson Bay and their impacts on the receiving marine system. Hudson Bay as a whole acted as a CO2 sink with an average flux of about –5 mmol CO2 m–2 day–1 (–3.3 Tg C) during the study, with an estimated uncertainty of 38% that corresponds to about ±1.9 mmol CO2 m–2 day–1. The calculated sink in this study is significantly larger than previous estimates for Hudson Bay in late summer and autumn, underlining the significance of properly accounting for seasonal variability in Arctic shelf seas. Nevertheless, on average, the bay was a weaker CO2 sink than Arctic shelves due to the strong influence of high-pCO2 river water. We estimated the total CO2 uptake in Hudson Bay during open water season to be –7.2 Tg C, contributing about 4%–6% of the annual CO2 flux estimated for the Arctic Ocean. Future work will attempt to use remotely sensed data and machine learning techniques to extend our observations and generate a more comprehensive estimate of regional sea–air CO2 fluxes in Hudson Bay.

The underway data presented in this article were uploaded to the Surface Ocean CO2 Atlas (SOCAT) and will be available to the public in SOCAT Version 2021.

Figure S1. Sea-ice concentration from daily AMSR2 and weekly CISDA sea-ice products and ship-recorded observations. Sea-ice concentration (%) from daily AMSR2 products (blue), weekly CISDA sea-ice products (black), and in situ observations recorded by the ship officers (red circles). AMSR2 = Advanced Microwave Scanning Radiometer 2; CISDA = Canadian Ice Service Digital Archive.

Figure S2. Regional maps of changes in sea-ice cover over time. Sea-ice concentration (%, color-coded) at approximately the start (May 28), middle (June 18), and end (July 16) of our cruise in Hudson Bay. Sea-ice concentrations were provided by the Canadian Ice Service. Dates are year-month-day.

Figure S3. Air temperature and sea–air CO2 gradient in Hudson Bay. Spatial variability of (a) NARR air temperature and (b) the measured sea–air gradient of CO2 partial pressure (ΔpCO2) during our study time (May 25 to July 13, 2018).

Figure S4. Nitrate and SPAR variability in Hudson Bay. Surface distribution of (a) nitrate (NO3) concentrations and (b) surface photosynthetic active radiation (SPAR) across Hudson Bay during our study. The dots indicate locations of CTD casts, and the white area represents sea-ice concentration greater than 90%, as of July 9, 2018, based on weekly ice charts provided by the Canadian Ice Service. SPAR = surface photosynthetic active radiation; CTD = conductivity-temperature-depth.

Figure S5. pCO2 versus FDOM. Scatterplot of pCO2 versus fluorescent dissolved organic matter (FDOM) in Hudson Bay during our study, with the measurements from the SW coast in red. pCO2 = CO2 partial pressure.

Figure S6. CO2 fluxes in the Nelson Estuary and NW polynya. Spatial variability of sea–air CO2 fluxes in (a) the Nelson and Churchill Estuaries with the sea-ice cover as of July 2, 2018, and (b) the NW polynya with the sea-ice cover as of June 16, 2018. NW polynya = northwest polynya.

Table S1. Shipboard and NARR match-up measurements of wind speed, air temperature, and sea–air CO2 fluxes from May 31 to July 2, 2018.a NARR = National Centers for Environmental Prediction North American Regional Reanalysis.

Text S1. Thermal and nonthermal pCO2 calculations. pCO2 = CO2 partial pressure.

Text S2. Uncertainties of CO2 flux estimation.

The authors thank the Captain and the crew of the Canadian research icebreaker CCGS Amundsen for their hard work and help during the BaySys cruise in 2018. They thank Rachel Mandryk and Yekaterina (Kate) Yezhova for their help during sampling; Tonya Burgers and Sebastian Luque for logistics and help during ship mobilization. They also thank Dr Elizabeth Jones and the second anonymous reviewer for their valuable comments and suggestions to improve the article. The ship thermosalinograph and conductivity-temperature-depth data used in this study were collected by systems onboard the Amundsen and made available by the Amundsen Science Program.

This study was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) Collaborative Research and Development (CRD) Program, Manitoba Hydro, and Fisheries and Oceans Canada.

The authors have no competing interests to declare. LAM serves as Elementa Associate Editor. She was not involved in the review of this article.

Contributed to the conception and design: MMMA, BGTE, TP.

Contributed to the acquisition of data: MMMA, DWC, BB, TP.

Contributed to the analysis and making the figures: MMMA.

Contributed to the interpretation of results: All authors.

Drafted the article: MMMA.

Revised the article: All authors.

Approved the submitted version for publication: All authors.

Ahmed
,
M
,
Else
,
BGT
.
2019
.
The ocean CO2 sink in the Canadian Arctic Archipelago: A present day budget and past trends due to climate change
.
Geophysical Research Letters
46
:
9777
9785
. DOI: http://dx.doi.org/10.1029/2019gl083547.
Ahmed
,
M
,
Else
,
BGT
,
Burgers
,
TM
,
Papakyriakou
,
T
.
2019
.
Variability of surface water pCO2 in the Canadian Arctic Archipelago from 2010 to 2016
.
Journal of Geophysical Research: Oceans
124
(
3
):
1876
1896
. DOI: http://dx.doi.org/10.1029/2018jc014639.
Ahmed
,
MMM
,
Else
,
BGT
,
Capelle
,
DW
,
Miller
,
LA
,
Papakyriakou
,
T
.
2020
.
Underestimation of surface pCO2 and air-sea CO2 fluxes due to freshwater stratification in an Arctic Shelf Sea, Hudson Bay
.
Elementa: Science of the Anthropocene
8
(
1
):
084
. DOI: http://dx.doi.org/10.1525/elementa.084.
Amundsen Science Data Collection
.
2019
.
CTD data collected by the CCGS Amundsen in the Canadian Arctic
.
ArcticNet Inc.
,
.
. DOI: http://dx.doi.org/10.5884/12713.
Available at
www.polardata.ca.
Accessed 2 November 2019
.
Arctic Climate Impact Assessment.
2005
.
Arctic climate impact assessment
.
Cambridge, UK
:
Cambridge University Press:
1042
.
Arrigo
,
KR
,
van Dijken
,
GL
.
2015
.
Continued increases in Arctic Ocean primary production
.
Progress in Oceanography
136
:
60
70
. DOI: http://dx.doi.org/10.1016/j.pocean.2015.05.002.
Bates
,
NR
.
2006
.
Air–sea CO2 fluxes and the continental shelf pump of carbon in the Chukchi Sea adjacent to the Arctic Ocean
.
Journal of Geophysical Research: Oceans
111
(
C10
). DOI: http://dx.doi.org/10.1029/2005jc003083.
Bates
,
NR
,
Mathis
,
J
.
2009
.
The Arctic Ocean marine carbon cycle: Evaluation of air-sea CO2 exchanges, ocean acidification impacts and potential feedbacks
.
Biogeosciences
6
(
11
):
2433
2459
. DOI: http://dx.doi.org/10.5194/bg-6-2433-2009.
Bauer
,
JE
,
Cai
,
WJ
,
Raymond
,
PA
,
Bianchi
,
TS
,
Hopkinson
,
CS
,
Regnier
,
PA
.
2013
.
The changing carbon cycle of the coastal ocean
.
Nature
504
(
7478
):
61
70
. DOI: http://dx.doi.org/10.1038/nature12857.
Borges
,
A
,
Delille
,
B
,
Frankignoulle
,
M
.
2005
.
Budgeting sinks and sources of CO2 in the coastal ocean: Diversity of ecosystems counts
.
Geophysical Research Letters
32
(
14
):
L1. 4601
. DOI: http://dx.doi.org/10.1029/2005gl023053.
Bourgeois
,
T
,
Orr
,
JC
,
Resplandy
,
L
,
Terhaar
,
J
,
Ethé
,
C
,
Gehlen
,
M
,
Bopp
,
L
.
2016
.
Coastal-ocean uptake of anthropogenic carbon
.
Biogeosciences
13
(
14
):
4167
4185
. DOI: http://dx.doi.org/10.5194/bg-13-4167-2016.
Breitburg
,
D
,
Levin
,
LA
,
Oschlies
,
A
,
Grégoire
,
M
,
Chavez
,
FP
,
Conley
,
DJ
,
Garçon
,
V
,
Gilbert
,
D
,
Gutiérrez
,
D
,
Isensee
,
K
,
Jacinto
,
GS
,
Limburg
,
KE
,
Montes
,
I
,
Naqvi
,
SWAAA
,
Pitcher
,
GC
,
Rabalais
,
NN
,
Roman
,
MR
,
Rose
,
KA
,
Seibel
,
BA
,
Telszewski
,
M
,
Yasuhara
,
M
,
Zhang
,
J
.
2018
.
Declining oxygen in the global ocean and coastal waters
.
Science
359
(
6371
). DOI: http://dx.doi.org/10.1126/science.aam7240.
Burt
,
WJ
,
Thomas
,
H
,
Miller
,
LA
,
Granskog
,
MA
,
Papakyriakou
,
TN
,
Pengelly
,
L
.
2016
.
Inorganic carbon cycling and biogeochemical processes in an Arctic inland sea (Hudson Bay)
.
Biogeosciences
13
(
16
):
4659
4671
. DOI: http://dx.doi.org/10.5194/bg-13-4659-2016.
Butterworth
,
BJ
,
Else
,
BGT
.
2018
.
Dried, closed-path eddy covariance method for measuring carbon dioxide flux over sea ice
.
Atmospheric Measurement Techniques
11
:
6075
6090
. DOI: http://dx.doi.org/10.5194/amt-11-6075-2018.
Butterworth
,
BJ
,
Miller
,
SD
.
2016
.
Air-sea exchange of carbon dioxide in the Southern Ocean and Antarctic marginal ice zone
.
Geophysical Research Letters
43
(
13
):
7223
7230
. DOI: http://dx.doi.org/10.1002/2016GL069581.
Cai
,
WJ
,
Dai
,
M
,
Wang
,
Y
.
2006
.
Air–sea exchange of carbon dioxide in ocean margins: A province-based synthesis
.
Geophysical Research Letters
33
(
12
). DOI: http://dx.doi.org/10.1029/2006gl026219.
Capelle
,
DW
,
Kuzyk
,
ZZA
,
Papakyriakou
,
T
,
Guéguen
,
C
,
Miller
,
LA
,
Macdonald
,
RW
.
2020
.
Effect of terrestrial organic matter on ocean acidification and CO2 flux in an Arctic Shelf Sea
.
Progress in Oceanography
185
(
2020
):
102319
. DOI: http://dx.doi.org/10.1016/j.pocean.2020.102319.
Carmack
,
EC
,
Yamamoto-Kawai
,
M
,
Haine
,
T
,
Bacon
,
S
,
Bluhm
,
B
,
Lique
,
C
,
Melling
,
H
,
Polyakov
,
IV
,
Straneo
,
F
,
Timmermans
,
ML
,
Williams
,
WJ
.
2016
.
Fresh water and its role in the Arctic Marine System: Sources, disposition, storage, export, and physical and biogeochemical consequences in the Arctic and global oceans
.
Journal of Geophysical Research:
Biogeosciences
121
(
3
):
675
717
. DOI: http://dx.doi.org/10.1002/2015jg003140.
Chavez
,
F
,
Strutton
,
P
,
Friederich
,
G
,
Feely
,
R
,
Feldman
,
G
,
Foley
,
DG
,
,
MJ
.
1999
.
Biological and chemical response of the equatorial Pacific Ocean to the 1997–98 El Niño
.
Science
286
(
5447
):
2126
2131
. DOI: http://dx.doi.org/10.1126/science.286.5447.2126.
Chen
,
C-TA
,
Borges
,
AV
.
2009
.
Reconciling opposing views on carbon cycling in the coastal ocean: Continental shelves as sinks and near-shore ecosystems as sources of atmospheric CO2
.
Deep Sea Research Part II: Topical Studies in Oceanography
56
(8–10):
578
590
. DOI: http://dx.doi.org/10.1016/j.dsr2.2009.01.001.
Chen
,
C-TA
,
Huang
,
T-H
,
Chen
,
Y-C
,
Bai
,
Y
,
He
,
X
,
Kang
,
Y
.
2013
.
Air-sea exchanges of CO2 in the world’s coastal seas
.
Biogeosciences
10
(
10
). DOI: http://dx.doi.org/10.5194/bg-10-6509-2013.
Dai
,
M
,
Yin
,
Z
,
Meng
,
F
,
Liu
,
Q
,
Cai
,
W-J
.
2012
.
Spatial distribution of riverine DOC inputs to the ocean: An updated global synthesis
.
Current Opinion in Environmental Sustainability
4
(
2
):
170
178
. DOI: http://dx.doi.org/10.1016/j.cosust.2012.03.003.
Dainard
,
PG
,
Guéguen
,
C
,
McDonald
,
N
,
Williams
,
WJ
.
2015
.
Photobleaching of fluorescent dissolved organic matter in Beaufort Sea and North Atlantic Subtropical Gyre
.
Marine Chemistry
177
:
630
637
. DOI: http://dx.doi.org/10.1016/j.marchem.2015.10.004.
Demarty
,
M
,
Bastien
,
J
,
Tremblay
,
A
,
Hesslein
,
RH
,
Gill
,
R
.
2009
.
Greenhouse gas emissions from boreal reservoirs in Manitoba and Québec, Canada, measured with automated systems
.
Environmental Science & Technology
43
(
23
):
8908
8915
. DOI: http://dx.doi.org/10.1021/es8035658.
Denfeld
,
BA
,
Frey
,
KE
,
Sobczak
,
WV
,
Mann
,
PJ
,
Holmes
,
RM
.
2013
.
Summer CO2 evasion from streams and rivers in the Kolyma River basin, north-east Siberia
.
Polar Research
32
(
1
):
19704
. DOI: http://dx.doi.org/10.3402/polar.v32i0.19704.
Déry
,
SJ
,
Mlynowski
,
TJ
,
Hernández-Henríquez
,
MA
,
Straneo
,
F
.
2011
.
Interannual variability and interdecadal trends in Hudson Bay streamflow
.
Journal of Marine Systems
88
(
3
):
341
351
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2010.12.002.
Déry
,
SJ
,
,
TA
,
MacDonald
,
MK
,
Gauli-Sharma
,
B
.
2016
.
Recent trends and variability in river discharge across northern Canada
.
Hydrology and Earth System Sciences
20
(
12
):
4801
4818
. DOI: http://dx.doi.org/10.5194/hess-20-4801-2016.
Déry
,
SJ
,
,
TA
,
MacDonald
,
MK
,
Koenig
,
KA
,
Guay
,
C
.
2018
.
Flow alteration impacts on Hudson Bay river discharge
.
Hydrological Process
32
(
24
):
3576
3587
. DOI: http://dx.doi.org/10.1002/hyp.13285.
Dickson
,
AG
,
Sabine
,
CL
,
Christian
,
JR
.
2007
.
Guide to best practices for ocean CO2 measurements. PICES special publication
.
Sidney, BC
:
North Pacific Marine Science Organization
:
191
,
Vol. 3
.
Dieckmann
,
GS
,
Nehrke
,
G
,
,
S
,
Göttlicher
,
J
,
Steininger
,
R
,
Kennedy
,
H
,
,
D
,
Thomas
,
D
.
2008
.
Calcium carbonate as ikaite crystals in Antarctic sea ice
.
Geophysical Research Letters
35
(
8
). DOI: http://dx.doi.org/10.1029/2008gl033540.
Drinkwater
,
K
,
Jones
,
E
.
1987
.
Density stratification, nutrient and chlorophyll distributions in the Hudson Strait region during summer and their relation to tidal mixing
.
Continental Shelf Research
7
(
6
):
599
607
. DOI: http://dx.doi.org/10.1016/0278-4343(87)90025-2.
Ducklow
,
HW
,
Steinberg
,
DK
,
Buesseler
,
KO
.
2001
.
Upper ocean carbon export and the biological pump
.
Oceanography
14
(
4
):
50
58
. DOI: http://dx.doi.org/10.5670/oceanog.2001.06.
Else
,
BGT
,
Papakyriakou
,
TN
,
Galley
,
RJ
,
Mucci
,
A
,
Gosselin
,
M
,
Miller
,
LA
,
,
EH
,
Thomas
,
H
.
2012
.
Annual cycles of pCO2sw in the southeastern Beaufort Sea: New understandings of air-sea CO2 exchange in arctic polynya regions
.
Journal of Geophysical Research: Oceans
117
(
C9
):
C00G13
. DOI: http://dx.doi.org/10.1029/2011jc007346.
Else
,
BGT
,
Papakyriakou
,
TN
,
Granskog
,
MA
,
Yackel
,
JJ
.
2008
a.
Observations of sea surface fCO2 distributions and estimated air–sea CO2 fluxes in the Hudson Bay region (Canada) during the open water season
.
Journal of Geophysical Research: Oceans
113
(
C8
):
C08026
. DOI: http://dx.doi.org/10.1029/2007jc004389.
Else
,
BGT
,
Yackel
,
JJ
,
Papakyriakou
,
TN
.
2008
b.
Application of satellite remote sensing techniques for estimating air-sea CO2 fluxes in Hudson Bay, Canada during the ice-free season
.
Remote Sensing of Environment
112
(
9
):
3550
3562
. DOI: http://dx.doi.org/10.1016/j.rse.2008.04.013.
Ferland
,
J
,
Gosselin
,
M
,
Starr
,
M
.
2011
.
Environmental control of summer primary production in the Hudson Bay system: The role of stratification
.
Journal of Marine Systems
88
(
3
):
385
400
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2011.03.015.
Friedlingstein
,
P
,
Jones
,
M
,
O’sullivan
,
M
,
Andrew
,
R
,
Hauck
,
J
,
Peters
,
GP
,
Peters
,
W
,
Pongratz
,
J
,
Sitch
,
S
,
Le Quéré
,
C
,
Bakker
,
DCE
,
,
JG
,
Ciais
,
P
,
Jackson
,
RB
,
Anthoni
,
P
,
Barbero
,
L
,
Bastos
,
A
,
Bastrikov
,
V
,
Becker
,
M
,
Bopp
,
L
,
Buitenhuis
,
E
,
Chandra
,
N
,
Chevallier
,
F
,
Chini
,
LP
,
Currie
,
KI
,
Feely
,
RA
,
Gehlen
,
M
,
Gilfillan
,
D
,
Gkritzalis
,
T
,
Goll
,
DS
,
Gruber
,
N
,
Gutekunst
,
S
,
Harris
,
I
,
Haverd
,
V
,
Houghton
,
RA
,
Hurtt
,
G
,
Ilyina
,
T
,
Jain
,
AK
,
Joetzjer
,
E
,
Kaplan
,
JO
,
Kato
,
E
,
Goldewijk
,
KK
,
Korsbakken
,
JI
,
Landschützer
,
P
,
Lauvset
,
SK
,
Lefèvre
,
N
,
Lenton
,
A
,
Lienert
,
S
,
Lombardozzi
,
D
,
Marland
,
G
,
McGuire
,
PC
,
Melton
,
JR
,
Metzl
,
N
,
Munro
,
DR
,
Nabel
,
JEMS
,
Nakaoka
,
S-I
,
Neill
,
C
,
Omar
,
AM
,
Ono
,
T
,
Peregon
,
A
,
Pierrot
,
D
,
Poulter
,
B
,
Rehder
,
G
,
Resplandy
,
L
,
Robertson
,
E
,
Rödenbeck
,
C
,
Séférian
,
R
,
Schwinger
,
J
,
Smith
,
N
,
Tans
,
PP
,
Tian
,
H
,
Tilbrook
,
B
,
Tubiello
,
FN
,
van der Werf
,
GR
,
Wiltshire
,
AJ
,
Zaehle
,
S
.
2019
.
Global carbon budget 2019
.
Earth System Science Data
11
(
4
):
1783
1838
. DOI: http://dx.doi.org/10.5194/essd-11-1783-2019.
Gattuso
J-P
,
Frankignoulle M
,
Wollast R
.
1998
.
Carbon and carbonate metabolism in coastal aquatic ecosystems
.
Annual Review of Ecology and Systematics
29
:
405
434
. DOI: http://dx.doi.org/10.1146/annurev.ecolsys.29.1.405.
Gattuso
,
J-P
,
Hansson
L
.
2011
.
Ocean acidification
.
Oxford, UK
:
Oxford University Press
:
352
.
Geilfus
,
N-X
,
Galley
,
R
,
Crabeck
,
O
,
Papakyriakou
,
T
,
Landy
,
J
,
Tison
,
JL
,
Rysgaard
,
S
.
2015
.
Inorganic carbon dynamics of melt-pond-covered first-year sea ice in the Canadian Arctic
.
Biogeosciences
12
(
6
):
2047
2061
. DOI: http://dx.doi.org/10.5194/bg-12-2047-2015.
Godin
,
P
,
Macdonald
,
RW
,
Kuzyk
,
ZZA
,
Goñi
,
MA
,
Stern
,
GA
.
2017
.
Organic matter compositions of rivers draining into Hudson Bay: Present-day trends and potential as recorders of future climate change
.
Journal of Geophysical Research: Biogeosciences
122
(
7
):
1848
1869
.
Gough
,
WA
,
Gagnon
,
AS
,
Lau
,
HP
.
2004
.
Interannual variability of Hudson Bay ice thickness
.
Polar Geography
28
(
3
):
222
238
. DOI: http://dx.doi.org/10.1080/789610188.
Gough
,
WA
,
Leung
,
A
.
2002
.
Nature and fate of Hudson Bay permafrost
.
Regional Environmental Change
2
(
4
):
177
184
. DOI: http://dx.doi.org/10.1007/s10113-002-0048-x.
Granskog
,
MA
,
Kuzyk
,
ZZA
,
Azetsu-Scott
,
K
,
Macdonald
,
RW
.
2011
.
Distributions of runoff, sea-ice melt and brine using δ18O and salinity data: A new view on freshwater cycling in Hudson Bay
.
Journal of Marine Systems
88
(
3
):
362
374
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2011.03.011.
Granskog
,
MA
,
Macdonald
,
RW
,
Mundy
,
C-J
,
Barber
,
DG
.
2007
.
Distribution, characteristics and potential impacts of chromophoric dissolved organic matter (CDOM) in Hudson Strait and Hudson Bay, Canada
.
Continental Shelf Research
27
(
15
):
2032
2050
. DOI: http://dx.doi.org/10.1016/j.csr.2007.05.001.
Griffiths
,
D
,
Pingree
,
R
,
Sinclair
,
M
.
1981
.
Summer tidal fronts in the near-arctic regions of Foxe Basin and Hudson Bay
.
Deep Sea Research Part A: Oceanographic Research Papers
28
(
8
):
865
873
. DOI: http://dx.doi.org/10.1016/s0198-0149(81)80006-4.
Guéguen
,
C
,
Granskog
,
MA
,
McCullough
,
G
,
Barber
,
DG
.
2011
.
Characterisation of colored dissolved organic matter in Hudson Bay and Hudson Strait using parallel factor analysis
.
Journal of Marine Systems
88
(
3
):
423
433
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2010.12.001.
Guéguen
,
C
,
Itoh
,
M
,
Kikuchi
,
T
,
Eert
,
J
,
Williams
,
WJ
.
2015
.
Variability in dissolved organic matter optical properties in surface waters in the Amerasian Basin
.
Frontiers in Marine Science
2
:
78
. DOI: http://dx.doi.org/10.3389/fmars.2015.00078.
Hochheim
,
KP
,
Lukovich
,
JV
,
Barber
,
D
.
2011
.
Atmospheric forcing of sea ice in Hudson Bay during the spring period, 1980–2005
.
Journal of Marine Systems
88
(
3
):
476
487
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2011.05.003.
Ingram
,
RG
,
Osler
,
J
,
Legendre
,
L
.
1989
.
Influence of internal wave induced vertical mixing on ice algal production in a highly stratified sound
.
Estuarine, Coastal and Shelf Science
29
(
5
):
435
446
. DOI: http://dx.doi.org/10.1016/0272-7714(89)90078-4.
Ingram
,
RG
,
Prinsenberg
,
S
.
1998
. Coastal oceanography of Hudson Bay and surrounding eastern Canadian Arctic waters, in
Robinson
,
AR
,
Brink
,
KH
eds.,
The sea, the global coastal ocean, processes and methods
.
New York, NY
:
Wiley
:
835
861
.
Jiang
,
LQ
,
Cai
,
WJ
,
Wanninkhof
,
R
,
Wang
,
Y
,
Lüger
,
H
.
2008
.
Air–sea CO2 fluxes on the US South Atlantic Bight: Spatial and seasonal variability
.
Journal of Geophysical Research: Oceans
113
(
C7
):
C07019
. DOI: http://dx.doi.org/10.1029/2007jc004366.
Kuzyk
,
ZZA
,
Macdonald
,
RW
,
Tremblay
,
J-É
,
Stern
,
GA
.
2010
.
Elemental and stable isotopic constraints on river influence and patterns of nitrogen cycling and biological productivity in Hudson Bay
.
Continental Shelf Research
30
(
2
):
163
176
. DOI: http://dx.doi.org/10.1016/j.csr.2009.10.014.
Landy
,
JC
,
Ehn
,
JK
,
Babb
,
DG
,
Thériault
,
N
,
Barber
,
DG
.
2017
.
Sea ice thickness in the Eastern Canadian Arctic: Hudson Bay Complex and Baffin Bay
.
Remote Sensing of Environment
200
:
281
294
. DOI: http://dx.doi.org/10.1016/j.rse.2017.08.019.
Lapoussière
,
A
,
Michel
,
C
,
Gosselin
,
M
,
Poulin
,
M
,
Martin
,
J
,
Tremblay
,
J-E
.
2013
.
Primary production and sinking export during fall in the Hudson Bay system, Canada
.
Continental Shelf Research
52
:
62
72
. DOI: http://dx.doi.org/10.1016/j.csr.2012.10.013.
Laruelle
,
GG
,
Lauerwald
,
R
,
Pfeil
,
B
,
Regnier
,
P
.
2014
.
Regionalized global budget of the CO2 exchange at the air-water interface in continental shelf seas
.
Global Biogeochemical Cycles
28
(
11
):
1199
1214
. DOI: http://dx.doi.org/10.1002/2014gb004832.
Lauvset
,
SK
,
Chierici
,
M
,
Counillon
,
F
,
Omar
,
A
,
Nondal
,
G
,
Johannessen
,
T
,
Olsen
,
A
.
2013
.
Annual and seasonal fCO2 and air-sea CO2 fluxes in the Barents Sea
.
Journal of Marine Systems
113–114
:
62
74
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2012.12.011.
Legendre
,
L
,
Robineau
,
B
,
Gosselin
,
M
,
Michel
,
C
,
Ingram
,
R
,
Fortier
,
L
,
Therriault
,
J
,
Demers
,
S
,
Monti
,
D
.
1996
.
Impact of freshwater on a subarctic coastal ecosystem under seasonal sea ice (southeastern Hudson Bay, Canada) II. Production and export of microalgae
.
Journal of Marine Systems
7
(
2
):
233
250
. DOI: http://dx.doi.org/10.1016/0924-7963(95)00007-0.
Liao
,
E
,
Lu
,
W
,
Yan
,
X-H
,
Jiang
,
Y
,
Kidwell
,
A
.
2015
.
The coastal ocean response to the global warming acceleration and hiatus
.
Scientific Reports
5
:
16630
. DOI: http://dx.doi.org/10.1038/srep16630.
Liu
,
K-K
,
Atkinson
,
L
,
Quiñones
,
R
,
Talaue-McManus
,
L
.
2010
.
Carbon and nutrient fluxes in continental margins: A global synthesis
.
. DOI: http://dx.doi.org/:10.1007/978-3-540-92735-8.
Mallory
,
M
,
Loseto
,
L
,
Ferguson
,
S
.
2010
. The future of Hudson Bay: New directions and research needs, in
Ferguson
,
SH
,
Loseto
,
H
,
Mallory
,
ML
eds.,
A little less Arctic: Top predators in the world’s largest northern inland sea, Hudson Bay
.
Dordrecht, the Netherlands
:
Springer
:
291
304
. DOI: http://dx.doi.org/10.1007/978-90-481-9121-5_14.
Manizza
,
M
,
Menemenlis
,
D
,
Zhang
,
H
,
Miller
,
CE
.
2019
.
Modeling the recent changes in the Arctic Ocean CO2 sink (2006–2013)
.
Global Biogeochemical Cycles
33
(
3
):
420
438
. DOI: http://dx.doi.org/10.1029/2018GB006070.
Markus
,
T
,
Stroeve
,
JC
,
Miller
,
J
.
2009
.
Recent changes in Arctic sea ice melt onset, freezeup, and melt season length
.
Journal of Geophysical Research: Oceans
114
(
C12
). DOI: http://dx.doi.org/10.1029/2009jc005436.
Matthes
,
L
,
Dalman
,
LA
,
Ehn
,
JK
,
Tremblay
,
J-E
,
Lee
,
J
,
Belanger
,
S
,
Peeken
,
I
,
Mundy
,
CJ
.
2020
.
First measurements of spring primary production in central Hudson Bay
.
Paper presented at the Ocean Sciences Meeting 2020
;
San Diego, CA
.
Mesinger
,
F
,
DiMego
,
G
,
Kalnay
,
E
,
Mitchell
,
K
,
Shafran
,
PC
,
Ebisuzaki
,
W
,
Jovic
,
D
,
Woollen
,
J
,
Rogers
,
E
,
Berbery
,
EH
,
Ek
,
MB
,
Fan
,
Y
,
Grumbine
,
R
,
Higgins
,
W
,
Li
,
H
,
Lin
,
Y
,
Manikin
,
G
,
Parrish
,
D
,
Shi
,
W
.
2006
.
North American regional reanalysis
.
Bulletin of the American Meteorological Society
87
(
3
):
343
360
. DOI: http://dx.doi.org/10.1175/bams-87-3-343.
Monti
,
D
,
Legendre
,
L
,
Therriault
,
J-C
,
Demers
,
S
.
1996
.
Horizontal distribution of sea-ice microalgae: environmental control and spatial processes (southeastern Hudson Bay, Canada)
.
Oceanographic Literature Review
11
(
43
):
1125
. DOI: http://dx.doi.org/10.3354/meps133229.
Mundy
,
C
,
Gosselin
,
M
,
Starr
,
M
,
Michel
,
C
.
2010
.
Riverine export and the effects of circulation on dissolved organic carbon in the Hudson Bay system, Canada
.
Limnology and Oceanography
55
(
1
):
315
323
.
Murata
,
A
,
Takizawa
,
T
.
2003
.
Summertime CO2 sinks in shelf and slope waters of the western Arctic Ocean
.
Continental Shelf Research
23
(
8
):
753
776
. DOI: http://dx.doi.org/10.1016/s0278-4343(03)00046-3.
Notz
,
D
,
Stroeve
,
J
.
2016
.
Observed Arctic sea-ice loss directly follows anthropogenic CO2 emission
.
Science
354
(
6313
):
747
750
. DOI: http://dx.doi.org/10.1126/science.aag2345.
Orr
,
JC
,
Fabry
,
VJ
,
Aumont
,
O
,
Bopp
,
L
,
Doney
,
SC
,
Feely
,
RA
,
,
A
,
Gruber
,
N
,
Ishida
,
A
,
Joos
,
F
,
Key
,
RM
,
Lindsay
,
K
,
Maier-Reimer
,
E
,
Matear
,
R
,
Monfray
,
P
,
Mouchet
,
A
,
Najjar
,
RG
,
Plattner
,
GK
,
Rodgers
,
KB
,
Sabine
,
CL
,
Sarmiento
,
JL
,
Schlitzer
,
R
,
Slater
,
RD
,
Totterdell
,
IJ
,
Weirig
,
MF
,
Yamanaka
,
Y
,
Yool
,
A
.
2005
.
Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms
.
Nature
437
(
7059
):
681
686
. DOI: http://dx.doi.org/10.1038/nature04095.
Pierre
,
KAS
,
Louis
,
VLS
,
Schiff
,
SL
,
Lehnherr
,
I
,
Dainard
,
PG
,
Gardner
,
AS
,
Aukes
,
PJ
,
Sharp
,
MJ
.
2019
.
Proglacial freshwaters are significant and previously unrecognized sinks of atmospheric CO2
.
Proceedings of the National Academy of Sciences
116
(
36
):
17690
17695
. DOI: http://dx.doi.org/10.1073/pnas.1904241116.
Pierrot
,
D
,
Neill
,
C
,
Sullivan
,
K
,
Castle
,
R
,
Wanninkhof
,
R
,
Lüger
,
H
,
Johannessen
,
T
,
Olsen
,
A
,
Feely
,
RA
,
Cosca
,
CE
.
2009
.
Recommendations for autonomous underway pCO2 measuring systems and data-reduction routines
.
Deep Sea Research Part II: Topical Studies in Oceanography
56
(
8–10
):
512
522
. DOI: http://dx.doi.org/10.1016/j.dsr2.2008.12.005.
Pipko
,
II
,
Pugach
,
SP
,
Semiletov
,
IP
,
Anderson
,
LG
,
Shakhova
,
NE
,
Gustafsson
,
Ö
,
Repina
,
IA
,
Spivak
,
EA
,
Charkin
,
AN
,
Salyuk
,
AN
.
2017
.
The spatial and interannual dynamics of the surface water carbonate system and air–sea CO2 fluxes in the outer shelf and slope of the Eurasian Arctic Ocean
.
Ocean Science
13
(
6
): 997. DOI: http://dx.doi.org/10.5194/os-13-997-2017.
Pipko
,
II
,
Semiletov
,
I
,
Pugach
,
S
,
Wåhlström
,
I
,
Anderson
,
L
.
2011
.
Interannual variability of air-sea CO2 fluxes and carbon system in the East Siberian Sea
.
Biogeosciences
8
(
7
):
1987
2007
. DOI: http://dx.doi.org/10.5194/bg-8-1987-2011.
Prinsenberg
,
S
.
1986
. The circulation pattern and current structure of Hudson Bay, in
Martini
,
IP
ed.,
.
Amsterdam, the Netherlands
:
Elsevier (Eslevier oceanography series
; vol.
44
):
187
204
. DOI: http://dx.doi.org/10.1016/s0422-9894(08)70903-6.
Prinsenberg
,
S
.
1988
.
Ice-cover and ice-ridge contributions to the freshwater contents of Hudson Bay and Foxe Basin
.
Arctic
41
(
1
):
6
11
. DOI: http://dx.doi.org/10.14430/arctic1686.
Prytherch
,
J
,
Brooks
,
IM
,
Crill
,
PM
,
Thornton
,
BF
,
Salisbury
,
DJ
,
Tjernström
,
M
,
Anderson
,
LG
,
Geibel
,
MC
,
Humborg
,
C
.
2017
.
Direct determination of the air-sea CO2 gas transfer velocity in Arctic sea ice regions
.
Geophysical Research Letters
44
(
8
):
3770
3778
. DOI: http://dx.doi.org/10.1002/2017gl073593.
Regnier
,
P
,
Friedlingstein
,
P
,
Ciais
,
P
,
Mackenzie
,
FT
,
Gruber
,
N
,
Janssens
,
IA
,
Laruelle
,
GG
,
Lauerwald
,
R
,
Luyssaert
,
S
,
,
AJ
,
Arndt
,
S
,
Arnosti
,
C
,
Borges
,
AV
,
Dale
,
AW
,
Gallego-Sala
,
A
,
Goddéris
,
Y
,
Goossens
,
N
,
Hartmann
,
J
,
Heinze
,
C
,
Ilyina
,
T
,
Joos
,
F
,
LaRowe
,
DE
,
Leifeld
,
J
,
Meysman
,
FJR
,
Munhoven
,
G
,
Raymond
,
PA
,
Spahni
,
R
,
Suntharalingam
,
P
,
Thullner
,
M
.
2013
.
Anthropogenic perturbation of the carbon fluxes from land to ocean
.
Nature Geoscience
6
(
8
):
597
607
. DOI: http://dx.doi.org/10.1038/ngeo1830.
Roff
,
J
,
Legendre
,
L
.
1986
. Physico-chemical and biological oceanography of Hudson Bay, in
Martini
,
IP
ed.,
.
Amsterdam, the Netherlands
:
Elsevier (Eslevier oceanography series
; vol.
44
):
265
292
. DOI: http://dx.doi.org/10.1016/s0422-9894(08)70907-3.
Roobaert
,
A
,
Laruelle
,
GG
,
Landschützer
,
P
,
Gruber
,
N
,
Chou
,
L
,
Regnier
,
P
.
2019
.
The spatiotemporal dynamics of the sources and sinks of CO2 in the global coastal ocean
.
Global Biogeochemical Cycles
33
(
12
):
1469
1783
. DOI: http://dx.doi.org/10.1029/2019GB006239.
Rosenberg
,
DM
,
Chambers
,
PA
,
Culp
,
JM
,
Franzin
,
WG
,
Nelson
,
PA
,
Salki
,
AG
,
Stainton
,
MP
,
Bodaly
,
R
,
Newbury
,
RW
.
2005
. Nelson and Churchill river basins, in
Cushing
,
CE
ed.,
Rivers of North America
.
Benke, AC
:
Elsevier
:
852
901
.
Ch. 19
. DOI: http://dx.doi.org/10.1016/b978-012088253-3/50022-5.
Rysgaard
,
S
,
Glud
,
RN
,
Sejr
,
M
,
Bendtsen
,
J
,
Christensen
,
P
.
2007
.
Inorganic carbon transport during sea ice growth and decay: A carbon pump in polar seas
.
Journal of Geophysical Research: Oceans
112
(
C3
). DOI: http://dx.doi.org/10.1029/2006jc003572.
Rysgaard
,
S
,
Mortensen
,
J
,
Juul-Pedersen
,
T
,
Sørensen
,
LL
,
Lennert
,
K
,
Søgaard
,
DH
,
Arendt
,
KE
,
Blicher
,
ME
,
Sejr
,
MK
,
Bendtsen
,
J
.
2012
.
High air-sea CO2 uptake rates in nearshore and shelf areas of Southern Greenland: Temporal and spatial variability
.
Marine Chemistry
128
:
26
33
:
C03016
. DOI: http://dx.doi.org/10.1016/j.marchem.2011.11.002.
Sakshaug
,
E
.
2004
. Primary and secondary production in the Arctic Seas, in
Stein
,
R
,
MacDonald
,
RW
eds.,
The organic carbon cycle in the Arctic Ocean
.
Berlin, Heidelberg
:
Springer
:
57
81
. DOI: http://dx.doi.org/10.1007/978-3-642-18912-8_3.
Saucier
,
F
,
Senneville
,
S
,
Prinsenberg
,
S
,
Roy
,
F
,
Smith
,
G
,
Gachon
,
P
,
Caya
,
D
,
Laprise
,
R
.
2004
.
Modelling the sea ice-ocean seasonal cycle in Hudson Bay, Foxe Basin and Hudson Strait, Canada
.
Climate Dynamics
23
(
3–4
):
303
326
. DOI: http://dx.doi.org/10.1007/s00382-004-0445-6.
Semiletov
,
IP
,
Pipko
,
II
,
Shakhova
,
NE
,
Dudarev
,
OV
,
Pugach
,
SP
,
Charkin
,
AN
,
McRoy
,
CP
,
Kosmach
,
D
,
Gustafsson
,
O
.
2011
.
Carbon transport by the Lena River from its headwaters to the Arctic Ocean, with emphasis on fluvial input of terrestrial particulate organic carbon vs. carbon transport by coastal erosion
.
Biogeosciences
8
(
9
):
2407
. DOI: http://dx.doi.org/10.5194/bg-8-2407-2011.
,
EH
,
Thomas
,
H
,
Chierici
,
M
,
Else
,
B
,
Fransson
,
A
,
Michel
,
C
,
Miller
,
LA
,
Mucci
,
A
,
Niemi
,
A
,
Papakyriakou
,
TN
,
Tremblay
,
JE
.
2011
.
Seasonal variability of the inorganic carbon system in the Amundsen Gulf region of the southeastern Beaufort Sea
.
Limnology and Oceanography
56
(
1
):
303
322
. DOI: http://dx.doi.org/10.4319/lo.2011.56.1.0303.
Shepherd
,
TG
.
2016
.
Effects of a warming Arctic
.
Science
353
(
6303
):
989
990
. DOI: http://dx.doi.org/10.1126/science.aag2349.
Sibert
,
V
,
Zakardjian
,
B
,
Gosselin
,
M
,
Starr
,
M
,
Senneville
,
S
,
LeClainche
,
Y.
2011
.
3D bio-physical model of the sympagic and planktonic productions in the Hudson Bay System
.
Journal of Marine Systems
88
(
3
):
401
422
. DOI: http://dx.doi.org/10.1016/j.jmarsys.2011.03.014.
Spreen
,
G
,
Kaleschke
,
L
,
Heygster
,
G
.
2008
.
Sea ice remote sensing using AMSR-E 89-GHz channels
.
Journal of Geophysical Research: Oceans
113
(
C2
):
C02S03
. DOI: http://dx.doi.org/10.1029/2005jc003384.
Stull
,
RB
.
1988
.
An introduction to boundary layer meteorology
.
Dordrecht, the Netherlands
:
Springer
:
666
. DOI: http://dx.doi.org/10.1007/978-94-009-3027-8_1.
Takahashi
,
T
,
Olafsson
,
J
,
Goddard
,
JG
,
Chipman
,
DW
,
Sutherland
,
S
.
1993
.
Seasonal variation of CO2 and nutrients in the high-latitude surface oceans: A comparative study
.
Global
Biogeochemical Cycles
7
(
4
):
843
878
. DOI: http://dx.doi.org/10.1029/93gb02263.
Takahashi
,
T
,
Sutherland
,
SC
,
Sweeney
,
C
,
Poisson
,
A
,
Metzl
,
N
,
Tilbrook
,
B
,
Bates
,
N
,
Wanninkhof
,
R
,
Feely
,
RA
,
Sabine
,
C
.
2002
.
Global sea-air CO2 flux based on climatological surface ocean pCO2, and seasonal biological and temperature effects
.
Deep Sea Research Part II: Topical Studies in Oceanography
49
(
9
):
1601
1622
. DOI: http://dx.doi.org/10.1016/s0967-0645(02)00003-6.
Takahashi
,
T
,
Sutherland
,
SC
,
Wanninkhof
,
R
,
Sweeney
,
C
,
Feely
,
RA
,
Chipman
,
DW
,
Hales
,
B
,
Friederich
,
G
,
Chavez
,
F
,
Sabine
,
C
.
2009
.
Climatological mean and decadal change in surface ocean pCO2, and net sea–air CO2 flux over the global oceans
.
Deep Sea Research Part II: Topical Studies in Oceanography
56
(8–10):
554
577
. DOI: http://dx.doi.org/10.1016/s0967-0645(02)00003-6.
Tivy
,
A
,
Howell
,
SE
,
Alt
,
B
,
McCourt
,
S
,
Chagnon
,
R
,
Crocker
,
G
,
Carrieres
,
T
,
Yackel
,
JJ
.
2011
.
Trends and variability in summer sea ice cover in the Canadian Arctic based on the Canadian Ice Service Digital Archive, 1960–2008 and 1968–2008
.
Journal of Geophysical Research: Oceans
116
(
C3
):
C03007
. DOI: http://dx.doi.org/10.1029/2009jc005855.
Tremblay
,
J-É
,
Anderson
,
LG
,
Matrai
,
P
,
Coupel
,
P
,
Bélanger
,
S
,
Michel
,
C
,
,
M
.
2015
.
Global and regional drivers of nutrient supply, primary production and CO2 drawdown in the changing Arctic Ocean
.
Progress in Oceanography
139
:
171
196
. DOI: http://dx.doi.org/10.1016/j.pocean.2015.08.009.
Tsunogai
,
S
,
Watanabe
,
S
,
Sato
,
T
.
1999
.
Is there a “continental shelf pump” for the absorption of atmospheric CO2?
Tellus B
51
(
3
):
701
712
. DOI: http://dx.doi.org/10.1034/j.1600-0889.1999.t01-2-00010.x.
Walsh
,
JJ
.
1991
.
Importance of continental margins in the marine biogeochemical cycling of carbon and nitrogen
.
Nature
350
(
6313
):
53
55
. DOI: http://dx.doi.org/10.1038/350053a0.
Wanninkhof
,
R
.
2014
.
Relationship between wind speed and gas exchange over the ocean revisited
.
Limnology and Oceanography: Methods
12
(
6
):
351
362
. DOI: http://dx.doi.org/10.4319/lom.2014.12.351.
Weiss
,
RF
.
1974
.
Carbon dioxide in water and seawater: The solubility of a non-ideal gas
.
Marine Chemistry
2
(
3
):
203
215
. DOI: http://dx.doi.org/10.1016/0304-4203(74)90015-2.
Wollast
,
R
,
1998
. Evaluation and comparison of the global carbon cycle in the coastal zone and in the open ocean, in
Brink
,
KH
,
Robinson
,
AR
eds.,
The sea, the global coastal ocean, processes and methods
.
New York, NY
:
Wiley
:
213
252
.
Yager
,
PL
,
Wallace
,
DW
,
Johnson
,
KM
,
Smith Jr
,
WO
,
Minnett
,
PJ
,
Deming
,
JW
.
1995
.
The Northeast Water Polynya as an atmospheric CO2 sink: A seasonal rectification hypothesis
.
Journal of Geophysical Research: Oceans
100
(
C3
):
4389
4398
. DOI: http://dx.doi.org/10.1029/94jc01962.

How to cite this article: Ahmed, MMM, Else, BGT, Butterworth, B, Capelle, DW, Guéguen, C, Miller, LA, Meilleur, C, Papakyriakou, T. 2021. Widespread surface water pCO2 undersaturation during ice-melt season in an Arctic continental shelf sea (Hudson Bay, Canada). Elementa: Science of Anthropocene 9(1). DOI: https://doi.org/10.1525/elementa.2020.00130

Domain Editor-in-Chief: Jody W. Deming, University of Washington, Seattle, WA, USA

Associate Editor: Jean-Éric Tremblay, Department of Biology, Université Laval, Quebec, Canada

Knowledge Domain: Ocean Science

Part of an Elementa Special Feature: BaySys

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.