Search

Methane emissions from plant foliage may play an important role in the global methane cycle, but their size and the underlying source processes remain poorly understood. Here, we quantify methane fluxes from the shoots of Scots pine trees, a dominant tree species in boreal forests, to identify source processes and environmental drivers, and we evaluate whether these fluxes can be constrained at the ecosystem-level by eddy covariance flux measurements. We show that shoot-level measurements conducted in forest, garden, or greenhouse settings; on mature trees and saplings; manually and with an automated CO 2 -, temperature-, and water-controlled chamber system; and with multiple methane analyzers all resulted in comparable daytime fluxes (0.144 ± 0.019 to 0.375 ± 0.074 nmol CH 4 g -1 foliar d.w. h -1 ). We further find that these emissions exhibit a pronounced diurnal cycle that closely follows photosynthetically active radiation and is further modulated by temperature. These diurnal patterns indicate that methane production is associated with diurnal cycle of sunlight, indicating that this production is either a byproduct of photosynthesis-associated biochemical reactions (e.g., the methionine cycle) or produced through nonenzymatic photochemical reactions in plant biomass. Moreover, we identified a light-dependent component in stand-level methane fluxes, which showed order-of-magnitude agreement with shoot-level measurements (0.968 ± 0.031 nmol CH 4 g -1 h -1 ) and which provides an upper limit for shoot methane emissions., Competing Interests: Competing interests statement:The authors declare no competing interest.

Wetlands are the largest natural source of methane (CH 4 ) to the atmosphere. The eddy covariance method provides robust measurements of net ecosystem exchange of CH 4 , but interpreting its spatiotemporal variations is challenging due to the co-occurrence of CH 4 production, oxidation, and transport dynamics. Here, we estimate these three processes using a data-model fusion approach across 25 wetlands in temperate, boreal, and Arctic regions. Our data-constrained model-iPEACE-reasonably reproduced CH 4 emissions at 19 of the 25 sites with normalized root mean square error of 0.59, correlation coefficient of 0.82, and normalized standard deviation of 0.87. Among the three processes, CH 4 production appeared to be the most important process, followed by oxidation in explaining inter-site variations in CH 4 emissions. Based on a sensitivity analysis, CH 4 emissions were generally more sensitive to decreased water table than to increased gross primary productivity or soil temperature. For periods with leaf area index (LAI) of ≥20% of its annual peak, plant-mediated transport appeared to be the major pathway for CH 4 transport. Contributions from ebullition and diffusion were relatively high during low LAI (<20%) periods. The lag time between CH 4 production and CH 4 emissions tended to be short in fen sites (3 ± 2 days) and long in bog sites (13 ± 10 days). Based on a principal component analysis, we found that parameters for CH 4 production, plant-mediated transport, and diffusion through water explained 77% of the variance in the parameters across the 19 sites, highlighting the importance of these parameters for predicting wetland CH 4 emissions across biomes. These processes and associated parameters for CH 4 emissions among and within the wetlands provide useful insights for interpreting observed net CH 4 fluxes, estimating sensitivities to biophysical variables, and modeling global CH 4 fluxes., (© 2023 John Wiley & Sons Ltd.)

