Person:

Richardson, Andrew

Climate change affects the phenology of many species. As temperature and precipitation are thought to control autumn color change in temperate deciduous trees, it is possible that climate change might also affect the phenology of autumn colors. Using long-term data for eight tree species in a New England hardwood forest, we show that the timing and cumulative amount of autumn color are correlated with variation in temperature and precipitation at specific times of the year. A phenological model driven by accumulated cold degree-days and photoperiod reproduces most of the interspecific and interannual variability in the timing of autumn colors. We use this process-oriented model to predict changes in the phenology of autumn colors to 2099, showing that, while responses vary among species, climate change under standard IPCC projections will lead to an overall increase in the amount of autumn colors for most species.

Primarily driven by concern about rising levels of atmospheric CO2, ecologists and earth system scientists are collecting vast amounts of data related to the carbon cycle. These measurements are generally time-consuming and expensive to make, and, unfortunately, we live in an era where research funding is increasingly hard to come by. Thus, important questions are: ‘Which data streams provide the most valuable information? ’ and, ‘How much data do we need? ’ These questions are relevant not only for model developers, who need observational data to improve, constrain and test their models, but also for experimentalists and those designing ecological observation networks. Here we address these questions using a model-data fusion approach. We constrain a process-oriented, forest ecosystem C cycle model with seventeen different data streams from the Harvard Forest. We iteratively rank each data source according to its contribution to reducing model uncertainty. Results show the importance of some measurements commonly unavailable to carbon cycle modelers, such as estimates of turnover times from different carbon pools. Surprisingly, many data sources are relatively redundant in the presence of others, and do not lead to a significant improvement in model performance. A few select data sources lead to the largest reduction in parameter based model uncertainty. Projections of future carbon cycling were poorly constrained when only hourly net ecosystem exchange measurements were used to inform the model. They were well constrained, however, with only five of the seventeen data streams, even though many individual parameters are not constrained. The approach taken here should stimulate further cooperation between modelers and measurement teams, and may be useful in the context of setting research priorities and allocating research funds.

Interannual variability in biosphere-atmosphere exchange of CO2 is driven by a diverse range of biotic and abiotic factors. Replicating this variability thus represents the ‘acid test’ for terrestrial biosphere models. Although such models are commonly used to project responses to both normal and anomalous variability in climate, they are rarely tested explicitly against inter-annual variability in observations. Herein, using standardized data from the North American Carbon Program, we assess the performance of 16 terrestrial biosphere models and 3 remote sensing products against long-term measurements of biosphere-atmosphere CO2 exchange made with eddy-covariance flux towers at 11 forested sites in North America. Instead of focusing on model-data agreement we take a systematic, variability-oriented approach and show that although the models tend to reproduce the mean magnitude of the observed annual flux variability, they fail to reproduce the timing. Large biases in modeled annual means are evident for all models. Observed interannual variability is found to commonly be on the order of magnitude of the mean fluxes. None of the models consistently reproduce observed interannual variability within measurement uncertainty. Underrepresentation of variability in spring phenology, soil thaw and snowpack melting, and difficulties in reproducing the lagged response to extreme climatic events are identified as systematic errors, common to all models included in this study.

