Response of Net Ecosystem Productivity of Three Boreal Forest Stands to Drought

Daylength predominates the bud growth initiation of winter deciduous forest trees in the monsoon region of China

Climate warming has induced significant shifts in spring phenology over both temperate and boreal forests. The timing of bud growth resuming from dormancy is crucial for predicting spring phenology. However, the mechanisms by which environmental cues, other than chilling accumulation, initiate bud growth remains unclear. By constructing a revised process-based spring phenology model incorporating photoperiod and temperature triggers of bud growth, we simulated the first leaf unfolding and first flowering dates of four deciduous forest trees during 1981-2014 at 102 stations across China's monsoon regions. Then, we revealed spatial patterns of the two triggers. Moreover, we compared fitting precision and robustness of the revised model with three mainstream models. Results show that the revised models can effectively simulate all spring phenology time series. Growth initiation of foliar and floral buds was induced by photoperiod lengthening in 80.8% and 77.7% of time series, and by temperature increasing in remaining 19.2% and 22.3% of time series, respectively. The proportions of time series with photoperiod- and temperature-initiated bud growth significantly increase and decrease from northern to southern climatic zones, respectively. Chilling exposure controls the predominant bud growth triggers in different climate zones. Specifically, in regions with long and severe winters where chilling requirement is easily fulfilled, rising temperature in spring alleviates the cold constraint and initiate bud growth. Conversely, in regions with short and mild winters, prolonged daylength in spring compensates the lack of chilling exposure to initiate bud growth. These findings suggest that photoperiod may limit spring phenology response to temperature in low-latitudes. Overall, our model slightly outperforms other models in terms of efficiency, accuracy, and robustness in modeling leaf unfolding and flowering dates. Therefore, this study deepens our understanding of the mechanisms of spring phenology, and improves the predicting capability of spring phenology models in the face of ongoing global warming.

Monitoring nature's calendar from space: Emerging topics in land surface phenology and associated opportunities for science applications

Vegetation phenology has been viewed as the nature's calendar and an integrative indicator of plant-climate interactions. The correct representation of vegetation phenology is important for models to accurately simulate the exchange of carbon, water, and energy between the vegetated land surface and the atmosphere. Remote sensing has advanced the monitoring of vegetation phenology by providing spatially and temporally continuous data that together with conventional ground observations offers a unique contribution to our knowledge about the environmental impact on ecosystems as well as the ecological adaptations and feedback to global climate change. Land surface phenology (LSP) is defined as the use of satellites to monitor seasonal dynamics in vegetated land surfaces and to estimate phenological transition dates. LSP, as an interdisciplinary subject among remote sensing, ecology, and biometeorology, has undergone rapid development over the past few decades. Recent advances in sensor technologies, as well as data fusion techniques, have enabled novel phenology retrieval algorithms that refine phenology details at even higher spatiotemporal resolutions, providing new insights into ecosystem dynamics. As such, here we summarize the recent advances in LSP and the associated opportunities for science applications. We focus on the remaining challenges, promising techniques, and emerging topics that together we believe will truly form the very frontier of the global LSP research field.

Spatial heterogeneity of soil carbon exchanges and their drivers in a boreal forest

Boreal forests have a large impact on the global greenhouse gas balance and their soils constitute an important carbon (C) reservoir. Mature boreal forests are typically a net CO₂ sink, but there are also examples of boreal forests that are persistent CO₂ sources. The reasons remain often unknown, presumably due to a lack of understanding of how biotic and abiotic drivers interact to determine the microbial respiration of soil organic matter (SOM). This study aimed at identifying the main drivers of microbial SOM respiration and CO₂ and CH₄ soil chamber-fluxes within dry and wet sampling areas at the mature boreal forest of Norunda, Sweden, a persistent net CO₂ source. The spatial heterogeneity of the drivers was assessed with a geostatistical approach combined with stepwise multiple regression. We found that heterotrophic soil respiration increased with SOM content and nitrogen (N) availability, while the SOM reactivity, i.e., SOM specific respiration, was determined by soil moisture and N availability. The latter suggests that microbial activity was N rather than C limited and that microbial N mining might be driving old-SOM decomposition, which was observed through a positive correlation between soil respiration and its δ¹³C values. SOM specific heterotrophic respiration was lower in wet than in dry areas, while no such dependencies were found for chamber-based soil CO₂ fluxes, implying that oxygen depletion resulted in lower SOM reactivity. The chamber-based soil CH₄ flux differed significantly between the wet and dry areas. In the wet area, we observed net CH₄ emission that was positively related to soil moisture and NH₄⁺-N content. Taken together, our findings suggest that N availability has a strong regulatory effect on soil CO₂ and CH₄ emissions at Norunda, and that microbial decomposition of old-SOM to release bioavailable N might be partly responsible for the net CO₂ emission at the site.

