Ecosystem warming extends vegetation activity but heightens vulnerability to cold temperatures.

Shifts in vegetation phenology are a key example of the biological effects of climate change1-3. However, there is substantial uncertainty about whether these temperature-driven trends will continue, or whether other factors-for example, photoperiod-will become more important as warming exceeds the bounds of historical variability4,5. Here we use phenological transition dates derived from digital repeat photography6 to show that experimental whole-ecosystem warming treatments7 of up to +9 °C linearly correlate with a delayed autumn green-down and advanced spring green-up of the dominant woody species in a boreal Picea-Sphagnum bog. Results were confirmed by direct observation of both vegetative and reproductive phenology of these and other bog plant species, and by multiple years of observations. There was little evidence that the observed responses were constrained by photoperiod. Our results indicate a likely extension of the period of vegetation activity by 1-2 weeks under a 'CO2 stabilization' climate scenario (+2.6 ± 0.7 °C), and 3-6 weeks under a 'high-CO2 emission' scenario (+5.9 ± 1.1 °C), by the end of the twenty-first century. We also observed severe tissue mortality in the warmest enclosures after a severe spring frost event. Failure to cue to photoperiod resulted in precocious green-up and a premature loss of frost hardiness8, which suggests that vulnerability to spring frost damage will increase in a warmer world9,10. Vegetation strategies that have evolved to balance tradeoffs associated with phenological temperature tracking may be optimal under historical climates, but these strategies may not be optimized for future climate regimes. These in situ experimental results are of particular importance because boreal forests have both a circumpolar distribution and a key role in the global carbon cycle11.

Evaluating the effects of moonlight on the vertical flight profiles of three western palaearctic swifts.

Recent studies have suggested the presence of moonlight mediated behaviour in avian aerial insectivores, such as swifts. Here, we use the combined analysis of state-of-the-art activity logger data across three swift species, the common, pallid and alpine swifts, to quantify flight height and activity in responses to moonlight-driven crepuscular and nocturnal light conditions. Our results show a significant response in flight heights to moonlight illuminance for common and pallid swifts, i.e. when moon illuminance increased flight height also increased, while a moonlight-driven response is absent in alpine swifts. We show a weak relationship between night-time illuminance-driven responses and twilight ascending behaviour, suggesting a decoupling of both crepuscular and night-time behaviour. We suggest that swifts optimize their flight behaviour to adapt to favourable night-time light conditions, driven by light-responsive and size-dependent vertical insect stratification and weather conditions.

Quantifying effects of cold acclimation and delayed springtime photosynthesis resumption in northern ecosystems.

Land carbon dynamics in temperate and boreal ecosystems are sensitive to environmental change. Accurately simulating gross primary productivity (GPP) and its seasonality is key for reliable carbon cycle projections. However, significant biases have been found in early spring GPP simulations of northern forests, where observations often suggest a later resumption of photosynthetic activity than predicted by models. Here, we used eddy covariance-based GPP estimates from 39 forest sites that differ by their climate and dominant plant functional types. We used a mechanistic and an empirical light use efficiency (LUE) model to investigate the magnitude and environmental controls of delayed springtime photosynthesis resumption (DSPR) across sites. We found DSPR reduced ecosystem LUE by 30-70% at many, but not all site-years during spring. A significant depression of LUE was found not only in coniferous but also at deciduous forests and was related to combined high radiation and low minimum temperatures. By embedding cold-acclimation effects on LUE that considers the delayed effects of minimum temperatures, initial model bias in simulated springtime GPP was effectively resolved. This provides an approach to improve GPP estimates by considering physiological acclimation and enables more reliable simulations of photosynthesis in northern forests and projections in a warming climate.

Pan-Arctic soil moisture control on tundra carbon sequestration and plant productivity.

