Most North American temperate forests are plantation or regrowth forests, which are actively managed. These forests are in different stages of their growth cycles and their ability to sequester atmospheric carbon is affected by extreme weather events. In this study, the impact of heat and drought events on carbon sequestration in an age-sequence (80, 45, and 17 years as of 2019) of eastern white pine (Pinus strobus L.) forests in southern Ontario, Canada was examined using eddy covariance flux measurements from 2003 to 2019.Over the 17-year study period, the mean annual values of net ecosystem productivity (NEP) were 180 ± 96, 538 ± 177 and 64 ± 165 g C m-2 yr-1 in the 80-, 45- and 17-year-old stands, respectively, with the highest annual carbon sequestration rate observed in the 45-year-old stand. We found that air temperature (Ta) was the dominant control on NEP in all three different-aged stands and drought, which was a limiting factor for both gross ecosystem productivity (GEP) and ecosystems respiration (RE), had a smaller impact on NEP. However, the simultaneous occurrence of heat and drought events during the early growing seasons or over the consecutive years had a significant negative impact on annual NEP in all three forests. We observed a similar trend of NEP decline in all three stands over three consecutive years that experienced extreme weather events, with 2016 being a hot and dry, 2017 being a dry, and 2018 being a hot year. The youngest stand became a net source of carbon for all three of these years and the oldest stand became a small source of carbon for the first time in 2018 since observations started in 2003. However, in 2019, all three stands reverted to annual net carbon sinks.Our study results indicate that the timing, frequency and concurrent or consecutive occurrence of extreme weather events may have significant implications for carbon sequestration in temperate conifer forests in Eastern North America. This study is one of few globally available to provide long-term observational data on carbon exchanges in different-aged temperate plantation forests. It highlights interannual variability in carbon fluxes and enhances our understanding of the responses of these forest ecosystems to extreme weather events. Study results will help in developing climate resilient and sustainable forestry practices to offset atmospheric greenhouse gas emissions and improving simulation of carbon exchange processes in terrestrial ecosystem models.

• A new PRI based LUE estimation method was proposed. • This method separated the half hourly observed PRI into PRI0 and ΔPRI. • PRI0 indicated daily maximal light use efficiency (LUE max ). • The ΔPRI linked PRI to different diurnal meteorological stress conditions. • This new method significantly improves LUE accuracy (R 2 from 0.1 to 0.7). Photosynthetic light use efficiency (LUE) determines the ability of a plant to assimilate atmospheric carbon dioxide to biomass and is known to be controlled by environmental conditions, light regimes and forest age. The photochemical reflectance index (PRI), derived from leaf or canopy remotely sensed spectra, has been shown to be an effective and accurate estimator of LUE. In this study, we propose a new LUE estimation method that separates the PRI into daily maximal PRI (PRI0) for indicating daily maximal light use efficiency (LUE max ) and ΔPRI, defined as the difference between PRI0 and instantaneous PRI, for estimating the diurnal physiological stress ( fstress) . We develop and apply the method across three temperate pine stands and a deciduous stand of different ages, in Southern Ontario, Canada. Half hourly canopy level spectra were acquired from a tower-based spectro-radiometer system (AMSPEC-III) over the growing season at the four stands. Results show that the PRI0 predicted well LUE max (R 2 > 0.6, p < 0.05) in both coniferous and deciduous stands and was able to track seasonal changes in pigment pools sizes. The ΔPRI was sensitive to short-term meteorological conditions, specifically temperature, vapor pressure deficit (VPD), and light variations resulting in strong correlations ( p < 0.05 ) with fstress and half hourly LUE. This new method significantly improves the estimation accuracy (R 2 increases from 0.1 to around 0.7) for PRI-based LUE estimation across all four stands of varying age and species composition and suggests that PRI-based LUE estimation has the ability to inform on both the effects of seasonal and diurnal change in photosynthetic efficiency under different meteorological conditions.

A Correction to this paper has been published: https://doi.org/10.1038/s41597-021-00851-9.

