Arctic rivers deliver a large amount of organic carbon to Arctic coastal waters each year. The Arctic Ocean receives approximately 11% of global river discharge (Lammers et al., 2001) and 33 Tg y–1 of riverine dissolved organic carbon (DOC; McGuire et al., 2009), resulting in an Arctic “Riverine Coastal Domain” (Carmack et al., 2015) that contains the highest concentrations of land-derived DOC among the world’s oceans (Benner et al., 2005). In recent decades, climate change and human activity have led to increased river discharge (McClelland et al., 2006; McGuire et al., 2009; Déry et al., 2011), accelerated coastal erosion (Mars and Houseknech, 2007), and permafrost thaw (Vonk et al., 2015), all leading to increased release and transport of riverine DOC to Arctic coastal waters. Amon et al. (2012) predict that DOC fluxes from Arctic rivers could increase in a warmer climate even if conditions are otherwise undisturbed.

In this study, we conducted 45-day incubation experiments to examine the biodegradation of DOC during late winter (prefreshet) conditions in riverine and coastal waters of Hudson Bay. Located at the southern margin of the Canadian Arctic, Hudson Bay is even more strongly river-influenced than the Arctic Ocean, with an annual runoff yield of approximately 0.9 m (Granskog et al., 2011). Winter river discharge to Hudson Bay is a significant fraction of the annual total (ca. one-fifth) and increasing during recent decades due in part to ongoing climate change and river regulation for hydropower production (Déry et al., 2016). We use the results of the incubation experiments to assess the biodegradation potential of DOC in late winter and to compare two contrasting Hudson Bay riverine coastal areas differing in amount and composition of DOC and in marine characteristics (e.g., salinity, stratification). The results provide some of the first information regarding biodegradability of wintertime riverine DOC within the Arctic watershed and point to basin properties including river regulation and dissolved nutrient concentrations as variables warranting further exploration in terms of their influence on the biodegradability of arctic riverine DOC.

2.1. Study area

Hudson Bay (including James Bay) is a large (approximately 1.2 × 106 km2) shelf sea that receives and processes Arctic Ocean outflow and then discharges it via Hudson Strait (Figure 1) to the Labrador Sea. Hudson Bay is ice-covered for up to 9 months of each year, with the exception of leads and polynyas, and the surrounding watershed contains the most southerly extent of permafrost in Canada (Figure 1). The location of Hudson Bay at the southern margin of the Arctic makes it more vulnerable to climate change including shortening of the season of sea ice cover (Hochheim and Barber, 2014) and hydrologic changes in the watershed (see also Macdonald and Kuzyk, 2011; Déry et al., 2016). Hudson Bay has a drainage basin of about 3.9 × 106 km2 (Figure 1) and receives river inflow totaling more than 760 km3 yr–1 or about 30% of the total Canadian river runoff (Déry et al., 2011). Winter river discharge to southwestern Hudson and James Bays increased significantly between 1964 and 2013 due in part to hydroelectric development on several large rivers, which has altered natural hydrological cycles so that river flow rates remain high in winter when demand for hydropower production peaks (Déry et al., 2016).

Figure 1.

A map of Hudson Bay and its overall drainage basin and permafrost cover. The individual basins of the Nelson, Hayes, and Great Whale Rivers are shown relative to vegetation zones. Inset photos show the lower Nelson River (A) and the Great Whale River (B). Hudson Bay’s general surface circulation patterns are shown by arrows. The coastal study areas, southwestern Hudson Bay (SWHB) and southeastern Hudson Bay (SEHB), are labeled. DOI: https://doi.org/10.1525/elementa.2020.00123.f1

Figure 1.

A map of Hudson Bay and its overall drainage basin and permafrost cover. The individual basins of the Nelson, Hayes, and Great Whale Rivers are shown relative to vegetation zones. Inset photos show the lower Nelson River (A) and the Great Whale River (B). Hudson Bay’s general surface circulation patterns are shown by arrows. The coastal study areas, southwestern Hudson Bay (SWHB) and southeastern Hudson Bay (SEHB), are labeled. DOI: https://doi.org/10.1525/elementa.2020.00123.f1

Close modal

With an abundance of wetlands, peatlands, and soil carbon in the Hudson Bay watershed (Godin et al., 2017), riverine inputs of DOC to Hudson Bay are very high relative to the surface area of the bay or its seawater volume (Granskog et al., 2011). James Bay and Hudson Bay river waters have been sampled opportunistically during ArcticNet-funded CCGS Amundsen expeditions and during irregular monitoring programs involving industry, government, and community partners. DOC concentrations in the sampled river waters range between 0.24 and 116 mg L–1 (Granskog et al., 2007; Kirk and St. Louis, 2009; Rosa et al., 2012; Godin et al., 2017), with DOC fluxes and composition varying among individual watersheds in relation to climate, vegetation, permafrost distribution, and hydrology (Rosa et al., 2012; Godin et al., 2017). The annual riverine DOC inputs are estimated from the scarce data at approximately 5.5 Tg C (Mundy et al., 2010), which is an order of magnitude more than the inputs of particulate organic carbon (POC; 0.46 ± 0.33 Tg C; Kuzyk et al., 2009). Dominance of DOC over its particulate counterpart in river discharge is expected as many of the large rivers drain the massive permafrost-bound wetland/peatland of the Hudson Bay Lowlands (Godin et al., 2017). Within Hudson Bay, river runoff is introduced to a strong, anticlockwise, surface coastal current system that largely confines the river water to a coastal corridor approximately 100 km wide (Granskog et al., 2009, 2011; St-Laurent et al., 2011).

This study focuses on two subareas within the Hudson Bay riverine coastal domain: the southwestern Hudson Bay (SWHB) and southeastern Hudson Bay (SEHB) study areas (Figure 1). The coastal waters of SWHB receive discharge from the Nelson and Hayes Rivers. These rivers drain mainly flat-lying sedimentary rocks covered by unconsolidated Quaternary sediments (Rosa et al., 2012). The Nelson River has an extensive watershed (>1.1 × 106 km2; Figure 1) that contains many dams and other control structures. Its flow is augmented by partial diversion of the Churchill River resulting in a mean annual discharge of 109 km3 yr–1 (2004–2013) or one-seventh of the total discharge to Hudson Bay and the single largest river discharging to Hudson Bay (Déry et al., 2016). The regulation of the Nelson River results in a relatively flat hydrograph with high winter flows compared to unregulated rivers (Figure 2A). Most of the Nelson River basin is covered by boreal forest (Boreal Shield and Boreal Plains), but Prairies with agricultural land lie in the southwestern part of the basin (Godin et al., 2017; Figure 1). The lower 150 km of the Nelson River crosses the Hudson Plains in a deep valley incised into the flat Hudson Bay Lowlands, with near-vertical, permanently frozen river banks (Rosenberg et al., 2005; see inset photo in Figure 1). The annual DOC flux of the Nelson River represents more than one-fifth of the Bay’s total input at approximately 1.2 Tg C yr–1 (Rosa et al., 2012). DOC concentrations in Nelson River water remain relatively high throughout winter, possibly because of abundant large lakes and wetlands in the catchment that continue to supply DOC when soils are frozen (Rosa et al., 2012). The Hayes River, which contributes to the same estuary as the Nelson River, is a smaller, unregulated river with a mean annual discharge of 21.3 km3 yr–1 (2004–2013; Déry et al., 2016). Its basin is relatively small and, similar to the neighboring Nelson River basin, spans the Boreal Shield and Hudson Plain ecozones, with continuous and discontinuous permafrost covering the whole basin and boreal forest being the dominant type of vegetation (Godin et al., 2017; Figure 1).

Figure 2.

Daily river discharges of the Nelson, Hayes, and Great Whale Rivers. Data are shown throughout a year (A) and in late winter–spring (B) for the Nelson River (daily averages for the period 2000–2016; 56° 23′ 51″N, 94° 22′ 10″W), Hayes River (2017; 56° 25′ 00″N, 92° 48′ 33″W), and Great Whale River (daily averages for the period 2000–2013; 55° 14′ 15″N, 76° 59′ 02″W). The sampling period of this study is highlighted by a yellow bar. Data derived from Environment and Natural Resources, Government of Canada. DOI: https://doi.org/10.1525/elementa.2020.00123.f2

Figure 2.

Daily river discharges of the Nelson, Hayes, and Great Whale Rivers. Data are shown throughout a year (A) and in late winter–spring (B) for the Nelson River (daily averages for the period 2000–2016; 56° 23′ 51″N, 94° 22′ 10″W), Hayes River (2017; 56° 25′ 00″N, 92° 48′ 33″W), and Great Whale River (daily averages for the period 2000–2013; 55° 14′ 15″N, 76° 59′ 02″W). The sampling period of this study is highlighted by a yellow bar. Data derived from Environment and Natural Resources, Government of Canada. DOI: https://doi.org/10.1525/elementa.2020.00123.f2

Close modal

The Nelson and Hayes Rivers discharge together into a large, mostly shallow (approximately 6 m), funnel-shaped estuary, which subsequently discharges onto a shallow, gently sloping marine shelf. The tide in the estuary is semidiurnal and the tidal range varies from 2 to 5 m for neap and spring tides, respectively (Wang et al., 2012). Winter ice cover reduces tidal amplitude and range due to under-ice friction while increasing tidal velocity due to reduced channel cross-section (Wang et al., 2012). Hydraulic residence times for the estuary were estimated by Manitoba Hydro utilizing a calibrated MIKE3 HD model (RSW-Environnement Illimité Inc., 2014) coupled with a 3D transport module from the Danish Hydraulic Institute. The model results indicate that under an average flow of 3,500 m3 s–1, the mean length of time a parcel spends in the estuary varies from <2 days at the Nelson River mouth to 40 days in the outer portion of the estuary (K Sydor, Manitoba Hydro, personal communication, 2019).

The watershed and coastal morphology of SEHB differs markedly from SWHB. Precambrian rocks create a more rolling topography, and small nearshore basins reach water depths of 100 m or more (Hudon et al., 1996). SEHB receives discharge directly from the Great Whale River, which is unregulated and similar in size to the Hayes River, having mean annual discharge of 19 km3 yr–1 for the period of 2004–2013 (Déry et al., 2016). SEHB is also in the path of James Bay outflow (Figure 1), which remains relatively fresh (salinity < 26) throughout the year, due to river inflow, and causes shallow stratification to be maintained throughout winter in SEHB (Eastwood et al., 2020). Tidal amplitudes are smaller in SEHB than in SWHB averaging <2 m. The Great Whale River drains a transitioning mixed forest-tundra zone (Taiga Shield; Figure 1), and the basin is covered by patchy and sporadic discontinuous permafrost (Godin et al., 2017). The annual DOC flux of the Great Whale River is approximately 0.12 Tg C yr–1, an order of magnitude lower than the Nelson River (Rosa et al., 2012).

2.2. Sampling locations and environmental conditions during sampling

