Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Contribution of deep soil layers to the transpiration of a temperate deciduous forest: Implications for the modelling of productivity

Climate change is imposing drier atmospheric and edaphic conditions on temperate forests. Here, we investigated how deep soil (down to 300 cm) water extraction contributed to the provision of water in the Fontainebleau-Barbeau temperate oak forest over two years, including the 2018 record drought. Deep water provision was key to sustain canopy transpiration during drought, with layers below 150 cm contributing up to 60% of the transpired water in August 2018, despite their very low density of fine roots. We further showed that soil databases used to parameterize ecosystem models largely underestimated the amount of water extractable from the soil by trees, due to a considerable underestimation of the tree rooting depth. The consensus database established for France gave an estimate of 207 mm for the soil water holding capacity (SWHC) at Fontainebleau-Barbeau, when our estimate based on the analysis of soil water content measurements was 1.9 times as high, reaching 390 ± 17 mm. Running the CASTANEA forest model with the database-derived SWHC yielded a 185 gC m^-2 y^-1 average underestimation of annual gross primary productivity under current climate, reaching up to 687 ± 117 gC m^-2 y^-1 under climate change scenario RCP8.5. It is likely that the strong underestimation of SWHC that we show at our site is not a special case, and concerns a large number of forest sites. Thus, we argue for a generalisation of deep soil water content measurements in forests, in order to improve the estimation of SWHC and the simulation of the forest carbon cycle in the current context of climate change.

Limits to reproduction and seed size-number trade-offs that shape forest dominance and future recovery

The relationships that control seed production in trees are fundamental to understanding the evolution of forest species and their capacity to recover from increasing losses to drought, fire, and harvest. A synthesis of fecundity data from 714 species worldwide allowed us to examine hypotheses that are central to quantifying reproduction, a foundation for assessing fitness in forest trees. Four major findings emerged. First, seed production is not constrained by a strict trade-off between seed size and numbers. Instead, seed numbers vary over ten orders of magnitude, with species that invest in large seeds producing more seeds than expected from the 1:1 trade-off. Second, gymnosperms have lower seed production than angiosperms, potentially due to their extra investments in protective woody cones. Third, nutrient-demanding species, indicated by high foliar phosphorus concentrations, have low seed production. Finally, sensitivity of individual species to soil fertility varies widely, limiting the response of community seed production to fertility gradients. In combination, these findings can inform models of forest response that need to incorporate reproductive potential.

Globally, tree fecundity exceeds productivity gradients

Lack of tree fecundity data across climatic gradients precludes the analysis of how seed supply contributes to global variation in forest regeneration and biotic interactions responsible for biodiversity. A global synthesis of raw seedproduction data shows a 250-fold increase in seed abundance from cold-dry to warm-wet climates, driven primarily by a 100-fold increase in seed production for a given tree size. The modest (threefold) increase in forest productivity across the same climate gradient cannot explain the magnitudes of these trends. The increase in seeds per tree can arise from adaptive evolution driven by intense species interactions or from the direct effects of a warm, moist climate on tree fecundity. Either way, the massive differences in seed supply ramify through food webs potentially explaining a disproportionate role for species interactions in the wet tropics.

The 2018 European heatwave led to stem dehydration but not to consistent growth reductions in forests

Heatwaves exert disproportionately strong and sometimes irreversible impacts on forest ecosystems. These impacts remain poorly understood at the tree and species level and across large spatial scales. Here, we investigate the effects of the record-breaking 2018 European heatwave on tree growth and tree water status using a collection of high-temporal resolution dendrometer data from 21 species across 53 sites. Relative to the two preceding years, annual stem growth was not consistently reduced by the 2018 heatwave but stems experienced twice the temporary shrinkage due to depletion of water reserves. Conifer species were less capable of rehydrating overnight than broadleaves across gradients of soil and atmospheric drought, suggesting less resilience toward transient stress. In particular, Norway spruce and Scots pine experienced extensive stem dehydration. Our high-resolution dendrometer network was suitable to disentangle the effects of a severe heatwave on tree growth and desiccation at large-spatial scales in situ, and provided insights on which species may be more vulnerable to climate extremes.

The within-population variability of leaf spring and autumn phenology is influenced by temperature in temperate deciduous trees