Current and future impacts of drought and ozone stress on Northern Hemisphere forests

Rising ozone (O₃ ) concentrations, coupled with an increase in drought frequency due to climate change, pose a threat to plant growth and productivity which could negatively affect carbon sequestration capacity of Northern Hemisphere (NH) forests. Using long-term observations of O₃ mixing ratios and soil water content (SWC), we implemented empirical drought and O₃ stress parameterizations in a coupled stomatal conductance-photosynthesis model to assess their impacts on plant gas exchange at three FLUXNET sites: Castelporziano, Blodgett and Hyytiälä. Model performance was evaluated by comparing model estimates of gross primary productivity (GPP) and latent heat fluxes (LE) against present-day observations. CMIP5 GCM model output data were then used to investigate the potential impact of the two stressors on forests by the middle (2041-2050) and end (2091-2100) of the 21st century. We found drought stress was the more significant as it reduced model overestimation of GPP and LE by ~11%-25% compared to 1%-11% from O₃ stress. However, the best model fit to observations at all the study sites was obtained with O₃ and drought stress combined, such that the two stressors counteract the impact of each other. With the inclusion of drought and O₃ stress, GPP at CPZ, BLO and HYY is projected to increase by 7%, 5% and 8%, respectively, by mid-century and by 14%, 11% and 14% by 2091-2100 as atmospheric CO₂ increases. Estimates were up to 21% and 4% higher when drought and O₃ stress were neglected respectively. Drought stress will have a substantial impact on plant gas exchange and productivity, off-setting and possibly negating CO₂ fertilization gains in future, suggesting projected increases in the frequency and severity of droughts in the NH will play a significant role in forest productivity and carbon budgets in future.

Canopy photosynthetic capacity drives contrasting age dynamics of resource use efficiencies between mature temperate evergreen and deciduous forests

Forest resource use efficiencies (RUEs) can vary with tree age, but the nature of these trends and their underlying mechanisms are not well understood. Understanding the age dynamics of forest RUEs and their drivers is vital for assessing the trade-offs between forest functions and resource consumption, making rational management policy, and projecting ecosystem carbon dynamics. Here we used the FLUXNET2015 and AmeriFlux datasets and published literature to explore the age-dependent variability of forest light use efficiency (LUE) and inherent water use efficiency as well as their main regulatory variables in temperate regions. Our results showed that evergreen forest RUEs initially increased before reaching the mature stage (i.e., around 90 years old), and then gradually declined; in contrast, RUEs continuously increased with age for mature deciduous forests. Changing canopy photosynthetic capacity (A_max) was the primary cause of age-related changes in RUEs across temperate forest sites. More importantly, soil nitrogen (N) increased in mature deciduous forests through time but decreased in older evergreen forests. The age-dependent changes in soil N were closely linked with the age dynamics of A_max for mature temperate forests. Additionally, soil nutrient conditions played a greater role in deciduous forest RUEs than evergreen forest RUEs. This study highlights the importance of stand age and forest type on temperate forest RUEs over the long term.

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.

Divergent long-term trends and interannual variation in ecosystem resource use efficiencies of a southern boreal old black spruce forest 1999-2017

Long-term trends in ecosystem resource use efficiencies (RUEs) and their controlling factors are key pieces of information for understanding how an ecosystem responds to climate change. We used continuous eddy covariance and microclimate data over the period 1999-2017 from a 120-year-old black spruce stand in central Saskatchewan, Canada, to assess interannual variability, long-term trends, and key controlling factors of gross ecosystem production (GEP) and the RUEs of carbon (CUE = net primary production [NPP]/GEP), light (LUE = GEP/absorbed photosynthetic radiation [APAR]), and water (WUE = GEP/evapotranspiration [E]). At this site, annual GEP has shown an increasing trend over the 19 years (p < 0.01), which may be attributed to rising atmospheric CO₂ concentration. Interannual variability in GEP, aside from its increasing trend, was most strongly related to spring temperatures. Associated with the significant increase in annual GEP were relatively small changes in NPP, APAR, and E, so that annual CUE showed a decreasing trend and annual LUE and WUE showed increasing trends over the 19 years. The long-term trends in the RUEs were related to the increasing CO₂ concentration. Further analysis of detrended RUEs showed that their interannual variation was impacted most strongly by air temperature. Two-factor linear models combining CO₂ concentration and air temperature performed well (R² ~0.60) in simulating annual RUEs. LUE and WUE were positively correlated both annually and seasonally, while LUE and CUE were mostly negatively correlated. Our results showed divergent long-term trends among CUE, LUE, and WUE and highlighted the need to account for the combined effects of climatic controls and the 'CO₂ fertilization effect' on long-term variations in RUEs. Since most RUE-based models rely primarily on one resource limitation, the observed patterns of relative change among the three RUEs may have important implications for RUE-based modeling of C fluxes.

