Journal topic
Biogeosciences, 16, 4535–4553, 2019
https://doi.org/10.5194/bg-16-4535-2019
Biogeosciences, 16, 4535–4553, 2019
https://doi.org/10.5194/bg-16-4535-2019

Research article 29 Nov 2019

Research article | 29 Nov 2019

Comparisons of dissolved organic matter and its optical characteristics in small low and high Arctic catchments

Comparisons of dissolved organic matter and its optical characteristics in small low and high Arctic catchments
Caroline Coch1,2,3, Bennet Juhls4, Scott F. Lamoureux5, Melissa J. Lafrenière5, Michael Fritz1, Birgit Heim1, and Hugues Lantuit1,2 Caroline Coch et al.
• 1Alfred-Wegener-Institute Helmholtz Centre for Polar and Marine Research, Telegrafenberg A45, 14473 Potsdam, Germany
• 2Institute of Earth and Environmental Science, University of Potsdam, Karl-Liebknecht-Straße 24/25, 14476 Potsdam, Germany
• 3World Wildlife Fund, The Living Planet Centre, Rufford House, Brewery Road, Woking, Surrey, GU214LL, UK
• 4Institute for Space Sciences, Department of Earth Sciences, Freie Universität Berlin, 12165 Berlin, Germany
• 5Department of Geography and Planning, Mackintosh-Corry Hall, Queen's University, Kingston, Ontario K7L 3N6, Canada

Correspondence: Caroline Coch (coch.caroline@gmail.com)

Abstract

Climate change is affecting the rate of carbon cycling, particularly in the Arctic. Permafrost degradation through deeper thaw and physical disturbances results in the release of carbon dioxide and methane to the atmosphere and to an increase in lateral dissolved organic matter (DOM) fluxes. Whereas riverine DOM fluxes of the large Arctic rivers are well assessed, knowledge is limited with regard to small catchments that cover more than 40 % of the Arctic drainage basin. Here, we use absorption measurements to characterize changes in DOM quantity and quality in a low Arctic (Herschel Island, Yukon, Canada) and a high Arctic (Cape Bounty, Melville Island, Nunavut, Canada) setting with regard to geographical differences, impacts of permafrost degradation, and rainfall events. We find that DOM quantity and quality is controlled by differences in vegetation cover and soil organic carbon content (SOCC). The low Arctic site has higher SOCC and greater abundance of plant material resulting in higher chromophoric dissolved organic matter (cDOM) and dissolved organic carbon (DOC) than in the high Arctic. DOC concentration and cDOM in surface waters at both sites show strong linear relationships similar to the one for the great Arctic rivers. We used the optical characteristics of DOM such as cDOM absorption, specific ultraviolet absorbance (SUVA), ultraviolet (UV) spectral slopes (S275–295), and slope ratio (SR) for assessing quality changes downstream, at base flow and storm flow conditions, and in relation to permafrost disturbance. DOM in streams at both sites demonstrated optical signatures indicative of photodegradation downstream processes, even over short distances of 2000 m. Flow pathways and the connected hydrological residence time control DOM quality. Deeper flow pathways allow the export of permafrost-derived DOM (i.e. from deeper in the active layer), whereas shallow pathways with shorter residence times lead to the export of fresh surface- and near-surface-derived DOM. Compared to the large Arctic rivers, DOM quality exported from the small catchments studied here is much fresher and therefore prone to degradation. Assessing optical properties of DOM and linking them to catchment properties will be a useful tool for understanding changing DOM fluxes and quality at a pan-Arctic scale.

1 Introduction

