Utilization of ancient permafrost carbon in headwaters of Arctic fluvial networks.

Northern high-latitude rivers are major conduits of carbon from land to coastal seas and the Arctic Ocean. Arctic warming is promoting terrestrial permafrost thaw and shifting hydrologic flowpaths, leading to fluvial mobilization of ancient carbon stores. Here we describe (14)C and (13)C characteristics of dissolved organic carbon from fluvial networks across the Kolyma River Basin (Siberia), and isotopic changes during bioincubation experiments. Microbial communities utilized ancient carbon (11,300 to >50,000 (14)C years) in permafrost thaw waters and millennial-aged carbon (up to 10,000 (14)C years) across headwater streams. Microbial demand was supported by progressively younger ((14)C-enriched) carbon downstream through the network, with predominantly modern carbon pools subsidizing microorganisms in large rivers and main-stem waters. Permafrost acts as a significant and preferentially degradable source of bioavailable carbon in Arctic freshwaters, which is likely to increase as permafrost thaw intensifies causing positive climate feedbacks in response to on-going climate change.

Climate-induced Arctic warming has led to increased soil temperatures, causing a succession of changes associated with permafrost degradation and widespread ground collapse or thermokarst, including deepening of the seasonally thawed surface active layer12 and alterations to watershed hydrology34. These changes threaten to destabilize ancient (Pleistocene-aged) northern circumpolar terrestrial permafrost, which contains vast quantities of organic carbon (OC) comprising as much as 50% of global below-ground soil carbon stocks5. Consequently, soil OC losses of up to 220 PgC by 2100 are predicted across Arctic regions67, leading to large-scale mobilization of permafrost-derived OC into Arctic inland and coastal waters8.

The great majority of OC (>80%)9 supplied to coastal seas from Arctic rivers is in the form of dissolved OC (DOC) and thus readily available to support microbial carbon demand during transit through fluvial networks10. Growing evidence exists from temperate and sub-tropical regions111213 that older aged terrestrial carbon can be an important subsidy to microbial communities, yet the processing and fate of ancient Arctic permafrost OC throughout fluvial networks is currently unknown. Understanding the response of aquatic systems to future increases in permafrost-derived DOC is necessary for determining its effects upon global carbon cycling and predicting its export via high-latitude rivers.

Studies from major Arctic rivers show that higher-order fluvial systems contain predominantly young (14C-enriched) DOC14151617, suggesting that limited mobilization of permafrost has occurred to date or that it has been removed upstream. A number of studies focusing on lower-order and headwater streams provide evidence for ancient DOC mobilization1819. Elucidating the fate of permafrost-derived DOC is crucial for establishing if on-going thaw is currently impacting Arctic fluvial networks and in determining whether its turnover may result in a positive carbon feedback to climate change7.

Here, we examine the quantity, age and source of the OC present, and constrain the nature of OC supporting microbial carbon demand with increasing water residence time through fluvial networks (higher order in the stream network). We base this work in the Kolyma River Basin, Siberia, which represents the largest global watershed (ca. 650,000 km2) completely underlain by continuous permafrost. The majority of this permafrost is comprised of frozen Yedoma, organic-rich (1–5% carbon by mass) Pleistocene-aged deposits that can contain 50–90% ice content520. Yedoma–ice complexes measured in the Kolyma Basin suggest that permafrost carbon by mass concentrations are on the order of 1.5±1.4%, with an average carbon inventory amounting to 14±8 kg m−3 (refs 21, 22, 23). Thawing Yedoma-derived OC is thought to be highly susceptible to biological degradation in fluvial networks1924, and has been proposed as a key feedback upon global climate20.

We first assess the initial age and bioavailability of DOC in streams and rivers throughout the basin using measurements of biologically degradable DOC (using DOC loss) and 14C-DOC measurements. We then determine the source and age of DOC subsidizing microbial demand throughout the network using simultaneous measurements of stable and radiocarbon signatures of DOC (13C and 14C) during short-term (28 days) bioincubation experiments. We show that a preferential loss of permafrost-derived DOC occurs with increasing water residence and inland water transit time, despite similar initially modern DOC ages. Our results demonstrate that microbial communities selectively remove ancient OC during water transit throughout the fluvial network, with a potential to rapidly act as a positive feedback upon climate change with increasing permafrost thaw.

Results

Initial 14C-DOC age and proportion bioavailable DOC