Spatiotemporal Simulation of Net Ecosystem Productivity and Its Response to Climate Change in Subtropical Forests

Subtropical forests have great potential as carbon sinks; however, the relationship between net ecosystem productivity (NEP) and climate change is still unclear. This study took Zhejiang Province, a subtropical region, as an example. Based on remote sensing classification data of forest resources, the integrated terrestrial ecosystem carbon cycle (InTEC) model was used to simulate the spatiotemporal dynamics of the forest NEP in Zhejiang Province during 1985–2015 and analyze its response to meteorological factors such as temperature, precipitation, relative humidity, and radiation. Three patterns emerged: (1) The optimized InTEC model can better simulate the forest NEP in Zhejiang Province, and the correlation coefficient between the simulated NEP and observed NEP was up to 0.75. (2) From 1985 to 2015, the increase in the total NEP was rapid, with an average annual growth rate of 1.52 Tg·C·yr−1. During 1985–1988, the forests in Zhejiang Province were carbon sources. After 1988, the forests turned into carbon sinks and this continued to increase. During 2000–2015, more than 97% of the forests in Zhejiang Province were carbon sinks. The total NEP reached 32.02 Tg·C·yr−1, and the annual mean NEP increased to 441.91 gC·m−2·yr−1. The carbon sequestration capacity of forests in the east and southwest of Zhejiang Province is higher than that in the northeast of Zhejiang Province. (3) From 2000 to 2015, there was an extremely significant correlation between forest NEP and precipitation, with a correlation coefficient of 0.85. Simultaneously, the forest NEP showed a negative correlation with temperature and radiation, with a correlation coefficient of −0.56 for both, and the forest NEP was slightly negatively correlated with relative humidity. The relative contribution rates of temperature, precipitation, relative humidity, and radiation data to NEP showed that the contribution of precipitation to NEP is the largest, reaching 61%, followed by temperature and radiation at 18% and 17%, respectively. The relative contribution rate of relative humidity is the smallest at only 4%. During the period of 1985–1999, due to significant man-made disturbances, the NEP had a weak correlation with temperature, precipitation, relative humidity, and radiation. The results of this study are important for addressing climate change and illustrating the response mechanism between subtropical forest NEP and climate change.

Changes of vegetation carbon sequestration in the tableland of Loess Plateau and its influencing factors

The variations of vegetation carbon sequestration have become a gauge for evaluating the ecological effect of vegetation restoration. In this study, the spatiotemporal patterns of the net ecosystem production (NEP) were simulated using an improved CASA model and GSMSR model. It showed that the NEP markedly increased in the tableland of Loess Plateau during 2003-2012, with an annual average growth of 3.65 g C·m^-2 a^-1. The mixed broadleaf-conifer forest ranked first (127.23 g C·m^-2 a^-1) while the bare land and sparse vegetation presented the lowest carbon sequestration (14.64 g C·m^-2 a^-1). The NEP manifested a significantly uneven overall spatial distribution: high in the southwest and low in the northeast. The spatial variations of NEP resulted from the combined effects of geographic position, terrain, meteorology, and soil and vegetation, respectively. Quantitative isolation revealed that the most dominant factor of vegetation carbon sequestration was soil and vegetation, while terrain exerted insignificant impacts on the NEP.

The impact of rising temperatures on water balance and phenology of European beech (Fagus sylvatica L.) stands