The terrestrial net ecosystem productivity (NEP) has increased during the past three decades, but the mechanisms responsible are still unclear. We analyzed 17 years (2001-2017) of eddy-covariance measurements of NEP, evapotranspiration (ET) and light and water use efficiency from a boreal coniferous forest in Southern Finland for trends and inter-annual variability (IAV). The forest was a mean annual carbon sink (252 [ ± 42] gC m - 2 a - 1 ), and NEP increased at rate +6.4-7.0 gC m - 2 a - 1 (or ca. +2.5% a - 1 ) during the period. This was attributed to the increasing gross-primary productivity GPP and occurred without detectable change in ET. The start of annual carbon uptake period was advanced by 0.7 d a - 1 , and increase in GPP and NEP outside the main growing season contributed ca. one-third and one-fourth of the annual trend, respectively. Meteorological factors were responsible for the IAV of fluxes but did not explain the long-term trends. The growing season GPP trend was strongest in ample light during the peak growing season. Using a multi-layer ecosystem model, we showed that direct CO 2 fertilization effect diminishes when moving from leaf to ecosystem, and only 30-40% of the observed ecosystem GPP increase could be attributed to CO 2 . The increasing trend in leaf-area index (LAI), stimulated by forest thinning in 2002, was the main driver of the enhanced GPP and NEP of the mid-rotation managed forest. It also compensated for the decrease of mean leaf stomatal conductance with increasing CO 2 and LAI, explaining the apparent proportionality between observed GPP and CO 2 trends. The results emphasize that attributing trends to their physical and physiological drivers is challenged by strong IAV, and uncertainty of LAI and species composition changes due to the dynamic flux footprint. The results enlighten the underlying mechanisms responsible for the increasing terrestrial carbon uptake in the boreal zone., (© 2022 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)

Understanding the critical soil moisture (SM) threshold (θ crit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θ crit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θ crit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θ crit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θ crit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θ crit . We further found that soil texture dominates the spatial variability of θ crit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θ crit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models., (© 2021 John Wiley & Sons Ltd.)

The leaf economics spectrum 1,2 and the global spectrum of plant forms and functions 3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species 2 . Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities 4 . However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability 4,5 . Here we derive a set of ecosystem functions 6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems 7,8 ., (© 2021. The Author(s).)

The regional variability in tundra and boreal carbon dioxide (CO 2 ) 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) CO 2 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 CO 2 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 CO 2 fluxes and test the accuracy and uncertainty of different statistical models. CO 2 fluxes were upscaled at relatively high spatial resolution (1 km 2 ) 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 CO 2 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 CO 2 budget, estimated using the annual NEE ensemble prediction, suggests the high-latitude region was on average an annual CO 2 sink during 1990-2015, although uncertainty remains high., (© 2021 John Wiley & Sons Ltd.)

Climate change will impact forest productivity worldwide. Forecasting the magnitude of such impact, with multiple environmental stressors changing simultaneously, is only possible with the help of process-based models. In order to assess their performance, such models require careful evaluation against measurements. However, direct comparison of model outputs against observational data is often not reliable, as models may provide the right answers due to the wrong reasons. This would severely hinder forecasting abilities under unprecedented climate conditions. Here, we present a methodology for model assessment, which supplements the traditional output-to-observation model validation. It evaluates model performance through its ability to reproduce observed seasonal changes of the most limiting environmental driver (MLED) for a given process, here daily gross primary productivity (GPP). We analyzed seasonal changes of the MLED for GPP in two contrasting pine forests, the Mediterranean Pinus halepensis Mill. Yatir (Israel) and the boreal Pinus sylvestris L. Hyytiälä (Finland) from three years of eddy-covariance flux data. Then, we simulated the same period with a state-of-the-art process-based simulation model (LandscapeDNDC). Finally, we assessed if the model was able to reproduce both GPP observations and MLED seasonality. We found that the model reproduced the seasonality of GPP in both stands, but it was slightly overestimated without site-specific fine-tuning. Interestingly, although LandscapeDNDC properly captured the main MLED in Hyytiälä (temperature) and in Yatir (soil water availability), it failed to reproduce high-temperature and high-vapor pressure limitations of GPP in Yatir during spring and summer. We deduced that the most likely reason for this divergence is an incomplete description of stomatal behavior. In summary, this study validates the MLED approach as a model evaluation tool, and opens up new possibilities for model improvement., (© 2021 The Authors. Ecological Applications published by Wiley Periodicals LLC on behalf of Ecological Society of America.)