We conducted 81 individual incubation experiments on samples from 19 streams and rivers throughout the Kolyma fluvial network, grouping by stream size and type (Fig. 1; Supplementary Table 1). Thaw streams, directly draining Yedoma outcrops, contained extremely high concentrations of DOC (10,939±1,278 μM) that are highly depleted in 14C (Δ14C−883±41‰; Fig. 2a). A large proportion of thaw stream DOC was utilized by microorganisms (DOCloss) during short-term (28 days, 20 °C) incubations (47±8%), confirming the high bioavailability of permafrost-derived carbon in fluvial systems1924 (Table 1; Supplementary Table 2). Mean DOC concentrations declined moving downstream through the fluvial network, yet proportions of DOCloss were similar (15–21%) among all other stream and river waters (analysis of variance P>0.05, Table 1; Supplementary Table 2). Erosion-impacted streams that receive greater amounts of permafrost soil inputs via mechanical erosion of the stream banks contain 14C-depleted DOC (−214±145‰; Fig. 2a), highlighting the strong potential for thermokarst processes to sporadically deliver aged OC to fluvial networks. By contrast, all other stream and river waters contained modern initial bulk DOC pools (Δ14C>−50‰; Fig. 2a), consistent with the majority of available radiocarbon data for DOC from major Arctic Rivers141518.

Figure 1

Study site and sample locations.

Location and type of waters collected throughout the fluvial network. Cherskiy and site of the Northeast science station is marked in white. Individual site latitude and longitude information is provided in Supplementary Table 1.

Figure 2

Downstream Δ14C DOC values.

Δ14C values of (a) initial DOC, and (b) utilized DOC supporting microbial demand (DOCloss) during bioincubations. Δ14C values of utilized DOC were calculated as outlined in equation (1) (see Methods). Table 1 displays mean values for initial and utilized DOC across sites and sample numbers for each site type.

Table 1

DOC and isotope characteristics throughout the fluvial network.

	n	Initial DOC (μM)	SUVA254 (l mgC−1 m−1)	DOCloss (%)	Init. Δ14C (‰)	OCloss Δ14C (‰)	Mean OCloss age (14C years)	Permafrost (%)	Contemporary (%)	In situ (%)
Yedoma thaw	5	10,939±1,278	1.63±0.19	47.2±7.6	−884±41	−873±48	16,576	97.2±0.8	1.7±0.5	1.1±0.5
Erosion streams	3	2,503±518	3.25±0.44	17.5±4.9	−214±145	−405±185	4,171	42.8±20.7	56.4±20.7	0.8±0.4
Streams	25	1,691±87.6	3.44±0.11	13.3±2.0	51±6	−21±50	170	13.0±4.0	85.9±4.0	1.1±0.2
Minor tributaries	7	766±110	3.56±0.19	21.0±5.5	25±13	26±88	Modern	5.7±3.5	93.7±3.6	0.6±0.2
Major tributaries	8	511±68	4.03±0.08	12.2±1.8	40±3	68±32	Modern	0.6±0.1	73.1±8.4	26.3±9.9
Main stem	6	425±46	3.67±0.10	14.6±1.8	22±9	74±42	Modern	0.7±0.1	88.6±2.8	10.7±4.4

DOC, dissolved organic carbon; OC, organic carbon.

Mean initial DOC concentration, specific ultraviolet absorbance at 254 nm (SUVA254), percent DOC loss over 28 days (DOCloss), calculated Δ14C value and mean radiocarbon age of DOCloss (using equation (1) in Methods). The sample size at each site type is provided (n) and individual data are detailed in Supplementary Table 2. Mean percentage (±s.e.m.) contribution of permafrost, contemporary and in situ-derived DOC to DOCloss determined using the dual-carbon isotope mixing model.

Radiocarbon age of DOC supporting microbial demand

Microbial communities in Yedoma thaw stream waters consistently utilized highly 14C-depleted DOC (−750 to −1,000‰), indicating that Pleistocene-age OC (11,300 to >50,000 14C years) is intrinsically biologically labile and immediately available for microbial metabolism after mobilization to freshwaters (Fig. 2b). Microbial demand in erosion-impacted streams was also supported by highly aged DOC (Δ14C, −405±185‰, Table 1), demonstrating that microbial communities in waters containing mixed modern and ancient OC pools can preferentially utilize an older OC fraction (up to 10,000 14C years). Furthermore, despite similar modern initial bulk DOC ages among other stream and river waters, Δ14C values of the DOC utilized by microorganisms across these sites (that is, the fraction consumed during bioincubations) varied considerably (Δ14C+300 to −630‰; Fig. 2b), indicating a disconnect between bulk DOC age and the age of DOC fuelling microbial demand1113. Modern (14C-enriched) DOC periodically utilized by microorganisms through the river network may have been delivered from soil pore waters that can contain high values of Δ14C-DOC25 or may be caused by the upstream preferential loss of aged OC.

