# Canadian permafrost stores large pools of ammonium and optically distinct dissolved organic matter

Julien Fouché, C. T Christiansen, M. J Lafrenière, P. Grogan, S. F Lamoureux

► To cite this version:

Julien Fouché, C. T Christiansen, M. J Lafrenière, P. Grogan, S. F Lamoureux. Canadian permafrost stores large pools of ammonium and optically distinct dissolved organic matter. Nature Communications, 2020, 11 (1), 10.1038/s41467-020-18331-w. hal-02944592

HAL Id: hal-02944592

https://hal.inrae.fr/hal-02944592

Submitted on 21 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

https://doi.org/10.1038/s41467-020-18331-w

OPEN

<figure><img src="image_2.png" /></figure>

Check for updates

# Canadian permafrost stores large pools of ammonium and optically distinct dissolved organic matter

J. Fouché<fnref n="1" />,²✉, C. T. Christiansen³,⁴, M. J. Lafrenière², P. Grogan⁵,⁴ & S. F. Lamoureux⁵

Permafrost degradation may lead to mobilization of carbon and nutrients and enhance microbial processing rates of previously frozen organic matter. Although the pool size and chemical composition of dissolved organic matter (DOM) are fundamental determinants of the carbon cycle in Arctic watersheds, its source within the seasonally thawing active layer and the underlying permafrost remains largely uncharacterized. Here, we used 25 soil cores that extended down into the permafrost from nine sites across Arctic Canada to quantify dissolved organic carbon (DOC) and nitrogen stocks, and to characterize DOM optical properties. Organic permafrost stores 5-7 times more DOC and ammonium than the active layer and mineral permafrost. Furthermore, the permafrost layers contain substantial low molecular weight DOM with low aromaticity suggesting high biodegradability. We conclude that soil organic matter stoichiometry and cryogenic processes determine permafrost DOM chemistry, and that thawing will mobilize large amounts of labile DOC and ammonium into Arctic watersheds.

<footnote n="1">¹LISAH, Univ Montpellier, INRAE, IRD, Institut Agro, 34060 Montpellier, France. ²Department of Geography and Planning, Queen's University, Kingston, ON K7L 3N6, Canada. ³Center for Permafrost (CENPERM), Department of Geoscience and Natural Resource Management, University of Copenhagen, DK-1350 Copenhagen, Denmark. ⁴Department of Biology, Queen's University, Kingston, ON K7L 3N6, Canada. ✉email: julien.fouche@supagro.fr</footnote>
Global warming will lead to massive declines in the extent of northern permafrost soils during this century¹. Permafrost degradation occurs through either a gradual thickening of the seasonally thawed active layer or the development of thermokarst landforms as land surface subsidence follows thawing of ice-rich permafrost². Permafrost-affected soils store an estimated 1000 ± 150 Pg organic carbon at 0–3 m depth³ and climate warming-induced thaw therefore has the potential to increase C emissions to the atmosphere through enhanced microbial decomposition of organic matter². In addition to being a reactive and easily mobilized pool of soil organic matter (DOM), dissolved organic matter (DOM) is an important substrate for microorganisms and a crucial component of overall net ecosystem carbon balance in Arctic watersheds⁴. Yet, research on the permafrost DOM pool size and composition remains scarce⁵⁻⁸ despite a clear need to evaluate permafrost thaw as a potential positive feedback on climate change. This study is the first to characterize the permafrost DOM pool systematically and to estimate the dissolved nitrogen pool across Northern Canada, a severely underrepresented region in climate-change research⁹.

Permafrost degradation alters ground surface hydrology and strongly modifies the amount and composition of DOM and nutrients which are mobilized from terrestrial to aquatic ecosystems¹⁰. Warming and permafrost thaw increase dissolved inorganic nitrogen (DIN) pool sizes through organic matter mineralization and additional permafrost DIN release¹¹, promoting deep-rooted plant growth¹², changing microbial and plant community composition¹³ while increasing microbial metabolic efficiency¹⁴, and significantly affecting carbon and nitrogen cycling in northern ecosystems¹⁵.

In addition to the amount of DOM exported, watershed carbon emissions partly depend on DOM biodegradability, which is determined by its chemical composition¹⁶, its sorption to soil minerals¹⁷, stage of degradation and photo-oxidation¹⁸, abiotic conditions¹⁹, nutrient availability¹⁸, and microbial community functionality¹⁶. Current evidence suggests that a substantial portion of permafrost DOM is relatively quickly mineralized by microbes and/or photo-oxidized after permafrost thaw during its transport from soils to rivers and towards the Arctic Ocean⁵. However, to date, the characterization of the composition and lability of permafrost DOM remains limited to only a few Arctic locations⁵,⁶,¹⁰,¹⁹.

In order to characterize and quantify DOM and DIN pools stored in Canadian permafrost soils, we utilized 25 active layer and permafrost (upper 3 m) cores from a wide range of permafrost-affected settings that varied in formation histories, ages, surficial geology, and climate conditions. For all cores, we grouped water extracts from the active layer or permafrost samples into either organic (soil C ≥ 12%) or mineral (soil C < 12%) layers. Soil samples were gently leached during a short period by using milliQ water in order to simulate the leaching that would occur immediately upon permafrost thaw. We analyzed each sample for dissolved organic carbon (DOC), total dissolved nitrogen (TDN), and DIN, and DOM optical properties using absorbance and fluorescence spectroscopy. In addition to using optical indices to characterize DOM aromaticity and molecular weight²⁰, we performed PARAFAC modeling to identify the main chemical components explaining the fluorescent DOM pool²¹. By quantifying and characterizing DOM in permafrost soils from a wide range of formation histories, ages, surficial geology, and climate conditions, we are able to highlight the distinct signature of DOM in permafrost layers compare to the active layer and to investigate the controls of DOM character. In addition, we estimate that large of pools of DOM and dissolved nitrogen have the potential to be released in surface waters with active layer thickening or the thaw induced exposure of permafrost.

## Results and discussion

**Properties of permafrost soils.** All sites contained syngenetic permafrost in which the active layer and the uppermost permafrost have experienced numerous freeze-thaw cycles since formation during the Holocene²². Permafrost organic matter radiocarbon ages ranged from 7850 ± 30 to 830 ± 20 y B.P. (Supplementary Data), with the western sites containing the oldest SOM and the northern Hudson Bay peatlands containing the youngest SOM.

For organic layers, permafrost soil C:N atomic ratios (14.50 [12.61–19.67], median [25th–75th]) were lower and H:C atomic ratios (0.13 [0.13–0.14]) were greater relative to the active layer, 23.89 [19.33–29.50] and 0.14 [0.12–0.15], respectively (Supplementary Fig. 1). Similarly for mineral layers, permafrost C:N ratios were lower (12.0 [3.23–18.09]) and H:C ratios higher (0.24 [0.15–0.86]) compared to the active layer, 16.21 [13.67–18.87] and 0.17 [0.15–0.21], respectively. For both thermal layers, in organic layers C:N ratios were higher and H:C ratios were lower than in mineral layers.

These stoichiometric properties are typical of boreal and tundra soils (Supplementary Fig. 1)²³. The higher C:N and very low H:C ratios of the organic layers relative to mineral layers suggest higher contents of condensed aromatic structures originating from peat²⁴. Permafrost layers displayed lower C:N properties suggesting different SOM composition (e.g., lignin, tannins, lipids, sugars or amino acids) and an enrichment in microbial biomass relative to the active layer²⁴. The absence of a downward trend of C:N and H:C within the permafrost (Supplementary Fig. 1), except at Daring Lake, indicates that soil development and microbial processing were effectively halted soon after permafrost aggradation²³.

