# Impact of Permafrost Thaw and Climate Warming on Riverine Export Fluxes of Carbon, Nutrients and Metals in Western Siberia

Oleg S Pokrovsky, Rinat M Manasypov, Sergey G Kopysov, Ivan V Krickov, Liudmila S Shirokova, Sergey V Loiko, Artem G Lim, Larisa G Kolesnichenko, Sergey N Vorobeyev, Sergey N Kirpotin

► To cite this version:

Oleg S Pokrovsky, Rinat M Manasypov, Sergey G Kopysov, Ivan V Krickov, Liudmila S Shirokova, et al.: Impact of Permafrost Thaw and Climate Warming on Riverine Export Fluxes of Carbon, Nutrients and Metals in Western Siberia. Water, 2020, 12 (6), pp.1817. 10.3390/w12061817. hal-04849119

HAL Id: hal-04849119

https://ut3-toulouseinp.hal.science/hal-04849119v1

Submitted on 19 Dec 2024

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

# Impact of Permafrost Thaw and Climate Warming on Riverine Export Fluxes of Carbon, Nutrients and Metals in Western Siberia

Oleg S. Pokrovsky<sup>1,2,3,*</sup>, Rinat M. Manasypov<sup>1 ID</sup>, Sergey G. Kopysov<sup>4</sup>, Ivan V. Krickov<sup>1 ID</sup>, Liudmila S. Shirokova<sup>2,3</sup>, Sergey V. Loiko<sup>1 ID</sup>, Artem G. Lim<sup>1 ID</sup>, Larisa G. Kolesnichenko<sup>1</sup>, Sergey N. Vorobyev<sup>1</sup> and Sergey N. Kirpotin<sup>1</sup>

Received: 30 April 2020; Accepted: 22 June 2020; Published: 24 June 2020

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

check for
updates

**Abstract:** The assessment of riverine fluxes of carbon, nutrients, and metals in surface waters of permafrost-affected regions is crucially important for constraining adequate models of ecosystem functioning under various climate change scenarios. In this regard, the largest permafrost peatland territory on the Earth, the Western Siberian Lowland (WSL) presents a unique opportunity of studying possible future changes in biogeochemical cycles because it lies within a south-north gradient of climate, vegetation, and permafrost that ranges from the permafrost-free boreal to the Arctic tundra with continuous permafrost at otherwise similar relief and bedrocks. By applying a "substituting space for time" scenario, the WSL south-north gradient may serve as a model for future changes due to permafrost boundary shift and climate warming. Here we measured export fluxes (yields) of dissolved organic carbon (DOC), major cations, macro- and micro-nutrients, and trace elements in 32 rivers, draining the WSL across a latitudinal transect from the permafrost-free to the continuous permafrost zone. We aimed at quantifying the impact of climate warming (water temperature rise and permafrost boundary shift) on DOC, nutrient and metal in rivers using a "substituting space for time" approach. We demonstrate that, contrary to common expectations, the climate warming and permafrost thaw in the WSL will likely decrease the riverine export of organic C and many elements. Based on the latitudinal pattern of riverine export, in the case of a northward shift in the permafrost zones, the DOC, P, N, Si, Fe, divalent heavy metals, trivalent and tetravalent hydrolysates are likely to decrease the yields by a factor of 2-5. The DIC, Ca, SO₄, Sr, Ba, Mo, and U are likely to increase their yields by a factor of 2-3. Moreover, B, Li, K, Rb, Cs, N-NO₃, Mg, Zn, As, Sb, Rb, and Cs may be weakly affected by the permafrost boundary migration (change of yield by a factor of 1.5 to 2.0). We conclude that modeling of C and element cycle in the Arctic and subarctic should be region-specific and that neglecting huge areas of permafrost peatlands might produce sizeable bias in our predictions of climate change impact.

**Keywords:** river flux; weathering; organic matter; permafrost; trace element; river

<footnote n="1">¹ BIO-GEO-CLIM Laboratory, Tomsk State University, 36 Lenina ave, 634050 Tomsk, Russia; rmmanassypov@gmail.com (R.M.M.); krickov_ivan@mail.ru (I.V.K.); s.loyko@yandex.ru (S.V.L.); lim_artyom@mail.ru (A.G.L.); klg77777@mail.ru (L.G.K.); soil@green.tsu.ru (S.N.V.); kirp@mail.tsu.ru (S.N.K.)</footnote>

<footnote n="2">² Geosciences and Environment Toulouse, UMR 5563 CNRS, 14 Avenue Edouard Belin, 31400 Toulouse, France; liudmila.shirokova@get.omp.eu</footnote>

<footnote n="3">³ N. Laverov Federal Center for Integrated Arctic Research, Russian Academy of Sciences, 23 Naberezhnaya Sev. Dviny, 163000 Arkhangelsk, Russia</footnote>

<footnote n="4">⁴ Institute of Monitoring of Climatic and Ecological Systems SB RAS, 10/3 Academichesky ave, 634055 Tomsk, Russia; wosypok@mail.ru</footnote>

<footnote>* Correspondence: oleg.pokrovsky@get.omp.eu</footnote>
# 1. Introduction

Arctic warming is anticipated to result in massive carbon (C) mobilization from permafrost peat to atmosphere, rivers and lakes, thereby potentially worsening global warming via greenhouse gases (GHG) emissions [1,2]. Permafrost peatlands cover roughly 2.8 million km² or 14% of permafrost-affected areas, mostly in Northern Eurasia (Bolshezemelskaya Tundra, $0.25 \times 10^6$ km²; western Siberia, $1.05 \times 10^6$ km²; Northern Siberia and Eastern Siberia lowlands, $0.5 \times 10^6$ km²) and contain a huge amount of highly vulnerable carbon in soil and surface waters [3]. Except for several regional studies of peatland lakes and small streams in Canada [4–6] and northern Sweden [7], the control factors, timing and reality of C and related elements release from soil and sediments are largely unknown, in part because element export by rivers across the permafrost peatlands is still poorly quantified.

One of the largest peatlands in the world is the Western Siberia Lowland (WSL) which exhibits a disproportionally high contribution to C storage and exchange fluxes [8] and presents a prominent exception to well-studied aquatic and soil ecosystems in mountainous territories of Northern America and Scandinavia. Specific features of the WSL are: (1) developed on flat, low runoff terrain with long water residence time in both lentic and lotic systems; (2) acts as important C stock in the form of organic-rich histosols (peat soils); (3) emits substantial amount of CO₂ and CH₄ to the atmosphere from the surface of inland waters; and (4) exhibits high dissolved organic carbon (DOC) concentration and low pH in waters in contact with peat soils. These factors determine rather unique and still poorly known aquatic communities in dystrophic to mesotrophic peatland waters. Due to the lack of nutrients, shallow photic layer and high sensitivity to water warming, aquatic ecosystems of frozen peatlands are highly vulnerable to the permafrost thaw and can respond unexpectedly to ongoing climate warming in terms in their C, nutrient, and metal storage and export fluxes as well biodiversity and organisms adaptation strategies.

Western Siberia is a key region for biogeochemical studies in the Arctic [2] as this region includes the continuous-discontinuous permafrost transition and experiences substantial thermokarst. As a result, this high-priority region is most susceptible to thaw-induced change in solute transport by rivers and its export to the ocean and atmosphere, notably via activation of deep underground flow [9] and exposure of large volumes of previously frozen peat and mineral soils [10–12]. Raised by incontestable Arctic amplification of overall climate change, the fate of carbon, nutrient and metal in Arctic rivers is at the forefront of field and modeling studies [13–20]. In particular, the Arctic Great Rivers Observatory (GRO) program assessed concentrations and export fluxes of 6 largest Arctic Rivers including the Ob River, draining sizeable part of the WSL [21,22]. However, widespread loss of hydrological monitoring in the beginning of 21st century was especially pronounced for small and medium-size rivers of the permafrost-affected part of the WSL [23]. At the same time, the Ob River is not suitable for modeling of possible changes in western Siberia because (1) it is strongly influenced by its largest permafrost-free tributary, the Irtysh River; and (2) it cannot be used for assessing the fluxes of northern (permafrost) zones as it integrates vast territory of variable permafrost coverage (20% in average) and landscape parameters. For these reasons, the study of small rivers at the WSL territory is more suitable for assessing both mechanisms of flux formation and its possible future changes.

Currently, a dominant paradigm is that riverine fluxes of C and inorganic nutrients are increasing in virtually all permafrost-affected rivers [14,24,25]. In particular, the increase in suspended versus dissolved transport of elements may be due to abrasion of lake shores and riverbanks. At the same time, enhanced groundwater input will lead to an increase in the transport of truly dissolved forms. Recently, following the pioneering work of Frey et al. [26–28], the concentrations of dissolved, particulate and colloidal carbon, nutrient and metals in WSL rivers, lakes and soil waters have been studied over a sizable latitudinal gradient [29–38]. These results allowed a first-order assessment of the consequences of climate warming and permafrost thaw on river water concentrations of dissolved and particulate forms of elements from western Siberia to the Arctic Ocean and C emission to the atmosphere. For this, a "substituting space for time" approach was employed, which postulates, in a broad context, that spatial phenomena which are observed today can be used to describe past and future events [39].
Thus, a concentration pattern of major and trace elements (TE) in WSL river suspended matter implies that, upon a progressive shift of the permafrost boundary northward, there will be a sizeable decrease in concentrations of alkalis, alkaline-earths, divalent heavy metals, and trivalent and tetravalent hydrolysates in northern rivers in currently discontinuous and continuous permafrost zones [31]. This decrease may be by a factor of 2–5 from the position of the minimum elemental concentration in sporadic to discontinuous permafrost zones relative to the continuous permafrost zone. Concerning the dissolved element concentration and potential transport in the WSL rivers, the following results were achieved implying the substituting space for time scenario. From the permafrost thaw perspective, the increase in depth of the active layer and connectivity of a river with underground water reservoirs may decrease the colloidal fraction (1 kDa–0.45 µm) of OC, Fe, Al and number of divalent metals as well as low-soluble trivalent and tetravalent hydrolysates in the sporadic and isolated permafrost zone as it becomes permafrost-free [30]. The forestation of wetlands and lake drainage may slightly diminish colloidal transport of DOC and metals at the expense of low molecular weight forms. Overall, given the significant role of seasonal and forestation effects on colloidal forms of OC and TE in WSL rivers, major changes in the speciation of riverine C and nutrients in the WSL may occur due to changes in vegetation rather than in temperature and precipitation [30]. A northward permafrost boundary shift with increase in air and water temperature may decrease or maintain, rather than increase of major nutrient (K, P, N, Si) and DOC concentrations in rivers draining through continuous permafrost zone [37].

All these predictions described above were made based solely on evolution of concentrations of soluble and particulate C, metals and nutrients in WSL rivers, without taking into account the hydrological flux. The latter could not be addressed until now due to lack of reliable information on river discharge over different seasons. The present work aims to fill this gap by quantifying the export fluxes of ~30 WSL rivers of various size, combined with new hydrological modeling of region and season-specific river runoff of the WSL territory. The chosen rivers encompass a large gradient of climate, vegetation and permafrost distribution and thus enables a substituting space for time scenario to provide a tentative prediction of possible changes in riverine fluxes for short- and long-term prospects. The present work is built on seasonally resolved sampling performed in 2015 and 2016, and incorporates thorough hydrological modeling to calculate the seasonal discharge. Building upon our previous studies of the WSL rivers [29–31,33,34,37], here we assess for the first time, element export fluxes (yields) across a large permafrost-affected territory and large number of rivers. The obtained elementary yields are essential for further modeling of biogeochemical cycles in the permafrost regions.

