Search

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic-boreal zone (ABZ) have provided valuable information but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic-boreal CO2 fluxes (ABCflux) that aggregates in situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures), and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19% of the monthly observations), snow diffusion (3%) and eddy covariance (78%) techniques. The largest number of observations were collected during the climatological summer (June-August; 32%), and fewer observations were available for autumn (September-October; 25%), winter (December-February; 18%), and spring (March-May; 25%). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (Virkkala et al., 2021b, 10.3334/ORNLDAAC/1934). Copyright:

The continued warming of the Arctic could release vast stores of carbon into the atmosphere from high-latitude ecosystems, especially from thawing permafrost. Increasing uptake of carbon dioxide (CO2) by vegetation during longer growing seasons may partially offset such release of carbon. However, evidence of significant net annual release of carbon from site-level observations and model simulations across tundra ecosystems has been inconclusive. To address this knowledge gap, we combined top-down observations of atmospheric CO2 concentration enhancements from aircraft and a tall tower, which integrate ecosystem exchange over large regions, with bottom-up observed CO2 fluxes from tundra environments and found that the Alaska North Slope is not a consistent net source nor net sink of CO2 to the atmosphere (ranging from -6 to +6TgCyr-1 for 2012-2017). Our analysis suggests that significant biogenic CO2 fluxes from unfrozen terrestrial soils, and likely inland waters, during the early cold season (September-December) are major factors in determining the net annual carbon balance of the North Slope, implying strong sensitivity to the rapidly warming freeze-up period. At the regional level, we find no evidence of the previously reported large late-cold-season (January-April) CO2 emissions to the atmosphere during the study period. Despite the importance of the cold-season CO2 emissions to the annual total, the interannual variability in the net CO2 flux is driven by the variability in growing season fluxes. During the growing season, the regional net CO2 flux is also highly sensitive to the distribution of tundra vegetation types throughout the North Slope. This study shows that quantification and characterization of year-round CO2 fluxes from the heterogeneous terrestrial and aquatic ecosystems in the Arctic using both site-level and atmospheric observations are important to accurately project the Earth system response to future warming.

Arctic warming is affecting snow cover and soil hydrology, with consequences for carbon sequestration in tundra ecosystems. The scarcity of observations in the Arctic has limited our understanding of the impact of covarying environmental drivers on the carbon balance of tundra ecosystems. In this study, we address some of these uncertainties through a novel record of 119 site-years of summer data from eddy covariance towers representing dominant tundra vegetation types located on continuous permafrost in the Arctic. Here we found that earlier snowmelt was associated with more tundra net CO2 sequestration and higher gross primary productivity (GPP) only in June and July, but with lower net carbon sequestration and lower GPP in August. Although higher evapotranspiration (ET) can result in soil drying with the progression of the summer, we did not find significantly lower soil moisture with earlier snowmelt, nor evidence that water stress affected GPP in the late growing season. Our results suggest that the expected increased CO2 sequestration arising from Arctic warming and the associated increase in growing season length may not materialize if tundra ecosystems are not able to continue sequestering CO2 later in the season.

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.

Retrogressive thaw slumps (RTS) are a form of abrupt permafrost thaw that can rapidly mobilize ancient frozen soil carbon, magnifying the permafrost carbon feedback. However, the magnitude of this effect is uncertain, largely due to limited information about the distribution and extent of RTS across the circumpolar region. Although deep learning methods such as Convolutional Neural Networks (CNN) have shown the ability to map RTS from high-resolution satellite imagery (≤10 m), challenges remain in deploying these models across large areas. Imagery selection and procurement remain one of the largest challenges to upscaling RTS mapping projects, as the user must balance cost with resolution and sensor quality. In this study, we compared the performance of three satellite imagery sources that differed in terms of sensor quality and cost in predicting RTS using a Unet3+ CNN model and identified RTS characteristics that impact detectability. Maxar WorldView imagery was the most expensive option, with a ground sample distance of 1.85 m in the multispectral bands (downloaded at 4 m resolution). Planet Labs PlanetScope imagery was a less expensive option with a ground sample distance of approximately 3.0–4.2 m (downloaded at 3 m resolution). Although PlanetScope imagery was downloaded at a higher resolution than WorldView, the radiometric footprint is around 10–12 m, resulting in less crisp imagery. Finally, Sentinel-2 imagery is freely available and has a 10 m resolution. We used 756 RTS polygons from seven sites across Arctic Canada and Siberia in model training and 63 RTS polygons in model testing. The mean IoU of the validation and testing data sets were 0.69 and 0.75 for the WorldView model, 0.70 and 0.71 for the PlanetScope model, and 0.66 and 0.68 for the Sentinel-2 model, respectively. The IoU of the RTS class was nonlinearly related to the RTS Area, showing a strong positive correlation that attenuated as the RTS Area increased. The models were better able to predict RTS that appeared bright on a dark background and were less able to predict RTS that had higher plant cover, indicating that bare ground was a primary way the models detected RTS. Additionally, the models performed less well in wet areas or areas with patchy ground cover. These results indicate that all imagery sources tested here were able to predict larger RTS, but higher-quality imagery allows more accurate detection of smaller RTS. [ABSTRACT FROM AUTHOR]

