Does growing atmospheric CO2 explain increasing carbon sink in a boreal coniferous forest?

Impact of vegetation composition and seasonality on sensitivity of modelled CO2 exchange in temperate raised bogs

Encroachment of vascular plants (VP) in temperate raised bogs, as a consequence of altered hydrological conditions and nutrient input, is widely observed. Effects of such vegetation shift on water and carbon cycles are, however, largely unknown and identification of responsible plant physiological traits is challenging. Process-based modelling offers the opportunity of gaining insights into ecosystem functioning beyond observations, and to infer decisive trait shifts of plant functional groups. We adapted the Soil-Vegetation-Atmosphere Transfer model pyAPES to a temperate raised bog site by calibration against measured peat temperature, water table and surface CO2 fluxes. We identified the most important traits determining CO2 fluxes by conducting Morris sensitivity analysis (MSA) under changing conditions throughout the year and simulated VP encroachment. We further investigated transferability of results to other sites by extending MSA to parameter ranges derived from literature review. We found highly variable intra-annual plant traits importance determining ecosystem CO2 fluxes, but only a partial shift of importance of photosynthetic processes from moss to VP during encroachment. Ecosystem respiration was dominated by peat respiration. Overall, carboxylation rate, base respiration rate and temperature sensitivity (Q10) were most important for determining bog CO2 balance and parameter ranking was robust even under the extended MSA.

Vegetation canopy structure mediates the response of gross primary production to environmental drivers across multiple temporal scales

Gross primary production (GPP) is a critical component of the global carbon cycle and plays a significant role in the terrestrial carbon budget. The impact of environmental factors on GPP can occur through both direct (by influencing photosynthetic efficiency) and indirect (through the modulation of vegetation structure) pathways, but the extent to which these mechanisms contribute has been seldom quantified. In this study, we used structural equation modeling and observations from the FLUXNET network to investigate the direct and indirect effects of environmental factors on terrestrial ecosystem GPP at multiple temporal scales. We found that canopy structure, represented by leaf area index (LAI), is a crucial intermediate factor in the GPP response to environmental drivers. Environmental factors affect GPP indirectly by altering canopy structure, and the relative proportion of indirect effects decreased with increasing LAI. The study also identified different effects of environmental factors on GPP across time scales. At the half-hourly time scale, radiation was the primary driver of GPP. In contrast, the influences of temperature and vapor pressure deficit took on greater prominence at longer time scales. About half of the total effect of temperature on GPP was indirect through the regulation of canopy structure, and the indirect effect increased with increasing time scale (GPPNT-based models: 0.135 (half-hourly) vs. 0.171 (daily) vs. 0.189 (weekly) vs. 0.217 (monthly); GPPDT-based models: 0.139 vs. 0.170 vs. 0.187 vs. 0.215; all values were reported in gC m-2 d-1 °C-1, P < 0.001); while the indirect effect of radiation on GPP was comparatively lower, accounting for less than a quarter of the total effect. Furthermore, we observed a direct, negative-to-positive impact of precipitation on GPP across timescales. These findings provide crucial information on the interplay between environmental factors and LAI on GPP and enable a deeper understanding of the driving mechanisms of GPP.

Modeling the response of Norway spruce tree-ring carbon and oxygen isotopes to selection harvest on a drained peatland forest

Continuous cover forestry (CCF) has gained interest as an alternative to even-aged management particularly on drained peatland forests. However, relatively little is known about the physiological response of suppressed trees when larger trees are removed as a part of CCF practices. Consequently, studies concentrating on process-level modeling of the response of trees to selection harvesting are also rare. Here, we compared, modeled and measured harvest response of previously suppressed Norway spruce (Picea abies) trees to a selection harvest. We quantified the harvest response by collecting Norway spruce tree-ring samples in a drained peatland forest site and measuring the change in stable carbon and oxygen isotopic ratios of wood formed during 2010-20, including five post-harvest years. The measured isotopic ratios were compared with ecosystem-level process model predictions for ${\kern0em }^{13}$C discrimination and ${\kern0em }^{18}$O leaf water enrichment. We found that the model predicted similar but lower harvest response than the measurements. Furthermore, accounting for mesophyll conductance was important for capturing the variation in ${\kern0em }^{13}$C discrimination. In addition, we performed sensitivity analysis on the model, which suggests that the modeled ${\kern0em }^{13}$C discrimination is sensitive to parameters related to CO2 transport through stomata to the mesophyll.