## 2. Study Site and Methods

### 2.1. *Rivers of Western Siberian Lowland, their Sampling and Analyses*

The 32 rivers of the Ob, Pur and Taz watersheds in the WSL were sampled (Figure 1). Detailed climatic, lithological and physio-geographic characteristics for the WSL are presented elsewhere [29,33,34]. The dominant lithology is Quaternary deposits (silt, clay, sand) overlaid by peat. The climate gradient of sampled rivers presents a decrease in mean annual air temperature (MAAT) from –0.5 °C in the south (Tomsk region) to –9.5 °C in the north (Arctic coast). Annual precipitation is fairly constant ranging between 550 mm in the south and 600 mm at the lower reaches of the Taz River. The river runoff ranges from 190 ± 30 mm y⁻¹ in the south to 300 ± 20 mm y⁻¹ in the north [40]. The distribution of the permafrost reflects the south-north MAAT gradient and changes from isolated and sporadic in the south to discontinuous and continuous in the north. In 2015, we sampled rivers in spring (18 May–25 June) and summer (25 July–21 August). In 2016, sampling was performed in spring (17 May–15 June 2016), summer (1–29 August 2016), and autumn (24 September–13 October 2016). We followed the change of seasons during our sampling campaign and moved from the south to the north in spring and from the north to the south in autumn thus collecting the river water at approximately the same period of the discharge. Note that more frequent sampling would be desirable to accurately evaluate the annual export flux, but rather harsh environment and logistical issues
constrained sampling. In contrast, route sampling is a common way to assess chemical weathering in extreme environments [41,42], and it is accepted that single sampling during high flow season provides the best agreement with time-series estimates [43].

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

**Figure 1.** The WSL river sampling points. The numbers on the map correspond to the following rivers: (1) Ob' (Pobeda), (2) Maliy Tatosh, (3) Chaya, (4) Vyalovka, (5) Ob' (Aleksandrovskoye), (6) Vakh, (7) Agan, (8) Tromyegan. (9) Vach-Yagun, (10) Vachinguriyagun, (11) Pintyr'yagun, (12) Kamgayakha, (13) Khatytayakha, (14) Pyakupur, (15) Lymbyd'yakha, (16) Apoku-Yakha, (17) Etu-Yakha, (18) Seryareyakha, (19) Purpe, (20) Aivasedapur, (21) Tydylyakha, (22) Yamsovey, (23) Pur, (24) Ngarka Khadyta-Yakha, (25) M. Kheyakha, (26) Nuny-Yakha, (27) Taz. The upper left insert is from Google Maps®.
The two sampled years contrasted in mean monthly and annual temperature and precipitation. The summer 2015 was warm (+0.3 °C above normal for July), while summer 2016 was cool (−3.1° below normal for July). The summer 2015 was wet (180% of normal precipitation for July), and the summer 2016 was less wet (but still 130% of normal precipitation for July). The winter precipitations of 2014–2015 and 2015–2016 were rather similar and sizably higher (ca. 140%) than the normal winter values. The normal values here are defined as between 1970 to 2000 based on the Roshydromet archives [44].

The surface water was collected in a polypropylene 1-L container. Samples were filtered through 0.45 µm cellulose acetate filters using a Nalgene 250 mL filter unit combined with a vacuum pump. All filtrations were run on site, in a protected environment, within 2 h of river water collection. Immediately after filtration, samples for DOC, DIC, major and trace elements were stored in the refrigerator during 1–2 months prior to the analyses, while the samples for nutrients were kept frozen. Dissolved (<0.45 µm) concentrations of nutrients, major and trace elements in WSL river waters were analyzed as described elsewhere [33,34,37,45]. Major anion concentrations (Cl⁻ and SO₄²⁻) were measured by ion chromatography (HPLC, Dionex ICS 2000) with an uncertainty of 2%. DOC and DIC were analyzed using a Shimadzu TOC 6000 with an uncertainty of 3–5% [46]. For all major and most trace elements, analyzed by ICP MS, the concentrations in the blanks were below analytical detection limits (≤0.1–1 ng/L for Cd, Ba, Y, Zr, Nb, rare earth elements (REEs), Hf, Pb, Th, U; 1 ng/L for Ga, Ge, Rb, Sr, Sb; ~10 ng/L for Ti, V, total P (P<sub>tot</sub>), Cr, Mn, Fe, Co, Ni, Cu, Zn, As). The international certified reference material SLRS-5 (Riverine Water Reference Material for Trace Metals) was used to validate the analysis. Further details of analytical uncertainties and detection limits for TE are provided elsewhere [30,34,45].

## 2.2. Hydrological Parameters, River Discharge Modeling, and Element Export Flux Calculations

The runoff was longitudinally modelled using average altitude of the watershed for gauged rivers [47]. For small and medium rivers draining palsa and polygonal bogs of the permafrost zone, we used empirical equations accounting for hydrological characteristics of these watersheds [48]. In this work, only open water period (May to October) was considered. The contribution of winter time export in large rivers typically does not exceed 10% of mean annual runoff [49]; moreover, small (<1000–10,000 km² watershed) WSL rivers freeze solid in winter [33]. We interpolated discharge by using watershed area change along the river mainstem. The river hydrographs were modelled using HBV Light program package (https://www.geo.uzh.ch/en/units/h2k/Services/HBV-Model.html) as described in refs [50–52]. For the gauged rivers, we modelled daily runoff based using mean water temperature and precipitation from adjacent meteostations. The data for reproducing the water discharge at key sites were taken from hydrological yearbooks and automatic information system of State monitoring of water bodies [53]. The daily air temperature and precipitation were taken from Russian Hydrometeocenter (URL: http://aisori.meteo.ru/ClimateR; [44]). The calibration of model parameters demonstrated that the HBV model adequately reproduces the seasonal dynamics of runoff (quality criteria of 0.75 to 0.90). The uncertainties of modelled discharge at ungauged rivers stem from several factors. First, there is uncertainties in characterizing watershed boundaries for the very flat, but poorly resolved, WSL territory. For small rivers, this uncertainty can amount to 30% due to lack of precise topographic information. For large rivers, it can be up to 7% due to difficulties in determination of exact position of the watershed divide on flat bog-lake landscapes. The second cause of uncertainty is a discrepancy between the model parameters of analogous rivers due to intrinsic differences in the conditions of runoff formation (up to 25%). Finally, the amount of precipitation obtained from the nearest meteostations might substantially differ from the actual precipitation at a given watershed, which is reflected, in particular, by low quality of the rain-flood modeling. The maximal overall uncertainty of daily discharges is determined by the most probable value of 30%. However, because we used two-month averaged discharges for calculating elementary yields, the real uncertainty of export fluxes is determined by the uncertainty of watershed delineations and ranges from 7 to 30%.
In order to calculate element export fluxes, we defined six latitudinal belts (56–58° N, 58–60° N, 60–62° N, 62–64° N, 64–66° N and 66–68° N). These latitude belts were selected based on: i) the permafrost map of the WSL, where the permafrost zones roughly follow latitude, and ii) necessity to integrate sufficient number of rivers in each permafrost zone and for statistical comparison. For statistical comparisons, we separated the years 2015 and 2016, because the autumn was sampled only in 2016. The number of rivers used for the latitudinal-average flux calculation ranged from 2 to 10 in 2016 and from 2 to 16 in 2015 for each latitudinal range.

# 3. Results

## 3.1. *Impact of the Watershed Size and Season on Element Export Fluxes*

The yields of DOC, representative major solutes (DIC, Ca, Mg) and nutrients (Si, P<sub>tot</sub>, K, Fe) are illustrated in Supplementary Figure S1. Flux magnitudes exhibited sizable variation, ranging over one (DOC, Si and Mg) to two (DIC, Ca, K, Fe and P<sub>tot</sub>) orders of magnitude. Generally, the variations were the highest for rivers having watersheds between 100 and 1000 km<sup>2</sup>, for both permafrost-affected and permafrost-free. The river watershed area (S<sub>watershed</sub>) exhibited rather weak control on elementary yields, as also illustrated by statistical treatment of individual seasons and full data set for both years (Supplementary Table S1). Overall, the watershed size exhibited more pronounced correlations with element yield in the permafrost-free zone compared to the permafrost zone. The correlations were quite low during spring flood but become pronounced during summer and autumn baseflow. The highest correlations were observed in summer, when Cl, SO<sub>4</sub>, Fe, Cu, Y, Mo, Sb, REEs and Th were positively linked to S<sub>watershed</sub> (R<sub>Spearman</sub> > 0.58, p < 0.01) in the permafrost-free zone.

Partial contribution of spring, summer, and autumn (2 months each) to the overall export of elements during the six-month open-water period (May to October) in 2016 demonstrated the dominant role of spring for DOC, Al and Fe (>60%), for both permafrost-free and permafrost-affected rivers (Figure 2). In contrast, soluble highly mobile elements (DIC, Na, K, Ca, Mg, Si) had less than 40% and 50% contribution of spring in the permafrost-free and permafrost zone, respectively, and a 2–3 times higher contribution during summer baseflow period compared to during the spring flood period. Nutrients (e.g., K, P<sub>tot</sub>) and metals (e.g., Mn) presented an intermediate case with a half of yield occurring during the spring flood period and with negligible (<10%) role of autumn period in permafrost-free zone but comparable contribution of summer and autumn in the permafrost-affected part of the WSL.