We apply and compare three widely applicable methods for estimating ecosystem transpiration (T) from eddy covariance (EC) data across 251 FLUXNET sites globally. All three methods are based on the coupled water and carbon relationship, but they differ in assumptions and parameterizations. Intercomparison of the three daily T estimates shows high correlation among methods (R between .89 and .94), but a spread in magnitudes of T/ET (evapotranspiration) from 45% to 77%. When compared at six sites with concurrent EC and sap flow measurements, all three EC-based T estimates show higher correlation to sap flow-based T than EC-based ET. The partitioning methods show expected tendencies of T/ET increasing with dryness (vapor pressure deficit and days since rain) and with leaf area index (LAI). Analysis of 140 sites with high-quality estimates for at least two continuous years shows that T/ET variability was 1.6 times higher across sites than across years. Spatial variability of T/ET was primarily driven by vegetation and soil characteristics (e.g., crop or grass designation, minimum annual LAI, soil coarse fragment volume) rather than climatic variables such as mean/standard deviation of temperature or precipitation. Overall, T and T/ET patterns are plausible and qualitatively consistent among the different water flux partitioning methods implying a significant advance made for estimating and understanding T globally, while the magnitudes remain uncertain. Our results represent the first extensive EC data-based estimates of ecosystem T permitting a data-driven perspective on the role of plants' water use for global water and carbon cycling in a changing climate., (© 2020 The Authors. Global Change Biology published by John Wiley & Sons Ltd.)

During the summer of 2018, a widespread drought developed over Northern and Central Europe. The increase in temperature and the reduction of soil moisture have influenced carbon dioxide (CO 2 ) exchange between the atmosphere and terrestrial ecosystems in various ways, such as a reduction of photosynthesis, changes in ecosystem respiration, or allowing more frequent fires. In this study, we characterize the resulting perturbation of the atmospheric CO 2 seasonal cycles. 2018 has a good coverage of European regions affected by drought, allowing the investigation of how ecosystem flux anomalies impacted spatial CO 2 gradients between stations. This density of stations is unprecedented compared to previous drought events in 2003 and 2015, particularly thanks to the deployment of the Integrated Carbon Observation System (ICOS) network of atmospheric greenhouse gas monitoring stations in recent years. Seasonal CO 2 cycles from 48 European stations were available for 2017 and 2018. Earlier data were retrieved for comparison from international databases or national networks. Here, we show that the usual summer minimum in CO 2 due to the surface carbon uptake was reduced by 1.4 ppm in 2018 for the 10 stations located in the area most affected by the temperature anomaly, mostly in Northern Europe. Notwithstanding, the CO 2 transition phases before and after July were slower in 2018 compared to 2017, suggesting an extension of the growing season, with either continued CO 2 uptake by photosynthesis and/or a reduction in respiration driven by the depletion of substrate for respiration inherited from the previous months due to the drought. For stations with sufficiently long time series, the CO 2 anomaly observed in 2018 was compared to previous European droughts in 2003 and 2015. Considering the areas most affected by the temperature anomalies, we found a higher CO 2 anomaly in 2003 (+3 ppm averaged over 4 sites), and a smaller anomaly in 2015 (+1 ppm averaged over 11 sites) compared to 2018. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

The role of plant phenology as a regulator for gross ecosystem productivity (GEP) in peatlands is empirically not well constrained. This is because proxies to track vegetation development with daily coverage at the ecosystem scale have only recently become available and the lack of such data has hampered the disentangling of biotic and abiotic effects. This study aimed at unraveling the mechanisms that regulate the seasonal variation in GEP across a network of eight European peatlands. Therefore, we described phenology with canopy greenness derived from digital repeat photography and disentangled the effects of radiation, temperature and phenology on GEP with commonality analysis and structural equation modeling. The resulting relational network could not only delineate direct effects but also accounted for possible effect combinations such as interdependencies (mediation) and interactions (moderation). We found that peatland GEP was controlled by the same mechanisms across all sites: phenology constituted a key predictor for the seasonal variation in GEP and further acted as a distinct mediator for temperature and radiation effects on GEP. In particular, the effect of air temperature on GEP was fully mediated through phenology, implying that direct temperature effects representing the thermoregulation of photosynthesis were negligible. The tight coupling between temperature, phenology and GEP applied especially to high latitude and high altitude peatlands and during phenological transition phases. Our study highlights the importance of phenological effects when evaluating the future response of peatland GEP to climate change. Climate change will affect peatland GEP especially through changing temperature patterns during plant phenologically sensitive phases in high latitude and high altitude regions., (© 2019 John Wiley & Sons Ltd.)