Source and radiocarbon age of DOC supporting microorganisms

We used a dual-isotope approach26 to examine the source and age of OC subsidizing DOCloss throughout the stream network. We partitioned source contributions that most likely explained utilized Δ14C and δ13C values of DOCloss, assuming three dominant sources of DOC (contemporary terrestrial OC, permafrost-derived terrestrial OC and internally produced in situ OC) (Fig. 3a).

Figure 3

Distribution of 14C and 13C isotopes and end-member contributions across the fluvial network.

(a) Initial (open squares) and utilized (black circles) Δ14C and δ13C values of DOC. Isotopic ranges of three contributing carbon sources are shown in red boxes (see Methods). (b) Mean percent (±s.e.m.) contribution of each carbon source to the DOC utilized over bioincubations (DOCloss) in waters from differing site types.

Microbial demand in Yedoma thaw waters was almost entirely (97.2±0.8%) supported by permafrost-derived OC, demonstrating that ancient DOC was highly bioavailable when present. Downstream, permafrost subsidized greater proportions of DOCloss in higher-order erosion-impacted (42.8±20.7%) and other small streams (13.0±4.0%) as compared with minor tributary (5.7±3.5%), major tributary (0.6±0.1%) and main-stem waters (0.7±0.1%) (Table 1). This landscape pattern of declining permafrost contribution to DOCloss through the fluvial network proved consistent across both years of investigation (2012 and 2013), showing that microbial demand in thaw streams and headwaters is consistently subsidized by significant proportions of permafrost OC. The fact that DOC character and age altered with increasing water residence time does not simply reflect catchment variability, as evidenced by identical patterns of DOC modernization in sub-basins within the larger Kolyma watershed. Furthermore, most sites drained similar soil and vegetation types, and spanned little elevation change within the Kolyma lowlands. Utilized Δ14C of DOCloss increased exponentially (r2=0.72, P<0.001, Supplementary Fig. 1) relative to initial Δ14C values, suggesting that older OC fractions were rapidly and preferentially mineralized when present in bulk DOC (initial Δ14C<−20‰). Terrestrially derived contemporary OC contributed towards a minor proportion of DOCloss in thaw waters (1.7±0.5%), but comprised the main OC source (>55%; Table 1) across stream and river waters, emphasizing the dominance of surface soil–water dynamics in determining bulk DOC character and thus biogeochemical processing and fate (Fig. 3b).

Environmental and compositional effects on 14C utilization

Seasonal and spatial trends in OC bioavailability across fluvial networks can be controlled by variations in water temperature, inorganic nutrient availability and OC character2427282930. To examine whether downstream variations in water temperature and inorganic nutrient concentrations influenced the fraction and age of DOC utilized, we conducted parallel experiments incubating waters at 4 °C and under nutrient-replete conditions. Total DOCloss was lower at colder incubation temperatures (mean Q10=1.21±0.07), but was not significantly different with respect to the mean Δ14C value of DOCloss (paired t-test, P>0.1; Fig. 4a). Bioenergetic constraints from lower water temperatures therefore did not consistently alter the DOC age or fraction utilized, implying that selective DOC loss continues during colder fall and winter periods, when greater contributions of older DOC could be expected1718. Similarly, inorganic nutrient concentration did not significantly influence Δ14C values of OCloss (paired t-test, P>0.1; Fig. 4b), indicating that stoichiometric differences cannot account for the observed downstream trends. Individual water samples displaying differences in the age of DOC utilized at colder temperatures, or under nutrient-replete conditions, all resulted in the loss of an older (14C-depleted) DOC fraction. This implies that our results represent conservative estimates of permafrost OC loss along the aquatic continuum under our experimental conditions (Fig. 3a,b).

Figure 4

Environmental and OC compositional controls on the age of carbon supporting freshwater microbial demand.