Leaf phenology is a major driver of ecosystem functioning in temperate forests and a robust indicator of climate change. Both the inter-annual and inter-population variability of leaf phenology have received much attention in the literature; in contrast, the within-population variability of leaf phenology has been far less studied. Beyond its impact on individual tree physiological processes, the within-population variability of leaf phenology can affect the estimation of the average budburst or leaf senescence dates at the population scale. Here, we monitored the progress of spring and autumn leaf phenology over 14 tree populations (9 tree species) in six European forests over the period of 2011 to 2018 (yielding 16 site-years of data for spring, 14 for autumn). We monitored 27 to 512 (with a median of 62) individuals per population. We quantified the within-population variability of leaf phenology as the standard deviation of the distribution of individual dates of budburst or leaf senescence (SD_BBi and SD_LSi, respectively). Given the natural variability of phenological dates occurring in our tree populations, we estimated from the data that a minimum sample size of 28 (resp. 23) individuals, are required to estimate SD_BBi (resp. SD_LSi) with a precision of 3 (resp. 7) days. The within-population of leaf senescence (average SD_LSi = 8.5 days) was on average two times larger than for budburst (average SD_BBi = 4.0 days). We evidenced that warmer temperature during the budburst period and a late average budburst date were associated with a lower SD_BBi, as a result of a quicker spread of budburst in tree populations, with a strong species effect. Regarding autumn phenology, we observed that later senescence and warm temperatures during the senescence period were linked with a high SD_LSi, with a strong species effect. The shares of variance explained by our models were modest suggesting that other factors likely influence the within-population variation in leaf phenology. For instance, a detailed analysis revealed that summer temperatures were negatively correlated with a lower SD_LSi.

Ecosystem transpiration and evaporation: Insights from three water flux partitioning methods across FLUXNET sites

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.

The FLUXNET2015 dataset and the ONEFlux processing pipeline for eddy covariance data

The FLUXNET2015 dataset provides ecosystem-scale data on CO₂, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Measurement(s)	net ecosystem exchange • carbon dioxide • water • energy
Technology Type(s)	eddy covariance • measurement device
Sample Characteristic - Environment	terrestrial biome • atmosphere
Sample Characteristic - Location	Earth (planet)

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12295910

"Green pointillism": detecting the within-population variability of budburst in temperate deciduous trees with phenological cameras

Phenological cameras have been used over a decade for identifying plant phenological markers (budburst, leaf senescence) and more generally the greenness dynamics of forest canopies. The analysis is usually carried out over the full camera field of view, with no particular analysis of the variability of phenological markers among trees. Here we show that images produced by phenological cameras can be used to quantify the within-population variability of budburst (WPVbb) in temperate deciduous forests. Using seven site-years of image analyses, we report a strong correlation (r² = 0.97) between the WPVbb determined with a phenological camera and its quantification through ground observation. We show that WPVbb varies strongly (by a factor of 4) from year to year in a given population and that those variations are linked with temperature conditions during the budburst period, with colder springs associated to a higher differentiation of budburst (higher WPVbb) among trees. Deploying our approach at the continental scale, i.e., throughout phenological cameras networks, would improve the understanding of the spatial (across populations) and temporal (across years) variations of WPVbb, which have strong implications on forest functioning, tree fitness and phenological modelling.

Carbon-use strategies in stem radial growth of two oak species, one Temperate deciduous and one Mediterranean evergreen: what can be inferred from seasonal variations in the δ13C of the current year ring?