In this article, we outline the set-up and the application of an eco-hydrological box model, with the aim to describe the water balance of deciduous (Fagus Sylvatica L.) forest stands. The water balance model (WBM) uses standard meteorological parameters as input variables and runs on a daily time step. It consists of two modules. The aboveground module (1) comprises routines for fog precipitation generation, precipitation interception and snowfall/snowmelt dynamics. Covered belowground processes (2) are bypass flow, percolation, soil evaporation and transpiration, where the latter two processes are considered separately. Preceding to the WBM, a routine is introduced, specifying the intra-annual foliage dynamics of beech. Emphasis is also laid on the inter-annual variation of beech phenology. Leaf sprouting and leaf senescence are calculated as functions of day-length and air temperature. The WBM was applied to four European beech dominated forest stands in the northeastern part of Austria. They are located on a gradient of declining annual precipitation (from west to east). The two easterly sites are located close to the (dry) limit of the natural distribution of beech. Records of soil moisture were used for the adjustment of 26 parameters. On all sites the calibration process (simulated annealing) delivered good predictions of soil moisture (Nash–Sutcliffe efficiency≥ 0.925). Then, the obtained parameterization was used to apply different scenarios of global warming. The temperature was increased step-wisely up to 4 °C. All scenarios were run (1) with present phenological conditions and (2) with phenology responding to higher temperatures. This way, we wanted to assign the effect of higher temperatures and longer growing seasons on the water dynamics of the forest stands. A warming of 1 °C corresponded roughly to an elongation of the growing season of 4.5 days, where the start of the growing season was affected more strongly than the end. Apparently, higher temperatures led to drier soils. The strongest change was observed in early summer, also amplified by an earlier start of the growing season. Rising temperatures led to lower export fluxes of liquid water, simultaneously increasing evapotranspiration (ET). The gain in ET was almost entirely assignable to increased soil evaporation. Drier soils led to a sharp depression of transpiration during summer months. This decline was compensated by the effect of elongated growing seasons. The risk of severe drought was increased by higher temperatures, but here the contribution of growing season length was negligible. Drier soils seem to hamper the stands’ productivity. For all warming scenarios, the estimated increase of the gross primary production, caused by longer periods of assimilation, is nullified by the effect of soil water deficit in mid-summer.

Direct and indirect climate change effects on carbon dioxide fluxes in a thawing boreal forest-wetland landscape

In the sporadic permafrost zone of northwestern Canada, boreal forest carbon dioxide (CO₂ ) fluxes will be altered directly by climate change through changing meteorological forcing and indirectly through changes in landscape functioning associated with thaw-induced collapse-scar bog ('wetland') expansion. However, their combined effect on landscape-scale net ecosystem CO₂ exchange (NEE_LAND ), resulting from changing gross primary productivity (GPP) and ecosystem respiration (ER), remains unknown. Here, we quantify indirect land cover change impacts on NEE_LAND and direct climate change impacts on modeled temperature- and light-limited NEE_LAND of a boreal forest-wetland landscape. Using nested eddy covariance flux towers, we find both GPP and ER to be larger at the landscape compared to the wetland level. However, annual NEE_LAND (-20 g C m^-2 ) and wetland NEE (-24 g C m^-2 ) were similar, suggesting negligible wetland expansion effects on NEE_LAND . In contrast, we find non-negligible direct climate change impacts when modeling NEE_LAND using projected air temperature and incoming shortwave radiation. At the end of the 21st century, modeled GPP mainly increases in spring and fall due to reduced temperature limitation, but becomes more frequently light-limited in fall. In a warmer climate, ER increases year-round in the absence of moisture stress resulting in net CO₂ uptake increases in the shoulder seasons and decreases during the summer. Annually, landscape net CO₂ uptake is projected to decline by 25 ± 14 g C m^-2 for a moderate and 103 ± 38 g C m^-2 for a high warming scenario, potentially reversing recently observed positive net CO₂ uptake trends across the boreal biome. Thus, even without moisture stress, net CO₂ uptake of boreal forest-wetland landscapes may decline, and ultimately, these landscapes may turn into net CO₂ sources under continued anthropogenic CO₂ emissions. We conclude that NEE_LAND changes are more likely to be driven by direct climate change rather than by indirect land cover change impacts.

Temperature and snowfall trigger alpine vegetation green-up on the world's roof