Element export as a function of watershed area (Supplementary Figure S1) and season (Supplementary Figure S2) represent two main groups, namely (1) soluble highly mobile elements (Cl, SO<sub>4</sub>, Sr, Rb, Ba and As, Sb, Mo, U (permafrost-free zone only) and (2) DOC and low soluble elements, typically present in the form of organic- and organo-mineral colloids especially in the permafrost zone [30], such as micronutrients (Fe, Co, Ni, Cu, Zn, V), toxic metals (Cr, Cd, Pb), Nb, and trivalent (Al, Ga, Y, REE) and tetravalent (Ti, Hf, Th) hydrolysates. Generally, the share of spring flood contribution in overall open-water export of elements was 20 ± 10% higher in permafrost-affected rivers than in rivers of the permafrost-free zone.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Partial contribution of spring, summer and autumn 2016 to overall open-water period export of elements by WSL rivers, located in the permafrost-free (A) and permafrost-affected (B) regions.</figcaption></figure>

## 3.2. Riverine Element Export Fluxes Across the Latitudinal Profile (Permafrost and Climate Gradient)

Because both concentrations [33,34,37] and seasonal fluxes (Supplementary Figure S1 and Table S1) of elements were not strongly impacted by the river size, the watershed-averaged fluxes can be calculated as a function of the watershed latitude. Furthermore, we used ternary molar diagrams (Ca – Mg – (Na + K) and Cl – SO₄ – HCO₃) to reveal the role of season, river size and permafrost coverage of major elementary composition (Supplementary Figure S3), following traditional geochemical classifications [54]. Regardless of the season, water chemical composition from permafrost-free rivers was distinctly different from water from permafrost-affected rivers and strongly enriched in Ca and HCO₃⁻. As a result, the latitude (which corresponds with the permafrost zonation) can be considered as the main factor controlling major cation concentrations in the WSL rivers. The average (±2 s.d.) fluxes of dissolved (<0.45 µm) major element export from the WSL watersheds for each latitudinal belt in 2016 and 2015 are shown in Figure 3 and Supplementary Figure S4, respectively. For the sake of scientific rigor, we illustrate the yields separately for 2015 and 2016, but the two-year average values of element export during spring and summer are given in Table 1.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Yields (watershed-area normalized export fluxes) of DOC (A), DIC (B), Ca (C), Mg (D), K (E), Si (F), Ptot (G) and Fe (H) during May-October 2016 in 33 WSL rivers across the latitudinal gradient. The vertical uncertainties represent the 2 s.d. of several rivers belonging to the same latitudinal belt. Thick horizontal bars represent the 2° latitudinal belts. The year is represented by open-water period (May to October), neglecting winter (November to April), when small rivers freeze solid in the north.</figcaption></figure>
<table><caption>**Table 1.** Mean (±SD) export fluxes of major and trace elements by WSL rivers in the permafrost-free and the permafrost zone.</caption><thead><tr><th rowspan="2">Element<br>kg/km²/y</th><th colspan="2">Biennial Flux, 2015 + 2016<br>(May–August)</th><th colspan="2">Annual Flux, 2016 (May–October)</th></tr><tr><th>Permafrost-Free</th><th>Permafrost</th><th>Permafrost-Free</th><th>Permafrost</th></tr></thead><tbody><tr><td>DOC</td><td>1192 ± 310</td><td>6340 ± 1858</td><td>1165 ± 202</td><td>3399 ± 1279</td></tr><tr><td colspan="5"><em>Major anions</em></td></tr><tr><td>DIC</td><td>2122 ± 389</td><td>1150 ± 316</td><td>2278 ± 403</td><td>839 ± 280</td></tr><tr><td>Cl</td><td>132 ± 50</td><td>2092 ± 1052</td><td>160 ± 84</td><td>1150 ± 692</td></tr><tr><td>SO₄</td><td>409 ± 201</td><td>308 ± 112</td><td>486 ± 319</td><td>234 ± 116</td></tr><tr><td colspan="5"><em>Major cations</em></td></tr><tr><td>Na</td><td>452 ± 115</td><td>1784 ± 769</td><td>528 ± 163</td><td>1116 ± 592</td></tr><tr><td>Mg</td><td>554 ± 103</td><td>615 ± 174</td><td>602 ± 107</td><td>414 ± 125</td></tr><tr><td>Ca</td><td>3045 ± 545</td><td>1585 ± 510</td><td>3125 ± 557</td><td>1004 ± 368</td></tr><tr><td>K</td><td>136 ± 21</td><td>175 ± 56</td><td>130 ± 20</td><td>108 ± 36</td></tr><tr><td colspan="5"><em>Macro-nutrients</em></td></tr><tr><td>P<sub>tot</sub></td><td>6.41 ± 2.4</td><td>22.4 ± 8.2</td><td>3.66 ± 1</td><td>5.71 ± 2.3</td></tr><tr><td>P-PO₄</td><td>1.99 ± 0.67</td><td>10.7 ± 3.8</td><td>-</td><td>-</td></tr><tr><td>N-NO₃⁻</td><td>7.4 ± 1.8</td><td>8.12 ± 3.4</td><td>-</td><td>-</td></tr><tr><td>N-NH₄⁺</td><td>3.57 ± 0.8</td><td>14.6 ± 6.1</td><td>-</td><td>-</td></tr><tr><td>Si</td><td>409 ± 73</td><td>1528 ± 415</td><td>426 ± 112</td><td>1076 ± 430</td></tr><tr><td colspan="5"><em>Micro-nutrients</em></td></tr><tr><td>B</td><td>2.26 ± 0.43</td><td>5.63 ± 2.1</td><td>2.23 ± 0.39</td><td>3.82 ± 1.8</td></tr><tr><td>Mn</td><td>3.46 ± 1.8</td><td>36.7 ± 13</td><td>1.4 ± 0.69</td><td>16.9 ± 11</td></tr><tr><td>Fe</td><td>32.9 ± 14.3</td><td>292.7 ± 93.8</td><td>17.3 ± 7.5</td><td>92.2 ± 42</td></tr><tr><td>Co</td><td>0.012 ± 0.005</td><td>0.15 ± 0.06</td><td>0.009 ± 0.002</td><td>0.066 ± 0.04</td></tr><tr><td>Ni</td><td>0.126 ± 0.03</td><td>0.599 ± 0.22</td><td>0.118 ± 0.018</td><td>0.322 ± 0.15</td></tr><tr><td>Cu</td><td>0.159 ± 0.06</td><td>0.283 ± 0.12</td><td>0.128 ± 0.041</td><td>0.177 ± 0.087</td></tr><tr><td>Zn</td><td>1.31 ± 1</td><td>4.55 ± 3</td><td>0.203 ± 0.091</td><td>1.42 ± 0.89</td></tr><tr><td>Rb</td><td>0.092 ± 0.016</td><td>0.261 ± 0.098</td><td>0.078 ± 0.014</td><td>0.159 ± 0.059</td></tr><tr><td>Mo</td><td>0.023 ± 0.01</td><td>0.011 ± 0.004</td><td>0.029 ± 0.017</td><td>0.006 ± 0.002</td></tr><tr><td>Ba</td><td>3.32 ± 0.52</td><td>7.89 ± 3.4</td><td>2.87 ± 0.37</td><td>2.95 ± 1.4</td></tr><tr><td colspan="5"><em>Toxicants</em></td></tr><tr><td>As</td><td>0.098 ± 0.021</td><td>0.337 ± 0.11</td><td>0.093 ± 0.026</td><td>0.143 ± 0.047</td></tr><tr><td>Cd</td><td>0.0009 ± 0.0004</td><td>0.005 ± 0.002</td><td>0.0008 ± 0.0003</td><td>0.003 ± 0.001</td></tr><tr><td>Sb</td><td>0.008 ± 0.003</td><td>0.016 ± 0.005</td><td>0.008 ± 0.003</td><td>0.009 ± 0.004</td></tr><tr><td>Pb</td><td>0.023 ± 0.009</td><td>0.065 ± 0.031</td><td>0.005 ± 0.002</td><td>0.024 ± 0.012</td></tr><tr><td colspan="5"><em>Geochemical traces</em></td></tr><tr><td>Li</td><td>0.193 ± 0.059</td><td>0.57 ± 0.17</td><td>0.272 ± 0.094</td><td>0.455 ± 0.18</td></tr><tr><td>Al</td><td>5.37 ± 2.8</td><td>41.5 ± 18</td><td>5.86 ± 2.6</td><td>20.4 ± 11</td></tr><tr><td>Ti</td><td>0.223 ± 0.075</td><td>1.58 ± 0.56</td><td>0.512 ± 0.1</td><td>1.4 ± 0.52</td></tr><tr><td>Cr</td><td>0.041 ± 0.01</td><td>0.379 ± 0.12</td><td>0.024 ± 0.009</td><td>0.204 ± 0.08</td></tr><tr><td>Ga</td><td>0.002 ± 0.001</td><td>0.007 ± 0.002</td><td>0.0009 ± 0.0004</td><td>0.004 ± 0.002</td></tr><tr><td>Sr</td><td>15.7 ± 3</td><td>17.2 ± 8.8</td><td>16.6 ± 3.4</td><td>9.12 ± 4.8</td></tr><tr><td>Y</td><td>0.015 ± 0.006</td><td>0.11 ± 0.05</td><td>0.013 ± 0.003</td><td>0.056 ± 0.03</td></tr><tr><td>La</td><td>0.036 ± 0.029</td><td>0.085 ± 0.035</td><td>0.012 ± 0.003</td><td>0.05 ± 0.026</td></tr><tr><td>Ce</td><td>0.059 ± 0.048</td><td>0.177 ± 0.082</td><td>0.02 ± 0.006</td><td>0.091 ± 0.051</td></tr><tr><td>Pr</td><td>0.003 ± 0.001</td><td>0.023 ± 0.01</td><td>0.003 ± 0.0008</td><td>0.012 ± 0.007</td></tr><tr><td>Nd</td><td>0.013 ± 0.006</td><td>0.096 ± 0.044</td><td>0.012 ± 0.003</td><td>0.049 ± 0.028</td></tr><tr><td>Sm</td><td>0.003 ± 0.001</td><td>0.022 ± 0.01</td><td>0.003 ± 0.0008</td><td>0.012 ± 0.006</td></tr><tr><td>Eu</td><td>0.0009 ± 0.0003</td><td>0.006 ± 0.002</td><td>0.001 ± 0.0002</td><td>0.003 ± 0.002</td></tr><tr><td>Tb</td><td>0.0006 ± 0.0005</td><td>0.012 ± 0.006</td><td>0.003 ± 0.0007</td><td>0.012 ± 0.007</td></tr><tr><td>Gd</td><td>0.003 ± 0.001</td><td>0.014 ± 0.006</td><td>0.0007 ± 0.0002</td><td>0.002 ± 0.0009</td></tr><tr><td>Dy</td><td>0.003 ± 0.001</td><td>0.02 ± 0.009</td><td>0.003 ± 0.0007</td><td>0.01 ± 0.006</td></tr><tr><td>Ho</td><td>0.0005 ± 0.0002</td><td>0.004 ± 0.002</td><td>0.0005 ± 0.0001</td><td>0.002 ± 0.001</td></tr><tr><td>Er</td><td>0.002 ± 0.0007</td><td>0.012 ± 0.005</td><td>0.001 ± 0.0004</td><td>0.006 ± 0.003</td></tr><tr><td>Tm</td><td>0.0002 ± 0.0001</td><td>0.002 ± 0.0007</td><td>0.0002 ± 0.00005</td><td>0.0009 ± 0.0005</td></tr><tr><td>Yb</td><td>0.001 ± 0.0006</td><td>0.011 ± 0.005</td><td>0.001 ± 0.0003</td><td>0.006 ± 0.003</td></tr><tr><td>Lu</td><td>0.0002 ± 0.0001</td><td>0.002 ± 0.0007</td><td>0.0002 ± 0.00005</td><td>0.0009 ± 0.0005</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th rowspan="2">Element<br>kg/km²/y</th><th colspan="2">Biennial Flux, 2015 + 2016<br>(May—August)</th><th colspan="2">Annual Flux, 2016 (May—October)</th></tr><tr><th>Permafrost-Free</th><th>Permafrost</th><th>Permafrost-Free</th><th>Permafrost</th></tr></thead><tbody><tr><td>Zr</td><td>0.901 ± 0.33</td><td>1.47 ± 0.31</td><td>1.27 ± 0.46</td><td>1.57 ± 0.6</td></tr><tr><td>Cs</td><td>0.0003 ± 0.0001</td><td>0.0021 ± 0.0009</td><td>0.0003 ± 0.0001</td><td>0.001 ± 0.0005</td></tr><tr><td>Nb</td><td>0.0007 ± 0.0002</td><td>0.004 ± 0.002</td><td>0.0009 ± 0.0003</td><td>0.002 ± 0.001</td></tr><tr><td>Hf</td><td>0.0007 ± 0.0004</td><td>0.008 ± 0.003</td><td>0.001 ± 0.0004</td><td>0.005 ± 0.003</td></tr><tr><td>Th</td><td>0.002 ± 0.0009</td><td>0.009 ± 0.004</td><td>0.002 ± 0.0005</td><td>0.005 ± 0.003</td></tr><tr><td>U</td><td>0.021 ± 0.01</td><td>0.005 ± 0.002</td><td>0.03 ± 0.019</td><td>0.003 ± 0.001</td></tr></tbody></table>

Although the maximal uncertainties on discharge do not exceed 30% (see Section 2.2), the regional element yield assessments are subjected to high variation among rivers (Figure S1) with s.d. often at 50% (see Figure 3, Supplementary Figure S4). Among possible causes of these sizable uncertainties are non-linear and hysteretic relationships between concentrations and discharge which could not be resolved due to low frequency of sampling. Another reason of this variability could be highly dynamic behavior of element concentration during summer and autumn baseflow and spring flood, reflecting source limitation, chemostatis, or transport limitation (i.e., refs. [55–57]).

The DOC yield was minimal in rivers south to 59°N and exhibited a maximum at 61–65° N and 63–65° N in 2015 and 2016, respectively, with overall magnitude of variation by a factor of 4 (Figure S4A and Figure 3A, respectively). The yields of DIC and Ca decreased more than five-fold from south to north and achieved minimal values in the 63–65° N belt, for both years of observation (Figure 3B,C and Supplementary Figure S4B,C). Magnesium showed a weak minimum of yield in the 63–65° N which was however pronounced only in 2016; the overall variations were less than a factor of 2 (Figure 3D and Supplementary Figure S4D). Potassium yield remained fairly constant across the latitudinal profile with overall variations less than a factor of 1.5 to 2.0 between various latitudinal belts (Figure 3E and Supplementary Figure S4E). Silicon showed a three-fold increase in export fluxes from the south to the north, quite similar for two years of monitoring (Figure 3F and Supplementary Figure S4F).

The latitudinal pattern of other major and trace element export fluxes followed these three main types of behavior described above: (1) minimal in the south and a northward increase by a factor of 3 to 5 (Si, P-PO₄; N-NH₄, Ni, REEs, Zr, Nb, Hf, and Th), with a local maximum at 63–65° N (DOC, Al, Ti, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Rb, Cd, Cs, and Pb); (2) northward decrease by a factor of 3 to 5 (DIC, Ca, SO₄, Sr, Mo, W, and U), and (3) overall independence of yield on latitude, or the latitudinal variations were less than a factor of 2 (Cl, N-NO₃, Li, B, Mg, Ptot, K, As, Zr, Sb, and Ba). Note that some elements (Mg, Zn, As, Rb) may belong to one or another group depending on the year of observation.

The yields of all 51 dissolved elements are listed in Table 1, which presents the export fluxes of two-year average (May to August) and open-water period (May to October 2016) in permafrost-free and permafrost-affected zones. The ratio of element flux in the non-permafrost zone to that in the permafrost zone ($R_{absent/permafrost}$, Figure 4) demonstrates a distinct order of elements whose export occurs preferentially in southern or northern part of the WSL. There are three main groups of elements - those exhibiting the highest (a factor of 2 to 10) yield in the permafrost-free zone (DIC, SO₄, Ca, Mo and U) and those showing a maximum in the permafrost-affected zone (As, Ptot, Li, B, Rb, Na, Si with $0.5 \le R_{absent/permafrost} \le 0.25$, and Ni, Th, Cd, DOC, Cs, Nb, Ti, Al, Fe, Co, Cr, Mn, Hf, Co, and REEs with $R_{absent/permafrost} < 0.25$). The elements of an intermediate group showed comparable (±30%) yields in permafrost-free and permafrost-affected parts of the WSL (Sr, Mg, K, Ba, Sb, Zr, Cu).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** A histogram of the ratio of element average yield in the permafrost-free rivers to that in the permafrost-affected rivers for 6 open-water months (May to October) of 2016 (based on the data listed in Table 1).</figcaption></figure>

# 4. Discussion

## 4.1. *Effect of Seasons, Watershed Size, Latitude and Permafrost Coverage on Elementary Yields*

There are several distinct groups of elements defined according to the latitudinal patterns of their yields. These groups reflect the relative mobility of elements, consistent with general knowledge of river hydrochemistry in the high latitude regions [58]. The DOC, organically-bound metals (V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Cd and Pb) and many low soluble TE-geochemical traces (Al, REEs, Nb, Ti, Zr, Hf, and Th) exhibited similar 1st-type latitudinal pattern (minimum in the south and maximum in the permafrost-affected zone). The overall open-water period transport of these elements was dominated by spring flood, as it was also observed in other boreal and subarctic rivers (i.e., the Severnaya Dvina River, [59]). The nutrients (K, Si, P-PO₄, Ptot, and N-NH₄⁺) and TE—analogous of macronutrients (Rb)—exhibited 50 ± 10% share of annual export during spring flood and a northward increase. Finally, the yield of geochemically mobile elements (DIC, Cl⁻, SO₄, Li, B, Mg, Ca, Sr, Mo, and U) decreased northward, with less than 40–50% of overall contribution provided by spring flood period.

The division of elements into these groups depends on (1) the share of different seasons in overall 6-month open water period export, (2) the shape of the latitudinal pattern (Figure 3 and Figure S4), and (3) the difference in total 2-year averaged riverine yields between permafrost-free and permafrost-affected zones (Table 1, Figure 4). Such a distinction is consistent with two main factors controlling the element export (transport) in the WSL rivers. First, this is source-limitation, when the input of elements from soils to the river is controlled by river connectivity to deep and shallow groundwater reservoirs [27]. Due to the presence of carbonate mineral concretions in clay-silt bedrocks, especially in the south of the WSL [33], the groundwater is enriched in soluble elements such as alkalis and alkaline-earth metals, oxyanions and U(VI). The second is transport-limitation, when the export of an element is controlled by the availability of its carrier such as organic and organo-ferric, organo-aluminum colloids [30,34,45]. This factor reflects the superposition of surface source (topsoil, vegetation) providing organic colloids, and deep soil (mineral horizons), together with Fe(II)-rich groundwaters, providing soluble TE. In addition to these two main groups of elements, the mineral nutrients are not limited by transport and exhibit quite complex pattern which reflects their deep (groundwater, bedrock lithology) and surface (plant litter, atmospheric deposition) sources. These elements are strongly impacted by their biotic uptake in the river channel or seasonal release from decaying aquatic macrophytes, plankton and periphyton [37,60]. Various internal factors, operating in soils and riparian zone of the river, are capable of modifying the export of nutrients. For example,
low export fluxes of both total P (phosphate and organic P) and P-PO₄ in southern (56–60° N) latitudinal belt compared to northern, permafrost-affected rivers (see Figure 3G, Figure S4G) may be due P retention via adsorption onto and coprecipitation with Al, Fe hydroxides in the deeper part of soil profile [61,62] as well as P removal in the form of Ca phosphate minerals [63] occurring in the unfrozen mineral soils, exposed to surface fluids. Furthermore, P uptake by abundant terrestrial and aquatic vegetation is most pronounced in the south. In the north, the mineral soils are essentially frozen, and the requirements of aquatic and terrestrial biota for this nutrient are much lower [64].

The effect of watershed size on DOC, DIC, cations, and Si export fluxes was rather minor. This is at odds with strong river size control on element yields as it is known in temperate and mountainous regions [65,66]. In particular, small catchments of wetlands exhibit generally lower runoff than the medium and large rivers, and the runoff is one of the major controlling factors of chemical weathering and element export [46,67–70]. Season also played a secondary role in determining element yield pattern. However, sizable correlations of element fluxes with S<sub>watershed</sub> in summer and autumn, observed solely in the permafrost-free zone (Supplementary Table S1) can be explained by more pronounced impact of deep underground waters in large rivers compared to small ones. These waters typically contain a high concentration of soluble, labile elements (e.g., Cl, SO₄, Na, Mo, As) [33,34] but also Fe(II), those oxidation in the riparian and hyporheic zones create large amount of organo-ferric colloids [45] capable to provide enhanced concentrations of typically insoluble low mobile elements such as trivalent and tetravalent hydrolysates.

The latitude was revealed to be the primary governing parameter of elementary yields and clearly marked the difference between permafrost-free, discontinuous and continuous permafrost regions. The most northern regions of the WSL exhibited rather high yields of DOC, DIC, Si and cations. It is possible that in continuous permafrost zone of frozen peat bogs, the underlining mineral layer is protected by the permafrost. As a result, the active (seasonally unfrozen) layer is located within the organic layer which comprises live vegetation, plant litter, and upper peat layer. This organic matrix is extremely reactive, and capable of releasing sizable amount of DOC, P, Si, Ca and nutrients over very short periods of time during contact with surface waters [71,72]. Water temperature (between 4 and 25 °C) has only minor effect on C and element release from both thawed and frozen peat [72]. Therefore, even short-term contact of water with surface peat and vegetation is capable mobilizing sizable amount of DOC and nutrients, despite low temperature in the northern regions. In this regard, the elementary export by the WSL rivers is strongly controlled by dynamics of peat formation/decay across the territory. The particularity of the WSL is that, currently, this region is recovering from the last glaciation. As a result, the ecosystems are highly non-stationary: the peat actively accumulates in the south [73], while in the north, the frozen peat is subjected to thawing and degradation [74–76]. The uptake of elements from groundwater, river and forest tree litter by growing peat in southern mires counteracts with DOC and element release from thawing/degrading peat in the northern palsas. The elements affected by these processes are those that exhibit the highest concentration in peat relative to undelaying mineral (silt, sand) horizons. According to a previous assessment of elementary peat composition in the WSL [77], the peat is sizably enriched in C, (V, Cr), trivalent (TE³⁺) and tetravalent (TE⁴⁺) trace metals (Al, Y, REEs, Ti, Zr, Hf, Th) and U, Zn, Pb and depleted in highly mobile alkalis and alkaline-earths metals, As, and Mo. Therefore, enhanced riverine yield of low-soluble TE³⁺, TE⁴⁺, U, and some divalent metals in the north relative to the south is possibly due to these elements being tightly linked to peatland evolutionary pattern across the WSL. Similarly, depletion of peat relative to underlying mineral horizons in soluble, highly mobile elements is consistent with enhanced export of these elements by southern rivers, where the peat accumulation occurs.

## 4.2. *Comparison of Major Cation, DIC, Si and DOC Export Fluxes in the WSL with Other Boreal Regions*

The average total dissolved cation flux (TDS_c = Na + K + Ca + Mg) over May-October 2016 ranged from 4.40 ± 0.55 t km⁻² y⁻¹ in the permafrost-free zone to 2.64 ± 0.59 t km⁻² y⁻¹ in the permafrost-affected zone, which is lower than the fluxes of Central Siberian rivers of the same latitude,
draining basaltic rocks (5 to 8 t km⁻² y⁻¹, [78]), large Siberian rivers such as Yenisey and Lena (6.2 and 6.8 t km⁻² y⁻¹, respectively, [79]), the Ob River in its middle course (6.0 t km⁻² y⁻¹, [80]), and the permafrost-free Eurasian Arctic rivers draining sedimentary rocks (Sev. Dvina, 9.5 t km⁻² y⁻¹, [59]; Pechora, 6.6 t km⁻² y⁻¹, [79]). The TDS_c yield of permafrost-affected WSL rivers is, however, comparable with previous estimations of that in the middle-size Siberian rivers (2.8, 2.5, and 2.3 t km⁻² y⁻¹ for Kolyma, Indigirka and Anabar, respectively [79]). The main reason for relatively low cationic fluxes of small WSL rivers compared to other large rivers of Northern Eurasia of similar runoff are i) low connectivity of WSL rivers with underlying bedrocks and groundwater, due to thick peat layer and permafrost, and ii) essentially weathered character of silicate rocks (sands, clays) in western Siberia. Note also that the transport of major cations in permafrost-free rivers is strongly pronounced during winter (e.g., 35–40% of total annual yield in the Severnaya Dvina River [59]), so it is possible that cationic fluxes of WSL rivers during May-October are somewhat underestimated relative to annual export. At the same time, the majority of small (<4000 km² watershed area) rivers in the northern part of the WSL freeze solid during the winter [33], so the winter flux may be non-negligible only for southern, permafrost-free rivers.