Tree ring synthesis is a key process in wood production; however, little is known of the origin and fate of the carbon involved. We used natural 13C abundance to investigate the carbon-use process for the ring development in a temperate deciduous (Quercus petraea (Matt.) Liebl.) and a Mediterranean evergreen (Quercus ilex L.) oak. The sapwood carbon reserves, phloem sucrose contents, stem respired CO2 efflux and their respective carbon isotope compositions (δ13C) were recorded over 1 year, in the native area of each species. The seasonal δ13C variation of the current year ring was determined in the total ring throughout the seasons, as well as in slices from the fully mature ring after the growth season (intra-ring pattern). Although the budburst dates of the two oaks were similar, the growth of Quercus ilex began 50 days later. Both species exhibited growth cessation during the hot and dry summer but only Q. ilex resumed in the autumn. In the deciduous oak, xylem starch storage showed clear variations during the radial growth. The intra-ring δ13C variations of the two species exhibited similar ranges, but contrasting patterns, with an early increase for Q. petraea. Comparison between δ13C of starch and total ring suggested that Q. petraea (but not Q. ilex) builds its rings using reserves during the first month of growth. Shifts in ring and soluble sugars δ13C suggested an interspecific difference in either the phloem unloading or the use of fresh assimilate inside the ring. A decrease in ring δ13C for both oaks between the end of the radial growth and the winter is attributed to a lignification of ring cell walls after stem increment. This study highlighted the differences in carbon-use during ring growth for evergreen and deciduous oaks, as well as the benefits of exploring the process using natural 13C abundance.

Wood phenology, not carbon input, controls the interannual variability of wood growth in a temperate oak forest

Although the analysis of flux data has increased our understanding of the interannual variability of carbon inputs into forest ecosystems, we still know little about the determinants of wood growth. Here, we aimed to identify which drivers control the interannual variability of wood growth in a mesic temperate deciduous forest. We analysed a 9-yr time series of carbon fluxes and aboveground wood growth (AWG), reconstructed at a weekly time-scale through the combination of dendrometer and wood density data. Carbon inputs and AWG anomalies appeared to be uncorrelated from the seasonal to interannual scales. More than 90% of the interannual variability of AWG was explained by a combination of the growth intensity during a first 'critical period' of the wood growing season, occurring close to the seasonal maximum, and the timing of the first summer growth halt. Both atmospheric and soil water stress exerted a strong control on the interannual variability of AWG at the study site, despite its mesic conditions, whilst not affecting carbon inputs. Carbon sink activity, not carbon inputs, determined the interannual variations in wood growth at the study site. Our results provide a functional understanding of the dependence of radial growth on precipitation observed in dendrological studies.

In situ assessment of the velocity of carbon transfer by tracing 13 C in trunk CO2 efflux after pulse labelling: variations among tree species and seasons

Phloem is the main pathway for transferring photosynthates belowground. In situ(13) C pulse labelling of trees 8-10 m tall was conducted in the field on 10 beech (Fagus sylvatica) trees, six sessile oak (Quercus petraea) trees and 10 maritime pine (Pinus pinaster) trees throughout the growing season. Respired (13) CO2 from trunks was tracked at different heights using tunable diode laser absorption spectrometry to determine time lags and the velocity of carbon transfer (V). The isotope composition of phloem extracts was measured on several occasions after labelling and used to estimate the rate constant of phloem sap outflux (kP ). Pulse labelling together with high-frequency measurement of the isotope composition of trunk CO2 efflux is a promising tool for studying phloem transport in the field. Seasonal variability in V was predicted in pine and oak by bivariate linear regressions with air temperature and soil water content. V differed among the three species consistently with known differences in phloem anatomy between broadleaf and coniferous trees. V increased with tree diameter in oak and beech, reflecting a nonlinear increase in volumetric flow with increasing bark cross-sectional area, which suggests changes in allocation pattern with tree diameter in broadleaf species. Discrepancies between V and kP indicate vertical changes in functional phloem properties.

Impact of carbohydrate supply on stem growth, wood and respired CO2 delta13C: assessment by experimental girdling

The present study examines the impact of the C source (reserves vs current assimilates) on tree C isotope signals and stem growth, using experimental girdling to stop the supply of C from leaves to stem. Two-year-old sessile oaks (Quercus petraea) were girdled at three different phenological periods during the leafy period: during early wood growth (Girdling Period 1), during late wood growth (Girdling Period 2) and just after growth cessation (Girdling Period 3). The measured variables included stem respiration rates, stem radial increment, delta(13)C of respired CO(2) and contents of starch and water-soluble fraction in stems (below the girdle) and leaves. Girdling stopped growth, even early in the growing season, leading to a decrease in stem CO(2) efflux (CO(2R)). Shift in substrate use from recently fixed carbohydrate to reserves (i.e., starch) induced (13)C enrichment of CO(2) respired by stem. However, change in substrate type was insufficient to explain alone all the observed CO(2R) delta(13)C variations, especially at the period corresponding to large growth rate of control trees. The below-girdle mass balance suggested that, during girdling periods, stem C was invested in metabolic pathways other than respiration and stem growth. After Girdling Period 1, the girdle healed and the effects of girdling on stem respiration were reversed. Stem growth restarted and total radial increment was similar to the control one, indicating that growth can be delayed when a stress event occurs early in the growth period. Concerning tree ring, seasonal shift in substrate use from reserves (i.e., starch) to recently fixed carbohydrate is sufficient to explain the observed (13)C depletion of tree ring during the early wood growth. However, the inter-tree intra-ring delta(13)C variability needs to be resolved in order to improve the interpretation of intra-seasonal ring signals in terms of climatic or ecophysiological information. This study highlighted, via carbohydrate availability effects, the importance of the characterization of stem metabolic pathways for a complete understanding of the delta(13)C signals.

