3.1 C, N and P concentrations in RSM and their link to seasons and watershed size

Mean bulk RSM concentration in the WSL river waters did not depend on the season of open-water period of the year and was equal to 7.1±3.9, 8.1±4.1 and 7.0±3.7 mg L⁻¹ in spring, summer and autumn, respectively (Table 1). The RSM concentrations weakly depended on the size of the watersheds (S_watershed) with a negative relationship in autumn (R²=0.33, p<0.05; Fig. S1a in the Supplement). Further, the RSM concentrations increased with permafrost coverage and latitude (R²=0.56 and 0.41), although this was visible only in autumn (Fig. S1b, c; Table S2). The sporadic permafrost zone exhibited the highest RSM concentration in summer (Fig. S1d). Finally, there was no correlation (p>0.05) between lake, bog or forest coverage and the RSM concentrations (R²<0.2; see also Table S2). For RSM concentrations, statistically significant differences between different permafrost zones, notably between permafrost-free and permafrost-bearing regions, were evidenced in summer and autumn using the Kruskal–Wallis and Mann–Whitney tests (Table S3).

The concentrations of C, N and P in WSL rivers averaged over three seasons were equal to 15.3±9.7 %, 1.2±0.9 % and 0.49±0.42 % in mass of RSM (1.05±0.805, 0.083±0.066 and 0.035±0.036 mg L⁻¹ in the river water). The watershed size sizably affected the C concentration; there was a power-law decrease of C with the size of watershed (R²=0.28, 0.47 and 0.25 in spring, summer and autumn, respectively Fig. 2a), but there was no relationship with the N and P concentrations in RSM (R²<0.2, Fig. 2b, c). Generally, a 2- to 3-fold decrease in C_org, from ca. 20 %–30 % in rivers with S_watershed<100 km² to C_org=5 %–10 % in rivers with S_watershed>10 000 km² was observed. The C:N ratio of RSM was independent of the watershed size in spring but decreased 2–3 times with an increase in S_watershed (R²=0.4) in summer and autumn (Fig. 2d).

Figure 2Particulate (>0.45 µm) C (a), N (b) and P (c) concentrations in the RSM (%), and C:N ratio (d) in RSM as a function of river watershed size. The solid lines represent a power-law fitting of the data with regression coefficients shown for each season in corresponding panels. Only the curves with R²>0.3 are depicted.

Download

The inter-annual variations of suspended nutrient concentration in WSL rivers were of secondary order importance when compared to season and watershed size control. We did not find any inter-annual differences (at p<0.05) in RSM concentrations and P concentrations in RSM collected in June and August in 2014, 2015, and 2016 for the same eight rivers (Agan, Tromyogan, Pyakupur, Ayvasedapur, Purpe, Yamsovey, Pur and Taz; Table S1).

3.2 Role of permafrost distribution and landscape parameters for C, N and P concentrations and fraction of particulate nutrients

There was a local maximum of C and N concentrations in isolated and sporadic permafrost zone (Fig. 3a, b, d, e), which was not seen for P (Fig. 4c, f). Overall, the differences in C and N concentrations in RSM among different permafrost zones were significant, as verified by the non-parametric Kruskal–Wallis H test ($\mathrm{0.005}<p<\mathrm{0.05}$), while the difference in P concentrations between permafrost zones was not significant (p>0.05, see Table S3c, d). Specifically, the C demonstrated a maximum concentration (significant at p<0.02 during all three seasons) at 62–64^∘ N (Fig. S2a). The latitude did not impact N and P concentrations in RSM per se (Fig. S2b, c). However, significant differences between adjacent permafrost zones were evidenced by C and N in summer and autumn (Table S3d).

Figure 3Box plot of first and third quartiles (25 % and 75 %) of C (a), N (b) and P (c) concentrations in RSM (%) in five permafrost zones over three seasons. The C, N and P concentrations in the river water are shown in panels (d), (e) and (f), respectively, and a C:N ratio is shown in (g).

Download

