Increasing wildfires threaten historic carbon sink of boreal forest soils.

Boreal forest fires emit large amounts of carbon into the atmosphere primarily through the combustion of soil organic matter1-3. During each fire, a portion of this soil beneath the burned layer can escape combustion, leading to a net accumulation of carbon in forests over multiple fire events4. Climate warming and drying has led to more severe and frequent forest fires5-7, which threaten to shift the carbon balance of the boreal ecosystem from net accumulation to net loss1, resulting in a positive climate feedback8. This feedback will occur if organic-soil carbon that escaped burning in previous fires, termed 'legacy carbon', combusts. Here we use soil radiocarbon dating to quantitatively assess legacy carbon loss in the 2014 wildfires in the Northwest Territories of Canada2. We found no evidence for the combustion of legacy carbon in forests that were older than the historic fire-return interval of northwestern boreal forests9. In forests that were in dry landscapes and less than 60 years old at the time of the fire, legacy carbon that had escaped burning in the previous fire cycle was combusted. We estimate that 0.34 million hectares of young forests (<60 years) that burned in the 2014 fires could have experienced legacy carbon combustion. This implies a shift to a domain of carbon cycling in which these forests become a net source-instead of a sink-of carbon to the atmosphere over consecutive fires. As boreal wildfires continue to increase in size, frequency and intensity7, the area of young forests that experience legacy carbon combustion will probably increase and have a key role in shifting the boreal carbon balance.

Diminishing warming effects on plant phenology over time.

Plant phenology, the timing of recurrent biological events, shows key and complex response to climate warming, with consequences for ecosystem functions and services. A key challenge for predicting plant phenology under future climates is to determine whether the phenological changes will persist with more intensive and long-term warming. Here, we conducted a meta-analysis of 103 experimental warming studies around the globe to investigate the responses of four phenophases - leaf-out, first flowering, last flowering, and leaf coloring. We showed that warming advanced leaf-out and flowering but delayed leaf coloring across herbaceous and woody plants. As the magnitude of warming increased, the response of most plant phenophases gradually leveled off for herbaceous plants, while phenology responded in proportion to warming in woody plants. We also found that the experimental effects of warming on plant phenology diminished over time across all phenophases. Specifically, the rate of changes in first flowering for herbaceous species, as well as leaf-out and leaf coloring for woody species, decreased as the experimental duration extended. Together, these results suggest that the real-world impact of global warming on plant phenology will diminish over time as temperatures continue to increase.

Abrupt permafrost thaw drives spatially heterogeneous soil moisture and carbon dioxide fluxes in upland tundra.

Permafrost thaw causes the seasonally thawed active layer to deepen, causing the Arctic to shift toward carbon release as soil organic matter becomes susceptible to decomposition. Ground subsidence initiated by ice loss can cause these soils to collapse abruptly, rapidly shifting soil moisture as microtopography changes and also accelerating carbon and nutrient mobilization. The uncertainty of soil moisture trajectories during thaw makes it difficult to predict the role of abrupt thaw in suppressing or exacerbating carbon losses. In this study, we investigated the role of shifting soil moisture conditions on carbon dioxide fluxes during a 13-year permafrost warming experiment that exhibited abrupt thaw. Warming deepened the active layer differentially across treatments, leading to variable rates of subsidence and formation of thermokarst depressions. In turn, differential subsidence caused a gradient of moisture conditions, with some plots becoming consistently inundated with water within thermokarst depressions and others exhibiting generally dry, but more variable soil moisture conditions outside of thermokarst depressions. Experimentally induced permafrost thaw initially drove increasing rates of growing season gross primary productivity (GPP), ecosystem respiration (Reco ), and net ecosystem exchange (NEE) (higher carbon uptake), but the formation of thermokarst depressions began to reverse this trend with a high level of spatial heterogeneity. Plots that subsided at the slowest rate stayed relatively dry and supported higher CO2 fluxes throughout the 13-year experiment, while plots that subsided very rapidly into the center of a thermokarst feature became consistently wet and experienced a rapid decline in growing season GPP, Reco , and NEE (lower carbon uptake or carbon release). These findings indicate that Earth system models, which do not simulate subsidence and often predict drier active layer conditions, likely overestimate net growing season carbon uptake in abruptly thawing landscapes.