In contrast to major cations, no difference in the export fluxes between small WSL rivers and large and medium size Eurasian rivers was detected for the DIC. The DIC export of small WSL rivers in the permafrost-free zone (2.2 ± 0.4 t km⁻² y⁻¹) is in agreement with recent estimations of DIC export in the middle course of the Ob River (2.9 t km⁻² y⁻¹, [80]) and with the mean riverine DIC yield of the entire Eurasian Arctic basin (2.2 t km⁻² y⁻¹, [79]). However, DIC export by permafrost-affected WSL rivers is somewhat lower (1–2 t km⁻² y⁻¹) and comparable to medium size rivers of Central and Eastern Siberia (0.6–2.2 t km⁻² y⁻¹, [79]). A tentative explanation of elevated DIC (but not TDS_c) flux in the WSL rivers is CO₂ and HCO₃⁻ generation by microbial (and photolytic) processing of peat soil organic carbon (both DOC and POC), which is the main cause of very high CO₂ emission from WSL inland waters [81,82]. The light isotopic composition of DIC in the WSL rivers (−25 ≤ δ¹³C_DIC < −10 ‰, [33]) is consistent with this hypothesis.

The riverine Si fluxes in the southern part (<61° N) of the WSL (0.5–1.0 t km⁻² y⁻¹) are comparable to those of the Ob River (0.62 t km⁻² y⁻¹, [80] and small rivers of the northern Sweden (0.9 t km⁻² y⁻¹, [83]). In contrast, the fluxes of the permafrost-affected rivers (1.0–1.5 t km⁻² y⁻¹, Figure 3F and Figure S4F and Table 1) contradict the expected trend of decreasing flux with the increasing latitude and decreasing temperature, given that the chemical weathering of silicate rocks is much slower in colder climates [67–70]. In fact, we noted that the northern fluxes of WSL are comparable with those of the temperate rivers such as Mississippi and Yangtze [84]. Moreover, if the silicate rock weathering significantly controls element delivery from the soil to the river, such a northward increase would occur for cations (Ca, Mg, Na) as it is known for typical silicate terrains of the boreal zone [85], but this is not observed in the WSL territory. As such, we hypothesize that two- to three-fold increase in Si yield of northern rivers relatively to southern rivers is linked to combination of (1) enhanced mobilization from plant litter via supapermafrost flow in the permafrost zone, (2) limited silicate secondary mineral formation in shallow, essentially frozen northern soils compared to southern soils, and (3) strong uptake of Si by both terrestrial plants and aquatic macrophytes and periphyton in the southern rivers.

The DOC fluxes in the permafrost-affected regions of the WSL territory were quite high (2 to 6 t km⁻² y⁻¹) compared to the middle and lower reaches of the Ob River (0.64 and 1.2 t km⁻² y⁻¹, respectively, [79,80]). At the same time, these fluxes are comparable with those in large boreal river draining permafrost-free wetlands (the Severnaya Dvina River, 4.2 ± 0.8 t Cₐrg km⁻² y⁻¹, [59]) and are the highest among all known rivers flowing to the Arctic Ocean. Indeed, the Ob, Yukon, Lena, Yenisey, and Mackenzie rivers exhibit a DOC yield from 0.5 to 2.5 t km⁻² y⁻¹ ([79,86]). We suggest that the main factor responsible for such high DOC yield in small WSL rivers is high proportion of peatlands on their watershed (i.e., typically from 40 to 60%, according to the GIS data [34]). The peatland-draining Taz (watershed = 150,000 km²), Pur (112,000 km²), and Nadym (64,000 km²) rivers, located entirely
in the discontinuous permafrost zone, also have a DOC yield of 1.9, 2.1, and 4.4 t C<sub>org</sub> km<sup>-2</sup> y<sup>-1</sup>, respectively ([79] and calculations based on data from the RHS).

The northward increase in DOC flux possibly reflects strong leaching of OM from the plant litter and organic-rich topsoil (Histosol). In the north, the adsorption of DOM on underlying mineral horizons is minimal because these horizons are frozen. As a result, the riverine DOC export in the permafrost zone of the WSL is controlled by water travel time through the peat layer and underlying mineral horizons and the water residence time necessary for DOC leaching from upper vegetation layer (moss, lichen, litter). However, quantitative modeling of DOM and element reactive transport in the WSL peatlands, on the scale of a small watershed or large river, is beyond the scope of this study.

### 4.3. Possible Evolution of Western Siberia Rivers Elementary Yields Under Climate Change Scenario

The space for time approach employed in the present work provides some future projections for riverine element behavior. However, it exhibits a number of shortcomings whose analysis goes beyond the scope of this work. In particular, this approach does not address the time scale, necessary for the northern ecosystem to reach the new “more southern” state; it ignores possible shift in the structure of vegetation and soil microbial community that can indirectly impact the carbon and nutrient, in terms of landscape storage and removal via rivers, and it does not include information on the altered seasonality such as extended hydrologic seasons, earlier snowmelt, higher precipitation that will likely occur as a result of climate change. Nevertheless, as a first order empirical assessment, the following predictions can be made. The first consequence of climate warming in western Siberia is thawing of frozen peat and underlying mineral horizons. The thickness of the active layer (ALT) is projected to increase more than 30% during this century across the tundra area in the Northern Hemisphere [87–89]. In the WSL, this increase will be most dramatic in the north, where the peat deposits are thinner than those in the discontinuous permafrost zone [48,90–92]. The main consequences of the ALT increase may be the involvement of mineral (clay) horizons into water infiltration within the soil profile [90]. As a result, the DOC originated from the leaching of the upper peat layers and plant litter degradation will be retained on mineral surfaces via adsorption onto and incorporation into clay interlayers [93–97]. For inorganic solutes, the effect of ALT increase will be lower than that of DOC, given much lower affinity of HCO<sub>3</sub><sup>-</sup> cations and Si to clay surfaces and the lack of unweathered (primary) silicate rocks underneath the peat soil column. However, if the thawing will open new water paths between deep groundwater reservoirs and the river, this may increase the riverine export of major cations and DIC [9,33].

The second consequence of climate change in western Siberia is a shift of the permafrost zone boundaries further north [98–100]. Within the substituting space for time scenario, such a permafrost boundary change can be considered equivalent to the northward shift of the river latitudes as shown in Figure 3. Based on latitudinal pattern of major and TE (as illustrated in Figure 3 and Figure S4), the following groups of the river water components likely to change their riverine yield over open-water period in case of anticipated shift in the permafrost zones: (1) Elements those yields likely to increase by a factor of 2 to 3: DIC, Ca, SO<sub>4</sub>, Sr, Ba, Mo, U; (2) Elements likely to decrease the yields by a factor of 2 to 5: DOC, Fe, Si, Ptot, P-PO<sub>4</sub>, N-NH<sub>4</sub>, divalent heavy metals, trivalent and tetravalent hydrolysates, and (3) Elements weakly affected by the permafrost boundary change and the ALT increase (change of yield by a factor of 1.5 to 2.0): B, Li, K, Rb, Cs, N-NO<sub>3</sub>, Mg, Zn, As, Sb, Rb, Cs).