The landscape parameters of the watershed (bogs, lakes and forest coverage) sizably affected (p<0.05) suspended C and N. Bogs and lakes in the watershed increased the concentrations of C and N in RSM, whereas forest generally decreased C in RSM (Fig. 4a–c for C, and Fig. S3a–c for N). This increase in C and N (%) with bog and lake coverage and a decrease in forest presence was mostly visible in summer and autumn. The increase in lake coverage of the watershed led to a decrease in P concentrations in RSM in summer and autumn (R²=0.31 and 0.22, respectively; Fig. S3d-e-f) that was especially visible in autumn in the permafrost-free zone ($R=-\mathrm{0.88}$; Table S2). During this period, the P concentrations in RSM positively correlated with the presence of forest in the permafrost zone (R=0.60; Table S2).

Figure 4The dependence of C concentrations in RSM (%) on the coverage of watershed by bogs (a), lakes (b) and forest (c). The solid lines represent second-degree polynomial fitting of the data with regression coefficients shown for each season in corresponding panels.

Download

The Mann–Whitney U test for the impact of watershed parameters demonstrated significant differences in C and N concentrations (all seasons) and P concentrations (summer baseflow) between watersheds having <10 % and >10 % lake coverage (see Table S3e). The differences were also observed among watersheds with <50 % and >50 % of bogs for C (all seasons) and N (summer and autumn; Table S3f). Finally, the forest coverage (<30 % and >30 %) exhibited a significant effect on C and N (all seasons) and P (autumn baseflow; see Table S3g).

The share of particulate carbon versus total carbon (dissolved + particulate C) did not demonstrate any significant dependence on S_watershed, bogs, forest and permafrost proportions on the watershed (R²<0.3, not shown). However, there was a localized maximum of the particulate carbon fraction around 64^∘ N within the isolated-to-sporadic permafrost zone (Fig. 5a and c). The presence of lakes sizably increased the particulate over total transport of C in rivers (R²=0.52 and 0.32 in spring and summer, respectively, Fig. 5b). The share of particulate phosphorus versus total phosphorus ranged from 10 % to 90 %. It did not demonstrate any link to the size of river watershed, percentage of forest and bogs, and type of permafrost distribution (not shown).

Figure 5Fraction of particulate OC of total (dissolved + particulate) form plotted as a function of latitude (a), lake fraction on the watershed (b) and a box plot of fractions for five permafrost zones (c). The solid lines in (a) and (b) represent a second-degree polynomial (a: autumn) and linear (b: spring) fitting of the data, with regression coefficients equal to 0.32 and 0.52, respectively.

Download

3.3 C, N, P and RSM export fluxes by WSL rivers

Based on available hydrological data, we calculated open-water-period fluxes of C, N and P in WSL rivers. This analysis takes into account the spatial and temporal variability of river discharge, performed using hydrological approaches elaborated for the dissolved (<0.45 µm) fraction of the river water (Pokrovsky et al., 2015, 2016). The seasonal fluxes of C, N, P and RSM export by WSL rivers were calculated separately for the spring (May and June), summer (July, August and September) and autumn periods (September–October) for each 2^∘-wide latitudinal belt of the full WSL territory, following the approach developed for dissolved C and major and trace elements in the river water (Fig. S4). These three seasons of the open-water period represent by far the largest contribution to overall annual element and RSM yield, following the results for other Arctic rivers (McClelland et al., 2016). Thus, six ice-covered months (November to April) represent only 12 % of annual POC export flux by the Ob River. Based on the results of three main seasons, open-water-period export fluxes of C, N, P and RSM were calculated (Fig. 6). There is a clear maximum of C and N export fluxes at the beginning of permafrost appearance, in the isolated-to-sporadic permafrost zone. The obtained particulate C and N yields are comparable with other Siberian rivers. For the two largest WSL rivers, Pur and Taz, we found May-to-October export fluxes of 69 and 80 kg C km² yr⁻¹, which are lower than the annual POC yield of the Ob River (191 kg C km⁻² yr⁻¹) but similar to that of the Yenisey River (103 kg C km⁻² yr⁻¹; McClelland et al., 2016).

Figure 6Total open-water-season (May to October) fluxes of particulate C (a), N (b) and P (c), and suspended matter (d) in five permafrost-free and four distinct permafrost zones of WSL (box plots of first and third quartiles). There is a clear maximum of C and N export at the beginning of permafrost appearance, in the isolated-to-sporadic permafrost zone.

Download

4.1 Concentrations of C, N and P in the RSM and impact of the watershed size