Water sample collection and measurement of temperature and salinity for this study were conducted from the river ice and landfast ice (approximately 1 m thick) in SWHB and SEHB between April 9 and April 13, 2017 (Table 1), prior to the start of spring freshet (Figure 2). The discharge of the Nelson River was about 3,500 m3 s–1 during the sampling period (Figure 2). The Nelson River was partially ice-covered between Split Lake and the coast (Sentinel-1 hosted on SIKU.org; not shown), and the last 40 km stretch of river was completely ice-covered (Figure 3A). The Hayes River, which was flowing at approximately 400 m3 s–1 (Figure 2), was ice-covered throughout its whole length (Figure 3A). Due to logistical constraints, surface samples were collected near the mouth of the Nelson and Hayes Rivers (Figure 3A). The Nelson/Hayes Estuary was ice-covered at the time of sampling except for an intermittent area of open water (flaw lead; Figure 3A), which is typical of this estuary throughout winter (Wang et al., 2012). The SWHB coastal sampling sites were evenly distributed over a distance of about 80 km along the landfast ice, paralleling the shoreline in an E-NE direction, and within 5 km of shore (Figure 3A). The water depth at these sampling sites varied from 3 to 5.5 m (Table 1). At the time of sampling, SWHB was receiving only approximately 1,400 µmol m–2 s–1 of photosynthetically active radiation, of which less than 4% was transmitted through the snow and ice cover at our sampling sites (L Dalman, personal communication, 2020). This small fraction implies that at the ice-covered portion of SWHB, riverine and coastal waters were receiving little light and had low photobleaching potential.

Table 1.

Stations sampled in the two study areas, southwestern Hudson Bay and southeastern Hudson Bay, as shown in Figure 3. DOI: https://doi.org/10.1525/elementa.2020.00123.t1

StationLatitudeLongitudeBottom Depth (m)Depth Sampled (m)Date Sampled (2017)
NR 57.033 –92.521 N/A 1.0 April 9
HR 57.008 –92.260 N/A 1.0 April 9
W-1 57.098 –91.972 5.0 1.0 April 9
W-1 57.098 –91.972 5.0 3.0 April 9
W-2 57.158 –91.711 4.0 1.0 April 12
W-3 57.162 –91.714 4.5 1.0 April 12
W-3 57.162 –91.714 4.5 2.0 April 12
W-4 57.223 –91.402 3.0 1.0 April 9
W-4 57.223 –91.402 3.0 3.0 April 9
W-5 57.246 –91.405 5.5 1.0 April 9
W-5 57.246 –91.405 5.5 6.0 April 9
GWR 55.280 –77.642 N/A 1.0 April 12
S-1 55.288 –77.808 41.0 1.0 April 11
S-1 55.288 –77.808 41.0 5.0 April 11
S-1 55.288 –77.808 41.0 10.0 April 11
S-1 55.288 –77.808 41.0 40.0 April 11
S-2 55.294 –77.821 61.0 1.0 April 12
S-3 55.303 –77.824 60.0 1.0 April 12
S-3 55.303 –77.824 60.0 5.0 April 12
S-3 55.303 –77.824 60.0 20.0 April 12
S-3 55.303 –77.824 60.0 50.0 April 12
S-4 55.430 –77.504 36.5 1.0 April 13
S-4 55.430 –77.504 36.5 10.0 April 13
S-4 55.430 –77.504 36.5 35.0 April 13
S-5 55.475 –77.595 31.0 1.0 April 13
S-5 55.475 –77.595 31.0 5.0 April 13
S-5 55.475 –77.595 31.0 30.0 April 13
StationLatitudeLongitudeBottom Depth (m)Depth Sampled (m)Date Sampled (2017)
NR 57.033 –92.521 N/A 1.0 April 9
HR 57.008 –92.260 N/A 1.0 April 9
W-1 57.098 –91.972 5.0 1.0 April 9
W-1 57.098 –91.972 5.0 3.0 April 9
W-2 57.158 –91.711 4.0 1.0 April 12
W-3 57.162 –91.714 4.5 1.0 April 12
W-3 57.162 –91.714 4.5 2.0 April 12
W-4 57.223 –91.402 3.0 1.0 April 9
W-4 57.223 –91.402 3.0 3.0 April 9
W-5 57.246 –91.405 5.5 1.0 April 9
W-5 57.246 –91.405 5.5 6.0 April 9
GWR 55.280 –77.642 N/A 1.0 April 12
S-1 55.288 –77.808 41.0 1.0 April 11
S-1 55.288 –77.808 41.0 5.0 April 11
S-1 55.288 –77.808 41.0 10.0 April 11
S-1 55.288 –77.808 41.0 40.0 April 11
S-2 55.294 –77.821 61.0 1.0 April 12
S-3 55.303 –77.824 60.0 1.0 April 12
S-3 55.303 –77.824 60.0 5.0 April 12
S-3 55.303 –77.824 60.0 20.0 April 12
S-3 55.303 –77.824 60.0 50.0 April 12
S-4 55.430 –77.504 36.5 1.0 April 13
S-4 55.430 –77.504 36.5 10.0 April 13
S-4 55.430 –77.504 36.5 35.0 April 13
S-5 55.475 –77.595 31.0 1.0 April 13
S-5 55.475 –77.595 31.0 5.0 April 13
S-5 55.475 –77.595 31.0 30.0 April 13

NR = Nelson River; HR = Hayes River; W series = coastal waters near HR outflow; GWR = Great Whale River; S series = coastal waters near GWR outflow; N/A = not available.

Figure 3.

The two riverine coastal study areas on April 11, 2017, with locations of sampling sites. The southwestern Hudson Bay area (A) and southeastern Hudson Bay area (B) are shown by satellite imagery (Sentinel-1, Google Earth, and SIKU.org by Arctic Eider Society). In Panel A, the dark area is open water, which is surrounded by mobile pack ice, and the landfast ice is the band that lies parallel to shore with an irregular outer edge. The landfast ice is ridged near the shoreline because it goes up and down with the large tides and interacts with the seabed. In Panel B, the snow-covered land is light-toned and the landfast sea ice is darker. DOI: https://doi.org/10.1525/elementa.2020.00123.f3

Figure 3.

The two riverine coastal study areas on April 11, 2017, with locations of sampling sites. The southwestern Hudson Bay area (A) and southeastern Hudson Bay area (B) are shown by satellite imagery (Sentinel-1, Google Earth, and SIKU.org by Arctic Eider Society). In Panel A, the dark area is open water, which is surrounded by mobile pack ice, and the landfast ice is the band that lies parallel to shore with an irregular outer edge. The landfast ice is ridged near the shoreline because it goes up and down with the large tides and interacts with the seabed. In Panel B, the snow-covered land is light-toned and the landfast sea ice is darker. DOI: https://doi.org/10.1525/elementa.2020.00123.f3

Close modal

In SEHB, the river was completely ice-covered throughout its entire length in April 2017 and the Great Whale River estuary had a continuous cover of landfast sea ice (Figure 3B). Sampling sites S-1, S-2, and S-3 were close to the river mouth (Figure 3B), where previous work shows the river plume spreading out in a thin layer under the landfast ice (Ingram and Larouche, 1987). The water depth at these sampling sites ranged from 41 to 61 m (Table 1). The remaining SEHB sampling sites (S-4 and S-5) were located approximately 25 km away in a NE direction within 5 km of shore with water depths of 36.5 and 31 m, respectively (Table 1; Figure 3B).

2.3. Temperature and salinity measurements

Water samples were collected from depths of 1–6 m in SWHB and 1–50 m in SEHB using a Kemmerer sampler lowered through a hole cut in the landfast ice. The water was placed into an acid-cleaned amber Nalgene bottle and stored in a cooler to avoid freezing during transport by snowmobile back to the field lab. Subsequently, samples for salinity analysis were subsampled into precleaned glass bottles rinsed three times with sample water before filling. The in situ temperature and salinity measurements were performed using an Idronaut Ocean Seven CTD sensor with the temperature accuracy of ±0.002 °C, conductivity accuracy of 0.003 mS cm–1, and pressure (depth) accuracy of ±0.05% full scale.

Water samples were shipped to the University of Manitoba for analyses. Salinity was measured on bottle samples using a Guildline Autosal 8400 salinometer with a precision of 0.002 or greater. Samples were standardized against the International Association for the Physical Sciences of the Ocean Standard Seawater. The salinity sample from the Great Whale River was lost, so that in this work we assume a near-zero salinity (0.05) as previously reported for Great Whale River water near the river mouth (Granskog et al., 2011).

2.4. Incubation experiments

DOC biodegradability and lability were assessed using incubation experiments that were started directly after sample collection and followed an incubation protocol modified from Vonk et al. (2015). Within a few hours of field sampling, samples were filtered through precombusted (4 h, 450 °C) 0.7-µm glass fiber filters. The filtrate was decanted into duplicate incubation sets of eight precleaned glass vials (20 mL volume). In total, there were 54 incubation sets: 11 duplicate incubation sets for samples collected in SWHB and 16 duplicate incubation sets for samples collected in SEHB. About 2 ml of headspace was left in each vial in order to avoid oxygen depletion and allow CO2 outgassing during DOC degradation. The vials were tightly sealed and stored in the dark at 4 °C and shaken daily throughout the incubation period. The vials were acidified (200 µL of 2 M of hydrochloric acid (HCI) added) at different time points (T = 0, T = 0.3, T = 1, T = 3, T = 8, T = 15, T = 25, and T = 45 days) in order to terminate microbial activity and lower pH of the water to remove inorganic carbon. The time series was designed to capture the dynamics of both rapid DOC degradation in the first few days and slower DOC degradation over a span of a few weeks to few months.

All plastic equipment used in the experiment, such as syringes, forceps, filter holders, and vial caps, were rinsed with 10% HCl, while glass equipment, such as vials, were precombusted at 550 °C for 6–12 h. Special care was taken at all times to avoid organic carbon contamination.

2.5. Analysis of DOC content

DOC analysis was conducted using a Thermalox™ TOC-TN analyzer employing the thermal catalytic oxidation technique. Each incubation set was analyzed fully during the same DOC analysis session so that the same calibration curve could be applied to the whole set. Because DOC concentrations spanned a wide range between the river and marine samples, a new calibration curve specific for the expected DOC range of the samples to be analyzed was generated for each DOC analysis session. Deep seawater reference material from the Hansell Laboratory (Rosenstiel School of Marine and Atmospheric Science, University of Miami) with a DOC range of 0.492–0.528 mg L–1 was used as an additional calibration point when generating calibration curves for the low-DOC SEHB samples.

Five to eight injections were used for each sample depending on the coefficient of variation, which was considered valid when the coefficient was ≤2%. The lower the DOC concentration, the greater the number of injections that were needed to achieve a valid coefficient of variation. Coefficients of variation among duplicate samples were typically ≤5%.

The DOC degraded during the total length of the incubation experiment was defined as biodegradable DOC (BDOC). The results of the incubation experiment were reported primarily as %BDOC (fraction of total DOC that was biodegraded) at time point x according to the following equation:

$%BDOC=([DOCT= 0−DOCT=x]/DOCT= 0)×100,$
1

where DOCT = 0 is the initial DOC concentration, and DOCT = x is the DOC concentration at time x days.

The lability of BDOC describes how fast the degradation of DOC is happening. In accordance with Kirchman et al. (1993), the BDOC degraded by Day 3 was defined as labile (L-BDOC) and the BDOC degraded from Day 3 to Day 45 was defined as semilabile (SL-BDOC). Hence, %L-BDOC is the fraction of the total DOC degraded by Day 3.

Additionally, the results of the incubation experiment were reported with biodegradation rate constants, k′ and k″, which were computed from the following equations describing the microbial degradation rates of L-BDOC and SL-BDOC, respectively:

2.1
2.2

where DOCT = 0 is the initial DOC concentration, and DOCT = x is the DOC concentration at time x days.

