The northern permafrost region has been projected to shift from a net sink to a net source of carbon under global warming. However, estimates of the contemporary net greenhouse gas (GHG) balance and budgets of the permafrost region remain highly uncertain. Here we construct the first comprehensive bottom-up budgets of CO2, CH4, and N2O across the terrestrial permafrost region using databases of more than 1000 in-situ flux measurements and a land cover-based ecosystem flux upscaling approach for the period 2000-2020. Estimates indicate that the permafrost region emitted a mean annual flux of 0.36 (-620, 652) Tg CO2-C y-1, 38 (21, 53) Tg CH4-C y-1, and 0.62 (0.03, 1.2) Tg N2O-N y-1 to the atmosphere throughout the period. While the region was a net source of CH4 and N2O, the CO2 budget was near neutral with large uncertainties. Terrestrial ecosystems remained a CO2 sink, but emissions from fire disturbances and inland waters largely offset the sink in vegetated ecosystems. Including lateral fluxes, the permafrost region was a net source of C and N, releasing 136 (-517, 821) Tg C y-1 and 3.2 (1.9, 4.8) Tg N y-1. Large uncertainty ranges in these estimates point to a need for further expansion of monitoring networks, continued data synthesis efforts, and better integration of field observations, remote sensing data, and ecosystem models to constrain the contemporary net GHG budgets of the permafrost region and track their future trajectory.

The long-term net sink of carbon (C), nitrogen (N) and greenhouse gases (GHGs) in the northern permafrost region is projected to weaken or shift under climate change. But large uncertainties remain, even on present-day GHG budgets. We compare bottom-up (data-driven upscaling, process-based models) and top-down budgets (atmospheric inversion models) of the main GHGs (CO2, CH4, and N2O) and lateral fluxes of C and N across the region over 2000-2020. Bottom-up approaches estimate higher land to atmosphere fluxes for all GHGs compared to top-down atmospheric inversions. Both bottom-up and top-down approaches respectively show a net sink of CO2 in natural ecosystems (-31 (-667, 559) and -587 (-862, -312), respectively) but sources of CH4 (38 (23, 53) and 15 (11, 18) Tg CH4-C yr-1) and N2O (0.6 (0.03, 1.2) and 0.09 (-0.19, 0.37) Tg N2O-N yr-1) in natural ecosystems. Assuming equal weight to bottom-up and top-down budgets and including anthropogenic emissions, the combined GHG budget is a source of 147 (-492, 759) Tg CO2-Ceq yr-1 (GWP100). A net CO2 sink in boreal forests and wetlands is offset by CO2 emissions from inland waters and CH4 emissions from wetlands and inland waters, with a smaller additional warming from N2O emissions. Priorities for future research include representation of inland waters in process-based models and compilation of process-model ensembles for CH4 and N2O. Discrepancies between bottom-up and top-down methods call for analyses of how prior flux ensembles impact inversion budgets, more in-situ flux observations and improved resolution in upscaling.

High-latitude tundra ecosystems are increasingly affected by climate warming. As an important fraction of soil microorganisms, fungi play essential roles in carbon (C) degradation, especially the old, chemically recalcitrant C. However, it remains obscure how fungi respond to climate warming and whether fungi, in turn, affect C stability of tundra. In a two-year winter soil warming experiment of 2 °C by snow fences, we investigated responses of fungal communities to warming in the active layer of the Alaskan tundra. Although fungal community composition, revealed by 28S rRNA gene amplicon sequencing, remained unchanged (P > 0.05), fungal functional gene composition, revealed by a microarray named GeoChip, was altered (P < 0.05). Changes in functional gene composition were linked to winter soil temperature, thaw depth, soil moisture, and gross primary productivity (Canonical Correlation Analysis, P < 0.05). Specifically, relative abundances of fungal genes encoding invertase, xylose reductase, and vanillin dehydrogenase significantly increased (P < 0.05), indicating higher C degradation capacities of fungal communities under warming. Accordingly, we detected changes of fungal gene networks under warming, including higher average path distance, lower average clustering coefficient, and lower percentage of negative links, indicating that warming potentially changed fungal interactions. Together, our study revealed higher C degradation capacities of fungal communities under short-term warming and highlights the potential impacts of fungal communities on mediating tundra ecosystem respiration, and consequently future C stability of high-latitude tundra.