Uncertainties in model projections of carbon cycling in terrestrial ecosystems stem from inaccurate parameterization of incorporated processes (endogenous uncertainties) and processes or drivers that are not accounted for by the model (exogenous uncertainties). Here, we assess endogenous and exogenous uncertainties using a model-data fusion framework benchmarked with an artificial neural network (ANN). We used 18 years of eddy-covariance carbon flux data from the Harvard forest, where ecosystem carbon uptake has doubled over the measurement period, along with 15 ancillary ecological data sets relative to the carbon cycle. We test the ability of combinations of diverse data to constrain projections of a process-based carbon cycle model, both against the measured decadal trend and under future long-term climate change. The use of high-frequency eddy-covariance data alone is shown to be insufficient to constrain model projections at the annual or longer time step. Future projections of carbon cycling under climate change in particular are shown to be highly dependent on the data used to constrain the model. Endogenous uncertainties in long-term model projections of future carbon stocks and fluxes were greatly reduced by the use of aggregated flux budgets in conjunction with ancillary data sets. The data-informed model, however, poorly reproduced interannual variability in net ecosystem carbon exchange and biomass increments and did not reproduce the long-term trend. Furthermore, we use the model-data fusion framework, and the ANN, to show that the long-term doubling of the rate of carbon uptake at Harvard forest cannot be explained by meteorological drivers, and is driven by changes during the growing season. By integrating all available data with the model-data fusion framework, we show that the observed trend can only be reproduced with temporal changes in model parameters. Together, the results show that exogenous uncertainty dominates uncertainty in future projections from a data-informed process-based model.

Forests dominate the global carbon cycle, but their role in methane (CH4) biogeochemistry remains uncertain. We analyzed whole-ecosystem CH4 fluxes from 2 years, obtained over a lowland evergreen forest in Maine, USA. Gross primary productivity provided the strongest correlation with the CH4 flux in both years, with an additional significant effect of soil moisture in the second, drier year. This forest was a neutral to net source of CH4 in 2011 and a small net sink in 2012. Interannual variability in the summer hydrologic cycle apparently shifts the ecosystem from being a net source to a sink for CH4. The small magnitude of the CH4 fluxes and observed control or CH4 fluxes by forest productivity and summer precipitation provide novel insight into the CH4 cycle in this globally important forest ecosystem.

The quantification of carbon fluxes between the terrestrial biosphere and the atmosphere is of scientific importance and also relevant to climate-policy making. Eddy covariance flux towers provide continuous measurements of ecosystem-level exchange of carbon dioxide spanning diurnal, synoptic, seasonal, and interannual time scales. However, these measurements only represent the fluxes at the scale of the tower footprint. Here we used remotely sensed data from the Moderate Resolution Imaging Spectroradiometer (MODIS) to upscale gross primary productivity (GPP) data from eddy covariance flux towers to the continental scale. We first combined GPP and MODIS data for 42 AmeriFlux towers encompassing a wide range of ecosystem and climate types to develop a predictive GPP model using a regression tree approach. The predictive model was trained using observed GPP over the period 2000–2004, and was validated using observed GPP over the period 2005–2006 and leave-one-out cross-validation. Our model predicted GPP fairly well at the site level. We then used the model to estimate GPP for each 1 km × 1 km cell across the U.S. for each 8-day interval over the period from February 2000 to December 2006 using MODIS data. Our GPP estimates provide a spatially and temporally continuous measure of gross primary production for the U.S. that is a highly constrained by eddy covariance flux data. Our study demonstrated that our empirical approach is effective for upscaling eddy flux GPP data to the continental scale and producing continuous GPP estimates across multiple biomes. With these estimates, we then examined the patterns, magnitude, and interannual variability of GPP. We estimated a gross carbon uptake between 6.91 and 7.33 Pg C yr− 1 for the conterminous U.S. Drought, fires, and hurricanes reduced annual GPP at regional scales and could have a significant impact on the U.S. net ecosystem carbon exchange. The sources of the interannual variability of U.S. GPP were dominated by these extreme climate events and disturbances.

Digital repeat photography is becoming widely used for near-surface remote sensing of vegetation. Canopy greenness, which has been used extensively for phenological applications, can be readily quantified from camera images. Important questions remain, however, as to whether the observed changes in canopy greenness are directly related to changes in leaf-level traits, changes in canopy structure, or some combination thereof.

We investigated relationships between canopy greenness and various metrics of canopy structure and function, using five years (2008–2012) of automated digital imagery, ground observations of phenological transitions, leaf area index (LAI) measurements, and eddy covariance estimates of gross ecosystem photosynthesis from the Harvard Forest, a temperate deciduous forest in the northeastern United States. Additionally, we sampled canopy sunlit leaves on a weekly basis throughout the growing season of 2011. We measured physiological and morphological traits including leaf size, mass (wet/dry), nitrogen content, chlorophyll fluorescence, and spectral reflectance and characterized individual leaf color with flatbed scanner imagery.