3.1. Physical characteristics of the riverine and coastal waters

The temperatures of the riverine and coastal waters were close to their respective freezing points in mid-April 2017, except for subsurface waters below depths of 15–25 m in SEHB, which were slightly warmer (–0.8 °C; Figure 4). In SWHB, the shallow water column (<6 m) was well mixed (W-2) to weakly stratified (remaining stations) with surface salinity ranging from 1 to 12 and bottom salinity ranging from 1 to 17 (Figure 4A). Surface salinity in SWHB increased with increased distance from the river mouths (compare stations W-2 and W-3, located approximately 35 km northeast of the Hayes River mouth, to stations W-4 and W-5, located about 60 km northeast of the Hayes River mouth). However, Site W-1, closest to the Hayes River mouth, had a surface salinity of approximately 10, likely because the sample was collected at high tide and the discharge from the Hayes is relatively weak at this time of year. Overall, the range of salinity in SWHB samples (0.06–16.51) was much narrower than that in SEHB samples (0.11–29.3) with the latter group having mostly very low or very high salinities (Figure 4C). The vertical profiles of salinity and temperature in SEHB (Figure 4B) indicate a highly stratified under-ice plume, which is associated with the Great Whale River (Ingram and Larouche, 1987). Low-salinity (<4) river plume waters extend to a depth of 3–5 m and overlie cold, saline waters (25–29) at intermediate depths (Figure 4C). Warmer, slightly saltier waters underlie the cold, saline waters extending from depths of 20–30 m to the bottom (Figure 4B). A warm, salty subsurface water mass has been seen previously in southeast Hudson Bay, presumably a remnant from the previous summer (Petrusevich et al., 2018; Eastwood et al., 2020). Despite this regional feature, a strong linear relationship exists between temperature and salinity (p << .001, R2 = .86) across all the collected samples, which indicates mixing between fresh, cold (approximately 0 °C) river waters and a common salty (approximately 26), colder (approximately –1.3 °C) marine water mass (Figure 4C).

Figure 4.

Temperature and salinity depth profiles and the T-S relationship across all the samples (C). Panel A shows profiles for the southwestern Hudson Bay (SWHB) coastal stations, and Panel B shows the profiles for the southeastern Hudson Bay (SEHB) coastal stations. The thin lines on the temperature-depth profiles represent corresponding freezing temperatures at each station. The empty circles show the TS values of the DOC samples. On the T-S plot, the dashed linear regression line encompasses all samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f4

Figure 4.

Temperature and salinity depth profiles and the T-S relationship across all the samples (C). Panel A shows profiles for the southwestern Hudson Bay (SWHB) coastal stations, and Panel B shows the profiles for the southeastern Hudson Bay (SEHB) coastal stations. The thin lines on the temperature-depth profiles represent corresponding freezing temperatures at each station. The empty circles show the TS values of the DOC samples. On the T-S plot, the dashed linear regression line encompasses all samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f4

Close modal

3.2. In situ concentrations of DOC

In situ DOC concentrations in the Nelson and Hayes Rivers and in SWHB were higher than in the Great Whale River and SEHB (Figure 5). The Nelson, Hayes, and Great Whale Rivers had DOC concentrations of 23.6, 14.2, and 8.2 mg L–1, respectively. SWHB coastal samples had DOC concentrations ranging between 10.2 and 17.3 mg L–1, and SEHB coastal samples had DOC concentrations ranging between 1.9 and 6.9 mg L–1.

In SWHB, there was no statistically significant, linear relationship between DOC concentration and salinity, reflecting nonconservative behavior of riverine DOC across the Nelson/Hayes Estuary and/or influence of both high-DOC (Nelson) and low-DOC (Hayes) riverine sources (Figure 5). After a drop in DOC concentrations as salinity increased from 0 to 3, DOC concentrations increased slightly around a salinity of 10 (Figure 5). In SEHB, the data distribution appeared to be bimodal in terms of salinity, which may be attributed to the high stratification of the coastal waters. If, however, we assume a linear relationship between DOC and salinity (R2 = .54, but the assumption of data normality is not met), then DOC concentration decreased by 0.12 mg L–1 per unit of salinity increase (Figure 5).

Figure 5.

Relationships between in situ DOC concentrations and salinity in the two study areas. Data from the Nelson and Hayes Rivers pertain to the southwestern Hudson Bay (SWHB) coastal samples; data from the Great Whale River pertain to the southeastern Hudson Bay (SEHB) coastal samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00173.f5

Figure 5.

Relationships between in situ DOC concentrations and salinity in the two study areas. Data from the Nelson and Hayes Rivers pertain to the southwestern Hudson Bay (SWHB) coastal samples; data from the Great Whale River pertain to the southeastern Hudson Bay (SEHB) coastal samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00173.f5

Close modal

Table 2 summarizes the results of the microbial degradation experiments. The Nelson, Hayes, and Great Whale rivers contained 60%, 33%, and 24% BDOC, respectively, while BDOC in coastal waters was 16%–65% (averaging 23% with standard error of the mean, SEM, ± 4%). The k′ values for the Nelson, Hayes, and Great Whale rivers calculated according to Equation 2.1 were 0.236 d–1 (R2 = .75), 0.059 d–1 (R2 = .99), and 0.027 d–1 (R2 = .67), respectively, while the k” values (Equation 2.2) ranged between 0.003 and 0.005 d–1 (R2 > .55). The k′ values for the coastal ocean reached 0.262 d–1 with the mean = 0.049 (SEM ± 0.014; each R2 > .50), while the k” values for the coastal ocean reached 0.008 d–1 with the mean = 0.002 (SEM ± 0.001; each R2 > .55).

Table 2.

Results of the microbial degradation experiments performed for the two study areas: southwestern Hudson Bay and southeastern Hudson Bay. DOI: https://doi.org/10.1525/elementa.2020.00123.t2

StationDepth (m)SalinityTemperature (°C)DOC (mg L–1)BDOCa (%)L-BDOCa (%)k′ (d–1)k″ (d–1)
NR 0.15 23.6 60 53 0.236 0.003
HR 0.06 0.05 14.2 33 16 0.059 0.005
W-1 10.9 –0.61 16.6 25 15 0.038 0.003
W-1 11.4 –0.64 14.2 24 15 0.041 0.002
W-2 0.97 –0.03 12.3 16 0.023 0.002
W-3 2.94 –0.16 17.3 40 35 0.099 0.001
W-3 3.05 –0.16 10.2
W-4 8.46 –0.46 13.0 22 16 0.052 0.002
W-4 9.37 –0.52 13.4 26 0.022 0.004
W-5 11.76 –0.62 15.1 37 0.017 0.008
W-5 16.51 –0.86 12.5 42 20 0.061 0.007
GWR N/A 0.05 8.2 24 0.027 0.005
S-1 0.11 0.05 6.0 16 16 0.043
S-1 14.75 –0.59 2.4
S-1 10 26.28 –1.33 5.0 55 51 0.244 0.002
S-1 40 28.84 –0.91 1.9
S-2 0.72 –0.04 5.2
S-3 2.93 0.06 6.9 21 10 0.028 0.003
S-3 23.64 –1.10 3.2 16 0.026 0.002
S-3 20 26.89 –1.38 3.3 28 15 0.043 0.004
S-3 50 29.30 –0.73 2.0
S-4 1.35 –0.06 6.6 23 23 0.052
S-4 10 25.57 –1.22 5.7 59 35 0.133 0.008
S-4 35 27.71 –0.83 2.9
S-5 24.60 –1.33 2.5
S-5 25.61 –1.36 2.3
S-5 30 27.04 –1.34 5.4 65 57 0.262 0.004
StationDepth (m)SalinityTemperature (°C)DOC (mg L–1)BDOCa (%)L-BDOCa (%)k′ (d–1)k″ (d–1)
NR 0.15 23.6 60 53 0.236 0.003
HR 0.06 0.05 14.2 33 16 0.059 0.005
W-1 10.9 –0.61 16.6 25 15 0.038 0.003
W-1 11.4 –0.64 14.2 24 15 0.041 0.002
W-2 0.97 –0.03 12.3 16 0.023 0.002
W-3 2.94 –0.16 17.3 40 35 0.099 0.001
W-3 3.05 –0.16 10.2
W-4 8.46 –0.46 13.0 22 16 0.052 0.002
W-4 9.37 –0.52 13.4 26 0.022 0.004
W-5 11.76 –0.62 15.1 37 0.017 0.008
W-5 16.51 –0.86 12.5 42 20 0.061 0.007
GWR N/A 0.05 8.2 24 0.027 0.005
S-1 0.11 0.05 6.0 16 16 0.043
S-1 14.75 –0.59 2.4
S-1 10 26.28 –1.33 5.0 55 51 0.244 0.002
S-1 40 28.84 –0.91 1.9
S-2 0.72 –0.04 5.2
S-3 2.93 0.06 6.9 21 10 0.028 0.003
S-3 23.64 –1.10 3.2 16 0.026 0.002
S-3 20 26.89 –1.38 3.3 28 15 0.043 0.004
S-3 50 29.30 –0.73 2.0
S-4 1.35 –0.06 6.6 23 23 0.052
S-4 10 25.57 –1.22 5.7 59 35 0.133 0.008
S-4 35 27.71 –0.83 2.9
S-5 24.60 –1.33 2.5
S-5 25.61 –1.36 2.3
S-5 30 27.04 –1.34 5.4 65 57 0.262 0.004

NR = Nelson River; HR = Hayes River; W series = coastal waters near HR outflow; GWR = Great Whale River; S series = coastal waters near GWR outflow; N/A = not available; DOC = dissolved organic carbon.

aFraction of total DOC that is biodegraded (BDOC); L-BDOC is the fraction of BDOC degraded by Day 3.

The DOC concentrations measured at intervals over the incubation period (45 days) are shown in Figure 6. For both SWHB and SEHB samples, DOC concentrations fell within the first 3 days. For SWHB, DOC concentrations continued to decline more slowly throughout the first one-third (15 days) of the experiment (Figure 6). For SEHB, little change occurred between 3 and 15 days. Despite the large early decreases in SWHB DOC concentrations, they remained well above SEHB DOC concentrations throughout the incubation period.

Figure 6.

Change in DOC concentration with time during 45-day incubation experiments. The blue and red borders separate labile (L) and semilabile (SL) fractions of BDOC based on degradation time. Data from the Nelson and Hayes Rivers pertain to the southwestern Hudson Bay (SWHB) coastal samples; data from the Great Whale River pertain to the southeastern Hudson Bay (SEHB) coastal samples. DOC = dissolved organic carbon; BDOC = biodegradable DOC. DOI: https://doi.org/10.1525/elementa.2020.00123.f6

Figure 6.

Change in DOC concentration with time during 45-day incubation experiments. The blue and red borders separate labile (L) and semilabile (SL) fractions of BDOC based on degradation time. Data from the Nelson and Hayes Rivers pertain to the southwestern Hudson Bay (SWHB) coastal samples; data from the Great Whale River pertain to the southeastern Hudson Bay (SEHB) coastal samples. DOC = dissolved organic carbon; BDOC = biodegradable DOC. DOI: https://doi.org/10.1525/elementa.2020.00123.f6

Close modal

