Abstract The greenhouse gas (GHG) balance of boreal peatlands in permafrost regions will be affected by climate change through disturbances such as permafrost thaw and wildfire. Although the future GHG balance of boreal peatlands including ponds is dominated by the exchange of both carbon dioxide (CO 2 ) and methane (CH 4 ), disturbance impacts on fluxes of the potent GHG nitrous oxide (N 2 O) could contribute to shifts in the net radiative balance. Here, we measured monthly (April to October) fluxes of N 2 O, CH 4 , and CO 2 from three sites located across the sporadic and discontinuous permafrost zones of western Canada. Undisturbed permafrost peat plateaus acted as N 2 O sinks (−0.025 mg N 2 O m −2 d −1 ), but N 2 O uptake was lower from burned plateaus (−0.003 mg N 2 O m −2 d −1 ) and higher following permafrost thaw in the thermokarst bogs (−0.054 mg N 2 O m −2 d −1 ). The thermokarst bogs had below‐ambient N 2 O soil gas concentrations, suggesting that denitrification consumed atmospheric N 2 O during reduction to dinitrogen. Atmospheric uptake of N 2 O in peat plateaus and thermokarst bogs increased with soil temperature and soil moisture, suggesting sensitivity of N 2 O consumption to further climate change. Four of five peatland ponds acted as N 2 O sinks (−0.018 mg N 2 O m −2 d −1 ), with no influence of thermokarst expansion. One pond with high nitrate concentrations had high N 2 O emissions (0.30 mg N 2 O m −2 d −1 ). Overall, our study suggests that the future net radiative balance of boreal peatlands will be dominated by impacts of wildfire and permafrost thaw on CH 4 and CO 2 fluxes, while the influence from N 2 O is minor.

Abstract Climate change is rapidly altering composition, structure, and functioning of the boreal biome, across North America often broadly categorized into ecoregions. The resulting complex changes in different ecoregions present a challenge for efforts to accurately simulate carbon dioxide (CO 2 ) and energy exchanges between boreal forests and the atmosphere with terrestrial ecosystem models (TEMs). Eddy covariance measurements provide valuable information for evaluating the performance of TEMs and guiding their development. Here, we compiled a boreal forest model benchmarking dataset for North America by harmonizing eddy covariance and supporting measurements from eight black spruce ( Picea mariana )-dominated, mature forest stands. The eight forest stands, located in six boreal ecoregions of North America, differ in stand characteristics, disturbance history, climate, permafrost conditions and soil properties. By compiling various data streams, the benchmarking dataset comprises data to parameterize, force, and evaluate TEMs. Specifically, it includes half-hourly, gap-filled meteorological forcing data, ancillary data essential for model parameterization, and half-hourly, gap-filled or partitioned component flux data on CO 2 (net ecosystem production, gross primary production [GPP], and ecosystem respiration [ER]) and energy (latent [LE] and sensible heat [H]) and their daily aggregates screened based on half-hourly gap-filling quality criteria. We present a case study with the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC) to: (1) demonstrate the utility of our dataset to benchmark TEMs and (2) provide guidance for model development and refinement. Model skill was evaluated using several statistical metrics and further examined through the flux responses to their environmental controls. Our results suggest that CLASSIC tended to overestimate GPP and ER among all stands. Model performance regarding the energy fluxes (i.e., LE and H) varied greatly among the stands and exhibited a moderate correlation with latitude. We identified strong relationships between simulated fluxes and their environmental controls except for H, thus highlighting current strengths and limitations of CLASSIC.