**Active layer and permafrost yields of DOC and nitrogen.** DOC content correlated with soil C content in both active layer ($r^2$ [log-log] = 0.748, *P* < 0.001) and permafrost ($r^2$ [log-log] = 0.834, *P* < 0.001, Supplementary Fig. 2). DOC and TDN pools were larger in organic layers relative to mineral soil layers (Fig. 1a, b). Organic permafrost yielded greater amounts of DOC and TDN than the active layer and mineral permafrost (Fig. 1a, b). We estimate that organic permafrost layers contained 0.396 [0.178–0.644] kg DOC m⁻³, and 37.6 [15.7–59.8] 10⁻³ kg TDN m⁻³, potentially releasing 5.3 times more DOC and 7.2 times more TDN upon thaw than either of the organic active layer, the mineral active layer, and the mineral permafrost layer (DOC: Kruskal–Wallis chi-squared = 103.22, *P* < 0.001, TDN: Kruskal–Wallis chi-squared = 85.3, *P* < 0.001, Fig. 1a, b), respectively. DOC pools in both active layer and permafrost were similar to values reported for mineral layers in Pleistocene Yedoma permafrost⁵ and in Orthels and Turbels from interior Alaska⁶, as well as for organic layers in Histels from Alaska and Siberia⁶,²⁵.

Within the active layer, organic and mineral soil samples displayed similar DOC (W = 48, *P* = 0.941) and TDN yields (W = 35, *P* = 0.963). Both DOC and TDN yields increased in the permafrost relative to the active layer (Fig. 1c, d), with mineral layers displaying greater DOC yields than organic layers in the permafrost (W = 1077, *P* = 0.019), extracting 1.3% [0.6–3.1] of soil C and 0.9% [0.3–1.4] of soil N, respectively. Our yields were in the range of values reported for permafrost soils in North America²⁵,²⁶.

In the active layer, dissolved organic nitrogen (DON = TDN – DIN) dominated the TDN pool representing ~78% of TDN in both mineral and organic layers. The DIN contribution increased with depth and ammonium-N dominated the DIN pool, with its concentration reaching 0.410 ± 0.326 mg N-NH₄⁺ g⁻¹ soil at
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Dissolved organic carbon and total dissolved nitrogen pools and yields in the active layer and permafrost. Samples are grouped by layer type with the following classification: organic layers (soil C ≥ 12%) and mineral layers (soil C < 12%). Letters indicate significant differences between thermal layers and soil types indicated by Mann-Whitney U test. The box plots summarize the distribution of **a** the amounts of dissolved organic carbon (DOC, kg C m⁻³) and **b** total dissolved nitrogen (TDN, kg N m⁻³) and the **c** water extractable yields in DOC (g DOC 100 g soil C⁻¹) and **d** TDN (g TDN 100 g soil N⁻¹) for each layer type and in the active layer and underlying permafrost (calculated from multiple 25 cm increments). In each box plot, the diamond represents the mean, the horizontal line represents the median, the end of the box the 25th and 75th percentiles, and the lines extending from the box are 1.5 interquartile ranges from the median. Data points outside of the 1.5 interquartile ranges are represented as dots.</figcaption></figure>

1–2 m depth in organic permafrost layers (Supplementary Fig. 3). In contrast, nitrate concentrations remained very low in both mineral and organic soils with 67 of 220 samples having NO₃⁻ concentrations below our instrument detection limit (0.006 mg N L⁻¹, Fig. 2, Supplementary Fig. 3). The TDN, NH₄⁺, and NO₃⁻ concentrations measured in this study were within the range of values reported for permafrost soils6,27–29.

While the active layer and mineral permafrost layers contained relatively little DIN, the pool of NH₄⁺ in the organic permafrost layers (32.8 [11.3–66.2] g N-NH₄⁺ m⁻³, Fig. 2) was very large compared to the general low nitrogen availability in Arctic environments¹³. Similar enrichment in NH₄⁺ were reported in various permafrost settings in Sweden, Siberia, and Greenland²⁸⁻³⁰ as being a particular feature of permafrost soils in the Arctic. Saturated and oxygen-limited conditions in peatlands and syngenetic permafrost aggradation likely prevented microbial processing and/or leaching of NH₄⁺ out of the soil profile. Seasonal NH₄⁺ leaching down through the active layer at the end of the growing season, leading to NH₄⁺ accumulation at the permafrost table during the winter freeze-up, likely contributed to the accumulation of a large pool of NH₄⁺²⁸ in the solute-enriched permafrost transition zone³¹ on decadal to centennial time scales. In addition, as microorganisms in permafrost soils are active under freezing conditions³², and some pore water remains within the permafrost, at least some of the large NH₄⁺ pool likely originates from in-situ mineralization below 0 °C over thousands of years6,28. A thaw-induced release of permafrost NH₄⁺ of this magnitude will likely affect Arctic terrestrial and aquatic ecosystems changing

trophic structures and providing new habitat niches¹³ while altering landscape-wide carbon cycling processes¹⁵.

**Optical properties of DOM in permafrost soils.** The PARAFAC modeling for DOM fluorescence was conducted in all permafrost soils, spanning a wide range of surficial geology and climate conditions. We related six identified fluorescent components to fractions previously reported in surface waters in northern latitudes (Fig. 3, Table 1). Components C2, C3, C4, and C5 were previously related to humic-like or fulvic-like components associated with high molecular weight (HMW), aromatic organic compounds originating from terrestrial sources such as plant inputs and SOM³³⁻³⁵. Components C1 and C6 are ubiquitous across a wide range of terrestrial and marine environments and correspond, respectively, to tyrosine-like (C1) and tryptophan-like (C6) components, representing proteinaceous compounds from microbial activity such as amino acids, peptide materials and free or bound proteins³⁴.

The optical properties of DOM differed between thermal layers and between organic and mineral soil leachates (Fig. 4). Permafrost DOM exhibited a distinctive signature with a greater proportion of low molecular weight (LMW) proteinaceous organic compounds than found in the active layer (Fig. 4). In the active layer, the HMW aromatic fluorophores derived from terrestrial sources were the most abundant in both mineral and organic layers (94.7% [88.3–96.7] of the total fluorescence). These fluorophores contributed less to the fluorescent DOM,
<figure><img src="image_4.png" /></figure>

**Fig. 2 Pools of ammonium and nitrate in the active layer and permafrost.** Samples are grouped by layer type with the following classification: organic layers (soil C ≥ 12%) and mineral layers (soil C < 12%). Letters indicate significant differences between thermal layers and soil types indicated by Mann-Whitney U test. The box plots summarize the distribution of **a** the amounts of ammonium (NH₄⁺) (kg N m⁻³) and **b** nitrate (NO₃⁻) (g N m⁻³) for each layer nature and in the active layer and permafrost. In each box plot, the diamond represents the mean, the horizontal line represents the median, the end of the box the 25th and 75th percentiles, and the lines extending from the box are 1.5 interquartile ranges from the median. Data points outside of the 1.5 interquartile ranges are represented as dots.

<figure><img src="image_5.png" /></figure>

**Fig. 3 Six fluorescent components identified using PARAFAC analysis.** Excitation and emission peak positions of the three-dimensional excitation-emission matrices of the independent components are indicated alongside descriptions in Table 1.