Wetlands are the largest natural source of methane (CH4) emissions globally. Northern wetlands (>45° N), accounting for 42 % of global wetland area, are increasingly vulnerable to carbon loss, especially as CH4 emissions may accelerate under intensified high-latitude warming. However, the magnitude and spatial patterns of high-latitude CH4 emissions remain relatively uncertain. Here we present estimates of daily CH4 fluxes obtained using a new machine learning-based wetland CH4 upscaling framework (WetCH4) that applies the most complete database of eddy covariance (EC) observations available to date, and satellite remote sensing informed observations of environmental conditions at 10-km resolution. The most important predictor variables included near-surface soil temperatures (top 40 cm), vegetation reflectance, and soil moisture. Our results, modeled from 138 site-years across 26 sites, had relatively strong predictive skill with a mean R2 of 0.46 and 0.62 and a mean absolute error (MAE) of 23 nmol m-2 s-1 and 21 nmol m-2 s-1 for daily and monthly fluxes, respectively. Based on the model results, we estimated an annual average of 20.8 ±2.1 Tg CH4 yr-1 for the northern wetland region (2016–2022) and total budgets ranged from 13.7–44.1 Tg CH4 yr-1, depending on wetland map extents. Although 86 % of the estimated CH4 budget occurred during the May–October period, a considerable amount (1.4 ±0.2 Tg CH4) occurred during winter. Regionally, the West Siberian wetlands accounted for a majority (51 %) of the interannual variation in domain CH4 emissions. Significant issues with data coverage remain, with only 23 % of the sites observing year-round and most of the data from 11 wetland sites in Alaska and 10 bog/fen sites in Canada and Fennoscandia, and in general, Western Siberian Lowlands are underrepresented by EC CH4 sites. Our results provide high spatiotemporal information on the wetland emissions in the high-latitude carbon cycle and possible responses to climate change. Continued, all-season tower observations and improved soil moisture products are needed for future improvement of CH4 upscaling. The dataset can be found at https://doi.org/10.5281/zenodo.10802154 (Ying et al., 2024). [ABSTRACT FROM AUTHOR]

Significant progress in permafrost carbon science made over the past decades include the identification of vast permafrost carbon stocks, the development of new pan‐Arctic permafrost maps, an increase in terrestrial measurement sites for CO2 and methane fluxes, and important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Process‐based modeling studies now include key elements of permafrost carbon cycling and advances in statistical modeling and inverse modeling enhance understanding of permafrost region C budgets. By combining existing data syntheses and model outputs, the permafrost region is likely a wetland methane source and small terrestrial ecosystem CO2 sink with lower net CO2 uptake toward higher latitudes, excluding wildfire emissions. For 2002–2014, the strongest CO2 sink was located in western Canada (median: −52 g C m−2 y−1) and smallest sinks in Alaska, Canadian tundra, and Siberian tundra (medians: −5 to −9 g C m−2 y−1). Eurasian regions had the largest median wetland methane fluxes (16–18 g CH4 m−2 y−1). Quantifying the regional scale carbon balance remains challenging because of high spatial and temporal variability and relatively low density of observations. More accurate permafrost region carbon fluxes require: (a) the development of better maps characterizing wetlands and dynamics of vegetation and disturbances, including abrupt permafrost thaw; (b) the establishment of new year‐round CO2 and methane flux sites in underrepresented areas; and (c) improved models that better represent important permafrost carbon cycle dynamics, including non‐growing season emissions and disturbance effects. Plain Language Summary: Climate change and the consequent thawing of permafrost threatens to transform the permafrost region from a carbon sink into a carbon source, posing a challenge to global climate goals. Numerous studies over the past decades have identified important factors affecting carbon cycling, including vegetation changes, periods of soil freezing and thawing, wildfire, and other disturbance events. Overall, studies show high wetland methane emissions and a small net carbon dioxide sink strength over the terrestrial permafrost region but results differ among modeling and upscaling approaches. Continued and coordinated efforts among field, modeling, and remote sensing communities are needed to integrate new knowledge from observations to modeling and predictions and finally to policy. Key Points: Rapid warming of northern permafrost region threatens ecosystems, soil carbon stocks, and global climate targetsLong‐term observations show importance of disturbance and cold season periods but are unable to detect spatiotemporal trends in C fluxCombined modeling and syntheses show the permafrost region is a small terrestrial CO2 sink with large spatial variability and net CH4 source [ABSTRACT FROM AUTHOR]

Boreal‐Arctic regions are key stores of organic carbon (C) and play a major role in the greenhouse gas balance of high‐latitude ecosystems. The carbon‐climate (C‐climate) feedback potential of northern high‐latitude ecosystems remains poorly understood due to uncertainty in temperature and precipitation controls on carbon dioxide (CO2) uptake and the decomposition of soil C into CO2 and methane (CH4) fluxes. While CH4 fluxes account for a smaller component of the C balance, the climatic impact of CH4 outweighs CO2 (28–34 times larger global warming potential on a 100‐year scale), highlighting the need to jointly resolve the climatic sensitivities of both CO2 and CH4. Here, we jointly constrain a terrestrial biosphere model with in situ CO2 and CH4 flux observations at seven eddy covariance sites using a data‐model integration approach to resolve the integrated environmental controls on land‐atmosphere CO2 and CH4 exchanges in Alaska. Based on the combined CO2 and CH4 flux responses to climate variables, we find that 1970‐present climate trends will induce positive C‐climate feedback at all tundra sites, and negative C‐climate feedback at the boreal and shrub fen sites. The positive C‐climate feedback at the tundra sites is predominantly driven by increased CH4 emissions while the negative C‐climate feedback at the boreal site is predominantly driven by increased CO2 uptake (80% from decreased heterotrophic respiration, and 20% from increased photosynthesis). Our study demonstrates the need for joint observational constraints on CO2 and CH4 biogeochemical processes—and their associated climatic sensitivities—for resolving the sign and magnitude of high‐latitude ecosystem C‐climate feedback in the coming decades. Key Points: We use in situ measurements to constrain the modeled joint climatic sensitivity of land‐atmosphere CH4 and CO2 exchangesA continued 1970‐present climate trend leads to positive C‐climate feedback in wet tundra sites but negative feedback in boreal and shrub fen sitesCH4 respiration dominates the positive tundra site feedback, CO2 respiration dominates the negative boreal and shrub fen sites feedback [ABSTRACT FROM AUTHOR]