Figure 7 expresses the results for each river and estuary in terms of the proportion of DOC that was biodegradable overall (%BDOC) and the relative proportions of L-BDOC and SL-BDOC to that BDOC total. Approximately one-half of the BDOC in SWHB coastal waters and about three-quarters of the BDOC in SEHB coastal waters was degraded by Day 3 (by definition L-BDOC; Figure 7). The BDOC in Nelson River water was dominated by L-BDOC, whereas the Hayes River water and the SWHB coastal water had equivalent contributions of L-BDOC and SL-BDOC. Interestingly, the Great Whale River BDOC was mostly semilabile, but the SEHB coastal waters contained an even representation of L-BDOC and SL-BDOC. In other words, the relative contribution of L-BDOC to total BDOC increased from the Great Whale River water to the SEHB coastal waters (Figure 7). The relative contribution of L-BDOC in SEHB coastal waters was nearly on par with that of the Nelson River.

Figure 7.

Biodegradable fraction of total DOC with labile and semilabile proportions in riverine and coastal samples. Labile (L-BDOC) versus semilabile (SL-BDOC) proportions of the biodegradable fraction of total DOC (%BDOC) are shown for river samples and coastal samples from southwestern Hudson Bay (SWHB) and southeastern Hudson Bay (SEHB). Error bars indicate standard error of the mean (n = 9 for SWHB; n = 15 for SEHB). DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f7

Figure 7.

Biodegradable fraction of total DOC with labile and semilabile proportions in riverine and coastal samples. Labile (L-BDOC) versus semilabile (SL-BDOC) proportions of the biodegradable fraction of total DOC (%BDOC) are shown for river samples and coastal samples from southwestern Hudson Bay (SWHB) and southeastern Hudson Bay (SEHB). Error bars indicate standard error of the mean (n = 9 for SWHB; n = 15 for SEHB). DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f7

Close modal

The %BDOC in relation to the initial DOC concentration and salinity is presented in Figure 8. Overall, %BDOC increased with DOC concentration, but at different rates depending on the location and salinity. In SEHB, %BDOC and DOC concentrations were related linearly to the high salinity (14.6–29.2) and low DOC waters (p < .001, R2 = .94), and the increase in %BDOC with increasing DOC is so steep that %BDOC reaches >50% for DOC concentrations ≥5 mg L–1 (Figure 8). In low-salinity SEHB waters (0.05–2.93), the %BDOC-DOC relationship has a very similar slope to that of the high salinity waters at low DOC concentrations but levels off at DOC concentrations above 7 mg L–1 and BDOC of approximately 23%. In SWHB, which has only low-salinity waters, %BDOC increases with initial DOC concentrations but at a much slower rate than in SEHB. The %BDOC-DOC relationship in SWHB waters seems to level off when %BDOC reaches 60% (Figure 8).

Figure 8.

Relationships between biodegradable DOC and initial DOC concentration and color visualized with salinity. The solid black line represents a significant linear trend for the high-salinity, southeastern Hudson Bay (SEHB) samples, while the dashed lines represent apparent linear trends that level off at high DOC values for the low-salinity SEHB and southwestern Hudson Bay (SWHB) samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f8

Figure 8.

Relationships between biodegradable DOC and initial DOC concentration and color visualized with salinity. The solid black line represents a significant linear trend for the high-salinity, southeastern Hudson Bay (SEHB) samples, while the dashed lines represent apparent linear trends that level off at high DOC values for the low-salinity SEHB and southwestern Hudson Bay (SWHB) samples. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f8

Close modal

4.1. DOC concentrations

4.1.1. DOC concentrations in rivers

The Nelson and Hayes Rivers had higher DOC concentrations than the Great Whale River in late winter, consistent with previous reports that the Nelson and Hayes Rivers on average contain two to six times greater DOC concentration than other Hudson Bay rivers (Mundy et al., 2010). The boreal forest, wetlands/peatlands, and thawing permafrost covering the Nelson and Hayes River drainage basins (Smith et al., 2015; Godin et al., 2017) are expected to leach a large amount of DOC to the rivers (Benner et al., 2004; Rontani et al., 2014). The much higher DOC of the Nelson River (23.6 mg L1) versus Hayes River (14.2 mg L1) in April was unexpected, given that similar concentrations are found in the two rivers during summer (Granskog et al., 2011; Godin et al., 2017). However, DOC concentrations in the Nelson have been noted as remaining high during winter (Rosa et al., 2012), perhaps due to sustained supplies from lakes, reservoirs, and wetlands throughout the winter period when soils are frozen. To our knowledge, there are no previous winter DOC data for the Hayes.