representing 77.0% [57.9–88.6] and 40.7% [23.9–70.8] of the fluorescence in permafrost mineral and organic layers, respectively (Supplementary Fig. 4). The tyrosine-like component (C1) was the most abundant fluorophore in the organic permafrost layers accounting for 51.2% [22.4–62.9] of the total fluorescence, and the tryptophan-like component (C6) representing 7.2% [5.0–10.1]. In the mineral permafrost, the tyrosine-like and tryptophan-like components represented 12.8% [2.7–34.4] and 7.6% [5.5–11.0] of the total fluorescence, respectively (Supplementary Fig. 4). These components contributed less in the active

layer than in the permafrost representing 5.3% [4.1–22.8] in organic and 5.3% [3.3–10.7] in mineral layers of the fluorescent DOM (Supplementary Fig. 4).

The multivariate statistical analyses based on absorbance and fluorescence indices demonstrated that permafrost DOM exhibited a distinct optical signature shared by both mineral and organic layers (Fig. 4a, b). The permafrost leachates were enriched in LMW proteinaceous compounds with low aromaticity while the active layer DOM was comprised of HMW aromatic condensed organic compounds. The fluorescence index
<table><caption>Table 1 Description of the six components identified by PARAFAC of 226 excitation-emission matrices.</caption><thead><tr><th>Component</th><th>Excitation max (nm)</th><th>Emission max (nm)</th><th>Description and likely structure</th><th>Cory and McKnight (2005) Polar lakes and other freshwaters</th><th>Fellman et al. (2010) Review</th><th>Mann et al. (2016) in main Arctic rivers</th><th>Walker et al. (2013) in main Arctic rivers</th><th>Openfluor comparison<sup>a</sup> (component with TCC &gt; 0.95)</th><th>Source of DOM</th></tr></thead><tbody><tr><td>1</td><td>275</td><td>306</td><td>Tyrosine-like— Amino acid, soluble bound microbial DOM, may indicate more degraded peptide material</td><td>Comp 13</td><td>C1</td><td>-</td><td>-</td><td>15</td><td>Terrestrial, Autochthonous, Marine</td></tr><tr><td>2</td><td>310</td><td>415</td><td>UVC humic-like—High molecular weight, syringaldehyde-like (byproducts of lignin breakdown), associated with high DOM contents</td><td>Comp 10</td><td>C7</td><td>AG3</td><td>C1</td><td>23</td><td>Terrestrial</td></tr><tr><td>3</td><td>&lt;255</td><td>450–500</td><td>UVA humic-like or fulvic like— High molecular weight, syringaldehyde-like (byproducts of lignin breakdown), associated with high DOM contents</td><td>Q2 (Comp 2)</td><td>C5</td><td>AG1</td><td>C4</td><td>12</td><td>Terrestrial</td></tr><tr><td>4</td><td>360 (&lt;255)</td><td>460</td><td>UVC humic-like—High molecular weight and aromatic humic acid</td><td>SQ2 (Comp 7)</td><td>C10</td><td>AG4</td><td>C2</td><td>22</td><td>Terrestrial</td></tr><tr><td>5</td><td>295–300</td><td>520–540</td><td>Semiquinone-like—High molecular weight and highly aromatic humic acid</td><td>SQ1 (Comp 5)</td><td>-</td><td>-</td><td>C3</td><td>1b</td><td>Terrestrial</td></tr><tr><td>6</td><td>280</td><td>345</td><td>Tryptophan-like—Amino acid, soluble free or bound microbial DOM, may indicate intact proteins or less degraded peptide material</td><td>Comp 8</td><td>C2</td><td>AG7</td><td>C5</td><td>24</td><td>Terrestrial, Autochthonous, Marine</td></tr></tbody></table>

<footnote>a number of studies in the open fluorescence database which previously identified components displaying similar optical properties (tucker congruence coefficient, TCC > 0.95).</footnote>

<footnote>b Excitation has a TCC = 0.95.</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Principal component analysis of DOM optical properties. The analyses takes into account the relative contribution of PARAFAC components and optical indices: the absorption coefficient at 350 (E<sub>350</sub>), the specific UV absorbance (SUV<sub>A254</sub>, L mg C<sup>-1</sup> m<sup>-1</sup>), the absorption spectral slope over the spectral band 275-295 nm (S<sub>295</sub>), the slope ratio between S<sub>295</sub> and the spectral slope over 350-400 nm (S<sub>R</sub>), the fluorescence index (FI), the freshness index (BIX), the humification index (HIX). **a** Explanatory variable loadings are shown as black dots and scores across the first and second principal components for all samples labeled by layer type—organic layers are shown as circles and mineral as diamonds—and colored by location in the active layer or permafrost. **b** Box plot of the scores of the first principal component for organic (soil C ≥ 12%) and mineral layers (soil C < 12%) in the active layer and in the permafrost. **c** Box plot of the scores of the second principal component for organic (soil C ≥ 12%) and mineral layers (soil C < 12%) in the active layer and in the permafrost. In each box plot, the diamond represents the mean, the horizontal line represents the median, the end of the box the 25th and 75th percentiles, and the lines extending from the box are 1.5 interquartile ranges from the median. Data points outside of the 1.5 interquartile ranges are represented as dots. Letters indicate significant differences between thermal layers and soil types indicated by Mann-Whitney U test.</figcaption></figure>

FI (a common proxy for the relative contribution between microbial and terrestrial DOM) and the ratio of absorption spectral slopes S<sub>R</sub>, which decreases with molecular weight and aromaticity, increased with depth (Fig. 4a, b, Supplementary Fig. 5) indicating more protein-derived substrate with a lower average molecular weight in the permafrost<sup>34</sup>. HIX and SUV<sub>A254</sub> values, which are known to increase with bulk aromaticity and condensing, decreased with depth and in the permafrost (Supplementary Fig. 5) suggesting that permafrost stores less condensed and aromatic DOM than in the active layer<sup>20,34</sup>. Similar values of HIX and SUV<sub>A254</sub> in active layer mineral and organic layers illustrate the similar long-term microbial processing of organic matter in the seasonally thawed zone Supplementary Fig. 5). In both thermal layers, DOM in organic and mineral layers only differed by the absorption at 350 nm, with organic layers displaying higher absorption values due to greater concentrations of aromatic DOM<sup>33</sup> (Fig. 4a, c).

Our values of optical indices are consistent with the few available characterizations of DOM in the active layer and permafrost<sup>5-8,25,36</sup>. Due to different organic matter sources, DOM in mineral permafrost layers has been previously characterized by lower SUV<sub>A254</sub> values ([0.6-1.2] L mg C<sup>-1</sup> m<sup>-1</sup>)<sup>5,8</sup> than organic active layer materials ([1-4.5] L mg C<sup>-1</sup> m<sup>-1</sup>)<sup>6,8,25,37</sup> due to a relative enrichment in microbial exometabolites and root exudates<sup>37</sup> in mineral soils compared to peaty organic layers. Last, molecular techniques also support our findings by showing enrichment of LMW compounds such as carbohydrates<sup>8</sup>, acetate and butyrate<sup>5</sup> in the permafrost in Alaska and Siberia.

**Controls of the pan-Canadian permafrost DOM pool.** Even though permafrost DOM displays contrasting optical properties from active layer DOM, the signature of permafrost DOM is weakly detected within Arctic catchments. This is mainly due to the high reactivity of DOM, which experiences sorption and desorption when percolating in the soil<sup>25</sup>, photo-oxidation of HMW aromatic components in streams and rivers, and microbial

processing of less aromatic compounds<sup>10,18,25</sup>. Indeed, while our results show the high contribution of tyrosine-like component to the fluorescent DOM in both mineral and organic permafrost in the Canadian Arctic region, studies that have conducted PARAFAC modeling on DOM in Arctic surface waters have only detected a low contribution of this component in headwater catchments<sup>38,39</sup>, with no contribution at all in Arctic rivers<sup>33,35</sup>. The major tyrosine contribution to the permafrost fluorescent DOM and its disappearance within the fluvial continuum provides additional evidence of its high reactivity<sup>5,10,19</sup>.