Abstract. Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land–atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80 % of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50 % of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56 % of the datasets. Many datasets contain data for species that make up 90 % or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (https://doi.org/10.5281/zenodo.3971689; Poyatos et al., 2020a). The “sapfluxnetr” R package – designed to access, visualize, and process SAPFLUXNET data – is available from CRAN.

• Phenological controls over aerodynamic resistance ( R ah ) were investigated. • R ah exhibits significant seasonal variability across a wide range of sites. • These shifts in R ah were caused by phenology in some ecosystems. • Accounting for variation in kB −1 is important for improving predictions of H . Surface roughness – a key control on land-atmosphere exchanges of heat and momentum – differs between dormant and growing seasons. However, how surface roughness shifts seasonally at fine time scales (e.g., days) in response to changing canopy conditions is not well understood. This study: (1) explores how aerodynamic resistance changes seasonally; (2) investigates what drives these seasonal shifts, including the role of vegetation phenology; and (3) quantifies the importance of including seasonal changes of aerodynamic resistance in “big leaf” models of sensible heat flux ( H ). We evaluated aerodynamic resistance and surface roughness lengths for momentum ( z 0m ) and heat ( z 0h ) using the kB −1 parameter (ln( z 0m / z 0h )). We used AmeriFlux data to obtain surface-roughness estimates, and PhenoCam greenness data for phenology. This analysis included 23 sites and ∼190 site years from deciduous broadleaf, evergreen needleleaf, woody savanna, cropland, grassland, and shrubland plant-functional types (PFTs). Results indicated clear seasonal patterns in aerodynamic resistance to sensible heat transfer ( R ah ). This seasonality tracked PhenoCam-derived start-of-season green-up transitions in PFTs displaying the most significant seasonal changes in canopy structure, with R ah decreasing near green-up transitions. Conversely, in woody savanna sites and evergreen needleleaf forests, patterns in R ah were not linked to green-up. Our findings highlight that decreases in kB −1 are an important control over R ah , explaining > 50% of seasonal variation in R ah across most sites. Decreases in kB −1 during green-up are likely caused by increasing z 0h in response to higher leaf area index. Accounting for seasonal variation in kB −1 is key for predicting H as well; assuming kB −1 to be constant resulted in significant biases that also exhibited strong seasonal patterns. Overall, we found that aerodynamic resistance can be sensitive to phenology in ecosystems having strong seasonality in leaf area, and this linkage is critical for understanding land-atmosphere interactions at seasonal time scales.

In eastern North America, many deciduous forest ecosystems grow at the northernmost extent of their geographical ranges, where climate change could aid or impede their growth. This region experiences frequent extreme weather conditions, allowing us to study the response of these forests to environmental conditions, reflective of future climates. Here we determined the impact of seasonal and annual climate variations and extreme weather events on the carbon (C) uptake capacity of an oak-dominated forest in southern Ontario, Canada, from 2012 to 2016. We found that changes in meteorology during late May to mid-July were key in determining the C sink strength of the forest, impacting the seasonal and annual variability of net ecosystem productivity (NEP). Overall, higher temperatures and dry conditions reduced ecosystem respiration (RE) much more than gross ecosystem productivity (GEP), leading to higher NEP. Variability in NEP was primarily driven by changes in RE, rather than GEP. The mean annual GEP, RE, and NEP values at our site during the study were 1,343 ± 85, 1,171 ± 139, and 206 ± 92 g C m-2 yr-1, respectively. The forest was a C sink even in years that experienced heat and water stresses. Mean annual NEP at our site was within the range of NEP (69-459 g C m-2 yr-1) observed in similar North American forests from 2012 to 2016. The growth and C sequestration capabilities of our oak-dominated forest were not adversely impacted by changes in environmental conditions and extreme weather events experienced over the study period.

Abstract The FLUXNET2015 dataset provides ecosystem-scale data on CO 2 , 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.

Abstract. The annual carbon and water dynamics of two eastern North American temperate forests were compared over a 6-year period from 2012 to 2017. The geographic location, forest age, soil, and climate were similar between the two stands; however, stand composition varied in terms of tree leaf-retention and shape strategy: one stand was a deciduous broadleaf forest, while the other was an evergreen needleleaf forest. The 6-year mean annual net ecosystem productivity (NEP) of the coniferous forest was slightly higher and more variable (218±109 g C m−2 yr−1) compared to that of the deciduous forest NEP (200±83 g C m−2 yr−1). Similarly, the 6-year mean annual evapotranspiration (ET) of the coniferous forest was higher (442±33 mm yr−1) than that of the deciduous forest (388±34 mm yr−1), but with similar interannual variability. Summer meteorology greatly impacted the carbon and water fluxes in both stands; however, the degree of response varied among the two stands. In general, warm temperatures caused higher ecosystem respiration (RE), resulting in reduced annual NEP values – an impact that was more pronounced at the deciduous broadleaf forest compared to the evergreen needleleaf forest. However, during warm and dry years, the evergreen forest had largely reduced annual NEP values compared to the deciduous forest. Variability in annual ET at both forests was related most to the variability in annual air temperature (Ta), with the largest annual ET observed in the warmest years in the deciduous forest. Additionally, ET was sensitive to prolonged dry periods that reduced ET at both stands, although the reduction at the coniferous forest was relatively larger than that of the deciduous forest. If prolonged periods (weeks to months) of increased Ta and reduced precipitation are to be expected under future climates during summer months in the study region, our findings suggest that the deciduous broadleaf forest will likely remain an annual carbon sink, while the carbon sink–source status of the coniferous forest remains uncertain.

Abstract This study presents the energy, water, and carbon (C) flux dynamics of a young afforested temperate white pine (Pinus strobus L.) forest in southern Ontario, Canada during the initial fourteen years (2003–2016) of establishment. Energy fluxes, namely, net radiation (Rn), latent heat (LE), and sensible heat (H) flux increased over time, due to canopy development. Annual values of ground heat flux (G) peaked in 2007 and then gradually declined in response to canopy closure. The forest became a consistent C-sink only 5 years after establishment owing in part to low respiratory fluxes from the former agricultural, sandy soils with low residual soil organic matter. Mean annual values of gross ecosystem productivity (GEP), ecosystem respiration (RE), and net ecosystem productivity (NEP) ranged from 494 to 1913, 515 to 1774 and −126 to 216 g C m−2 year−1 respectively, over the study period. Annual evapotranspiration (ET) values ranged from 328 to 429 mm year−1 over the same period. Water use efficiency (WUE) increased with stand age with a mean WUE value of 3.92 g C kg−1 H2O from 2008 to 2016. Multivariable linear regression analysis conducted using observed data suggested that the overall, C and water dynamics of the stand were primarily driven by radiation and temperature, both of which explained 77%, 48%, 28%, and 76% of the variability in GEP, RE, NEP, and ET, respectively. However, late summer droughts, which were prevalent in the region, reduced NEP. The reduction in NEP was enhanced when summer drought events were accompanied by increased heat such as those in 2005, 2012 and 2016. This study contributes to our understanding of the energy, water and C dynamics of afforested temperate conifer plantations and how these forests may respond to changing climate conditions during the crucial initial stage of their life cycle. Our findings also demonstrate the potential of pine plantation stands to sequester atmospheric CO2 in eastern North America.