Δ14C values of DOC utilized (±propagated s.e.m.) in parallel incubations conducted (a) at 20 and 4 °C (n=19) and (b) under ambient and nutrient-enriched conditions (n=11). Dashed line represents 1:1 line. (c) Δ14C of mean DOC utilized across sites in relation to initial DOC aromaticity (inferred from SUVA254, see text for details). Decreasing SUVA254 indicates declining DOC aromaticity34.

Discussion

Together, our findings indicate that terrestrial permafrost OC export is already underway and actively contributing to DOC turnover in Arctic fluvial networks. The presence of permafrost OC has previously gone undetected in fluvial networks due to the preferential loss of permafrost OC removing headwater signatures dominated by permafrost inputs31, combined with the historical dominance of field sampling on larger main-stem rivers1532. In addition, the persistence of modern DOC sourced from vegetation and surface soils causes permafrost OC dynamics to be effectively masked when examining bulk initial OC alone33.

The composition of DOC in Arctic stream and river networks is expected to shift under future climate change scenarios12. Permafrost thaw is promoting greater groundwater and sub-surface flow from Arctic and sub-arctic watersheds, likely resulting in the export of DOC characterized by lower-molecular weight compounds and a lower degree of aromaticity34172731. Across our sites, we observed a linear increase in the proportion of bioavailable carbon (DOCloss) with decreasing DOC aromaticity as evidenced by lower SUVA254 values3435 (r2=0.50, P<0.001; Supplementary Fig. 2), suggesting that increasing future permafrost contributions may promote enhanced microbial metabolism within Arctic fluvial networks. Furthermore, the microbial utilization of older DOC fractions was observed with decreasing OC aromaticity (Fig. 4c) and increasing DOCloss (r2=0.37, P<0.001), suggesting that DOC compositional changes will result in greater ancient carbon turnover. Heterotrophic metabolism and respiration is dependent on low-molecular weight DOC from thawing permafrost soils2428. The lower aromatic content and molecular weight of permafrost-derived DOC may efficiently support bacterial production and favour higher carbon-use efficiencies36.

The microbial uptake of highly biologically available permafrost carbon can result in the preferential loss of ancient OC in aquatic networks, and the selective export of modern DOC to the ocean. Similar trends in the processing of permafrost-derived DOC among fluvial networks of other major arctic rivers appear likely. All arctic watersheds are warming6 and exhibit similar trends and controls upon DOC export32. Freshwater discharge is also strengthening from at least the Eurasian arctic watersheds3738. Yedoma permafrost in Alaska has also recently been shown to be highly bioavailable in good agreement with studies in Siberia19243940. Additional mechanisms, such as the selective burial of older aged OC, changing microbial community composition or preferential photodegradation of ancient OC may also play an important or complementary role in controlling permafrost turnover in northern high-latitude systems101241. Photodegradation has been shown to influence surface water DOC processing in the Arctic41, yet it is currently unclear how sunlight affects the bulk age of DOC. Recognizing that modern bulk DOC signatures in Arctic rivers may disguise a rapidly cycling old C fraction, future efforts to detect the impacts of permafrost thaw should focus on headwater systems where the rapid removal and preferential biological utilization of ancient permafrost-derived DOC appears prevalent. As permafrost thaw accelerates168, it seems apparent that this will be accompanied by an increase in the amount of bioavailable DOC in aquatic networks, with aged DOC accounting for increased subsidies to aquatic microorganisms ultimately fuelling a positive feedback on climate.

Methods

Study area and sample collection

Surface water samples were collected directly from streams and rivers on several expeditions in 2012 and 2013 (Supplementary Table 1). Sampling concentrated during the months of August to October, when surface active layers are deepest and thus maximum export of permafrost-derived OC to inland waters is expected. Additional sampling soon after the spring ice-out period (June) provided insights during a period of minimal surface active layer thaw. We grouped sites by stream width, separating streams (<1–5 m), minor (>5 to ≤200 m) and major tributaries (>200 m), and the Kolyma main stem (Fig. 1). A small proportion of highly turbid (>100 nephelometric turbidity units (NTUs)) upland streams visibly affected by active mechanical soil erosion were examined independently (erosion streams) from unaffected streams (<100 NTU), thus providing more robust source determination26 (see the Isotope mixing model section). These streams were always small (<3 m width) and situated within active thermokarst zones with mixed vegetation covered and exposed yedoma banks.