Model of methane transport in tree stems: Case study of sap flow and radial diffusion

The transport processes of methane (CH4 ) in tree stems remain largely unknown, although they are critical in assessing the whole-forest CH4 dynamics. We used a physically based dynamic model to study the spatial and diurnal dynamics of stem CH4 transport and fluxes. We parameterised the model using data from laboratory experiments with Pinus sylvestris and Betula pendula and compared the model to experimental data from a field study. Stem CH4 flux in laboratory and field conditions were explained by the axial advective CH4 transport from soil with xylem sap flow and the radial CH4 diffusion through the stem conditions. Diffusion resistance caused by the bark permeability did not significantly affect gas transport or stem CH4 flux in the laboratory experiments. The role of axial diffusion of CH4 in trees was unresolved and requires further studies. Due to the transit time of CH4 in the stem, the diurnal dynamics of stem CH4 fluxes can deviate markedly from the diurnal dynamics of sap flow.

Persistent net release of carbon dioxide and methane from an Alaskan lowland boreal peatland complex

Permafrost degradation in peatlands is altering vegetation and soil properties and impacting net carbon storage. We studied four adjacent sites in Alaska with varied permafrost regimes, including a black spruce forest on a peat plateau with permafrost, two collapse scar bogs of different ages formed following thermokarst, and a rich fen without permafrost. Measurements included year-round eddy covariance estimates of net carbon dioxide (CO2 ), mid-April to October methane (CH4 ) emissions, and environmental variables. From 2011 to 2022, annual rainfall was above the historical average, snow water equivalent increased, and snow-season duration shortened due to later snow return. Seasonally thawed active layer depths also increased. During this period, all ecosystems acted as slight annual sources of CO2 (13-59 g C m-2  year-1 ) and stronger sources of CH4 (11-14 g CH4  m-2 from ~April to October). The interannual variability of net ecosystem exchange was high, approximately ±100 g C m-2  year-1 , or twice what has been previously reported across other boreal sites. Net CO2 release was positively related to increased summer rainfall and winter snow water equivalent and later snow return. Controls over CH4 emissions were related to increased soil moisture and inundation status. The dominant emitter of carbon was the rich fen, which, in addition to being a source of CO2 , was also the largest CH4 emitter. These results suggest that the future carbon-source strength of boreal lowlands in Interior Alaska may be determined by the area occupied by minerotrophic fens, which are expected to become more abundant as permafrost thaw increases hydrologic connectivity. Since our measurements occur within close proximity of each other (≤1 km2 ), this study also has implications for the spatial scale and data used in benchmarking carbon cycle models and emphasizes the necessity of long-term measurements to identify carbon cycle process changes in a warming climate.

Estimating intraseasonal intrinsic water-use efficiency from high-resolution tree-ring δ13 C data in boreal Scots pine forests