The impact of climate warming on riverine fluxes in the WSL is not restricted to the shift of permafrost zones and the change of water flow path (deep versus surface). It has to be placed in the context of changing precipitation, plant biomass productivity and modification of the seasonality [101]. Complex evaluation of these factors goes beyond the scope of this study and it requires ecosystem-level regional modeling. Overall, the WSL is likely to have increased lateral export of DIC but decreased export of DOC. However, the changes of both fluxes (+2 and -3 to -4 t C km<sup>-2</sup> y<sup>-1</sup> for DIC and DOC, respectively) are dwarfed compared to possible magnitude of C emission from WSL inland waters (rivers and lakes) to the atmosphere: 10–20 t km<sup>-2</sup> y<sup>-1</sup> in permafrost-free and isolated zone and 20–40 t C km<sup>-2</sup> y<sup>-1</sup> in discontinuous and continuous permafrost zone [102].
**5. Conclusions**

Based on a two-year seasonal sampling of 32 western Siberian rivers of various size (from 10 to 10⁵ km² in watershed area) draining through a sizable permafrost gradient, we measured riverine export fluxes of dissolved (<0.45 µm) C, nutrients, major and TE over a six-month open-water period (May to October). The export fluxes of DOC, DIC, major cations, macro- and micro-nutrients, toxicants, and geochemical tracers were weakly dependent on the size of the river. The primary parameter of export fluxes control was latitude, which marked the position of the permafrost zones.

There are several distinct groups of elements defined according to the latitudinal patterns of their yields. These groups reflect the relative mobility of elements, consistent with general knowledge of river hydrochemistry in high latitude regions. The DOC, organically-bound metals (V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Cd and Pb) and many low soluble TE (Al, Ti, Zr, Hf, Th and REEs) exhibited similar latitudinal pattern with a minimum in the south and a maximum in the permafrost zone). A northward increase in Si export flux may be due to a decrease in Si uptake by plants in the north and strong Si retaining by mire and forest vegetation in the south.

An increase in DOC export fluxes from the south to the north (by a factor of 3 to 4 depending on the years of observation) could be due to leaching of OM from the plant litter and organic-rich topsoil (Histosol). The removal of DOC by adsorption on mineral horizons was hypothesized to be very low in the north. As a result, the riverine DOC export in the permafrost zone of the WSL may be strongly controlled by the water residence time necessary for DOC leaching from upper vegetation layer (moss, lichen, litter). This calls a need for quantitative modeling of DOM and element reactive transport in WSL peatlands, both at the scale of small watershed and large rivers. Furthermore, the peculiarity of western Siberia is that the elementary export by WSL rivers is strongly controlled by dynamics of peat formation and decay across the territory. Because the WSL is currently recovering from the last glaciation, this territory is dominated by non-stationary ecosystems with strong latitudinal contrast: the fresh peat is accumulating in the south whereas the old frozen peat is thawing and degrading in the north. As a result, uptake of elements from groundwater, river and forest tree litter by growing peat in southern bog competes with DOC and element release from thawing/degrading peat in the northern palsas. Note that, unlike many other permafrost-affected regions in the world whose CO₂ uptake rate during weathering is likely to increase under climate warming, the WSL may increase its riverine export of DIC but decrease the export of DOC.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2073-4441/12/6/1817/s1, Figure S1: Correlation between elementary fluxes and watershed area for permafrost-free and permafrost-affected regions; Figure S2: Partial contribution of spring, summer and autumn 2016 to overall open-water period export of anions and trace elements by WSL rivers; Figure S3: Ternary molar diagrams of major cations and anions in the WSL rivers. Figure S4: Yields of DOC (A), DIC (B), Ca (C), Mg (D), K (E), Si (F), Ptot (G) and Fe (H) during May-August 2015 in 33 WSL rivers across the latitudinal gradient. Table S1: Spearman correlation coefficients (p < 0.05) between element export flux (yield) and watershed area.

