Umeå University's logo

Global warming causes permafrost thawing, transferring large amounts of soil carbon into rivers, which inevitably accelerates riverine CO2 release. However, temporally and spatially explicit variations of riverine CO2 emissions remain unclear, limiting the assessment of land carbon-climate feedback. Using new and published 5685 riverine CO2 partial pressure data in the Arctic and Tibetan Plateau, we show that current riverine CO2 emission across the Northern Hemisphere permafrost zone is 200 ± 15 Tg C yr⁻1. The emission offsets 28.1 ± 2.1% of the land carbon uptake in the Northern Hemisphere permafrost zone, with large regional variability of 13.1 to 63.1%. Our findings suggest that CO2 emissions increased at a rate of 0.42 ± 0.16 Tg C yr⁻1 during 2000 to 2020, and this is primarily driven by increased precipitation and accelerated permafrost thawing under climate change. This study highlights increased riverine carbon emission and strengthening of the permafrost carbon feedback to climate after incorporating carbon release from rivers.

Carbon emission from Arctic rivers constitutes a positive feedback between the climate warming and C cycle. However, in case of rivers with extensive floodplains, the impacts of temporary water bodies and secondary channels on CO₂ exchange with atmosphere, compared to the main stem and tributaries, remain strongly understudied. In order to quantify the relative role of various water bodies of the Arctic river basin in the C cycle, the hydrochemical variables and greenhouse gases GHG concentrations and fluxes were measured within the floodplain of the largest Arctic River, Ob, in its low reaches located in the permafrost zone. These included the main stem, secondary channels, tributaries and floodplain lakes sampled over a 900 km north-south transect (25,736 km² of the main stem and adjacent floodplain area; 7893 km² water surface) during peak of spring flood (May 2023). In addition to main stem and tributaries, providing less than a half of overall C flux, floodplain lakes and secondary channels acted as important factor of C emission from the floodplain water surfaces. Multi-parametric statistical treatment of the data suggested two main processes of C emission from the Ob River floodplain waters: terrestrial organic matter-rich flooded wetlands (fens) provided elevated pCO₂, whereas the sites of possible groundwater discharge in the secondary channels decreased the CO₂ fluxes due to more alkaline environments, rich in labile metals and anionic elements. Based on available high-resolution Landsat-8 images, which matched the period of field work, it was found that the total water coverage of the floodplain during spring 2023 was 30 % of overall territory, compared to 18 % during the baseflow. Based on chamber-measured CO₂ fluxes (1.56 ± 0.47 g C-CO2 m−2 d−1), overall CO₂ emissions during 2 months of the spring flood from the entire Lower Ob River floodplain water surfaces including the main stem amounted to 0.73 ± 0.25 Tg C. Diffuse CH4 flux represented <1 % of total C flux. The main stem of the Ob River accounted for 34 % and 18 % of CO₂ and CH₄ emissions, respectively, whereas the floodplain lakes provided 59 % and 50 % of CO₂ and CH₄ emission, respectively. Considering that the low reaches of the Ob River represent >70 % of total river basin floodplain, and that during some years, the entire floodplain can be covered by water, emissions from the river – if assessed solely from summer (July–August) measurements – can be at least 3 times underestimated. It is therefore important to account for extended water surface during high water levels on Arctic rivers when assessing global riverine C emissions.

The CO₂ flux (FCO₂) from lakes to the atmosphere is a large component of the global carbon cycle and depends on the air–water CO₂ concentration gradient (ΔCO₂) and the gas transfer velocity (k). Both ΔCO₂ and k can vary on multiple timescales and understanding their contributions to FCO₂ is important for explaining variability in fluxes and developing optimal sampling designs. We measured FCO₂ and ΔCO₂ and derived k for one full ice-free period in 18 lakes using floating chambers and estimated the contributions of ΔCO₂ and k to FCO₂ variability. Generally, k contributed more than ΔCO₂ to short-term (1–9 d) FCO₂ variability. With increased temporal period, the contribution of k to FCO₂ variability decreased, and in some lakes resulted in ΔCO₂ contributing more than k to FCO₂ variability over the full ice-free period. Increased contribution of ΔCO₂ to FCO₂ variability over time occurred across all lakes but was most apparent in large-volume southern-boreal lakes and in deeper (> 2 m) parts of lakes, whereas k was linked to FCO₂ variability in shallow waters. Accordingly, knowing the variability of both k and ΔCO₂ over time and space is needed for accurate modeling of FCO₂ from these variables. We conclude that priority in FCO₂ assessments should be given to direct measurements of FCO₂ at multiple sites when possible, or otherwise from spatially distributed measurements of ΔCO₂ combined with k-models that incorporate spatial variability of lake thermal structure and meteorology.

