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.

Environmental Controls on Multi-Scale Dynamics of Net Carbon Dioxide Exchange From an Alpine Peatland on the Eastern Qinghai-Tibet Plateau.

Peatlands are characterized by their large carbon storage capacity and play an essential role in the global carbon cycle. However, the future of the carbon stored in peatland ecosystems under a changing climate remains unclear. In this study, based on the eddy covariance technique, we investigated the net ecosystem CO2 exchange (NEE) and its controlling factors of the Hongyuan peatland, which is a part of the Ruoergai peatland on the eastern Qinghai-Tibet Plateau (QTP). Our results show that the Hongyuan alpine peatland was a CO2 sink with an annual NEE of -226.61 and -185.35 g C m-2 in 2014 and 2015, respectively. While, the non-growing season NEE was 53.35 and 75.08 g C m-2 in 2014 and 2015, suggesting that non-growing seasons carbon emissions should not be neglected. Clear diurnal variation in NEE was observed during the observation period, with the maximum CO2 uptake appearing at 12:30 (Beijing time, UTC+8). The Q10 value of the non-growing season in 2014 and 2015 was significantly higher than that in the growing season, which suggested that the CO2 flux in the non-growing season was more sensitive to warming than that in the growing season. We investigated the multi-scale temporal variations in NEE during the growing season using wavelet analysis. On daily timescales, photosynthetically active radiation was the primary driver of NEE. Seasonal variation in NEE was mainly driven by soil temperature. The amount of precipitation was more responsible for annual variation of NEE. The increasing number of precipitation event was associated with increasing annual carbon uptake. This study highlights the need for continuous eddy covariance measurements and time series analysis approaches to deepen our understanding of the temporal variability in NEE and multi-scale correlation between NEE and environmental factors.

Effects of drought and meteorological forcing on carbon and water fluxes in Nordic forests during the dry summer of 2018.

The Nordic region was subjected to severe drought in 2018 with a particularly long-lasting and large soil water deficit in Denmark, Southern Sweden and Estonia. Here, we analyse the impact of the drought on carbon and water fluxes in 11 forest ecosystems of different composition: spruce, pine, mixed and deciduous. We assess the impact of drought on fluxes by estimating the difference (anomaly) between year 2018 and a reference year without drought. Unexpectedly, the evaporation was only slightly reduced during 2018 compared to the reference year at two sites while it increased or was nearly unchanged at all other sites. This occurred under a 40 to 60% reduction in mean surface conductance and the concurrent increase in evaporative demand due to the warm and dry weather. The anomaly in the net ecosystem productivity (NEP) was 93% explained by a multilinear regression with the anomaly in heterotrophic respiration and the relative precipitation deficit as independent variables. Most of the variation (77%) was explained by the heterotrophic component. Six out of 11 forests reduced their annual NEP with more than 50 g C m-2 yr-1 during 2018 as compared to the reference year. The NEP anomaly ranged between -389 and +74 g C m-2 yr-1 with a median value of -59 g C m-2 yr-1. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018.

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO2 exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO2 uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Estimating canopy gross primary production by combining phloem stable isotopes with canopy and mesophyll conductances.

Gross primary production (GPP) is a key component of the forest carbon cycle. However, our knowledge of GPP at the stand scale remains uncertain, because estimates derived from eddy covariance (EC) rely on semi-empirical modelling and the assumptions of the EC technique are sometimes not fully met. We propose using the sap flux/isotope method as an alternative way to estimate canopy GPP, termed GPPiso/SF , at the stand scale and at daily resolution. It is based on canopy conductance inferred from sap flux and intrinsic water-use efficiency estimated from the stable carbon isotope composition of phloem contents. The GPPiso/SF estimate was further corrected for seasonal variations in photosynthetic capacity and mesophyll conductance. We compared our estimate of GPPiso/SF to the GPP derived from PRELES, a model parameterized with EC data. The comparisons were performed in a highly instrumented, boreal Scots pine forest in northern Sweden, including a nitrogen fertilized and a reference plot. The resulting annual and daily GPPiso/SF estimates agreed well with PRELES, in the fertilized plot and the reference plot. We discuss the GPPiso/SF method as an alternative which can be widely applied without terrain restrictions, where the assumptions of EC are not met.