Following permafrost thaw, three main processes lead to the decrease in the proteinaceous contribution during lateral flow downslope along the permafrost table, prior to reaching water bodies. First, soil microbes could preferentially mineralize tyrosine-like DOM components, such as amino acids and oligopeptides<sup>40</sup>, as they utilize aliphatic compounds<sup>41</sup> almost immediately following permafrost thaw<sup>5</sup>. Recent studies in Arctic permafrost environments highlight that low SUV<sub>A254</sub><sup>19</sup>, high S<sub>R</sub><sup>42</sup> values and a high contribution of tyrosine-like<sup>43</sup> correlated with greater lability of bulk DOM<sup>6,36,43,44</sup>. Therefore, we suggest that Canadian permafrost likely contains a high proportion of biodegradable DOM. Secondly, in northern ecosystems plants take up amino acids and oligopeptides directly as a source of N, leading to a decrease in proteinaceous compounds in exported waters<sup>45</sup>. Finally, the retention of less aromatic compounds and protein-like fluorophores in mineral subsoils and the preferential release of HMW aromatic fluorophores lead to the depletion of tyrosine-like components in riverine fluorescent DOM and the relative enrichment in aromatic compounds<sup>25</sup>.

Our results demonstrate robust relationships between SOM stoichiometry (C:N and H:C ratios) and chromophoric and fluorescent DOM properties and DOC yield highlighting the control of permafrost organic matter on the DOM pool (Fig. 5). E<sub>350</sub> decreased with increasing H:C ratios (Fig. 5a), suggesting that chromophoric DOM abundance was greater for C-rich SOM<sup>24</sup>. The C:N ratio was negatively correlated with S<sub>R</sub> and FI (Fig. 5b, c) demonstrating that N-rich SOM (i.e., enriched in
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Relationships between indices of DOM character and soil properties. Relationships between E<sub>350</sub> (a), S<sub>R</sub> (b), FI (c, e), SUVA<sub>254</sub> (L mg C<sup>-1</sup> m<sup>-1</sup>) (e)) and soil properties (H:C (a), C:N (b, c, f), and DOC yield (g DOC 100 g soil C<sup>-1</sup>) (d, e, f)) for both active layer and permafrost samples. Each point represents one sample. The linear regression lines and equations are for active layer and permafrost samples.</figcaption></figure>

microbial biomass) yields LMW DOM with low aromaticity. In addition, SOM properties (i.e., C:N ratio), the degree of aromaticity (i.e., SUVA<sub>254</sub>) and the microbial origin of fluorescent DOM (i.e., FI) was correlated with the water extractable DOC yield (Fig. 5d–f). Therefore, permafrost samples with lower C:N ratio yield proportionally more DOM that contains a higher proportion of less aromatic and more proteinaceous compounds, which are likely labile<sup>44</sup> relative to the active layer. As soil C:N and H:C ratios are widely measured in permafrost soils<sup>23</sup>, these findings allow for upscaling of site or plot level analyses to a circumpolar characterization of the DOM pool.

Our results further indicate that the uppermost permafrost degradation has the potential to release a larger amount of DOC and TDN (mostly as DON and NH<sub>4</sub><sup>+</sup>) than the active layer currently releases to high-latitude terrestrial and aquatic ecosystems. These findings are in agreement with observations of thawing permafrost leading to shifts in the supply of carbon and nutrients to surface waters<sup>46,47</sup>. Using the NCSCDv2 dataset<sup>3</sup> for the pan-Canadian permafrost area, we estimated the DOC pool available right upon thaw increasing from 0.2 [0.2–0.5] Pg C in 0–1 m, to 1.1 [0.7–2.1] Pg C in 1–2 m and decreasing to 0.9 [0.7–1.5] Pg C in 2–3 m depth interval. Although DOC is a relatively small component of the total permafrost carbon pool, its depletion through hydrological mobilization will likely contribute relatively fast to soil carbon loss in Arctic landscapes<sup>48</sup>.

In addition to an increase in DOC availability in the large areas of peatlands located in southern boundaries of the permafrost region (i.e., the Mackenzie River region, the Hudson Bay lowlands), which have been subject to rapid thaw<sup>49</sup> (Fig. 6a–c), the DOC pool increased with volumetric water content (Fig. 6d). While being more prone to thermokarst processes<sup>2</sup>,

the degradation of ice-rich permafrost could potentially export more carbon to surface waters. Our data further show that the DOC pool was positively correlated with the NH<sub>4</sub><sup>+</sup> pool and the large DOM pools are characterized by lower aromaticity (i.e., low SUVA<sub>254</sub>) (Fig. 6e, f). This finding supports reported observations of the accumulation of NH<sub>4</sub><sup>+</sup> and biodegradable DOC in Arctic permafrost soils<sup>47</sup>. While proteins compose most of the DON that are thereafter degraded to amino acids and mineralized to NH<sub>4</sub><sup>+</sup>, we demonstrate here the buildup of both NH<sub>4</sub><sup>+</sup> and tyrosine-like fluorophores in the permafrost, which then halted their decomposition and export. As reported in glacier ice, the dominant contribution of the tyrosine-like component to fluorescent DOM is evidence of microbial metabolism in permafrost<sup>50</sup>. The distribution of ground ice has been proven to be the main factor controlling the fate of permafrost carbon (i.e., thermokarst activity and methane production)<sup>2</sup>. Our data demonstrate that ice-rich permafrost soils have the potential to release great amount of NH<sub>4</sub><sup>+</sup> and DOM characterized by a specific low aromaticity protein-rich optical signature that suggests rapid mineralization potential.

The uppermost layers of Canadian permafrost, found immediately below the seasonally thawing active layer, currently store substantial DOC and NH<sub>4</sub><sup>+</sup> pools, particularly in organic permafrost soils. This suggests that permafrost thaw—in addition to the release of reactive DOM—is likely to enhance nutrient inputs and availability into terrestrial, as well as aquatic systems. This nutrient priming mechanism will almost certainly interact with and affect microbial processing of organic matter and primary productivity, altering the net carbon balance of Arctic ecosystems<sup>15</sup>. Here, we show that permafrost DOM exhibits a highly consistent and distinctive chemical composition across a
<figure><img src="image_8.png" /><figcaption>**Fig. 6 Maps and controls of the DOC pool in the pan-Canadian permafrost area.** a–c Maps of pool of DOC (kg C m⁻³) that would be available upon permafrost thaw for three soil layers: 0–1 m, 1–2 m, 2–3 m. The DOC pool in permafrost soils was estimated by multiplying soil carbon pool of the three different gelisol types from the NCSCDv2 by the median of water extractable DOC yield of mineral and organic soils. Blue dots represent the locations of the coring sites. The dark gray area represents the non-permafrost soils within the Canadian permafrost distribution zone (isolated, sporadic, discontinuous and continuous). **d** Volumetric water content (%), **e** pool of NH₄⁺ (g N m⁻³), and **f** SUVA₂₅₄ (L mg C⁻¹ m⁻¹) as a function of the pool of DOC (kg C m⁻³) for both active layer and permafrost samples. Each point represents one sample. The regression lines (dashed lines) and equations are for active layer and permafrost samples combined.</figcaption></figure>

very large region of northern Canada and spanning numerous distinct environments and climates. In conclusion, our results support the premise that degradation of organic rich permafrost will increase carbon release and potentially affect aquatic systems through carbon and nutrient additions.

## Methods