The RSM values in WSL rivers (2 to 18 mg L⁻¹) are similar to other boreal rivers of low runoff, which drain peatlands such as Severnaya Dvina (2.3 to 16 mg L⁻¹; Pokrovsky et al., 2010), but are lower than the Ob River itself (around 30 mg L⁻¹; Gebhardt et al., 2004) and other big rivers of the Kara Sea basin (average 22 mg L⁻¹; Gordeev et al., 1996). The POC values of the WSL rivers (0.5 to 3.0 mg L⁻¹ POC) are consistent with recent data on WSL river transects sampled in 2015 (Vorobyev et al., 2017) and are in agreement with those of the Ob–Taz River confluence measured in June (1.3 mg L⁻¹; Gebhardt et al., 2004), the Ob River at Salekhard in May through October (0.8 to 2.4 mg POC L⁻¹; Le Fouest et al., 2013), the low reaches of the Ob River (1.2 to 2.4 mg POC L⁻¹; McClelland et al., 2016), the mean multi-annual values of POC in subarctic rivers of northern Eurasia draining peatlands (3.2, 0.3, 0.9 mg POC L⁻¹ for Severnaya Dvina, Pechora and Ob; as compiled in Gordeev et al., 1996) and the Lena River basin (0.5 mg L⁻¹; Kutscher et al., 2017).

However, the C_org concentrations in RSM of WSL rivers (5 % to 40 %), notably in small- and medium-sized (<10 000–100 000 km²) ones, are an order of magnitude higher than those in other world rivers, which drain mineral substrates (typically 1 % C_org in RSM; Meybeck, 1993) and are significantly higher than the values of the Siberian rivers (2.3 %, 3.6 %, 5.8 % and 3.0 % for Ob, Yenisey, Lena and Kolyma, respectively; Gordeev and Kravchishina, 2009). For example, the typical concentrations of C_org in RSM of large (S_watershed>100 000 km²) central Siberian rivers that drain larch forest are only 0.4 % to 0.5 % (Pokrovsky et al., 2005). The C_org concentrations in the RSM of Severnaya Dvina River (which has a sizeable proportion of bogs and lakes within its watershed compared to WSL rivers) are 2.7±0.7 % in May and 4.8±1.1 % in August (Savenko et al., 2004). The N_org content in RSM ranges from 0.3 % to 1.8  % (0.05 to 0.2 mg particulate N_tot L⁻¹), which is much higher than that in sedimentary rocks (0.05 % to 0.06 %; Houlton et al., 2018) but is comparable with the value reported for the freshwater part of Ob River estuary (0.16 mg N L⁻¹; Gebhardt et al., 2004), the Ob River at Salekhard in May to October (0.1 to 0.3 mg PON L⁻¹; PON – particulate organic nitrogen; Le Fouest et al., 2013), the Yukon River (0.14±0.09 mg particulate N L⁻¹; Guo and MacDonald, 2006) and small rivers in the North Slope, Alaska (0.05 to 0.6 mg PON L⁻¹; McClelland et al., 2014).

High concentrations of C (and N) in the RSM of WSL rivers may stem from the organic nature of soils that prevail on river watersheds. The Histosols, one of the dominant soil groups of WSL, are capable of providing a sizeable amount of organic particles given the higher susceptibility of peat to physical disintegration compared to mineral soils. The enrichment of the river water in C-rich particles may occur at both the river bank (especially in small rivers flowing through the wetlands) and within the extensive floodplains via the remobilization of organic-rich sediments during high flow periods.

The concentrations of C and N in RSM decreased with increase in S_watershed, thereby illustrating the importance of organic particles in small rivers draining the peatlands and the role of mineral matter from bank abrasion in larger rivers. The impact of watershed size is more significant for C than for N. Presumably this is because N is more affected by autochthonous processes, and because particulate N may partly be generated from phytoplankton and macrophytes in the river. Small rivers (S_watershed<100–1000 km²) exhibited the largest scatter in particulate C, N and P concentrations. This is probably due to multiple sources of POM and the very short transit time in the watershed that results in fast responses of river particulate load to minor variations in surface hydrology, including high sensitivity to local storm events.