Small water bodies (i.e., ponds; <0.01 km2) play an important role in Earth System processes, including carbon cycling and emissions of methane. Detection and monitoring of ponds using satellite imagery has been extremely difficult and many water maps are biased toward lakes (>0.01 km2). We leverage high‐resolution (3 m) optical satellite imagery from Planet Labs and deep learning methods to map seasonal changes in pond and lake areal extent across four regions in Alaska. Our water maps indicate that changes in open water extent over the snow‐free season are especially pronounced in ponds. To investigate potential impacts of seasonal changes in pond area on carbon emissions, we provide a case study of open water methane emission budgets using the new water maps. Our approach has widespread applications for water resources, habitat and land cover change assessments, wildlife management, risk assessments, and other biogeochemical modeling efforts. Plain Language Summary: Small water bodies (<0.01 km2) are an important driver of many Earth system processes. Despite their importance, many existing water mapping products have difficulty detecting these small water features and their seasonal changes in surface area. We used deep learning and high‐resolution (3 m) satellite imagery to map and monitor seasonal changes in the areal extent of lakes and small ponds across four regions in Alaska. The resulting water maps accounted for considerably more water coverage than existing products. The maps also effectively tracked widespread seasonal changes in pond and lake area that were not previously identified. This demonstrates the importance of monitoring surface water at high spatial resolutions and across seasons. Key Points: Deep learning and 3 m resolution satellite imagery from Planet Labs can detect and track ponds and lakes >0.0001 km2Total surface area for ponds (<0.01 km2) in boreal forest and tundra environments can vary by 20%–40% throughout an individual seasonPonds can contribute to a broad range (8%–37%) of total methane emissions from lakes and ponds in northern boreal forest and tundra [ABSTRACT FROM AUTHOR]

Arctic‐boreal landscapes are experiencing profound warming, along with changes in ecosystem moisture status and disturbance from fire. This region is of global importance in terms of carbon feedbacks to climate, yet the sign (sink or source) and magnitude of the Arctic‐boreal carbon budget within recent years remains highly uncertain. Here, we provide new estimates of recent (2003–2015) vegetation gross primary productivity (GPP), ecosystem respiration (Reco), net ecosystem CO2 exchange (NEE; Reco − GPP), and terrestrial methane (CH4) emissions for the Arctic‐boreal zone using a satellite data‐driven process‐model for northern ecosystems (TCFM‐Arctic), calibrated and evaluated using measurements from >60 tower eddy covariance (EC) sites. We used TCFM‐Arctic to obtain daily 1‐km2 flux estimates and annual carbon budgets for the pan‐Arctic‐boreal region. Across the domain, the model indicated an overall average NEE sink of −850 Tg CO2‐C year−1. Eurasian boreal zones, especially those in Siberia, contributed to a majority of the net sink. In contrast, the tundra biome was relatively carbon neutral (ranging from small sink to source). Regional CH4 emissions from tundra and boreal wetlands (not accounting for aquatic CH4) were estimated at 35 Tg CH4‐C year−1. Accounting for additional emissions from open water aquatic bodies and from fire, using available estimates from the literature, reduced the total regional NEE sink by 21% and shifted many far northern tundra landscapes, and some boreal forests, to a net carbon source. This assessment, based on in situ observations and models, improves our understanding of the high‐latitude carbon status and also indicates a continued need for integrated site‐to‐regional assessments to monitor the vulnerability of these ecosystems to climate change. [ABSTRACT FROM AUTHOR]

Long‐term atmospheric CO2 concentration records have suggested a reduction in the positive effect of warming on high‐latitude carbon uptake since the 1990s. A variety of mechanisms have been proposed to explain the reduced net carbon sink of northern ecosystems with increased air temperature, including water stress on vegetation and increased respiration over recent decades. However, the lack of consistent long‐term carbon flux and in situ soil moisture data has severely limited our ability to identify the mechanisms responsible for the recent reduced carbon sink strength. In this study, we used a record of nearly 100 site‐years of eddy covariance data from 11 continuous permafrost tundra sites distributed across the circumpolar Arctic to test the temperature (expressed as growing degree days, GDD) responses of gross primary production (GPP), net ecosystem exchange (NEE), and ecosystem respiration (ER) at different periods of the summer (early, peak, and late summer) including dominant tundra vegetation classes (graminoids and mosses, and shrubs). We further tested GPP, NEE, and ER relationships with soil moisture and vapor pressure deficit to identify potential moisture limitations on plant productivity and net carbon exchange. Our results show a decrease in GPP with rising GDD during the peak summer (July) for both vegetation classes, and a significant relationship between the peak summer GPP and soil moisture after statistically controlling for GDD in a partial correlation analysis. These results suggest that tundra ecosystems might not benefit from increased temperature as much as suggested by several terrestrial biosphere models, if decreased soil moisture limits the peak summer plant productivity, reducing the ability of these ecosystems to sequester carbon during the summer. [ABSTRACT FROM AUTHOR]

