Next Article in Journal
Measuring Social Vulnerability to Climate Change at the Coast: Embracing Complexity and Context for More Accurate and Equitable Analysis
Next Article in Special Issue
The Influence of Seasonal Variability of Eutrophication Indicators on Carbon Dioxide and Methane Diffusive Emissions in the Largest Shallow Urban Lake in China
Previous Article in Journal
Divergent Geochemical Pathways of Carbonate Aquifer Evolution in a Classic Karst Terrain: (1) Polygenetic Cave Development Identified Using Longitudinal Groundwater Geochemistry
Previous Article in Special Issue
Large Uncertainties in CO2 Water–Air Outgassing Estimation with Gas Exchange Coefficient KT for a Large Lowland River
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissolved Carbon Concentrations and Emission Fluxes in Rivers and Lakes of Central Asia (Sayan–Altai Mountain Region, Tyva)

by
Arisiya A. Byzaakay
1,2,
Larisa G. Kolesnichenko
1,
Iury Ia. Kolesnichenko
1,
Aldynay O. Khovalyg
2,
Tatyana V. Raudina
1,
Anatoly S. Prokushkin
3,
Inna V. Lushchaeva
1,
Zoia N. Kvasnikova
4,
Sergey N. Vorobyev
1,
Oleg S. Pokrovsky
5,6,* and
Sergey Kirpotin
1
1
Bio-Clim-Land Centre of Excellence, Tomsk State University, Lenin Avenue, 36, 634050 Tomsk, Russia
2
Research Organization Department, Tuvan State University, Lenina Street, 36, 667000 Kyzyl, Russia
3
V.N. Sukachev Institute of Forest SB RAS, Akademgorodok 50/28, 660036 Krasnoyarsk, Russia
4
Geology and Geography Faculty, Tomsk State University, Lenin Avenue, 36, 634050 Tomsk, Russia
5
Geosciences and Environment Toulouse, UMR 5563 CNRS, 31400 Toulouse, France
6
N. Laverov Federal Center for Integrated Arctic Research of the Ural Branch of the Russian Academy of Sciences, Northern Dvina Emb. 23, 163000 Arkhangelsk, Russia
*
Author to whom correspondence should be addressed.
Water 2023, 15(19), 3411; https://doi.org/10.3390/w15193411
Submission received: 3 September 2023 / Revised: 20 September 2023 / Accepted: 23 September 2023 / Published: 28 September 2023
(This article belongs to the Special Issue Recent Progress in CO2 Emission from the World’s Rivers)

Abstract

:
The carbon (C) cycle in inland waters, including carbon concentrations in and carbon dioxide (CO2) emissions from water surfaces, are at the forefront of biogeochemical studies, especially in regions strongly impacted by ongoing climate change. Towards a better understanding of C storage, transport and emission in Central Asian mountain regions, an area of knowledge that has been extremely poorly studied until now, here, we carried out systematic measurements of dissolved C and CO2 emissions in rivers and lakes located along a macrotransect of various natural landscapes in the Sayan–Altai mountain region, from the high mountains of the Western Sayan in the northwest of Tyva to the arid (dry) steppes and semideserts in the intermountain basins in the southeast of Tyva on the border with Mongolia. New data on major hydrochemical parameters and CO2 fluxes (fCO2) gathered by floating chambers and dissolved organic and inorganic carbon (DOC and DIC, respectively) concentrations collected over the four main hydrological seasons allowed us to assess the current C biogeochemical status of these water bodies in order to judge possible future changes under climate warming. We further tested the impact of permafrost, river watershed size, lake area and climate parameters as well as ‘internal’ biogeochemical drivers (pH, mineralization, organic matter quality and bacterial population) on CO2 concentration and emissions in lakes and rivers of this region and compared them with available data from other subarctic and mountain settings. We found strong environmental control of the CO2 pattern in the studied water bodies, with thermokarst lakes being drastically different from other lakes. In freshwater lakes, pCO2 negatively correlated with O2, whereas the water temperature exerted a positive impact on pCO2 in large rivers. Overall, the large complexity of counteracting external and internal drivers of CO2 exchange between the water surfaces and the atmosphere (CO2-rich underground DIC influx and lateral soil and subsurface water; CO2 production in the water column due to dissolved and particulate OC biodegradation; CO2 uptake by aquatic biota) precluded establishing simple causalities between a single environmental parameter and the fCO2 of rivers and lakes. The season-averaged CO2 emission flux from the rivers of Tyva measured in this study was comparable, with some uncertainty, to the C uptake fluxes from terrestrial ecosystems of the region, which were assessed in other works.

1. Introduction

The carbon (C) cycle in inland waters, including the dissolved and particulate concentrations of carbon and the carbon dioxide (CO2) emissions, are at the forefront of biogeochemical studies, especially in the regions most sensitive to ongoing climate change such as boreal and subarctic zones [1,2,3,4,5,6,7,8,9,10] or the tropical/equatorial belt [11,12,13,14]. In contrast to these numerous works, C storage, transport and emission in central continental, mountain and arid regions remain strongly understudied, due either to limited access and logistics or to the still underestimated potential role of these remote territories in C cycling in inland waters. This is especially true for the Central Asian mountain system encompassing the Tibet, Himalaya, Pamir, Altai and Sayan regions. Exceptions are the thorough works on DOC, DIC and POC fluxes in the Himalayan rivers [15,16,17] and C concentrations and fluxes in thermokarst lakes of the Tibetan Plateau [18,19,20,21,22,23]. However, the northern part of the Central Asian Mountain System remains virtually unexplored from the viewpoint of hydrochemistry and carbon balance in its rivers and lakes.
The Altai–Sayan Mountain system is a specific inland region of the northern part of Central Asia which covers the territory of four countries: Russia, Mongolia, China and Kazakhstan [24]. In Russia, it is located within the Tyva Republic, Altai Republic, Buryatia, Khakassia, and Krasnoyarsk Krai. It is characterized by a high degree of continentality and aridity as well as the highest level of endemicity under a huge variety of ecosystems and landscapes, many of which are vulnerable to climatic changes [25]. The water bodies such as rivers and lakes are of particular importance for the sustainability, ecosystem services and conservation of the biodiversity of this arid region [26,27,28,29,30,31]. The particular interest in the Altai–Sayan region is that currently occurring climate changes in this territory contradict existing world models and forecasts. In particular, in the Tyva Republic, our group reported the first ever observed occurrence of the natural phenomenon of ‘greening’ or afforestation of steppes and bare sands [32,33]. This finding contrasts with a number of prediction models that stipulated further progressive drying of the arid Altai and Sayan regions (ASR), such as Khakassia, Tyva and Mongolia [34,35]. The other existing assessments of future climates in arid regions of Eurasia also predict an increase in aridity and even propose the propagation of desertification in the steppe regions [36,37,38]. For instance, in the Tyva Republic, the areal extension of steppe ecosystems, including dry steppes, is predicted to increase by 20–65%, whereas in semideserts, the projected increase is predicted to be several hundred percent more compared to their current extent [24].
In contrast to the above-mentioned predictions, some other studies suggest an increase in the climatic instability, such as drastic transformation of the main atmospheric circulation, in the south Siberia regions and at the border with Central Asia, which could lead to a decrease in the role of atmospheric transport from Atlantic regions, resulting in the blocking of anticyclones and an increase in meridional atmospheric transport [39,40]. This process, in turn, could cause many catastrophic weather events, such as the dramatic appearance of rainstorms, hurricanes and water floods, together leading to the progressive humidification of arid, previously dry regions [41]. Such catastrophic weather events have been observed over the last few years in Khakassia, Tuva and Mongolia [42,43,44].
It is clear that the ongoing climate changes can drastically impact the current hydrochemical, biogeochemical and hydrological status of the inland rivers and lakes of the region. For example, it is known that lakes at high altitudes and cold climates are particularly sensitive to global change [45,46] and modifications in carbon biogeochemistry, including the status of dissolved organic matter (DOM), and C emission fluxes may strongly alter the role of these lakes in the global C cycle [47,48,49]. Among the consequences of the progressive humidification of Central Asian arid regions under ongoing climate change, the treeline shift (greening of upland and advancement of treeline) can strongly affect the biogeochemical functioning of lakes, given that 15% of all lakes globally are located at elevations 500 m or more above sea level [50]. Further, existing paleo-reconstructions suggest that lake productivity [51] and biogeochemistry [52] are sensitive to changes in DOM input linked to treeline position. Therefore, thorough assessment of today’s status of the C biogeochemical cycle (concentration, emission from the water surfaces) is necessary to be able to judge the possible future changes in Central Asian regions induced by climate instability.
Over the past decade, a few studies in this area addressed the hydrochemical [53,54,55,56,57,58], hydrological and hydrographic [59,60,61] status of the water bodies, including the balneological aspect of mineral springs [62,63,64,65]. There are also some studies conducted on the hydrochemistry of water bodies in adjacent regions of Mongolia [66]. However, the status of aquatic C and CO2 emissions from the water surfaces of the Altai–Sayan region remain unknown. Towards filling this gap in knowledge, here, we assessed the concentration of dissolved organic and inorganic carbon and C emissions in lakes and rivers of the Tyva region by selecting large and small water bodies affected by permafrost to different degrees. A novel physicogeographical and climatic transect of inland water bodies, which we implemented in this study, extends from the northwest (the highlands of the western Sayan) to the southeast (the semideserts of the Ubsunur Depression on the border with Mongolia). This transect comprises a large variety of natural ecosystems and landscapes of the region, from glacial–nival high mountain belts to foothill taiga forests, intermountain basins with mixed herb and steppe ecosystems and dry semideserts. As a working hypothesis, we anticipated strong environmental control of C biogeochemical parameters of lakes and rivers, including climate, altitude, permafrost extent and size of the watershed as the main ‘external’ drivers of C concentration and emission. We also assessed the link between carbon parameters of the water bodies and possible ‘internal’ biogeochemical drivers (hydrochemical parameters) such as pH, mineralization, quality of dissolved organic matter and bacterial abundance. We tested these controls across the four main hydrological seasons (spring, summer, autumn and winter) in five lakes and 15 rivers of different sizes and landscape contexts.

2. Materials and Methods

2.1. Water Bodies of the Tyva Republic