Interspecific variability of stem photosynthesis among tree species

The photosynthetic characteristics of current-year stems of six deciduous tree species, two evergreen tree species and ginkgo (Ginkgo biloba L.) were compared. Gas exchange, chlorophyll concentration, nitrogen concentration and maximum quantum yield of PSII were measured in stems in summer and winter. A light-induced decrease in stem CO(2) efflux was observed in all species. The apparent gross photosynthetic rate in saturating light ranged from 0.72 micromol m(-2 )s(-1) (ginkgo, in winter) to 3.73 micromol m(-2) s(-1) (Alnus glutinosa (L.) Gaertn., in summer). Despite this variability, a unique correlation (slope = 0.75), based on our results and those reported in the literature, was found between gross photosynthetic rate and dark respiration rate. Mass-based gross photosynthetic rate decreased with stem mass per area and correlated to chlorophyll concentration and nitrogen concentration, both in summer and winter. The radial distribution of stem chlorophyll differed among species, but all species except ginkgo had chlorophyll as deep as the pith. In summer, the maximum quantum yield of stem PSII (estimated from the ratio of variable to maximal fluorescence; F (v)/F (m)) of all species was near the optimal value found for leaves. By contrast, the values were highly variable in winter, suggesting large differences in sensitivity to low-temperature photoinhibition. The winter values of Fv/Fm were only 31-60% of summer values for the deciduous species, whereas the evergreen conifer species maintained high F (v)/F (m) in winter. The results highlight the interspecific variability of gross photosynthesis in the stem and its correlation with structural traits like those found for leaves. The structural correlations suggest that the selection of photosynthetic traits has operated under similar constraints in stems and leaves.

Tree stem phosphoenolpyruvate carboxylase (PEPC): lack of biochemical and localization evidence for a C4-like photosynthesis system

Here, the kinetic properties and immunolocalization of phosphoenolpyruvate carboxylase (PEPC) and ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) in young stems of Fagus sylvatica were investigated. The aim of the study was to test the hypothesis that there is a C4-like photosynthesis system in the stems of this C3 tree species. The activity, optimal pH and L-malate sensitivity of PEPC, and the Michaelis-Menten constant (Km) for phosphoenolpyruvate (PEP), were measured in protein extracts from current-year stems and leaves. A gel blot experiment and immunolocalization studies were performed to examine the isozyme complexity of PEPC and the tissue distribution of PEPC and Rubisco in stems. Leaf and stem PEPCs exhibited similar, classical values characteristic of C3 PEPCs, with an optimal pH of c. 7.8, a Km for PEP of c. 0.3 mM and a IC50 for L-malate (the L-malate concentration that inhibits 50% of PEPC activity at the Km for PEP) of c. 0.1 mM. Western blot analysis showed the presence of two PEPC subunits (molecular mass c. 110 kDa) both in leaves and in stems. Immunogold labelling did not reveal any differential localization of PEPC and Rubisco, neither between nor inside cells. This study suggests that C4-type photosynthesis does not occur in stems of F. sylvatica and underlines the importance of PEPC in nonphotosynthetic carbon fixation by most stem tissues (fixation of respired CO2 and fixation via the anaplerotic pathway).

Seasonal, daily and diurnal variations in the stable carbon isotope composition of carbon dioxide respired by tree trunks in a deciduous oak forest

