Abstract Photosynthesis of terrestrial ecosystems in the Arctic-Boreal region is a critical part of the global carbon cycle. Solar-induced chlorophyll Fluorescence (SIF), a promising proxy for photosynthesis with physiological insight, has been used to track gross primary production (GPP) at regional scales. Recent studies have constructed empirical relationships between SIF and eddy covariance-derived GPP as a first step to predicting global GPP. However, high latitudes pose two specific challenges: (a) Unique plant species and land cover types in the Arctic–Boreal region are not included in the generalized SIF-GPP relationship from lower latitudes, and (b) the complex terrain and sub-pixel land cover further complicate the interpretation of the SIF-GPP relationship. In this study, we focused on the Arctic-Boreal vulnerability experiment (ABoVE) domain and evaluated the empirical relationships between SIF for high latitudes from the TROPOspheric Monitoring Instrument (TROPOMI) and a state-of-the-art machine learning GPP product (FluxCom). For the first time, we report the regression slope, linear correlation coefficient, and the goodness of the fit of SIF-GPP relationships for Arctic-Boreal land cover types with extensive spatial coverage. We found several potential issues specific to the Arctic-Boreal region that should be considered: (a) unrealistically high FluxCom GPP due to the presence of snow and water at the subpixel scale; (b) changing biomass distribution and SIF-GPP relationship along elevational gradients, and (c) limited perspective and misrepresentation of heterogeneous land cover across spatial resolutions. Taken together, our results will help improve the estimation of GPP using SIF in terrestrial biosphere models and cope with model-data uncertainties in the Arctic-Boreal region.

Abstract Photosynthesis of terrestrial ecosystems in the Arctic-Boreal region is a critical part of the global carbon cycle. Solar-induced chlorophyll Fluorescence (SIF), a promising proxy for photosynthesis with physiological insight, has been used to track gross primary production (GPP) at regional scales. Recent studies have constructed empirical relationships between SIF and eddy covariance-derived GPP as a first step to predicting global GPP. However, high latitudes pose two specific challenges: (a) Unique plant species and land cover types in the Arctic–Boreal region are not included in the generalized SIF-GPP relationship from lower latitudes, and (b) the complex terrain and sub-pixel land cover further complicate the interpretation of the SIF-GPP relationship. In this study, we focused on the Arctic-Boreal vulnerability experiment (ABoVE) domain and evaluated the empirical relationships between SIF for high latitudes from the TROPOspheric Monitoring Instrument (TROPOMI) and a state-of-the-art machine learning GPP product (FluxCom). For the first time, we report the regression slope, linear correlation coefficient, and the goodness of the fit of SIF-GPP relationships for Arctic-Boreal land cover types with extensive spatial coverage. We found several potential issues specific to the Arctic-Boreal region that should be considered: (a) unrealistically high FluxCom GPP due to the presence of snow and water at the subpixel scale; (b) changing biomass distribution and SIF-GPP relationship along elevational gradients, and (c) limited perspective and misrepresentation of heterogeneous land cover across spatial resolutions. Taken together, our results will help improve the estimation of GPP using SIF in terrestrial biosphere models and cope with model-data uncertainties in the Arctic-Boreal region.

Forests play a crucial role in the global carbon (C) cycle by storing and sequestering a substantial amount of C in the terrestrial biosphere. Due to temporal dynamics in climate and vegetation activity, there are significant regional variations in carbon dioxide (CO2) fluxes between the biosphere and atmosphere in forests that are affecting the global C cycle. Current forest CO2 flux dynamics are controlled by instantaneous climate, soil, and vegetation conditions, which carry legacy effects from disturbances and extreme climate events. Our level of understanding from the legacies of these processes on net CO2 fluxes is still limited due to their complexities and their long-term effects. Here, we combined remote sensing, climate, and eddy-covariance flux data to study net ecosystem CO2 exchange (NEE) at 185 forest sites globally. Instead of commonly used non-dynamic statistical methods, we employed a type of recurrent neural network (RNN), called Long Short-Term Memory network (LSTM) that captures information from the vegetation and climate's temporal dynamics. The resulting data-driven model integrates interannual and seasonal variations of climate and vegetation by using Landsat and climate data at each site. The presented LSTM algorithm was able to effectively describe the overall seasonal variability (Nash-Sutcliffe efficiency, NSE = 0.66) and across-site (NSE = 0.42) variations in NEE, while it had less success in predicting specific seasonal and interannual anomalies (NSE = 0.07). This analysis demonstrated that an LSTM approach with embedded climate and vegetation memory effects outperformed a non-dynamic statistical model (i.e. Random Forest) for estimating NEE. Additionally, it is shown that the vegetation mean seasonal cycle embeds most of the information content to realistically explain the spatial and seasonal variations in NEE. These findings show the relevance of capturing memory effects from both climate and vegetation in quantifying spatio-temporal variations in forest NEE.