Our results show that observed spring and autumn phenological transition dates are well captured by information extracted from digital repeat photography. However, spring development of both LAI and the measured physiological and morphological traits are shown to lag behind spring increases in canopy greenness, which rises very quickly to its maximum value before leaves are even half their final size. Based on the hypothesis that changes in canopy greenness represent the aggregate effect of changes in both leaf-level properties (specifically, leaf color) and changes in canopy structure (specifically, LAI), we developed a two end-member mixing model. With just a single free parameter, the model was able to reproduce the observed seasonal trajectory of canopy greenness. This analysis shows that canopy greenness is relatively insensitive to changes in LAI at high LAI levels, which we further demonstrate by assessing the impact of an ice storm on both LAI and canopy greenness.

Our study provides new insights into the mechanisms driving seasonal changes in canopy greenness retrieved from digital camera imagery. The nonlinear relationship between canopy greenness and canopy LAI has important implications both for phenological research applications and for assessing responses of vegetation to disturbances.

Sensor-based monitoring of vegetation phenology is being widely used to quantify phenological responses to climate variability and change. Digital repeat photography, in particular, can characterize the seasonality of canopy greenness. However, these data cannot be directly compared to satellite vegetation indices (e.g. NDVI, the normalized difference vegetation index) that require information about vegetation properties at near-infrared (NIR) wavelengths. Here, we develop a new method, using an inexpensive, NIR-enabled camera originally designed for security monitoring, to calculate a “camera NDVI” from sequential visible and visible + NIR photographs. We use a lab experiment for proof-of-concept, and then test the method using a year of data from an ongoing field campaign in a mixed temperate forest. Our analysis shows that the seasonal cycle of camera NDVI is almost identical to that of NDVI measured using narrow-band radiometric instruments, or as observed from space by the MODIS platform. This camera NDVI thus provides different information about the state of the canopy than can be obtained using only visible-wavelength imagery. In addition to phenological monitoring, our method should be useful for a variety of applications, including continuous monitoring of plant stress and quantifying vegetation responses to manipulative treatments in large field experiments.

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 needle leaf 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 needle leaf 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.

Phenology, by controlling the seasonal activity of vegetation on the land surface, plays a fundamental role in regulating photosynthesis and other ecosystem processes, as well as competitive interactions and feedbacks to the climate system. We conducted an analysis to evaluate the representation of phenology, and the associated seasonality of ecosystem-scale CO2 exchange, in 14 models participating in the North American Carbon Program Site Synthesis. Model predictions were evaluated using long-term measurements (emphasizing the period 2000-2006) from 10 forested sites within the AmeriFlux and Fluxnet-Canada networks. In deciduous forests, almost all models consistently predicted that the growing season started earlier, and ended later, than was actually observed; biases of 2 weeks or more were typical. For these sites, most models were also unable to explain more than a small fraction of the observed interannual variability in phenological transition dates. Finally, for deciduous forests, misrepresentation of the seasonal cycle resulted in over-prediction of gross ecosystem photosynthesis by +160 ± 145 g C m-2 y-1 during the spring transition period, and +75 ± 130 g C m-2 y-1 during the autumn transition period (13% and 8% annual productivity, respectively) compensating for the tendency of most models to under-predict the magnitude of peak summertime photosynthetic rates. Models did a better job of predicting the seasonality of CO2 exchange for evergreen forests. These results highlight the need for improved understanding of the environmental controls on vegetation phenology, and incorporation of this knowledge into better phenological models. Existing models are unlikely to predict future responses of phenology to climate change accurately, and therefore will misrepresent the seasonality and interannual variability of key biosphere-atmosphere feedbacks and interactions in coupled global climate models.