This study initiated from work conducted as part of the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Frontiers project Arctic Development and Adaptation to Permafrost in Transition (ADAPT), which brought together 15 laboratories from across Canada to develop an integrated Earth systems science framework on diverse aspects related to thawing permafrost conditions in the Canadian Arctic²². Data produced with the standard protocols for cryostratigraphy and soil analyses from the ADAPT project are published through the Nordicana D database⁵¹ and briefly described below.

**Permafrost cores.** We investigated 25 soil cores that were collected by members of the ADAPT project teams²² in nine distinct regions across the Canadian Arctic (Table 1). The cored areas represent a broad range of Canadian permafrost conditions, soil types, and landscapes covered by diverse ecosystems such as boreal forest, palsa peatlands, and various tundra vegetation, including dry heath, mesic shrub tundra, and wet sedge. Each soil core was obtained in 2012 and 2013 and consisted of the seasonally thawing active layer and uppermost permafrost soil layers. Five sites lack active layer samples (Arviat 1A, SAS A and B, VDT A and B). After retrieval, the soil cores were kept frozen in the field and shipped frozen to Laval University, Quebec City, Quebec, Canada. Here, the soil cores were stored frozen until subsampled into 5 cm segments in 2014. In the following analyses, each meter of a soil core contributes 3–10 subsamples. In each sampled region, we obtained duplicate soil cores from at least two distinct sites differing in soil and/or vegetation type, soil moisture, and topography. Measurements of soil moisture and bulk density and the water extracts were performed at Queen's University on frozen soil core samples. All soil samples were classified as located in either active layer or permafrost soil layers according to stratigraphic information based on CT scans and site descriptions and classified as organic or mineral layers following the USDA Soil Taxonomy⁵². Our active layer thickness, i.e., maximum thaw depth during summer, is indicative of sampling year and it does not take decadal and local spatial variability into account.

**Soil properties.** After being oven dried at 105 °C, all subsamples were analyzed for total mass content of carbon (C), nitrogen (N), and hydrogen (H) by combustion using a LECO CHN628 Elemental Analyzer (LECO Corporation, St. Joseph, MI) at the Centre d'Etudes Nordiques (CEN) Radiocarbon Dating Laboratory at Laval University, Quebec⁵¹. At Laval University, prior to combustion, we determined bulk density and gravimetric water content and homogenized the subsample, using a mortar and pestle. For additional subsamples, conventional ¹⁴C dating was done at the CEN Radiocarbon Laboratory, and samples for ¹⁴C-AMS dating were prepared at CEN and dated at the University of California Riverside radiocarbon facility (selected sites only; Table S1). Further information about CHN and radiocarbon analyses are available http://www.cen.ulaval.ca/en/page.php?lien=labradio.

The total carbon content comprises organic and inorganic carbon, soil samples were not acidified prior to analyses. Based on site descriptions, peat properties (i.e., pH, parental material), no studied soils were developed from carbonate bedrock, and soil C and N contents were strongly correlated in mineral soils (r = 0.98, P < 0.001), and previous studies in the same sites²²,⁵³⁻⁵⁵, we assume very little inorganic carbon (~<1%) contribute to the total carbon content in both organic and mineral layers.

Due to the lack of clay-rich layers with more than 2% C, we consider soil materials as mineral layers when the C content was below 12% (by weight) and as organic layers when C content was above 12% C (by weight)⁵².

**Water extracts and analyses.** For all active layer and permafrost soil samples, 20 g of frozen soil from a 5 cm segment was transferred to a glass beaker and allowed to thaw at room temperature for 30 min before 100 mL of deionized water was added (mass to volume ratio 1:5). Each soil and water slurry was homogenized for 1 h, using a shaker table, after which slurries were allowed to settle for 30 min before vacuum-filtration through glass fiber filters (Fisher G4; 1.2 mm pore-size). Water extracts were stored frozen in HDPE vials until further analyses. We analyzed water extracts of these soils for dissolved organic carbon (DOC) and nitrogen (DTN), ammonium (NH₄⁺) and nitrate (NO₃⁻), and we characterized the composition of DOM using UV-Vis absorbance and fluorescence spectroscopy. All analyses were conducted in the Department of Geography and Planning at Queen's University. DOC and TDN measurements were performed simultaneously using high-temperature combustion on a Shimadzu TOC-V Analyzer with a TMN-1 chemiluminescence detection unit (Shimadzu North America, Columbia, MD). Milli-Q blanks were run at the beginning and throughout every run to ensure consistency. A new calibration was generated for each DOC and TDN run and
based on replicate analyses of standards, which were analyzed within every run. DOC and TDN were calculated as the mean of between three and five injections with the coefficient of variance always <2%. For DOC analyses, the detection limit was 0.14 mg L⁻¹ and analytical errors were 2.9%. For TDN analyses, the detection limit was 0.02 mg L⁻¹ and analytical errors were 4.0%. Concentrations of NH₄⁺ and NO₂⁻-NO₃⁻ were measured by colorimetry using a flow Analyzer (Astoria-Pacific 2-channel Flow Analyzer, EPA Method 353.2). For NH₄⁺, the detection limit was 0.006 ppm N and the analytical error was 2.6%. For NO₃⁻, the detection limit was 0.006 ppm N and the analytical error was 1.6%. Method Detection Limit was determined by running a low-concentration standard when samples were run, and calculated as concentration of DOC/TDN/NH₄⁺/NO₃⁻ in the blanks over all the runs plus thrice standard deviation of the concentrations of DOC/TDN/NH₄⁺/NO₃⁻ in the low-concentration standard over all the runs. The analytical error was calculated as the ratio of the standard deviation of the concentrations of DOC/TDN/NH₄⁺/NO₃⁻ in the 'method check' standard multiple times over all runs by the average concentrations of DOC/TDN/NH₄⁺/NO₃⁻ in the 'method check' standard multiplied by 100. All calculations of water extracts sample elemental contents and pools were corrected for dilution associated with the deionized water addition and sample-specific gravimetric soil moisture content.

We applied seven different indices to investigate the DOM optical properties: the absorption coefficient at 350 (E₃₅₀), the specific UV absorbance (SUVA₂₅₄, L mg C⁻¹ m⁻¹), the absorption spectral slope over the spectral band 275-295 nm (S₂₉₅), the slope ratio between S₂₉₅ and the spectral slope over 350-400 nm (S<sub>R</sub>), the fluorescence index (FI), the freshness index (BIX), the humification index (HIX). Excitation-emission matrices (EEMs) were collected and decomposed into individual fluorescent components using a parallel factor analysis (PARAFAC).

UV-visible absorbance and fluorescence measurements and Emission Excitation Matrices were collected at room temperature using a Horiba Aqualog (Horiba-Jobin Yvone Scientific Edison, NJ). Each sample was run individually in 1-cm path length quartz cuvette. Absorbance spectra were blank-corrected using Milli-Q water. UV-visible absorbance spectra and excitation values were collected between 240 and 600 nm at 3 nm intervals and emission wavelengths were scanned from 214.16 to 621.03 nm at 3.15 nm intervals. The samples were run at varied integration time and gain settings to optimize the quality of the results for each sample. EEMs files were corrected for blanks and inner filter effects, 1st and 2nd order Rayleigh scatter effects were removed using the manufacturer correction procedure. The signal was normalized using a 1 ppm Quinone Sulfate standard, made up in 0.05 M H₂SO₄, and run daily during measurements, to account for fluctuations in the instrument light source characteristics. All intensities were reported in terms of Quinone Sulfate Units (QSU).