Abstract Wetlands are responsible for 20%–31% of global methane (CH 4 ) emissions and account for a large source of uncertainty in the global CH 4 budget. Data‐driven upscaling of CH 4 fluxes from eddy covariance measurements can provide new and independent bottom‐up estimates of wetland CH 4 emissions. Here, we develop a six‐predictor random forest upscaling model (UpCH4), trained on 119 site‐years of eddy covariance CH 4 flux data from 43 freshwater wetland sites in the FLUXNET‐CH4 Community Product. Network patterns in site‐level annual means and mean seasonal cycles of CH 4 fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash‐Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH 4 emissions of 146 ± 43 TgCH 4 y −1 for 2001–2018 which agrees closely with current bottom‐up land surface models (102–181 TgCH 4 y −1 ) and overlaps with top‐down atmospheric inversion models (155–200 TgCH 4 y −1 ). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH 4 fluxes has the potential to produce realistic extra‐tropical wetland CH 4 emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid‐to‐arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC ( https://doi.org/10.3334/ORNLDAAC/2253 ).

Abstract Climate change and permafrost thaw may impact the mobilization of terrestrial dissolved organic carbon (DOC), mercury (Hg), and neurotoxic methylmercury (MeHg) into aquatic ecosystems; thus, understanding processes that control analyte export in northern catchments is needed. We monitored water chemistry for 3 years (2019–2021) at a peatland catchment (Scotty Creek) and a mixed catchment (Smith Creek) in the Dehcho (Northwest Territories), within the discontinuous permafrost zone of boreal western Canada. The peatland catchment had higher DOC and dissolved MeHg, but lower total Hg concentrations (mean ± standard deviation; 19 ± 2.6 mg DOC L −1 ; 0.08 ± 0.04 ng DMeHg L −1 ; 1.1 ± 0.3 ng THg L −1 ) than the mixed catchment (12 ± 4.4 mg DOC L −1 ; 0.05 ± 0.01 ng DMeHg L −1 ; 3.1 ± 2.2 ng THg L −1 ). Analyte concentrations increased with discharge at the mixed catchment, suggesting transport limitation and the flushing of near‐surface, organic‐rich flow paths during wet periods. In contrast, analyte concentrations in the peatland catchment were not primarily associated with discharge. MeHg concentrations, MeHg:THg, and MeHg:DOC increased with water temperature, suggesting enhanced Hg methylation during warmer periods. Mean open water season DOC and total MeHg yields were greater and more variable from the peatland than the mixed catchment (1.1–6.6 vs. 1.4–2.4 g DOC m −2 ; 5.2–36 vs. 6.1–10 ng MeHg m −2 ). Crucial storage thresholds controlling runoff generation likely drove greater inter‐annual variability in analyte yields from the peatland catchment. Our results suggest climate change may influence the production and transport of MeHg from boreal‐Arctic catchments as temperatures increase, peatlands thaw, and runoff generation is altered.

Abstract. 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, https://doi.org/10.3334/ORNLDAAC/1934).

Peatlands and forests cover large areas of the boreal biome and are critical for global climate regulation. They also regulate regional climate through heat and water vapour exchange with the atmosphere. Understanding how land-atmosphere interactions in peatlands differ from forests may therefore be crucial for modelling boreal climate system dynamics and for assessing climate benefits of peatland conservation and restoration. To assess the biophysical impacts of peatlands and forests on peak growing season air temperature and humidity, we analysed surface energy fluxes and albedo from 35 peatlands and 37 evergreen needleleaf forests - the dominant boreal forest type - and simulated air temperature and vapour pressure deficit (VPD) over hypothetical homogeneous peatland and forest landscapes. We ran an evapotranspiration model using land surface parameters derived from energy flux observations and coupled an analytical solution for the surface energy balance to an atmospheric boundary layer (ABL) model. We found that peatlands, compared to forests, are characterized by higher growing season albedo, lower aerodynamic conductance, and higher surface conductance for an equivalent VPD. This combination of peatland surface properties results in a ∼20% decrease in afternoon ABL height, a cooling (from 1.7 to 2.5 °C) in afternoon air temperatures, and a decrease in afternoon VPD (from 0.4 to 0.7 kPa) for peatland landscapes compared to forest landscapes. These biophysical climate impacts of peatlands are most pronounced at lower latitudes (∼45°N) and decrease toward the northern limit of the boreal biome (∼70°N). Thus, boreal peatlands have the potential to mitigate the effect of regional climate warming during the growing season. The biophysical climate mitigation potential of peatlands needs to be accounted for when projecting the future climate of the boreal biome, when assessing the climate benefits of conserving pristine boreal peatlands, and when restoring peatlands that have experienced peatland drainage and mining. © 2020 The Author(s). Published by IOP Publishing Ltd. (Less)