Stream temperature was measured in the field (YSI multi-meter) and sample turbidity immediately upon return to the laboratory (Hach 2100Q). Waters were filtered through pre-combusted glass fibre filters (0.7 μm) on the day of collection using a combusted glass filtration system. After filtration, waters were immediately transferred into pre-combusted amber glass vials (40 ml vials, ∼35 ml water) leaving a headspace. Triplicate aliquots for initial DOC concentration were immediately acidified (trace element grade HCl, ≤pH 2) and stored until analysis. Aliquots for δ13C and Δ14C were immediately stored frozen (−20 °C) in pre-leached, acid-cleaned HDPE bottles.

DOC and bioincubations

DOC was calculated as the mean of between three and five injections on a Shimadzu TOC-V where the coefficient of variance across measurements was <2%. Triplicate bioincubation DOC measurements were typically within <5% of each other. The bioavailability of DOC (DOCloss) was determined as the difference in triplicate measures of DOC before and after a 28-day laboratory incubation at 20±0.5 °C or 4±1 °C containing the indigenous microbial community from each site. The absorbance of water at 254 nm was measured (Shimadzu UV-1800) and the specific ultraviolet absorbance (SUVA254) determined as a proxy for aromaticity34. After bioincubation, waters were refiltered as above, and samples collected as before for DOC, δ13C and Δ14C analyses. Parallel bioavailability incubations conducted with nutrient-amended waters were supplemented with inorganic nitrate and phosphate (KNO3/KH2PO4) according to Redfield ratio, resulting in nutrient-replete incubations in relation to ambient site nutrient concentrations. We used analysis of variance to test between site means and paired t-tests to compare parallel incubations results. All tests including Univariate regression and exponential fits were conducted in SPSS21 (IBM).

Stable and radiocarbon isotope analyses