Long-term atmospheric CO2 concentration records have suggested a reduction in the positive effect of warming on high-latitude carbon uptake since the 1990s. A variety of mechanisms have been proposed to explain the reduced net carbon sink of northern ecosystems with increased air temperature, including water stress on vegetation and increased respiration over recent decades. However, the lack of consistent long-term carbon flux and in situ soil moisture data has severely limited our ability to identify the mechanisms responsible for the recent reduced carbon sink strength. In this study, we used a record of nearly 100 site-years of eddy covariance data from 11 continuous permafrost tundra sites distributed across the circumpolar Arctic to test the temperature (expressed as growing degree days, GDD) responses of gross primary production (GPP), net ecosystem exchange (NEE), and ecosystem respiration (ER) at different periods of the summer (early, peak, and late summer) including dominant tundra vegetation classes (graminoids and mosses, and shrubs). We further tested GPP, NEE, and ER relationships with soil moisture and vapor pressure deficit to identify potential moisture limitations on plant productivity and net carbon exchange. Our results show a decrease in GPP with rising GDD during the peak summer (July) for both vegetation classes, and a significant relationship between the peak summer GPP and soil moisture after statistically controlling for GDD in a partial correlation analysis. These results suggest that tundra ecosystems might not benefit from increased temperature as much as suggested by several terrestrial biosphere models, if decreased soil moisture limits the peak summer plant productivity, reducing the ability of these ecosystems to sequester carbon during the summer.

Early spring onset increases carbon uptake more than late fall senescence: modeling future phenological change in a US northern deciduous forest.

In deciduous forests, spring leaf development and fall leaf senescence regulate the timing and duration of photosynthesis and transpiration. Being able to model these dates is therefore critical to accurately representing ecosystem processes in biogeochemical models. Despite this, there has been relatively little effort to improve internal phenology predictions in widely used biogeochemical models. Here, we optimized the phenology algorithms in a regionally developed biogeochemical model (PnET-CN) using phenology data from eight mid-latitude PhenoCam sites in eastern North America. We then performed a sensitivity analysis to determine how the optimization affected future predictions of carbon, water, and nitrogen cycling at Bartlett Experimental Forest, New Hampshire. Compared to the original PnET-CN phenology models, our new spring and fall models resulted in shorter season lengths and more abrupt transitions that were more representative of observations. The new phenology models affected daily estimates and interannual variability of modeled carbon exchange, but they did not have a large influence on the magnitude or long-term trends of annual totals. Under future climate projections, our new phenology models predict larger shifts in season length in the fall (1.1-3.2 days decade-1) compared to the spring (0.9-1.5 days decade-1). However, for every day the season was longer, spring had twice the effect on annual carbon and water exchange totals compared to the fall. These findings highlight the importance of accurately modeling season length for future projections of carbon and water cycling.

Developmental changes in the reflectance spectra of temperate deciduous tree leaves and implications for thermal emissivity and leaf temperature.

Leaf optical properties impact leaf energy balance and thus leaf temperature. The effect of leaf development on mid-infrared (MIR) reflectance, and hence thermal emissivity, has not been investigated in detail. We measured a suite of morphological characteristics, as well as directional-hemispherical reflectance from ultraviolet to thermal infrared wavelengths (250 nm to 20 µm) of leaves from five temperate deciduous tree species over the 8 wk following spring leaf emergence. By contrast to reflectance at shorter wavelengths, the shape and magnitude of MIR reflectance spectra changed markedly with development. MIR spectral differences among species became more pronounced and unique as leaves matured. Comparison of reflectance spectra of intact vs dried and ground leaves points to cuticular development - and not internal structural or biochemical changes - as the main driving factor. Accompanying the observed spectral changes was a drop in thermal emissivity from about 0.99 to 0.95 over the 8 wk following leaf emergence. Emissivity changes were not large enough to substantially influence leaf temperature, but they could potentially lead to a bias in radiometrically measured temperatures of up to 3 K. Our results also pointed to the potential for using MIR spectroscopy to better understand species-level differences in cuticular development and composition.

Seasonal variation in the canopy color of temperate evergreen conifer forests.