Ecosystem and soil respiration radiocarbon detects old carbon release as a fingerprint of warming and permafrost destabilization with climate change.

The permafrost region has accumulated organic carbon in cold and waterlogged soils over thousands of years and now contains three times as much carbon as the atmosphere. Global warming is degrading permafrost with the potential to accelerate climate change as increased microbial decomposition releases soil carbon as greenhouse gases. A 19-year time series of soil and ecosystem respiration radiocarbon from Alaska provides long-term insight into changing permafrost soil carbon dynamics in a warmer world. Nine per cent of ecosystem respiration and 23% of soil respiration observations had radiocarbon values more than 50‰ lower than the atmospheric value. Furthermore, the overall trend of ecosystem and soil respiration radiocarbon values through time decreased more than atmospheric radiocarbon values did, indicating that old carbon degradation was enhanced. Boosted regression tree analyses showed that temperature and moisture environmental variables had the largest relative influence on lower radiocarbon values. This suggested that old carbon degradation was controlled by warming/permafrost thaw and soil drying together, as waterlogged soil conditions could protect soil carbon from microbial decomposition even when thawed. Overall, changing conditions increasingly favoured the release of old carbon, which is a definitive fingerprint of an accelerating feedback to climate change as a consequence of warming and permafrost destabilization. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Making the case for an International Decade of Radiocarbon.

Radiocarbon (14C) is a critical tool for understanding the global carbon cycle. During the Anthropocene, two new processes influenced 14C in atmospheric, land and ocean carbon reservoirs. First, 14C-free carbon derived from fossil fuel burning has diluted 14C, at rates that have accelerated with time. Second, 'bomb' 14C produced by atmospheric nuclear weapon tests in the mid-twentieth century provided a global isotope tracer that is used to constrain rates of air-sea gas exchange, carbon turnover, large-scale atmospheric and ocean transport, and other key C cycle processes. As we write, the 14C/12C ratio of atmospheric CO2 is dropping below pre-industrial levels, and the rate of decline in the future will depend on global fossil fuel use and net exchange of bomb 14C between the atmosphere, ocean and land. This milestone coincides with a rapid increase in 14C measurement capacity worldwide. Leveraging future 14C measurements to understand processes and test models requires coordinated international effort-a 'decade of radiocarbon' with multiple goals: (i) filling observational gaps using archives, (ii) building and sustaining observation networks to increase measurement density across carbon reservoirs, (iii) developing databases, synthesis and modelling tools and (iv) establishing metrics for identifying and verifying changes in carbon sources and sinks. This article is part of the Theo Murphy meeting issue 'Radiocarbon in the Anthropocene'.

Permafrost thaw with warming reduces microbial metabolic capacities in subsurface soils.

Microorganisms are major constituents of the total biomass in permafrost regions, whose underlain soils are frozen for at least two consecutive years. To understand potential microbial responses to climate change, here we examined microbial community compositions and functional capacities across four soil depths in an Alaska tundra site. We showed that a 5-year warming treatment increased soil thaw depth by 25.7% (p = .011) within the deep organic layer (15-25 cm). Concurrently, warming reduced 37% of bacterial abundance and 64% of fungal abundances in the deep organic layer, while it did not affect microbial abundance in other soil layers (i.e., 0-5, 5-15, and 45-55 cm). Warming treatment altered fungal community composition and microbial functional structure (p < .050), but not bacterial community composition. Using a functional gene array, we found that the relative abundances of a variety of carbon (C)-decomposing, iron-reducing, and sulphate-reducing genes in the deep organic layer were decreased, which was not observed by the shotgun sequencing-based metagenomics analysis of those samples. To explain the reduced metabolic capacities, we found that warming treatment elicited higher deterministic environmental filtering, which could be linked to water-saturated time, soil moisture, and soil thaw duration. In contrast, plant factors showed little influence on microbial communities in subsurface soils below 15 cm, despite a 25.2% higher (p < .05) aboveground plant biomass by warming treatment. Collectively, we demonstrate that microbial metabolic capacities in subsurface soils are reduced, probably arising from enhanced thaw by warming.