Little is known about the sources of DOC in the Nelson River. Kirk and St. Louis (2009) found wide variation in DOC concentrations (15.1 ± 19.7 mg L–1, n = 66) when they conducted biweekly sampling in the northern part of the Nelson River at the Limestone Generating Station (120 km from the river mouth). Particularly high concentrations (>20 mg L–1) occurred irregularly but always during periods of relatively high flow (>3,000 m3 s–1; see Figure S1 in Kirk and St. Louis, 2009). Manitoba Hydro’s monitoring data associated with the Keeyask development (https://keeyask.com/wp-content/uploads/2020/06/PEMP-2020-01-Physical-Environment-Monitoring.pdf) indicate a typical DOC concentration of 8–10 mg L–1 for the open-water season. However, DOC concentrations in the Nelson River also vary from place to place. Data (n = 851 samples) collected throughout the Nelson River watershed as part of the Coordinated Aquatic Monitoring Program (2014) show clustering of high DOC concentrations in the northern part of the basin (Figure 9). Small tributaries and backbays monitored during the open-water period have high concentrations of DOC (>20 mg L–1; T Papakyriakou, personal communication, 2020).

Figure 9.

DOC concentrations in samples along a latitudinal gradient in the Nelson River watershed. Data are from the Coordinated Aquatic Monitoring Program 2004–2013 (n = 850; one high outlier excluded near 51.68°N), spanning all seasons. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f9

Figure 9.

DOC concentrations in samples along a latitudinal gradient in the Nelson River watershed. Data are from the Coordinated Aquatic Monitoring Program 2004–2013 (n = 850; one high outlier excluded near 51.68°N), spanning all seasons. DOC = dissolved organic carbon. DOI: https://doi.org/10.1525/elementa.2020.00123.f9

Close modal

In comparison to other major northern rivers (Holmes et al., 2018), the Nelson clearly emerges as a DOC-rich river (Figure 10). DOC concentrations measured across this and previous studies (Kirk and St. Louis, 2009; Rosa et al., 2012; Godin et al., 2017) place the Nelson River at the upper end of the range for major northern rivers, particularly during winter (Figure 10). Presumably, DOC supply from the Nelson watershed is stronger especially during winter due to its sub-Arctic setting and more extensive wetlands. With these high DOC concentrations and relatively high flows maintained during winter by regulation, the Nelson River winter DOC fluxes would likely exceed those of the Mackenzie, Yenisei, and Ob Rivers despite their larger size.

Figure 10.

DOC concentrations in the major Arctic rivers compared to those measured in Hudson Bay rivers. Mean monthly DOC concentrations for the major Arctic rivers (PARTNERS rivers; Holmes et al., 2018) and April DOC concentrations for the Nelson, Hayes, and Great Whale Rivers of Hudson Bay (A), with a more detailed comparison of April discharge versus DOC concentrations for the same set of rivers (B). Error bars represent standard error of the mean; n values for means in Panels A and B are provided in Table S1 and Table S2, respectively. DOI: https://doi.org/10.1525/elementa.2020.00123.f10

Figure 10.

DOC concentrations in the major Arctic rivers compared to those measured in Hudson Bay rivers. Mean monthly DOC concentrations for the major Arctic rivers (PARTNERS rivers; Holmes et al., 2018) and April DOC concentrations for the Nelson, Hayes, and Great Whale Rivers of Hudson Bay (A), with a more detailed comparison of April discharge versus DOC concentrations for the same set of rivers (B). Error bars represent standard error of the mean; n values for means in Panels A and B are provided in Table S1 and Table S2, respectively. DOI: https://doi.org/10.1525/elementa.2020.00123.f10

Close modal

4.1.2. DOC concentrations in the coastal waters

The high DOC flux from the Nelson River in April (6,730 t day–1, assuming 3,300 m3 s–1 and DOC concentration of 23.6 mg L–1), together with the flux from the Hayes (440 t day–1, assuming 360 m3 s–1 and 14.2 mg L–1, respectively), produces a low-salinity, high-DOC estuarine environment in SWHB coastal waters (Figure 5). These conditions are not unusual in the Arctic, where high concentrations of DOC in surface coastal waters are sourced from river discharge containing higher concentrations of terrigenous DOC (Amon and Benner, 2003), but they are rarely described in winter. In contrast, SEHB coastal waters are characterized in April by very fresh (approximately 0) river plume waters overlying salty (25–30) waters with DOC concentrations nowhere exceeding 7 mg L–1 (Table 1; Figure 5). In both environments, DOC concentrations show a tendency to decrease with increasing salinity (Figure 5), as expected from the low DOC concentrations (1.3 mg L–1) in marine waters of Hudson Bay (Mundy et al., 2010). Conservative mixing of riverine and marine water masses would yield a linear salinity–DOC relationship with a steeper slope for the SWHB coastal area where DOC concentrations in river water were higher (cf. Hansell et al., 2004). The nonlinear salinity–DOC relationship that we observed in SWHB (Figure 5) thus provides a first indication that the DOC supplied by the Nelson and Hayes Rivers is more labile and/or degradable than that supplied by the Great Whale River.

4.2.1. %BDOC and lability in rivers

Comparison of %BDOC values for the three Hudson Bay rivers and other Arctic rivers (Table 3) indicates relatively high biodegradability of Hudson Bay riverine DOC in late winter. The values in the table vary widely but in general fall within the range of 10%–40% BDOC after 1–3 months of incubation; the available data are biased toward the summer season and incubation temperatures of 15 °C and 20 °C. In view of the low incubation temperature, we used 4 °C, which might be expected to limit microbial degradation, the %BDOC results for the Hayes and Great Whale Rivers are notable for their comparability to those reported for other Arctic rivers. The 60% BDOC of the Nelson River is exceptional (Table 3). We speculate that continued DOC supply during winter (e.g., from reservoirs, wetlands/peatlands, and thawing permafrost) together with better preservation during transport compared to open-water conditions account for the high biodegradability of the Hudson Bay riverine DOC and the Nelson River DOC, specifically.

Table 3.

Summary of recent incubation experiments examining microbial degradation of DOC in Arctic rivers. DOI: https://doi.org/10.1525/elementa.2020.00123.t3

ReferenceStudy SiteSampling TimeInitial DOC (mg L–1)Incubation Time (Days)Incubation Temperature (°C)Parameter Analyzed%BDOCa
Kawahigashi et al. (2004)  Yenisei River August 4–30 97 20 CO2 4–28
Holmes et al. (2008)  Alaskan rivers Spring freshet 2–15 90 20 DOC 14–33
July 2–8 90 20 DOC Up to 9
Balcarczyk et al. (2009)  Alaska streams June 2–7 40 DOC 4–17
20 DOC 4–35
Mann et al. (2012)  Kolyma River Preflush (May) 2–3 28 20 DOC Up to 10
Freshet (May to June) 8–14 28 20 DOC 1–20
Wickland et al. (2012)  Yukon River Winter 28 CO2 27–35
Spring freshet 3–17 28 CO2 2–29
Summer 3–7 28 CO2 4–19
Yukon River and its tributaries Winter 2–4 28 15 CO2 9–53
Spring freshet 2–17 28 15 CO2 5–47
Summer 3–14 28 15 CO2 3–36
Herlemann et al. (2014)  Kalix River June 28 10 DOC 6–8
Shirokova et al. (2017)  Severnaya Dvina River June (spring flood) 10–15 Approximately 25 20 DOC Negligible
August 12–15 Approximately 25 20 DOC 15–20
This study Nelson River April 24 45 DOC 60
Hayes River April 14 45 DOC 33
Great Whale River April 45 DOC 24
ReferenceStudy SiteSampling TimeInitial DOC (mg L–1)Incubation Time (Days)Incubation Temperature (°C)Parameter Analyzed%BDOCa
Kawahigashi et al. (2004)  Yenisei River August 4–30 97 20 CO2 4–28
Holmes et al. (2008)  Alaskan rivers Spring freshet 2–15 90 20 DOC 14–33
July 2–8 90 20 DOC Up to 9
Balcarczyk et al. (2009)  Alaska streams June 2–7 40 DOC 4–17
20 DOC 4–35
Mann et al. (2012)  Kolyma River Preflush (May) 2–3 28 20 DOC Up to 10
Freshet (May to June) 8–14 28 20 DOC 1–20
Wickland et al. (2012)  Yukon River Winter 28 CO2 27–35
Spring freshet 3–17 28 CO2 2–29
Summer 3–7 28 CO2 4–19
Yukon River and its tributaries Winter 2–4 28 15 CO2 9–53
Spring freshet 2–17 28 15 CO2 5–47
Summer 3–14 28 15 CO2 3–36
Herlemann et al. (2014)  Kalix River June 28 10 DOC 6–8
Shirokova et al. (2017)  Severnaya Dvina River June (spring flood) 10–15 Approximately 25 20 DOC Negligible
August 12–15 Approximately 25 20 DOC 15–20
This study Nelson River April 24 45 DOC 60
Hayes River April 14 45 DOC 33
Great Whale River April 45 DOC 24

DOC = dissolved organic carbon.

aPercentage of total DOC that is biodegradable.

Permafrost is extensive near the mouth of the Nelson and Hayes Rivers, and its degradation is associated with increasing river bank instability (Stainton, 2019; K Skaftfeld, personal communication, 2020). Permafrost DOC is not only biodegradable but also degraded relatively rapidly (approximately 50% within 1 week) once released into streams (Drake et al., 2015; Mann et al., 2015; Spencer et al., 2015; Ward and Cory, 2015; Müller et al., 2018). Both the Nelson and Hayes had relatively large proportions of L-BDOC (DOC degraded within 3 days; Figure 7), consistent with expectations for permafrost DOC. We do not know the transport time from DOC source to the river mouths where sampling occurred, but faster transport in the Nelson River with its high (regulated) winter flows perhaps contributed to its higher BDOC and L-BDOC compared to the Hayes.

Another potential source of degradable DOC along the Nelson is the series of lakes and reservoirs formed by damming, which support primary production during summer yet have limited storage relative to the larger and deeper reservoirs of, for example, La Grande Rivière system. The possible transport of labile aquatic DOC from the reservoirs during winter might explain the measured L-BDOC (Figure 7). Water residence times in the Nelson reservoirs are too short for DOC buildup in deep waters, but labile DOC could be sourced from reservoir sediments. Higher preservation of riverine DOC during winter versus summer transport certainly may be expected because of the darkness (under the snow and ice cover), which limits photodegradation, and low temperatures, which limit microbial degradation. The small areas of open water below each dam may have allowed some minor photobleaching that, while not affecting DOC concentrations, facilitated photo-enhanced biomineralization of DOC (Fichot and Benner, 2014). Guéguen et al. (2016) found that during summer, photobleaching rates for Nelson River dissolved organic matter (DOM; i.e., decrease in DOM signal resulting from exposure to solar radiation) ranged from 0.005 to 0.030 h–1, corresponding to half-lives of 4.9–9.9 days. However, there was so little open water during April along the lower Nelson River that we may assume photobleaching had negligible effects on our results, at least until the riverine DOC arrived at the Nelson/Hayes Estuary and may have been recirculated in the area of the flaw lead (Figure 3A).

Previous work suggests that %BDOC in northern rivers is strongly related to nutrient availability (cf. Wickland et al. 2012). Indeed, Wickland et al. (2012) found that the molar ratio between DOC and dissolved inorganic nitrogen (DIN; typically NO3 + NH4) was a good predictor of %BDOC in the Yukon River and tributaries and speculated that higher %BDOC in winter arose from greater nutrient availability at that time. The Nelson River does appear to be a nitrate-rich river relative to other Hudson Bay rivers (see Table 1 in Kuzyk et al., 2010). Lacking ammonium data for our samples, our best estimates for the DIN concentrations in late winter are previously collected Coordinated Aquatic Monitoring Program data for the Nelson and Hayes Rivers in February to March (8.8 and 6.3 µM, respectively) and previously published April DIN concentrations for the Great Whale River (4 µM; Hudon et al., 1996). Substituting these values with our DOC data into the Wickland et al.’s (2012) equation yields predictions of %BDOC no more than 10% for any of the samples. The Wickland et al.’s (2012) equation may underestimate the late winter %BDOC in the Hudson Bay rivers, particularly the Nelson, because we are underestimating DIN or the differences in DOM composition or microbial community among rivers may significantly alter the relationship.

4.2.2. %BDOC and lability in the coastal waters

A marked decline in %BDOC, and especially the proportion of labile DOC, was observed between the Nelson River and SWHB coastal waters (Figure 7), which we attribute to the dynamics of the Nelson/Hayes Estuary and specifically the accumulation of older, semilabile DOC along the southern shore where winter sampling was conducted. Manitoba Hydro’s MIKE model simulates river water residence times of 30–40 days for that portion of the estuary (K Sydor, Manitoba Hydro, personal communication, 2019). Presumably, the long residence times mean that the in situ DOC was not all freshly introduced riverine DOC, but a mixture including older residual, much less labile DOC. If the entire estuary had been covered with landfast ice allowing sampling in the central (thalweg) area, where the residence time is much shorter, presumably we would have seen higher biodegradability and lability of the DOC. During summer, Guéguen et al. (2016) found a strong, negative linear correlation between salinity and DOC (i.e., conservative behavior) throughout the estuary even in the thalweg area. However, the low lability of DOC in the coastal waters in summer might simply reflect poor preservation of labile Nelson River DOC during transport to the estuary. During open-water conditions, DOC may be degraded upstream, closer to its point of release, whereas during the cold, dark season, the site of degradation shifts downstream to the river mouth and inner estuary.

In contrast to SWHB, SEHB coastal waters had, on average, similar %BDOC to the Great Whale River and, surprisingly, a higher proportion of L-BDOC (Figure 7). Upon closer inspection, and by separating the very fresh surface (river plume) waters from the saltier coastal waters in Figure 8, the higher-salinity SEHB coastal waters (and not the very fresh plume waters) can be seen to have the higher %BDOC and %L-BDOC than the river (Figure 8). The high %BDOC (>60% in some cases) occurred despite low DOC concentrations in the salty SEHB coastal waters (Figure 8). A possible explanation is the priming effect or co-metabolism (Bianchi, 2011; Ward et al., 2016) accompanied by increased nutrient (dissolved inorganic nitrogen or phosphorus) availability in the more marine-dominated waters (cf. McClelland et al., 2007; Holmes et al., 2008). Marine organic matter (e.g., algae) produced in the coastal ocean is much more labile than the terrestrial organic matter supplied by rivers (Bauer et al., 2013), so that its presence is believed to improve the efficiency with which refractory tDOC is degraded (also known as priming effect or co-metabolism). Sea-ice microalgae begin growing in southeast Hudson Bay around April 10 (Gosselin et al., 1985), possibly providing a marine organic matter source. However, the priming effect cannot be confirmed or rejected with the available data. The DOC concentrations were very low and the sources are not known. Remnant old(er) marine DOC or disintegrated POC produced in the previous summer may account for some or all of the observed %BDOC.

In terms of nutrients, very low nitrate concentrations are characteristic of Hudson Bay surface waters, while river waters are notoriously low in phosphate. Thus, %BDOC may have increased in the more saline estuary waters due to better nutrient availability, possibly resulting from the release of adsorbed NH4+ from riverine sediments as they encountered the saline conditions of the estuary (cf. Rysgaard et al., 1999) combined with seawater supply of phosphate, which has 5–6-fold higher phosphate concentrations than river water (Z Kuzyk, personal communication, 2020). The %BDOC in the very fresh plume waters of the SEHB coast increased in direct relation to concentrations of DOC, but the apparent “leveling off” at DOC concentrations above 7 mg L–1 (Figure 8) is consistent with the notion of nutrient limitation in that water mass. The under-ice plume of the Great Whale River apparently also retains river water (and associated DOC) for months at a time (Ingram and Larouche, 1987), which may mean that the plume becomes increasingly phosphate-poor over time as river-derived phosphate is consumed and little crosses the sharp halocline from the deep marine waters. Thus, we speculate that the river plume waters lacked the nutrients for efficient biodegradation of DOC, whereas the more saline coastal surface waters provided a more favorable environment for degradation.

Interestingly, the only saline SEHB samples showing no biodegradation were those drawn from the warmer water mass below 30-m depth (Table 2). The cold saline waters found between depths of 0 and 20 m at stations S1–S4 were the waters that showed high %BDOC (16%–59%; Table 2). Presumably, the warm water mass, which has been isolated from the surface and thus riverine DOC inputs since the previous summer (Petrusevich et al., 2018; Eastwood et al., 2020), contained only residual recalcitrant DOC.

Biodegradation rate constants complement %BDOC results as they provide additional insight regarding biodegradation dynamics. However, methodological differences in how L-BDOC and SL-BDOC are defined and in curve-fitting approaches make directly comparing results difficult in most cases. For example, Kawahigashi et al. (2004) reported k′ values for the Yenisei River (1.833–2.514 d–1), nearly an order of magnitude greater than the k′ values for Hudson Bay rivers (0.027–0.236 d–1) and k” values that were much lower (<0.001 d–1) than the k” values in our study (0.003–0.005 d–1). These differences arise from Kawahigashi et al. (2004) having defined the degradation of L-BDOC as the first 15–90 min of the incubation period (rather than 3 days), having used a model that forced the best fit of both L-BDOC and SL-BDOC simultaneously, and having a longer total incubation time. These sorts of methodological differences are what motivated Vonk et al. (2015) to propose a standardized BDOC protocol, which we adopted with few modifications for the present work.

A potential limitation in our incubation protocol due to the remote setting and lack of facilities at our field stations is that samples were shaken only once per day rather than constantly agitated (e.g., in a rotating chamber). Ward et al. (2019) suggested that DOC biodegradation rate in rivers may depend on water mixing as a result of interactions among suspended particles, dissolved constituents, and free-living and particle-bound microbes. Had we been able to implement constant mixing of our river waters, we expect it would have increased the biodegradation rates observed.

Compared to rivers, coastal ocean biodegradation rate constants have been reported much more frequently and compiled in some recent review articles (e.g., Lønborg and Álvarez-Salgado, 2012; Kaiser et al., 2017). DOC biodegradation rate constants for the Hudson Bay coastal waters (k′ = 0.049 ± 0.014 d–1, n = 24; mean k” = 0.002 ± 0.001 d–1, n = 24) are comparable to the general biodegradation rate constants (k combining k′ and k", incubation time > 40 days) reported for the global coastal ocean (0.066 ± 0.065 d–1, n = 127; Lønborg and Álvarez-Salgado, 2012). This agreement is not surprising as the coastal ocean DOC is supplied mostly by rivers. Moreover, if we apply the formula for the correlation of tDOC decay constants with incubation time (n = 8) derived in the review study of Kaiser et al. (2017), the resulting 3-day biodegradation rate (k′) of 0.053 ± 0.024 d–1 and 45-day biodegradation rate (k combining k′ and k") of 0.008 ± 0.002 d–1 are very similar to our findings for the Hudson Bay coastal waters.

4.2.4. Implications of DOC degradation for carbon cycling in Hudson Bay coastal waters

For more than a decade, coastal waters of Hudson Bay have been recognized as supersaturated in CO2 with respect to the atmosphere during late summer, while offshore waters are undersaturated, and the distribution of CO2 has been recognized as strongly related to river discharge (Else et al., 2008). A recent carbon budget integrating inorganic and organic forms of carbon demonstrated that on an annual basis, the remineralization of terrestrial DOC is the main driver of CO2 accumulation in the bay’s coastal surface waters (Capelle et al., 2020). During the month of July, Burt et al. (2016) observed elevated dissolved inorganic carbon (DIC) concentrations in the Nelson/Hayes Estuary (approximately 100 uM excess DIC; see Figure 6B in Burt et al., 2016, and associated discussion) as compared to the remaining coastal waters of Hudson Bay, which implied that biodegradation of the DOC delivered by these rivers plays a disproportionate role in the CO2 patterns.

Our biodegradation results may be used to provide a first test of the hypothesis of Burt et al. (2016) that the Nelson/Hayes Estuary contains approximately 100 µM excess DIC relative to the remaining riverine coastal domain of Hudson Bay. Our DOC concentrations and %BDOC values for the Nelson and Hayes Rivers (Table 2), in combination with the April river discharge data (Figure 2), imply that 1,105 µM DIC is produced as a result of biodegradation of the riverine DOC: [(Riverine DOC × %BDOC × Discharge)Nelson + (Riverine DOC × %BDOC × Discharge)Hayes] / (Combined discharge of the Nelson and Hayes Rivers) = [(1,963 µM × 0.600 × 3,646 m3 s–1) + (1,183 µM × 0.325 × 366 m3 s–1)] / 4,012 m3 s–1 = 1,105 µM of DIC. Assuming a 40-day residence time of river water in the Nelson/Hayes Estuary (based on Manitoba Hydro’s MIKE model), the addition of river water with the 1,105 µM DIC excess concentration would result in about 1.5 × 1016 µmol DIC added to the estuary in excess, yielding a DIC excess concentration of 136 µmol L–1 in the estuary (assuming a total volume of water in the estuary is 1.13 × 1011 m3 based on 4,500 km2 with an average depth of 25 m). Thus, the DIC excess of approximately 100 µM observed by Burt et al. (2016) in the Nelson/Hayes Estuary is reasonable or even low in comparison to what appears to be the potential DIC enrichment in the Nelson/Hayes Estuary resulting from the biodegradation of riverine DOC. A 30-day residence time would yield a DIC excess concentration of 102 µM, in near-perfect agreement with the observations of Burt et al. (2016). Also likely is that some of the excess DIC produced by biodegradation during April and the subsequent ice-covered months (May and part of June) would be lost from the estuary by CO2 evasion prior to an observation period in July. The only other estimate of Nelson/Hayes freshwater residence time that we can find is the one given by Granskog et al. (2009), which was in the order of months. Our biodegradation results imply that a DIC buildup much stronger than 100 µM would be expected if this longer residence time is accurate.

The Arctic Ocean receives the greatest riverine DOC input in the global ocean on a per volume basis, and much progress has been made in recent years on quantifying this input (Holmes et al., 2011; Tank et al., 2012). The purpose of this study was to contribute to the scarce but accumulating knowledge of DOC degradation in Arctic rivers and the Arctic riverine coastal domain in winter, which has been mostly ignored because of the assumption that low base flow and low DOC concentration result in insignificant DOC supply. However, recent observations show an increase in winter freshwater discharge in the Arctic (Déry et al., 2016; St. Jacques et al., 2009; Walvoord and Striegl, 2007), and this study reports relatively high DOC concentrations for the Nelson, Hayes, and Great Whale rivers and adjacent Hudson Bay coastal waters during winter.

The winter DOC concentration in the Nelson River is high relative to other major Arctic rivers, probably due to relatively high water levels mobilizing DOC from permafrost in river banks or from wetlands, lakes, or reservoirs. The Nelson River winter DOC is also highly biodegradable. The high biodegradability and lability of the Nelson River DOC may relate to high contributions of fresh permafrost DOC leaching from sources along the lower Nelson River and its tributaries and/or DOC production in the system’s reservoirs. The composition of Nelson River DOC warrants further study.

The receiving coastal waters of SWHB have lower %BDOC than the source rivers due to DOC degradation. In each study area, %BDOC increased with DOC concentration, which is likely due to heterotrophic microbial communities performing better at a higher concentration of organic matter that includes more biodegradable components. Likely, in the estuarine settings, there were differences in the extent of DOC degradation that had occurred prior to the sampling and lowered both %BDOC and DOC concentration. The long residence time of the river water in the Nelson/Hayes Estuary (ca. 40 days) causes an accumulation of less degradable DOC in the estuary, which, in turn, leads to a decrease in DOC lability and a decrease in the slope of the %BDOC-DOC relationship in the SWHB coastal waters (Figure 7). In the SEHB coastal waters, the lability of the DOC is higher in the salty waters than in the source river, presumably due to the priming effect and higher phosphate concentrations in the coastal ocean.

We acknowledge that the explanations for the biodegradability and lability dynamics described above are mostly hypothetical and need to be verified by analyzing heterotrophic bacterial communities and assessing DOC chemical structure/properties. Thus, future studies need to complement biodegradation experiments with such measurements as well as expand on the types of rivers and seasons studied. Models of organic carbon budget, water–atmosphere CO2 exchange, and food web dynamics in the Arctic watersheds and the Arctic Ocean would greatly benefit from this research.

Data sets used within this study have been archived on the Canadian Watershed Information Network (CanWIN; http://lwbin-datahub.ad.umanitoba.ca). River discharge data for the Nelson, Hayes, and Great Whale rivers can be accessed at the website of Environment and Natural Resources, Government of Canada (https://wateroffice.ec.gc.ca). DOC concentration and discharge data for the major Arctic rivers (PARTNERS) can be accessed at the website of the Arctic Great Rivers Observatory (https://arcticgreatrivers.org/data/).

Table S1. Number of samplesa (n) used to produce Figure 10A.

Table S2. Number of samplesa (n) used to produce Figure 10B.

The authors thank the Arctic Eider Society’s Community-Driven Research Network and the Hunters and Trappers Organization of Kuujjuarapik for support of the fieldwork in SE Hudson Bay and Churchill Wild for facilitating fieldwork in SW Hudson Bay. They thank all the collaborators of the BaySys Program who assisted in the field (D. Babb, J. Liu, K. Munson, L. Dalman, N. Theriault, and V. Petrusevich) or supported lab analyses (D. Armstrong and Dr F. Wang). They also thank K. Wong and K. Sydor from Manitoba Hydro for providing information on residence times for the Nelson Estuary. J. Kirk kindly shared dissolved organic carbon data for the Nelson River from her past work; Manitoba Hydro and Manitoba partners kindly shared the Coordinated Aquatic Monitoring Program data. P. Kimber assisted with figures. This study was much improved by comments from Dr C. Guéguen (Université de Sherbrooke, Sherbrooke, QC), B. Gill (Manitoba Hydro), and Dr K. Kaiser (Texas A&M University, Galveston, TX) as well as from two anonymous reviewers and the Editor-in-Chief Dr J. Deming.

This study is a contribution to the Arctic Science Partnership (ASP, asp-net.org) and ArcticNet. Funding was provided by the Natural Sciences and Engineering Council of Canada (NSERC) Collaborative Research and Development project: BaySys (CRDPJ 470028-14) led by Dr D. Barber (Academic PI) and K. Sydor (Manitoba Hydro, Industry PI) and by the NCE ArcticNet (Hudson Bay Freshwater-Marine Coupling project; Dr Z. Kuzyk PI). Individual support from NSERC was provided to Dr Z. A. Kuzyk, Dr T. Papakyriakou, Dr D. W. Capelle, and Dr S. Rysgaard.

The authors declare that they have no conflict of interest.

Contributed to the conception and design: ZVK, ZAK, TP.

Contributed to the acquisition of data: ZVK, CMK.

Contributed to the analysis and interpretation of data: ZVK, ZAK, DWC, SR, TP.

Drafted the article or revised it critically for important intellectual content: ZVK, ZAK, DWC, CMK, SR, TP.

Approved the submitted version for publication: ZVK, DWC, CMK, SR, TP, ZAK.

Amon
,
RMW
,
Benner
,
R
.
2003
.
Combined neutral sugars as indicators of the diagenetic state of dissolved organic matter in the Arctic Ocean
.
Deep Sea Research I
50
:
151
169
.
Amon
,
RMW
,
Meon
,
B
.
2004
.
The biogeochemistry of dissolved organic matter and nutrients in two large Arctic estuaries and potential implications for our understanding of the Arctic Ocean system
.
Marine Chemistry
92
:
311
330
.
Amon
,
RMW
,
Rinehart
,
AJ
,
Duan
,
S
,
Louchouarn
,
P
,
Prokushkin
,
A
,
Guggenberger
,
G
,
Bauch
,
D
,
Stedmon
,
C
,
Raymond
,
PA
,
Holmes
,
RM
,
McClelland
,
JW
,
Peterson
,
BJ
,
Walker
,
SA
,
Zhulidov
,
AV
.
2012
.
Dissolved organic matter sources in large Arctic rivers
.
Geochimica et Cosmochimica Acta
94
:
217
237
.
Balcarczyk
,
KL
,
Jones
Jr.,
JB
,
Jaffe´
,
R
,
Maie
,
N.
2009
.
Stream dissolved organic matter bioavailability and composition in watersheds underlain with discontinuous permafrost
.
Biogeochemistry
94
:
255
270
.
Bauer
,
JE
,
Cai
,
WJ
,
Raymond
,
PA
,
Bianchi
,
TS
,
Hopkinson
,
CS
,
Regnier
,
PA
.
2013
.
The changing carbon cycle of the coastal ocean
.
Nature
504
:
61
70
.
Bélanger
,
S
,
Xie
,
H
,
Krotkov
,
N
,
Larouche
,
P
,
Vincent
,
WF
,
Babin
,
M
.
2006
.
Photomineralization of terrigenous dissolved organic matter in Arctic coastal waters from 1979 to 2003: Interannual variability and implications of climate change
.
Global Biogeochemical Cycles
20
:
GB4005
.
Bendtsen
,
J
,
Hilligsøe
,
KM
,
Hansen
,
JLS
,
Richardson
,
K
.
2015
.
Analysis of remineralisation, lability, temperature sensitivity and structural composition of organic matter from the upper ocean
.
Progress in Oceanography
130
:
125
145
.
Benner
,
R
,
Benitez-Nelson
,
B
,
Kaiser
,
K
,
Amon
,
RMW
.
2004
.
Export of young terrigenous dissolved organic carbon from rivers to the Arctic Ocean
.
Geophysical Research Letters
31
(
5
):
L05305
. DOI: http://dx.doi.org/10.1029/2003GL019251.
Benner
,
R
,
Louchouarn
,
P
,
Amon
,
RMW
.
2005
.
Terrigenous dissolved organic matter in the Arctic Ocean and its transport to surface and deep waters of the North Atlantic
.
Global Biochemical Cycles
19
:
GB2025
.
Bianchi
,
TS
.
2011
.
The role of terrestrially derived organic carbon in the coastal ocean: A changing paradigm and the priming effect
.
Proceedings of the National Academy of Sciences of the United States of America
.
108
(
49
):
19473
19481
.
Burt
,
WJ
,
Miller
,
LA
,
Granskog
,
MA
,
Papakyriakou
,
TN
.
2016
.
Inorganic carbon cycling and biogeochemical processes in an Arctic inland sea (Hudson Bay)
.
Biogeosciences
13
(
16
):
4659
4671
.
Capelle
,
DW
,
Kuzyk
,
ZA
,
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
:
102319
.
Carmack
,
E
,
Winsor
,
P
,
Williams
,
W
.
2015
.
The contiguous Panarctic Riverine Coastal Domain: A unifying concept
.
Progress in Oceanography
139
:
13
23
.
Cauwet
,
C
,
Sidorov
,
I
.
1996
.
The biogeochemistry of Lena River: Organic carbon and nutrients distribution
.
Marine Chemistry
53
:
211
217
.
Coordinated Aquatic Monitoring Program
.
2014
.
Three year summary report (2008–2010)
.
Report prepared for the Manitoba/Manitoba Hydro MOU working group
.
:
North/South Consultants Inc
..
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
:
341
351
.
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
:
4801
4818
.
Drake
,
TW
,
Wickland
,
KP
,
Spencer
,
RGM
,
McKnight
,
DM
,
Striegl
,
RG
.
2015
.
Ancient low–molecular-weight organic acids in permafrost fuel rapid carbon dioxide production upon thaw
.
Environmental Sciences
112
:
13946
13951
.
Eastwood
,
RA
,
Macdonald
,
RW
,
Ehn
,
JK
,
Heath
,
J
,
Arragutainaq
,
L
,
Myers
,
PG
,
Barber
,
DG
,
Kuzyk
,
ZA
.
2020
.
Role of river runoff and sea ice brine rejection in controlling stratification throughout winter in southeast Hudson Bay
.
Estuaries and Coasts
43
:
756
786
.
Else
,
BGT
,
Papakyriakou
,
TN
,
Granskog
,
MA
,
Yackel
,
JJ
.
2008
.
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
113
:
C08026
.
Ertel
,
JR
,
Hedges
,
JI
,
Devol
,
AH
,
Richey
,
JE
,
Ribeiro
,
M
.
1986
.
Dissolved humic substances in the Amazon River system
.
Limnology and Oceanography
31
:
739
754
.
Fichot
,
CG
,
Benner
,
R
.
2014
.
The fate of terrigenous dissolved organic carbon in a river-influenced ocean margin
.
Global Biogeochemical Cycles
28
:
300
318
.
Godin
,
P
,
Macdonald
,
RW
,
Kuzyk
,
ZA
,
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
:
1848
1869
.
Gosselin
,
M
,
Legendre
,
L
,
Demers
,
S
,
Ingram
,
RG
.
1985
.
Responses of sea-ice microalgae to climatic and fortnightly tidal energy inputs (Manitounuk Sound, Hudson Bay)
.
Canadian Journal of Fisheries and Aquatic Sciences
42
:
999
1006
.
Granskog
,
MA
,
Kuzyk
,
ZA
,
Azetsu-Scott
,
K.
,
Macdonald
,
R. W
.
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
.
Granskog
,
MA
,
Macdonald
,
RW
,
Kuzyk
,
ZA
,
Senneville
,
S
,
Mundy
,
C
,
Barber
,
DG
,
Stern
,
GA
,
Saucier
,
F
.
2009
.
Coastal conduit in southwestern Hudson Bay (Canada) in summer: Rapid transit of freshwater and significant loss of colored dissolved organic matter
.
Journal of Geophysical Research
114
:
C08012
.
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
:
2032
2050
.
Guéguen
,
C
,
Mokhtar
,
M
,
Perroud
,
A
,
McCullough
,
G
,
Papakyriakou
,
T
.
2016
.
Mixing and photoreactivity of dissolved organic matter in the Nelson/Hayes estuarine system (Hudson Bay, Canada)
.
Journal of Marine Systems
161
:
42
48
.
Guillemette
,
F
,
Bianchi
,
TS
,
Spencer
,
RGM
.
2017
.
Old before your time: Ancient carbon incorporation in contemporary aquatic foodwebs
.
Limnology and Oceanography
62
:
1682
1700
.
Hansell
,
DA
,
,
D
,
Bates
,
NR
.
2004
.
Degradation of terrigenous dissolved organic carbon in the western Arctic Ocean
.
Science
304
:
858
861
.
Hedges
,
JI
,
Keil
,
RG
,
Benner
,
R
.
1997
.
What happens to terrestrial organic matter in the ocean?
Organic Geochemistry
27
(
5
6
):
195
212
.
Herlemann
,
DPR
,
Manecki
,
M
,
Meeske
,
C
,
Pollehne
,
F
,
Labrenz
,
M
,
Schulz-Bull
,
D
,
Dittmar
,
T
,
Jurgens
K.
2014
.
Uncoupling of bacterial and terrigenous dissolved organic matter dynamics in decomposition experiments
.
PLOS ONE
9
(
4
):
e93945
.
Hernes
,
PJ
,
Benner
,
R
.
2002
.
Transport and diagenesis of dissolved and particulate terrigenous organic matter in the North Pacific Ocean
.
Deep Sea Research I
49
(
12
):
2119
2132
.
Hochheim
,
KP
,
Barber
,
DG
.
2014
.
An update on the ice climatology of the Hudson Bay system
.
Arctic Antarctic and Alpine Research
46
:
66
83
.
Holmes
,
RM
,
McClelland
,
JW
,
Peterson
,
BJ
,
Tank
,
SE
,
Bulygina
,
E
,
Eglinton
,
TI
,
Gordeev
,
VV
,
Gurtovaya
,
TY
,
Raymond
,
PA
,
Repeta
,
DJ
,
Staples
,
R
,
Striegl
,
RG
,
Zhulidov
,
AV
,
Zimov
,
SA
.
2011
.
Seasonal and annual fluxes of nutrients and organic matter from large rivers to the Arctic ocean and surrounding seas
.
Estuaries and Coasts
35
:
369
382
.
Holmes
,
RM
,
McClelland
,
JW
,
Raymond
,
PA
,
Frazer
,
BB
,
Peterson
,
BJ
,
Stieglitz
,
M
.
2008
.
Lability of DOC transported by Alaskan rivers to the Arctic Ocean
.
Geophysical Research Letters
35
(
3
):
L03402
.
Holmes
,
RM
,
McClelland
,
JW
,
Tank
,
SE
,
Spencer
,
RGM
,
Shiklomanov
,
AI
.
2018
.
Arctic Great Rivers Observatory
.
Water Quality Dataset, Version 20200106
.
Available at
https://www.arcticgreatrivers.org/data.
Hudon
,
C
,
Morin
,
R
,
Bunch
,
J
,
Harland
,
R
.
1996
.
Carbon and nutrient output from the Great Whale River (Hudson Bay) and a comparison with other rivers around Quebec
.
Canadian Journal of Fisheries and Aquatic Sciences
53
:
1513
1525
.
Hutchins
,
RHS
,
Prairie
,
YT
,
del Giorgio
,
PA
.
2019
.
Large-scale landscape drivers of CO2, CH4, DOC, and DIC in boreal river networks
.
Global Biogeochemical Cycles
33
:
125
142
.
Ingram
,
RG
,
Larouche
,
P
.
1987
.
Variability of an under-ice river plume in Hudson Bay
.
Journal of Geophysical Research
92
:
9541
9547
.
Ittekkot
,
V
.
1988
.
Global trends in the nature of organic matter in river suspensions
.
Nature
332
(
6163
):
436
.
Kaiser
,
K
,
Benner
,
R
,
Amon
,
R
.
2017
.
The fate of terrigenous dissolved organic carbon on the Eurasian shelves and export to the North Atlantic
.
Journal of Geophysical Research: Oceans
122
:
4
22
.
Kawahigashi
,
M
,
Kaiser
,
K
,
Kalbitz
,
K
,
Rodionov
,
A
,
Guggenberger
,
G
.
2004
.
Dissolved organic matter in small streams along a gradient from discontinuous to continuous permafrost
.
Global Change Biology
10
:
1576
1586
.
Kirchman
,
DL
,
Lancelot
,
C
,
Fasham
,
M
,
Legendre
,
L
,
,
B
,
Scott
,
M
.
1993
. Dissolved organic matter in biogeochemical models of the ocean, in
Evans
,
GT
,
Fasham
,
MJR
eds.,
Towards a model of ocean biogeochemical processes
.
Berlin, Germany
:
Springer-Verlag
:
209
225
.
Kirk
,
JL
, St
.
Louis
,
VL
.
2009
.
Multiyear total and methyl mercury exports from two major sub-Arctic rivers draining into Hudson Bay, Canada
.
Environmental Science and Technology
43
(
7
):
2254
2261
.
Kohler
,
H
,
Meon
,
B
,
Gordeev
,
VV
,
Spitzy
,
A
,
Amon
,
RMW
.
2003
. Dissolved organic matter (DOM) in the estuaries of Ob and Yenisei and the adjacent Kara Sea, Russia, in
Stein
,
R
,
Fahl
,
K
,
Fiitterer
,
DK
,
Galimov
,
EM
,
Stepanets
,
OV
eds.,
Siberian river run-off in the kara sea: Characterization, quantification, variability, and environmental significance.Proc. Mar. Sci
. Vol.
6
.
New York
:
Elsevier
:
281
308
.
Kuzyk
,
ZA
,
Macdonald
,
RW
,
Johannessen
,
SC
,
Gobeil
,
C
,
Stern
,
GA
.
2009
.
Towards a sediment and organic carbon budget for Hudson Bay
.
Marine Geology
264
:
190
208
.
Kuzyk
,
ZA
,
Macdonald
,
RW
,
Tremblay
,
JE
,
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
:
163
176
.
Lammers
,
RB
,
Shiklomanov
,
AI
,
Vorosmarty
,
CJ
,
Fekete
,
BM
,
Peterson
,
BJ
.
2001
.
Assessment of contemporary Arctic river runoff based on observational discharge records
.
Journal of Geophysical Research
10
:
63321
63334
.
Lapierre
,
J-F
,
Guillemette
,
F
,
Berggren
,
M
,
Giorgio
,
PA
.
2013
.
Increases in terrestrially derived carbon stimulate organic carbon processing and CO2 emissions in boreal aquatic ecosystems
.
Nature
4
:
3972
.
Lønborg
,
C
,
,
XA
.
2012
.
Recycling versus export of bioavailable dissolved organic matter in the coastal ocean and efficiency of the continental shelf pump
.
Global Biogeochemical Cycles
26
:
GB3018
.
Lu
,
C-J
,
Benner
,
R
,
Fichot
,
CG
,
Fukuda
,
H
,
Yamashita
,
Y
,
Ogawa
,
H
.
2016
.
Sources and transformations of dissolved lignin phenols and chromophoric dissolved organic matter in Otsuchi Bay, Japan
.
Frontiers in Marine Science
3
:
85
.
Macdonald
,
RW
,
Kuzyk
,
ZA
.
2011
.
The Hudson Bay system: A northern inland sea in transition. Preface to the Journal of Marine Systems special issue: The Hudson Bay system
.
Journal of Marine Systems
88
:
337
340
.
MacGilchrist
,
GA
,
Naveira Garabato
,
AC
,
Tsubouchi
,
T
,
Bacon
,
S
,
Torres-Valdés
,
S
,
Azetsu-Scott
,
K
.
2014
.
The arctic ocean carbon sink
.
Deep Sea Research Part I: Oceanographic Research Papers
86
:
39
55
.
Mann
,
PJ
,
Davydova
,
A
,
Zimov
,
N
,
Spencer
,
RGM
,
Davydov
,
S
,
Bulygina
,
E
,
Zimov
,
S
,
Holmes
,
RM
.
2012
.
Controls on the composition and lability of dissolved organic matter in Siberia’s Kolyma River basin
.
Journal of Geophysical Research
117
:
G01028
.
Mann
,
PJ
,
Eglinton
,
TI
,
McIntyre
,
CP
,
Zimov
,
N
,
Davydova
,
A
,
Vonk
,
JE
,
Holmes
,
RM
,
Spencer
,
RGM
.
2015
.
Utilization of ancient permafrost carbon in headwaters of Arctic fluvial networks
.
Nature Communications
6
:
7856
.
Mars
,
JC
,
Houseknech
,
DW
.
2007
.
Quantitative remote sensing study indicates doubling of coastal erosion rate in past 50 years along a segment of the Arctic coast of Alaska
.
Geology
35
:
583
586
.
McClelland
,
JW
,
Déry
SJ
,
Peterson
BJ
,
Holmes
RM
,
Wood
EF
.
2006
.
A pan-arctic evaluation of changes in river discharge during the latter half of the 20th century
.
Geophysical Research
Letters
33
:
L06715
.
McClelland
,
JW
,
Stieglitz
,
M
,
Pan
,
F
,
Holmes
,
RM
,
Peterson
,
BJ
.
2007
.
Recent changes in nitrate and dissolved organic carbon export from the upper Kuparuk River, North Slope, Alaska
.
Journal of Geophysical Research
112
:
G04S60
.
McGuire
,
,
Anderson
,
LG
,
Christensen
,
TR
,
Dallimore
,
S
,
Guo
,
L
,
Hayes
,
DJ
,
Heimann
,
M
,
Lorenson
,
TD
,
MacDonald
,
RW
,
Roulet
,
N
.
2009
.
Sensitivity of the carbon cycle in the Arctic to climate change
.
Ecological Monographs
79
(
4
):
523
555
.
Müller
,
O
,
Seuthe
,
L
,
Bratbak
,
G
,
Paulsen
,
ML
.
2018
.
Bacterial response to permafrost derived organic matter input in an Arctic Fjord
.
Frontiers in Marine Science
5
:
263
.
Mundy
,
CJ
,
Gosselin
,
M
,
Starr
,
M
,
Michel
,
M
.
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
.
Osburn
,
CL
,
Retamal
,
L
,
Vincent
,
WF
.
2009
.
Photoreactivity of chromophoric dissolved organic matter transported by the Mackenzie river to the Beaufort sea
.
Marine Chemistry
115
(
1
2
):
10
20
.
Petrusevich
,
VY
,
Dmitrenko
,
IA
,
Kozlov
,
IE
,
Kirillov
,
SA
,
Kuzyk
,
ZA
,
Komarov
,
AS
,
Heath
,
JP
,
Barber
,
DG
,
Ehn
,
JK
.
2018
.
Tidally-generated internal waves in Southeast Hudson Bay
.
Continental Shelf Research
167
:
65
76
.
Rontani
,
JF
,
Charrière
,
B
,
Sempéré
,
R
,
Doxaran
,
D
,
Vaultier
,
F
,
Vonk
,
JE
,
Volkman
,
JK
.
2014
.
Degradation of sterols and terrigenous organic matter in waters of the Mackenzie Shelf
.
75
:
61
73
.
Rosa
,
E
,
Gaillardet
,
J
,
Hillaire-Marcel
,
C
,
Hélie
,
J-F
,
Richard
,
L-F
.
2012
.
Rock denudation rates and organic carbon exports along a latitudinal gradient in the Hudson, James, and Ungava Bays watershed
.
49
:
742
757
.
Rosenberg
,
D
,
Chambers
,
P
,
Culp
,
J
,
Franzin
,
W
,
Nelson
,
P
,
Salki
,
AG
,
Stainton
,
MP
,
Bodaly
,
RA
,
Newbury
,
RW
.
2005
.
Nelson and Churchill river basins
, in
Benke
AC
,
Cushing
,
CE
eds.,
Rivers of North America
.
Amsterdam, the Netherlands
:
Elsevier
:
853
901
.
RSW-Environnement Illimité Inc
.
2014
.
Conawapa generating station—Stage IV studies—Physical environment of the Nelson River Estuary, hydrodynamic modeling component study, Revision 1. WBS: 1.4.3.3.2.6.2275, Manitoba Hydro File: 00192-11100-0062
.
Report prepared for Manitoba Hydro, Water Resources Engineering Department, Power Planning Division, Power Supply
,
.
Rysgaard
,
S
,
Thastum
,
P
,
Dalsgaard
,
T
,
Christensen
,
PB
,
Sloth
,
NP
.
1999
.
Effects of salinity on NH4+ adsorption capacity, nitrification, and denitrification in Danish estuarine sediments
.
Estuaries
22
(
1
):
21
30
.
Shirokova
,
LS
,
Chupakova
,
AA
,
Chupakov
,
AV
,
Pokrovsky
,
OS
.
2017
.
Transformation of dissolved organic matter and related trace elements in the mouth zone of the largest European Arctic river: Experimental modeling
.
Inland Waters
7
(
3
):
272
282
.
Smith
,
A
,
Delavau
,
C
,
,
T
.
2015
.
Identification of geographical influences and flow regime characteristics using regional water isotope surveys in the lower Nelson River, Canada
.
40
(
1
):
23
35
.
Spencer
,
RGM
,
Mann
,
PJ
,
Dittmar
,
T
,
Eglinton
,
TI
,
McIntyre
,
C
,
Holmes
,
RM
,
Zimov
,
N
,
Stubbins
,
A
.
2015
.
Detecting the signature of permafrost thaw in Arctic rivers
.
Geophysical Research Letters
42
:
2830
2835
.
Stainton
,
T
.
2018
. Sources and transport of particulate matter in the Nelson River system, Manitoba [M.Sc. thesis].
Winnipeg, Manitoba
:
Department of Geological Sciences, University of Manitoba
:
154 (pages) + 44 (Appendices)
.
St. Jacques
,
JM
,
Sauchyn
,
D
.
2009
.
Increasing winter baseflow and mean annual streamflow from possible permafrost thawing in the Northwest Territories, Canada
.
Geophysical Research Letters
36
:
L01401
.
St. Laurent
,
P
,
Straneo
,
F
,
Dumais
,
J-F
,
Barber
,
DG
.
2011
.
What is the fate of the river waters of Hudson Bay?
Journal of Marine Systems
88
:
352
361
.
Tank
,
SE
,
Frey
,
KE
,
Striegl
,
RG
,
Raymond
,
PA
,
Holmes
,
RM
,
McClelland
,
JW
,
Peterson
,
BJ
.
2012
.
Landscape-level controls on dissolved carbon flux from diverse catchments of the circumboreal
.
Global Biogeochemical Cycles
26
:
GB0E02
.
Tranvik
,
LJ
,
Cole
,
JJ
,
Prairie
,
YT
.
2018
.
The study of carbon in inland waters—From isolated ecosystems to players in the global carbon cycle
.
Limnology and Oceanography Letters
3
:
41
48
.
Vonk
,
JE
,
Tank
,
SE
,
Mann
,
PJ
,
Spencer
,
RGM
,
Treat
,
CC
,
Striegl
,
RG
,
Abbott
,
BW
,
Wickland
,
KP
.
2015
.
Biodegradability of dissolved organic carbon in permafrost soils and aquatic systems: A meta-analysis
.
Biogeosciences
12
:
6915
6930
.
Walvoord
,
MA
,
Striegl
,
RG
.
2007
.
Increased groundwater to stream discharge from permafrost thawing in the Yukon River basin: Potential impacts on lateral export of carbon and nitrogen
.
Geophysical Research Letters
34
:
L12402
.
Wang
,
R
,
McCullough
,
GK
,
Gunn
,
GG
,
Hochheim
,
KP
,
Dorostkar
,
A
,
Sydor
,
K
,
Barber
,
DG
.
2012
.
An observational study of ice effects on Nelson River estuarine variability, Hudson Bay, Canada
.
Continental Shelf Research
47
:
68
77
.
Ward
,
CP
,
Cory
,
RM
.
2015
.
Chemical composition of dissolved organic matter draining permafrost soils
.
Geochimica et Cosmochimica Acta
167
:
63
79
.
Ward
,
ND
,
Bianchi
,
TS
,
Sawakuchi
,
HO
,
Gagne-Maynard
,
W
,
Cunha
,
AC
,
Brito
,
DO
,
Neu
,
V
,
Valerio
,
AM
,
Silva
,
R
,
Krusche
,
AV
,
Richey
,
JE
,
Keil
,
RG
.
2016
.
The reactivity of plant-derived organic matter and the potential importance of priming effects along the lower Amazon River
.
Journal of Geophysical Research Biogeosciences
121
:
1522
1539
.
Ward
,
ND
,
Sawakuchi
,
HO
,
Richey
,
JE
,
Keil
,
RG
,
Bianchi
,
TS
.
2019
.
Enhanced aquatic respiration associated with mixing of clearwater tributary and turbid Amazon River waters
.
Frontiers in Earth Science
7
:
101
.
Wickland
,
KP
,
Aiken
,
GR
,
Butler
,
K
,
Dornblaser
,
MM
,
Spencer
,
RGM
,
Striegl
,
RG
.
2012
.
Biodegradability of dissolved organic carbon in the Yukon River and its tributaries: Seasonality and importance of inorganic nitrogen
.
Global Biogeochemical Cycles
26
:
GB0E03
.
Yasunaka
,
S
,
Siswanto
,
E
,
Olsen
,
A
,
Hoppema
,
M
,
Watanabe
,
E
,
Fransson
,
A
,
Chierici
,
M
,
Murata
,
A
,
Lauvset
,
SK
,
Wanninkhof
,
R
,
Takahashi
,
T
.
2018
.
Arctic Ocean CO2 uptake: An improved multiyear estimate of the air–sea CO2 flux incorporating chlorophyll a concentrations
.
Biogeosciences
15
:
1643
1661
.

How to cite this article: Kazmiruk, ZV, Capelle, DW, Kamula, CM, Rysgaard, S, Papakyriakou, T, Kuzyk, ZA. 2021. High biodegradability of riverine dissolved organic carbon in late winter in Hudson Bay, Canada. Elementa: Science of Anthropocene 9(1). DOI: https://doi.org/10.1525/elementa.2020.00123.

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

Associate Editor: Jeff Bowman, Integrative Oceanography Division, Scripps Institution of Oceanography, CA, USA

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/.