The decrease of C:N in the RSM from small to large rivers likely reflected a shift in the main origin of suspended matter, from peat in small rivers to more lithogenic (deep soil) in large rivers. This was mostly visible in summer and autumn; in spring the rivers exhibited a very homogeneous C:N signature that may be linked to a dominant source of RSM from bank abrasion and sediment transport as well as deposition within the riparian zone. In fact, the flood plain of the Ob River and other rivers of the WSL extend to more than 10 times the width of the main channel (Vorobyev et al., 2015). Note that the C:N ratio in large rivers (>100 000 km²) approaches that of average sedimentary rocks (8.1; Houlton et al., 2018). In this regard, highly homogeneous C:N ratios in the particulate load of Arctic rivers (7 to 18 for Mackenzie, Yukon, Kolyma, Lena, Yenisey and Ob regardless of season; McClelland et al., 2016) are interpreted as the mixture of deep soil sources where $\mathrm{C}:\mathrm{N}<\mathrm{10}$ (Schädel et al., 2014) and upper organic-rich horizons of soils have elevated C:N (Gentsch et al., 2015). The Ob River demonstrates the youngest POC of all Arctic rivers (−203 ‰ to −220 ‰ Δ¹⁴C; McClelland et al., 2016), which certainly indicates a relatively fresh (ca. 1000–2000 years old) origin of particulate carbon that is presumably from intermediate peat horizons.

We believe that variations in C:N in RSM reflect different sources of organic material feeding the river depending on seasons and latitudes. A compilation of C:N ratios in peat and mineral horizons as well as in thermokarst lake sediments for four main sites of latitudinal transect considered in this study is given in Fig. S5 of Supplement. The range of C:N values in RSM rivers (10 to 20) is closer to that in sediments of thermokarst lakes (20 to 30). Note that the resuspension of sediments may be an important source of water column POC (Yang et al., 2016).

The minerotrophic bogs, which are mostly linked to rivers via hydrological networks, have a C:N ratio in upper peat horizons ranging from 24 to 28. In mineral soils of the region, the C:N range is between 10 and 15, regardless of latitude, from the tundra-situated Taz River riparian zone to the taiga-situated middle channel of the Ob River. For upper organic horizons the C:N is always higher than the bottom mineral horizons. The old alluvial deposits of the Pyakupur River (discontinuous permafrost zone) had only 0.2 % of POC, with C:N equal to 6. Overall, there is an enrichment in N relative to C in the course of water transport of organic and organic-mineral solid particles from soils and riparian deposits to the river water.

Another important observation following from the consistently low C:N ratios of RSM across rivers of various size and climatic zones is that the flocculation and aggregation of riverine dissolved organic matter (DOM) in lotic waters of Siberian lowlands may be quite low. Further, the absence of a significant relationship between the lake proportion at the watershed and the C:N ratio implies the negligible impact of DOM coagulation due to photolysis (von Wachenfeldt et al., 2008; von Wachenfeldt and Tranvik, 2008) or bacterial activity (von Wachenfeldt et al., 2009), with the subsequent transformation of coagulation products (Kortelainen et al., 2006b, 2013) as it is known in European humic lakes. Note that, because the range of C:N in RSM of WSL is far from that reported for DOM in soil solutions of boreal taiga (ca. 100, Ilina et al., 2014; 40 to 80, Dymov et al., 2013) and humic (peatland) lakes (>50, Chupakov et al., 2017), the coagulation of DOM from soil waters producing particles in the rivers is also unlikely.

4.2 A maximum of C and N in the isolated and sporadic permafrost zone and the impact of river watershed characteristics