Lowering water table reduces carbon sink strength and carbon stocks in northern peatlands.

Peatlands at high latitudes have accumulated >400 Pg carbon (C) because saturated soil and cold temperatures suppress C decomposition. This substantial amount of C in Arctic and Boreal peatlands is potentially subject to increased decomposition if the water table (WT) decreases due to climate change, including permafrost thaw-related drying. Here, we optimize a version of the Organizing Carbon and Hydrology In Dynamic Ecosystems model (ORCHIDEE-PCH4) using site-specific observations to investigate changes in CO2 and CH4 fluxes as well as C stock responses to an experimentally manipulated decrease of WT at six northern peatlands. The unmanipulated control peatlands, with the WT <20 cm on average (seasonal max up to 45 cm) below the surface, currently act as C sinks in most years (58 ± 34 g C m-2  year-1 ; including 6 ± 7 g C-CH4 m-2  year-1 emission). We found, however, that lowering the WT by 10 cm reduced the CO2 sink by 13 ± 15 g C m-2  year-1 and decreased CH4 emission by 4 ± 4 g CH4 m-2  year-1 , thus accumulating less C over 100 years (0.2 ± 0.2 kg C m-2 ). Yet, the reduced emission of CH4 , which has a larger greenhouse warming potential, resulted in a net decrease in greenhouse gas balance by 310 ± 360 g CO2-eq  m-2  year-1 . Peatlands with the initial WT close to the soil surface were more vulnerable to C loss: Non-permafrost peatlands lost >2 kg C m-2 over 100 years when WT is lowered by 50 cm, while permafrost peatlands temporally switched from C sinks to sources. These results highlight that reductions in C storage capacity in response to drying of northern peatlands are offset in part by reduced CH4 emissions, thus slightly reducing the positive carbon climate feedbacks of peatlands under a warmer and drier future climate scenario.

Responses of tundra soil microbial communities to half a decade of experimental warming at two critical depths.

Northern-latitude tundra soils harbor substantial carbon (C) stocks that are highly susceptible to microbial degradation with rising global temperatures. Understanding the magnitude and direction (e.g., C release or sequestration) of the microbial responses to warming is necessary to accurately model climate change. In this study, Alaskan tundra soils were subjected to experimental in situ warming by ∼1.1 °C above ambient temperature, and the microbial communities were evaluated using metagenomics after 4.5 years, at 2 depths: 15 to 25 cm (active layer at outset of the experiment) and 45 to 55 cm (transition zone at the permafrost/active layer boundary at the outset of the experiment). In contrast to small or insignificant shifts after 1.5 years of warming, 4.5 years of warming resulted in significant changes to the abundances of functional traits and the corresponding taxa relative to control plots (no warming), and microbial shifts differed qualitatively between the two soil depths. At 15 to 25 cm, increased abundances of carbohydrate utilization genes were observed that correlated with (increased) measured ecosystem carbon respiration. At the 45- to 55-cm layer, increased methanogenesis potential was observed, which corresponded with a 3-fold increase in abundance of a single archaeal clade of the Methanosarcinales order, increased annual thaw duration (45.3 vs. 79.3 days), and increased CH4 emissions. Collectively, these data demonstrate that the microbial responses to warming in tundra soil are rapid and markedly different between the 2 critical soil layers evaluated, and identify potential biomarkers for the corresponding microbial processes that could be important in modeling.

Winter warming rapidly increases carbon degradation capacities of fungal communities in tundra soil: Potential consequences on carbon stability.