Rapid temperature increase and its impacts on alpine ecosystems in the Qinghai-Tibetan Plateau, the world's highest and largest plateau, are a matter of global concern. Satellite observations have revealed distinctly different trend changes and contradicting temperature responses of vegetation green-up dates, leading to broad debate about the Plateau's spring phenology and its climatic attribution. Large uncertainties in remote-sensing estimates of phenology significantly limit efforts to predict the impacts of climate change on vegetation growth and carbon balance in the Qinghai-Tibetan Plateau, which are further exacerbated by a lack of detailed ground observation calibration. Here, we revealed the spatiotemporal variations and climate drivers of ground-based herbaceous plant green-up dates using 72 green-up datasets for 22 herbaceous plant species at 23 phenological stations, and corresponding daily mean air temperature and daily precipitation data from 19 climate stations across eastern and southern parts of the Qinghai-Tibetan Plateau from 1981 to 2011. Results show that neither the continuously advancing trend from 1982 to 2011, nor a turning point in the mid to late 1990s as reported by remote-sensing studies can be verified by most of the green-up time series, and no robust evidence for a warmer winter-induced later green-up dates can be detected. Thus, chilling requirements may not be an important driver influencing green-up responses to spring warming. Moreover, temperature-only control of green-up dates appears mainly at stations with relatively scarce preseason snowfall and lower elevation, while coupled temperature and precipitation controls of green-up dates occur mostly at stations with relatively abundant preseason snowfall and higher elevation. The diversified interactions between snowfall and temperature during late winter to early spring likely determine the spatiotemporal variations of green-up dates. Therefore, prediction of vegetation growth and carbon balance responses to global climate change on the world's roof should integrate both temperature and snowfall variations.

Climate-driven uncertainties in modeling terrestrial gross primary production: a site level to global-scale analysis

We used a land surface model to quantify the causes and extents of biases in terrestrial gross primary production (GPP) due to the use of meteorological reanalysis datasets. We first calibrated the model using meteorology and eddy covariance data from 25 flux tower sites ranging from the tropics to the northern high latitudes and subsequently repeated the site simulations using two reanalysis datasets: NCEP/NCAR and CRUNCEP. The results show that at most sites, the reanalysis-driven GPP bias was significantly positive with respect to the observed meteorology-driven simulations. Notably, the absolute GPP bias was highest at the tropical evergreen tree sites, averaging up to ca. 0.45 kg C m(-2)  yr(-1) across sites (ca. 15% of site level GPP). At the northern mid-/high-latitude broadleaf deciduous and the needleleaf evergreen tree sites, the corresponding annual GPP biases were up to 20%. For the nontree sites, average annual biases of up to ca. 20-30% were simulated within savanna, grassland, and shrubland vegetation types. At the tree sites, the biases in short-wave radiation and humidity strongly influenced the GPP biases, while the nontree sites were more affected by biases in factors controlling water stress (precipitation, humidity, and air temperature). In this study, we also discuss the influence of seasonal patterns of meteorological biases on GPP. Finally, using model simulations for the global land surface, we discuss the potential impacts of site-level reanalysis-driven biases on the global estimates of GPP. In a broader context, our results can have important consequences on other terrestrial ecosystem fluxes (e.g., net primary production, net ecosystem production, energy/water fluxes) and reservoirs (e.g., soil carbon stocks). In a complementary study (Barman et al., ), we extend the present analysis for latent and sensible heat fluxes, thus consistently integrating the analysis of climate-driven uncertainties in carbon, energy, and water fluxes using a single modeling framework.

Autumn warming reduces the CO2 sink of a black spruce forest in interior Alaska based on a nine-year eddy covariance measurement

Nine years (2003-2011) of carbon dioxide (CO2) flux were measured at a black spruce forest in interior Alaska using the eddy covariance method. Seasonal and interannual variations in the gross primary productivity (GPP) and ecosystem respiration (RE) were associated primarily with air temperature: warmer conditions enhanced GPP and RE. Meanwhile, interannual variation in annual CO2 balance was controlled predominantly by RE, and not GPP. During these 9 years of measurement, the annual CO2 balance shifted from a CO2 sink to a CO2 source, with a 9-year average near zero. The increase in autumn RE was associated with autumn warming and was mostly attributed to a shift in the annual CO2 balance. The increase in autumn air temperature (0.22 °C yr(-1)) during the 9 years of study was 15 times greater than the long-term warming trend between 1905 and 2011 (0.015 °C yr(-1)) due to decadal climate oscillation. This result indicates that most of the shifts in observed CO2 fluxes were associated with decadal climate variability. Because the natural climate varies in a cycle of 10-30 years, a long-term study covering at least one full cycle of decadal climate oscillation is important to quantify the CO2 balance and its interaction with the climate.