Background: High-resolution soil moisture estimates are critical for planning water management and assessing environmental quality. In-situ measurements alone are too costly to support the spatial and temporal resolutions needed for water management. Recent efforts have combined calibration data with machine learning algorithms to fill the gap where high resolution moisture estimates are lacking at the field scale. This study aimed to provide calibrated soil moisture models and methodology for generating gridded estimates of soil moisture at multiple depths, according to user-defined temporal periods, spatial resolution and extent. Methods: We applied nearly one million national library soil moisture records from over 100 sites, spanning the U.S. Midwest and West, to build Quantile Random Forest (QRF) calibration models. The QRF models were built on covariates including soil moisture estimates from North American Land Data Assimilation System (NLDAS), soil properties, climate variables, digital elevation models, and remote sensing-derived indices. We also explored an alternative approach that adopted a regionalized calibration dataset for the Western U.S. The broad-scale QRF models were independently validated according to sampling depths, land cover type, and observation period. We then explored the model performance improved with local samples used for spiking. Finally, the QRF models were applied to estimate soil moisture at the field scale where evaluation was carried out to check estimated temporal and spatial patterns. Results: The broad-scale QRF model showed moderate performance (R2 = 0.53, RMSE = 0.078 m3/m3) when data points from all depth layers (up to 100 cm) were considered for an independent validation. Elevation, NLDAS-derived moisture, soil properties, and sampling depth were ranked as the most important covariates. The best model performance was observed for forest and pasture sites (R2 > 0.5; RMSE < 0.09 m3/m3), followed by grassland and cropland (R2 > 0.4; RMSE < 0.11 m3/m3). Model performance decreased with sampling depths and was slightly lower during the winter months. Spiking the national QRF model with local samples improved model performance by reducing the RMSE to less than 0.05 m3/m3 for grassland sites. At the field scale, model estimates illustrated more accurate temporal trends for surface than subsurface soil layers. Model estimated spatial patterns need to be further improved and validated with management data. Conclusions: The model accuracy for top 0-20 cm soil depth (R2 > 0.5, RMSE < 0.08 m3/m3) showed promise for adopting the methodology for soil moisture monitoring. The success of spiking the national model with local samples showed the need to collect multi-year high frequency (e.g., hourly) sensor-based field measurements to improve estimates of soil moisture for a longer time period. Future work should improve model performance for deeper depths with additional hydraulic properties and use of locally-selected calibration datasets. [ABSTRACT FROM AUTHOR]