High-latitude tundra ecosystems are increasingly affected by climate warming. As an important fraction of soil microorganisms, fungi play essential roles in carbon degradation, especially the old, chemically recalcitrant carbon. However, it remains obscure how fungi respond to climate warming and whether fungi, in turn, affect carbon stability of tundra. In a 2-year winter soil warming experiment of 2°C by snow fences, we investigated responses of fungal communities to warming in the active layer of an Alaskan tundra. Although fungal community composition, revealed by the 28S rRNA gene amplicon sequencing, remained unchanged (p > .05), fungal functional gene composition, revealed by a microarray named GeoChip, was altered (p < .05). Changes in functional gene composition were linked to winter soil temperature, thaw depth, soil moisture, and gross primary productivity (canonical correlation analysis, p < .05). Specifically, relative abundances of fungal genes encoding invertase, xylose reductase and vanillin dehydrogenase significantly increased (p < .05), indicating higher carbon degradation capacities of fungal communities under warming. Accordingly, we detected changes in 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 reveals higher carbon degradation capacities of fungal communities under short-term warming and highlights the potential impacts of fungal communities on tundra ecosystem respiration, and consequently future carbon stability of high-latitude tundra.

Lower soil moisture and deep soil temperatures in thermokarst features increase old soil carbon loss after 10 years of experimental permafrost warming.

Almost half of the global terrestrial soil carbon (C) is stored in the northern circumpolar permafrost region, where air temperatures are increasing two times faster than the global average. As climate warms, permafrost thaws and soil organic matter becomes vulnerable to greater microbial decomposition. Long-term soil warming of ice-rich permafrost can result in thermokarst formation that creates variability in environmental conditions. Consequently, plant and microbial proportional contributions to ecosystem respiration may change in response to long-term soil warming. Natural abundance δ13 C and Δ14 C of aboveground and belowground plant material, and of young and old soil respiration were used to inform a mixing model to partition the contribution of each source to ecosystem respiration fluxes. We employed a hierarchical Bayesian approach that incorporated gross primary productivity and environmental drivers to constrain source contributions. We found that long-term experimental permafrost warming introduced a soil hydrology component that interacted with temperature to affect old soil C respiration. Old soil C loss was suppressed in plots with warmer deep soil temperatures because they tended to be wetter. When soil volumetric water content significantly decreased in 2018 relative to 2016 and 2017, the dominant respiration sources shifted from plant aboveground and young soil respiration to old soil respiration. The proportion of ecosystem respiration from old soil C accounted for up to 39% of ecosystem respiration and represented a 30-fold increase compared to the wet-year average. Our findings show that thermokarst formation may act to moderate microbial decomposition of old soil C when soil is highly saturated. However, when soil moisture decreases, a higher proportion of old soil C is vulnerable to decomposition and can become a large flux to the atmosphere. As permafrost systems continue to change with climate, we must understand the thresholds that may propel these systems from a C sink to a source.

Statistical upscaling of ecosystem CO2 fluxes across the terrestrial tundra and boreal domain: Regional patterns and uncertainties.

The regional variability in tundra and boreal carbon dioxide (CO2 ) fluxes can be high, complicating efforts to quantify sink-source patterns across the entire region. Statistical models are increasingly used to predict (i.e., upscale) CO2 fluxes across large spatial domains, but the reliability of different modeling techniques, each with different specifications and assumptions, has not been assessed in detail. Here, we compile eddy covariance and chamber measurements of annual and growing season CO2 fluxes of gross primary productivity (GPP), ecosystem respiration (ER), and net ecosystem exchange (NEE) during 1990-2015 from 148 terrestrial high-latitude (i.e., tundra and boreal) sites to analyze the spatial patterns and drivers of CO2 fluxes and test the accuracy and uncertainty of different statistical models. CO2 fluxes were upscaled at relatively high spatial resolution (1 km2 ) across the high-latitude region using five commonly used statistical models and their ensemble, that is, the median of all five models, using climatic, vegetation, and soil predictors. We found the performance of machine learning and ensemble predictions to outperform traditional regression methods. We also found the predictive performance of NEE-focused models to be low, relative to models predicting GPP and ER. Our data compilation and ensemble predictions showed that CO2 sink strength was larger in the boreal biome (observed and predicted average annual NEE -46 and -29 g C m-2  yr-1 , respectively) compared to tundra (average annual NEE +10 and -2 g C m-2  yr-1 ). This pattern was associated with large spatial variability, reflecting local heterogeneity in soil organic carbon stocks, climate, and vegetation productivity. The terrestrial ecosystem CO2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO2 sink during 1990-2015, although uncertainty remains high.