Evergreen conifer forests are the most prevalent land cover type in North America. Seasonal changes in the color of evergreen forest canopies have been documented with near-surface remote sensing, but the physiological mechanisms underlying these changes, and the implications for photosynthetic uptake, have not been fully elucidated. Here, we integrate on-the-ground phenological observations, leaf-level physiological measurements, near surface hyperspectral remote sensing and digital camera imagery, tower-based CO2 flux measurements, and a predictive model to simulate seasonal canopy color dynamics. We show that seasonal changes in canopy color occur independently of new leaf production, but track changes in chlorophyll fluorescence, the photochemical reflectance index, and leaf pigmentation. We demonstrate that at winter-dormant sites, seasonal changes in canopy color can be used to predict the onset of canopy-level photosynthesis in spring, and its cessation in autumn. Finally, we parameterize a simple temperature-based model to predict the seasonal cycle of canopy greenness, and we show that the model successfully simulates interannual variation in the timing of changes in canopy color. These results provide mechanistic insight into the factors driving seasonal changes in evergreen canopy color and provide opportunities to monitor and model seasonal variation in photosynthetic activity using color-based vegetation indices.

On quantifying the apparent temperature sensitivity of plant phenology.

Many plant phenological events are sensitive to temperature, leading to changes in the seasonal cycle of ecosystem function as the climate warms. To evaluate the current and future implications of temperature changes for plant phenology, researchers commonly use a metric of temperature sensitivity, which quantifies the change in phenology per degree change in temperature. Here, we examine the temperature sensitivity of phenology, and highlight conditions under which the widely used days-per-degree sensitivity approach is subject to methodological issues that can generate misleading results. We identify several factors, in particular the length of the period over which temperature is integrated, and changes in the statistical characteristics of the integrated temperature, that can affect the estimated apparent sensitivity to temperature. We show how the resulting artifacts can lead to spurious differences in apparent temperature sensitivity and artificial spatial gradients. Such issues are rarely considered in analyses of the temperature sensitivity of phenology. Given the issues identified, we advocate for process-oriented modelling approaches, informed by observations and with fully characterised uncertainties, as a more robust alternative to the simple days-per-degree temperature sensitivity metric. We also suggest approaches to minimise and assess spurious influences in the days-per-degree metric.

Asymmetric responses of ecosystem productivity to rainfall anomalies vary inversely with mean annual rainfall over the conterminous United States.

The CONterminous United States (CONUS) presents a large range of climate conditions and biomes where terrestrial primary productivity and its inter-annual variability are controlled regionally by rainfall and/or temperature. Here, the response of ecosystem productivity to those climate variables was investigated across different biomes from 2010 to 2018 using three climate datasets of precipitation, air temperature or drought severity, combined with several proxies of ecosystem productivity: a remote sensing product of aboveground biomass, an net primary productivity (NPP) remote sensing product, an NPP model-based product and four gross primary productivity products. We used an asymmetry index (AI) where positive AI indicates a greater increase of ecosystem productivity in wet years compared to the decline in dry years, and negative AI indicates a greater decline of ecosystem productivity in dry years compared to the increase in wet years. We found consistent spatial patterns of AI across the CONUS for the different products, with negative asymmetries over the Great Plains and positive asymmetries over the southwestern CONUS. Shrubs and, to a lesser extent, evergreen forests show a persistent positive asymmetry, whilst (natural) grasslands appear to have transitioned from positive to negative anomalies during the last decade. The general tendency of dominant negative asymmetry response for ecosystem productivity across the CONUS appears to be influenced by the negative asymmetry of precipitation anomalies. AI was found to be a function of mean rainfall: more positive AIs were found in dry areas where plants are adapted to drought and take advantage of rainfall pulses, and more negative AIs were found in wet areas, with a threshold delineating the two regimes corresponding to a mean annual rainfall of 200-400 mm/year.

Later springs green-up faster: the relation between onset and completion of green-up in deciduous forests of North America.