The stable C isotope composition (delta13C) of CO2 respired by trunks was examined in a mature temperate deciduous oak forest (Quercus petraea). Month-to-month, day-to-day and diurnal, measurements were made to determine the range of variations at different temporal scales. Trunk growth and respiration rates were assessed. Phloem tissue was sampled and was analysed for total organic matter and soluble sugar 13C composition. The CO2 respired by trunk was always enriched in 13C relative to the total organic matter, sometimes by as much as 5 per thousand. The delta13C of respired CO2 exhibited a large seasonal variation (3.3 per thousand), with a relative maximum at the beginning of the growth period. The lowest values occurred in summer when the respiration rates were maximal. After the cessation of radial trunk growth, the respired CO2 delta13C values showed a progressive increase, which was linked to a parallel increase in soluble sugar content in the phloem tissue (R=0.95; P<0.01). At the same time, the respiration rates declined. This limited use of the substrate pool might allow the discrimination during respiration to be more strongly expressed. The late-season increase in CO2 delta13C might also be linked to a shift from recently assimilated C to reserves. At the seasonal scale, CO2 delta13C was negatively correlated with air temperature (R=-0.80; P<0.01). The diurnal variation sometimes reached 3 per thousand, but the range and the pattern depended on the period within the growing season. Contrary to expectations, diurnal variations were maximal in winter and spring when the leaves were missing or not totally functional. By contrast to the seasonal scale, these diurnal variations were not related to air temperature or sugar content. Our study shows that seasonal and diurnal variations of respired 13C exhibited a similar large range but were probably explained by different mechanisms.

Mixotrophy in orchids: insights from a comparative study of green individuals and nonphotosynthetic individuals of Cephalanthera damasonium

Some green orchids obtain carbon (C) from their mycorrhizal fungi and photosynthesis. This mixotrophy may represent an evolutionary step towards mycoheterotrophic plants fully feeding on fungal C. Here, we report on nonphotosynthetic individuals (albinos) of the green Cephalanthera damasonium that likely represent another evolutionary step. Albino and green individuals from a French population were compared for morphology and fertility, photosynthetic abilities, fungal partners (using microscopy and molecular tools), and nutrient sources (as characterized by 15N and 13C abundances). Albinos did not differ significantly from green individuals in morphology and fertility, but tended to be smaller. They harboured similar fungi, with Thelephoraceae and Cortinariaceae as mycorrhizal partners and few rhizoctonias. Albinos were nonphotosynthetic, fully mycoheterotrophic. Green individuals carried out photosynthesis at compensation point and received almost 50% of their C from fungi. Orchid fungi also colonized surrounding tree roots, likely to be the ultimate C source. Transition to mycoheterotrophy may require several simultaneous adaptations; albinos, by lacking some of them, may have reduced ecological success. This may limit the appearance of cheaters in mycorrhizal networks.

The carbon isotope composition of CO2 respired by trunks: comparison of four sampling methods

The (13)C natural abundance of CO(2) respired by plants has been used in the laboratory to examine the discrimination processes that occur during respiration. Currently, field measurements are being expanded to interpret the respiration delta(13)C signature measured at ecosystem and global levels. In this context, forests are particularly important to consider as they represent 80% of the continental biomass. The objective of this investigation was to compare four methods of sampling the CO(2) respired by trunks for the determination of its carbon isotope composition: three in situ methods using chambers placed on the trunk, and one destructive method using cores of woody tissues. The in situ methods were based either on a Keeling plot approach applied at the tissue level or on an initial flush of the chamber with nitrogen or with CO(2)-free air. In parallel, we investigated the possibility of an apparent discrimination during tissue respiration by comparing the delta(13)C signature of the respired CO(2) and that of the organic matter. The study was performed on six tree species widely distributed in temperate and mediterranean areas. The four methods were not significantly different when overall means were considered. However, considering the individual data, the Keeling plot approach and the nitrogen flush methods gave fairly homogeneous results, whereas the CO(2)-free air method produced more variable results. The core method was not correlated with any of the chamber methods. Regardless of the methodology, the respired CO(2) generally was enriched in (13)C relative to the total organic matter. This apparent enrichment during respiration was variable, reaching as much as 3-5 per thousand. This study showed that, on the whole, the different sampling techniques gave similar results, but one should be aware of the variability associated with each method.