Dependence of the evolution of carbon dynamics in the northern permafrost region on the trajectory of climate change.

We conducted a model-based assessment of changes in permafrost area and carbon storage for simulations driven by RCP4.5 and RCP8.5 projections between 2010 and 2299 for the northern permafrost region. All models simulating carbon represented soil with depth, a critical structural feature needed to represent the permafrost carbon-climate feedback, but that is not a universal feature of all climate models. Between 2010 and 2299, simulations indicated losses of permafrost between 3 and 5 million km2 for the RCP4.5 climate and between 6 and 16 million km2 for the RCP8.5 climate. For the RCP4.5 projection, cumulative change in soil carbon varied between 66-Pg C (1015-g carbon) loss to 70-Pg C gain. For the RCP8.5 projection, losses in soil carbon varied between 74 and 652 Pg C (mean loss, 341 Pg C). For the RCP4.5 projection, gains in vegetation carbon were largely responsible for the overall projected net gains in ecosystem carbon by 2299 (8- to 244-Pg C gains). In contrast, for the RCP8.5 projection, gains in vegetation carbon were not great enough to compensate for the losses of carbon projected by four of the five models; changes in ecosystem carbon ranged from a 641-Pg C loss to a 167-Pg C gain (mean, 208-Pg C loss). The models indicate that substantial net losses of ecosystem carbon would not occur until after 2100. This assessment suggests that effective mitigation efforts during the remainder of this century could attenuate the negative consequences of the permafrost carbon-climate feedback.

Model parameterization to represent processes at unresolved scales and changing properties of evolving systems.

Modeling has become an indispensable tool for scientific research. However, models generate great uncertainty when they are used to predict or forecast ecosystem responses to global change. This uncertainty is partly due to parameterization, which is an essential procedure for model specification via defining parameter values for a model. The classic doctrine of parameterization is that a parameter is constant. However, it is commonly known from modeling practice that a model that is well calibrated for its parameters at one site may not simulate well at another site unless its parameters are tuned again. This common practice implies that parameter values have to vary with sites. Indeed, parameter values that are estimated using a statistically rigorous approach, that is, data assimilation, vary with time, space, and treatments in global change experiments. This paper illustrates that varying parameters is to account for both processes at unresolved scales and changing properties of evolving systems. A model, no matter how complex it is, could not represent all the processes of one system at resolved scales. Interactions of processes at unresolved scales with those at resolved scales should be reflected in model parameters. Meanwhile, it is pervasively observed that properties of ecosystems change over time, space, and environmental conditions. Parameters, which represent properties of a system under study, should change as well. Tuning has been practiced for many decades to change parameter values. Yet this activity, unfortunately, did not contribute to our knowledge on model parameterization at all. Data assimilation makes it possible to rigorously estimate parameter values and, consequently, offers an approach to understand which, how, how much, and why parameters vary. To fully understand those issues, extensive research is required. Nonetheless, it is clear that changes in parameter values lead to different model predictions even if the model structure is the same.

Drainage enhances modern soil carbon contribution but reduces old soil carbon contribution to ecosystem respiration in tundra ecosystems.

Warming temperatures are likely to accelerate permafrost thaw in the Arctic, potentially leading to the release of old carbon previously stored in deep frozen soil layers. Deeper thaw depths in combination with geomorphological changes due to the loss of ice structures in permafrost, may modify soil water distribution, creating wetter or drier soil conditions. Previous studies revealed higher ecosystem respiration rates under drier conditions, and this study investigated the cause of the increased ecosystem respiration rates using radiocarbon signatures of respired CO2 from two drying manipulation experiments: one in moist and the other in wet tundra. We demonstrate that higher contributions of CO2 from shallow soil layers (0-15 cm; modern soil carbon) drive the increased ecosystem respiration rates, while contributions from deeper soil (below 15 cm from surface and down to the permafrost table; old soil carbon) decreased. These changes can be attributed to more aerobic conditions in shallow soil layers, but also the soil temperature increases in shallow layers but decreases in deep layers, due to the altered thermal properties of organic soils. Decreased abundance of aerenchymatous plant species following drainage in wet tundra reduced old carbon release but increased aboveground plant biomass elevated contributions of autotrophic respiration to ecosystem respiration. The results of this study suggest that drier soils following drainage may accelerate decomposition of modern soil carbon in shallow layers but slow down decomposition of old soil carbon in deep layers, which may offset some of the old soil carbon loss from thawing permafrost.

Microbial functional diversity covaries with permafrost thaw-induced environmental heterogeneity in tundra soil.

Permafrost soil in high latitude tundra is one of the largest terrestrial carbon (C) stocks and is highly sensitive to climate warming. Understanding microbial responses to warming-induced environmental changes is critical to evaluating their influences on soil biogeochemical cycles. In this study, a functional gene array (i.e., geochip 4.2) was used to analyze the functional capacities of soil microbial communities collected from a naturally degrading permafrost region in Central Alaska. Varied thaw history was reported to be the main driver of soil and plant differences across a gradient of minimally, moderately, and extensively thawed sites. Compared with the minimally thawed site, the number of detected functional gene probes across the 15-65 cm depth profile at the moderately and extensively thawed sites decreased by 25% and 5%, while the community functional gene ß-diversity increased by 34% and 45%, respectively, revealing decreased functional gene richness but increased community heterogeneity along the thaw progression. Particularly, the moderately thawed site contained microbial communities with the highest abundances of many genes involved in prokaryotic C degradation, ammonification, and nitrification processes, but lower abundances of fungal C decomposition and anaerobic-related genes. Significant correlations were observed between functional gene abundance and vascular plant primary productivity, suggesting that plant growth and species composition could be co-evolving traits together with microbial community composition. Altogether, this study reveals the complex responses of microbial functional potentials to thaw-related soil and plant changes and provides information on potential microbially mediated biogeochemical cycles in tundra ecosystems.

Cross-scale controls on carbon emissions from boreal forest megafires.

Climate warming and drying is associated with increased wildfire disturbance and the emergence of megafires in North American boreal forests. Changes to the fire regime are expected to strongly increase combustion emissions of carbon (C) which could alter regional C balance and positively feedback to climate warming. In order to accurately estimate C emissions and thereby better predict future climate feedbacks, there is a need to understand the major sources of heterogeneity that impact C emissions at different scales. Here, we examined 211 field plots in boreal forests dominated by black spruce (Picea mariana) or jack pine (Pinus banksiana) of the Northwest Territories (NWT), Canada after an unprecedentedly large area burned in 2014. We assessed both aboveground and soil organic layer (SOL) combustion, with the goal of determining the major drivers in total C emissions, as well as to develop a high spatial resolution model to scale emissions in a relatively understudied region of the boreal forest. On average, 3.35 kg C m-2 was combusted and almost 90% of this was from SOL combustion. Our results indicate that black spruce stands located at landscape positions with intermediate drainage contribute the most to C emissions. Indices associated with fire weather and date of burn did not impact emissions, which we attribute to the extreme fire weather over a short period of time. Using these results, we estimated a total of 94.3 Tg C emitted from 2.85 Mha of burned area across the entire 2014 NWT fire complex, which offsets almost 50% of mean annual net ecosystem production in terrestrial ecosystems of Canada. Our study also highlights the need for fine-scale estimates of burned area that represent small water bodies and regionally specific calibrations of combustion that account for spatial heterogeneity in order to accurately model emissions at the continental scale.

Biotic responses buffer warming-induced soil organic carbon loss in Arctic tundra.

Climate warming can result in both abiotic (e.g., permafrost thaw) and biotic (e.g., microbial functional genes) changes in Arctic tundra. Recent research has incorporated dynamic permafrost thaw in Earth system models (ESMs) and indicates that Arctic tundra could be a significant future carbon (C) source due to the enhanced decomposition of thawed deep soil C. However, warming-induced biotic changes may influence biologically related parameters and the consequent projections in ESMs. How model parameters associated with biotic responses will change under warming and to what extent these changes affect projected C budgets have not been carefully examined. In this study, we synthesized six data sets over 5 years from a soil warming experiment at the Eight Mile Lake, Alaska, into the Terrestrial ECOsystem (TECO) model with a probabilistic inversion approach. The TECO model used multiple soil layers to track dynamics of thawed soil under different treatments. Our results show that warming increased light use efficiency of vegetation photosynthesis but decreased baseline (i.e., environment-corrected) turnover rates of SOC in both the fast and slow pools in comparison with those under control. Moreover, the parameter changes generally amplified over time, suggesting processes of gradual physiological acclimation and functional gene shifts of both plants and microbes. The TECO model predicted that field warming from 2009 to 2013 resulted in cumulative C losses of 224 or 87 g/m2 , respectively, without or with changes in those parameters. Thus, warming-induced parameter changes reduced predicted soil C loss by 61%. Our study suggests that it is critical to incorporate biotic changes in ESMs to improve the model performance in predicting C dynamics in permafrost regions.

Fuel-reduction management alters plant composition, carbon and nitrogen pools, and soil thaw in Alaskan boreal forest.

Increasing wildfire activity in Alaska's boreal forests has led to greater fuel-reduction management. Management has been implemented to reduce wildfire spread, but the ecological impacts of these practices are poorly known. We quantified the effects of hand-thinning and shearblading on above- and belowground stand characteristics, plant species composition, carbon (C) and nitrogen (N) pools, and soil thaw across 19 sites dominated by black spruce (Picea mariana) in interior Alaska treated 2-12 years prior to sampling. The density of deciduous tree seedlings was significantly higher in shearbladed areas compared to unmanaged forest (6.4 vs. 0.1 stems/m2 ), and unmanaged stands exhibited the highest mean density of conifer seedlings and layers (1.4 stems/m2 ). Understory plant community composition was most similar between unmanaged and thinned stands. Shearblading resulted in a near complete loss of aboveground tree biomass C pools while thinning approximately halved the C pool size (1.2 kg C/m2 compared to 3.1 kg C/m2 in unmanaged forest). Significantly smaller soil organic layer (SOL) C and N pools were observed in shearbladed stands (3.2 kg C/m2 and 116.8 g N/m2 ) relative to thinned (6.0 kg C/m2 and 192.2 g N/m2 ) and unmanaged (5.9 kg C/m2 and 178.7 g N/m2 ) stands. No difference in C and N pool sizes in the uppermost 10 cm of mineral soil was observed among stand types. Total C stocks for measured pools was 2.6 kg C/m2 smaller in thinned stands and 5.8 kg C/m2 smaller in shearbladed stands when compared to unmanaged forest. Soil thaw depth averaged 13 cm deeper in thinned areas and 46 cm deeper in shearbladed areas relative to adjacent unmanaged stands, although variability was high across sites. Deeper soil thaw was linked to shallower SOL depth for unmanaged stands and both management types, however for any given SOL depth, thaw tended to be deeper in shearbladed areas compared to unmanaged forest. These findings indicate that fuel-reduction management alters plant community composition, C and N pools, and soil thaw depth, with consequences for ecosystem structure and function beyond those intended for fire management.