The importance of floodplains in carbon (C) evasion from the lotic systems is especially important in continental plains of low runoff such as the organic-rich Western Siberian Lowland (WSL). To quantify the relative importance of the floodplain compared to main stem CO2 emissions, we monitored a large region of the Ob River's middle course (permafrost-free zone) over 3 months from spring to summer. We calculated seasonal water coverage using remote sensing, GIS and hydrologically-based approaches and measured CO2 emissions using floating chambers. There was a strongly pronounced seasonality in the water area's extent of the floodplain with water covering > 40 % of land during the ∼ 30 days of the most intensive spring flood (May – June) and subsequently declining to ≤ 10 % during summer (July-August). Maximal CO2 emissions were recorded in most shallow water bodies of the floodplain, notably in temporary flooded fens and birch forests. The CO2 emissions during the study period ranged from 0.2 ± 0.2 to 0.9 ± 0.2 g Cm−2 d−1 for the floodplain and 0.03 ± 0.34 g C m−2 d−1 for the Ob's main channel.

CO2 emissions from the floodplain were ∼ 163 ± 20 t C per km for the river's main stem during the 95 day study period. The partial contributions of temporary flooded zones, main stem, and permanent lakes / secondary channels to total emissions (1820 km² area) were 70, 16, and 14 %, respectively. Over spring and summer seasons, contributions from flooded zones ranged from 43 to 99 % of total CO2 emissions from water surfaces of the Ob River's middle course. Extrapolation of obtained results to the entire territory of the Ob River floodplain indicates that not accounting for the floodplain emissions may sizably—up to an order of magnitude—underestimate the CO2 emissions from riverine systems in Western Siberia during open water period. Future work on the Ob River floodplain in the permafrost-bearing zone should be prioritized and would allow adequate upscaling of C emission from this environmentally important territory.

Current estimates of carbon dioxide (CO2) evasion from Arctic lakes are highly uncertain because few studies integrate seasonal variability, specifically evasion during spring ice-melt. We quantified annual CO2 evasion for 14 clear-water Arctic lakes in Northern Sweden through mass balance (ice-melt period) and high-frequency loggers (open-water period). On average, 80% (SD: ± 18) of annual CO2 evasion occurred within 10 d following ice-melt. The contribution of the ice-melt period to annual CO2 evasion was high compared to earlier studies of Arctic lakes (47% ± 32%). Across all lakes, the proportion of ice-melt : annual CO2 evasion was negatively related to the dissolved organic carbon concentration and positively related to the mean depth of the lakes. The results emphasize the need for measurements of CO2 exchange at ice-melt to accurately quantify CO2 evasion from Arctic lakes.

Climate-sensitive northern cryosphere inland waters emit greenhouse gases (GHGs) into the atmosphere, yet their total emissions remain poorly constrained. We present a data-driven synthesis of GHG emissions from northern cryosphere inland waters considering water body types, cryosphere zones, and seasonality. We find that annual GHG emissions are dominated by carbon dioxide ([Formula: see text] teragrams of CO2; [Formula: see text]) and methane ([Formula: see text] teragrams of CH4), while the nitrous oxide emission ([Formula: see text] gigagrams of N2O) is minor. The annual CO2-equivalent (CO2e) GHG emissions from northern cryosphere inland waters total [Formula: see text] or [Formula: see text] petagrams of CO2e using the 100- or 20-year global warming potentials, respectively. Rivers emit 64% more CO2e GHGs than lakes, despite having only one-fifth of their surface area. The continuous permafrost zone contributed half of the inland water GHG emissions. Annual CO2e emissions from northern cryosphere inland waters exceed the region's terrestrial net ecosystem exchange, highlighting the important role of inland waters in the cryospheric land-aquatic continuum under a warming climate.

Global change impacts important environmental drivers for pelagic gross primary production (GPP) in northern lakes, such as temperature, light, nutrient, and inorganic carbon availability. Separate and/or synergistic impacts of these environmental drivers on pelagic GPP remain largely unresolved. Here, we assess key drivers of pelagic GPP by combining detailed depth profiles of summer pelagic GPP with environmental and climatic data across 45 small and shallow lakes across northern Sweden (20 boreal, 6 subarctic, and 19 arctic lakes). We found that across lakes summer pelagic GPP was strongest associated with lake water temperatures, lake carbon dioxide (CO2) concentrations impacted by lake water pH, and further moderated by dissolved organic carbon (DOC) concentrations influencing light and nutrient conditions. We further used this dataset to assess the extent of additional DOC-induced warming of epilimnia (here named internal warming), which was especially pronounced in shallow lakes (decreasing 0.96°C for every decreasing m in average lake depth) and increased with higher concentrations of DOC. Additionally, the total pools and relative proportion of dissolved inorganic carbon and DOC, further influenced pelagic GPP with drivers differing slightly among the boreal, subarctic and Arctic biomes. Our study provides novel insights in that global change affects pelagic GPP in northern lakes not only by modifying the organic carbon cycle and light and nutrient conditions, but also through modifications of inorganic carbon supply and temperature. Considering the large-scale impacts and similarities of global warming, browning and recovery from acidification of lakes at higher latitudes throughout the northern hemisphere, these changes are likely to operate on a global scale.