In deciduous forests, spring leaf phenology controls the onset of numerous ecosystem functions. While most studies have focused on a single annual spring event, such as budburst, ecosystem functions like photosynthesis and transpiration increase gradually after budburst, as leaves grow to their mature size. Here, we examine the "velocity of green-up," or duration between budburst and leaf maturity, in deciduous forest ecosystems of eastern North America. We use a diverse data set that includes 301 site-years of phenocam data across a range of sites, as well as 22 years of direct ground observations of individual trees and 3 years of fine-scale high-frequency aerial photography, both from Harvard Forest. We find a significant association between later start of spring and faster green-up: - 0.47 ± 0.04 (slope ± 1 SE) days change in length of green-up for every day later start of spring within phenocam sites, - 0.31 ± 0.06 days/day for trees under direct observation, and - 1.61 ± 0.08 days/day spatially across fine-scale landscape units. To explore the climatic drivers of spring leaf development, we fit degree-day models to the observational data from Harvard Forest. We find that the default phenology parameters of the ecosystem model PnET make biased predictions of leaf initiation (39 days early) and maturity (13 days late) for red oak, while the optimized model has biases of 1 day or less. Springtime productivity predictions using optimized parameters are closer to results driven by observational data (within 1%) than those of the default parameterization (17% difference). Our study advances empirical understanding of the link between early and late spring phenophases and demonstrates that accurately modeling these transitions is important for simulating seasonal variation in ecosystem productivity.

Greenness indices from digital cameras predict the timing and seasonal dynamics of canopy-scale photosynthesis.

The proliferation of digital cameras co-located with eddy covariance instrumentation provides new opportunities to better understand the relationship between canopy phenology and the seasonality of canopy photosynthesis. In this paper we analyze the abilities and limitations of canopy color metrics measured by digital repeat photography to track seasonal canopy development and photosynthesis, determine phenological transition dates, and estimate intra-annual and interannual variability in canopy photosynthesis. We used 59 site-years of camera imagery and net ecosystem exchange measurements from 17 towers spanning three plant functional types (deciduous broadleaf forest, evergreen needleleaf forest, and grassland/crops) to derive color indices and estimate gross primary productivity (GPP). GPP was strongly correlated with greenness derived from camera imagery in all three plant functional types. Specifically, the beginning of the photosynthetic period in deciduous broadleaf forest and grassland/crops and the end of the photosynthetic period in grassland/crops were both correlated with changes in greenness; changes in redness were correlated with the end of the photosynthetic period in deciduous broadleaf forest. However, it was not possible to accurately identify the beginning or ending of the photosynthetic period using camera greenness in evergreen needleleaf forest. At deciduous broadleaf sites, anomalies in integrated greenness and total GPP were significantly correlated up to 60 days after the mean onset date for the start of spring. More generally, results from this work demonstrate that digital repeat photography can be used to quantify both the duration of the photosynthetically active period as well as total GPP in deciduous broadleaf forest and grassland/crops, but that new and different approaches are required before comparable results can be achieved in evergreen needleleaf forest.

Historical tree phenology data reveal the seasonal rhythms of the Congo Basin rainforest.

Tropical forest phenology directly affects regional carbon cycles, but the relation between species-specific and whole-canopy phenology remains largely uncharacterized. We present a unique analysis of historical tropical tree phenology collected in the central Congo Basin, before large-scale impacts of human-induced climate change. Ground-based long-term (1937-1956) phenological observations of 140 tropical tree species are recovered, species-specific phenological patterns analyzed and related to historical meteorological records, and scaled to characterize stand-level canopy dynamics. High phenological variability within and across species and in climate-phenology relationships is observed. The onset of leaf phenophases in deciduous species was triggered by drought and light availability for a subset of species and showed a species-specific decoupling in time along a bi-modal seasonality. The majority of the species remain evergreen, although central African forests experience relatively low rainfall. Annually a maximum of 1.5% of the canopy is in leaf senescence or leaf turnover, with overall phenological variability dominated by a few deciduous species, while substantial variability is attributed to asynchronous events of large and/or abundant trees. Our results underscore the importance of accounting for constituent signals in canopy-wide scaling and the interpretation of remotely sensed phenology signals.

Acclimation of phenology relieves leaf longevity constraints in deciduous forests.