**Author Contributions:** O.S.P. designed the study and wrote the paper; R.M.M., I.V.K., and S.V.L. performed sampling, analysis of cations and anions, and their interpretation; S.N.V. and S.N.K. were responsible for the choice of sampling objects and statistical treatment; S.V.L. provided the background information on soil, peat, and permafrost active layer; L.S.S. was in charge of DOC, DIC, and anion measurements and their interpretation; L.G.K. provided GIS-based interpretation, mapping, and identification of river watersheds; S.G.K. performed all primary hydrological data collection, and their analysis and interpretation. All 10 authors participated in field expeditions. All authors have read and agreed to the published version of the manuscript.

**Funding:** Russian Science Foundation: No 18-17-00237 and 18-77-10045. Russian Fund for Fundaental Research: 19-55-15002, 20-05-00729_a.

**Acknowledgments:** We acknowledge main support from RSF grant No 18-17-00237 and RFBR grants No 19-55-15002, 20-05-00729_a, and RSF grant No 18-77-10045 for field work.

**Conflicts of Interest:** The authors declare no conflict of interest.
## References

1. Turetsky, M.R.; Abbott, B.W.; Jones, M.C.; Anthony, K.W.; Olefeldt, D.; Schuur, E.A.G.; Grosse, G.; Kuhry, P.; Hugelius, G.; Koven, C.; et al. Carbon release through abrupt permafrost thaw. *Nat. Geosci.* **2020**, *13*, 138–143.
[CrossRef]

2. Vonk, J.E.; Tank, S.E.; Walvoord, M.A. Integrating hydrology and biogeochemistry across frozen landscapes. *Nat. Commun.* **2019**, *10*, 1–4. [CrossRef] [PubMed]

3. Smith, L.C.; Sheng, Y.; MacDonald, G.M. A first Pan-Arctic assessment of the influence of glaciation, permafrost, topography and peatlands on Northern hemisphere lake distribution. *Permafr. Periglac. Proc.* **2007**, *18*, 201–208. [CrossRef]

4. McEnroe, N.A.; Roulet, N.T.; Moore, T.R.; Garneau, M. Do pool surface area and depth control CO₂ and CH₄ fluxes from an ombrotrophic raised bog, James Bay, Canada? *J. Geophys. Res. Biogeosci.* **2009**, *114*. [CrossRef]

5. Olefeldt, D.; Roulet, N.T. Effects of permafrost and permafrost thaw on dissolved organic carbon in a subarctic peatland complex. *J. Geophys. Res. Biogeosci.* **2012**, *117*. [CrossRef]

6. Olefeldt, D.; Roulet, N.; Giesler, R.; Persson, A. Total waterborne carbon export and DOC composition from ten nested subarctic peatland catchments—Importance of peatland cover, groundwater influence, and inter-annual variability of precipitation patterns. *Hydrol. Process.* **2013**, *27*, 2280–2294. [CrossRef]

7. Roehm, C.L.; Giesler, R.; Karlsson, J. Bioavailability of terrestrial organic carbon to lake bacteria: The case of a degrading subarctic permafrost mire complex. *J. Geophys. Res. Biogeosci.* **2009**, *114*. [CrossRef]

8. Amesbury, M.J.; Gallego-Sala, A.; Loisel, J. Peatlands as prolific carbon sinks. *Nat. Geosci.* **2019**, *12*, 880–881. [CrossRef]

9. Frey, K.E.; McClelland, J.W. Impacts of permafrost degradation on arctic river biogeochemistry. *Hydrol. Process.* **2009**, *23*, 169–182. [CrossRef]

10. Bartsch, A.; Leibman, M.; Strozzi, T.; Khomutov, A.; Widhalm, B.; Babkina, E.; Mullanurov, D.; Ermokhina, K.; Kroisleitner, C.; Bergstedt, H. Seasonal Progression of Ground Displacement Identified with Satellite Radar Interferometry and the Impact of Unusually Warm Conditions on Permafrost at the Yamal Peninsula in 2016. *Remote Sens.* **2019**, *11*, 1865. [CrossRef]

11. Beilman, D.W.; MacDonald, G.M.; Smith, L.C.; Reimer, P.J. Carbon accumulation in peatlands of West Siberia over the last 2000 years. *Glob. Biogeochem. Cycles* **2009**, *23*. [CrossRef]

12. Liljedahl, A.K.; Boike, J.; Daanen, R.P.; Fedorov, A.N.; Frost, G.V.; Grosse, G.; Hinzman, L.D.; Iijma, Y.; Jorgenson, J.C.; Matveyeva, N.; et al. Pan-Arctic ice-wedge degradation in warming permafrost and its influence on tundra hydrology. *Nat. Geosci.* **2016**, *9*, 312–318. [CrossRef]

13. Colombo, M.; Brown, K.A.; De Vera, J.; Bergquist, B.A.; Orians, K.J. Trace metal geochemistry of remote rivers in the Canadian Arctic Archipelago. *Chem. Geol.* **2019**, *525*, 479–491. [CrossRef]

14. Drake, T.W.; Tank, S.E.; Zhulidov, A.V.; Holmes, R.M.; Gurtovaya, T.; Spencer, R.G.M. Increasing alkalinity export from large russian Arctic rivers. *Environ. Sci. Technol.* **2018**, *52*, 8302–8308. [CrossRef]

15. Fabre, C.; Sauvage, S.; Tananaev, N.; Noël, G.E.; Teisserenc, R.; Probst, J.L.; Pérez, J.M.S. Assessment of sediment and organic carbon exports into the Arctic ocean: The case of the Yenisei River basin. *Water Res.* **2019**, *158*, 118–135. [CrossRef]

16. Griffin, C.G.; McClelland, J.W.; Frey, K.E.; Fiske, G.; Holmes, R.M. Quantifying CDOM and DOC in major Arctic rivers during ice-free conditions using Landsat TM and ETM+ data. *Remote Sens. Environ.* **2018**, *209*, 395–409. [CrossRef]

17. Guo, L.; Ping, C.-L.; Macdonald, R.W. Mobilization pathways of organic carbon from permafrost to arctic rivers in a changing climate. *Geophys. Res. Lett.* **2007**, *34*. [CrossRef]

18. Hindshaw, R.S.; Teisserenc, R.; Le Dantec, T.; Tananaev, N. Seasonal change of geochemical sources and processes in the Yenisei River: A Sr, Mg and Li isotope study. *Geochim. Cosmochim. Acta* **2019**, *255*, 222–236. [CrossRef]

19. Li, M.; Peng, C.; Zhou, X.; Yang, Y.; Guo, Y.; Shi, G.; Zhu, Q. Modeling global riverine DOC flux dynamics from 1951 to 2015. *J. Adv. Modeling Earth Syst.* **2019**, *11*, 514–530. [CrossRef]

20. Murphy, M.J.; Porcelli, D.; Pogge von Strandmann, P.A.E.; Hirst, C.A.; Kutscher, L.; Katchinoff, J.A.; Mörth, C.-M.; Maximov, T.; Andersson, P.S. Tracing silicate weathering processes in the permafrost-dominated Lena River watershed using lithium isotopes. *Geochim. Cosmochim. Acta* **2019**, *245*, 154–171. [CrossRef]
21. Holmes, R.M.; Coe, M.T.; Fiske, G.J.; Gurtovaya, T.; McClelland, J.W.; Shiklomanov, A.I.; Spencer, R.G.M.; Tank, S.E.; Zhulidov, A.V. Climate change impacts on the hydrology and biogeochemistry of Arctic Rivers. In *Climatic Changes and Global warming of Inland Waters; Impacts and Mitigation for Ecosystems and Societies*; Goldman, C.R., Kumagi, M., Robarts, R.D., Eds.; John Wiley and Sons: Hoboken, NJ, USA, 2013; pp. 1–26.

22. McClelland, J.W.; Déry, S.J.; Peterson, B.J.; Holmes, R.M.; Wood, E.F. A pan-arctic evaluation of changes in river discharge during the latter half of the 20th century. *Geophys. Res. Lett.* **2006**, *33*. [CrossRef]

23. Shiklomanov, A.I.; Lammers, R.B.; Vörösmarty, C.J. Widespread decline in hydrological monitoring threatens Pan-Arctic research. *Eos Trans. Am. Geophys. Union* **2002**, *83*, 13–20. [CrossRef]

24. Shogren, A.J.; Zarnetske, J.P.; Abbott, B.W.; Iannucci, F.; Frei, R.J.; Griffin, N.A.; Bowden, W.B. Revealing biogeochemical signatures of Arctic landscapes with river chemistry. *Sci. Rep.* **2019**, *9*, 1–11. [CrossRef] [PubMed]

25. Toohey, R.C.; Herman-Mercer, N.M.; Schuster, P.F.; Mutter, E.A.; Koch, J.C. Multidecadal increases in the Yukon River Basin of chemical fluxes as indicators of changing flowpaths, groundwater, and permafrost. *Geophys. Res. Lett.* **2016**, *43*, 12–120. [CrossRef]

26. Frey, K.E.; McClelland, J.W.; Holmes, R.M.; Smith, L.C. Impacts of climate warming and permafrost thaw on the riverine transport of nitrogen and phosphorus to the Kara Sea. *J. Geophys. Res.* **2007**, *112*. [CrossRef]

27. Frey, K.E.; Siegel, D.I.; Smith, L.C. Geochemistry of west Siberian streams and their potential response to permafrost degradation. *Water Resour. Res.* **2007**, *43*. [CrossRef]

28. Frey, K.E.; Smith, L.C. Amplified carbon release from vast West Siberian peatlands by 2100. *Geophys. Res. Lett.* **2005**, *32*. [CrossRef]

29. Krickov, I.; Lim, A.; Manasypov, R.M.; Loiko, S.V.; Shirokova, L.S.; Kirpotin, S.N.; Karlsson, J.; Pokrovsky, O.S. Riverine particulate C and N generated at the permafrost thaw front: Case study of western Siberian rivers across a 1700-km latitudinal transect. *Biogeosciences* **2018**, *15*, 6867–6884. [CrossRef]

30. Krickov, I.V.; Pokrovsky, O.S.; Manasypov, R.M.; Lim, A.G.; Shirokova, L.S.; Viers, J. Colloidal transport of carbon and metals by western Siberian rivers during different seasons across a permafrost gradient. *Geochim. Cosmochim. Acta* **2019**, *265*, 221–241. [CrossRef]

31. Krickov, I.V.; Lim, A.G.; Manasypov, R.M.; Loiko, S.V.; Vorobyev, S.N.; Shevchenko, V.P.; Dara, O.M.; Gordeev, V.V.; Pokrovsky, O.S. Major and trace elements in suspended matter of western Siberian rivers: First assessment across permafrost zones and landscape parameters of watersheds. *Geochim. Cosmochim. Acta* **2020**, *269*, 429–450. [CrossRef]