The Net Landscape Carbon Balance-Integrating terrestrial and aquatic carbon fluxes in a managed boreal forest landscape in Sweden.

The boreal biome exchanges large amounts of carbon (C) and greenhouse gases (GHGs) with the atmosphere and thus significantly affects the global climate. A managed boreal landscape consists of various sinks and sources of carbon dioxide (CO2 ), methane (CH4 ), and dissolved organic and inorganic carbon (DOC and DIC) across forests, mires, lakes, and streams. Due to the spatial heterogeneity, large uncertainties exist regarding the net landscape carbon balance (NLCB). In this study, we compiled terrestrial and aquatic fluxes of CO2 , CH4 , DOC, DIC, and harvested C obtained from tall-tower eddy covariance measurements, stream monitoring, and remote sensing of biomass stocks for an entire boreal catchment (~68 km2 ) in Sweden to estimate the NLCB across the land-water-atmosphere continuum. Our results showed that this managed boreal forest landscape was a net C sink (NLCB = 39 g C m-2  year-1 ) with the landscape-atmosphere CO2 exchange being the dominant component, followed by the C export via harvest and streams. Accounting for the global warming potential of CH4 , the landscape was a GHG sink of 237 g CO2 -eq m-2  year-1 , thus providing a climate-cooling effect. The CH4 flux contribution to the annual GHG budget increased from 0.6% during spring to 3.2% during winter. The aquatic C loss was most significant during spring contributing 8% to the annual NLCB. We further found that abiotic controls (e.g., air temperature and incoming radiation) regulated the temporal variability of the NLCB whereas land cover types (e.g., mire vs. forest) and management practices (e.g., clear-cutting) determined their spatial variability. Our study advocates the need for integrating terrestrial and aquatic fluxes at the landscape scale based on tall-tower eddy covariance measurements combined with biomass stock and stream monitoring to develop a holistic understanding of the NLCB of managed boreal forest landscapes and to better evaluate their potential for mitigating climate change.

The concentrations and characteristics of dissolved organic matter in high-latitude lakes determine its ambient reducing capacity.

The reducing capacity (RC) of natural organic matter plays an important role in the carbon cycle and biogeochemical fates of environmental contaminants in the aquatic system. However, the electron donation potentials of dissolved organic matter (DOM) from high-latitude lakes are still uncertain. In this study, we collected DOM samples from high-latitude lakes across the Arctic and boreal regions in Sweden and Norway to investigate the effects of the DOM concentration and characteristics on its ambient reducing capacity (ARC). Mercury (Hg(II)) abiotic reduction in darkness was used to determine the ARC. The results showed that the DOM in Arctic lakes is less terrestrial-dominant than in reference sites (i.e., forest lakes). Between the two categories of Arctic lakes, tundra lakes are more terrestrial-influenced compared to mountain lakes. Additionally, terrestrial-originated DOM is a main controlling factor for enhancing the ambient reducing capacity, whereas the DOM concentration, i.e., dissolved organic carbon (DOC), resulted in variations in the Hg/DOC ratios that also cause the variations of the observed ARC values. Thus, comparisons of the ARC values can be conducted while oxidant/DOC ratios are kept the same and reported through the method using heavy metals as a chemical probe. After correction for Hg/DOC ratio interference, the ambient reducing capacity of DOM followed the order: boreal forest lakes > Arctic tundra lakes > Arctic mountain lakes. This study highlights that the DOM concentration should also be considered when estimating the ARC as compared to the previous that mainly focusing on the properties of DOM such as its origins. As climate change is projected to be severe in high latitudes, this study demonstrates a significant connection between aquatic DOM geochemical reactivity and terrestrial inputs, which is crucial for a better prediction of the role of DOM in high-latitude lakes in the context of climate change.