Intrinsic water-use efficiency (iWUE), a key index for carbon and water balance, has been widely estimated from tree-ring δ¹³ C at annual resolution, but rarely at high-resolution intraseasonal scale. We estimated high-resolution iWUE from laser-ablation δ¹³ C analysis of tree-rings (iWUE_iso ) and compared it with iWUE derived from gas exchange (iWUE_gas ) and eddy covariance (iWUE_EC ) data for two Pinus sylvestris forests from 2002 to 2019. By carefully timing iWUE_iso via modeled tree-ring growth, iWUE_iso aligned well with iWUE_gas and iWUE_EC at intraseasonal scale. However, year-to-year patterns of iWUE_gas , iWUE_iso , and iWUE_EC were different, possibly due to distinct environmental drivers on iWUE across leaf, tree, and ecosystem scales. We quantified the modification of iWUE_iso by postphotosynthetic δ¹³ C enrichment from leaf sucrose to tree rings and by nonexplicit inclusion of mesophyll and photorespiration terms in photosynthetic discrimination model, which resulted in overestimation of iWUE_iso by up to 11% and 14%, respectively. We thus extended the application of tree-ring δ¹³ C for iWUE estimates to high-resolution intraseasonal scale. The comparison of iWUE_gas , iWUE_iso , and iWUE_EC provides important insights into physiological acclimation of trees across leaf, tree, and ecosystem scales under climate change and improves the upscaling of ecological models.

Explicitly accounting for needle sugar pool size crucial for predicting intra-seasonal dynamics of needle carbohydrates δ18 O and δ13 C

We explore needle sugar isotopic compositions (δ¹⁸ O and δ¹³ C) in boreal Scots pine (Pinus sylvestris) over two growing seasons. A leaf-level dynamic model driven by environmental conditions and based on current understanding of isotope fractionation processes was built to predict δ¹⁸ O and δ¹³ C of two hierarchical needle carbohydrate pools, accounting for the needle sugar pool size and the presence of an invariant pinitol pool. Model results agreed well with observed needle water δ¹⁸ O, δ¹⁸ O and δ¹³ C of needle water-soluble carbohydrates (sugars + pinitol), and needle sugar δ¹³ C (R²  = 0.95, 0.84, 0.60, 0.73, respectively). Relative humidity (RH) and intercellular to ambient CO₂ concentration ratio (C_i /C_a ) were the dominant drivers of δ¹⁸ O and δ¹³ C variability, respectively. However, the variability of needle sugar δ¹⁸ O and δ¹³ C was reduced on diel and intra-seasonal timescales, compared to predictions based on instantaneous RH and C_i /C_a , due to the large needle sugar pool, which caused the signal formation period to vary seasonally from 2 d to more than 5 d. Furthermore, accounting for a temperature-sensitive biochemical ¹⁸ O-fractionation factor and mesophyll resistance in ¹³ C-discrimination were critical. Interpreting leaf-level isotopic signals requires understanding on time integration caused by mixing in the needle sugar pool.

Drivers of Decadal Carbon Fluxes Across Temperate Ecosystems

Long-running eddy covariance flux towers provide insights into how the terrestrial carbon cycle operates over multiple timescales. Here, we evaluated variation in net ecosystem exchange (NEE) of carbon dioxide (CO₂) across the Chequamegon Ecosystem-Atmosphere Study AmeriFlux core site cluster in the upper Great Lakes region of the USA from 1997 to 2020. The tower network included two mature hardwood forests with differing management regimes (US-WCr and US-Syv), two fen wetlands with varying levels of canopy sheltering and vegetation (US-Los and US-ALQ), and a very tall (400 m) landscape-level tower (US-PFa). Together, they provided over 70 site-years of observations. The 19-tower Chequamegon Heterogenous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 campaign centered around US-PFa provided additional information on the spatial variation of NEE. Decadal variability was present in all long-term sites, but cross-site coherence in interannual NEE in the earlier part of the record became weaker with time as non-climatic factors such as local disturbances likely dominated flux time series. Average decadal NEE at the tall tower transitioned from carbon source to sink to near neutral over 24 years. Respiration had a greater effect than photosynthesis on driving variations in NEE at all sites. Declining snowfall offset potential increases in assimilation from warmer springs, as less-insulated soils delayed start of spring green-up. Higher CO₂ increased maximum net assimilation parameters but not total gross primary productivity. Stand-scale sites were larger net sinks than the landscape tower. Clustered, long-term carbon flux observations provide value for understanding the diverse links between carbon and climate and the challenges of upscaling these responses across space.