Leaf phenology is key for regulating total growing-season mass and energy fluxes. Long-term temporal trends towards earlier leaf unfolding are observed across Northern Hemisphere forests. Phenological dates also vary between years, whereby end-of-season (EOS) dates correlate positively with start-of-season (SOS) dates and negatively with growing-season total net CO2 assimilation (Anet). These associations have been interpreted as the effect of a constrained leaf longevity or of premature carbon (C) sink saturation-with far-reaching consequences for long-term phenology projections under climate change and rising CO2. Here, we use multidecadal ground and remote-sensing observations to show that the relationships between Anet and EOS are opposite at the interannual and the decadal time scales. A decadal trend towards later EOS persists in parallel with a trend towards increasing Anet-in spite of the negative Anet-EOS relationship at the interannual scale. This finding is robust against the use of diverse observations and models. Results indicate that acclimation of phenology has enabled plants to transcend a constrained leaf longevity or premature C sink saturation over the course of several decades, leading to a more effective use of available light and a sustained extension of the vegetation CO2 uptake season over time.

Using long-term data from a whole ecosystem warming experiment to identify best spring and autumn phenology models.

Predicting vegetation phenology in response to changing environmental factors is key in understanding feedbacks between the biosphere and the climate system. Experimental approaches extending the temperature range beyond historic climate variability provide a unique opportunity to identify model structures that are best suited to predicting phenological changes under future climate scenarios. Here, we model spring and autumn phenological transition dates obtained from digital repeat photography in a boreal Picea-Sphagnum bog in response to a gradient of whole ecosystem warming manipulations of up to +9°C, using five years of observational data. In spring, seven equally best-performing models for Larix utilized the accumulation of growing degree days as a common driver for temperature forcing. For Picea, the best two models were sequential models requiring winter chilling before spring forcing temperature is accumulated. In shrub, parallel models with chilling and forcing requirements occurring simultaneously were identified as the best models. Autumn models were substantially improved when a CO2 parameter was included. Overall, the combination of experimental manipulations and multiple years of observations combined with variation in weather provided the framework to rule out a large number of candidate models and to identify best spring and autumn models for each plant functional type.

Earlier snowmelt may lead to late season declines in plant productivity and carbon sequestration in Arctic tundra ecosystems.

Arctic warming is affecting snow cover and soil hydrology, with consequences for carbon sequestration in tundra ecosystems. The scarcity of observations in the Arctic has limited our understanding of the impact of covarying environmental drivers on the carbon balance of tundra ecosystems. In this study, we address some of these uncertainties through a novel record of 119 site-years of summer data from eddy covariance towers representing dominant tundra vegetation types located on continuous permafrost in the Arctic. Here we found that earlier snowmelt was associated with more tundra net CO2 sequestration and higher gross primary productivity (GPP) only in June and July, but with lower net carbon sequestration and lower GPP in August. Although higher evapotranspiration (ET) can result in soil drying with the progression of the summer, we did not find significantly lower soil moisture with earlier snowmelt, nor evidence that water stress affected GPP in the late growing season. Our results suggest that the expected increased CO2 sequestration arising from Arctic warming and the associated increase in growing season length may not materialize if tundra ecosystems are not able to continue sequestering CO2 later in the season.

Testing Hopkins' Bioclimatic Law with PhenoCam data.

PREMISE OF THE STUDY: We investigated the spatial and temporal patterns of vegetation phenology with phenometrics derived from PhenoCam imagery. Specifically, we evaluated the Bioclimatic Law proposed by Hopkins, which relates phenological transitions to latitude, longitude, and elevation. METHODS: "Green-up" and "green-down" dates-representing the start and end of the annual cycles of vegetation activity-were estimated from measures of canopy greenness calculated from digital repeat photography. We used data from 65 deciduous broadleaf (DB) forest sites, 18 evergreen needleleaf (EN) forest sites, and 21 grassland (GR) sites. RESULTS: DB green-up dates were well correlated with mean annual temperature and varied along spatial gradients consistent with the Bioclimatic Law. Interannual variation in DB phenology was most strongly associated with temperature anomalies during a relatively narrow window of time. EN phenology was not well correlated with either climatic factors or spatial gradients, but similar to DB phenology, interannual variation was most closely associated with temperature anomalies. For GR sites, mean annual precipitation explained most of the spatial variation in the duration of vegetation activity, whereas both temperature and precipitation anomalies explained interannual variation in phenology. DISCUSSION: PhenoCam data provide an objective and consistent means by which spatial and temporal patterns in vegetation phenology can be investigated.