Climate change has important impacts on carbon cycling, particularly in the Arctic. Approximately 1300 Gt of organic carbon is stored in permafrost soils in the Northern Hemisphere (Hugelius et al., 2014), which is 40 % more than currently circulating in the atmosphere. Thawing permafrost and deepening of the active layer lead to the mobilization of this carbon (Osterkamp, 2007; Woo et al., 2008), the release of carbon dioxide (CO2) and methane (CH4) to the atmosphere (Schaefer et al., 2014), and an increase in riverine dissolved organic carbon (DOC) fluxes (Frey and Smith, 2005; Le Fouest et al., 2018). Also associated with warming is the development of surface (physical) disturbances such as active layer detachments or retrogressive thaw slumps (Lacelle et al., 2010; Lamoureux and Lafrenière, 2009; Lewkowicz, 2007; Ramage et al., 2018) and thermal perturbation of the subsurface (Lafrenière and Lamoureux, 2013). As these processes influence freshwater systems, they ultimately have impacts on the biological production and the biogeochemistry of the Arctic Ocean. The six largest Arctic rivers (Mackenzie, Yukon, Ob', Yenisey, Lena, Kolyma) drain 53 % of the Arctic Ocean drainage basin (Holmes et al., 2012) and transport huge amounts of nutrients and dissolved organic matter (DOM) to the ocean. However, there are limited flux estimates and information on DOM quality available for the remaining 47 %, which is sourced by smaller watersheds. “Small” in this context refers to smaller than the large Arctic rivers, as the actual size distribution of these watersheds remains unknown.

Terrigenous DOM is an important source of DOC originating from allochthonous (terrestrial such as soil and plants) and autochthonous (in situ production) sources (Aiken, 2014), and it is modified by biotic and abiotic processes during its lateral transport to the ocean (Tank et al., 2018; Vonk et al., 2015a, b). Yet little is known about the transformation of DOM along short distances in small catchments. The composition and the vulnerability of riverine DOM to transformation is influenced by several factors such as soil organic matter and vegetation, sorption processes in the mineral layers, and biodegradation and photodegradation processes (Cory et al., 2014; Mann et al., 2012; Vonk et al., 2015b; Ward and Cory, 2015; Ward et al., 2017). Chromophoric or coloured dissolved organic matter (cDOM) is a fraction of DOM, which absorbs light in the ultraviolet and visible wavelengths (Green and Blough, 1994). Optical characteristics of cDOM such as absorption coefficients and spectral slopes can serve as proxies for DOM molecular weight and aromaticity, which in turn can help to characterize the lability of DOM (Helms et al., 2008; Neff et al., 2006; Spencer et al., 2009; Striegl et al., 2005; Weishaar et al., 2003). High SUVA values (UV specific absorbance at 254 nm) in combination with low S275–295 (spectral cDOM slope between 275 and 295 nm) values indicate “fresh” DOM, or systems of shorter residence time receiving a greater input of fresh DOM from the catchment area. In contrast, low SUVA and high S275–295 are considered indicators of limited inputs of fresh DOM, a higher relative contribution of autochthonous DOM, greater exposure to photobleaching, and longer residence time (Anderson and Stedmon, 2007; Fichot and Benner, 2012; Fichot et al., 2013; Helms et al., 2008; Whitehead et al., 2000).

Previous studies have focused on characterizing the cDOM–DOC relationship for the large Arctic rivers and coastal shelf areas, which exhibit a strong seasonality (Spencer et al., 2008; Stedmon et al., 2011; Walker et al., 2013). Some studies have investigated cDOM–DOC relationships in smaller Arctic catchments: Dvornikov et al. (2018) examined cDOM characteristics in surface waters of the Yamal Peninsula, and cDOM–DOC relationships were examined for subarctic catchments (Balcarczyk et al., 2009; Cory et al., 2015; Larouche et al., 2015; O'Donnell et al., 2014) and the high Arctic (Fouché et al., 2017; Wang et al., 2018). Optical parameters have also been used to assess the impact of permafrost disturbance on stream geochemistry in Alaska (Abbott et al., 2014; Larouche et al., 2015) and the Northwest Territories (NWT), Canada (Littlefair et al., 2017). As most studies focused on downstream reaches, knowledge on the spatial variability across catchments is limited. To our knowledge, no study has examined this relationship in a low Arctic setting or attempted to resolve geographic differences between the low and high Arctic.

Here, we investigate cDOM and DOC in surface waters in the low Arctic (Herschel Island, Yukon, Canada) and the high Arctic (Cape Bounty, Melville Island, Nunavut, Canada). The aim of this study is to (1) compare the variability and relation of DOC concentration and cDOM in low and high Arctic surface water environments, (2) to investigate changes in DOM composition along longitudinal stream profiles and with regard to permafrost disturbances, and (3) examine changes in DOM concentration and composition throughout the summer season with occasional rainfall events.

2 Study area

This study was carried out in two Arctic locations, Herschel Island in the low Arctic and at the Cape Bounty Arctic Watershed Observatory (CBAWO), Melville Island, in the high Arctic (Fig. 1a). Herschel Island (Yukon, Canada) is located at 6935 N and 13905 W in the Beaufort Sea off the Yukon coast. The island is composed of unconsolidated and fine-grained marine and glaciogenic sediments as it was formed by the Laurentide Ice Sheet (Rampton, 1982). The island is situated in the zone of continuous permafrost with ground ice content between 30 % and 60 % for the entire island. Physical permafrost degradation typically occurs in the form of retrogressive thaw slumps (Lantuit and Pollard, 2008) and active layer detachments (Coch et al., 2020). Ramage et al. (2019) reported mean active layer depths of 52.2±20.2 cm. Soil organic carbon content (SOCC) for valleys on the eastern side of Herschel Island was estimated to be 11.4±3.7 kg m2 at 0–30 cm depth and 26.4±8.9 kg m2 at 0–100 cm depth with a C:N ratio of 12.9±2.2 at 0–100 cm depth (Ramage et al., 2019). The dominant vegetation type is lowland tundra (Myers-Smith et al., 2011; Smith et al., 1989) and can be classified into subzone E (CAVM, 2003), which corresponds to the low Arctic. The mean annual air temperature and yearly precipitation between 1971 and 2000 at Komakuk Beach, the nearest long-term meteorological station ∼40 km away from our study site, are −11C and 161.3 mm respectively. The mean July temperature is 7.8 C and average precipitation is 27.3 mm (Environment and Climate Change Canada, 2018). Snowmelt is the largest hydrological event of the year occurring in May to early June. Summer base flow after mid-June is controlled by rainfall events (Coch et al., 2018). The active layer freezes up by mid-November (Burn, 2012). The studied catchments unofficially named Ice Creek West (1.4 km2) and Ice Creek East (1.6 km2) are adjacent to each other and merge before draining into the Beaufort Sea (Fig. 1b, Table 1). Both sampled ponds in Ice Creek West are <1 ha, and there are degrading ice-wedge polygons present in the headwaters of Ice Creek West (Coch et al., 2020).

Figure 1Maps of the study area showing (a) the location of Herschel Island and Cape Bounty in the Canadian Arctic including the Circumpolar Arctic Vegetation Map (CAVM) bioclimatic zones (CAVM, 2003), (b) the studied catchments Ice Creek West and Ice Creek East on Herschel Island, and (c) the studied catchments Boundary River with its subcatchment Robin Creek (dashed watershed), West River, and East River. The watershed names are indicated with numbers, and the general flow direction is southwards towards the ocean. Note that samples from flowing water (rivers and streams) are indicated by circles, whereas samples from standing water (ponds and lakes) are indicated by triangles. Yellow marks locations where DOC concentration and cDOM measurements are available, while only DOC concentrations are available at red locations. The basemap for (a) is the World Hillshade layer by Esri (created 2015, last updated 2019). The background images for (b) and (c) are true colour mosaics (Herschel Island: WorldView-3 quasi-true colour RGB composite, acquired on 8 August 2015; Cape Bounty: Sentinel-2 quasi-true colour RGB composite, acquired on 7 August 2016).

Table 1Characteristics of studied watersheds on Herschel Island (low Arctic) and Cape Bounty (high Arctic) showing catchment size (km2), channel length (km), circumarctic vegetation map (CAVM) bioclimatic zone (CAVM, 2003), soil organic carbon content (SOCC) (Hugelius et al., 2013; Ramage et al., 2019), and maximum catchment elevation above sea level (m a.s.l.).

The CBAWO is situated on the south coast of Melville Island (Nunavut, Canada) at 7455 N and 10935 W. The geology is characterized by Devonian sandstone and siltstone bedrock overlain by Quaternary marine and glacial sediments (Hodgson et al., 1984). The soils are categorized as Cryosols with a thin organic horizon. The site is situated in the zone of continuous permafrost, and active layer depths typically range from 50 to 70 cm (Lafrenière et al., 2013). Permafrost degradation such as deep thaw and physical disturbances have altered hydrochemical fluxes of the rivers (Lamoureux and Lafrenière, 2017). The vegetation cover is patchy with polar semi-desert, mesic tundra, and wet sedge meadows (Edwards and Treitz, 2018), and it falls into subzones B and C (CAVM, 2003). Soil organic carbon is estimated to be 3.0 kg m2 at 0–30 cm depth and 10.2 kg m2 at 0–100 cm depth (Hugelius et al., 2013), with a C:N ratio of 10.0 at 0–100 cm depth (ADAPT, 2014). The nearest long-term meteorological station is located ∼300 km away, at Mould Bay (NWT). Between 1971 and 2000, the mean annual air temperature and precipitation were −17.5C and 111 mm respectively. The mean July temperature is 4.0 C and precipitation is 13.5 mm. Snowmelt and nival runoff start in early to mid-June with base flow establishing around early to mid-July. Refreezing of the active layer starts late August or early September (Lamoureux and Lafrenière, 2017; Lewis et al., 2012). Samples were taken downstream in Boundary River (152.5 km2), its sub-catchment Robin Creek (14.8 km2), and the neighbouring watersheds West River (8.6 km2) and East River (12.0 km2) (all unofficial names). There is an active retrogressive thaw slump in the Robin Creek watershed and a number of recent (since 2007) active layer detachments and other disturbances in the other watersheds. The sampled lakes and ponds cover a range of sizes from below 1 ha in West River to the larger downstream West and East lakes ( 140–160 ha).

3 Methods

3.1 Field methods and hydrochemistry

To explore downstream changes in DOM across regions, we used a transect approach in this study. Samples were taken along longitudinal stream profiles in catchments, and additionally samples were collected from standing water bodies (ponds and lakes) (Fig. 1). We also obtained discharge records and water samples from the outflow of both catchments at Herschel Island over the course of the 2016 season as detailed below.

Field work on Herschel was carried out in July–August 2016. We measured discharge using a cutthroat flume equipped with a U20 Onset Hobo level logger in Ice Creek West. Discharge data at 30 min intervals are available from 15 May 2016 and at 5 min intervals after 22 July 2016 (Coch et al., 2018, 2020). In Ice Creek East, discharge was determined using the area velocity method in combination with a U20 Onset Hobo level logger (see Coch et al., 2020, for a detailed description). Data in Ice Creek East are available at 5 min intervals after 25 July 2016. Weather data are available from the local Environment and Climate Change Canada station and from an additional station deployed in Ice Creek West during the summer. Water samples were collected in bottles triple rinsed with sample at the outflow of both streams between 20 July and 10 August. At the outflow of Ice Creek West, water samples were collected using an automatic water sampler (ISCO 3700) at a 12 h interval between 25 July and 10 August and more frequently during rainfall events (between 1 and 3 h). Prior to the automatic sampling, and also in Ice Creek East, water samples were taken manually once per day. We collected water samples along longitudinal profiles of the channels (11 in Ice Creek West, 12 in Ice Creek East) starting in the headwaters (∼2000 m distance from the outflow) and following the channel downstream (Fig. 1). Longitudinal profiles were sampled three times in Ice Creek West (20, 25, and 30 July) and once in Ice Creek East (30 July). Samples of flowing water are available from Ice Creek West (n=90), Ice Creek East (n=32), and the alluvial fan (n=8). Standing water samples (n=4) were collected from two ponds in the Ice Creek West catchment.

The field work at CBAWO took place during August 2017. All samples were collected manually after triple rinsing the bottle with sample water. Similar to the Herschel field work, we collected samples along longitudinal stream profiles. Robin Creek is a subcatchment of Boundary River (Fig. 1), where stream samples were collected at six locations downstream of a retrogressive thaw slump. Three lakes were also sampled in the Boundary River catchment, and two samples are from the main channel of the Boundary River. A total of 21 stream samples and nine samples from lakes and ponds are available from the West River catchment, some of which were collected after the rainfall event on 12 August 2017. In East River, four samples are available from the stream and eight samples from standing water bodies.

Within 24 h of sampling, electrical conductivity (EC) and pH were measured in the field laboratory. After collection, water samples were filtered through pre-rinsed 0.7 µm GF/F syringe filters and were then stored cool and dark for transport to the Alfred Wegener Institute, University of Hamburg and Geoscience Research Centre GFZ, Germany, where analysis for DOC and cDOM was carried out. Water samples for absorbance measurements were kept in brown glass bottles without acidifying. Samples for DOC analyses were acidified with HCl (30 % Suprapur) prior to the measurements. In 2016, DOC measurements were performed on a Shimadzu TOC-L analyzer with a TNM-L module (University of Hamburg), whereas a Shimadzu TOC-VCPH analyzer was used in 2017 (AWI). The error for these measurements is below 10 %.

Inorganic carbon was sparged out using synthetic air prior to the measurement. As we had a shortage of HCl in the field in 2016, 82 of the samples were frozen and acidified upon return to Germany. After new acid was acquired later in the summer, sample duplicates (n=47) were processed directly in the field and also frozen. The frozen duplicate was thawed and acidified upon return to determine the effect of different sample treatment (Coch et al., 2018). There is a significant linear relationship (p<0.05, n=47, R2=0.87) between DOC concentrations of unfrozen and frozen sample duplicates. Samples that were frozen in the field, and subsequently thawed and acidified upon return to Germany, showed lower DOC concentrations (by 13 %) than samples that were acidified directly in the field and kept unfrozen. We corrected the frozen samples for this offset (Supplement, Fig. S1). In both years, deionized water used in the field was also analyzed as a blank following the same procedure.

The absorbance was measured for the wavelengths from 200 to 800 nm with a 1 nm increment using a LAMBDA 950 UV-vis spectrophotometer (GFZ Potsdam). The measurements were performed in duplicates using a 5 cm cuvette and Milli-Q water as a blank. Some of the water samples showed fine particles precipitated in the sample bottle. They appeared in the form of small thin flakes, which partly remained in suspension or accumulated at the bottom of the flask. This precipitation occurred after the samples were filtered through 0.7 µm glass fibre filters, transported to the laboratory for storage of about 4 weeks. This information was documented, and although absorbance spectra were measured, they were not further analyzed in this study as interference of the spectral characteristics by the particles might have occurred. This was the case for 25 (out of 55) samples at Cape Bounty and for eight samples (out of 134) at Herschel Island.

The Napierian spectral absorption coefficient of cDOM (acDOM(λ)) was calculated out of the mean of each duplicate with

$\begin{array}{}\text{(1)}& {a}_{\mathrm{cDOM}}\left(\mathit{\lambda }\right)\phantom{\rule{0.25em}{0ex}}\left({\mathrm{m}}^{-\mathrm{1}}\right)=\mathrm{2.303}\frac{\left({A}_{\mathrm{sample}}\left(\mathit{\lambda }\right)-{A}_{\mathrm{reference}}\left(\mathit{\lambda }\right)\right)}{L},\end{array}$

where Asample is the absorbance of the sample, Areference the absorbance of the Milli-Q reference, and L is the optical path length of the used cuvette in the spectrophotometer (L=0.05 m). The decadal absorption is multiplied by 2.303, which is the conversion factor of base 10 to base e logarithm to derive the Napierian absorption coefficient that is used for cDOM.

$\begin{array}{}\text{(2)}& \mathrm{cDOM}\left(\mathit{\lambda }\right)\phantom{\rule{0.25em}{0ex}}\left({\mathrm{m}}^{-\mathrm{1}}\right)\frac{\mathrm{2.303}\cdot \left({A}_{\mathrm{sample}}\left(\mathit{\lambda }\right)-{A}_{\mathrm{reference}}\left(\mathit{\lambda }\right)\right)}{L}\end{array}$

The absorption was corrected for scatter using a baseline correction by subtracting acDOM(700) (Hancke et al., 2014; Helms et al., 2008). At that wavelength, absorption by cDOM is assumed negligible (Mitchell et al., 2002). Spectral slope values of acDOM for wavelength ranges from 275 to 295 nm (S275–295) and 350 to 400 nm (S350–400) were calculated using Eq. (2) and a non-linear fit. These spectral slope values indicate photochemical or microbial alteration of DOM (Helms et al., 2008). The ratio of both spectral slopes (S275–295 : S350–400) defines the slope ratio (SR). The SUVA (L mg−1 m−1) was calculated by dividing the decadal absorption (A254L) at 254 nm (m−1), where A254 is the absorbance at 254 nm and L the optical path length of the used cuvette in the spectrophotometer by DOC (mg L−1). Both parameters, SR and SUVA, have been related to the relative molecular weight and aromaticity of DOM (Helms et al., 2008; Weishaar et al., 2003).

A cDOM absorption spectrum, acDOM(λ), is generally expressed as an exponential function:

$\begin{array}{}\text{(3)}& {a}_{\mathrm{cDOM}}\left(\mathit{\lambda }\right)\left(\mathit{\lambda }\right)={a}_{\mathrm{cDOM}}\left(\mathit{\lambda }\right)\cdot {e}^{-S\left(\mathit{\lambda }-{\mathit{\lambda }}_{\mathrm{0}}\right)},\end{array}$

where λ0 is the absorption coefficient at the reference wavelength and S is the spectral slope of acDOM(λ) for the chosen wavelength range. To compare our data with different studies, we converted absorption coefficient values reported in various studies to acDOM350 using an interpolation method developed by Massicotte et al. (2017). Throughout the paper all data are reported as mean ± standard deviation.

3.2 Statistical analyses

We used RStudio (version 1.0.153) to perform statistical tests (RStudio Team, 2016). Normality was tested using the Shapiro–Wilk normality test. To determine the difference in means of two populations, we applied the Welch's two-sample t test if the data were normally distributed with unequal variances. In the case where data were not normally distributed, we used the Wilcoxon–Mann–Whitney test. To measure the relationship between two variables, we used the Pearson correlation coefficient for normally distributed data and the Spearman rank correlation if the data were not normally distributed.

4 Results

4.1 DOM characteristics and relationships in and across low Arctic (Herschel Island) and high Arctic (Cape Bounty) catchments

Comparing DOC concentrations and cDOM absorption between both study sites, significant differences can be found. The cDOM absorption is significantly higher (p<0.05) in samples from Herschel Island compared with Cape Bounty across the entire spectrum (Fig. 2a) with acDOM350 of 14.5±5.1 and 5.5±4.9 m−1 respectively (Fig. 2b). DOC concentrations show a similar pattern with significantly higher values (p<0.05) on Herschel Island (10.0±1.6 mg L−1) compared to Cape Bounty (2.5±2.0 mg L−1) (Fig. 2c, Table 2).

Table 2Descriptive statistics (mean ± standard deviation) of dissolved organic carbon, DOC (mg L−1), specific ultraviolet absorbance, SUVA (L mg−1 m−1), coloured dissolved organic matter (cDOM) absorption at 350 nm, acDOM350 (m−1), cDOM slope S275–295 (10−3 nm−1), slope ratio SR, electrical conductivity EC (µS cm−1), pH, and the number (n) of all samples/samples with cDOM absorption measurements. The statistics are given for specific streams, samples from flowing waters, standing waters, and all samples on Herschel Island (HE) and Cape Bounty (CB) respectively. The symbols “>” and “<” indicate significant inter-group differences at the alpha = 0.95 level. When the inter-group differences are significantly different at the alpha = 0.99 level, then they are bold. When the difference is not significant, “” is used.

NA – not available; n/a – not applicable.

Comparing the streams on Herschel Island (Fig. 3b), the highest DOC and acDOM350 values are found in the headwaters of Ice Creek West. Ice Creek West has significantly higher (p<0.05) values in DOC (10.4±1.5 mg L−1) and acDOM350 (16.1±5.4 m−1) than Ice Creek East, which are 8.7±1.1 mg L−1 and 11.1±1.8 m−1 respectively. At Cape Bounty, West River shows the highest DOC (2.5±1.7 mg L−1) and acDOM350 (8.5±5.2 m−1) compared to other sampled streams. The highest values of DOC and acDOM350 values of flowing are recorded in West River after the 8 August rainfall event.

Figure 2Dissolved organic matter (DOM) quality and quantity for the sites from Herschel Island (HE, in orange) and Cape Bounty (CB, in blue). (a) Average cDOM absorption (m−1) for the wavelengths (λ) between 250 and 700 nm. The coloured shaded areas represent the standard deviation from the mean (solid line). Boxplots of (b) coloured dissolved organic matter (cDOM) absorption at 350 nm acDOM350 (m−1), (c) dissolved organic carbon, DOC (mg L−1), (d) specific ultraviolet absorbance, SUVA (L mg−1 m−1), (e) cDOM slope S275–295 (10−3 nm−1), and slope ratio SR. The plots depict distributions for all samples (all) and subsets of flowing water (F) and standing water (S) for each of the site.

Figure 3Absorption of coloured dissolved organic matter (cDOM) at 350 nm (m−1) versus dissolved organic carbon (DOC) concentration (mg L−1) for (a) all sites; (b) sites on Herschel Island depicting the sampling locations Ice Creek West (ICW) upstream, downstream, and ponds; Ice Creek East (ICE) upstream, downstream, and alluvial fan; and (c) sites at Cape Bounty (West River, East River, Boundary River, Robin Creek). Note that flowing water is indicated by a circle while standing water such as lakes or ponds is indicated by a triangle. The cDOM-to-DOC relationships are divided into two different groups (c).

For both study areas, DOC concentrations are substantially higher in standing waters compared to flowing water (Fig. 2c). Furthermore, on Herschel Island, generally upstream waters show higher DOC concentrations and acDOM350 compared to downstream waters. At Cape Bounty, samples from standing waters in the East River catchment show the highest DOC concentrations and acDOM350 values.

Mean S275–295 on Herschel Island is generally higher ($\mathrm{16.4}±\mathrm{1.5}×{\mathrm{10}}^{-\mathrm{1}}$ nm−1) compared to Cape Bounty ($\mathrm{14.8}±\mathrm{3.2}×{\mathrm{10}}^{-\mathrm{1}}$ nm−1), whereas SUVA values show a broader range at Cape Bounty (from 1.35 to 5.16 mg L−1 m−1) compared to Herschel Island (from 2.0 to 4.3 mg L−1 m−1) (Fig. 2d, Table 2). At Cape Bounty, standing water samples show significantly larger spectral slopes and slope ratios (p<0.05) and mostly smaller SUVA (p<0.05) than flowing water samples (Fig. 2d–f). We observed a significant positive relationship (ρ=0.78, p<0.05) between acDOM350 and DOC concentration for all samples at both sites (Fig. 3a). Whereas on Herschel Island the relationship between acDOM350 and DOC follows one linear trend (ρ=0.72, p<0.05), the relationship at Cape Bounty is broadly separated into two groups, which we schematically indicate by solid ellipsoids that correspond to flowing and standing water. Correlations for both groups are significant (<0.05) and show different spectral slopes. One sample identified as standing water falls into the group of flowing water. We identify this sample as an outlier likely affected by a local DOM source.

We found a strong negative relationship between SUVA and S275–295 among all water samples from both locations ($\mathit{\rho }=-\mathrm{0.64}$, p<0.05, Fig. 4a). This relationship is even stronger when only samples from Herschel are considered ($\mathit{\rho }=-\mathrm{0.72}$, p<0.05), whereas samples from Cape Bounty show higher deviations. Reported SUVA and S275–295 from Walker et al. (2013) indicate a similar relationship among the large Arctic rivers (Ob', Lena, Yenisei, Kolyma and Mackenzie). SUVA and S275 of the large Arctic rivers are located in the centre of our reported values from Herschel Island and Cape Bounty (highlighted in Fig. 4a).

Figure 4Slope of coloured dissolved organic matter ultraviolet cDOM UV absorption 275–295 (10−3 nm−1) versus specific ultraviolet absorbance SUVA (L mg−1 m−1) for (a) all sites, (b) sites on Herschel Island depicting the sampling locations Ice Creek West (ICW) upstream, downstream, and ponds; Ice Creek East (ICE) upstream, downstream, and alluvial fan; and (c) sites at Cape Bounty (West River, East River, Boundary River, Robin Creek). Note that flowing water is indicated by a circle while standing water such as lakes or ponds is indicated by a triangle. Midflow SUVA and S275–295 values for the large Arctic rivers (Walker et al., 2013) are added to (a).

Whereas no difference was found between standing and flowing water on Herschel Island, pH and EC values were significantly higher in standing water at Cape Bounty than in flowing water (Table 2). Robin Creek showed the highest EC values and the largest variability (145±213µS cm−1) of the Cape Bounty rivers, whereas West River showed the overall lowest EC values (60±17µS cm−1). On Herschel Island, both adjacent streams show EC and pH values in the same order of magnitude with a slight decrease at the alluvial fan outflows.

4.2 Hydrochemical and DOM patterns along longitudinal stream transects

The studied streams followed different hydrochemical patterns from upstream to downstream (Fig. 5). EC and pH are significantly higher (p<0.05) in surface waters on Herschel Island (1050±370 and 8.2±0.2µS cm−1) than at Cape Bounty (137±136 and 7.2±0.5µS cm−1). On Herschel Island EC increased from upstream to downstream in Ice Creek West, whereas it varied less in Ice Creek East (Fig. 5e). At Cape Bounty, river samples exhibit no clear visible trends, except for Robin Creek, where an active retrogressive thaw slump is hydrologically connected to the stream. Here, we observed a substantial downstream decrease in EC.

Figure 5Stream transects showing values of (a) dissolved organic carbon (DOC) concentration (mg L−1), (b) absorption of coloured dissolved organic matter (cDOM) at 350 nm, acDOM350 (m−1), (c) specific ultraviolet absorbance SUVA (L mg−1 m−1), (d) cDOM Slope S275–295 (10−3 nm−1), (e) electrical conductivity (EC) for streams on Herschel Island (left) and Cape Bounty (right). Note that Ice Creek West (ICW) and Ice Creek East (ICE) on Herschel Island were sampled at different dates as indicated in the legend.

At Herschel Island, overall, DOC concentration and acDOM350 (Fig. 5a) decreased downstream at all sampling periods. However, we observed a stronger decrease in Ice Creek West compared to Ice Creek East. Ice Creek West shows an increase in DOC and acDOM350 concentration at ∼1300 m, where a tributary joins the main stem. On 30 July 2016, when both streams on Herschel Island were sampled simultaneously, Ice Creek East showed significantly lower (p<0.05) DOC concentrations and acDOM350 than Ice Creek West throughout the entire profile. At Cape Bounty, DOC concentrations and acDOM350 do not show clear downstream trends but are at a rather low <2 mg L−1 in all streams. One exception is West River after the rainfall event on 12 August 2017, where we found a slight downstream increase in DOC and acDOM350 with generally higher levels of DOC and acDOM350 compared to the period before rainfall. In Robin Creek, we observed an increase in DOC from 1.3 to 1.7 mg L−1 as the stream is impacted by a retrogressive thaw slump, and then DOC decreases thereafter. Boundary River shows similar concentrations to Robin Creek. Generally, for Cape Bounty rivers other than West River the number of samples is likely too low to allow clear statements about downstream trends.

SUVA values on Herschel Island showed a similar, however weaker, decreasing downstream trend DOC and acDOM350 (Fig. 5c). Different from DOC and acDOM350, SUVA values of both Herschel Island streams are very similar. In contrast, at Cape Bounty, West River (sampled after rainfall) showed higher SUVA than the remaining rivers (sampled before the rainfall).

Spectral slope values (S275–295) at Herschel Island showed an increase downstream (Fig. 5d). When sampled on the same day, Ice Creek West showed only slightly smaller spectral slope values along the stream profile compared to Ice Creek East on 30 July 2016. Significant differences were observed between different sampling periods in Ice Creek West. Spectral slope values were smallest after the first rainfall event (19 July) and increase progressively over the course of the season. The Cape Bounty rivers showed the highest spectral slopes for East River ($\mathrm{16.1}±\mathrm{1.6}×{\mathrm{10}}^{-\mathrm{3}}$ nm−1) and the lowest for West River ($\mathrm{11.9}±\mathrm{0.8}×{\mathrm{10}}^{-\mathrm{3}}$ nm−1). A slight downstream increase in spectral slope was recorded in West River.

4.3 Temporal changes of DOM under different meteorological conditions

The mean annual air temperature on Herschel Island was −6.3C in 2016 with mean temperatures of 9.4 C in July and 7.7 C in August. During the monitoring period, rainfall events of 33.9 mm (19 July), 9.3 mm (30 July), and 12.7 mm (5 August) were recorded. Cape Bounty had a mean annual air temperature of −15.3C in 2017, with mean air temperatures of 4.5 C in July and 1.6 C in August. During the monitoring period, two rainfall events of 0.2 mm (4 August) and 1.2 mm (8 August) occurred.

Changes in discharge, DOM composition, and conductivity over the summer season were observed for both streams at Herschel Island. Rainfall response is direct with steep rising hydrographs and elongated falling limbs (Fig. 6a) in both streams (detailed presentation of rainfall response in Coch et al., 2018). In both streams, DOC, acDOM350, and SUVA were highest following the 33.9 mm rainfall event (Event 1). Following the rainfall event, declining DOC was accompanied by a decline in acDOM350 and SUVA and an increase in S275–295 (Fig. 6b–e). EC is steadily increasing after peak flow in both streams (Fig. 6f).

Figure 6Time series from Herschel Island in 2016 showing (a) discharge (m3 s−1) and hourly rainfall (mm) from Ice Creek West, (b) dissolved organic carbon (DOC) concentration (mg L−1), (c) coloured dissolved organic matter absorption at 350 nm acDOM350 (m−1), (d) specific ultraviolet absorbance SUVA (L mg−1 m−1), and (e) the cDOM slope S275–295 (10−3 nm−1) over the summer season in 2016 for Ice Creek West (magenta) and Ice Creek East (green) respectively. The onset of rainfall events is marked with vertical blue lines. As described in the methods, DOC concentrations were corrected between 30 July and 7 August. The scale depicts only S275–295 values below $\mathrm{18}×{\mathrm{10}}^{-\mathrm{3}}$ nm−1 to capture the variability; hence the two outliers in Ice Creek East (Fig. 2) are not displayed.

The subsequent rainfall event (Event 2, 9.3 mm) led to an increase in DOC, acDOM350, and S275–295 and a decrease in SUVA (Fig. 6b–e) in Ice Creek West on 30 July. This dynamic was only captured to some extent in Ice Creek East, which was sampled at a longer time interval. Base flow increased after this rainfall event (Fig. 6a). EC in both streams drops with peak flow and increases thereafter.

The hydrochemical response to rainfall Event 3 (12.7 mm) was different from the response to Event 2. An initial decrease in DOC, acDOM350, and S275–295 is followed by a sharp increase in these parameters in Ice Creek West. SUVA shows two peaks on 3 August followed by a general decreasing trend until the end of our sampling period 10 August. Ice Creek East had a different response, showing an increase in DOC and acDOM350 and a distinct decrease in SUVA.

5 Discussion

5.1 Catchment processes and DOM alteration

5.1.1 Regional catchment properties of DOM

Our study sites show strong differences in DOM quantity and quality related to their geographic location and environmental setting. DOM characteristics, such as SUVA and S275–295, provide insights into potential sources, degradation state, and properties of the water bodies. Herschel Island (low Arctic) shows on average significantly higher values in DOC, acDOM350, and SUVA than Cape Bounty (high Arctic). Catchment topography (Connolly et al., 2018), vegetation type, and soil characteristics (Harms et al., 2016) are important drivers of DOC concentrations in catchments. The greater abundance of vegetation in the low Arctic (Fig. 1) delivers more organic material resulting in high amounts of lignin introduced into the aquatic system (Sulzberger and Durisch-Kaiser, 2009) compared to the high Arctic. The Herschel Island tundra catchments have thick organic moss mats and a dense layer of vascular plant cover, which delivers plant detritus that is continuously decomposed resulting in DOC and acDOM350 values higher than in the high Arctic catchment at Cape Bounty.

In addition to higher DOM concentrations, we observed that the sampled surface waters in the Herschel Island catchments contain DOM with specific optical characteristics for high aromaticity and high molecular weight (Guéguen et al., 2007; Guo et al., 2007) that are indicative of fresh organic matter (Neff et al., 2006; Stedmon et al., 2011). High value ranges and variability of SUVA and S275–295 in the surface waters of the Cape Bounty catchments point towards a broad spectra of different DOM sources and quality. The sampled surface waters throughout the 2017 summer season contain DOM with high spectral slope values (S275–295) and low SUVA as well as surface waters with low S275–295 and high SUVA. This is due to different flow pathways, residence time, and permafrost disturbance, delivering DOM of different quality.

At Cape Bounty two different water types were identified based on the acDOM350-to-DOC ratios (Fig. 3c). The group of surface waters with lower acDOM350-to-DOC ratios is dominated by standing water bodies. High residence times in standing waters make photodegradation of DOM a dominant process (Vonk et al., 2015b) and result in an increase in S275–295 and decrease in the cDOM–DOC ratio. The group of surface waters with higher acDOM350-to-DOC ratios, higher SUVA, and lower S275–295 are, in contrast, dominated by flowing water. Higher turbidity in flowing waters potentially limits photodegradation processes (Cory et al., 2015, 2014) preserving low S275–295. Within the catchments, there may be also more import of fresh organic material to the flowing water bodies.

5.1.2 Downstream patterns of DOM and impact of permafrost disturbance

Transport and degradation of DOM is a dynamic process. Vonk et al. (2015b) showed that the microbial and photo-degradability decreased from small streams towards larger rivers within the continuous permafrost zone. The fate of DOM along lateral flow pathways from headwater streams through lakes and large rivers to the ocean is altered by photochemical and biological oxidation (Cory et al., 2015, 2014). Studies show the importance of headwater systems where photodegradation (Cory et al., 2014) and bacterial respiration of ancient permafrost-derived DOC are prevalent (Mann et al., 2015). Our sampling strategy along the rivers in combination with detailed mapping of the catchments with a focus on permafrost disturbances provide insights into upstream to downstream patterns in small coastal catchments in both the low and high Arctic.

At Herschel Island, we found a high variability of DOC, SUVA, and S275–295 in the headwaters of Ice Creek West. The locations at 2000 and 1300 m from the outflow show distinct high values of DOC and acDOM350 compared to the other locations downstream of them. These high concentrations are a result of degrading ice-wedge polygons, which heavily influence DOM in the headwaters of the stream (Coch et al., 2020). The location at 1300 m marks the inflow of another headwater tributary impacted by degrading ice-wedge polygons. Thus, the main expected sources for fresh mobilized DOM, from deeper permafrost soil horizons, are the headwaters and tributary water. This is supported by high SUVA and low S275–295 indicating high molecular weight.

DOC and acDOM350 are highest in the headwaters and decrease downstream. Combined with increasing S275–295 along both streams, our results are indicative of a progressive photochemical degradation of DOM. S275–295 has been found to be a good indicator for photodegradation of DOM (Fichot and Benner, 2012; Fichot et al., 2013; Helms et al., 2008) and also been observed along a flow path continuum of the Kolyma River (Frey et al., 2016). They found a relatively constant proportion of bioavailable DOM along the entire flow path, indicating an acclimatization of aquatic microorganisms to downstream DOM changes and/or the generation of labile DOM for microbial processing through photodegradation. Cory et al. (2014, 2015) show at a subarctic site that DOC in headwater streams, which are directly sourced by soil water, have low prior exposure to light and are therefore prone to photodegradation to CO2.

At Cape Bounty, optical data of downstream patterns are more limited (see Sect. 5.3). West River shows an increase in DOC concentration downstream (3 August 2017), which is also reflected in an increase in acDOM350. As discussed by Fouché et al. (2017) and Wang et al. (2018), the West River is characterized by a downstream increase in autochthonous DOM. SUVA and S275–295 do not show strong differences downstream in the West River, suggesting little modification of DOM through microbial and/or photodegradation processes. A retrogressive thaw slump at Robin Creek heavily impacted DOM quality. At ∼2100 m distance from the outflow, closest to the slump, we see the highest DOC, S275–295, and EC values and lowest SUVA, indicative of low aromaticity and lower molecular weight.

Abbott et al. (2014) found that DOM is most biodegradable during active disturbance at sites in the subarctic. SUVA values at thermokarst outflows in that study were half of the magnitude of the undisturbed reference waters, indicating less aromatic DOC. High S275–295 and SR were observed in conjunction with geomorphic disturbance in headwater streams of West River by Fouché et al. (2017). Impact of retrogressive thaw slumps on DOM quality was also studied in the subarctic Peel Plateau by Littlefair et al. (2017). They reported dynamics at modestly sized slumps similar to those we observed at Robin Creek: DOC concentration is highest directly at the slump outflow and is lower downstream compared to undisturbed upstream conditions. The authors attribute low SUVA and high S275–295 within the disturbed site to deep permafrost flow pathways. SUVA values of surface waters in Lake Bounty catchments within the slump and downstream are very similar to the ones reported by Littlefair et al. (2017). Overall, DOM characteristics in both study areas are affected by local permafrost disturbances. In sampling transects which are not affected by permafrost disturbances, gradual degradation was observed.

5.1.3 Rainfall event impacts on DOM

Rain magnitude, intensity, and antecedent conditions play an important role in mobilizing DOM in permafrost catchments. At Herschel Island, we captured the response to three different rainfall events through continuous sampling at the outflow and repeated sampling along the longitudinal stream profile in Ice Creek West.

Rainfall Event 1 (33.9 mm) was captured only at the receding hydrograph at the outflow (Fig. 6), but along the stream profile in Ice Creek West (Fig. 5). This event of high magnitude and intensity led to high SUVA and low S275–295 values indicating fresh plant-derived DOM that is prone to degradation processes, both microbial and photodegradation. After this event, the hydrograph recedes, and the DOM signature during the “post-rain” conditions suggests a source from deeper in the active layer that contains potentially older carbon (decreased SUVA and increased S275–295 at the outflow and throughout the profile) than surface soils with mostly recently fixed carbon from the vegetation cover. The contrasting response of Ice Creek West to rainfall Events 2 and 3 suggests different sourcing of DOM and controlling factors.

During the second rainfall event (9.3 mm), as DOC concentration increased, we found a decrease in SUVA accompanied by an increase in S275–295. This indicates a decrease in aromaticity and a lower molecular weight, indicative of more decomposed material. The following event of 12.7 mm led to a decrease in DOC concentration and S275–295 and an increase in SUVA, indicating an increase in aromaticity and a higher molecular weight – suggesting fresher and lignin-rich plant-derived DOM. A change in water sources for these two rainfall events was examined by Coch et al. (2020). Whereas runoff during the 9.3 mm rainfall event showed the signature of supra-permafrost water, which was forced out during that rainfall event, runoff during the subsequent 12.7 mm rainfall event reflected the isotopic signature of rain (Coch et al., 2020). Thus, the DOM was first sourced from the surface and through the entire active layer and had a longer residence time than the rain event after.

These results indicate that antecedent (pre-rainfall) soil water conditions play a crucial role in the sourcing of DOM. The second rainfall event (9.3 mm) occurred about 10 d later, whereas the time difference between the second and the third one was less than 4 d (i.e. the soil was saturated, mobilizing surface OM). In addition to the antecedent conditions, the magnitude and intensity of the rainfall event might also play an important role here. The 9.3 mm rainfall event occurred over a period of 3 d. Thus, the flow pathways during this event might be deeper in the active layer, mobilizing more decomposed OM (Marín-Spiotta et al., 2014). In contrast, the subsequent 12.7 mm event occurred within 1 d, which presumably led to increased overland flow and the mobilization of surface OM. Base flow in this catchment increases with summer rainfall and as the summer season progresses (Coch et al., 2018). The authors also reported a linear increase in DOC export with increasing runoff. Our dataset shows that the quality of exported DOC depends on the intensity of rainfall and the antecedent conditions, which in turn determine hydrological flow pathways and sourcing of DOM.

When sampling before and after rainfall on Cape Bounty (West River), we found a substantial increase in DOC concentration compared to the pre-rainfall concentrations. Fouché et al. (2017) conducted an extensive study of DOM quality in four headwater streams of West River (Cape Bounty) and also reported an increase in DOC concentrations and fluxes during storm flow. They observed a change in DOM quality: enrichment in fresh low-molecular-weight (LMW), microbially derived, components as indicated by an increase in S275–295 and a decrease in SUVA during rainfall. Although we do not have data on the optical properties for West River before the rainfall event, similar concentrations of DOC in the West and East rivers point towards similar optical characteristics at that time. Base flow in undisturbed high Arctic headwater streams seems therefore characterized by more high-molecular-weight (HMW) humic-like components with high aromaticity (low spectral slope and increase in SUVA) relative to storm flow DOM. In turn, storm flow leads to an export of DOM characterized by lower molecular weight and decreased aromaticity (high spectral slope, decreased SUVA). Fouché et al. (2017) explain this pattern by a change in flow pathways from shallow active layer soils (base flow) to subsurface runoff (rainfall), where soluble components from mineral soils deeper in the active layer are mobilized. Associated with the change in DOM quality, they also found an increase in total dissolved solids (TDSs) supporting this hypothesis. Impacts of changing flow pathways on DOM quality are also reported from a subarctic setting by Balcarczyk et al. (2009). The increased residence time of percolating water through the active layer leads to a selective sorption of compounds to mineral soil particles. The authors describe that hydrophobic compounds are absorbed, while hydrophilic compounds remain in the solution, and are therefore exported from the catchment (Balcarczyk et al., 2009). Further, an increased residence time and subsurface flow mobilizes DOC that is more microbially degraded (Striegl et al., 2005; Ward and Cory, 2015).

Several studies anticipate a shift towards deeper flow pathways as active layer depths increase with climate change (Drake et al., 2018; Liljedahl et al., 2016; Mann et al., 2015; O'Donnell et al., 2014; Ward and Cory, 2015). These studies found that permafrost-derived DOM is more labile and thus easily used by bacteria compared to surface (organic mat) DOM. As described above, we show that different flow pathways are activated during storm flow conditions at the low and high Arctic locations, which influences the quality of DOM exported. At the low Arctic setting our data suggest that more permafrost-derived DOM is exported with increasing base flow during the season and during a rainfall event of smaller magnitude and lower intensity. Based on the optical properties, this material shows low molecular weight and aromaticity (i.e. it is already altered). In contrast, high-magnitude and high-intensity rainfall events that act on saturated soil lead to shorter residence time in the flow path and thus export more fresh (less altered due to different degradation processes), near-surface-derived DOM (higher SUVA and lower S275–295). As summer rainfall is projected to increase across the Arctic (Bintanja, 2018; Bintanja and Andry, 2017), an increase in DOC export is expected (Coch et al., 2018). Small catchments in the subarctic Canadian Shield have already shifted towards a nival–pluvial flow regime, leading to substantial increases in organic matter fluxes during fall and winter (Spence et al., 2011, 2015). The DOM quality will depend on the residence time and thus flow pathways within the catchment, which in turn are controlled by the frequency and magnitude of the rainfall events and the thaw depth of the active layer.

5.2 DOM dynamics of small and large Arctic catchments

The knowledge of ecosystem responses to external disturbances is necessary for predictive models. Flux of DOM is a significant input into Arctic coastal oligotrophic marine environments. The major Arctic catchments cover approximately half of the Arctic drainage basin, whereas the remainder is covered by the complex network of smaller catchments. In the Arctic, most historical data and studies on riverine DOM dynamics are from the major Arctic catchments. Research on small catchments may yield different information to riverine DOM of the major Arctic catchments.

We linked DOC and acDOM350 from this study and the literature (Tables 3 and S1 and S2 in the Supplement) to latitude and the soil organic carbon content (SOCC) at 0–30 and 0–100 cm depth as retrieved from Hugelius et al. (2013). We found a positive correlation ($\mathit{\rho }=\mathrm{0.53}/\mathrm{0.51}$, p<0.05) between SOCC and DOC concentration, indicating that vegetation coverage and the connected SOCC are influencing DOM. The relationship between acDOM350 and SOCC is also significant, although weaker ($\mathit{\rho }=\mathrm{0.26}/\mathrm{0.34}$, p<0.05). It is important to bear in mind that the northern circumpolar soil carbon database is a product of upscaling and will most likely not cover the spatial variability reported in the studies. Nevertheless, the data show a decrease in SOCC at higher latitudes, influenced by climate and accordingly vegetation cover and related soil cover. DOM and SOCC are further influenced by watershed topography. Longer residence times in low-relief terrain and high hydrologic connectivity facilitate leaching and export of DOM from soil organic matter (Connolly et al., 2018; Harms et al., 2016). However, DOM characteristics are not always influenced by the subcatchments of the investigated waterbodies. Several of the studies contain data on waterbodies with small subcatchments but located within the large Arctic river floodplains with very large catchment sizes. For example, Dvornikov et al. (2018) report untypically high cDOM for tundra lakes in Yamal (western Siberia) and Dvornikov et al. (2018) for tundra lakes in the Lena delta (central Siberia), which are all influenced by the spring flood of the large rivers seasonally flashing high amounts of organic material into these lakes.

Table 3Correlation matrix using the Spearman's ρ correlation coefficient between latitude, dissolved organic carbon concentration (DOC), coloured dissolved organic carbon absorption at 350 nm (acDOM350), soil organic carbon content (SOCC) in 0–30 and 0–100 cm depth (Hugelius et al., 2014). Significance levels of p<0.05 and p≤0.01 (in bold) are indicated.

Strong positive correlations between DOC and acDOM350 found in this study (Fig. 3a) are also characteristic for riverine DOM of the large Arctic rivers (Walker et al., 2013), across the land–ocean continuum in the eastern Arctic (Juhls et al., 2019; Mann et al., 2016) and globally (Massicotte et al., 2017). We used DOC and cDOM data available from surface waters in northeastern Canada (Breton et al., 2009), Scandinavia (Forsström et al., 2015; Kellerman et al., 2015), and Alaska (Cory et al., 2015; Larouche et al., 2015). Comparing our data from the low Arctic and high Arctic sites to those found in the literature confirms the strong positive relationship (ρ=0.85, p<0.05) between DOC and acDOM350 (Fig. S3), indicating the robustness for using the optical parameter acDOM350 as a proxy for DOC concentration in terrestrial freshwater systems. Compared to DOC concentrations and cDOM magnitudes of other small Arctic catchments, some of the samples from Cape Bounty show extremely low values, which reflects the low supply of organic matter due to the low plant abundance in the high Arctic. DOC and cDOM from Herschel Island are within the range of most studies on Arctic catchments whereas studies at high Arctic sites with low vegetation cover are underrepresented. However, compared to other reported studies with DOC and cDOM in the same value range, acDOM350 is slightly depleted, which is visible in a lower cDOM-to-DOC ratio. This can be a result of stronger photodegradation compared to other sites with eventually more turbid water types, but information on turbidity or suspended matter is frequently not provided in those studies focusing on DOM dynamics and properties.

Due to snowmelt dynamics and active layer development throughout the summer season in Arctic catchments, there is also a strong seasonal influence on the DOC-to-cDOM relationship leading to variability throughout the season and regions (Mannino et al., 2008; Vantrepotte et al., 2015). Walker et al. (2013) report SUVA for three different flow regimes of the large Arctic rivers: peak flow (spring freshet), midflow (summer), and base flow (winter). The SUVA values reported in this study (2.9±0.4 L mg−1 m−1 for Herschel Island and 2.8±1.1 L mg−1 m−1 for Cape Bounty) are higher than the mean mid-flow SUVA for the five large Arctic rivers (2.4 L mg−1 m−1), which ranges between 2.0 L mg−1 m−1 in the Mackenzie River and 2.7 L mg−1 m−1 in the Ob'. This supports the model proposed by Vonk et al. (2015b), in which DOM exported from smaller rivers has a higher aromaticity, which suggests that the material is fresh and less altered by different degradation processes. However, our results also show that a broad range of SUVA as well as S275–295 can be found in small Arctic catchments. This highlights the importance of small rivers and streams for the magnitude of potential modification of DOM before waters are exported in the large rivers or directly to the Arctic Ocean. The ratio of SUVA versus S275–295 values from the large Arctic rivers (Fig. 4a) falls in the same value range as the ratio from the catchments of Cape Bounty and Herschel Island. However, SUVA values of the large Arctic rivers are in the low range only, reflecting the higher degradation status of the transported DOM. The results of this study suggest that small Arctic catchments potentially deliver material that is fresher and more prone to degradation compared to DOM of the large Arctic rivers.

5.3 Limitations of cDOM measurements from terrestrial sources

There are some constraints to optical DOM measurements and the nature of the samples themselves that we encountered in this study. As described in the methods section, some samples formed precipitates inside the bottles in the form of small thin flakes, which partly remained in suspension or accumulated at the bottom. All samples were filtered in the field through 0.7 µm glass fibre filters, and the precipitation occurred after filtration during storage. At Cape Bounty, these problematic samples had very high acDOM350 values of 13.9±13.8 m−1 with a maximum of 75.8 m−1 and SUVA values of 10.1±11.5 L mg−1 m−1 with a maximum of 59.5 L mg−1 m−1. Those values are significantly higher (p<0.05) than the mean values reported in Table 2 and are not realistic for natural surface waters. At Herschel Island, acDOM350 and SUVA did not differ significantly from the mean (11.8±0.8 m−1 and 3.5±0.4 L mg−1 m−1 respectively).

As described in the methods section (Sect. 3.1), samples showing precipitates in the laboratory were excluded from the study, even if the absorption values were plausible when compared to the corresponding DOC concentration (Fig. 7). At Cape Bounty, this was the case for 25 out of 55 samples. We assume that absorbance measurements are high because of scattering by newly formed colloid complexes and precipitates and absorbance from other absorbing dissolved constituents. Also, Hansen et al. (2016) and Weishaar et al. (2003) report that SUVA values above 6.0 L mg−1 m−1 are indicative for an optical disturbance due to other constituents in the sample (Hansen et al., 2016; Weishaar et al., 2003).

Figure 7Relationship between coloured dissolved organic matter absorption acDOM350 (m−1) and dissolved organic matter concentration DOC (mg L−1) at Herschel Island (HE) displayed as triangle and Cape Bounty (CB) shown as circle. Samples marked in orange were excluded from the study due to flocculation after filtration (Sect. 3.1). The samples circled in black show disproportionately high absorption values in relation to the DOC concentration.

We suggest that the optical interference could be due to polymeric iron (hydr)oxides or high concentrations of dissolved iron and changing pH conditions of the sample. Dissolved iron in terrestrially dominated waters is dominantly complexed with humic and fulvic acids. With changing temperature and changing pH of the sample filtrates, redox reaction can result in colloid formation and phase changes, which then strongly affect the optical properties of the sample filtrate by scattering. Poulin et al. (2014) describe how iron (Fe(II,III)) is known to interfere with the absorption of cDOM with a linear dependency of increasing acDOM with increasing Fe(III) concentration in the water. Poulin et al. (2014) suggest correcting cDOM absorption coefficients according to the iron concentrations using correction coefficients. Coch et al. (2018) report total aqueous dissolved iron concentrations from Herschel Island. High total iron concentration is found to occur in high acDOM350 (Fig. S2), which indicates a potential influence of iron concentration on the absorption. Fraction of Fe(II) and Fe(III) on the total iron concentration was not measured as a standard hydrochemistry measurement; thus the correction could not be performed. However, Fig. 7 clearly shows that the samples that were removed fell into the problematic group (circled), where cDOM was overestimated compared to DOC concentration. This conservative approach also removed other samples with reasonable cDOM-to-DOC ratios.

Poulin et al. (2014) also showed that in samples with low pH the dominant fraction of iron is Fe(II), which then potentially can precipitate as Fe(III) with increasing pH during transport and storage. The Cape Bounty samples that showed a substantially lower pH, likely caused by low vegetation, are therefore more prone to precipitate Fe(III) colloids that affect the optical absorption measurements and lead to the high absorption values at 700 nm (Fig. S4). Herschel Island samples originally already had a higher pH compared to Cape Bounty. Thus, we expect that the dominant fraction of iron on Herschel Island was Fe(III), which leads to a lower potential of Fe(III) precipitation compared to Cape Bounty. Catchment properties that influence riverine pH such as the local lithology may play an important role. In the case of alpine and high Arctic catchments with thin or no soil cover, a bedrock composition of acid rocks in the catchments will lead to lower pH values in surface waters as is the case for Cape Bounty. Conversely, surface waters from Herschel Island catchments on glacial moraines and marine sediments are characterized by higher alkalinity.

6 Conclusion

This study investigates DOM optical properties in low and high Arctic surface water environments and downstream patterns with regard to permafrost disturbance and rainfall events. We find that both Arctic locations exhibit a distinct signature of DOC concentration and acDOM350 linked to the differences in vegetation cover and SOCC content. Compared to the high Arctic (Cape Bounty), DOC and cDOM in the low Arctic (Herschel Island) are higher due to the greater abundance of plant material and higher SOCC. In both regions, the strong terrestrial signature of DOM is apparent in the optical properties, which is typical for small headwater catchments. The relationship between acDOM350 and DOC is very strong across both regions and including data from the literature, proving the applicability of cDOM as a tracer for DOC throughout different aquatic Arctic environments (rivers, streams, and lakes).

Comparing DOM optical characteristics (SUVA, S275–295, and SR) from large Arctic rivers to the surface waters in our study, we find that the low and high Arctic small catchments potentially deliver fresh, less altered DOM. However, the results also show that DOM characteristics indicate organic matter modification and degradation cover a broad spectrum and can be highly variable in small catchments in space (along longitudinal transects) and time (throughout the season). Degrading ice-wedge polygons and retrogressive thaw slumps impact DOM quantity and quality in the catchments. This underlines the importance of small catchments for potential DOM modification and degradation before waters enter bigger rivers or coastal waters of the Arctic Ocean.

The optical characteristics of DOM prove to be useful for assessing downstream patterns in the studied streams. The downstream increase in S275–295 is indicative for photodegradation processes, which is apparent in most of the streams. Although the temporal resolution of data at Cape Bounty is limited, we found a similar response to rainfall events as in the Herschel Island study. Rainfall leading to runoff with a short residence time (rainfall of high magnitude and intensity, dry antecedent conditions in the catchment) leads to the export of fresh near-surface-derived DOM (higher SUVA, lower S275–295). In contrast, base flow conditions and long residence times (including low-magnitude rainfall events and a saturated catchment) favours the export of permafrost-derived DOM that has undergone microbial processing in the soil. Examining flow pathways and residence time will be crucial to assess the impacts of projected increasing summer rainfall across the Arctic.

Monitoring optical properties of DOM in combination with a mapping of permafrost disturbances across river catchments will be a useful tool for assessing DOM fluxes and DOM quality changes at a pan-Arctic scale.

Data availability
Data availability.

Data have been made available through PANGAEA (https://doi.pangaea.de/10.1594/PANGAEA.897289; Coch et al., 2019).

Supplement
Supplement.

Author contributions
Author contributions.

CC, HL, and SFL developed the study design. Field work was conducted by CC in 2016 and by CC, SFL, and MJL in 2017. MF partly funded and supervised laboratory analyses for cDOM measurements. BJ processed the absorbance spectra and contributed to developing the paper. CC ran laboratory analyses and processed and interpreted the data with input from BH, MJL, SFL, and HL and prepared the paper with editorial contributions from all co-authors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

Thanks are due to the two anonymous reviewers, who helped to improve the paper. We are grateful to the Yukon Territorial government, Yukon Parks (Herschel Island Qikiqtaryuk Territorial Park), and the Aurora Research Institute for their support during this project. Caroline Coch received a scholarship and travel support from the Studienstiftung des deutschen Volkes. Research at CBAWO is supported by the Canadian Natural Sciences and Engineering Research Council (NSERC) and ArcticNet National Centres of Excellence. Polar Continental Shelf Program provided field logistical support. We thank the hamlet of Resolute and the Nunavut Research Institute for their support. The authors wish to thank Antje Eulenburg, Christian Knoblauch, Birgit Grabellus, Justus Gimsa, Marek Jaskólski, Jennifer Krutzke, Nicole Mätzing, Paul Overduin, Julian Schneider, and Samuel Stettner for their help in the instrumentation setup, data collection in the field, and laboratory analyses. We acknowledge Saskia Foerster and Sabine Chabrillat (GFZ Potsdam) for providing access to the laboratory spectrometer. A special thanks to Cameron Eckert, Richard Gordon, Ricky Joe, Paden Lennie, Edward McLeod, and Samuel McLeod for their support and helpful insights in the field. Many thanks also to the CBAWO 2017 field team.

Financial support
Financial support.

This research has been supported by Horizon 2020 (Nunataryuk (grant no. 773421)), the Geo.X Research Network for Geosciences in Berlin and Potsdam (grant no. SO_087_GeoX), and the Helmholtz Association (grant no. VH-NG-801 to Hugues Lantuit).

The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.

Review statement
Review statement.

This paper was edited by Tom J. Battin and reviewed by two anonymous referees.

References

Abbott, B. W., Larouche, J. R., Jones, J. B., Bowden, W. B., and Balser, A. W.: Elevated dissolved organic carbon biodegradability from thawing and collapsing permafrost, J. Geophys. Res.-Biogeo., 119, 2049–2063, https://doi.org/10.1002/2014jg002678, 2014.

ADAPT: Carbon, nitrogen and water content of the active layer from sites across the Canadian Arctic Nordicana D21, Nordicana, Québec, Canada, https://doi.org/10.5885/45327AD-5245D08606AB4F52, 2014.

Aiken, G. R.: Fluorescence and dissolved organic matter: A chemist's perspective: Chapter 2, 35–74, https://doi.org/10.1017/CBO9781139045452.005, Cambridge University Press, 2014.

Anderson, N. J. and Stedmon, C. A.: The effect of evapoconcentration on dissolved organic carbon concentration and quality in lakes of SW Greenland, Freshwater Biol., 52, 280–289, https://doi.org/10.1111/j.1365-2427.2006.01688.x, 2007.

Balcarczyk, K. L., Jones, J. B., Jaffé, R., and Maie, N.: Stream dissolved organic matter bioavailability and composition in watersheds underlain with discontinuous permafrost, Biogeochemistry, 94, 255–270, https://doi.org/10.1007/s10533-009-9324-x, 2009.

Bintanja, R.: The impact of Arctic warming on increased rainfall, Sci. Rep.-UK, 8, 16001, https://doi.org/10.1038/s41598-018-34450-3, 2018.

Bintanja, R. and Andry, O.: Towards a rain-dominated Arctic, Nat. Clim. Change, 7, 263–267, https://doi.org/10.1038/nclimate3240, 2017.

Breton, J., Prairie, Y., Vallières, C., and Laurion, I.: Limnological properties of permafrost thaw ponds in northeastern Canada, Can. J. Fish. Aquat. Sci., 66, 1635–1648, https://doi.org/10.1139/f09-108, 2009.

Burn, C. R.: Herschel Island Qikiqtaryuk: A Natural and Cultural History of Yukon's Arctic Island, edited by: Burn, C. R., University of Calgary Press, 242 pp., 2012.

CAVM: Circumpolar Arctic Vegetation Map, Conservation of Arctic Flora and Fauna (CAFF) Map, 2003.

Coch, C., Lamoureux, S. F., Knoblauch, C., Eischeid, I., Fritz, M., Obu, J., and Lantuit, H.: Summer rainfall DOC, solute and sediment fluxes in a small Arctic coastal catchment on Herschel Island (Yukon Territory, Canada), Arctic Science, 4, 750–780, https://doi.org/10.1139/as-2018-0010, 2018.

Coch, C., Juhls, B., Lamoureux, S., Lafrenière, M., Fritz, M., Heim, B., and Lantuit, H.: Colored dissolved organic matter (cDOM) absorption measurements in terrestrial waters on Herschel Island (Low Arctic) and Melville Island (High Arctic) in 2016 and 2017, PANGAEA, https://doi.org/10.1594/PANGAEA.897289, 2019.

Coch, C., Ramage, J. L., Lamoureux, S., Knoblauch, C., Meyer, H., and Lantuit, H.: Spatial variability of dissolved organic carbon (DOC), solutes and suspended sediment in disturbed Low Arctic coastal watersheds, J. Geophys. Res.-Biogeo., in review, 2020.

Connolly, C. T., Khosh, M. S., Burkart, G. A., Douglas, T. A., Holmes, R. M., Jacobson, A. D., Tank, S. E., and McClelland, J. W.: Watershed slope as a predictor of fluvial dissolved organic matter and nitrate concentrations across geographical space and catchment size in the Arctic, Environ. Res. Lett., 13, 104015, https://doi.org/10.1088/1748-9326/aae35d, 2018.

Cory, R. M., Ward, C. P., Crump, B. C., and Kling, G. W.: Carbon cycle. Sunlight controls water column processing of carbon in arctic fresh waters, Science, 345, 925–928, https://doi.org/10.1126/science.1253119, 2014.

Cory, R. M., Harrold, K. H., Neilson, B. T., and Kling, G. W.: Controls on dissolved organic matter (DOM) degradation in a headwater stream: the influence of photochemical and hydrological conditions in determining light-limitation or substrate-limitation of photo-degradation, Biogeosciences, 12, 6669–6685, https://doi.org/10.5194/bg-12-6669-2015, 2015.

Drake, T. W., Tank, S. E., Zhulidov, A. V., Holmes, R. M., Gurtovaya, T., and Spencer, R. G. M.: Increasing Alkalinity Export from Large Russian Arctic Rivers, Environ. Sci. Technol., 52, 8302–8308, https://doi.org/10.1021/acs.est.8b01051, 2018.

Dvornikov, Y., Leibman, M., Heim, B., Bartsch, A., Herzschuh, U., Skorospekhova, T., Fedorova, I., Khomutov, A., Widhalm, B., Gubarkov, A., and Rößler, S.: Terrestrial CDOM in Lakes of Yamal Peninsula: Connection to Lake and Lake Catchment Properties, Remote Sensing, 10, 167, https://doi.org/10.3390/rs10020167, 2018.

Edwards, R. and Treitz, P.: Vegetation Greening Trends at Two Sites in the Canadian Arctic: 1984–2015, Arct. Antarct. Alp. Res., 49, 601–619, https://doi.org/10.1657/aaar0016-075, 2018.

Environment and Climate Change Canada, Historical Data: available at: http://climate.weather.gc.ca/historical_data/search_historic_data_e.html, last access: 12 December 2018.

Fichot, C. G. and Benner, R.: The spectral slope coefficient of chromophoric dissolved organic matter (S275–295) as a tracer of terrigenous dissolved organic carbon in river-influenced ocean margins, Limnol. Oceanogr., 57, 1453–1466, https://doi.org/10.4319/lo.2012.57.5.1453, 2012.

Fichot, C. G., Kaiser, K., Hooker, S. B., Amon, R. M., Babin, M., Belanger, S., Walker, S. A., and Benner, R.: Pan-Arctic distributions of continental runoff in the Arctic Ocean, Sci. Rep.-UK, 3, 1053, https://doi.org/10.1038/srep01053, 2013.

Forsström, L., Rautio, M., Cusson, M., Sorvari, S., Albert, R.-L., Kumagai, M., and Korhola, A.: Dissolved organic matter concentration, optical parameters and attenuation of solar radiation in high-latitude lakes across three vegetation zones, Ecoscience, 22, 17–31, https://doi.org/10.1080/11956860.2015.1047137, 2015.

Fouché, J., Lafrenière, M. J., Rutherford, K., and Lamoureux, S.: Seasonal hydrology and permafrost disturbance impacts on dissolved organic matter composition in High Arctic headwater catchments, Arctic Science, 3, 378–405, https://doi.org/10.1139/as-2016-0031, 2017.

Frey, K. E. and Smith, L. C.: Amplified carbon release from vast West Siberian peatlands by 2100, Geophys. Res. Lett., 32, L09401, https://doi.org/10.1029/2004gl022025, 2005.

Frey, K. E., Sobczak, W. V., Mann, P. J., and Holmes, R. M.: Optical properties and bioavailability of dissolved organic matter along a flow-path continuum from soil pore waters to the Kolyma River mainstem, East Siberia, Biogeosciences, 13, 2279–2290, https://doi.org/10.5194/bg-13-2279-2016, 2016.

Green, S. A. and Blough, N. V.: Optical absorption and fluorescence properties of chromophoric dissolved organic matter in natural waters, Limnol. Oceanogr., 39, 1903–1916, https://doi.org/10.4319/lo.1994.39.8.1903, 1994.

Guéguen, C., Guo, L., Yamamoto-Kawai, M., and Tanaka, N.: Colored dissolved organic matter dynamics across the shelf-basin interface in the western Arctic Ocean, J. Geophys. Res., 112, C05038, https://doi.org/10.1029/2006JC003584, 2007.

Guo, L., Ping, C.-L., and Macdonald, R. W.: Mobilization pathways of organic carbon from permafrost to arctic rivers in a changing climate, Geophys. Res. Lett., 34, L13603, https://doi.org/10.1029/2007gl030689, 2007.

Hancke, K., Hovland, E. K., Volent, Z., Pettersen, R., Johnsen, G., Moline, M., and Sakshaug, E.: Optical properties of CDOM across the Polar Front in the Barents Sea: Origin, distribution and significance, J. Marine Syst., 130, 219–227, https://doi.org/10.1016/J.JMARSYS.2012.06.006, 2014.

Hansen, A. M., Kraus, T. E. C., Pellerin, B. A., Fleck, J. A., Downing, B. D., and Bergamaschi, B. A.: Optical properties of dissolved organic matter (DOM): Effects of biological and photolytic degradation, Limnol. Oceanogr., 61, 1015–1032, https://doi.org/10.1002/lno.10270, 2016.

Harms, T. K., Edmonds, J. W., Genet, H., Creed, I. F., Aldred, D., Balser, A., and Jones, J. B.: Catchment influence on nitrate and dissolved organic matter in Alaskan streams across a latitudinal gradient, J. Geophys. Res.-Biogeo., 121, 350–369, https://doi.org/10.1002/2015jg003201, 2016.

Helms, J. R., Stubbins, A., Ritchie, J. D., Minor, E. C., Kieber, D. J., and Mopper, K.: Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and photobleaching of chromophoric dissolved organic matter, Limnol. Oceanogr., 53, 955–969, https://doi.org/10.4319/lo.2008.53.3.0955, 2008.

Hodgson, D. A., Vincent, J. S., and Fyles, J. G.: Quaternary geology of central Melville Island, Northwest Territories, Geological Survey of Canada, Canada, https://doi.org/10.4095/119784, 1984.

Holmes, R. M., McClelland, J. W., Peterson, B. J., Tank, S. E., Bulygina, E., Eglinton, T. I., Gordeev, V. V., Gurtovaya, T. Y., Raymond, P. A., Repeta, D. J., Staples, R., Striegl, R. G., Zhulidov, A. V., and Zimov, S. A.: Seasonal and Annual Fluxes of Nutrients and Organic Matter from Large Rivers to the Arctic Ocean and Surrounding Seas, Estuar. Coast., 35, 369–382, https://doi.org/10.1007/s12237-011-9386-6, 2012.

Hugelius, G., Tarnocai, C., Broll, G., Canadell, J. G., Kuhry, P., and Swanson, D. K.: The Northern Circumpolar Soil Carbon Database: spatially distributed datasets of soil coverage and soil carbon storage in the northern permafrost regions, Earth Syst. Sci. Data, 5, 3–13, https://doi.org/10.5194/essd-5-3-2013, 2013.

Hugelius, G., Strauss, J., Zubrzycki, S., Harden, J. W., Schuur, E. A. G., Ping, C.-L., Schirrmeister, L., Grosse, G., Michaelson, G. J., Koven, C. D., O'Donnell, J. A., Elberling, B., Mishra, U., Camill, P., Yu, Z., Palmtag, J., and Kuhry, P.: Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps, Biogeosciences, 11, 6573–6593, https://doi.org/10.5194/bg-11-6573-2014, 2014.

Juhls, B., Overduin, P. P., Hölemann, J., Hieronymi, M., Matsuoka, A., Heim, B., and Fischer, J.: Dissolved organic matter at the fluvial–marine transition in the Laptev Sea using in situ data and ocean colour remote sensing, Biogeosciences, 16, 2693–2713, https://doi.org/10.5194/bg-16-2693-2019, 2019.

Kellerman, A. M., Kothawala, D. N., Dittmar, T., and Tranvik, L. J.: Persistence of dissolved organic matter in lakes related to its molecular characteristics, Nat. Geosci., 8, 454–457, https://doi.org/10.1038/ngeo2440, 2015.

Lacelle, D., Bjornson, J., and Lauriol, B.: Climatic and geomorphic factors affecting contemporary (1950–2004) activity of retrogressive thaw slumps on the Aklavik Plateau, Richardson Mountains, NWT, Canada, Permafrost Periglac., 21, 1–15, https://doi.org/10.1002/ppp.666, 2010.

Lafrenière, M. J. and Lamoureux, S. F.: Thermal Perturbation and Rainfall Runoff have Greater Impact on Seasonal Solute Loads than Physical Disturbance of the Active Layer, Permafrost Periglac., 24, 241–251, https://doi.org/10.1002/ppp.1784, 2013.

Lafrenière, M. J., Laurin, E., and Lamoureux, S. F.: The Impact of Snow Accumulation on the Active Layer Thermal Regime in High Arctic Soils, Vadose Zone J., 12, 1–13, https://doi.org/10.2136/vzj2012.0058, 2013.

Lamoureux, S. F. and Lafrenière, M. J.: Fluvial Impact of Extensive Active Layer Detachments, Cape Bounty, Melville Island, Canada, Arct. Antarct. Alp. Res., 41, 59–68, https://doi.org/10.1657/1938-4246(08-030)[lamoureux]2.0.co;2, 2009.

Lamoureux, S. F. and Lafrenière, M. J.: More than just snowmelt: integrated watershed science for changing climate and permafrost at the Cape Bounty Arctic Watershed Observatory, WIREs Water, 5, e1255, https://doi.org/10.1002/wat2.1255, 2017.

Lantuit, H. and Pollard, W. H.: Fifty years of coastal erosion and retrogressive thaw slump activity on Herschel Island, southern Beaufort Sea, Yukon Territory, Canada, Geomorphology, 95, 84–102, https://doi.org/10.1016/j.geomorph.2006.07.040, 2008.

Larouche, J. R., Abbott, B. W., Bowden, W. B., and Jones, J. B.: The role of watershed characteristics, permafrost thaw, and wildfire on dissolved organic carbon biodegradability and water chemistry in Arctic headwater streams, Biogeosciences, 12, 4221–4233, https://doi.org/10.5194/bg-12-4221-2015, 2015.

Le Fouest, V., Matsuoka, A., Manizza, M., Shernetsky, M., Tremblay, B., and Babin, M.: Towards an assessment of riverine dissolved organic carbon in surface waters of the western Arctic Ocean based on remote sensing and biogeochemical modeling, Biogeosciences, 15, 1335–1346, https://doi.org/10.5194/bg-15-1335-2018, 2018.

Lewis, T., Lafrenière, M. J., and Lamoureux, S. F.: Hydrochemical and sedimentary responses of paired High Arctic watersheds to unusual climate and permafrost disturbance, Cape Bounty, Melville Island, Canada, Hydrol. Process., 26, 2003–2018, https://doi.org/10.1002/hyp.8335, 2012.

Lewkowicz, A. G.: Dynamics of active-layer detachment failures, Fosheim Peninsula, Ellesmere Island, Nunavut, Canada, Permafrost Periglac., 18, 89–103, https://doi.org/10.1002/ppp.578, 2007.

Liljedahl, A. K., Boike, J., Daanen, R. P., Fedorov, A. N., Frost, G. V., Grosse, G., Hinzman, L. D., Iijma, Y., Jorgenson, J. C., Matveyeva, N., Necsoiu, M., Raynolds, M. K., Romanovsky, V. E., Schulla, J., Tape, K. D., Walker, D. A., Wilson, C. J., Yabuki, H., and Zona, D.: Pan-Arctic ice-wedge degradation in warming permafrost and its influence on tundra hydrology, Nat. Geosci., 9, 312–318, https://doi.org/10.1038/ngeo2674, 2016.

Littlefair, C. A., Tank, S. E., and Kokelj, S. V.: Retrogressive thaw slumps temper dissolved organic carbon delivery to streams of the Peel Plateau, NWT, Canada, Biogeosciences, 14, 5487–5505, https://doi.org/10.5194/bg-14-5487-2017, 2017.

Mann, P. J., Davydova, A., Zimov, N., Spencer, R. G. M., Davydov, S., Bulygina, E., Zimov, S., and Holmes, R. M.: Controls on the composition and lability of dissolved organic matter in Siberia's Kolyma River basin, J. Geophys. Res.-Biogeo., 117, G01028, https://doi.org/10.1029/2011jg001798, 2012.

Mann, P. J., Eglinton, T. I., McIntyre, C. P., Zimov, N., Davydova, A., Vonk, J. E., Holmes, R. M., and Spencer, R. G.: Utilization of ancient permafrost carbon in headwaters of Arctic fluvial networks, Nat. Commun., 6, 7856, https://doi.org/10.1038/ncomms8856, 2015.

Mann, P. J., Spencer, R. G. M., Hernes, P. J., Six, J., Aiken, G. R., Tank, S. E., McClelland, J. W., Butler, K. D., Dyda, R. Y., and Holmes, R. M.: Pan-Arctic Trends in Terrestrial Dissolved Organic Matter from Optical Measurements, Front. Earth Sci., 4, 25, https://doi.org/10.3389/feart.2016.00025, 2016.

Mannino, A., Russ, M. E., and Hooker, S. B.: Algorithm development and validation for satellite-derived distributions of DOC and CDOM in the U.S. Middle Atlantic Bight, J. Geophys. Res., 113, C07051, https://doi.org/10.1029/2007jc004493, 2008.

Marín-Spiotta, E., Gruley, K. E., Crawford, J., Atkinson, E. E., Miesel, J. R., Greene, S., Cardona-Correa, C., and Spencer, R. G. M.: Paradigm shifts in soil organic matter research affect interpretations of aquatic carbon cycling: transcending disciplinary and ecosystem boundaries, Biogeochemistry, 117, 279–297, https://doi.org/10.1007/s10533-013-9949-7, 2014.

Massicotte, P., Asmala, E., Stedmon, C., and Markager, S.: Global distribution of dissolved organic matter along the aquatic continuum: Across rivers, lakes and oceans, Sci. Total Environ., 609, 180–191, https://doi.org/10.1016/j.scitotenv.2017.07.076, 2017.

Mitchell, B., Kahru, M., Wieland, J., and Stramska, M.: Determination of spectral absorption coefficients of particles, dissolved material and phytoplankton for discrete water samples, NASA Technical Memorandum, 2002.

Myers-Smith, I. H., Hik, D. S., Kennedy, C., Cooley, D., Johnstone, J. F., Kenney, A. J., and Krebs, C. J.: Expansion of Canopy-Forming Willows Over the Twentieth Century on Herschel Island, Yukon Territory, Canada, Ambio, 40, 610–623, https://doi.org/10.1007/s13280-011-0168-y, 2011.

Neff, J. C., Finlay, J. C., Zimov, S. A., Davydov, S. P., Carrasco, J. J., Schuur, E. A. G., and Davydova, A. I.: Seasonal changes in the age and structure of dissolved organic carbon in Siberian rivers and streams, Geophys. Res. Lett., 33, L23401, https://doi.org/10.1029/2006gl028222, 2006.

O'Donnell, J. A., Aiken, G. R., Walvoord, M. A., Raymond, P. A., Butler, K. D., Dornblaser, M. M., and Heckman, K.: Using dissolved organic matter age and composition to detect permafrost thaw in boreal watersheds of interior Alaska, J. Geophys. Res.-Biogeo., 119, 2155–2170, https://doi.org/10.1002/2014jg002695, 2014.

Osterkamp, T. E.: Characteristics of the recent warming of permafrost in Alaska, J. Geophys. Res., 112, F02S02, https://doi.org/10.1029/2006jf000578, 2007.

Poulin, B. A., Ryan, J. N., and Aiken, G. R.: Effects of iron on optical properties of dissolved organic matter, Environ. Sci. Technol., 48, 10098–10106, https://doi.org/10.1021/es502670r, 2014.

Ramage, J. L., Irrgang, A. M., Morgenstern, A., and Lantuit, H.: Increasing coastal slump activity impacts the release of sediment and organic carbon into the Arctic Ocean, Biogeosciences, 15, 1483–1495, https://doi.org/10.5194/bg-15-1483-2018, 2018.

Ramage, J. L., Fortier, D., Hugelius, G., Lantuit, H., and Morgenstern, A.: Distribution of carbon and nitrogen along hillslopes in three valleys on Herschel Island, Yukon Territory, Canada, Catena, 178, 132–140, https://doi.org/10.1016/j.catena.2019.02.029, 2019.

Rampton, V. N.: Quaternary geology of the Yukon Coastal Plain, Geological Survey of Canada, Canada, 1982.

RStudio Team: RStudio: Integrated Development for R, Boston, MA, 2016.

Schaefer, K., Lantuit, H., Romanovsky, V. E., Schuur, E. A. G., and Witt, R.: The impact of the permafrost carbon feedback on global climate, Environ. Res. Lett., 9, 085003, https://doi.org/10.1088/1748-9326/9/8/085003, 2014.

Smith, C. A. S., Kennedy, C. E., Hargrave, A. E., and McKenna, K. M.: Soil and vegetation of Herschel Island, Yukon Territory, Ottawa, ON, 101 pp., 1989.

Spence, C., Kokelj, S., and Ehsanzadeh, E.: Cold Region Hydrology in a Changing Climate, in: Proceedings of symposium H02 held during IUGG2011 in Melbourne, Australia, July 2011, IAHS publication 346, 3–8, 2011

Spence, C., Kokelj, S. V., Kokelj, S. A., McCluskie, M., and Hedstrom, N.: Evidence of a change in water chemistry in Canada's subarctic associated with enhanced winter streamflow, J. Geophys. Res.-Biogeo., 120, 113–127, https://doi.org/10.1002/2014jg002809, 2015.

Spencer, R. G. M., Aiken, G. R., Wickland, K. P., Striegl, R. G., and Hernes, P. J.: Seasonal and spatial variability in dissolved organic matter quantity and composition from the Yukon River basin, Alaska, Global Biogeochem. Cy., 22, GB4002, https://doi.org/10.1029/2008gb003231, 2008.

Spencer, R. G. M., Aiken, G. R., Butler, K. D., Dornblaser, M. M., Striegl, R. G., and Hernes, P. J.: Utilizing chromophoric dissolved organic matter measurements to derive export and reactivity of dissolved organic carbon exported to the Arctic Ocean: A case study of the Yukon River, Alaska, Geophys. Res. Lett., 36, L06401, https://doi.org/10.1029/2008gl036831, 2009.

Stedmon, C. A., Amon, R. M. W., Rinehart, A. J., and Walker, S. A.: The supply and characteristics of colored dissolved organic matter (CDOM) in the Arctic Ocean: Pan Arctic trends and differences, Mar. Chem., 124, 108–118, https://doi.org/10.1016/j.marchem.2010.12.007, 2011.

Striegl, R. G., Aiken, G. R., Dornblaser, M. M., Raymond, P. A., and Wickland, K. P.: A decrease in discharge-normalized DOC export by the Yukon River during summer through autumn, Geophys. Res. Lett., 32, L21413, https://doi.org/10.1029/2005gl024413, 2005.

Sulzberger, B. and Durisch-Kaiser, E.: Chemical characterization of dissolved organic matter (DOM): A prerequisite for understanding UV-induced changes of DOM absorption properties and bioavailability, Aquat. Sci., 71, 104–126, https://doi.org/10.1007/s00027-008-8082-5, 2009.

Tank, S. E., Fellman, J. B., Hood, E., and Kritzberg, E. S.: Beyond respiration: Controls on lateral carbon fluxes across the terrestrial-aquatic interface, Limnol. Oceanogr. Lett., 3, 76–88, https://doi.org/10.1002/lol2.10065, 2018.

Vantrepotte, V., Danhiez, F. P., Loisel, H., Ouillon, S., Meriaux, X., Cauvin, A., and Dessailly, D.: CDOM-DOC relationship in contrasted coastal waters: implication for DOC retrieval from ocean color remote sensing observation, Opt. Express, 23, 33–54, https://doi.org/10.1364/OE.23.000033, 2015.

Vonk, J. E., Tank, S. E., Bowden, W. B., Laurion, I., Vincent, W. F., Alekseychik, P., Amyot, M., Billet, M. F., Canário, J., Cory, R. M., Deshpande, B. N., Helbig, M., Jammet, M., Karlsson, J., Larouche, J., MacMillan, G., Rautio, M., Walter Anthony, K. M., and Wickland, K. P.: Reviews and syntheses: Effects of permafrost thaw on Arctic aquatic ecosystems, Biogeosciences, 12, 7129–7167, https://doi.org/10.5194/bg-12-7129-2015, 2015a.

Vonk, J. E., Tank, S. E., Mann, P. J., Spencer, R. G. M., Treat, C. C., Striegl, R. G., Abbott, B. W., and Wickland, K. P.: Biodegradability of dissolved organic carbon in permafrost soils and aquatic systems: a meta-analysis, Biogeosciences, 12, 6915–6930, https://doi.org/10.5194/bg-12-6915-2015, 2015b.

Walker, S. A., Amon, R. M. W., and Stedmon, C. A.: Variations in high-latitude riverine fluorescent dissolved organic matter: A comparison of large Arctic rivers, J. Geophys. Res.-Biogeo., 118, 1689–1702, https://doi.org/10.1002/2013JG002320, 2013.

Wang, Y., Spencer, R. G. M., Podgorski, D. C., Kellerman, A. M., Rashid, H., Zito, P., Xiao, W., Wei, D., Yang, Y., and Xu, Y.: Spatiotemporal transformation of dissolved organic matter along an alpine stream flow path on the Qinghai–Tibet Plateau: importance of source and permafrost degradation, Biogeosciences, 15, 6637–6648, https://doi.org/10.5194/bg-15-6637-2018, 2018.

Ward, C. P. and Cory, R. M.: Chemical composition of dissolved organic matter draining permafrost soils, Geochim. Cosmochim. Ac., 167, 63–79, https://doi.org/10.1016/j.gca.2015.07.001, 2015.

Ward, C. P., Nalven, S. G., Crump, B. C., Kling, G. W., and Cory, R. M.: Photochemical alteration of organic carbon draining permafrost soils shifts microbial metabolic pathways and stimulates respiration, Nat. Commun., 8, 772, https://doi.org/10.1038/s41467-017-00759-2, 2017.

Weishaar, J. L., Aiken, G. R., Bergamaschi, B. A., Fram, M. S., Fujii, R., and Mopper, K.: Evaluation of Specific Ultraviolet Absorbance as an Indicator of the Chemical Composition and Reactivity of Dissolved Organic Carbon, Environ. Sci. Technol., 37, 4702–4708, https://doi.org/10.1021/es030360x, 2003.

Whitehead, R. F., de Mora, S., Demers, S., Gosselin, M., Monfort, P., and Mostajir, B.: Interactions of ultraviolet-B radiation, mixing, and biological activity on photobleaching of natural chromophoric dissolved organic matter: A mesocosm study, Limnol. Oceanogr., 45, 278–291, https://doi.org/10.4319/lo.2000.45.2.0278, 2000.

Woo, M.-K., Kane, D. L., Carey, S. K., and Yang, D.: Progress in permafrost hydrology in the new millennium, Permafrost Periglac., 19, 237–254, https://doi.org/10.1002/ppp.613, 2008.