Carbon uptake in Eurasian boreal forests dominates the high-latitude net ecosystem carbon budget.

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.

Pan-Arctic soil moisture control on tundra carbon sequestration and plant productivity.

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.

Machine learning based estimation of field-scale daily, high resolution, multi-depth soil moisture for the Western and Midwestern United States.

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.

Respiratory loss during late-growing season determines the net carbon dioxide sink in northern permafrost regions.

Warming of northern high latitude regions (NHL, > 50 °N) has increased both photosynthesis and respiration which results in considerable uncertainty regarding the net carbon dioxide (CO2) balance of NHL ecosystems. Using estimates constrained from atmospheric observations from 1980 to 2017, we find that the increasing trends of net CO2 uptake in the early-growing season are of similar magnitude across the tree cover gradient in the NHL. However, the trend of respiratory CO2 loss during late-growing season increases significantly with increasing tree cover, offsetting a larger fraction of photosynthetic CO2 uptake, and thus resulting in a slower rate of increasing annual net CO2 uptake in areas with higher tree cover, especially in central and southern boreal forest regions. The magnitude of this seasonal compensation effect explains the difference in net CO2 uptake trends along the NHL vegetation- permafrost gradient. Such seasonal compensation dynamics are not captured by dynamic global vegetation models, which simulate weaker respiration control on carbon exchange during the late-growing season, and thus calls into question projections of increasing net CO2 uptake as high latitude ecosystems respond to warming climate conditions.

Earlier snowmelt may lead to late season declines in plant productivity and carbon sequestration in Arctic tundra ecosystems.

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.

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.

Increased high-latitude photosynthetic carbon gain offset by respiration carbon loss during an anomalous warm winter to spring transition.

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.

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.

Global Satellite Retrievals of the Near-Surface Atmospheric Vapor Pressure Deficit from AMSR-E and AMSR2.

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.

Assessing global surface water inundation dynamics using combined satellite information from SMAP, AMSR2 and Landsat.

A method to assess global land surface water (fw) inundation dynamics was developed by exploiting the enhanced fw sensitivity of L-band (1.4 GHz) passive microwave observations from the Soil Moisture Active Passive (SMAP) mission. The L-band fw (fwLBand ) retrievals were derived using SMAP H-polarization brightness temperature (Tb ) observations and predefined L-band reference microwave emissivities for water and land endmembers. Potential soil moisture and vegetation contributions to the microwave signal were represented from overlapping higher frequency Tb observations from AMSR2. The resulting fwLBand global record has high temporal sampling (1-3 days) and 36-km spatial resolution. The fwLBand annual averages corresponded favourably (R=0.85, p-value<0.001) with a 250-m resolution static global water map (MOD44W) aggregated at the same spatial scale, while capturing significant inundation variations worldwide. The monthly fwLBand averages also showed seasonal inundation changes consistent with river discharge records within six major US river basins. An uncertainty analysis indicated generally reliable fwLBand performance for major land cover areas and under low to moderate vegetation cover, but with lower accuracy for detecting water bodies covered by dense vegetation. Finer resolution (30-m) fwLBand results were obtained for three sub-regions in North America using an empirical downscaling approach and ancillary global Water Occurrence Dataset (WOD) derived from the historical Landsat record. The resulting 30-m fwLBand retrievals showed favourable spatial accuracy for water (commission error 31.46%, omission error 30.20%) and land (commission error 0.87%, omission error 0.96%) classifications and seasonal wet and dry periods when compared to independent water maps derived from Landsat-8 imagery. The new fwLBand algorithms and continuing SMAP and AMSR2 operations provide for near real-time, multi-scale monitoring of global surface water inundation dynamics and potential flood risk.

Cold season emissions dominate the Arctic tundra methane budget.

Arctic terrestrial ecosystems are major global sources of methane (CH4); hence, it is important to understand the seasonal and climatic controls on CH4 emissions from these systems. Here, we report year-round CH4 emissions from Alaskan Arctic tundra eddy flux sites and regional fluxes derived from aircraft data. We find that emissions during the cold season (September to May) account for ≥ 50% of the annual CH4 flux, with the highest emissions from noninundated upland tundra. A major fraction of cold season emissions occur during the "zero curtain" period, when subsurface soil temperatures are poised near 0 °C. The zero curtain may persist longer than the growing season, and CH4 emissions are enhanced when the duration is extended by a deep thawed layer as can occur with thick snow cover. Regional scale fluxes of CH4 derived from aircraft data demonstrate the large spatial extent of late season CH4 emissions. Scaled to the circumpolar Arctic, cold season fluxes from tundra total 12 ± 5 (95% confidence interval) Tg CH4 y(-1), ∼ 25% of global emissions from extratropical wetlands, or ∼ 6% of total global wetland methane emissions. The dominance of late-season emissions, sensitivity to soil environmental conditions, and importance of dry tundra are not currently simulated in most global climate models. Because Arctic warming disproportionally impacts the cold season, our results suggest that higher cold-season CH4 emissions will result from observed and predicted increases in snow thickness, active layer depth, and soil temperature, representing important positive feedbacks on climate warming.