Large-sized rare tree species contribute disproportionately to functional diversity in resource acquisition in African tropical forest.

Increasing evidence is available for a positive effect of biodiversity on ecosystem productivity and standing biomass, also in highly diverse systems as tropical forests. Biodiversity conservation could therefore be a critical aspect of climate mitigation policies. There is, however, limited understanding of the role of individual species for this relationship, which could aid in focusing conservation efforts and forest management planning. This study characterizes the functional specialization and redundancy for 95% of all tree species (basal area weighted percentage) in a diverse tropical forest in the central Congo Basin and relates this to species' abundance, contribution to aboveground carbon, and maximum size. Functional characterization is based on a set of traits related to resource acquisition (wood density, specific leaf area, leaf carbon, nitrogen and phosphorus content, and leaf stable carbon isotope composition). We show that within both mixed and monodominant tropical forest ecosystems, the highest functional specialization and lowest functional redundancy are solely found in rare tree species and significantly more in rare species holding large-sized individuals. Rare species cover the entire range of low and high functional redundancy, contributing both unique and redundant functions. Loss of species supporting functional redundancy could be buffered by other species in the community, including more abundant species. This is not the case for species supporting high functional specialization and low functional redundancy, which would need specific conservation attention. In terms of tropical forest management planning, we argue that specific conservation of large-sized trees is imperative for long-term maintenance of ecosystem functioning.

Publisher Correction: Tracking vegetation phenology across diverse biomes using Version 2.0 of the PhenoCam Dataset.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Tracking vegetation phenology across diverse biomes using Version 2.0 of the PhenoCam Dataset.

Monitoring vegetation phenology is critical for quantifying climate change impacts on ecosystems. We present an extensive dataset of 1783 site-years of phenological data derived from PhenoCam network imagery from 393 digital cameras, situated from tropics to tundra across a wide range of plant functional types, biomes, and climates. Most cameras are located in North America. Every half hour, cameras upload images to the PhenoCam server. Images are displayed in near-real time and provisional data products, including timeseries of the Green Chromatic Coordinate (Gcc), are made publicly available through the project web page ( https://phenocam.sr.unh.edu/webcam/gallery/ ). Processing is conducted separately for each plant functional type in the camera field of view. The PhenoCam Dataset v2.0, described here, has been fully processed and curated, including outlier detection and expert inspection, to ensure high quality data. This dataset can be used to validate satellite data products, to evaluate predictions of land surface models, to interpret the seasonality of ecosystem-scale CO2 and H2O flux data, and to study climate change impacts on the terrestrial biosphere.

Satellite-observed pantropical carbon dynamics.

Changes in terrestrial tropical carbon stocks have an important role in the global carbon budget. However, current observational tools do not allow accurate and large-scale monitoring of the spatial distribution and dynamics of carbon stocks1. Here, we used low-frequency L-band passive microwave observations to compute a direct and spatially explicit quantification of annual aboveground carbon (AGC) fluxes and show that the tropical net AGC budget was approximately in balance during 2010 to 2017, the net budget being composed of gross losses of -2.86 PgC yr-1 offset by gross gains of -2.97 PgC yr-1 between continents. Large interannual and spatial fluctuations of tropical AGC were quantified during the wet 2011 La Niña year and throughout the extreme dry and warm 2015-2016 El Niño episode. These interannual fluctuations, controlled predominantly by semiarid biomes, were shown to be closely related to independent global atmospheric CO2 growth-rate anomalies (Pearson's r = 0.86), highlighting the pivotal role of tropical AGC in the global carbon budget.