32. Manasypov, R.M.; Pokrovsky, O.S.; Kirpotin, S.N.; Shirokova, L.S. Thermokarst lake waters across the permafrost zones of western Siberia. *Cryosphere* **2014**, *8*, 1177–1193. [CrossRef]

33. Pokrovsky, O.S.; Manasypov, R.M.; Shirokova, L.S.; Loiko, S.V.; Krickov, I.V.; Kopysov, S.; Zemtzov, V.A.; Kulizhsky, S.P.; Vorobyev, S.N.; Kirpotin, S.N. Permafrost coverage, watershed area and season control of dissolved carbon and major elements in western Siberia rivers. *Biogeosciences* **2015**, *12*, 6301–6320. [CrossRef]

34. Pokrovsky, O.S.; Manasypov, R.M.; Loiko, S.; Krickov, I.A.; Kopysov, S.G.; Kolesnichenko, L.G.; Vorobyev, S.N.; Kirpotin, S.N. Trace element transport in western Siberia rivers across a permafrost gradient. *Biogeosci. Discuss.* **2016**, *13*, 1877–1900. [CrossRef]

35. Polishchuk, Y.M.; Bogdanov, A.N.; Muratov, I.N.; Polishchuk, V.Y.; Lim, A.; Manasypov, R.M.; Shirokova, L.S.; Pokrovsky, O.S. Minor contribution of small thaw ponds to the pools of carbon and methane in the inland waters of the permafrost-affected part of the Western Siberian Lowland. *Environ. Res. Lett.* **2018**, *13*, 045002. [CrossRef]

36. Raudina, T.V.; Loiko, S.V.; Lim, A.G.; Krickov, I.V.; Shirokova, L.S.; Istigichev, G.I.; Kuzmina, D.M.; Kulizhsky, S.P.; Vorobyev, S.N.; Pokrovsky, O.S. Dissolved organic carbon and major and trace elements in peat porewater of sporadic, discontinuous, and continuous permafrost zones of western Siberia. *Biogeosciences* **2017**, *14*, 3561–3584. [CrossRef]

37. Vorobyev, S.N.; Pokrovsky, O.S.; Serikova, S.; Manasypov, R.M.; Krickov, I.V.; Shirokova, L.S.; Lim, A.; Kolesnichenko, L.G.; Kirpotin, S.N.; Karlsson, J. Permafrost boundary shift in Western Siberia may not modify dissolved nutrient concentrations in rivers. *Water* **2017**, *9*, 985. [CrossRef]

38. Raudina, T.V.; Loiko, S.V.; Lim, A.; Manasypov, R.M.; Shirokova, L.S.; Istigechev, G.I.; Kuzmina, D.M.; Kulizhsky, S.P.; Vorobyev, S.N.; Pokrovsky, O.S. Permafrost thaw and climate warming may decrease the CO₂, carbon, and metal concentration in peat soil waters of the Western Siberia Lowland. *Sci. Total Environ.* **2018**, *634*, 1004–1023. [CrossRef]
39. Blois, J.L.; Williams, J.W.; Fitzpatrick, M.C.; Jackson, S.T.; Ferrier, S. Space can substitute for time in predicting climate-change effects on biodiversity. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 9374–9379. [CrossRef]

40. Nikitin, S.P.; Zemtsov, V.A. *The Variability of Hydrological Parameters of Western Siberia*; Nauka: Novosibirsk, Russia, 1986; 204p. (In Russian)

41. Huh, Y.; Edmond, J.M. The fluvial geochemistry of the rivers of Eastern Siberia: III. Tributaries of the Lena and Anabar draining the basement terrain of the Siberian Craton and the Trans-Baikal Highlands. *Geochim. Cosmochim. Acta* **1999**, *63*, 967–981. [CrossRef]

42. Spence, J.; Telmer, K. The role of sulfur in chemical weathering and atmospheric CO₂ fluxes: Evidence from major ions, δ¹³C_DIC, and δ³⁴S SO₄ in rivers of the Canadian Cordillera. *Geochim. Cosmochim. Acta* **2005**, *69*, 5441–5458. [CrossRef]

43. Qin, J.; Huh, Y.; Edmond, J.M.; Du, G.; Ran, J. Chemical and physical weathering in the Min Jiang, a headwater tributary of the Yangtze River. *Chem. Geol.* **2006**, *227*, 53–69. [CrossRef]

44. Bulygina, O.N.; Razuvaev, V.N.; Alexsandrova, T.M. Description of Daily Air Temperature and Precipitation Data set from Russian Meteorological Stations and from some Meteorological Stations over the Former USSR Territory (TTTR) [Electronic Resource]. Available online: http://meteo.ru/english/climate/descrip11.htm (accessed on 17 March 2020).

45. Pokrovsky, O.S.; Manasypov, R.M.; Loiko, S.V.; Shirokova, L.S. Organic and organo-mineral colloids in discontinuous permafrost zone. *Geochim. Cosmochim. Acta* **2016**, *188*, 1–20. [CrossRef]

46. Prokushkin, A.S.; Pokrovsky, O.S.; Shirokova, L.S.; Korets, M.A.; Viers, J.; Prokushkin, S.G.; Amon, R.M.W.; Guggenberger, G.; McDowell, W.H. Sources and the flux pattern of dissolved carbon in rivers of the Yenisey basin draining the Central Siberian Plateau. *Environ. Res. Lett.* **2011**, *6*, 045212. [CrossRef]

47. Rozhdestvensky, A.V.; Buzin, V.A.; Dobroumov, B.M.; Lobanova, A.G.; Lobanov, V.A.; Plitkin, G.A.; Tumanovskaya, S.M. SP 33-101-2003. Svod Pravil. Opredelenie Osnovny Khraschetnykh Gidrologicheskikh Kharakteristik [Set of Rules. Determination of Basic Design Hydrological Characteristics]. Odobren dlya Primeneniya v Kachestve Normativnogo Dokumenta Postanovleniem Gosroya Rossii No 218 ot 26.12.2003. Available online: https://files.stroyinf.ru/Data2/1/4294815/4294815038.htm (accessed on 18 March 2020).

48. Novikov, S.M.; Moskvin, Y.P.; Trofimov, S.A.; Usova, L.I.; Batuev, V.I.; Tumanovskaya, S.M.; Smirnova, V.P.; Markov, M.L.; Korotkevicth, A.E.; Potapova, T.M. *Hydrology of Bog Territories of the Permafrost Zone of Western Siberia*; BBM Publishing House: St. Petersburg, Russia, 2009; 535p. (In Russian)

49. McClelland, J.W.; Holmes, R.M.; Peterson, B.J.; Raymond, P.A.; Striegl, R.G.; Zhulidov, A.V.; Zimov, S.A.; Zimov, N.; Tank, S.E.; Spencer, R.G.M.; et al. Particulate organic carbon and nitrogen export from major Arctic rivers. *Glob. Biogeochem. Cycles* **2016**, *30*, 629–643. [CrossRef]

50. Bergström, S. The HBV Model. In *Structure and Applications*; Swedish Meteorological and Hydrological Institute (SMHI): Norrköping, Sweden, 1992; p. 35.

51. Bergström, S. The HBV model. In *Computer Models of Watershed Hydrology*; Singh, V.P., Ed.; Water Resources Publications: Highlands Ranch, CO, USA, 1995; Chapter 13; pp. 443–476.

52. Seibert, J.; Vis, M.J.P. Teaching hydrological modeling with a user-friendly catchment-runoff-model software package. *Hydrol. Earth Syst. Sci.* **2012**, *16*, 3315–3325. [CrossRef]

53. AIS GMVO Avtomatizirovannaya Informatsionnaya Sistema Gosudarstvennogo Monitoring Vodnykh ob"Ektov [Automated Information System of State Monitoring of Water Bodies]. [Electronic Resource]. 2020. Available online: https://gmvo.skniivh.ru/ (accessed on 17 March 2020).

54. Critelli, T.; Vespasiano, G.; Apollaro, C.; Muto, F.; Marini, L.; De Rosa, R. Hydrogeochemical study of an ophiolitic aquifer: A case study of Lago (Southern Italy, Calabria). *Environ. Earth Sci.* **2015**, *74*, 533–543. [CrossRef]

55. Raymond, P.A. Temperature versus hydrologic controls of chemical weathering fluxes from United States forests. *Chem. Geol.* **2015**, *458*, 1–13. [CrossRef]

56. Raymond, P.A.; Hamilton, S.K. Anthropogenic influences on riverine fluxes of dissolved inorganic carbon to the oceans. *Limnol. Oceanogr. Lett.* **2018**, *3*, 143–155. [CrossRef]

57. Stackpoole, S.M.; Stets, E.G.; Striegl, R.G. The impact of climate and reservoirs on longitudinal riverine carbon fluxes from two major water-sheds in the Central and Intermontane West. *J. Geophys. Res. Biogeosci.* **2014**, *119*, 848–863. [CrossRef]
58. Pokrovsky, O.S.; Viers, J.; Dupré, B.; Chabaux, F.; Gaillardet, J.; Audry, S.; Prokushkin, A.S.; Shirokova, L.S.; Kirpotin, S.N.; Lapitsky, S.A.; et al. Biogeochemistry of carbon, major and trace elements in watersheds of Northern Eurasia drained to the Arctic Ocean: The change of fluxes, sources and mechanisms under the climate warming prospective. *Comptes Rendus Geosci.* **2012**, *344*, 663-677. [CrossRef]

59. Chupakov, A.V.; Pokrovsky, O.S.; Moreva, O.Y.; Shirokova, L.S.; Neverova, N.V.; Chupakova, A.A.; Kotova, E.I.; Vorobyeva, T.Y. High resolution multi-annual riverine fluxes of organic carbon, nutrient and trace element from the largest European Arctic river, Severnaya Dvina. *Chem. Geol.* **2020**, *538*. [CrossRef]

60. Newcomer Johnson, T.A.; Kaushal, S.S.; Mayer, P.M.; Smith, R.M.; Sivirichi, G.M. Nutrient retention in restored streams and rivers: A global review and synthesis. *Water* **2016**, *8*, 116. [CrossRef]

61. Bloom, P.R. Phosphorus adsorption by an aluminum-peat complex. *Soil Sci. Soc. Am. J.* **1981**, *45*, 267-272.
[CrossRef]

62. Richardson, C. Mechanisms controlling phosphorus retention capacity in fresh-water wetlands. *Science* **1985**, *228*, 1424-1427. [CrossRef] [PubMed]

63. Penn, C.J.; Camberato, J.J. A critical review on soil chemical processes that control how soil pH affects phosphorus availability to plants. *Agriculture* **2019**, *9*, 120. [CrossRef]

64. Prokushkin, A.S.; Hagedorn, F.; Pokrovsky, O.S.; Viers, J.; Kirdyanov, A.V.; Masyagina, O.V.; Prokushkina, M.P.; McDowell, W.H. Permafrost regime affects the nutritional status and productivity of larches in Central Siberia. *Forests* **2018**, *9*, 314. [CrossRef]

65. Stets, E.G.; Striegl, R.G. Carbon export by rivers draining the conterminous United States. *Inland Waters* **2012**, *2*, 177-184. [CrossRef]