A capability for mapping meter-level resolution soil moisture with frequent temporal sampling over large regions is essential for quantifying local-scale environmental heterogeneity and eco-hydrologic behavior. However, available surface soil moisture (SSM) products generally involve much coarser grain sizes ranging from 30 m to several 10 s of kilometers. Hence, a new method is proposed to estimate 3-m resolution SSM using a combination of multi-sensor fusion, machine-learning (ML), and Cumulative Distribution Function (CDF) matching approaches. This method established favorable SSM correspondence between 3-m pixels and overlying 9-km grid cells from overlapping Planet SuperDove (PSD) observations and NASA Soil Moisture Active-Passive (SMAP) mission products. The resulting 3-m SSM predictions showed improved accuracy by reducing absolute bias and RMSE by ~0.01 cm3/cm3 over the original SMAP data in relation to in situ soil moisture measurements for the Australian Yanco region while preserving the high sampling frequency (1–3 day global revisit) and sensitivity to surface wetness (R 0.865) from SMAP. Heterogeneous soil moisture distributions varying with vegetation biomass gradients and irrigation regimes were generally captured within a selected study area. Further algorithm refinement and implementation for regional applications will allow for improvement in water resources management, precision agriculture, and disaster forecasts and responses. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Past efforts to synthesize and quantify the magnitude and change in carbon dioxide (CO2) fluxes in terrestrial ecosystems across the rapidly warming Arctic-Boreal Zone (ABZ) have provided valuable information, but were limited in their geographical and temporal coverage. Furthermore, these efforts have been based on data aggregated over varying time periods, often with only minimal site ancillary data, thus limiting their potential to be used in large-scale carbon budget assessments. To bridge these gaps, we developed a standardized monthly database of Arctic-Boreal CO2 fluxes (ABCflux) that aggregates in-situ measurements of terrestrial net ecosystem CO2 exchange and its derived partitioned component fluxes: gross primary productivity and ecosystem respiration. The data span from 1989 to 2020 with over 70 supporting variables that describe key site conditions (e.g., vegetation and disturbance type), micrometeorological and environmental measurements (e.g., air and soil temperatures) and flux measurement techniques. Here, we describe these variables, the spatial and temporal distribution of observations, the main strengths and limitations of the database, and the potential research opportunities it enables. In total, ABCflux includes 244 sites and 6309 monthly observations; 136 sites and 2217 monthly observations represent tundra, and 108 sites and 4092 observations represent the boreal biome. The database includes fluxes estimated with chamber (19 % of the monthly observations), snow diffusion (3 %) and eddy covariance (78 %) techniques. The largest number of observations were collected during the climatological summer (June-August; 32 %), and fewer observations were available for autumn (September-October; 25 %), winter (December-February; 18 %), and spring (March-May; 25 %). ABCflux can be used in a wide array of empirical, remote sensing and modeling studies to improve understanding of the regional and temporal variability in CO2 fluxes, and to better estimate the terrestrial ABZ CO2 budget. ABCflux is openly and freely available online (https://doi.org/10.3334/ORNLDAAC/1934, Virkkala et al., 2021a). [ABSTRACT FROM AUTHOR]

The increase in wildfire occurrence and severity seen over the past decades in the boreal and Arctic biomes is expected to continue in the future in response to rapid climate change in this region. Recent studies documented positive trends in gross primary productivity (GPP) for Arctic boreal biomes driven by warming, but it is unclear how GPP trends are affected by wildfires. Here, we used satellite vegetation observations and environmental data with a diagnostic GPP model to analyze recovery from large fires in Alaska over the period 2000–2019. We confirmed earlier findings that warmer‐than‐average years provide favorable climate conditions for vegetation growth, leading to a GPP increase of 1 Tg C yr−1, contributed mainly from enhanced productivity in the early growing season. However, higher temperatures increase the risk of wildfire occurrence leading to direct carbon loss over a period of 1–3 years. While mortality related to severe wildfires reduce ecosystem productivity, post‐fire productivity in moderately burned areas shows a significant positive trend. The rapid GPP recovery following fires reported here might be favorable for maintaining the region's net carbon sink, but wildfires can indirectly promote the release of long‐term stored carbon in the permafrost. With the projected increase in severity and frequency of wildfires in the future, we expect a reduction of GPP and therefore amplification of climate warming in this region. Plain Language Summary: We analyzed the indirect impact of wildfires in Alaska on the amount of carbon fixed during photosynthesis, which is known as ecosystem gross primary productivity (GPP). Satellite observations and our remote sensing GPP model confirmed earlier findings of increasing vegetation growth in warmer‐than‐average years. On the other hand, risk of wildfire occurrence is higher during the warmer years, which leads to direct vegetation loss, and as result reduction in ecosystem carbon uptake over a period of 1–3 years. We also observed rapid GPP recovery following fires, but severe wildfires are expected to reduce the GPP and therefore amplification of climate warming in this region due to the release of long‐term stored carbon in the permafrost. Key Points: Gross primary productivity (GPP) data in Alaska show increasing trends in the last two decades due to increased carbon uptake in the early growing seasonTemperature enhances plant productivity in the early growing season and increases the risk of wildfire occurrencesEcosystem GPP is negatively affected by wildfires, but shows fast recovery following less severe wildfires [ABSTRACT FROM AUTHOR]

Soil organic carbon (SOC) is an important metric of soil health and the terrestrial carbon balance. Short‐term climate variations affect SOC through changes in temperature and moisture, which control vegetation growth and soil decomposition. We evaluated a satellite data‐driven carbon model, operating under the NASA Soil Moisture Active‐Passive (SMAP) mission, as a means of monitoring global surface SOC dynamics. The SMAP Level 4 Carbon (L4C) product estimates a daily global carbon budget including surface (0‐ to 5‐cm depth) SOC. We found that the L4C mean latitudinal SOC distribution is generally consistent with alternative assessments from static soil inventory records and dynamic global vegetation models (r ≥ 0.89). Within forest systems, based on inventory data, L4C SOC is most similar in magnitude to litterfall but is correlated with coarse woody debris (r=0.86) and total SOC (r=0.81). L4C SOC is sensitive to seasonal and annual climate variability, with mean residence times that range from 1.5 years in the wet tropics to 17 years in the cold tundra. Incorporating soil moisture retrievals from the SMAP L‐band (1.4 GHz) microwave radiometer within the L4C algorithm provides enhanced soil moisture sensitivity under low‐to‐moderate vegetation cover (<5 kg m−2 vegetation water content). The L‐band soil moisture had the greatest impact on the L4C carbon budget in semiarid regions, which span almost 60% of the globe and account for substantial variability in the terrestrial carbon sink. The L4C operational product enables prognostic investigations into effects of recent climate trends and anomalies (e.g., droughts and pluvials) on shallow soil carbon dynamics. Plain Language Summary: Healthy soils are important for agriculture, climate‐change adaptation, and biodiversity. The most common way that scientists describe soil health is by taking a sample and measuring the amount of soil organic carbon, which is carbon stored in living things like microbes and fungi or that came from once‐living things like trees and other plants and is now in the soil. Collecting these samples requires a lot of time and work, which makes it very difficult for scientists to describe the soil health of large areas like countries or continents. Earth‐orbiting satellites that measure surface conditions like the temperature and amount of moisture in the soil have been used to inform computer models of how plants and soil change in response to a changing climate. Because satellites can see the entire globe, the computer models they inform have global coverage. We used one such computer model to estimate the amount of soil organic carbon stored in the world's soils. We compared our estimates of soil organic carbon to estimates from other methods and found that they agree very well once we accounted for the different soil depths used. What is particularly new and exciting about using this computer model for this purpose is that we can estimate seasonal and annual changes in soil organic carbon. This allows us to estimate how much soil health is impacted by short‐term natural events like droughts or floods. Key Points: SMAP L4C provides an estimated daily, global terrestrial carbon budget including surface (0‐ to 5‐cm depth) soil organic carbon (SOC)L4C SOC is consistent with other global benchmarks including soil inventory records and dynamic global vegetation models (r > 0.89)SMAP L‐band soil moisture has greatest influence on SOC variability in semiarid lands, which encompass about 60% of the global domain [ABSTRACT FROM AUTHOR]

The contribution of soil heterotrophic respiration to the boreal–Arctic carbon (CO 2) cycle and its potential feedback to climate change remains poorly quantified. We developed a remote-sensing-driven permafrost carbon model at intermediate scale (∼1 km) to investigate how environmental factors affect the magnitude and seasonality of soil heterotrophic respiration in Alaska. The permafrost carbon model simulates snow and soil thermal dynamics and accounts for vertical soil carbon transport and decomposition at depths up to 3 m below the surface. Model outputs include soil temperature profiles and carbon fluxes at 1 km resolution spanning the recent satellite era (2001–2017) across Alaska. Comparisons with eddy covariance tower measurements show that the model captures the seasonality of carbon fluxes, with favorable accuracy in simulating net ecosystem CO 2 exchange (NEE) for both tundra (R>0.8 , root mean square error (RMSE – 0.34 g C m -2 d -1), and boreal forest (R>0.73 ; RMSE – 0.51 g C m -2 d -1). Benchmark assessments using two regional in situ data sets indicate that the model captures the complex influence of snow insulation on soil temperature and the temperature sensitivity of cold-season soil heterotrophic respiration. Across Alaska, we find that seasonal snow cover imposes strong controls on the contribution from different soil depths to total soil heterotrophic respiration. Earlier snowmelt in spring promotes deeper soil warming and enhances the contribution of deeper soils to total soil heterotrophic respiration during the later growing season, thereby reducing net ecosystem carbon uptake. Early cold-season soil heterotrophic respiration is closely linked to the number of snow-free days after the land surface freezes (R=-0.48 , p<0.1), i.e., the delay in snow onset relative to surface freeze onset. Recent trends toward earlier autumn snow onset in northern Alaska promote a longer zero-curtain period and enhanced cold-season respiration. In contrast, southwestern Alaska shows a strong reduction in the number of snow-free days after land surface freeze onset, leading to earlier soil freezing and a large reduction in cold-season soil heterotrophic respiration. Our results also show nonnegligible influences of subgrid variability in surface conditions on the model-simulated CO 2 seasonal cycle, especially during the early cold season at 10 km scale. Our results demonstrate the critical role of snow cover affecting the seasonality of soil temperature and respiration and highlight the challenges of incorporating these complex processes into future projections of the boreal–Arctic carbon cycle. [ABSTRACT FROM AUTHOR]

The contribution of soil heterotrophic respiration to the boreal-Arctic carbon (CO2) cycle and its potential feedback to climate change remain poorly quantified. We developed a remote sensing driven permafrost carbon model at intermediate scale (~ 1 km) to investigate how environmental factors affect the magnitude and seasonality of soil heterotrophic respiration in Alaska. The permafrost carbon model simulates snow and soil thermal dynamics, and accounts for vertical soil carbon transport and decomposition at depths up to 3 m below surface. Model outputs include soil temperature profiles and carbon fluxes at 1-km resolution spanning the recent satellite era (2001-2017) across Alaska. Comparisons with eddy covariance tower measurements show that the model captures the seasonality of carbon fluxes, with favorable accuracy in predicting net ecosystem CO2 exchange (NEE) in both tundra (R > 0.8, RMSE = 0.34 g C m-2 d-1) and boreal forest (R > 0.73, RMSE = 0.51 g C m-2 d-1). Benchmark assessments using two regional in-situ datasets indicate that the model captures the complex influence of snow insulation on soil temperature, and the temperature sensitivity of cold-season soil respiration. Across Alaska, we find that seasonal snow cover imposes strong controls on the contribution from different soil depths to total soil carbon emissions. Earlier snow melt in spring promotes deeper soil warming and enhances the contribution of deeper soils to total soil respiration during the later growing season, thereby reducing net ecosystem carbon uptake. Early cold-season soil respiration is closely linked to the number of snow-free days after land surface freezes (R = -0.48, p < 0.1), i.e. the delay in snow onset relative to surface freeze onset. Recent trends toward earlier autumn snow onset in northern Alaska promote a longer zero-curtain period and enhanced cold-season respiration. In contrast, southwestern Alaska shows a strong reduction in the number of snow-free days after land surface freeze onset, leading to earlier soil freezing and a large reduction in cold-season soil respiration. Our results also show non-negligible influences of sub-grid variability of surface conditions on the model simulated CO2 seasonal cycle, especially during the early cold season at 10-km scale. Our results demonstrate the critical role of snow cover affecting the seasonality of soil temperature and respiration and highlight the challenges of incorporating these complex processes into future projections of boreal-Arctic carbon cycle. [ABSTRACT FROM AUTHOR]

Arctic and boreal ecosystems play an important role in the global carbon (C) budget, and whether they act as a future net C sink or source depends on climate and environmental change. Here, we used complementary in situ measurements, model simulations, and satellite observations to investigate the net carbon dioxide (CO2) seasonal cycle and its climatic and environmental controls across Alaska and northwestern Canada during the anomalously warm winter to spring conditions of 2015 and 2016 (relative to 2010–2014). In the warm spring, we found that photosynthesis was enhanced more than respiration, leading to greater CO2 uptake. However, photosynthetic enhancement from spring warming was partially offset by greater ecosystem respiration during the preceding anomalously warm winter, resulting in nearly neutral effects on the annual net CO2 balance. Eddy covariance CO2 flux measurements showed that air temperature has a primary influence on net CO2 exchange in winter and spring, while soil moisture has a primary control on net CO2 exchange in the fall. The net CO2 exchange was generally more moisture limited in the boreal region than in the Arctic tundra. Our analysis indicates complex seasonal interactions of underlying C cycle processes in response to changing climate and hydrology that may not manifest in changes in net annual CO2 exchange. Therefore, a better understanding of the seasonal response of C cycle processes may provide important insights for predicting future carbon–climate feedbacks and their consequences on atmospheric CO2 dynamics in the northern high latitudes. [ABSTRACT FROM AUTHOR]

A new automated method enabling consistent satellite assessment of seasonal lake ice phenology at 5 km resolution was developed for all lake pixels (water coverage ≥90 %) in the Northern Hemisphere using 36.5 GHz H-polarized brightness temperature (Tb) observations from the Advanced Microwave Scanning Radiometer for EOS and Advanced Microwave Scanning Radiometer 2 (AMSR-E/2) sensors. The lake phenology metrics include seasonal timing and duration of annual ice cover. A moving t test (MTT) algorithm allows for automated lake ice retrievals with daily temporal fidelity and 5 km resolution gridding. The resulting ice phenology record shows strong agreement with available ground-based observations from the Global Lake and River Ice Phenology Database (95.4% temporal agreement) and favorable correlations (R) with alternative ice phenology records from the Interactive Multisensor Snow and Ice Mapping System (R = 0:84 for water clear of ice (WCI) dates; R = 0:41 for complete freeze over (CFO) dates) and Canadian Ice Service (R = 0:86 for WCI dates; R = 0:69 for CFO dates). Analysis of the resulting 12-year (2002-2015) AMSR-E/2 ice record indicates increasingly shorter ice cover duration for 43 out of 71 (60.6 %) Northern Hemisphere lakes examined, with significant (p <0.05) regional trends toward earlier ice melting for only five lakes. Higherlatitude lakes reveal more widespread and larger trends toward shorter ice cover duration than lower-latitude lakes, consistent with enhanced polar warming. This study documents a new satellite-based approach for rapid assessment and regional monitoring of seasonal ice cover changes over large lakes, with resulting accuracy suitable for global change studies. [ABSTRACT FROM AUTHOR]

The Arctic is currently undergoing intense changes in climate; vegetation composition and productivity are expected to respond to such changes. To understand the impacts of climate change on the function of Arctic tundra ecosystems within the global carbon cycle, it is crucial to improve the understanding of vegetation distribution and heterogeneity at multiple scales. Information detailing the fine-scale spatial distribution of tundra communities provided by high resolution vegetation mapping, is needed to understand the relative contributions of and relationships between single vegetation community measurements of greenhouse gas fluxes (e.g., ∼1 m chamber flux) and those encompassing multiple vegetation communities (e.g., ∼300 m eddy covariance measurements). The objectives of this study were: (1) to determine whether dominant Arctic tundra vegetation communities found in different locations are spectrally distinct and distinguishable using field spectroscopy methods; and (2) to test which combination of raw reflectance and vegetation indices retrieved from field and satellite data resulted in accurate vegetation maps and whether these were transferable across locations to develop a systematic method to map dominant vegetation communities within larger eddy covariance tower footprints distributed along a 300 km transect in northern Alaska. We showed vegetation community separability primarily in the 450-510 nm, 630-690 nm and 705-745 nm regions of the spectrum with the field spectroscopy data. This is line with the different traits of these arctic tundra communities, with the drier, often non-vascular plant dominated communities having much higher reflectance in the 450-510 nm and 630-690 nm regions due to the lack of photosynthetic material, whereas the low reflectance values of the vascular plant dominated communities highlight the strong light absorption found here. High classification accuracies of 92% to 96% were achieved using linear discriminant analysis with raw and rescaled spectroscopy reflectance data and derived vegetation indices. However, lower classification accuracies (∼70%) resulted when using the coarser 2.0 m WorldView-2 data inputs. The results from this study suggest that tundra vegetation communities are separable using plot-level spectroscopy with hand-held sensors. These results also show that tundra vegetation mapping can be scaled from the plot level (<1 m) to patch level (<500 m) using spectroscopy data rescaled to match the wavebands of the multispectral satellite remote sensing. We find that developing a consistent method for classification of vegetation communities across the flux tower sites is a challenging process, given the spatial variability in vegetation communities and the need for detailed vegetation survey data for training and validating classification algorithms. This study highlights the benefits of using fine-scale field spectroscopy measurements to obtain tundra vegetation classifications for landscape analyses and use in carbon flux scaling studies. Improved understanding of tundra vegetation distributions will also provide necessary insight into the ecological processes driving plant community assemblages in Arctic environments. [ABSTRACT FROM AUTHOR]

Abstract: Conservation tillage management has been advocated for carbon sequestration and soil quality preservation purposes. Past satellite image analyses have had difficulty in differentiating between no-till (NT) and minimal tillage (MT) conservation classes due to similarities in surface residues, and may have been restricted by the availability of cloud-free satellite imagery. This study hypothesized that the inclusion of high temporal data into the classification process would increase conservation tillage accuracy due to the added likelihood of capturing spectral changes in MT fields following a tillage disturbance. Classification accuracies were evaluated for Random Forest models based on 250-m and 500-m MODIS, 30-m Landsat, and 30-m synthetic reflectance values. Synthetic (30-m) data derived from the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) were evaluated because high frequency Landsat image sets are often unavailable within a cropping season due to cloud issues. Classification results from a five-date Landsat model were substantially better than those reported by previous classification tillage studies, with 94% total and ≥88% class producer''s accuracies. Landsat-derived models based on individual image scenes (May through August) yielded poor MT classifications, but a monthly increase in accuracy illustrated the importance of temporal sampling for capturing regional tillage disturbance signatures. MODIS-based model accuracies (90% total; ≥82% class) were lower than in the five-date Landsat model, but were higher than previous image-based and survey-based tillage classification results. Almost all the STARFM prediction-based models had classification accuracies higher than, or comparable to, the MODIS-based results (>90% total; ≥84% class) but the resulting model accuracies were dependent on the MODIS/Landsat base pairs used to generate the STARFM predictions. Also evident within the STARFM prediction-based models was the ability for high frequency data series to compensate for degraded synthetic spectral values when classifying field-based tillage. The decision to use MODIS or STARFM-based data within conservation tillage analysis is likely situation dependent. A MODIS-based approach requires little data processing and could be more efficient for large-area mapping; however a STARFM-based analysis might be more appropriate in mixed-pixel situations that could potentially compromise classification accuracy. [Copyright &y& Elsevier]

Abstract: We used an object-oriented approach in conjunction with the Random Forest algorithm to classify agricultural practices, including tillage (till or no-till (NT)), crop intensity, and grassland-based conservation reserve (CR). The object-oriented approach allowed for per-field classifications and the incorporation of contextual elements in addition to spectral features. Random Forest is a classification tree-based advanced classifier that avoids data over-fitting associated with many tree-based models and incorporates an unbiased internal classification accuracy assessment. Landsat satellite imagery was chosen for its continuous coverage, cost effectiveness, and image accessibility. Classification results for 2007 included producer''s accuracies of 91% for NT and 31% for tillage when applying Random Forest to image objects generated from a May Landsat image. Low classification accuracies likely were attributed to the misclassification of conservation-based tillage practices as NT. Results showed that the binary separation of tillage from NT management is likely not appropriate due to surface spectral and textural similarities between NT and conservation-type tillage practices. Crop and CR lands resulted in producer''s accuracies of 100% and 90%, respectively. Crop and fallow producer''s accuracies were 95% and 82% in the 2007 classification, despite post-senesced vegetation; misclassification within the fallow class was attributed to pixel-mixing problems in areas of narrow (<100 m) strip management. A between-date normalized difference vegetation index approach was successfully used to detect areas having “changed” in vegetation status between the 2007 and prior image dates; classified “changed” objects were then merged with “unchanged” objects to produce crop status maps. Field crop intensity was then determined from the multi-year analysis of generated crop status maps. [Copyright &y& Elsevier]

For assessing the performance of satellite products and algorithms, SIC data sets were derived from ship-borne photographic observations acquired along cruise paths and compared with six passive microwave remote sensing products. Frozen soil - One of the key issues in satellite microwave sensing of frozen soil is the determination of microwave radiation response depth (MRRD). Multi-source data fusion approaches, emerging techniques such as microsatellites and artificial intelligence, light detection and ranging (LIDAR) and structure from motion photogrammetry, and next generation satellite missions will enable unprecedented remote sensing performance in cold land studies [[10]]. [Extracted from the article]

Cold regions, including high-latitude and high-altitude landscapes, are experiencing profound environmental changes driven by global warming. With the advance of earth observation technology, remote sensing has become increasingly important for detecting, monitoring, and understanding environmental changes over vast and remote regions. This paper provides an overview of recent achievements, challenges, and opportunities for land remote sensing of cold regions by (a) summarizing the physical principles and methods in remote sensing of selected key variables related to ice, snow, permafrost, water bodies, and vegetation; (b) highlighting recent environmental nonstationarity occurring in the Arctic, Tibetan Plateau, and Antarctica as detected from satellite observations; (c) discussing the limits of available remote sensing data and approaches for regional monitoring; and (d) exploring new opportunities from next-generation satellite missions and emerging methods for accurate, timely, and multi-scale mapping of cold regions. [ABSTRACT FROM AUTHOR]

Near-surface atmospheric Vapor Pressure Deficit (VPD) is a key environmental variable affecting vegetation water stress, evapotranspiration, and atmospheric moisture demand. Although VPD is readily derived from in situ standard weather station measurements, more spatially continuous global observations for regional monitoring of VPD are lacking. Here, we document a new method to estimate daily (both a.m. and p.m.) global land surface VPD at a 25-km resolution using a satellite passive microwave remotely sensed Land Parameter Data Record (LPDR) derived from the Advanced Microwave Scanning Radiometer (AMSR) sensors. The AMSR-derived VPD record shows strong correspondence (correlation coefficient ≥ 0.80, p-value < 0.001) and overall good performance (0.48 kPa ≤ Root Mean Square Error ≤ 0.69 kPa) against independent VPD observations from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data. The estimated AMSR VPD retrieval uncertainties vary with land cover type, satellite observation time, and underlying LPDR data quality. These results provide new satellite capabilities for global mapping and monitoring of land surface VPD dynamics from ongoing AMSR2 operations. Overall good accuracy and similar observations from both AMSR2 and AMSR-E allow for the development of climate data records documenting recent (from 2002) VPD trends and potential impacts on vegetation, land surface evaporation, and energy budgets. [ABSTRACT FROM AUTHOR]

Background: High-resolution soil moisture estimates are critical for planning water management and assessing environmental quality. In-situ measurements alone are too costly to support the spatial and temporal resolutions needed for water management. Recent efforts have combined calibration data with machine learning algorithms to fill the gap where high resolution moisture estimates are lacking at the field scale. This study aimed to provide calibrated soil moisture models and methodology for generating gridded estimates of soil moisture at multiple depths, according to user-defined temporal periods, spatial resolution and extent., Methods: We applied nearly one million national library soil moisture records from over 100 sites, spanning the U.S. Midwest and West, to build Quantile Random Forest (QRF) calibration models. The QRF models were built on covariates including soil moisture estimates from North American Land Data Assimilation System (NLDAS), soil properties, climate variables, digital elevation models, and remote sensing-derived indices. We also explored an alternative approach that adopted a regionalized calibration dataset for the Western U.S. The broad-scale QRF models were independently validated according to sampling depths, land cover type, and observation period. We then explored the model performance improved with local samples used for spiking. Finally, the QRF models were applied to estimate soil moisture at the field scale where evaluation was carried out to check estimated temporal and spatial patterns., Results: The broad-scale QRF model showed moderate performance (R 2 = 0.53, RMSE = 0.078 m 3 /m 3 ) when data points from all depth layers (up to 100 cm) were considered for an independent validation. Elevation, NLDAS-derived moisture, soil properties, and sampling depth were ranked as the most important covariates. The best model performance was observed for forest and pasture sites (R 2 > 0.5; RMSE < 0.09 m 3 /m 3 ), followed by grassland and cropland (R 2 > 0.4; RMSE < 0.11 m 3 /m 3 ). Model performance decreased with sampling depths and was slightly lower during the winter months. Spiking the national QRF model with local samples improved model performance by reducing the RMSE to less than 0.05 m 3 /m 3 for grassland sites. At the field scale, model estimates illustrated more accurate temporal trends for surface than subsurface soil layers. Model estimated spatial patterns need to be further improved and validated with management data., Conclusions: The model accuracy for top 0-20 cm soil depth (R 2 > 0.5, RMSE < 0.08 m 3 /m 3 ) showed promise for adopting the methodology for soil moisture monitoring. The success of spiking the national model with local samples showed the need to collect multi-year high frequency ( e.g. , hourly) sensor-based field measurements to improve estimates of soil moisture for a longer time period. Future work should improve model performance for deeper depths with additional hydraulic properties and use of locally-selected calibration datasets., Competing Interests: The authors declare that they have no competing interests., (© 2022 Xia et al.)

The regional variability in tundra and boreal carbon dioxide (CO 2 ) 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) CO 2 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 CO 2 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 CO 2 fluxes and test the accuracy and uncertainty of different statistical models. CO 2 fluxes were upscaled at relatively high spatial resolution (1 km 2 ) 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 CO 2 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 CO 2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO 2 sink during 1990-2015, although uncertainty remains high., (© 2021 John Wiley & Sons Ltd.)

Recent warming in the Arctic, which has been amplified during the winter 1-3 , greatly enhances microbial decomposition of soil organic matter and subsequent release of carbon dioxide (CO 2 ) 4 . However, the amount of CO 2 released in winter is highly uncertain and has not been well represented by ecosystem models or by empirically-based estimates 5,6 . Here we synthesize regional in situ observations of CO 2 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 CO 2 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 CO 2 emissions from northern terrestrial regions and indicate that enhanced soil CO 2 loss due to winter warming may offset growing season carbon uptake under future climatic conditions.