Forests dominate carbon (C) exchanges between the terrestrial biosphere and the atmosphere on land. In the long term, the net carbon flux between forests and the atmosphere has been significantly impacted by changes in forest cover area and structure due to ecological disturbances and management activities. Current empirical approaches for estimating net ecosystem productivity (NEP) rarely consider forest age as a predictor, which represents variation in physiological processes that can respond differently to environmental drivers, and regrowth following disturbance. Here, we conduct an observational synthesis to empirically determine to what extent climate, soil properties, nitrogen deposition, forest age and management influence the spatial and interannual variability of forest NEP across 126 forest eddy-covariance flux sites worldwide. The empirical models explained up to 62% and 71% of spatio-temporal and across-site variability of annual NEP, respectively. An investigation of model structures revealed that forest age was a dominant factor of NEP spatio-temporal variability in both space and time at the global scale as compared to abiotic factors, such as nutrient availability, soil characteristics and climate. These findings emphasize the importance of forest age in quantifying spatio-temporal variation in NEP using empirical approaches.

Abstract. Extreme climatic events, such as droughts and heat stress, induce anomalies in ecosystem–atmosphere CO2 fluxes, such as gross primary production (GPP) and ecosystem respiration (Reco), and, hence, can change the net ecosystem carbon balance. However, despite our increasing understanding of the underlying mechanisms, the magnitudes of the impacts of different types of extremes on GPP and Reco within and between ecosystems remain poorly predicted. Here we aim to identify the major factors controlling the amplitude of extreme-event impacts on GPP, Reco, and the resulting net ecosystem production (NEP). We focus on the impacts of heat and drought and their combination. We identified hydrometeorological extreme events in consistently downscaled water availability and temperature measurements over a 30-year time period. We then used FLUXNET eddy covariance flux measurements to estimate the CO2 flux anomalies during these extreme events across dominant vegetation types and climate zones. Overall, our results indicate that short-term heat extremes increased respiration more strongly than they downregulated GPP, resulting in a moderate reduction in the ecosystem's carbon sink potential. In the absence of heat stress, droughts tended to have smaller and similarly dampening effects on both GPP and Reco and, hence, often resulted in neutral NEP responses. The combination of drought and heat typically led to a strong decrease in GPP, whereas heat and drought impacts on respiration partially offset each other. Taken together, compound heat and drought events led to the strongest C sink reduction compared to any single-factor extreme. A key insight of this paper, however, is that duration matters most: for heat stress during droughts, the magnitude of impacts systematically increased with duration, whereas under heat stress without drought, the response of Reco over time turned from an initial increase to a downregulation after about 2 weeks. This confirms earlier theories that not only the magnitude but also the duration of an extreme event determines its impact. Our study corroborates the results of several local site-level case studies but as a novelty generalizes these findings on the global scale. Specifically, we find that the different response functions of the two antipodal land–atmosphere fluxes GPP and Reco can also result in increasing NEP during certain extreme conditions. Apparently counterintuitive findings of this kind bear great potential for scrutinizing the mechanisms implemented in state-of-the-art terrestrial biosphere models and provide a benchmark for future model development and testing.

Abstract. Extreme climatic events, such as droughts and heat stress induce anomalies in ecosystem-atmosphere CO2 fluxes, such as gross primary production (GPP) and ecosystem respiration (Reco), and, hence, can change the net ecosystem carbon balance. However, despite our increasing understanding of the underlying mechanisms, the magnitudes of the impacts of different types of extremes on GPP and Reco within and between ecosystems remain poorly predicted. Here we aim to identify the major factors controlling the amplitude of extreme event impacts on GPP, Reco, and the resulting net ecosystem production (NEP). We focus on the impacts of heat and drought and their combination. We identified hydrometeorological extreme events in consistently downscaled water availability and temperature measurements over a 30 year time period. We then used FLUXNET eddy-covariance flux measurements to estimate the CO2 flux anomalies during these extreme events across dominant vegetation types and climate zones. Overall, our results indicate that short-term heat extremes increased respiration more strongly than they down-regulated GPP, resulting in a moderate reduction of the ecosystem’s carbon sink potential. In the absence of heat stress, droughts tended to have smaller and similarly dampening effects on both GPP and Reco, and, hence, often resulted in neutral NEP responses. The combination of drought and heat typically led to a strong decrease in GPP, whereas heat and drought impacts on respiration partially offset each other. Taken together, compound heat and drought events led to the strongest C sink reduction compared to any single-factor extreme. A key insight of this paper, however, is that duration matters most: for heat stress during droughts, the magnitude of impacts systematically increased with duration, whereas under heat stress without drought, the response of Reco over time turned from an initial increase to a down-regulation after about two weeks. This confirms earlier theories that not only the magnitude but also the duration of an extreme event determines its impact. Our study corroborates the results of several local site-level case studies, but as a novelty generalizes these findings at the global scale. Specifically, we find that the different response functions of the two antipodal land-atmosphere fluxes GPP and Reco can also result in increasing NEP during certain extreme conditions. Apparently counterintuitive findings of this kind bear great potential for scrutinizing the mechanisms implemented in state-of-the-art terrestrial biosphere models and provide a benchmark for future model development and testing.