Inland waters play a critical role in the carbon cycle by emitting significant amounts of land-exported carbon to the atmosphere. While carbon gas emissions from individual aquatic systems have been extensively studied, how networks of connected streams and lakes regulate integrated fluxes of organic and inorganic forms remain poorly understood. Here, we develop a process-based model to simulate the fate of terrestrial dissolved organic carbon (DOC) and carbon dioxide (CO₂) in artificial inland water networks with variable topology, hydrology, and DOC reactivity. While the role of lakes is highly dependent on DOC reactivity, we find that the mineralization of terrestrial DOC is more efficient in lake-rich networks. Regardless of typology and hydrology, terrestrial CO₂ is emitted almost entirely within the network boundary. Consequently, the proportion of exported terrestrial carbon emitted from inland water networks increases with the CO₂ versus DOC export ratio. Overall, our results suggest that CO₂ emissions from inland waters are governed by interactions between the relative amount and reactivity of terrestrial DOC and CO₂ inputs and the network configuration of recipient lakes and streams.

Accounting for temporal changes in carbon dioxide (CO₂) effluxes from freshwaters remains a challenge for global and regional carbon budgets. Here, we synthesize 171 site-months of flux measurements of CO₂ based on the eddy covariance method from 13 lakes and reservoirs in the Northern Hemisphere, and quantify dynamics at multiple temporal scales. We found pronounced sub-annual variability in CO₂ flux at all sites. By accounting for diel variation, only 11% of site-months were net daily sinks of CO₂. Annual CO₂ emissions had an average of 25% (range 3%-58%) interannual variation. Similar to studies on streams, nighttime emissions regularly exceeded daytime emissions. Biophysical regulations of CO₂ flux variability were delineated through mutual information analysis. Sample analysis of CO₂ fluxes indicate the importance of continuous measurements. Better characterization of short- and long-term variability is necessary to understand and improve detection of temporal changes of CO₂ fluxes in response to natural and anthropogenic drivers. Our results indicate that existing global lake carbon budgets relying primarily on daytime measurements yield underestimates of net emissions.

Despite the importance of small and medium size rivers of Siberian boreal zone in greenhouse gases (GHG) emission, major knowledge gaps exist regarding its temporal variability and controlling mechanisms. Here we sampled 11 pristine rivers of the southern taiga biome (western Siberia Lowland, WSL), ranging in watershed area from 0.8 to 119,000 km2, to reveal temporal pattern and examine main environmental controllers of GHG emissions from the river water surfaces. Floating chamber measurements demonstrated that CO2 emissions from water surface decreased by 2 to 4-folds from spring to summer and autumn, were independent of the size of the watershed and stream order and did not exhibit sizable (>30 %, regardless of season) variations between day and night. The CH4 concentrations and fluxes increased in the order “spring ≤ summer < autumn” and ranged from 1 to 15 μmol L−1 and 5 to 100 mmol m−2 d−1, respectively. The CO2 concentrations and fluxes (range from 100 to 400 μmol L−1 and 1 to 4 g C m−2 d−1, respectively) were positively correlated with dissolved and particulate organic carbon, total nitrogen and bacterial number of the water column. The CH4 concentrations and fluxes were positively correlated with phosphate and ammonia concentrations. Of the landscape parameters, positive correlations were detected between riparian vegetation biomass and CO2 and CH4 concentrations. Over the six-month open-water period, areal emissions of C (>99.5 % CO2; <0.5 % CH4) from the watersheds of 11 rivers were equal to the total downstream C export in this part of the WSL. Based on correlations between environmental controllers (watershed land cover and the water column parameters), we hypothesize that the fluxes are largely driven by riverine mineralization of terrestrial dissolved and particulate OC, coupled with respiration at the river bottom and riparian sediments. It follows that, under climate warming scenario, most significant changes in GHG regimes of western Siberian rivers located in permafrost-free zone may occur due to changes in the riparian zone vegetation and water coverage of the floodplains.