Absorbance data were converted to Napieran absorption coefficients (E in m⁻¹) by multiplying raw absorbance values by 2.303 and dividing by the cuvette path length (m)⁵⁶. E₃₅₀ was calculated and represents a quantitative measure of the chromophoric fraction of DOM that is positively correlated with DOC and lignin phenol content⁵⁷. SUVA₂₅₄ was calculated as the decadal UV absorbance at 254 nm divided by the DOC concentration²⁰. SUVA₂₅₄ is positively correlated to bulk DOM aromaticity because aromatic compounds absorb more light in the UV-visible spectra. The DOM derived from plant and soil organic matter dissolution, which is composed of large biopolymers, shows higher aromaticity than autochthonous and microbially derived DOM that comprised small biopolymers and monomers assimilable and produced by microorganisms²⁰,⁵⁸. The absorption spectral slope S₂₉₅ was calculated over the spectral band 275-295 nm. S<sub>R</sub> was calculated as the ratio between S₂₉₅ and the absorption spectral slope over 350-400 nm⁵⁹. S₂₉₅ and S<sub>R</sub> differ between DOM sources being indicators of molecular weight and aromaticity. Higher slope coefficients are associated with lower molecular weight and decreasing DOM aromaticity¹⁸,⁵⁹,⁶⁰.

FI, which was calculated as the ratio of emission at 470 to 520 nm at an excitation wavelength of 370 nm⁶¹,⁶², illustrates the source of DOM from terrestrially derived (plant and soil organic matter, ~1.2) to microbial (bacteria and algae by-products, ~1.8)³⁴. BIX is calculated as the ratio of emission at 380 nm divided by the emission intensity maximum observed between 420 and 435 nm at an excitation wavelength of 310 nm⁶³. BIX is an indicator of the relative freshness of the bulk DOM, increasing with more recently derived DOM³⁴. HIX was calculated as the area under the emission spectra 435-480 nm divided by the peak area over two spectral bands 300-345 nm + 435-480 nm, at an excitation of 254 nm³⁴,⁶⁴. An increase in HIX is associated with a shift of emission spectra toward longer wavelengths and illustrates the increased contribution of HMW aromatic fluorescing molecules with lower H:C ratios³⁴.

PARAFAC models decompose data that are arranged in three-dimensional arrays (sample × excitation wavelength × emission wavelength) into accurate spectra and relative concentrations of known fluorescent groups showing similar chemical composition in a complex DOM mixture²¹,⁶⁵. PARAFAC individual components refer to groups as 'humic-like, fluvic-like or protein-like' with various properties rather than pure organic molecules³⁴.

PARAFAC modeling was performed using the DOMFluor toolbox within MATLAB⁶⁶ in order to decompose EEMs into independent fluorescent components. Excitation wavelengths below 255 nm and emission wavelengths below 280.81 were removed from matrices in order to avoid deteriorating signal-to-noise ratios⁶⁷. Prior to modeling, regions of the spectra influenced by 1st and 2nd order scatter peaks and with no fluorescence (i.e., emission << excitation) were cut and replaced by missing values. A six-component model was generated

incorporating a total of 226 samples. Model parameters were constrained to be non-negative. The model, which identified components with unimodal emission maxima, was validated using residual analysis characterized by instrument noise and a lack of systematic pattern, split-half analysis and random initialization²¹,⁶⁶ and explained 99.9% of total the fluorescence signature.