δ13C analyses were conducted at the University of California, Davis Stable Isotope Facility and 14C analyses at the Eidgenössiche Technische Hochschule (ETH) Zürich. 13C-DOC samples were analysed using an O.I. Analytical Model 1010 TOC analyzer (precision of ±0.2‰) interfaced to a PDZ Europa 20–20 IRMS (Sercon Ltd). 13C-DOC measurements were calibrated against the δ13C values of KHP and IHSS Suwannee River humic acid in Milli-Q. Waters for 14C-DOC analyses were freeze-dried (Christ Alpha 2-4, LSC with a low-carbon vacuum hybrid pump, Vacubrand RC-6; Martin Christ, Labex Instrument AB, Sweden) directly in pre-combusted (850 °C/5 h) quartz tubes. Samples were acidified under HCl vapours to remove carbonates and flame sealed with pre-combusted CuO under vacuum. CO2 was cryogenically captured and quantified (∼30 μg carbon) before measurement using a microscale radiocarbon dating system (MICADAS; http://www.ams.ethz.ch/instruments/micadas) and gas-feeding system4243. Combusted NIST SRM 4990C oxalic acid was used as a standard for normalization, and blanks were determined using radiocarbon-free CO2 both at a concentration of 5% CO2 in He. The modern oxalic acid standard was measured to better than 1% relative error and until samples were fully consumed.

Radiocarbon contents are reported as Δ14C (‰) and 14C age. All radiocarbon values are corrected for procedural blanks with the errors propagated. The isotopic signature (Δ14C and δ13C) of the utilized DOC fraction (DOCloss) was calculated using triplicate measurements of initial DOC (DOCinitial) and final DOC concentrations (DOCfinal) alongside the associated change in isotopic composition using simple mass balance (equation (1)). Individual errors associated with multiple isotope measurements were propagated to assess error on δ13C and Δ14Cutilized.

Isotopic mixing model

Feasible contributions of three end-member sources (permafrost DOC, contemporary DOC and in situ autochthonous DOC) to DOCloss at sites were calculated using a dual-carbon isotope mixing model solved using a Monte Carlo simulation approach. Monte Carlo simulations were conducted using the MIXSIAR package26 v2.1.2 (https://github.com/brianstock/MixSIAR/releases) within the R programming environment (http://www.R-project.org). MixSIAR is a Bayesian mixing model44 that allows the contributions of source end-members to isotope data to be estimated, taking into account the uncertainty in source values45 and model error46 (residual and process errors). A dual-isotope approach (14C/13C) was used to increase the accuracy of contribution estimates and to solve for three potential end-member sources.

We defined two separate terrestrial end-members to differentiate between contemporary OC pools (vegetation, surface soils and surface active layer) and permafrost-derived OC that is delivered to freshwaters via ground collapse, or thermokarst processes, bank erosion or deep groundwater flow. An in situ-derived DOC source from aquatic primary producers (for example, periphytic algae) was also included. End-member values of each DOC source were determined as outlined below with final mean values and ranges defined as follows: permafrost δ13C −26.3±1.3‰, Δ14C −940±84‰; contemporary δ13C −28.5±2.0‰, Δ14C 71.6±170‰; in situ autochthonous δ13C −19.8±3.0‰, Δ14C −78.0±170‰. Error estimates on end-member sources incorporated the entire variance in source values, enabling MixSIAR to fit the most viable end-member source to the model fit26. To account for potential fractionation of δ13C, the MIXSIAR package was run including discrimination factors set to ±1‰. No fractionation was allowed for Δ14C values.

Each experiment was treated as an individual ‘random' event, and site type (yedoma thaw, erosion stream, stream, minor and major tributary and main-stem waters) was included as a ‘random' effect47, which improved model fit (determined as a reduction in the Akaike information criterion). Site type (σsite=7.1) explained a greater variation in source than individual information alone (σind=1.8), indicating that the total variation in the source of DOC fuelling microbial demand was driven by location along the fluvial network. Here we present posterior mean contribution estimates (±s.e.m) among sites (Table 1; Supplementary Table 2).

Isotopic end-member determination

The contemporary δ13C (−28.5±2.0‰; n=14) and Δ14C (71.6±150‰; n=21) was derived from literature18 and unpublished data (http://arcticgreatrivers.org/data, Vonk ;Supplementary Table 3) collected during the months of May and June from a range of small streams and rivers from the Kolyma River Basin. Data were limited to the early spring period, when water flow is restricted to surface horizons and thus most likely contains a topsoil and vegetation signature. A larger error estimate was applied to the Δ14C end-member values to incorporate soil pore water-derived DOC that can be highly 14C enriched25.

The δ13C value of the permafrost end-member (−26.3±1.3‰; n=374) is taken from a recent review of Yedoma deposits across Siberia48 and was also recently used in a similar dual-carbon isotope model to determine relative contributions of Yedoma to the East Siberian Arctic Shelf49. The value closely aligns with other δ13C DOC values (−25.7±0.3‰; n=7) measured in Yedoma thaw waters19 and soils (−25.1±0.3‰; n=2) from the Kolyma River Basin (>2 m)50. The Δ14C value for permafrost OC was taken from a synthesis of Yedoma data as detailed in Vonk et al. (2012; Supplementary Table S4 and references herein)49. This value (−940±84‰; n=300) closely aligns with Yedoma permafrost thaw water DOC measured in the Kolyma River Basin1933 (−944±23‰; n=6). By using a maximum estimate for the old carbon permafrost end-member, we ensure conservative estimates of old C supply to stream and river waters.

The δ13C value of aquatic primary producers (for example, periphytic algae) is highly variable and difficult to assess, with literature values varying from −40 to −10‰51. δ13C end-member values were estimated using dissolved inorganic carbon measurements from the Kolyma River (Raymond et al., unpublished) (−7.5 to −7.2‰) and Lena River (−4.0 to −3.5‰), and periphytic OC δ13C (−23.0 to −16.6‰) calculated using an established relationship between δ13C of DIC and of epilithic algae from across 70 different streams and rivers spanning arctic and temperate environments52. This linear relationship was derived from direct measurements of herbivore δ13C and DIC δ13C values52 and demonstrated that a consistent fractionation (for example, 19 or 20‰ assumed fixed photosynthetic fractionation) did not accurately constrain δ13C of algae from DIC measurements in lotic systems. This resulted in a relative fractionation for this study ranging from 13.1 to 15.5‰. Our calculated δ13C for algae closely align with values measured in a single plankton bloom on the Pantileikha River within the Kolyma Basin of −23.3‰ (ref. 53) and of the end-members calculated for marine OC in Vonk et al.49 of −24‰. Both values lie within our end-member estimates once the model discrimination factor of up to 1‰ has been applied for δ13C fractionation. Phytoplankton should reflect the current Δ14C of DIC in waters. Δ14C of DIC from the Kolyma River ranged from −77.0 to −79.0‰ (Raymond et al., unpublished), and was therefore adopted as the mean end-member value in the mixing model.

Additional information

How to cite this article: Mann, P. J. et al. Utilization of ancient permafrost carbon in headwaters of Arctic fluvial networks. Nat. Commun. 6:7856 doi: 10.1038/ncomms8856 (2015).