66. Giesler, R.; Lyon, S.W.; Morth, C.-M.; Karlsson, J.; Karlsson, E.M.; Jantze, E.J.; Destouni, G.; Humborg, C. Catchment-scale dissolved carbon concentrations and export estimates across six subarctic streams in northern Sweden. *Biogeosciences* **2014**, *11*, 525-537. [CrossRef]

67. Zakharova, E.A.; Pokrovsky, O.S.; Dupré, B.; Zaslavskaya, M.B. Chemical weathering of silicate rocks in Aldan Shield and Baikal Uplift: Insights from long-term seasonal measurements of solute fluxes in rivers. *Chem. Geol.* **2005**, *214*, 223-248. [CrossRef]

68. White, A.F.; Blum, A.E. Effects of climate on chemical weathering in watersheds. *Geochim. Cosmochim. Acta* **1995**, *59*, 1729-1747. [CrossRef]

69. Dessert, C.; Dupré, B.; Gaillardet, J.; Francois, L.M.; Allégre, C.J. Basalt weathering laws and the impact of basalt weathering on the global carbon cycle. *Chem. Geol.* **2003**, *202*, 257-273. [CrossRef]

70. Oliva, P.; Viers, J.; Dupré, B. Chemical weathering in granitic environments. *Chem. Geol.* **2003**, *202*, 225-256.
[CrossRef]

71. Payandi-Rolland, D.; Shirokova, L.S.; Nakhle, P.; Tesfa, M.; Abdou, A.; Causserand, C.; Lartiges, B.; Rols, J.-L.; Guérin, F.; Bénézeth, P.; et al. Aerobic release and biodegradation of dissolved organic matter from frozen peat: Effects of temperature and heterotrophic bacteria. *Chem. Geol.* **2020**, *536*, 119448. [CrossRef]

72. Fraysse, F.; Pokrovsky, O.S.; Meunier, J.-D. Experimental study of terrestrial plant litter interaction with aqueous solutions. *Geochim. Cosmochim. Acta* **2010**, *74*, 70-84. [CrossRef]

73. Antipina, T.G.; Preis, Y.I.; Zenin, V.N. Dynamics of forest vegetation and climate in the southern taiga of Western Siberia in the Late Holocene according to spore-pollen analysis and AMS dating of the peat bog. *Russ. J. Ecol.* **2019**, *50*, 445-452. [CrossRef]

74. Fotiev, S.M. Arctic peatlands of the Yamal-Gydan province of Western Siberia. *Earth's Cryosphere* **2017**, *XXI*, 3-15. [CrossRef]

75. Loiko, S.; Raudina, T.; Lim, A.; Kuzmina, D.; Kulizhskiy, S.; Pokrovsky, O. Microtopography Controls of Carbon and Related Elements Distribution in the West Siberian Frozen Bogs. *Geosciences* **2019**, *9*, 291.
[CrossRef]

76. Panova, N.K.; Trofimova, S.S.; Antipina, T.G.; Zinoviev, E.V.; Gilev, A.V.; Erokhin, N.G. Holocene dynamics of vegetation and ecological conditions in the southern Yamal Peninsula according to the results of comprehensive analysis of a relict peat bog deposit. *Russ. J. Ecol.* **2010**, *41*, 20-27. [CrossRef]

77. Stepanova, V.M.; Pokrovsky, O.S.; Viers, J.; Mironycheva-Tokareva, N.P.; Kosykh, N.P.; Vishnyakova, E.K. Major and trace elements in peat profiles in Western Siberia: Impact of the landscape context, latitude and permafrost coverage. *Appl. Geochem.* **2015**, *53*, 53-70. [CrossRef]

78. Pokrovsky, O.S.; Schott, J.; Kudryavtzev, D.I.; Dupré, B. Basalt weathering in Central Siberia under permafrost conditions. *Geochim. Cosmochim. Acta* **2005**, *69*, 5659-5680. [CrossRef]
79. Gordeev, V.V.; Martin, J.-M.; Sidorov, I.S.; Sidorova, M.V. A reassessment of the Eurasian river input of water, sediment, major elements, and nutrients to the Arctic Ocean. *Am. J. Sci.* **1996**, *296*, 664–691. [CrossRef]

80. Vorobyev, S.N.; Pokrovsky, O.S.; Kolesnichenko, L.G.; Manasypov, R.M.; Shirokova, L.S.; Karlsson, J.; Kirpotin, S.N. Biogeochemistry of dissolved carbon, major, and trace elements during spring flood periods on the Ob River. *Hydrol. Process.* **2019**, *33*, 1579–1594. [CrossRef]

81. Serikova, S.; Pokrovsky, O.S.; Ala-Aho, P.; Kazantsev, V.; Kirpotin, S.N.; Kopysov, S.G.; Krickov, I.V.; Laudon, H.; Manasypov, R.M.; Shirokova, L.S.; et al. High riverine CO₂ emissions at the permafrost boundary of Western Siberia. *Nat. Geosci.* **2018**, *11*, 825–829. [CrossRef]

82. Serikova, S.; Pokrovsky, O.S.; Laudon, H.; Krickov, I.V.; Lim, A.G.; Manasypov, R.M.; Karlsson, J. High carbon emissions from thermokarst lakes of Western Siberia. *Nat. Commun.* **2019**, *10*, 1–7. [CrossRef] [PubMed]

83. Smedberg, E.; Mörth, C.-M.; Swaney, D.P.; Humborg, C. Modeling hydrology and silicon-carbon interactions in taiga and tundra biomes from a landscape perspective: Implications for global warming feedbacks. *Glob. Biogeochem. Cycles* **2006**, *20*. [CrossRef]

84. Guo, L.; Zhang, J.-Z.; Guéguen, C. Speciation and fluxes of nutrients (N, P, Si) from the upper Yukon River. *Glob. Biogeochem. Cycles* **2004**, *18*. [CrossRef]

85. Zakharova, E.A.; Pokrovsky, O.S.; Dupré, B.; Gaillardet, J.; Efimova, L.E. Chemical weathering of silicate rocks in Karelia region and Kola peninsula, NW Russia: Assessing the effect of rock composition, wetlands and vegetation. *Chem. Geol.* **2007**, *242*, 255–277. [CrossRef]

86. Raymond, P.A.; McClelland, J.W.; Holmes, R.M.; Zhulidov, A.V.; Mull, K.; Peterson, B.J.; Striegl, R.G.; Aiken, G.R.; Gurtovaya, T.Y. Flux and age of dissolved organic carbon exported to the Arctic Ocean: A carbon isotopic study of the five largest arctic rivers. *Glob. Biogeochem. Cycles* **2007**, *21*. [CrossRef]

87. Anisimov, O.A.; Shiklomanov, N.I.; Nelson, F.E. Variability of seasonal thaw depth in permafrost regions: A stochastic modeling approach. *Ecol. Model.* **2002**, *153*, 217–227. [CrossRef]

88. Dankers, R.; Burke, E.J.; Price, J. Simulation of permafrost and seasonal thaw depth in the JULES land surface scheme. *Cryosphere* **2011**, *5*, 773–790. [CrossRef]

89. Stendel, M.; Christensen, J.H. Impact of global warming on permafrost conditions in a coupled GCM. *Geophys. Res. Lett.* **2002**, *29*, 1632. [CrossRef]

90. Babkina, E.A.; Leibman, M.O.; Dvornikov, Y.A.; Fakashchuk, N.Y.; Khairullin, R.R.; Khomutov, A.V. Activation of cryogenic processes in Central Yamal as a result of regional and local change in climate and thermal state of permafrost. *Russ. Meteorol. Hydrol.* **2019**, *44*, 283–290. [CrossRef]

91. Botch, M.S.; Kobak, K.I.; Vinson, T.S.; Kolchugina, T.P. Carbon pools and accumulation in peatlands of the former Soviet Union. *Glob. Biogeochem. Cycles* **1995**, *9*, 37–46. [CrossRef]

92. Liss, O.L.; Abramova, L.I.; Avetov, N.A.; Berezina, N.A.; Inisheva, L.I.; Kurnishnikova, T.V.; Sluka, Z.A.; Tolpysheva, T.Y.; Shvedchikova, N.K. *Wetland Systems of West Siberia and Their Importance for Nature Conservation*; Grifi K Publisher: Tula, Russia, 2001; 584p. (In Russian)

93. Gentsch, N.; Mikutta, R.; Alves, R.J.E.; Barta, J.; Capek, P.; Gitte, A.; Hugelius, G.; Kuhry, P.; Lashchinskiy, N.; Palmtag, J.; et al. Storage and transformation of organic matter fractions in cryoturbated permafrost soils across the Siberian Arctic. *Biogeosciences* **2015**, *12*, 4525–4542. [CrossRef]

94. Kaiser, C.; Meyer, H.; Biasi, C.; Rusalimova, O.; Barsukov, P.; Richter, A. Conservation of soil organic matter through cryoturbation in arctic soils in Siberia. *J. Geophys. Res.* **2007**, *112*, 9–17. [CrossRef]

95. Kawahigashi, M.; Kaiser, K.; Kalbitz, K.; Rodionov, A.; Guggenberger, G. Dissolved organic matter in small streams along a gradient from discontinuous to continuous permafrost. *Glob. Chang. Biol.* **2004**, *10*, 1576–1586. [CrossRef]

96. Mergelov, N.S.; Targulian, V.O. Accumulation of organic matter in the mineral layers of permafrost-affected soils of coastal lowlands in East Siberia. *Eurasian Soil Sci.* **2011**, *44*, 249–260. [CrossRef]

97. Oosterwoud, M.R.; Temminghoff, E.J.M.; van der Zee, S.E.A.T.M. Quantification of DOC concentrations in relation with soil properties of soils in tundra and taiga of Northern European Russia. *Biogeosci. Discuss.* **2010**, *7*, 3189–3226. [CrossRef]

98. Anisimov, O.; Kokorev, V.; Zhil’tsova, Y. Temporal and spatial patterns of modern climatic warming: Case study of Northern Eurasia. *Clim. Chang.* **2013**, *118*, 871–883. [CrossRef]

99. Romanovsky, V.E.; Drozdov, D.S.; Oberman, N.G.; Malkova, G.V.; Kholodov, A.L.; Marchenko, S.S.; Moskalenko, N.G.; Sergeev, D.O.; Ukrainteva, N.G.; Abramov, A.A.; et al. Thermal state of permafrost in Russia. *Permafr. Periglac. Proc.* **2010**, *21*, 136–155. [CrossRef]
100. Vasiliev, A.A.; Streletskaya, I.D.; Shirokov, R.S.; Oblogov, G.E. Evolution of cryolithozone of coastal zone of western Yamal during climate change. *Kriosf. Zemli* **2011**, *2*, 56–64.

101. Ernakovich, J.; Hopping, K.A.; Berdanier, A.B.; Simpson, R.T.; Kachergis, E.J.; Steltzer, H.; Wallenstein, M.D. Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. *Glob. Chang. Biol.* **2014**, *20*, 3256–3269. [CrossRef] [PubMed]

102. Karlsson, J.; Serikova, S.; Vorobyev, S.N.; Rocher-Ros, G.; Denfeld, B.; Pokrovsky, O.S. Carbon emission from Western Siberian inland waters. *Nature Comm.* (submitted).

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

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).