PARAFAC component intensities were normalized to the sum of component fluorescence intensities and expressed as their relative contribution to the total fluorescence. The identified component emission and excitation loadings were compared with the open fluorescence database⁶⁸ (https://openfluor.lablicate.com/) and with previous studies using PARAFAC analysis that are not included in the database. Components have been previously reported by 97 independent studies from openfluor database in various environments (Tucker Congruence Coefficient, TCC > 0.95) and in many studies in Arctic ecosystems (Table S4).

**Statistical analyses.** All statistical analysis were carried out using the R free software version 3.2.1⁶⁹. Bartlett's and Shapiro-Wilk tests were performed to, respectively, evaluate the variance homogeneity and the normal distribution of each variable of interest. Even log-transformed most variables were not normally distributed. A Principal Component Analysis (PCA) was performed using PARAFAC components and absorbance and fluorescence indices (biogeochemical data and soil properties were not included) was used to examine the differences in optical properties between thermal layers and layer natures. Non-parametric analyses of variance using Kruskal-Wallis (defined by its *chi-squared H and P value*) and the multiple pairwise-comparisons using pairwise Wilcoxon rank sum test (defined by its *p value*) were used to determine significant differences in concentrations and yields of DOC, DTN, NH₄⁺, NO₃⁻ and optical properties among thermal layers (active layer vs. permafrost), and between soil layer types (organic vs. mineral). Wilcoxon rank sum test tests (defined by *W and P value*) were used to assess the significance of difference between two specific groups. Before performing the PCA auto scaling was conducted on all variables with the FactoMineR package in order to decrease the leverage of high values. Regression analyses were employed to determine relationships between concentrations and yields of DOC, DTN, NH₄⁺, NO₃⁻, optical properties and soil properties.

**Canadian Permafrost landscape DOC pool estimates.** We estimated the DOC pool in permafrost soils at the pan-Canadian scale using the updated version of the Northern Circumpolar Soil Carbon Database (NCSCDv2)⁷⁰ at a spatial resolution of 1 km. NCSCDv2 provides information on SOC storage stratified both by depth layers and soil types. In the database, the class of Gelifol (permafrost soils) is subdivided in 3 suborders: Histel, Turbel, Orthel. We classified the soil type of each borehole as "organic permafrost soils" (Histels) or "mineral permafrost soils" (Orthels and Turbels). Soils were defined as "organic permafrost soils" if the organic surface layer was thicker than 40 cm. We defined soils as "mineral per-mafrost soils" for both Turbels (cryoturbuted soils) and Orthels (non-cryotur-bated), if the organic layer was less than 40 cm thick⁵². Non-permafrost soils were not taken into account in any calculations.

To create maps of soil DOC, we multiplied maps of soil carbon by the median yield of water-extractable DOC per soil type and at three depth layers (0-1 m, 1-2 m, 2-3 m). We did not calculate different yields between the active layer and permafrost layers. We examined the changes in yield between soils and depth layers using Kruskal-Wallis and Mann-Whitney test. We used quadrature to combine uncertainty of DOC yields with the 14.5% uncertainty in soil carbon from the NCSCD. We used the median and the 25th to 75th percentiles as the representative DOC yield and range of uncertainties, respectively.

## Data availability

The authors declare that all data supporting the findings of this study other than the excitation-emission matrices (EEMs) and the PARAFAC parameters are available within the paper and its supplementary information and supplementary data files. The EEMs and the PARAFAC code that generated the 6-component model are available from the corresponding author upon request. All data associated with this study will shortly be posted on the Polar Data Catalogue metadata website (https://www.polardata.ca) and Nordicana-D, the online data report series archived by the Centre for Northern Studies (http://www.cen.culaval.ca/nordicadan/en_index.aspx).

Received: 14 July 2020; Accepted: 18 August 2020;
Published online: 09 September 2020

## References

1. Chadburn, S. E. et al. An observation-based constraint on permafrost loss as a function of global warming. *Nat. Clim. Change* **7**, 340–344 (2017).

2. Turetsky, M. R. et al. Carbon release through abrupt permafrost thaw. *Nat. Geosci.* **13**, 138–143 (2020).

3. Hugelius, G. et al. Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps. *Biogeosciences* **11**, 6573–6593 (2014).
4. McGuire, A. D. et al. An analysis of the carbon balance of the Arctic Basin from 1997 to 2006. *Tellus B* **62**, 455–474 (2010).

5. Drake, T. W., Wickland, K. P., Spencer, R. G. M., McKnight, D. M. & Striegl, R. G. Ancient low-molecular-weight organic acids in permafrost fuel rapid carbon dioxide production upon thaw. *Proc. Natl Acad. Sci. USA* **112**, 13946–13951 (2015).

6. Wickland, K. P. et al. Dissolved organic carbon and nitrogen release from boreal Holocene permafrost and seasonally frozen soils of Alaska. *Environ. Res. Lett.* **13**, 065011 (2018).

7. Selvam, B. P. et al. Degradation potentials of dissolved organic carbon (DOC) from thawed permafrost soils. *Sci. Rep.* **7**, 45811 (2017).

8. Ward, C. P. & Cory, R. M. Chemical composition of dissolved organic matter draining permafrost thaw. *Geochim. Et Cosmochim. Acta* **167**, 63–79 (2015).

9. Metcalfe, D. B. et al. Patchy field sampling biases understanding of climate change impacts across the Arctic. *Nat. Ecol. Evol.* **2**, 1443–1448 (2018).

10. Mann, P. J. et al. Utilization of ancient permafrost carbon in headwaters of Arctic fluvial networks. *Nat. Commun.* **6**, 7856 (2015).

11. Harms, T. K., Abbott, B. W. & Jones, J. B. Thermo-erosion gullies increase nitrogen available for hydrologic export. *Biogeochemistry* **117**, 299–311 (2014).

12. Keuper, F. et al. Experimentally increased nutrient availability at the permafrost thaw front selectively enhances biomass production of deep-rooting subarctic peatland species. *Glob. Change Biol.* **23**, 4257–4266 (2017).

13. Gough, L. et al. Effects of long-term nutrient additions on Arctic tundra, stream, and lake ecosystems: beyond NPP. *Oecologia* **182**, 653–665 (2016).

14. Chen, L. et al. Nitrogen availability regulates topsoil carbon dynamics after permafrost thaw by altering microbial metabolic efficiency. *Nat. Commun.* **9**, 3951 (2018).

15. Mack, M. C., Schuur, E. A. G., Bret-Harte, M. S., Shaver, G. R. & Chapin, F. S. Ecosystem carbon storage in Arctic tundra reduced by long-term nutrient fertilization. *Nature* **431**, 440–443 (2004).

16. Marschner, B. & Kalbitz, K. Controls of bioavailability and biodegradability of dissolved organic matter in soils. *Geoderma* **113**, 211–235 (2003).

17. Kalbitz, K., Schwesig, D., Rethemeyer, J. & Matzner, E. Stabilization of dissolved organic matter by sorption to the mineral soil. *Soil Biol. Biochem.* **37**, 1319–1331 (2005).

18. Mann, P. J. et al. Controls on the composition and lability of dissolved organic matter in Siberia's Kolyma River basin. *J. Geophys. Res.* **117**, G01028 (2012).

19. Vonk, J. E. et al. Biodegradability of dissolved organic carbon in permafrost soils and aquatic systems: a meta-analysis. *Biogeoosciences* **12**, 6915–6930 (2015).

20. Weishaar, J. L. et al. Evaluation of specific ultraviolet absorbance as an indicator of the chemical composition and reactivity of dissolved organic carbon. *Environ. Sci. Technol.* **37**, 4702–4708 (2003).

21. Murphy, K. R., Stedmon, C. A., Graeber, D. & Bro, R. Fluorescence spectroscopy and multi-way techniques. *PARAFAC. Anal. Methods* **5**, 6557–6566 (2013).

22. Vincent, W. F., Lemay, M. & Allard, M. Arctic permafrost landscapes in transition: towards an integrated Earth system approach. *Arct. Sci.* **3**, 39–64 (2017).

23. Schädel, C. et al. Circumpolar assessment of permafrost C quality and its vulnerability over time using long-term incubation data. *Glob. Change Biol.* **20**, 641–652 (2014).

24. Hodgkins, S. B. et al. Changes in peat chemistry associated with permafrost thaw increase greenhouse gas production. *Proc. Natl Acad. Sci. USA* **111**, 5819–5824 (2014).

25. Zhang, X. et al. Importance of lateral flux and its percolation depth on organic carbon export in Arctic tundra soil: Implications from a soil leaching experiment. *J. Geophys. Res.* **122**, 796–810 (2017).

26. Neff, J. C. & Hooper, D. U. Vegetation and climate controls on potential CO2, DOC and DON production in northern latitude soils. *Glob. Change Biol.* **8**, 872–884 (2002).

27. Reyes, F. R. & Lougheed, V. L. Rapid nutrient release from permafrost thaw in arctic aquatic ecosystems. *Arct., Antarct. Alp. Res.* **47**, 35–48 (2015).

28. Beermann, F. et al. Permafrost thaw and liberation of inorganic nitrogen in Eastern Siberia. *Permafr. Periglac. Process.* **28**, 605–618 (2017).

29. Keuper, F. et al. A frozen feast: thawing permafrost increases plant-available nitrogen in subarctic peatlands. *Glob. Change Biol.* **18**, 1998–2007 (2012).

30. Elberling, B., Christiansen, H. H. & Hansen, B. U. High nitrous oxide production from thawing permafrost. *Nat. Geosci.* **3**, 332–335 (2010).

31. Shur, Y., Hinkel, K. M. & Nelson, F. E. The transient layer: implications for geocryology and climate-change science. *Permafr. Periglac. Process.* **16**, 5–17 (2005).

32. Jansson, J. K. & Hofmockel, K. S. Soil microbiomes and climate change. *Nat. Rev. Microbiol.* **18**, 35–46 (2020).

33. Mann, P. J. et al. Pan-arctic trends in terrestrial dissolved organic matter from optical measurements. *Front. Earth Sci.* **4**, 25 (2016).

34. Fellman, J. B., Hood, E. & Spencer, R. G. Fluorescence spectroscopy opens new windows into dissolved organic matter dynamics in freshwater ecosystems: a review. *Limnol. Oceanogr.* **55**, 2452–2462 (2010).

35. Walker, S. A., Amon, R. M. W. & Stedmon, C. A. Variations in high-latitude riverine fluorescent dissolved organic matter: a comparison of large Arctic rivers. *J. Geophys. Res.* **118**, 1689–1702 (2013).

36. Wickland, K. P., Neff, J. C. & Aiken, G. R. Dissolved organic carbon in alaskan boreal forest: sources, chemical characteristics, and biodegradability. *Ecosystems* **10**, 1323–1340 (2007).

37. Raudina, T. V. et al. Dissolved organic carbon and major and trace elements in peat powerwater of sporadic, discontinuous, and continuous permafrost zones of western Siberia. *Biogeoosciences* **14**, 3561–3584 (2017).

38. Fouché, J., Lafrenière, M. J., Rutherford, K. & Lamoureux, S. Seasonal hydrology and permafrost disturbance impacts on dissolved organic matter composition in High Arctic headwater catchments. *Arct. Sci.* **3**, 378–405 (2017).

39. Olefeldt, D., Persson, A. & Turetsky, M. R. Influence of the permafrost boundary on dissolved organic matter characteristics in rivers within the Boreal and Taiga plains of western Canada. *Environ. Res. Lett.* **9**, 035005 (2014).

40. Farrell, M. et al. Oligopeptides represent a preferred source of organic N uptake: a global phenomenon? *Ecosystems* **16**, 133–145 (2013).

41. Textor, S. R., Wickland, K. P., Podgorski, D. C., Johnston, S. E. & Spencer, R. G. M. Dissolved organic carbon turnover in permafrost-influenced watersheds of interior Alaska: molecular insights and the priming effect. *Front. Earth Sci.* **7**, 275 (2019).

42. McDowell, W. H. et al. A comparison of methods to determine the biodegradable dissolved organic carbon from different terrestrial sources. *Soil Biol. Biochem.* **38**, 1933–1942 (2006).

43. Cory, R. M. & Kaplan, L. A. Biological lability of streamwater fluorescent dissolved organic matter. *Limnol. Oceanogr.* **57**, 1347–1360 (2012).

44. Balcarczyk, K. L., Jeremy, B. J. J., Jaffé, R. & Maie, N. Stream dissolved organic matter bioavailability and composition in watersheds underlain with discontinuous permafrost. *Biogeochemistry* **94**, 255–270 (2009).

45. Näsholm, T. et al. Boreal forest plants take up organic nitrogen. *Nature* **392**, 914–916 (1998).

46. Kendrick, M. R. et al. Linking permafrost thaw to shifting biogeochemistry and food web resources in an arctic river. *Glob. Change Biol.* **24**, 5738–5750 (2018).

47. Abbott, B. W., Larouche, J. R., Jones, J. B., Bowden, W. B. & Balser, A. W. Elevated dissolved organic carbon biodegradability from thawing and collapsing permafrost. *J. Geophys. Res.:* **119**, 2049–2063 (2014).

48. Plaza, C. et al. Direct observation of permafrost degradation and rapid soil carbon loss in tundra. *Nat. Geosci.* **12**, 627–631 (2019).

49. Snow, Water, Ice and Permafrost in the Arctic (SWIPA) 2017 |AMAP. http://www.amap.no/documents/doc/Snow-Water-Ice-and-Permafrost-in-the-Arctic-SWIPA-2017/1610 (2017).

50. Pautler, B. G. et al. Molecular characterization of dissolved organic matter in glacial ice: coupling natural abundance 1H NMR and fluorescence spectroscopy. *Environ. Sci. Technol.* **46**, 3753–3761 (2012).

51. ADAPT 2016. Cryostratigraphy, carbon and nitrogen content and 14C dating of permafrost cores from sites across the Canadian Arctic, v. 1.0 (2013-2014). *Nordica D25, https://doi.org/10.5885/45427AD-06F05740704B4CA3.*

52. Soil Survey Staff. *Soil taxonomy: A basic system of soil classification for making and interpreting soil surveys.* 2nd edition. (Natural Resources Conservation Service. U.S. Department of Agriculture. Handbook, 1999).

53. Paré, M. C. & Bedard-Haughn†, A. Surface soil organic matter qualities of three distinct Canadian arctic sites. *Arct. Antarct. Alp. Res.* **45**, 88–98 (2013).

54. Grewer, D. M., Lafrenière, M. J., Lamoureux, S. F. & Simpson, M. J. Redistribution of soil organic matter by permafrost disturbance in the Canadian High Arctic. *Biogeochemistry* **128**, 397–415 (2016).

55. Hugelius, G., Kuhry, P., Tarnocai, C. & Virtanen, T. Soil organic carbon pools in a periglacial landscape: a case study from the central Canadian Arctic. *Permafr. Periglac. Process.* **21**, 16–29 (2010).

56. Hu, C., Muller-Karger, F. E. & Zepp, R. G. Absorbance, absorption coefficient, and apparent quantum yield: a comment on common ambiguity in the use of these optical concepts. *Limnol. Oceanogr.* **47**, 1261–1267 (2002).

57. Mann, P. J. et al. Pan-arctic trends in terrestrial dissolved organic matter from optical measurements. *Front. Earth Sci.* **4**, 25 (2016).

58. Lehmann, J. & Kleber, M. The contentious nature of soil organic matter. *Nature* **528**, 60–68 (2015).

59. Helms, J. R. et al. Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and photobleaching of chromophoric dissolved organic matter. *Limnol. Oceanogr.* **53**, 955–969 (2008).

60. Frey, K. E., Sobczak, W. V., Mann, P. J. & 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. *Biogeoosciences* **13**, 2279–2290 (2016).

61. Knight, D. M. et al. Spectrofluorometric characterization of dissolved organic matter for indication of precursor organic material and aromaticity. *Limnol. Oceanogr.* **46**, 38–48 (2001).
62. Cory, R. M. & McKnight, D. M. Fluorescence spectroscopy reveals ubiquitous presence of oxidized and reduced quinones in dissolved organic matter. *Environ. Sci. Technol.* **39**, 8142–8149 (2005).

63. Parlanti, E., Wörz, K., Geoffroy, L. & Lamotte, M. Dissolved organic matter fluorescence spectroscopy as a tool to estimate biological activity in a coastal zone submitted to anthropogenic inputs. *Org. Geochem.* **31**, 1765–1781 (2000).

64. Ohno, T. Fluorescence inner-filtering correction for determining the humification index of dissolved organic matter. *Environ. Sci. Technol.* **36**, 742–746 (2002).

65. Murphy, K. R., Stedmon, C. A., Waite, T. D. & Ruiz, G. M. Distinguishing between terrestrial and autochthonous organic matter sources in marine environments using fluorescence spectroscopy. *Mar. Chem.* **108**, 40–58 (2008).

66. Stedmon, C. A. & Bro, R. Characterizing dissolved organic matter fluorescence with parallel factor analysis: a tutorial. *Limnol. Oceanogr.* **6**, 572–579 (2008).

67. Stedmon, C. A., Markager, S. & Bro, R. Tracing dissolved organic matter in aquatic environments using a new approach to fluorescence spectroscopy. *Mar. Chem.* **82**, 239–254 (2003).

68. Murphy, K. R., Stedmon, C. A., Wenig, P. & Bro, R. OpenFluor-an online spectral library of auto-fluorescence by organic compounds in the environment. *Anal. Methods* **6**, 658–661 (2014).

69. R Core Team, R. F. for S. C. R: A Language and Environment for Statistical Computing. (2012).

70. Hugelius, G. et al. 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 (2013).

## Acknowledgements

We are grateful for field and lab assistance from Jonathan Roger, Steve Koziar, Daniel Lamhonwah, Megan Rueckwald, Yvette Chirinian, Tova Pinsky, Olivia Rodee, and Gil-lian Thiel. We thank Warwick Vincent, Mickaël Lemay, and all ADAPT contributors for providing research infrastructure and science support. This research was funded by the Natural Sciences and Engineering Research Council (NSERC) Discovery Frontiers ADAPT program, the ArcticNet Network of Centres of Excellence, and NSERC Dis-covery grants to M.J.L., S.F.L., and P.G. C.T.C. was financed by an Ontario Trillium Scholarship from the Ontario Ministry of Training, Colleges, and Universities. Logistics support was provided by the Polar Continental Shelf Program, Natural Resources Canada. We would also like to thank Dr Ashley Rudy, Dr Benjamin Amann, Dr Greg King, Celine Géran, and Romain Thouvenin for their support.

## Author contributions

M.J.L., S.F.L. and P.G. conceived of the research approach; C.T.C. collected soil cores with ADAPT collaborators, and C.T.C. and P.G. coordinated subsampling strategy and logistics; J.F. and C.T.C. performed all laboratory analyses; J.F. conducted the PARAFAC modeling; and J.F. drafted the manuscript with all co-authors contributing to the data interpretation and to the final manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41467-020-18331-w.

Correspondence and requests for materials should be addressed to J.F.

Peer review information *Nature Communication* thanks the anonymous reviewers for their contributions to the peer review of this work.

Reprints and permission information is available at http://www.nature.com/reprints

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020