We visited 20 water bodies during four hydrological seasons (autumn 2021–summer 2022), as shown in Figure 1 and Figure 2. The main physicogeographical parameters of the rivers and lakes are listed in Table 1 and Table 2 and are described in detail in the Appendix A. The Republic of Tyva is located in the Sayan–Altai Mountain region, between 50 and 53° N and between 88 and 99° E; the elevations range from 2000 to 500 m, creating a large variety of landscapes, from high-altitude belts to basins with steppes and semideserts. The climate of the region is continental, with cold, long winters and hot summers; the mean monthly temperatures range from −41 in January to +35 °C in July; the precipitation is low (115 to 350 mm y−1), and about 70% of it falls during the warm season of the year [67].
During the period of our study, the temperature ranged from −28 °C to 21 °C, and the precipitation was equal to 322 and 144 m during 2021 and 2022, respectively. (http://www.pogodaiklimat.ru/history/36096_2.htm, accessed on 23 September 2023). The water objects sampled in this work include numerous water bodies mainly belonging to the Yenisei River basin, whereas only a small number belonged to the closed basin of the Ubsunur basin.
The largest river flowing through the region is the Yenisei River, formed at the confluence of the Maly and Bolshoy Yenisei rivers. Most rivers of the Yenisei basin are of mountain origin, that is, they mainly have snow and groundwater feeding. The studied lakes are mostly drainless and mainly fed by groundwater, with the exception of Lake Chagytai [68,69]. We collected the water samples during four hydrological seasons, in autumn (24 October 2021–27 October 2021), winter (7 March 2022–11 March 2022), spring (18 May 2022–22 May 2022) and summer (19 August 2022–23 August 2022). During each season, we sampled 20 water bodies (15 rivers and five lakes), with 140 samples in total. These were selected according to the criterion of proximity to the weather stations throughout the climate/landscape macrotransect of Tyva.

2.2. Analytical Methods

The list of measured parameters included temperature, pH, electrical conductivity (E.C.) and concentration of dissolved gases (CO2 and O2), dissolved organic (DOC) and inorganic carbon (DIC), optical properties of organic matter, as well as CO2 emission flux from the water surface. Dissolved oxygen, pH, electrical conductivity and temperature were measured in situ using an EXO2 multiparameter probe and a WTW Multi 3320 multimeter. The measurement of pCO2 in water was carried out in situ using the GM70 data logger Vaisala®. The pCO2 was measured in situ by an infrared gas analyzer (IRGA, GMT222, Vaisala, Finland) [70]. The sensor was enclosed in a semipermeable membrane and placed directly into the surface water (30–50 cm depth), where it was allowed to equilibrate for approximately 30 min. The sensors were calibrated against standard gas mixtures (0, 800, 3000, 8000 ppm) before and after the sampling. Following calibration, results were corrected using water temperature and barometric pressure.
Carbon dioxide emissions from the water surface were measured by direct floating chamber method using SensAir sensors. We used a freely drifting chamber (30 cm diameter, covered with aluminum tape). The CO2 accumulation rate inside the chamber was recorded continuously at 5 sec intervals for 5–10 min and used to compute (by linear regression if R2 > 0.75) CO2 flux and kCO2 following Kuhn et al., 2018 [71]. For all calculations, the CO2 air–water equilibrium was calculated assuming an air concentration of 400 ppm. Further details of pCO2 and fCO2 measurements in rivers and lakes of adjacent regions are provided elsewhere [1,72,73,74,75].
River and lake water was collected from the surface (depth 0.5 m) in a precleaned polypropylene container with a capacity of 1 L and immediately filtered through nitrate-cellulose filters (<0.45 microns, Sartorius Minisart High Flow). DOC and DIC were measured in the BIO-GEO-CLIM Laboratory (TSU) using a total organic carbon analyzer of the TOC-LCSN series, Shimadzu, with an uncertainty of 2%. As an indicator of the quality of DOM, we measured visual and UV absorbance using spectrophotometry (Agilent Cary 300 spectrophotometer). These measurements allowed for the assessment of spectral slopes (S290–350, S275–295, S350–400); the slope ratio (SR) is useful for approximating the molecular weight of DOM and the degree of its aromaticity. The relative amount of aromatic compounds was assessed via specific UV absorbance (SUVA254). The weight-average molecular weight index (WAMW) of DOM, which reflects the degree of polymerization, was quantified using the absorbance at 280 nm. Finally, the humification index was characterized using the ratio of absorbance at 250 and 400 nm (the E2:E4 ratio).
Total bacterial cell (TBC) concentration was measured after sample fixation in glutaraldehyde (in the field immediately after collection) using flow cytometry (Guava® EasyCyteTM systems, Merck, Rahway, NJ, USA). Cells were stained using 1 µL of 10 times diluted SYBR GREEN solution (10,000×, Merck), added to 250 µL of each sample before analysis. Particles were identified as cells based on green fluorescence and forward scatter.
To build the maps, the DAICHI satellite (ALOS) survey with a resolution of 30 m2 was used together with the 3D Analysis module in the ArcGIS environment. To delineate the catchments of the studied rivers, the Flow Direction module was used in the ArcGIS environment, which calculates flow lines based on data on the heights of nearby points.
Carbon concentrations and fluxes measured in rivers and lakes were tested for normality using a Shapiro–Wilk test. In case the data were not normally distributed, we used nonparametric statistics. Comparisons of dissolved C parameters between seasons and groups of the water bodies were conducted using a nonparametric Mann–Whitney test at a significance level of 0.05. For comparison of unpaired data, a nonparametric H-criterion Kruskal–Wallis test was used to assess the differences between sites. The Pearson rank order correlation coefficient (p < 0.05) was used to determine the relationship between CO2 concentrations and flux and the main hydrochemical parameters of the water column, which are considered as potential drivers (pH, O2, water temperature, electrical conductivity, DOC, DIC and total bacterial number). Boxplots of median, IQR and outlies were constructed for visualization of the differences between different seasons and groups of rivers or lakes.

3. Results

3.1. Major Hydrochemical Parameters

The values of pH varied significantly during different seasons of the year (F = 11.31, p = 0.000, Figure 3a) and between the two major types of the water bodies—lakes and rivers (F = 10.67, p = 0.000; Figure 3b,c). The electrical conductivity (E.C., µS cm−1) in the studied water bodies varied widely and increased in the following order: thermokarst lakes (29 ± 22) < large rivers (172 ± 97) < small rivers (207 ± 121) < freshwater lakes (534 ± 252) < salt lake (54,600 ± 24,380).
The water bodies of the Alash plateau, both rivers and lakes, exhibited the lowest E.C. In winter and autumn during low-water season, the electrical conductivity was significantly higher than in spring and summer (F = 3.3, p = 0.03), see Supplementary Figures S1 and S2. For lakes, a significant relationship (p < 0.05) between electrical conductivity and the pH was established (Figure 4). No such relationship was revealed for rivers. A strong direct relationship between electrical conductivity and DIC content was also observed (Figure 5), which reflected the dominance of bicarbonate ions in the major salt composition of both rivers and lakes (Supplementary Table S1).

3.2. DIC vs. DOC Concentrations

Elevated DIC concentrations were observed in freshwater lakes, especially in Torehol Lake where they reached 120 mg L−1, likely due to the impact of carbonate-rich groundwaters, which are reported to occur within the lake watershed (Figure 6a). The minimal DIC values were recorded in thermokarst lake waters, from 2.8 to 6.6 mg L−1. The seasonal dynamics of DIC demonstrated rather low variations (within 30%) with minimal values observed in spring and maximal ones in winter, see Figure 6b.
The maximal DOC values were measured in the waters of Torehol Lake and the thermokarst lakes. Of all the rivers studied, the Dyttyg-Khem, Durgen and Biche-Bayan-Kol waters were highly enriched in DOM, whereas the majority of the water bodies ranged from 2 to 6 mg L−1 in DOC concentration. As in the case of DIC spatial and temporal patterns (Figure 6), the DOC concentrations demonstrated relative stability across seasons (Figure 7), with an exception of an anomalously high value (14.8 mg L−1) in Torehol Lake during winter. The highest SUVA254 values, which reflect the DOM aromaticity, were observed in the waters of thermokarst lakes, followed by rivers, whereas the minimal values were recorded in freshwater lakes (Figure 8a–c). The highest SUVA254 was recorded in spring and the lowest in winter (Figure 8b).
The ratio E254:E436 is known to indicate the relative role of allochthonous versus autochthonous organic substances in water bodies [76,77]. The DOM composition in the studied water objects was dominated by allochthonous substances, with the exception of the summer and part of the spring periods, when there was an increase in the E254:E436 ratio, indicating an active process of photosynthesis and destruction of detritus in the rivers and lakes themselves, leading to an increase in autochthonous organic matter [78,79].
There was no strong variation in the value of SUVA254 in the waters of both small and large rivers (Figure 8a). However, small river waters had a lower SR value, which indicated an increased degree of aromaticity compared to large rivers (Figure 8d). A high index of SUVA254 (3.8 ± 1 L mg−1 m−1) and a low index of SR (0.9 ± 0.1) in the waters of thermokarst lakes indicated the presence of high-molecular-weight aromatic compounds. In contrast, freshwater lakes, poor in DOC, exhibited much lower SUVA254 and higher SR values compared to thermokarst lakes. The ratio of optical densities E254:E436 in rivers followed the order ‘summer > spring > autumn > winter’ (Figure 9a). This ratio was much higher in thermokarst lakes compared to rivers and freshwater lakes (Figure 9b).

3.3. Spatial and Seasonal Pattern of pCО2 and fCO2 in Rivers and Lakes

The highest values of pCO2 were noted for a group of small rivers, with a maximum in the Adyr-Khem River, the smallest of the studied rivers draining the mound peat bog (Table 3). A high content of dissolved CO2 was also noted in the waters of the Chadan River, originating in the western Tannu-Ola, and in the waters of the Khule River (Torgalyg), whose source is located on the eastern Tannu-Ola. This region (Tannu-Ola) is known for its CO2-rich underground discharges at the earth’s surface. The average values of pCO2 in large rivers differ significantly from small ones; the minimum was marked for the Alash River (Figure 10, Table 3).
In lakes, the dissolved CO2 content was sizably lower than that in rivers; the minimum values were noted for the salt lake Cheder and Tore-Khol Lake, located within a sandy substrate.
In terms of seasonal variations, in both rivers and lakes, the pCO2 was significantly higher in winter due to CO2 accumulation under ice (Figure 11). Note that because shallow thermokarst lakes freeze solid in winter, measurements of the water column during this period in these water bodies were not possible.
Overall, CO2 emissions from the surface of rivers and lakes were quite low (median 0.05–0.15 g C m−2 d−1; Table 4) and did not differ significantly (p > 0.05) among water bodies, including the two main categories—rivers and lakes (Figure 12a,b). There was no significant difference in fCO2 between different seasons of the year (Figure 12c). Exceptionally high CO2 fluxes in water were encountered in the Chadan River and the Ak-Sug River (12 and 5.6 g C m−2 d−1, respectively). Despite their small size, these rivers do not freeze solid during the ice-on period, and they likely possess sizable bicarbonate/CO2-rich underground sources that are discharged at the river bed and hence are mostly pronounced during the winter baseflow.
Analysis of correlations between O2 and CO2 concentrations of all objects across seasons did not allow us to establish any relationship (r = −0.13, p > 0.05); however, when considering solely large rivers and thermokarst lakes, a negative relationship (r = −0.47 and r = −0.95) was observed (Table S1; Figure S4). In freshwater lakes, pCO2 also negatively correlated with O2 level (r = −0.50, p < 0.05). The water temperature exerted a significant (r = 0.64, p < 0.05) impact on pCO2 in large rivers (Table S1).

3.4. Testing Potential Drivers of CO2 Concentrations and Fluxes

During the summer period, the pCO2 of the water column in lakes correlated with the lake water surface area (p < 0.05), whereas during other seasons, this relationship was not significant (p < 0.05); see Figure 13а. The flux of CO2 positively correlated with lake surface area during the spring and winter periods (p < 0.01), as illustrated in Figure 13b. In thermokarst lakes, pCO2 and fCO2 increased with lake surface area in spring but decreased with Sarea in summer and autumn (Figure S5), although these qualitative trends could not be statistically supported due to too low a number of sampled lakes. The river watershed area positively correlated with pCO2 and fCO2 (p < 0.05) in the autumn period, whereas during other seasons, the correlation was absent (p > 0.05), as shown in Figure S6.
Among the ‘internal’ factors likely to control the CO2 pattern (concentration and fluxes) in the water bodies, we tested major hydrochemical parameters of the water column, including pH, conductivity, DIC, DOC, concentration, DOM quality and a microbiological parameters (i.e., total bacterial count), as presented in Table 5. Considering all water objects together, neither pCO2 nor fCO2 correlated significantly (r > 0.39 at p = 0.05) with any tested parameter. Note that this result does not contradict several significant relationships of pCO2 with lake area or O2, which were analyzed for individual seasons. Individual analyses of groups of rivers and lakes demonstrated strong positive control (p < 0.05) of E.C. and pH on fCO2 in large rivers (r = 0.5 and 0.64, respectively) and negative control of pH on pCO2 in lakes (r = −0.99); see Table S1 of the Supplement.
We found that TBC does not correlate with any parameter of the CO2 system in rivers except in weak positive correlations with DOC (RPearson = 0.32, p > 0.05). In lakes, there were weak but not significant (p > 0.05) correlations with DOC (RPearson = 0.31) and fCO2 (RPearson = 0.36, Table 5). This mostly occurred in thermokarst lakes (Table S1). It is therefore possible that some bacterial processing of DOC in the water column (biodegradation) generates excessive CO2 emissions. However, these processes are obscured by other external factors of CO2 delivery to the river or lake, which occur throughout soil and surface waters, respiration of sediments or biodegradation of particulate organic matter.
It is also important to note that wind speed might strongly affect the CO2 emission fluxes, especially in wind-unprotected large lakes. Unfortunately, because we did not use wind-speed-based flux calculations [75] in this study, some environmental controllers listed in Table 5 could have exhibited a lower impact on fCO2 than it occurs under field conditions.

4. Discussion

4.1. Major Solutes, Dissolved Organic and Inorganic Carbon

It is known that inland water systems are sizable sources of CO2, which is a highly important prerequisite for the retroactive link between the climate and C cycle in aquatic systems [80,81,82]. Although at present, the climate changes in the territory of Tyva are already strongly pronounced (e.g., ref. [83]), it is not yet possible to quantify this impact on water bodies. However, the present study can provide a solid background for assessing the current status of C cycle in order to be able to judge the future changes which are likely to occur over the next decades.
The level of DIC concentration was sizably higher than that of DOC, and hence the C export in the studied rivers and C storage in the lakes are dominated by inorganic carbon. Only during spring was the DOC concentration higher than DIC, which could be linked to the snow melt, dilution of the groundwater signal and the dominance of surface flow. In particular, enhanced DOC delivery to lakes and rivers occurs during melt-snow interaction with surface soil and organic-rich litter. This period contrasts with the shallow subsurface flow in summer and autumn, when the water hydrochemical signal is essentially controlled by rock weathering at the watershed. In this regard, the present study corroborates the knowledge of processes controlling river and lake water chemical composition across seasons as established in the boreal, high latitude and subarctic rivers [84,85,86] and subarctic lakes [87]. Note here that the connection between the groundwaters and the lake, especially in permafrost-covered regions such the Altai–Sayan Mountain system, is weaker than that for rivers. A comobilization of DIC and other ions from shallow subsurface/groundwaters, bearing the signal of chemical weathering, for all water objects is confirmed by the observation that the DIC values strongly correlated with the electrical conductivity index. In the Khemchik, Ak-Sug, Alash and Adyr-Khem mountain rivers flowing in the northwestern part of the Republic of Tyva in the Western Sayans, the DIC concentration ranged from 4 to 20 mg L−1 and in other rivers from 9~13 to 55 mg L−1. It is possible that a reason for such reduced DIC values in mountain rivers is an increased precipitation in this area [83], the presence of permafrost [88] and high runoff [89]. An opposite situation is observed in the water bodies of the more arid Ubsu-Nur basin: Tore-Khol Lake and the Tes-Khem River, where the DIC values are the highest. Here, evaporative concentration of major solutes, occurring both within the river watersheds and the lake water column, can be responsible for elevated DIC and S.C. values.
Increased DOC concentrations were observed in lakes compared to rivers, which may be explained by higher water residence time and autochthonous production of OM. This was especially pronounced in thermokarst lakes, which exhibited the highest SUVA (Figure 8b,c) and humification index (E254/E436 ratio) (Figure 9b). The delivery of allochthonous material to the water bodies occurs essentially during summer and spring following the shallow water paths through organic-rich soil layers (Figure 9a). According to the ratio of optical densities E254/E436, allochthonous (peat-originated) organic matter prevails in thermokarst lakes and significantly exceeds the values of the ratios in other water bodies. Low water residence time in thermokarst lakes, their small size and their location within wetland zones may lead to sizable accumulation of terrestrial humic-rich organic matter. Much lower values of the optical index E254/E436 in rivers and freshwater lakes likely correspond to a decrease in the humification index [90] due to biotically driven degradation of humiclike aromatic components in these water objects with high water residence time [91].

4.2. Dissolved C Pattern and CO2 Fluxes: Driving Factors and Comparison with Other Regions

We identified the oxygen regime as one of the main controlling factors of pCO2 level in both rivers and lakes across seasons. The highest CO2 concentrations were observed in O2-impoverished waters, likely due to the important impact of partially anoxic sediment respiration on the CO2 regime of the water column. Similar effects are reported in many boreal and subarctic waters across the world, notably in lakes and channels of the Ob River floodplain (e.g., ref. [73]).
The pCO2 was generally higher in lakes of a smaller surface area, which is a well-known phenomenon in various lakes of the boreal and subarctic zones [1,75]. This can be linked to enhanced terrestrial input of CO2-rich waters from the watershed in small lakes and sizable primary productivity of plankton and periphyton/macrophytes in large, mature lakes, leading to CO2 uptake in the latter, as it is known in thermokarst water bodies in the neighboring territory of permafrost peatlands [92].
The seasonal variations of CO2 concentration were characterized by a maximum pCO2 during winter in rivers and during spring in lakes. The former can be explained by enhanced underground discharge within the river’s main stem and CO2 accumulation under ice, whereas the spring maximum of CO2 in lakes is likely due to a lack of phytoplankton activity and the enhanced input of biologically labile terrestrial DOM from the watershed during freshet. This DOM can be subjected to intensive bacterial processing leading to a maximum of CO2 emissions [93]. In the summer period, lakes act as a CO2 sink due to biological production processes (macrophytes, cyanobacterial bloom) in the water column (e.g., ref. [94]).
The data on carbon (DOC, DIC, CO2) concentrations and CO2 fluxes in the Altai–Sayan mountain system were compared with the data reported for the Qinghai–Tibetan–Plateau [18,95,96], since this region has relatively similar climate and is also impacted by permafrost. To account for thermokarst lakes, we chose the largest compilation of lakes in the sporadic to continuous permafrost zone of Western Siberia, reported in Serikova et al. (2019) [1]. A synthesis of available information is provided in Table 6. The DOC levels of rivers and lakes of the Tuva Republic obtained in this study are quite similar within the range of natural seasonal and spatial variations to those in thermokarst lakes of the QTP as well as the Tibetan rivers. The DIC concentrations in lakes and rivers of the Sayan–Altai Mountain system are also similar to the values of the Qinghai–Tibetan plateau, with the exception of saline lakes, whose C cycle is strongly controlled by local evaporative processes.
In contrast to the DIC and DOC concentrations, for which the results of different regions are fairly comparable, the pCO2 in the Sayan–Altai Mountain region had generally lower values than those observed in the Qinghai–Tibetan Plateau and Western Siberia. It is possible that sizable autochthonous production and CO2 uptake by plankton in lakes and by macrophytes in rivers could be partially responsible for lower CO2 levels in the water bodies of the studied region during the summer period, when the water temperature becomes higher than that in subarctic Western Siberia or high-altitude QTP. At the same time, our results on CO2 flux are comparable to those of Jia et al. (2022) [95] for the lakes and rivers of southwestern and northeastern regions of the QTP (0.3 ± 0.2 (g C m−2 d−1)). However, unlike the maximal flux observed in the present work during winter, the latter authors reported the lowest flux during ice-covered seasons. The lowest CO2 emission is typically observed during cold periods, as also reported by other studies [96,97,98,99,100]. In our case, high wintertime pCO2 and fCO2 could indicate the existence of sizable underground sources of CO2 (discharge of CO2-rich fluids, especially pronounced during low-water-level periods). These effects are visible for large rivers, which demonstrated strong positive impact of E.C., pH and DIC on fCO2 (Table S1).
Overall, the differences of CO2 emissions among different works on northern and mountain aquatic systems may be due to the peculiarities inherent to our region as well as differences in the methodologies used. It is interesting that the carbon uptake flux in the terrestrial Tyva steppe regions is estimated as 184 ± 41 g C m−2 yr−1 [101]. This is fairly comparable with the CO2 emission flux from the rivers of the region, assessed in the present study (213 g C m−2 yr−1). Such a comparison is consistent with the importance of inland water bodies (rivers and lakes) in overall C balance of the terrestrial biomes, as also demonstrated in Western Siberia [102]. However, at the present status of the research, a straightforward comparison of the data on the C pattern revealed in the Sayan–Altai region and those reported in the adjacent subarctic and mountain territories requires more thorough assessment of possible controlling factors such as local climate variation, productivity of the terrestrial compartments, underground influx, respiration of sediments and primary productivity in the river and lakes. In particular, it was beyond the scope of the present study to address the soil type and chemical weathering within each individual catchment of a river or lake. Future research will focus on relating the GIS-based specific watershed parameters (including soils, vegetation, climate and lithology) and CO2 emissions pattern, as it was recently developed for a number of rivers in Siberia (e.g., refs. [2,72,74]).

5. Conclusions

A thorough, first-time assessment of the contemporary status of the C biogeochemical cycle (concentration in and emission from the water surfaces) was conducted in lakes and rivers of the Central Asian Sayan–Altai Mountain system in order to provide a background for judging possible future changes induced by particulate climate instability in this region. Using a physicogeographical and climatic transect of inland water bodies, which comprised a large variety of natural ecosystems and landscapes of the region, we found that permafrost exerts the largest impact on lake water hydrochemistry and C pattern including CO2 exchange with the atmosphere, whereas the size of the river watershed had relatively little impact on the CO2 pattern. In contrast, pCO2 decreased with an increase in lake size, which could be linked to a combination of factors, notably (i) enhanced input of terrestrial, biolabile DOM in small lakes which served as a source of CO2 production in the water column; (ii) terrestrial input of CO2-rich shallow groundwater and soil waters, more pronounced in small lakes; and (iii) CO2 uptake in large lakes due to macrophytes, periphyton and phytoplankton activity. The oxygen regime was found to be an important controlling factor of the pCO2 level in both rivers and lakes during specific seasons, likely due to sediment respiration processes. However, despite these qualitative features of possible CO2 regime control in the studied water bodies, a pairwise correlation analysis did not demonstrate statistically significant relationships between CO2 flux and physicogeographical parameters of river watersheds, river and lake size and internal parameters of the water column, including basic hydrochemical characteristics, DOM concentration and quality and bacterial concentration. A likely reason is the complexity of the ‘external’ and ‘internal’ factors controlling the CO2 exchange between water surfaces and the atmosphere and their frequent counteraction of CO2 production/uptake in the water column, sediments and the watershed.
We argue that further assessment of possible controlling factors such as local climate variation, productivity of the terrestrial compartments, underground influx, respiration of sediments and primary productivity in the rivers and lakes are necessary for comparison of the aquatic C pattern revealed in the Sayan–Altai region in this study to those reported in the adjacent subarctic and mountain territories.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15193411/s1, Figure S1. Box plot of median and IQR range (with outlies as dots) of electrical conductivity in rivers (a); lakes (b); DIC concentrations in rivers (c); and lakes (d), averaged across seasons. Figure S2. Boxplots (median and IQR range) of electrical conductivity during different seasons (lakes and rivers combined together). Figure S3. Plots of DIC concentration (mg L−1) as a function of electrical conductivity in large rivers (a), small rivers (b) and lakes (c). Figure S4. The relationship between pCO2 and oxygen saturation degree of lake and river waters in the Tyva region (all seasons together). Figure S5. The relationship between (a) pCO2 and (b) fCO2 and thermokarst lakes water surface area of the Tyva region. Figure S6. The relationship between pCO2 (a) and fCO2 (b) in river watershed area in the Tyva region (all seasons together). Table S1. Pearson correlation coefficients between major hydrochemical and microbiological parameters of the water column, including pH, O2, S.C., DOC, DIC and optical parameters of DOM reflecting DOM quality. Sarea, S.C., and TBC represent area of the lake surface or river watershed, Electrical Conductivity, and Total Bacterial Count, respectively.

Author Contributions

Conceptualization, A.A.B., L.G.K., A.O.K. and S.K.; methodology, I.I.K., T.V.R., A.S.P., I.V.L., Z.N.K. and S.N.V.; software, I.I.K.; validation, S.N.V. and O.S.P.; formal analysis, A.A.B. and L.G.K.; investigation, A.A.B., L.G.K. and I.I.K.; resources, A.O.K. and S.K.; data curation, A.A.B. and L.G.K.; writing—original draft preparation, A.A.B. and L.G.K.; writing—review and editing, O.S.P.; visualization, A.A.B. and L.G.K.; supervision, S.K.; project administration, A.O.K.; funding acquisition, S.K. and A.O.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Tomsk State University Development Program “Priority 2030” (A.A.B., L.G.K., I.I.K., T.V.R., A.S.P., S.N.V., O.S.P., S.K.) and grant RSCF № 23-14-20015 (A.O.K). Research was performed using the equipment of the unique research installation ‘System of experimental bases located along the latitudinal gradient’ with financial support Ministry of Science and Higher Education of Russia (RF-2296.61321Х0043, 13.UNU.21.0005, contract No. 075-15-2021-672).

Data Availability Statement

Full data set of this study is available at the Mendeley Data Portal: Byzaakay, Arisia; Лариса, Кoлесниченкo; Pokrovsky, Oleg (2023), “Dissolved carbon concentrations and emission fluxes in rivers and lakes of Central Asia (Sayan-Altai mountain region, Tyva)”, Mendeley Data, V1, https://data.mendeley.com/datasets/d35ctw7tb3/1, accessed 23 September 2023.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Detailed description of rivers and lakes of the Tyva Republic sampled in this work.
The Yenisei River is one of the longest and deepest rivers in the world and Russia, flowing into the Kara Sea of the Arctic Ocean. The length is 3487 km, the catchment area is 2,580,000 km2, and the annual flow is 624.41 km3. The flow velocity is from 0.3 to 5 m/s, depending on the location of the river [103]. The research was carried out in the upper reaches of the river. The river originates from the confluence of two sources—the Big Yenisei (Biy-Khem) and the Small Yenisei (Kaa-Khem). The main tributaries in the upper reaches are Elegest, Khemchik, Us, Kantegir. The city of Kyzyl and the Service-Khem kozhuun of the Republic of Tyva, together with the administrative center of Shagonar, is located in the river basin. The non-flooded channel of the Upper Yenisei (Service-Khema) on the territory of Tyva runs through the flat terrain of the Tuva basin with a pronounced steppe landscape, has many channels and islands [104].
The basin of the Bolshoy Yenisei River (Biy-Khem) is located in the wettest, taiga, teeming with lakes and wetlands of Tyva. In the center of it is the vast Todzhinsky basin. The Biy-Khem basin, in the area from the mouth of the Seiba to the mouth of the Uyuk, is characterized by the presence of high mountains with steep slopes and weak forest cover. Below, before the confluence with the Maly Yenisei River (Kaa-Khem), smoothed open landforms and a wide river valley are typical. The length of the river is about 560 km. The width of the riverbed The Yenisei varies from 20-80 m in the upper reaches to 120–290 m in the middle and lower reaches, depths, respectively, from 1–1.5 m to 1.5–4 m. The flow velocity varies from 1.4 to 2.4 m/s. The average long-term water consumption of the Yenisei River in the closing alignment (Kara—Haak village) is 594 m3/s [104].
The Maly Yenisei River (Kaa-Khem) is formed by the confluence of two rivers: Kyzyl-Khem and Balyktyg-Khem. The source of Kyzyl-Khem is located on the territory of Mongolia, whereas Balyktyg-Khem originates from the northern slopes of the Sengilen Highlands. In the upper reaches of the Yenisei River, it is a typical mountain river, its bed abounds with rocky ledges that form numerous rapids and shivers. After entering the Ulugh-Khem basin, it flows in low steppe shores. Compared to Biy-Khem, Kaa-Khem has fewer tributaries. The right tributaries of the river are Unzhey, Honga, Uzhep, Derzyg; the left ones are Shivey, Sizim, Buren. The river is quite full-flowing, fast, the speed is 1.8–2.3 m/s, the width of the channel is from 20 to 80 m. The water discharge is highly variable and depends on the amount of precipitation (Table 1). In May–early June, the spring flood comes, and in July-August, as a result of heavy rains, sometimes there is a summer flood. Groundwater takes part in feeding rivers all year round. The depth of groundwater is about 8-10 m, in floodplains of rivers and streams they come close to the surface [105].
The Tes-Khem River has a mountainous character, belongs to the Ubsu-Nur water system. The length of the river is 139 km, the width of the river is 10–100 m, the catchment area is 4390 km2. The river’s feeding is predominated by rain. In summer, there is a high probability of rain floods. In winter, the river freezes. The river originates on the Tannu-Ola mountain range. The mouth of the river is at an altitude of 1067 m. To the west of Tes-Khem, beyond the sands of Tsugeer-Els, there is the lake Tore-Khol, formed by the damming of the former tributary of Tes-Khem.
The Khemchik River originates on the eastern slope of the Kozer ridge from a peak of 3122 m, belonging to the Shapshalsky ridge system, on the border with the Altai Republic. The entire river is located on the territory of the Republic of Tyva. It flows between two mountain systems—the Western Sayan from the north and the Western Tannu-Ola from the south, collecting all its runoff from them. The length of the river is 320 km, the catchment area is 27 thousand km2. The average water discharge is 102 m3/s [104]. All sources and tributaries are fed by runoff from high ridges, and therefore the river system in the Khemchik basin is widely used for irrigation [106]. The valley is narrow, with steep banks, and there are many boulders in the riverbed. In the Khemchik basin, the river has a flat character with bends. Large stone remains are not uncommon. The flow velocity is low. The river is essentially fed from underground discharge. The Khemchik’s flood occurs in summer, from June to August inclusive (with a maximum in July), and is the result of summer precipitation in the form of rain. In September, water consumption decreases significantly, and this decrease continues until the beginning of winter, along with a decrease in precipitation [107].
The Alash River, a left tributary of the Khemchik River, is formed by the confluence of the Chulcha and Kara-Khol rivers; its length is 125 km and the slope is up to 3.5 m/km. The catchment area is 4630 km. The width of the riverbed is up to 60–90 m. The valley of the river is wide. The river acquires a mountain-steppe character. Absolute elevation marks in the area of the river flow range from 1000 to 1200 m [108].
The Ak-Sug River (in the upper reaches of the Ak-Khem), a left tributary of the Khemchik River, originates in one of the lakes of the Dashtyg-Khem ridge at an altitude of 2115 m above sea level. The length of the river is 160 km, the area of its catchment area is 3170 km2. The average annual water discharge is 14 m3/s. Agriculture with irrigation is developed in the valley. Annual water consumption is 1% of the annual flow of the river (15 million m3). The river has 46 tributaries less than 10 km, the total length of which is 173 km, there are also 82 lakes in the catchment area, the total area of which is 10.23 km2. The Kyzyl-Taiga Mountain is located at the headwaters of the Ak-Sug River [109]. The absolute heights of the terrain through which the Ak-Sug River flows vary significantly from 1320 to 2100 m. The width of the riverbed ranges from 2–5 m in winter, to 15–20 m in summer, the depth, respectively, from 0.2–0.3 m to 0.7–1.1 m. The river valley on the investigated section of the Ak-Sug river is forested, the lower parts of the slopes are relatively flat, the upper ones with rocky and shrubby tundra are steep and rocky. At the beginning of the lower third of this segment on the left bank of the Ak-Suga at absolute altitudes of 1350–1500 m there is an eponymous deposit of molybdenum-copper sulfide ores. The deposit was exposed by erosion and subjected to the movement of part of the ore material by a glacier [110].
The Chadan River is a right tributary of the Khemchik River. The source is located on the western ridge of Tannu-Ola. The length of the watercourse is 98 km, the catchment area is 2200 km2. The absolute height of the terrain is 800 m [111].
The Durgen River is a left tributary of the The Yenisei of the 3rd order: flows into the Mezhegey River, a tributary of the Elegest River. The source is located in the highlands of the East Tannu-Ola ridge, then the river flows along its northern macro slope. According to the State Water Register, the total length of the river is 93 km. The speed of the river flow changes in accordance with the change in terrain. The river is distinguished by a rare cascade of waterfalls. The river is located within the Mezhegeysky coal deposit (Khayan et al., 2014) [112]. In the middle course, the river flows through the territory of two villages in which there is no sewerage, the river is used for watering livestock. Part of the catchment area is occupied by a specially protected natural area—the Durgensky Nature Reserve.
The Chaa-Khol River, a right tributary of the The Yenisei River, flows through the territory of the Ulug-Khem and Chaa-Khol kozhuuns of Tyva. The total length of the river is 90 km. The catchment area has 1730 km2. It begins on the northern slope of the Western Tannu-Ola ridge between the Tyndy-Ula and Hule-Bozh mountains. It flows in a northerly direction through a valley overgrown with larch forest in the upper reaches and across the steppe in the lower reaches. In the middle course it is divided into several channels. It flows into the Yenisei River at a distance of 3315 km from the mouth. The height of the mouth is 540 m above sea level [103].
The Khule River (Torgalyg), a right tributary of the Shagonar River, originates on the northern macroscline of the Eastern Tannu-Ola, flows through the Central Tuva Basin. The length of the river is 40 km, the width varies from 0.5 to 1.5 m, the catchment area is 610 km². Altitude above sea level: 535 m [113].
Anyyak-Chyrgaki River, a tributary of the Chirgaki River, a right tributary of the Khemchik River. The length of the watercourse is 52 km, the catchment area is 410 km2. The bottom of the river is pebbly. During the study period, the water is transparent to the bottom, its color is bluish-green, the depth of the river ranges from 0.2 to 2.0 m. The absolute elevation of the terrain is 800 m. The source of the river begins on the north-western forested slopes of the Western Tannu-Ola [114].
The Biche-Bayan-Gol River is a right tributary of the Yenisei, into which it flows 12.5 km from the mouth, 14 km northwest of the city of Kyzyl. The length of the watercourse is 32 km [115].
The Adyr-Khem River flows through the territory of the Barun-Khemchik district on the Alash plateau. The width of the Adyr-Khem River is 1.5 m, the depth is from 0.25–1m, the flow velocity is about 0.25 m/s [116].
Lake Tore-Khol is the only large freshwater lake in the Ubsunur basin. The lake has no tributaries, receiving water from springs of sand dunes, which are located in a horseshoe-shaped bowl of the southern shore and drain into the lake in the form of a stream. The Torehole is located in a shallow depression with flat sides, the average depth is 7 m, the area is 42 km2. The greatest depth (up to 40 m) was noted at the isthmus in the southern part. The absolute height of the terrain is 1149 m. The lake is oligotrophic, the least productive of all the lakes of the Ubsunur basin [117]. It contains a relatively small percentage of organogens and bacteria. Therefore, the lake is transparent, clean, rich in underutilized oxygen by organisms. In the summer, the surface layers warm up well, in the middle of August in the afternoon the temperature rises to 21 °C at an air temperature of 32 °C [118].
Lake Chagytai is the largest freshwater lake of the Tuva basin, located in its southern part on the border of the settled hilly foothills of the Vostochny Tannu-Ola ridge, at an altitude of 1003 meters above sea level. This is the largest freshwater lake in the Tuva basin. It has the shape of an almost regular circle with a diameter of about 6 km, the length of the coastline is about 20 km. The area of the lake is 2860 hectares, its maximum depth reaches 17 meters [103]. In the west and southeast of the reservoir, the slopes of the mountains are covered with taiga. The bottom of the lake is sandy and pebbly. The shores are mostly flat, sometimes rocky, sometimes sandy. The south-eastern shore is swampy, overgrown with talus, birches, larches. The only river flowing from the lake, the Mazhalyk River, originates here.
Lake Cheder is located 45 km south of the city of Kyzyl, in the south of the Tuva basin in a drainless depression, on board of which the sandy-clay rocks of the Jurassic come out. The depression is surrounded by a hilly, treeless plain. The lake is drainless, located in a shallow depression, within the vast Central Tuva basin, at an absolute mark of 706 m. The lake is fed by a stream flowing into the lake in the southern part, as well as groundwater from quaternary lake sediments. The lake has a slightly elongated shape. The water surface area is 430 ha, length—4.5 km, width—0.8–1.5 km, depth—1.5–1.8 m [119]. The water in the lake is highly mineralized. The main component is sodium sulfate. The mineralization of water in the lake, depending on the level, ranges from 80 in April and May to 200 grams per liter in August. Potassium salts (0.250 g/L), fluorine (0.003 g/L), iodine (0.001 g/L), strontium (0.001 g/L) are present in small quantities. The shores and bottom of the lake are composed of silt mud with a thickness of up to 2 meters. The mud is gray or gray-black with the smell of hydrogen sulfide. On the territory of the lake there is 46-m deep well with mineral drinking water. The mineral water captured by the well is cold, low-mineralized chloride magnesium-sodium-calcium in composition and slightly alkaline in the nature of the reaction of the medium [120].
Thermokarst lakes are located in the Alash Highlands in the Barun-Khemchik district, 6.5 km from the border with Khakassia and 150 km from the city of Ak-Dovurak. Absolute marks in the area of lakes range from 1800 to 2000 m. Absolute heights in the lakes range from 1840 to 1860 m. There are about 200 thermokarst lakes in the study area. Basically, the lakes are small, most of the lakes have a rounded, oval-elongated shape, some of the shores are overgrown and swampy. Two round-shaped lakes were selected for the study.

References

  1. Serikova, S.; Pokrovsky, O.; Laudon, H.; Kritskov, I.; Lim, A.G.; Manasypov, R.M.; Karlsson, J. High carbon emissions from thermokarst lakes of Western Siberia. Nat. Commun. 2019, 10, 1552. [Google Scholar] [CrossRef] [PubMed]
  2. Krickov, I.; Lim, A.; Shirokova, L.; Korets, M.; Karlsson, J.; Pokrovsky, O. Environmental controllers for carbon emission and concentration patterns in Siberian rivers during different seasons. Sci. Total Environ. 2023, 859, 160202. [Google Scholar] [CrossRef] [PubMed]
  3. Kirpotin, S.; Antoshkina, O.; Berezin, A.; Elshehawi, S.; Feurdean, A.; Lapshina, E.D.; Pokrovsky, O.S.; Peregon, A.M.; Semenova, N.M.; Tanneberger, F.; et al. Great Vasyugan Mire: How the world’s largest peatland helps addressing the world’s largest problems. Ambio 2021, 50, 2038–2049. [Google Scholar] [CrossRef] [PubMed]
  4. Chupakov, A.; Pokrovsky, O.; Moreva, O.; Kotova, E.; Vorobyeva, T.Y.; Shirokova, L.S. Export of organic carbon, nutrients and metals by the mid-sized Pechora River to the Arctic Ocean. Chem. Geol. 2023, 632, 121524. [Google Scholar] [CrossRef]
  5. Stackpoole, S.; Butman, D.; Clow, D.; Verdin, K.; Gaglioti, B.V.; Genet, H.; Striegl, R.G. Inland waters and their role in the carbon cycle of Alaska. Ecol. Appl. 2017, 27, 1403–1420. [Google Scholar] [CrossRef]
  6. Hutchins, R.; Casas-Ruiz, J.; Prairie, Y.; Del Giorgio, P. Magnitude and drivers of integrated fluvial network greenhouse gas emissions across the boreal landscape in Québec. Water Res. 2020, 173, 115556. [Google Scholar] [CrossRef]
  7. Pokrovsky, O.; Shirokova, L.; Manasypov, R.; Kirpotin, S.; Kulizhsky, S.; Kolesnichenko, L.; Loiko, S.; Vorobyev, S. Thermokarst lakes of Western Siberia: A complex biogeochemical multidisciplinary approach. Int. J. Environ. Stud. 2014, 71, 733–748. [Google Scholar] [CrossRef]
  8. Dristi, A.; Xu, Y.J. Large Uncertainties in CO2 Water–Air Outgassing Estimation with Gas Exchange Coefficient KT for a Large Lowland River. Water 2023, 15, 2621. [Google Scholar] [CrossRef]
  9. Huang, Y.; Zhang, L.; Ran, L. Total Organic Carbon Concentration and Export in a Human-Dominated Urban River: A Case Study in the Shenzhen River and Bay Basin. Water 2022, 14, 2102. [Google Scholar] [CrossRef]
  10. Wang, C.; Xv, Y.; Li, S.; Li, X. Interconnected River–Lake Project Decreased CO2 and CH4 Emission from Urban Rivers. Water 2023, 15, 1986. [Google Scholar] [CrossRef]
  11. Hastie, A.; Lauerwald, R.; Ciais, P.; Regnier, P. Aquatic carbon fluxes dampen the overall variation of net ecosystem productivity in the Amazon basin: An analysis of the interannual variability in the boundless carbon cycle. Glob. Chang. Biol. 2019, 25, 2094–2111. [Google Scholar] [CrossRef] [PubMed]
  12. Moyer, R.; Powell, C.; Gordon, D.; Long, J.; Bliss, C. Abundance, distribution, and fluxes of dissolved organic carbon (DOC) in four small sub-tropical rivers of the Tampa Bay Estuary (Florida, USA). Appl. Geochem. 2015, 63, 550–562. [Google Scholar] [CrossRef]
  13. Soria Reinoso, I.; Alcocer, J.; Sánchez-Carrillo, S.; Garcia-Oliva, F.; Cuevas-Lara, D.; Cortes-Guzman, D.; Oseguera, L.A. The Seasonal Dynamics of Organic and Inorganic Carbon along the Tropical Usumacinta River Basin Mexico. Water 2022, 14, 2703. [Google Scholar] [CrossRef]
  14. Huang, T.H.; Fu, Y.H.; Pan, P.Y.; Chen, C.T.A. Fluvial carbon fluxes in tropical rivers. Curr. Opin. Environ. Sustain. 2012, 4, 162–169. [Google Scholar] [CrossRef]
  15. Galy, A.; Yang, Y.; Fang, X. Sequestration of carbon as carbonate in the critical zone: Insights from the Himalayas and Tibetan Plateau. Acta Geochim. 2017, 36, 389–391. [Google Scholar] [CrossRef]
  16. Emberson, R.; Galy, A.; Hovius, N. Weathering of Reactive Mineral Phases in Landslides Acts as a Source of Carbon Dioxide in Mountain Belts. J. Geophys. Res. Earth Surf. 2018, 123, 2695–2713. [Google Scholar] [CrossRef]
  17. Ruan, X.; Galy, A. On the significance of periglacial conditions in active mountain belts for chemical weathering processes: Insights from the Chayu area, SE Tibet. Chem. Geol. 2021, 585, 120581. [Google Scholar] [CrossRef]
  18. Song, C.; Wang, G.; Hu, Z.; Zhang, T.; Huang, K.; Chen, X.; Li, Y. Net ecosystem carbon budget of a grassland ecosystem in central Qinghai-Tibet Plateau: Integrating terrestrial and aquatic carbon fluxes at catchment scale. Agric. For. Meteorol. 2020, 290, 108021. [Google Scholar] [CrossRef]
  19. Song, C.; Wang, G.; Mao, T.; Huang, K.; Sun, X.; Hu, Z.; Chang, R.; Chen, X.; Raymond, P.A. Spatiotemporal variability and sources of DIC in permafrost catchments of the Yangtze River source region: Insights from stable carbon isotope and water chemistry. Water Resour. Res. 2020, 55, e2019WR025343. [Google Scholar] [CrossRef]
  20. Wang, L.; Du, Z.; Wei, Z.; Xu, Q.; Feng, Y.; Lin, P.; Lin, J.; Chen, S.; Qiao, Y.; Shi, J.; et al. High methane emissions from thermokarst lakes on the Tibetan Plateau are largely attributed to ebullition fluxes. Sci. Total Environ. 2021, 801, 149692. [Google Scholar] [CrossRef]
  21. Zhang, H.; Wang, F.; Shan, S.; Ren, P.; Luo, C.; Fu, W.; Sun, S.; Wang, X. Sources and cycling of dissolved organic and inorganic carbon on the northern Qinghai-Tibetan Plateau: Radiocarbon results from Qinghai Lake. Sci. Total Environ. 2022, 851, 158123. [Google Scholar] [CrossRef] [PubMed]
  22. Zhong, J.; Li, S.L.; Zhu, X.; Jing, L.; Liu, J.; Sen, X.; Sheng, X.; Cong-Qiang, L. Dynamics and fluxes of dissolved carbon under short-term climate variabilities in headwaters of the Changjiang River, draining the Qinghai-Tibet Plateau. J. Hydrol. 2021, 596, 126128. [Google Scholar] [CrossRef]
  23. Lin, P.; Du, Z.; Wang, L.; Liu, J.; Xu, Q.; Du, J.; Jiang, R. Hotspots of riverine greenhouse gas (CH4, CO2, N2O) emissions from Qinghai Lake Basin on the northeast Tibetan Plateau. Sci. Total Environ. 2023, 857 Pt 1, 159373. [Google Scholar] [CrossRef] [PubMed]
  24. Kokorin, A. (Ed.) Assessment Report: Climate Change and Its Impact on Ecosystems, Population and Economy of the Russian Portion of the Altai-Sayan Ecoregion; WWF Russia: Moscow, Russia, 2011; p. 152. [Google Scholar]
  25. Golubeva, E.; Kerzhentsev, A.; Koropachinsky, I.; Kurbatskaya, S.S.; Pivovarova, Z.F.; Sedelnikov, V.P.; Surov, A.V.; Tebleeva, U.T.; Khakimov, F.I.; Khudyakov, Y.S. Experiment “Ubsu-Nur”; Part 1; Intellect: Moscow, Russia, 1995; Volume 1, 336p. (In Russian) [Google Scholar]
  26. Anandhi, A.; Kannan, N. Vulnerability assessment of water resources—Translating a theoretical concept to an operational framework using systems thinking approach in a changing climate: Case study in Ogallala Aquifer. J. Hydrol. 2018, 557, 460–474. [Google Scholar] [CrossRef]
  27. Joseph, J.; Ghosh, S.; Pathak, A.; Sahai, A. Hydrologic impacts of climate change: Comparisons between hydrological parameter uncertainty and climate model uncertainty. J. Hydrol. 2018, 566, 1–22. [Google Scholar] [CrossRef]
  28. Chaturvedi, A.; Pandey, B.; Yadav, A.; Saroj, S. Chapter 5—An overview of the potential impacts of global climate change on water resources. In Water Conservation in the Era of Global Climate; Elsevier: Amsterdam, The Netherlands, 2021; pp. 99–120. [Google Scholar] [CrossRef]
  29. Chen, Y.; Li, B.; Li, Z.; Li, W. Water resource formation and conversion and water security in arid region of Northwest China. J. Geogr. Sci 2016, 26, 939–952. [Google Scholar] [CrossRef]
  30. Luo, M.; Liu, T.; Meng, F.; Duan, Y.; Bao, A.; Xing, W.; Feng, X.; Maeyer, P.; Frankl, A. Identifying climate change impacts on water resources in Xinjiang, China. Sci. Total Environ. 2019, 676, 613–626. [Google Scholar] [CrossRef]
  31. Chen, Z.; Chen, Y.; Li, B. Quantifying the effects of climate variability and human activities on runoff for Kaidu River Basin in arid region of northwest China. Theor. Appl. Clim. 2013, 111, 537–545. [Google Scholar] [CrossRef]
  32. Kirpotin, S.; Callaghan, T.V.; Peregon, A.; Babenko, A.S.; Berman, D.I.; Bulakhova, N.A.; Byzaakay, A.A.; Chernykh, M.T.; Chursin, V.; Interesova, E.A.; et al. Impacts of environmental change on biodiversity and vegetation dynamics in Siberia. Ambio 2021, 50, 1926–1952. [Google Scholar] [CrossRef]
  33. Callaghan, T.; Shaduyko, O.; Kirpotin, S.; Gordov, E. Siberian environmental change: Synthesis of recent studies and opportunities for networking. Ambio 2021, 50, 2104–2127. [Google Scholar] [CrossRef]
  34. Shukla, P.R.; Skea, J.; Buendía, E.; Masson-Delmotte, V.; Pörtner, H.; Roberts, D.C.; Zhai, P.; Slade, R.; Connors, S.C.; Diemen, S.V.; et al. IPCC, 2019: Climate Change and Land: An IPCC Special Report on Climate Change, Desertification, Land Degradation, Sustainable Land Management, Food Security, and Greenhouse Gas Fluxes in Terrestrial Ecosystems; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2019; p. 896. [Google Scholar] [CrossRef]
  35. Shukla, P.R.; Skea, J.; Reisinger, A.; Slade, R.; Fradera, R.; Pathak, M.; Khourdajie, A.A.; Belkacemi, M.; van Diemen, R.; Hasija, A.; et al. IPCC, 2022: Summary for Policymakers. In Climate Change 2022: Mitigation of Climate Change. Contribution of Working Group III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2023. [Google Scholar] [CrossRef]
  36. Zalibekov, Z.G. The arid regions of the world and their dynamics in conditions of modern climatic warming. Arid Ecosyst. 2011, 1, 1–7. [Google Scholar] [CrossRef]
  37. Kulikov, A.I.; Ubugunov, L.L.; Mangataev, A.T. Global climate change and its impact on ecosystems. Arid Ecosyst. 2014, 4, 135–141. [Google Scholar] [CrossRef]
  38. Berdugo, M.; Delgado-Baquerizo, M.; Soliveres, S.; Hernández-Clemente, R.; Zhao, Y.; Gaitán, J.J.; Gross, N.; Saiz, H.; Maire, V.; Lehmann, A.; et al. Global ecosystem thresholds driven by aridity. Science 2020, 367, 787–790. [Google Scholar] [CrossRef] [PubMed]
  39. Yao, Y.; Luo, D.; Dai, A.; Simmonds, I. Increased quasi stationarity and persistence of winter Ural blocking and Eurasian extreme cold events in response to Arctic warming. Part I: Insights from observational analyses. J. Clim. 2017, 30, 3549–3568. [Google Scholar] [CrossRef]
  40. Mokhov, I.I.; Akperov, M.G.; Prokofyeva, M.A.; Timazhev, A.V.; Lupo, A.R.; Treut, H.L. Blockings in the Northern hemisphere and Euro-Atlantic region: Estimates of changes from reanalysis data and model simulations. Dokl. Earth Sc. 2013, 449, 430–433. [Google Scholar] [CrossRef]
  41. Watanabe, T.; Matsuyama, H.; Kuzhevskaia, I.; Nechepurenko, O.; Chursin, V.; Zemtsov, V. Long-Term Trends of Extreme Climate Indexes in the Southern Part of Siberia in Comparison with Those of Surrounding Regions. Atmosphere 2023, 14, 1131. [Google Scholar] [CrossRef]
  42. Kattsov, V.M. (Ed.) The Third Assessment Report on Climate Change and Its Consequences on the Territory of the Russian Federation; Science-Intensive Technologies: St. Petersburg, Russia, 2022; 676p. (In Russian) [Google Scholar]
  43. Kuular, H.B.O.; Chuldum, A.F.; Hertek, S.B.O.; Namzin, S.A. Spatio-temporal distribution of the spring burnability of the vegetation cover of the Republic of Tyva in 2000–2020. Nat. Resour. Environ. Soc. 2021, 2, 18–22. (In Russian) [Google Scholar] [CrossRef]
  44. Climate Risk Country Profile: Mongolia (2021): The World Bank Group and the Asian. Available online: https://www.adb.org/sites/default/files/publication/709901/climate-risk-country-profile-mongolia.pdf (accessed on 23 July 2023).
  45. Adrian, R.; O’Reilly, C.M.; Zagarese, H.; Baines, S.B.; Hessen, D.O.; Keller, W.; Livingstone, D.M.; Sommaruga, R.; Straile, D.; Van Donk, E.; et al. Lakes as sentinels of climate change. Limnol. Oceanogr. 2009, 54, 2283–2297. [Google Scholar] [CrossRef]
  46. Mladenov, N.; Sommaruga, R.; Morales-Baquero, R.; Laurion, I.; Camarero, L.; Diéguez, M.C.; Camacho, A.; Delgado, A.; Torres, O.; Chen, Z.; et al. Dust inputs and bacteria influence dissolved organic matter in clear alpine lakes. Nat. Commun. 2011, 2, 405. [Google Scholar] [CrossRef]
  47. Moser, K.A.; Baron, J.S.; Brahney, J.; Oleksy, I.A.; Saros, J.E.; Hundey, E.J.; Sadro, S.A.; Kopáček, J.; Sommaruga, R.; Kainz, M.J.; et al. Mountain lakes: Eyes on global environmental change. Glob. Planet. Chang. 2019, 178, 77–95. [Google Scholar] [CrossRef]
  48. Rogora, M.; Frate, L.; Carranza, M.L.; Freppaz, M.; Stanisci, A.; Bertani, I.; Bottarin, R.; Brambilla, A.; Canullo, R.; Carbognani, M.; et al. Assessment of climate change effects on mountain ecosystems through a cross-site analysis in the Alps and Apennines. Sci. Total Environ. 2018, 624, 1429–1442. [Google Scholar] [CrossRef] [PubMed]
  49. Tranvik, L.J.; Downing, J.A.; Cotner, J.B.; Loiselle, S.A.; Striegl, R.G.; Ballatore, T.J.; Dillon, P.; Finlay, K.; Fortino, K.; Knoll, L.B.; et al. Lakes and Reservoirs as Regulators of Carbon Cycling and Climate. Limnol. Oceanogr. 2009, 54, 2298–2314. [Google Scholar] [CrossRef]
  50. Lehner, B.; Doll, P. Development and validation of a global database of lakes, reservoirs and wetlands. J. Hydrol. 2004, 296, 1–22. [Google Scholar] [CrossRef]
  51. MacDonald, G.; Edwards, T.; Moser, K.; Pienitz, R.; Smol, J. Rapid response of treeline vegetation and lakes to past climate warming. Nature 1993, 361, 243–246. [Google Scholar] [CrossRef]
  52. Bunting, L.; Leavitt, P.R.; Weidman, R.P.; Vinebrooke, R.D. Regulation of the Nitrogen Biogeochemistry of Mountain Lakes by Subsidies of Terrestrial Dissolved Organic Matter and the Implications for Climate Studies. Limnol. Oceanogr. 2010, 55, 333–345. Available online: http://www.jstor.org/stable/20622883 (accessed on 23 September 2023). [CrossRef]
  53. Oidup, C.K. Lithium-uranium mineralization of salt lakes and underground springs of Central Tuva. Geospheric. Res. 2018, 3, 22–33. [Google Scholar] [CrossRef]
  54. Kyzyl, O.M. Hydrogeochemistry of mineral water sources of the Republic of Tyva. Bull. Kurgan State Univ. 2015, 4, 2222–3371. Available online: https://cyberleninka.ru/article/n/gidrogeohimiya-mineralnyh-istochnikov-vod-respubliki-tyva (accessed on 18 March 2023). (In Russian).
  55. Yu, I.K.; Khvashchevskaya, A.A. Features of the chemical composition of groundwater springs of Western Tuva. Resort Base Nat. Health-Impr. Areas Tuva Adjac. Reg. 2015, 2, 167–173. (In Russian) [Google Scholar]
  56. Kalnaya, O.I.; Rychkova, O.I.; Ayunova, O.D.; Archimaeva, T.P. Ecological and hydrochemical state of watercourses in the area of Kyzyl (Republic of Tyva). Nat. Resour. Environ. Soc. 2020, 4, 52–60. (In Russian) [Google Scholar]
  57. Kalnaya, O.I.; Ainova, O.D. Retrospective analysis of the ecological and hydrochemical state of two-dimensional and auxiliary means in the region of Lake Duskhol (Svatikovo), central Tuva. Nat. Resour. Environ. Soc. 2022, 1, 37–47. (In Russian) [Google Scholar]
  58. Kirova, N.A.; Kalnaya, O.I.; Ayunova, O.D. On the Issue of Hydrochemistry and Biology of Lake Duskhol (Tuva). Izv. AO RGO 2018, 82–88. Available online: https://cyberleninka.ru/article/n/k-voprosu-o-gidrohimii-i-biologii-ozera-dus-hol-tuva. (accessed on 20 March 2023). (In Russian).
  59. Sambuu, A.; Ayunova, O.; Chupikova, S.; Dugerzhaa, C.S.; Biche-ool, O.E. Water resources of the republic of Tuva and their current state. Int. J. Hydrol. Sci. Technol. 2020, 7, 60–66. (In Russian) [Google Scholar] [CrossRef]
  60. State Report “On the State and Environmental Protection of the Republic of Tyva in 2019”—Text: Electronic//Ministry of Natural Resources and Ecology of the Republic of Tyva: Official Website. 2019. Available online: https://mpr.rtyva.ru/upload/files/f992a01d-bc33-488f-8f71-ea266cfa02ec.pdf (accessed on 14 February 2023). (In Russian).
  61. Chupikova, S.A.; Prudnikov, S.D.; Chuldum, A.F. Morphometric analysis of the water intake of relig-Hem (Tuva) using GIS and DDZ. Bull. SGUGiT (Sib. State Geosystem Univ. Technol. Univ.) 2023, 28, 76–88. (In Russian) [Google Scholar]
  62. Mineeva, L.A.; Arakchaa, K.K.D.; Kyzyl, O.M. Physical and chemical characteristics of mineral waters from the Shumak and Choigan deposits. Bull. Irkutsk. State Univ. Ser. Earth Sci. 2016, 17, 115–134. (In Russian) [Google Scholar]
  63. Kopylova, Y.G.; Guseva, N.V.; Arakchaa, K.D.; Khvashchevskaya, A.A. Geochemistry of carbon dioxide mineral waters of the Choigan natural complex (northeastern Tuva). Russ. Geol. Geophys. 2014, 55, 1295–1305. [Google Scholar] [CrossRef]
  64. Sodnam, N.I.; Kashkak, E.S.; Oorzhak, U.S. Physico-chemical review of arzhanov Ovyursky district. Bull. Tuva State Univ. Nat. Agric. Sci. 2019, 2, 52–58. (In Russian) [Google Scholar]
  65. Klopotova, N.G.; Sidorina, N.G. The current state and study of therapeutic hydromineral resources of Tuva. Resort Base Nat. Health-Impr. Areas Tuva Adjac. Reg. 2013, 1, 16–19. (In Russian) [Google Scholar]
  66. Szopińska, M.; Szumińska, D.; Polkowska, Ż.; Machowiak, K.; Lehmann, S.; Chmiel, S. The chemistry of river–lake systems in the context of permafrost occurrence (Mongolia, Valley of the Lakes) Part, I. Analysis of ion and trace metal concentrations. Sediment. Geol. 2016, 340, 74–83. [Google Scholar] [CrossRef]
  67. Kuular, H.B. Climate warming in the Republic of Tyva according to ground observations. Nat. Resour. Environ. Soc. 2021, 1, 62–67. (In Russian) [Google Scholar] [CrossRef]
  68. Arakchaa, K.D. Prospects for the development of balneoresources of lake Cheder. Nat. Resour. Environ. Soc. 2021, 1, 21–29. Available online: https://cyberleninka.ru/article/n/perspektivy-osvoeniya-balneoresursov-ozera-cheder (accessed on 23 September 2023). (In Russian).
  69. Dubovik, D.S.; Yakutin, M.V. The use of Earth remote sensing data in monitoring of small lakes of the Ubsunur basin (Tyva). Interexpo Geo-Sib. 2016, 4, 74–78. [Google Scholar]
  70. Johnson, M.S.; Billett, M.F.; Dinsmore, K.J.; Wallin, M.B.; Dyson, K.E.; Jassal, R.S. Direct and continuous measurement of dissolved carbon dioxide in freshwater aquatic systems—Method and applications. Ecohydrology 2010, 3, 68–78. [Google Scholar] [CrossRef]
  71. Kuhn, M.; Lundin, E.J.; Giesler, R.; Johansson, M.; Karlsson, J. Emissions from thaw ponds largely offset the carbon sink of northern permafrost wetlands. Sci. Rep. 2018, 8, 9535. [Google Scholar] [CrossRef] [PubMed]
  72. Lim, A.; Kritskov, I.; Vorobyev, S.; Korets, M.; Kopysov, S.; Shirokova, L.S.; Karlsson, J.; Pokrovsky, O.S. Carbon emission and export from the Ket River, western Siberia. Biogeosciences 2022, 19, 5859–5877. [Google Scholar] [CrossRef]
  73. Krickov, I.; Serikova, S.; Pokrovsky, O.; Vorobyev, S.N.; Lim, A.G.; Siewert, M.B.; Karlsson, J. Sizable carbon emission from the floodplain of Ob River. Ecol. Indic. 2021, 131, 108164. [Google Scholar] [CrossRef]
  74. Vorobyev, S.N.; Karlsson, J.; Kolesnichenko, Y.Y.; Korets, M.A.; Pokrovsky, O.S. Fluvial carbon dioxide emission from the Lena River basin during the spring flood. Biogeosciences 2021, 18, 4919–4936. [Google Scholar] [CrossRef]
  75. Zabelina, S.; Shirokova, L.; Klimov, S.; Chupakov, A.; Lim, A.; Polishchuk, Y.M.; Polishchuk, V.Y.; Bogdanov, A.N.; Muratov, I.N.; Guerin, F.; et al. Carbon emission from thermokarst lakes in NE European tundra. Limnol. Oceanogr. 2020, 66, S216–S230. [Google Scholar] [CrossRef]
  76. Drozdova, O.Y.; Ilina, S.M.; Anokhina, N.A.; Zavgorodnyaya, Y.A.; Demin, V.V.; Lapitskiy, S.A. Organic Matter Transformation in the Conjugate Series of Surface Water in Northern Karelia. Water Resour. 2019, 46, 52–58. [Google Scholar] [CrossRef]
  77. Ilina, S.M.; Drozdova, O.Y.; Lapitskiy, S.A.; Alekhin, Y.V.; Demin, V.; Zavgorodnyaya, Y.; Shirokova, L.; Viers, J.; Pokrovsky, O.S. Size fractionation and optical properties of dissolved organic matter in the continuum soil solution-bog-river and terminal lake of a boreal watershed. Org. Geochem. 2014, 66, 14–24. [Google Scholar] [CrossRef]
  78. Jaffé, R.; Boyer, J.N.; Lu, X.; Maie, N.; Yang, C.; Scully, N.M.; Mock, S. Source characterization of dissolved organic matter in a subtropical mangrove-dominated estuary by fluorescence analysis. Mar. Chem. 2004, 84, 195–210. [Google Scholar] [CrossRef]
  79. Battin, T.J. Dissolved organic materials and its optical properties in a blackwater tributary of the upper Orinoco River, Venezuela. Org. Geochem. 1998, 28, 561–569. [Google Scholar] [CrossRef]
  80. Li, M.; Peng, C.; Wang, M.; Xue, W.; Zhang, K.; Wang, K.; Shi, G.; Zhu, Q. The carbon flux of global rivers: A re-evaluation of amount and spatial patterns. Ecol. Indic. 2017, 80, 40–51. [Google Scholar] [CrossRef]
  81. Cole, J.J.; Prairie, Y.T.; Caraco, N.F.; McDowell, W.H.; Tranvik, L.J.; Striegl, R.G.; Duarte, C.M.; Kortelainen, P.; Downing, J.A.; Middelburg, J.J.; et al. Plumbing the Global Carbon Cycle: Integrating Inland Waters into the Terrestrial Carbon Budget. Ecosystems 2007, 10, 172–185. [Google Scholar] [CrossRef]
  82. Raymond, P.; Hartmann, J.; Lauerwald, R.; Lauerwald, R.; Sobek, S.; McDonald, C.; Hoover, M.; Butman, D.; Striegl, R.; Mayorga, E.; et al. Global carbon dioxide emissions from inland waters. Nature 2013, 503, 355–359. [Google Scholar] [CrossRef] [PubMed]
  83. Kirpotin, S.N.; Kvasnikova, Z.N.; Potapova, S.A.; Volkova, I.I.; Volkov, I.V.; Pyak, A.I.; Byzaakay, A.A.; Kolesnichenko, L.G.; Lushchaeva, I.V.; Khovalyg, A.O.; et al. Pilot Studies of the Unique Highland Palsa Mire in Western Sayan (Tuva Republic, Russian Federation). Atmosphere 2022, 13, 32. [Google Scholar] [CrossRef]
  84. Pokrovsky, O.S.; Viers, J.; Dupre, 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. [Google Scholar] [CrossRef]
  85. Krickov, I.; Lim, A.G.; 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. [Google Scholar] [CrossRef]
  86. Bagard, M.L.; Chabaux, F.; Stille, P.; Rihs, S.; Pokrovsky, O.S.; Viers, J.; Dupré, B.; Prokushkin, A.S.; Schmitt, A.D. Seasonal variability of element fluxes in two Central Siberian rivers draining high latitude permafrost dominated areas. Geochim. Cosmochim. Acta 2011, 75, 3335–3357. [Google Scholar] [CrossRef]
  87. Manasypov, R.M.; Vorobyev, S.N.; Loiko, S.V.; Kritzkov, I.V.; Kirpotin, S.N.; Kulizhsky, S.P.; Kolesnichenko, L.G.; Zemtsov, V.A.; Sinkinov, V.V.; Pokrovsky, O.S.; et al. Seasonal dynamics of organic carbon and metals in thermokarst lakes from the discontinuous permafrost zone of western Siberia. Biogeosciences 2015, 12, 3009–3028. [Google Scholar] [CrossRef]
  88. Frey, K.E.; McClelland, J.W. Impacts of permafrost degradation on arctic river biogeochemistry. Hydrol. Process. 2009, 23, 169–182. [Google Scholar] [CrossRef]
  89. Calabrese, S.; Parolari, A.J.; Porporato, A. Hydrologic Transport of Dissolved Inorganic Carbon and Its Control on Chemical Weathering. J. Geophys. Res. Earth Surf. 2017, 122, 2016–2032. [Google Scholar] [CrossRef]
  90. Park, S.; Joe, K.S.; Han, S.H.; Kim, H.S. Characteristics of dissolved organic carbon in the leachate from Moonam Sanitary Landfill. Environ. Technol. 1999, 20, 419–424. [Google Scholar] [CrossRef]
  91. Hutchins, R.H.S.; Aukes, P.; Schiff, S.L.; Dittmar, T.; Prairie, Y.T.; del Giorgio, P.A. The optical, chemical, and molecular dissolved organic matter succession along a boreal soil-stream-river continuum. J. Geophys. Res. Biogeosci. 2017, 122, 2892–2908. [Google Scholar] [CrossRef]
  92. Shirokova, L.S.; Pokrovsky, O.S.; Kirpotin, S.N.; Desmukh, C.; Pokrovsky, B.G.; Audry, S.; Viers, J. Biogeochemistry of organic carbon, CO2, CH4, and trace elements in thermokarst water bodies in discontinuous permafrost zones of Western Siberia. Biogeochemistry 2013, 113, 573–593. [Google Scholar] [CrossRef]
  93. Payandi-Rolland, D.; Shirokova, L.S.; Nakhle, P.; Tesfa, M.; Abdou, A.; Causserand, C.; Lartiges, B.; Rols, J.-L.; Guerin, F.; Benezeth, 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. [Google Scholar] [CrossRef]
  94. Payandi-Rolland, D.; Shirokova, L.S.; Lariux, J.; Bénézeth, P.; Pokrovsky, O.S. Laboratory growth capacity of an invasive cyanobacterium (Microcystis aeruginosa) on organic substrates from surface waters of permafrost peatlands. Environ. Sci. Process. Impacts 2023, 25, 659–669. [Google Scholar] [CrossRef]
  95. Jia, J.; Sun, K.; Lu, S.; Li, M.; Wang, Y.; Yu, G.; Gao, Y. Determining whether Qinghai–Tibet Plateau waterbodies have acted like carbon sinks or sources over the past 20 years. Sci. Bull. 2022, 67, 2345–2357. [Google Scholar] [CrossRef] [PubMed]
  96. Wang, L. In-situ measurement on air–water flux of CH4, CO2 and their carbon stable isotope in lakes of northeast Tibetan Plateau. Adv. Clim. Chang. Res. 2022, 13, 279–289. [Google Scholar] [CrossRef]
  97. Yu, C.; Li, Y.; Jin, H.; Ma, Q.; Yu, Z.; Shi, K.; Li, X.; Chen, G. Organic versus Inorganic Carbon Exports from Glacier and Permafrost Watersheds in Qinghai–Tibet Plateau. Aquat. Geochem. 2021, 27, 283–296. [Google Scholar] [CrossRef]
  98. Maurischat, P.; Lehnert, L.; Zerres, H.D.V.; Tran, T.V.; Kalbitz, K.; Rinnan, Å.; Li, X.G.; Dorji, T.; Guggenberger, G. The glacial–terrestrial–fluvial pathway: A multiparametrical analysis of spatiotemporal dissolved organic matter variation in three catchments of Lake Nam Co, Tibetan Plateau. Sci. Total Environ. 2022, 838, 156542. [Google Scholar] [CrossRef]
  99. Jin, Z.D.; You, C.F.; Wang, Y.; Shi, Y. Hydrological and solute budgets of Lake Qinghai, the largest lake on the Tibetan Plateau. Quat. Int. 2010, 218, 151–156. [Google Scholar] [CrossRef]
  100. Xie, J.X.; Li, Y.; Zhai, C.X.; Li, C.; Lan, Z. CO2 absorption by alkaline soils and its implication to the global carbon cycle. Environ. Geol. 2009, 56, 953–961. [Google Scholar] [CrossRef]
  101. Golubyatnikov, L.; Kurganova, I.N.; de Gerenyu, V.L. Estimation of C-CO2 balance of natural steppe ecosystems: Khakassia and Tuva (Eastern Siberia, Russia) case studies. IOP Conf. Ser. Earth Environ. Sci. 2020, 606, 012013. [Google Scholar] [CrossRef]
  102. Karlsson, J.; Serikova, S.; Rocher-Ros, G.; Denfeld, B.; Vorobyev, S.N.; Pokrovsky, O.S. Carbon emission from Western Siberian inland waters. Nat. Commun. 2021, 12, 825. [Google Scholar] [CrossRef] [PubMed]
  103. Karabaeva, L. Surface Water Resources of the USSR: Hydrological Knowledge. Angaro-Yenisei District. Question 1. Yenisei; Hydrometeoizdat: Leningrad, Russia, 1973; Volume 16, p. 823. [Google Scholar]
  104. State Report on the State and Environmental Protection of the Republic of Tyva in 2018. Available online: http://docs.cntd.ru/document/561543920 (accessed on 25 December 2022).
  105. Grebneva, V.A. Geography of the Tuva ASSR; Tuvan Book Publishing House: Gaithersburg, MD, USA, 1972; Volume 2, 132p. [Google Scholar]
  106. Ochur-ool, A.O. Ecological and Geochemical State of the Landscapes of the Khemchik Basin (Western Tyva). Ph.D. Thesis, Candidate of Geographical Sciences, Kyzyl-Tomsk, Russia, 2016; 175p. Available online: https://vital.lib.tsu.ru/vital/access/services/Download/vtls:000563256/SOURCE1 (accessed on 18 June 2023).
  107. Makarieva, O.M.; Nesterova, N.V.; Yampolsky, G.P.; Kudymova, E.Y.; Ostashov, A.A.; Kolupaeva, A.D. The Oenc is a maximum category of aqueducts of versatile observation and unexplored mountain river Khemchik (Republic of Tyva) based on mathematical modeling methods. Eng. Surv. 2019, XIII, 36–51. [Google Scholar] [CrossRef]
  108. State Water Register: Alash River. Available online: https://textual.ru/gvr/index.php?card=212139 (accessed on 15 August 2023).
  109. State Water Register: Ak-Sug River. Available online: https://textual.ru/gvr/index.php?card=211467 (accessed on 17 August 2023).
  110. Potapova, S.A. Analysis of the physical and geographical conditions of the Ak-Sug river basin as a stage of monitoring vulnerable ecosystems of the Republic of Tuva/Geographical studies of Siberia and the Altai-Sayan transboundary region. In Proceedings of the International Scientific and Practical Conference, Veliky Novgorod, Russia, 28–29 June 2021; pp. 423–434. [Google Scholar]
  111. Chadan River: According to the State Water Register. Available online: http://www.textual.ru/gvr/index.php?card=212241 (accessed on 20 August 2023).
  112. Khayan, A.B. Geoecological state of the Durgen river. Ecology of Southern Siberia and adjacent territories. In Proceedings of the Conference “Ecological Siberia and the Limit of Territories”; 2014; pp. 154–155. Available online: https://www.elibrary.ru/item.asp?id=23627503 (accessed on 4 February 2023).
  113. Nazyn, C.D. The first information about the algae of the Torgalyg River (Tuva, Russia). Resort Base Nat. Health-Improv. Areas Tuva Adjac. Reg. 2013, 189–192. (In Russian) [Google Scholar]
  114. State Water Register: Anyyak-Chyrgaki River. Available online: https://textual.ru/gvr/index.php?card=212216 (accessed on 14 December 2021).
  115. State Water Register: Biche-Bayan-Kol River. Available online: https://textual.ru/gvr/index.php?card=212200 (accessed on 14 December 2021).
  116. Alash Plateau: The Great Soviet Encyclopedia. Available online: https://bio.1sept.ru/article.php?ID=200203204 (accessed on 18 August 2023).
  117. Shilkrot, G.S.; Kretova, S.P.; Smirnova, E.V. Ecosystem of Lake Ubsu-Nur potential and water quality. Informational problems of studying the biosphere. Pushchino 1990, 236–303. [Google Scholar]
  118. Protected Areas of Russia: Lake Tore-Khol. Available online: http://oopt.aari.ru/oopt/%D0%9E%D0%B7%D0%B5%D1%80%D0%BE%D0%A2%D0%BE%D1%80%D0%B5-%D0%A5%D0%BE%D0%BB%D1%8C (accessed on 20 December 2021).
  119. Protected Areas of Russia: Lake Cheder. Available online: http://oopt.aari.ru/oopt/%D0%9E%D0%B7%D0%B5%D1%80%D0%BE-%D0%A7%D0%B5%D0%B4%D0%B5%D1%80 (accessed on 16 August 2023).
  120. Erdynieva, L.S. Use of Natural and Ecological Resources of Lake Cheder. The World of Science, Culture, Enlightenment 2010. pp. 312–315. Available online: https://cyberleninka.ru/article/n/ispolzovanie-prirodno-ekologicheskih-resursov-ozera-cheder (accessed on 20 May 2023).
Figure 1. Location of research objects within the Sayan–Altai Mountain system (Central Asia).
Figure 1. Location of research objects within the Sayan–Altai Mountain system (Central Asia).
Water 15 03411 g001
Figure 2. River catchments of the study region. The direction of water flow is shown by different colors.
Figure 2. River catchments of the study region. The direction of water flow is shown by different colors.
Water 15 03411 g002
Figure 3. Boxplot of median pH values and IQR range (with outlies as dots) (a) during different seasons of the year (both rivers and lakes), and separately in (b) rivers and (c) lakes.
Figure 3. Boxplot of median pH values and IQR range (with outlies as dots) (a) during different seasons of the year (both rivers and lakes), and separately in (b) rivers and (c) lakes.
Water 15 03411 g003
Figure 4. The relationship between electrical conductivity and DIC concentrations in the waters of the studied lakes (the salt lake Cheder is excluded).
Figure 4. The relationship between electrical conductivity and DIC concentrations in the waters of the studied lakes (the salt lake Cheder is excluded).
Water 15 03411 g004
Figure 5. Linear relationship between DIC concentration and electrical conductivity in rivers and lakes of the Tyva region (a) in different seasons (lakes and rivers together) and (b) in different types of objects (averaged across seasons).
Figure 5. Linear relationship between DIC concentration and electrical conductivity in rivers and lakes of the Tyva region (a) in different seasons (lakes and rivers together) and (b) in different types of objects (averaged across seasons).
Water 15 03411 g005
Figure 6. Median (and IQR) boxplots of DIC concentration (mg/L−1) (a) in different types of objects (averaged across seasons) and (b) in different seasons (lakes and rivers together).
Figure 6. Median (and IQR) boxplots of DIC concentration (mg/L−1) (a) in different types of objects (averaged across seasons) and (b) in different seasons (lakes and rivers together).
Water 15 03411 g006
Figure 7. Median (and IQR) boxplots of DOC concentration (mg L−1) in different types of objects (a) in different seasons (lakes and rivers together) and (b) in different types of objects (averaged across seasons).
Figure 7. Median (and IQR) boxplots of DOC concentration (mg L−1) in different types of objects (a) in different seasons (lakes and rivers together) and (b) in different types of objects (averaged across seasons).
Water 15 03411 g007
Figure 8. Median (and IQR) boxplots of SUVA254 (a) in large and small rivers (averaged across seasons); (b) in different seasons (lakes and rivers together); (c) in different types of objects (averaged across seasons) and SR; and (d) in large and small rivers (averaged across seasons).
Figure 8. Median (and IQR) boxplots of SUVA254 (a) in large and small rivers (averaged across seasons); (b) in different seasons (lakes and rivers together); (c) in different types of objects (averaged across seasons) and SR; and (d) in large and small rivers (averaged across seasons).
Water 15 03411 g008
Figure 9. The ratio of optical densities E254/E436 in the rivers and lakes of Tyva (a) during different seasons, and (b) for different types of objects, averaged across seasons.
Figure 9. The ratio of optical densities E254/E436 in the rivers and lakes of Tyva (a) during different seasons, and (b) for different types of objects, averaged across seasons.
Water 15 03411 g009
Figure 10. Median (and IQR) boxplots of CO2 partial pressure (ppmv) in (a) rivers and (b) lakes, averaged across seasons. Different groups of rivers (F = 5.42, p = 0.02) and lakes (F = 4.44, p = 0.03) are significantly different.
Figure 10. Median (and IQR) boxplots of CO2 partial pressure (ppmv) in (a) rivers and (b) lakes, averaged across seasons. Different groups of rivers (F = 5.42, p = 0.02) and lakes (F = 4.44, p = 0.03) are significantly different.
Water 15 03411 g010
Figure 11. Median (and IQR) boxplots of CO2 partial pressure (ppmv) in (a) large and (b) small rivers across seasons, which are statistically different (large rivers, F = 12.20, p = 0.001; small rivers, F = 7.801, p = 0.001, accordingly). (c) Median (and IQR) boxplots of CO2 partial pressure (ppmv) in lakes across seasons (F = 0.9223, p = 0.46).
Figure 11. Median (and IQR) boxplots of CO2 partial pressure (ppmv) in (a) large and (b) small rivers across seasons, which are statistically different (large rivers, F = 12.20, p = 0.001; small rivers, F = 7.801, p = 0.001, accordingly). (c) Median (and IQR) boxplots of CO2 partial pressure (ppmv) in lakes across seasons (F = 0.9223, p = 0.46).
Water 15 03411 g011
Figure 12. Boxplots (median and IQR range) of fCO2 in (a) rivers, where the difference is insignificant (F = 2.30, p = 0.13); (b) lakes (F = 0.30, p = 0.743); and (c) all types of objects across seasons (F = 2.54, p = 0.064).
Figure 12. Boxplots (median and IQR range) of fCO2 in (a) rivers, where the difference is insignificant (F = 2.30, p = 0.13); (b) lakes (F = 0.30, p = 0.743); and (c) all types of objects across seasons (F = 2.54, p = 0.064).
Water 15 03411 g012
Figure 13. The relationship between (a) pCO2 and (b) fCO2 and lake water surface area of the Tyva region (all seasons together).
Figure 13. The relationship between (a) pCO2 and (b) fCO2 and lake water surface area of the Tyva region (all seasons together).
Water 15 03411 g013
Table 1. General information about the rivers sampled in this work.
Table 1. General information about the rivers sampled in this work.
NameFlow Velocity m/sDepth (m)Length (km)Average Discharge, m3/cCatchment Area (km²)Height at the Samp-Ling Point, mAverage Catchment Height (m)Slope of the Riverbed (m/km)Location
Yenisei0.25–2.62–334871020102,806650119616.8Ulug-Khem basin
Big Yenisei1.4–2.4 1.5–460559457,76663014483.1Todzhin-skaya basin; Kyzyl basin
Small Yenisei1.8–2.31–2.456341136,39563615552.8Sangilen Highlands; Ulugh-Khem basin
Tes-Hem1.1–2.11–2.175755.618,430106718427.9East Tannu-Ola
Hem-chik-0.75–23201023268850192314.4Shapshal ridge; Khemchik basin
Alash0.430.30–2172 474192020639.2Alash Plateau
Ak—Sug0.310.25–116014997.41150196626.8Alash Plateau
Chadan--98-881.5800156728.8West Tannu-Ola
Durgen0.54 0.66–193-121.71200175142.3The northern slope of East Tannu-Ola
Chaa-Hol0.28 0.5–290-320.3540169443.1The northern slope of the western Tannu-Ola
Huule (Torgalyg)-0.4–253-1090535127330.9The northern slope of the Eastern Tannu-Ola; The Central Tuva basin
Anyyak-Chyrgaki0.1730.2–252-1859800151913.9West Tannu-Ola
Dyttyg-Hem-0.2–0.834-426.91250171036.3Southern slope of East Tannu-Ola
Biche-Bayan-Kol0.340.3–0.832-15.3750122226.0Uyuk Ridge
Adyr-khem0.170.5–28.25-8.251850207666.2Alash Plateau
Table 2. General information about the studied lakes [15].
Table 2. General information about the studied lakes [15].
NameDepth (m)Water Mirror Area, km²TypeHeight (m)Location
Tore-Khol6–8 (max. 40 m)68.8Freshwater1148Ubsunur basin
Chagytai1728.6Freshwater lake1005The foot of the northern slope of the Tannu-Ola ridge
Cheder1.5–24.3Salt lake706South of the Tuva basin, a drainless depression
Thermokarst 140.3Thermokarst lake1850Alash Plateau
Thermokarst 250.1Thermokarst lake1850Alash Plateau
Table 3. The values of pCO2 (ppmv) in rivers and lakes of the region.
Table 3. The values of pCO2 (ppmv) in rivers and lakes of the region.
RiversAverageMedian
Large rivers
Yenisei689578
Big Yenisei860587
Small Yenisei790571
Tes-Khem672609
Khemchik705715
Alash563502
Small rivers
Ak—Sug809552
Chadan14701112
Durgen739.5493
Chaa-Hol778729
Huule (Torgalyg)11331003
Anyyak-Chyrgaki909932
Dyttyg-Hem840660
Biche-Bayan-Kol743743
Adyr-Khem20432105
Lakes
Tore-Khol332337
Chagytai504535
Cheder321292
Thermokarst lake 1754705
Thermokarst lake 2694806
Table 4. Measured fCO2 fluxes (g C m−2 d−1) in rivers and lakes of the Tyva Republic (Sayan–Altai Mountain region).
Table 4. Measured fCO2 fluxes (g C m−2 d−1) in rivers and lakes of the Tyva Republic (Sayan–Altai Mountain region).
AutumnWinterSpringSummerAverage ValueMedianStandard Deviation
Large rivers
Yenisei0.100.390.060.080.160.090.15
Big Yenisei 0.220.110.090.140.110.07
Small Yenisei0.190.710.030.120.260.150.31
Tes-Khem0.420.160.100.100.190.130.15
Khemchik0.190.350.000.080.160.140.15
Alash0.070.100.00−0.020.040.030.06
Small rivers
Ak—Sug0.105.600.040.111.460.102.76
Chadan 12.110.010.134.080.136.95
Durgen0.200.250.000.020.120.110.13
Chaa-Hol0.31 0.040.150.170.150.14
Huule (Torgalyg)0.140.150.050.450.200.150.17
Anyyak-Chyrgaki1.89 0.060.330.760.330.99
Dyttyg-Hem 0.510.120.030.220.120.26
Biche-Bayan-Kol
Adyr-Khem3.580.480.16−0.061.040.321.71
Lakes
Tore-Khol0.030.020.070.110.060.050.04
Chagytai0.160.010.110.000.070.060.08
Cheder0.070.090.05−0.040.040.060.06
Thermokarst lake 10.09 0.060.040.060.060.02
Thermokarst lake 20.13 0.020.060.070.060.06
Table 5. Pearson correlation coefficients between major hydrochemical and microbiological parameters of the water column, including pH, O2, S.C., DOC, DIC and optical parameters of DOM reflecting DOM quality. Sarea, E.C. and TBC represent the area of the lake surface or the river watershed, electrical conductivityf and total bacterial count, respectively.
Table 5. Pearson correlation coefficients between major hydrochemical and microbiological parameters of the water column, including pH, O2, S.C., DOC, DIC and optical parameters of DOM reflecting DOM quality. Sarea, E.C. and TBC represent the area of the lake surface or the river watershed, electrical conductivityf and total bacterial count, respectively.
All Objects
AllE.C.TwaterpHO2pCO2TBCDICDOCSUVA254E2:E3 E254:E436SRfCO2
Sarea0.030.160.230.23−0.18−0.29−0.07−0.46−0.2−0.05−0.41−0.14−0.16
E.C. 0.050.760.25−0.24−0.010.950.23−0.750.72−0.480.85−0.06
Twater −0.26−0.69−0.120.020.190.110.01−0.370.12−0.11−0.33
pH 0.54−0.330.090.70−0.01−0.770.71−0.580.61−0.04
O2 −0.13−0.160.06−0.28−0.420.53−0.460.250.09
pCO2 0.00−0.24−0.280.29−0.070.01−0.150.85
TBC 0.140.330.09−0.030.22−0.160.05
DIC 0.26−0.720.64−0.440.75−0.08
DOC 0.250.030.690.21−0.23
SUVA254 −0.610.81−0.560.08
E2:E3 −0.510.750.21
E254:E436 −0.39−0.16
SR 0.06
Rivers
E.C.TwaterpHO2pCO2TBCDICDOCSUVA254E2:E3E254:E436SRfCO2
Sarea0.220.380.170.02−0.46−0.230.180.08−0.06−0.040.05−0.06−0.40
E.C. −0.160.790.54−0.05−0.110.88−0.21−0.770.42−0.530.34−0.04
Twater −0.37−0.71−0.190.160.060.110.18−0.720.28−0.46−0.36
pH 0.67−0.23−0.180.65−0.18−0.780.61−0.550.46−0.09
O2 −0.07−0.260.16−0.16−0.480.74−0.440.660.06
pCO2 0.22−0.05−0.300.160.17−0.200.170.93
TBC 0.020.320.17−0.170.28−0.400.12
DIC −0.11−0.680.14−0.390.08−0.08
DOC 0.45−0.270.89−0.54−0.24
SUVA254 −0.360.75−0.390.11
E2:E3 −0.490.650.33
E254:E436 −0.58−0.20
SR 0.28
Lakes
E.C.TwaterpHO2pCO2TBCDICDOCSUVA254E2:E3E254:E436SRfCO2
Sarea0.950.580.720.11−0.310.190.92−0.32−0.93−0.53−0.830.980.17
E.C. 0.560.750.17−0.380.410.99−0.15−0.87−0.25−0.770.920.07
Twater −0.07−0.510.19−0.380.57−0.29−0.45−0.14−0.410.46−0.53
pH 0.59−0.540.650.710.16−0.68−0.28−0.530.750.55
O2 −0.920.540.140.20−0.35−0.03−0.300.180.16
pCO2 −0.46−0.350.050.590.080.57−0.360.10
TBC 0.450.34−0.190.39−0.160.260.36
DIC −0.17−0.84−0.19−0.760.900.05
DOC 0.480.680.67−0.420.02
SUVA254 0.580.96−0.94−0.06
E2:E3 0.59−0.59−0.42
E254:E436 −0.89−0.02
SR 0.28
Table 6. Comparison of carbon concentrations in lakes and rivers and CO2 (mean and standard deviation) emissions from water surfaces in the Tyva region to other regions of the world that have similar climate and landscape parameters.
Table 6. Comparison of carbon concentrations in lakes and rivers and CO2 (mean and standard deviation) emissions from water surfaces in the Tyva region to other regions of the world that have similar climate and landscape parameters.
SitesPeriodpCO2 (ppmv)fCO2 (g C m−2 d−1)DOC (mg L−1)DIC (mg L−1)Reference
AllAnnual cycle1495 ± 5770.46 ± 1.63.7 ± 2.929 ± 24This study
RiversAnnual cycle929 ± 6110.60 ± 1.92.7 ± 1.724.6 ± 12This study
LakesAnnual cycle385 ± 1230.06 ± 0.057.2 ± 3.8 78 ± 35This study
Thermokarst lakesAnnual cycle724 ± 3680.067 ± 0.0388 ± 1.65 ± 1.6This study
Fenghuoshan catchment, ChinaAnnual cycle1260 ± 1456.3 ± 0.917.9 ± 5.533 ± 3[19]
Thermokarst lakes, ChinaJune–September 9.5 ± 5.738 ± 35[96]
The NamCo basin and source area of Yellow River (SAID)July 2015 (melting of glacier) 2.3 ± 1.311.3 ± 10.6[97]
Yangtze River source regionBiweekly from May to October1086 ± 275 2.8 ± 0.624.4 ± 1.9[18]
Qinghai Lake and the inflowing rivers23 May 2021 1.0 ± 0.613.3 ± 7.7[21]
Three catchments within the Nam Co watershedJune/July 2018, May 2019 and September 2019 2.4 ± 0.38.8 ± 5.5[98]
76 lakes, Western Siberia LowlandMay–June, August–October1044 ± 15401.7 ± 1.716 ± 100.7 ± 0.8[1]
Lena upstream of Kirenga29 May to 17 June 2016714 ± 220.85 ± 0.0613.9 ± 1.420.0 ± 1.2[74]
Southwestern and northeastern regions of the QTP (70 lakes, four rivers and one reservoir on the QTP)20-year period (i.e., from the 2000s to the 2020s) 0.3 ± 0.2 [95]
Two saline lakes (Qinghai Lake and Hala Lake) in the Tibetan PlateauContinuously measured on 20, 23 October 2018 13.1 ± 0.4 [20]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Byzaakay, A.A.; Kolesnichenko, L.G.; Kolesnichenko, I.I.; Khovalyg, A.O.; Raudina, T.V.; Prokushkin, A.S.; Lushchaeva, I.V.; Kvasnikova, Z.N.; Vorobyev, S.N.; Pokrovsky, O.S.; et al. Dissolved Carbon Concentrations and Emission Fluxes in Rivers and Lakes of Central Asia (Sayan–Altai Mountain Region, Tyva). Water 2023, 15, 3411. https://doi.org/10.3390/w15193411

AMA Style

Byzaakay AA, Kolesnichenko LG, Kolesnichenko II, Khovalyg AO, Raudina TV, Prokushkin AS, Lushchaeva IV, Kvasnikova ZN, Vorobyev SN, Pokrovsky OS, et al. Dissolved Carbon Concentrations and Emission Fluxes in Rivers and Lakes of Central Asia (Sayan–Altai Mountain Region, Tyva). Water. 2023; 15(19):3411. https://doi.org/10.3390/w15193411

Chicago/Turabian Style

Byzaakay, Arisiya A., Larisa G. Kolesnichenko, Iury Ia. Kolesnichenko, Aldynay O. Khovalyg, Tatyana V. Raudina, Anatoly S. Prokushkin, Inna V. Lushchaeva, Zoia N. Kvasnikova, Sergey N. Vorobyev, Oleg S. Pokrovsky, and et al. 2023. "Dissolved Carbon Concentrations and Emission Fluxes in Rivers and Lakes of Central Asia (Sayan–Altai Mountain Region, Tyva)" Water 15, no. 19: 3411. https://doi.org/10.3390/w15193411

APA Style

Byzaakay, A. A., Kolesnichenko, L. G., Kolesnichenko, I. I., Khovalyg, A. O., Raudina, T. V., Prokushkin, A. S., Lushchaeva, I. V., Kvasnikova, Z. N., Vorobyev, S. N., Pokrovsky, O. S., & Kirpotin, S. (2023). Dissolved Carbon Concentrations and Emission Fluxes in Rivers and Lakes of Central Asia (Sayan–Altai Mountain Region, Tyva). Water, 15(19), 3411. https://doi.org/10.3390/w15193411

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop