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.

Wildfire combustion and carbon stocks in the southern Canadian boreal forest: Implications for a warming world.

Boreal wildfires are increasing in intensity, extent, and frequency, potentially intensifying carbon emissions and transitioning the region from a globally significant carbon sink to a source. The productive southern boreal forests of central Canada already experience relatively high frequencies of fire, and as such may serve as an analog of future carbon dynamics for more northern forests. Fire-carbon dynamics in southern boreal systems are relatively understudied, with limited investigation into the drivers of pre-fire carbon stocks or subsequent combustion. As part of NASA's Arctic-Boreal Vulnerability Experiment, we sampled 79 stands (47 burned, 32 unburned) throughout central Saskatchewan to characterize above- and belowground carbon stocks and combustion rates in relation to historical land use, vegetation characteristics, and geophysical attributes. We found southern boreal forests emitted an average of 3.3 ± 1.1 kg C/m2 from field sites. The emissions from southern boreal stands varied as a function of stand age, fire weather conditions, ecozone, and soil moisture class. Sites affected by historical timber harvesting had greater combustion rates due to faster carbon stock recovery rates than sites recovering from wildfire events, indicating that different boreal forest land use practices can generate divergent carbon legacy effects. We estimate the 2015 fire season in Saskatchewan emitted a total of 36.3 ± 15.0 Tg C, emphasizing the importance of southern boreal fires for regional carbon budgets. Using the southern boreal as an analog, the northern boreal may undergo fundamental shifts in forest structure and carbon dynamics, becoming dominated by stands <70 years old that hold 2-7 kg C/m2 less than current mature northern boreal stands. Our latitudinal approach reinforces previous studies showing that northern boreal stands are at a high risk of holding less carbon under changing disturbance conditions.

Climate change decreases the cooling effect from postfire albedo in boreal North America.

Fire is a primary disturbance in boreal forests and generates both positive and negative climate forcings. The influence of fire on surface albedo is a predominantly negative forcing in boreal forests, and one of the strongest overall, due to increased snow exposure in the winter and spring months. Albedo forcings are spatially and temporally heterogeneous and depend on a variety of factors related to soils, topography, climate, land cover/vegetation type, successional dynamics, time since fire, season, and fire severity. However, how these variables interact to influence albedo is not well understood, and quantifying these relationships and predicting postfire albedo becomes increasingly important as the climate changes and management frameworks evolve to consider climate impacts. Here we developed a MODIS-derived 'blue sky' albedo product and a novel machine learning modeling framework to predict fire-driven changes in albedo under historical and future climate scenarios across boreal North America. Converted to radiative forcing (RF), we estimated that fires generate an annual mean cooling of -1.77 ± 1.35 W/m2 from albedo under historical climate conditions (1971-2000) integrated over 70 years postfire. Increasing postfire albedo along a south-north climatic gradient was offset by a nearly opposite gradient in solar insolation, such that large-scale spatial patterns in RF were minimal. Our models suggest that climate change will lead to decreases in mean annual postfire albedo, and hence a decreasing strength of the negative RF, a trend dominated by decreased snow cover in spring months. Considering the range of future climate scenarios and model uncertainties, we estimate that for fires burning in the current era (2016) the cooling effect from long-term postfire albedo will be reduced by 15%-28% due to climate change.

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.

Large loss of CO2 in winter observed across the northern permafrost region.

Recent warming in the Arctic, which has been amplified during the winter1-3, greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO2)4. However, the amount of CO2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates5,6. Here we synthesize regional in situ observations of CO2 flux from arctic and boreal soils to assess current and future winter carbon losses from the northern permafrost domain. We estimate a contemporary loss of 1662 Tg C yr-1 from the permafrost region during the winter season (October through April). This loss is greater than the average growing season carbon uptake for this region estimated from process models (-1032 Tg C yr-1). Extending model predictions to warmer conditions in 2100 indicates that winter CO2 emissions will increase 17% under a moderate mitigation scenario-Representative Concentration Pathway (RCP) 4.5-and 41% under business-as-usual emissions scenario-RCP 8.5. Our results provide a new baseline for winter CO2 emissions from northern terrestrial regions and indicate that enhanced soil CO2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.