Complementary to previous results on dissolved (<0.45 µm) C and N concentrations in WSL rivers acquired by Frey et al. (2007a) and Vorobyev et al. (2017) that demonstrated the weak or nonexistent impact of permafrost on DOC and dissolved organic nitrogen (DON), the particulate C and N were affected by the presence of permafrost in summer and autumn but were not affected by its presence in spring. Moreover, during freshet the permafrost distribution did not influence the bulk RSM concentration in WSL rivers. This strongly implies that the delivery of RSM in rivers, and its chemical composition, is tightly linked to the thickness of the active layer and are limited by the transport of soil particles over the suprapermafrost flow to the river channel. This thickness is highest in September at the end of the active season. In agreement with this, the C and N demonstrated a maximum concentration and export fluxes at 62–64^∘ N, in the isolated-to-sporadic permafrost zone, that were most visible during summer and autumn (Figs. 3a–b and 6a–b). This latitudinal belt can be considered as a large-scale thawing front for the frozen peat that corresponds to the southern boundary of permafrost persistence. It is important to note that WSL rivers exhibit maximum CO₂ emission fluxes at the isolated-to-sporadic permafrost belt (Serikova et al., 2018), which could be linked to the strong processing of POC and PON in the water column of WSL rivers. Interestingly, that rate of POC biodegradation, leading to potential CO₂ emissions, sizably exceeds that of DOC in boreal humic waters (Attermeyer et al., 2018). Furthermore, a maximum percentage of particulate C over total C (suspended + dissolved) was also in the isolated and sporadic permafrost zones in spring; this maximum shifted to the sporadic permafrost zone in summer and moved northward to the discontinuous permafrost zone in autumn (Fig. 5c). We believe that this corresponds to a progressive increase in the thickness of the active layer that controls the degree of peat and mineral particles leaching from the soil profile to the river. The thickness of this layer increases from spring to autumn, and more importantly, it moves northward during this period (Trofimova and Balybina, 2014). The enhanced mobilization of nutrients at the “hot spot” of permafrost thaw in frozen peat landscapes was recently demonstrated on a local scale in western Siberia (Loiko et al., 2017).

The impact of watershed characteristics on particulate C and N was clearly pronounced, with increased C and N concentrations in RSM where there were increased bog and lake proportions and decreased C and N concentrations where there was increasing forest coverage. The stronger impact of lakes compared to bogs on C concentration in RSM suggests that the generation of C-rich particles occurs more efficiently in large water bodies than in stagnant shallow water bodies. Given the very short transit time of water from the surrounding peat to the lakes via suprapermafrost flow (Ala-aho et al., 2018a, b; Raudina et al., 2018), the allochthonous chromophoric DOM-rich material from peat soil water that arrives to the lakes may be subjected to fast degradation and coagulation, as shown in Scandinavian lakes (Kortelainen et al., 2006b; von Wachenfeldt and Tranvik, 2008). Second, the peat abrasion at the border of the thermokarst lakes and thaw ponds, which are highly abundant in the territory (Polishchuk et al., 2017, 2018), occurs due to wave erosion and thermo-abrasion (Shirokova et al., 2013; Manasypov et al., 2015). The physical disintegration of peat at the lake coast likely generates a large amount of suspended organic-rich material that can be exported to hydrological networks during, for example, lake drainage or already existing connecting channels (Kirpotin et al., 2008, 2011). Note that the maximal lake coverage of the WSL territory is in the 63 to 64^∘ N latitudinal belt (Polishchuk et al., 2017), where maximum C and N concentrations and RSM export fluxes also occur. Because the majority of thermokarst lakes are isolated water bodies without inlets and outlets, this connectivity is achieved via water movement along the permafrost table in the thawed active layer (Raudina et al., 2018), in the form of the so-called suprapermafrost flow between peat bogs, lakes and rivers.

Finally, for particulate P, neither its concentration nor the particulate fraction were affected by permafrost distribution, probably due to the various processes of biological uptake and the mineral precipitation controlling the removal of P, both in the soil profile and in the river water. For example, lakes and bogs retained particulate P, similar to that of dissolved P, which is in agreement with global assessments (Bouwman et al., 2013), P behavior in European northern wetlands and lakes (Lidman et al., 2014), and recent results on dissolved P in the WSL rivers (Vorobyev et al., 2017).

4.3 Mechanisms of RSM generation and prospectives for climate warming in western Siberia