Solar-induced chlorophyll fluorescence (SIF) has been increasingly used as a proxy for terrestrial gross primary productivity (GPP). Previous work mainly evaluated the relationship between satellite-observed SIF and gridded GPP products both based on coarse spatial resolutions. Finer resolution SIF (1.3 km × 2.25 km) measured from the Orbiting Carbon Observatory-2 (OCO-2) provides the first opportunity to examine the SIF-GPP relationship at the ecosystem scale using flux tower GPP data. However, it remains unclear how strong the relationship is for each biome and whether a robust, universal relationship exists across a variety of biomes. Here we conducted the first global analysis of the relationship between OCO-2 SIF and tower GPP for a total of 64 flux sites across the globe encompassing eight major biomes. OCO-2 SIF showed strong correlations with tower GPP at both midday and daily timescales, with the strongest relationship observed for daily SIF at the 757 nm (R 2  = 0.72, p < 0.0001). Strong linear relationships between SIF and GPP were consistently found for all biomes (R 2  = 0.57-0.79, p < 0.0001) except evergreen broadleaf forests (R 2  = 0.16, p < 0.05) at the daily timescale. A higher slope was found for C 4 grasslands and croplands than for C 3 ecosystems. The generally consistent slope of the relationship among biomes suggests a nearly universal rather than biome-specific SIF-GPP relationship, and this finding is an important distinction and simplification compared to previous results. SIF was mainly driven by absorbed photosynthetically active radiation and was also influenced by environmental stresses (temperature and water stresses) that determine photosynthetic light use efficiency. OCO-2 SIF generally had a better performance for predicting GPP than satellite-derived vegetation indices and a light use efficiency model. The universal SIF-GPP relationship can potentially lead to more accurate GPP estimates regionally or globally. Our findings revealed the remarkable ability of finer resolution SIF observations from OCO-2 and other new or future missions (e.g., TROPOMI, FLEX) for estimating terrestrial photosynthesis across a wide variety of biomes and identified their potential and limitations for ecosystem functioning and carbon cycle studies., (© 2018 John Wiley & Sons Ltd.)

In this study latent heat flux (λE) measurements made at 65 boreal and arctic eddy-covariance (EC) sites were analyses by using the Penman-Monteith equation. Sites were stratified into nine different ecosystem types: harvested and burnt forest areas, pine forests, spruce or fir forests, Douglas-fir forests, broadleaf deciduous forests, larch forests, wetlands, tundra and natural grasslands. The Penman-Monteith equation was calibrated with variable surface resistances against half-hourly eddy-covariance data and clear differences between ecosystem types were observed. Based on the modeled behavior of surface and aerodynamic resistances, surface resistance tightly control λE in most mature forests, while it had less importance in ecosystems having shorter vegetation like young or recently harvested forests, grasslands, wetlands and tundra. The parameters of the Penman-Monteith equation were clearly different for winter and summer conditions, indicating that phenological effects on surface resistance are important. We also compared the simulated λE of different ecosystem types under meteorological conditions at one site. Values of λE varied between 15% and 38% of the net radiation in the simulations with mean ecosystem parameters. In general, the simulations suggest that λE is higher from forested ecosystems than from grasslands, wetlands or tundra-type ecosystems. Forests showed usually a tighter stomatal control of λE as indicated by a pronounced sensitivity of surface resistance to atmospheric vapor pressure deficit. Nevertheless, the surface resistance of forests was lower than for open vegetation types including wetlands. Tundra and wetlands had higher surface resistances, which were less sensitive to vapor pressure deficits. The results indicate that the variation in surface resistance within and between different vegetation types might play a significant role in energy exchange between terrestrial ecosystems and atmosphere. These results suggest the need to take into account vegetation type and phenology in energy exchange modeling., (© 2014 John Wiley & Sons Ltd.)