A framework of particulate C, N and P generation in WSL rivers across the permafrost gradient is shown in Fig. 7. We suggest that the concentration and export fluxes of suspended particles depend on both the supply and losses in the catchments. The sources of suspended particles in WSL rivers include the following: (i) vegetation litter, which is washed by surficial flow to the river, especially in spring; (ii) surface (peat) soil horizons, which are also most active in spring, especially in the north; (iii) deep peat and mineral horizons, which provide the particles via bank abrasion in spring and via suprapermafrost flow in summer and autumn; (iv) lake coastal abrasion due to wave erosion; and finally, (v) the autochthonous organic debris of macrophytes, periphyton and phytoplankton, whose contribution is maximal in summer and autumn. A non-steady-state physical erosion of peat soils in WSL provides maximum particulate nutrients within the most fragile zone of actively thawing permafrost between 62 and 64^∘ N of the isolated-to-sporadic permafrost zone. The maximal thickness of the active layer progressively moves north during the active season, thereby leading to maximal export of particulate C, N and P at the thawing front. However, we also suggest that part of the difference in mobilized particulates is masked by retention in recipient waters. The transit time of water and particles in the southern WSL rivers is much longer than that in northern rivers (Ala-aho et al., 2018a, b); hence the biological uptake mechanisms (Richardson et al., 2013; Attermeyer et al., 2018), together with physiochemical processes such as the photodegradation of POC (Mayer et al., 2006; Riggsbee et al., 2008) or cryo-coagulation (Pokrovsky et al., 2018), have sufficient time to act on suspended matter of soil and shallow subsurface waters and to remove the nutrients from the river water as well. In rivers of the continuous permafrost zone, a relatively small stock of nutrient-rich particles within the soil profile and on soil surface (as plant litter) is largely compensated for by a more rapid flushing and shorter travel time through soils and rivers as well as lower microbial and phytoplankton activity. As a result, the zone of isolated-to-sporadic permafrost exhibits both the maximal release of soil particles and minimal uptake by in-stream processes. Further to the north, shallow unfrozen peat depth and low biomass cannot supply sufficiently high suspended nutrients, and the particulate transport of C and N decreases. In contrast, for P, opposite gradients in supply and in-stream removal may cancel out the net effect of temperature and permafrost on suspended P in the river water.

Figure 7A cartoon of spatial and temporal partitioning of particulate nutrients in WSL rivers across the permafrost gradient. The panels (a), (b), (c) and (d) represent the main sources (a: lakes and bogs in summer, and b: alluvial deposits in spring) and sinks (c: photodegradation and bio-degradation, and d: uptake by taiga forest) of particulate nutrients in WSL rivers. The panel (e) depicts the spatial gradient of C and N in RSM occurring in spring (blue line) and autumn (red line). A non-steady-state physical erosion of peat soils in WSL provides the maximum of particulate nutrients within the zone of most “fragile” actively thawing permafrost. The maximal thickness of the active layer progressively moves to the north during the active season, thus leading to the maximal removal of particulate C, N and P at the thawing front.

Download

Based on these results we can speculate on the conditions following warming and permafrost thaw. The drainage of the lakes and the colonization of the bogs by forest is a very common scenario of landscape evolution in western Siberia under ongoing climate warming (Kirpotin et al., 2009, 2011). Scenarios of thermokarst lake evolution under climate warming and permafrost thaw in western Siberia include (1) the draining of large thermokarst lakes into hydrological network, which is especially pronounced in discontinuous permafrost zone (Smith et al., 2005; Polishchuk et al., 2014) and (2) the appearance of new depressions, subsidence and small thaw ponds (<100–1000 m²), which is evidenced across all permafrost zones in this region (Shirokova et al., 2013; Bryksina and Polishchuk, 2015). In terms of landscape change, the area of hollows and subsidence will increase, and the coverage of palsa by mounds and polygons will decrease (Moskalenko, 2012; Pastukhov and Kaverin, 2016; Pastukhov et al., 2016). From a short-term prospective (10–50 years), assuming a soil temperature rise of 0.15^∘ to 0.3^∘C per 10 years in WSL (Pavlov and Malkova, 2009; Anisimov et al., 2012), the northern part of the WSL (discontinuous and continuous permafrost zones) will transform into sporadic and isolated permafrost zones (Anisimov and Reneva, 2006). This will lead to an increase in C and N concentrations in the RSM, C and N particulate export yield of the watershed and an overall increase in the particulate transport versus dissolved transport of C and P. Given the contemporary maximum of C and N at the permafrost thawing front, this increase may be 2 fold. However, from a longer-term prospective (50–100 years), even the continuous permafrost zone may disappear (Romanovsky et al., 2008; Nadyozhina et al., 2008); this will decrease the particulate C and N concentrations in the northern rivers and, consequently, their export to the coastal zone of the Kara Sea. Judging from the actual difference in nutrient concentrations and fluxes among adjusting permafrost zones, this decrease may be of around a factor of 2 to 3. Furthermore, from the same long-term prospective, the drainage of lakes and disappearance of bogs due to the colonization of northern palsa by forests (Anisimov et al., 2011; Anisimov and Sherstiukov, 2016; Kirpotin et al., 2008, 2009, 2011) should lead to a further decrease in the particulate nutrient load of WSL rivers.