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 yr−1 during the spring transition period and +75 ± 130 g C m−2 yr−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.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Accurately simulating gross primary productivity (GPP) in terrestrial ecosystem models is critical because errors in simulated GPP propagate through the model to introduce additional errors in ...simulated biomass and other fluxes. We evaluated simulated, daily average GPP from 26 models against estimated GPP at 39 eddy covariance flux tower sites across the United States and Canada. None of the models in this study match estimated GPP within observed uncertainty. On average, models overestimate GPP in winter, spring, and fall, and underestimate GPP in summer. Models overpredicted GPP under dry conditions and for temperatures below 0°C. Improvements in simulated soil moisture and ecosystem response to drought or humidity stress will improve simulated GPP under dry conditions. Adding a low‐temperature response to shut down GPP for temperatures below 0°C will reduce the positive bias in winter, spring, and fall and improve simulated phenology. The negative bias in summer and poor overall performance resulted from mismatches between simulated and observed light use efficiency (LUE). Improving simulated GPP requires better leaf‐to‐canopy scaling and better values of model parameters that control the maximum potential GPP, such asεmax (LUE), Vcmax (unstressed Rubisco catalytic capacity) or Jmax (the maximum electron transport rate).
Key Points
Gross primary productivity (GPP) from 26 models tested at 39 flux tower sites
Simulated light use efficiency controls model performance
Models overpredict GPP under dry conditions
Terrestrial biosphere models can help identify physical processes that control carbon dynamics, including land-atmosphere CO2 fluxes, and have great potential to predict the terrestrial ecosystem ...response to changing climate. The skill of models that provide continental scale carbon flux estimates, however, remains largely untested. This paper evaluates the performance of continental-scale flux estimates from 17 models against observations from 36 North American flux towers. Fluxes extracted from regional model simulations are compared with co-located flux tower observations at monthly and annual time increments. Site-level model simulations are used to help interpret sources of the mismatch between the regional simulations and site-based observations. On average the regional model runs overestimate the annual gross primary productivity (5%) and total respiration (15%), and significantly underestimate the annual net carbon uptake (64%) during the time period 2000-2005. Comparison with site-level simulations implicate choices specific to regional model simulations as contributors to the gross flux biases, but not the net carbon uptake bias. The models perform the best at simulating carbon exchange at deciduous broadleaf sites; likely because a number of models use prescribed phenology to simulate seasonal fluxes. The models do not perform as well for crop, grass and evergreen sites. The regional models match the observations most closely in terms of seasonal correlation and seasonal magnitude of variation, but have very little skill at inter-annual correlation and minimal skill at inter-annual magnitude of variability. The comparison of site versus regional level model runs demonstrate that 1) the inter-annual correlation is higher for site-level model runs but the skill remains low, and 2) the underestimation of year-to-year variability for all fluxes is an inherent weakness of the models. The best performing regional models that do not use flux tower calibration are CLM-CN, CASA-GFEDv2 and SIB3. Two flux tower calibrated, empirical models, EC-MOD and MOD17+, perform as well as the best process-based models. This suggests that 1) empirical, calibrated models can perform as well as complex, process-based models, and 2) combining process-based model structure with relevant constraining data could significantly improve model performance.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, INZLJ, KILJ, NLZOH, NMLJ, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZRSKP
Most carbon stocks and fluxes in the western United States are found in mountainous terrain, where observations and modeling are difficult. Terrestrial biosphere models generally underestimate ...above-ground biomass (AGB) over this region. Here, we identify methods to reduce this underestimation by focusing upon 1) biases in meteorological datasets, 2) model representation of water stress, and 3) spatial resolution. We adopted the widely-used Community Land Model version 4.5 (CLM 4.5) with six different meteorological datasets and found a 6-fold variation in simulated AGB across Utah/Colorado. Simulations underestimated AGB because of warm and dry biases within the meteorological datasets that reduced water availability and restricted plant growth. To eliminate the AGB underestimation we adopted a meteorological dataset designed for complex terrain (gridMET), combined with a representation of plant hydraulic stress (CLM 5.0). Conversely, changes in spatial resolution (meteorological variables and land surface description) had negligible impact on simulated AGB.
•Default CLM 4.5 underestimated above-ground biomass (AGB) across Central Rockies.•Modeled AGB varied 6-fold with bias-corrected meteorological datasets.•GridMET meteorology and representing plant hydraulics removed the low bias in AGB.•Simulations implemented at fine and coarse spatial resolution provided the same AGB.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Droughts in the western United States are expected to intensify with climate change. Thus, an adequate representation of ecosystem response to water stress in land models is critical for predicting ...carbon dynamics. The goal of this study was to evaluate the performance of the Community Land Model (CLM) version 4.5 against observations at an old-growth coniferous forest site in the Pacific Northwest region of the United States (Wind River AmeriFlux site), characterized by a Mediterranean climate that subjects trees to water stress each summer. CLM was driven by site-observed meteorology and calibrated primarily using parameter values observed at the site or at similar stands in the region. Key model adjustments included parameters controlling specific leaf area and stomatal conductance. Default values of these parameters led to significant underestimation of gross primary production, overestimation of evapotranspiration, and consequently overestimation of photosynthetic 13C discrimination, reflected in reduced 13C : 12C ratios of carbon fluxes and pools. Adjustments in soil hydraulic parameters within CLM were also critical, preventing significant underestimation of soil water content and unrealistic soil moisture stress during summer. After calibration, CLM was able to simulate energy and carbon fluxes, leaf area index, biomass stocks, and carbon isotope ratios of carbon fluxes and pools in reasonable agreement with site observations. Overall, the calibrated CLM was able to simulate the observed response of canopy conductance to atmospheric vapor pressure deficit (VPD) and soil water content, reasonably capturing the impact of water stress on ecosystem functioning. Both simulations and observations indicate that stomatal response from water stress at Wind River was primarily driven by VPD and not soil moisture. The calibration of the Ball–Berry stomatal conductance slope (mbb) at Wind River aligned with findings from recent CLM experiments at sites characterized by the same plant functional type (needleleaf evergreen temperate forest), despite significant differences in stand composition and age and climatology, suggesting that CLM could benefit from a revised mbb value of 6, rather than the default value of 9, for this plant functional type. Conversely, Wind River required a unique calibration of the hydrology submodel to simulate soil moisture, suggesting that the default hydrology has a more limited applicability. This study demonstrates that carbon isotope data can be used to constrain stomatal conductance and intrinsic water use efficiency in CLM, as an alternative to eddy covariance flux measurements. It also demonstrates that carbon isotopes can expose structural weaknesses in the model and provide a key constraint that may guide future model development.
Most carbon stocks and fluxes in the western United States are found in mountainous terrain, where observations and modeling are difficult. Terrestrial biosphere models generally underestimate ...above-ground biomass (AGB) over this region. In this report, we identify methods to reduce this underestimation by focusing upon 1) biases in meteorological datasets, 2) model representation of water stress, and 3) spatial resolution. We adopted the widely-used Community Land Model version 4.5 (CLM 4.5) with six different meteorological datasets and found a 6-fold variation in simulated AGB across Utah/Colorado. Simulations underestimated AGB because of warm and dry biases within the meteorological datasets that reduced water availability and restricted plant growth. To eliminate the AGB underestimation we adopted a meteorological dataset designed for complex terrain (gridMET), combined with a representation of plant hydraulic stress (CLM 5.0). Conversely, changes in spatial resolution (meteorological variables and land surface description) had negligible impact on simulated AGB.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Abstract
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
CO
2
exchange, in 14 models participating in the
N
orth
A
merican
C
arbon
P
rogram
S
ite
S
ynthesis. Model predictions were evaluated using long‐term measurements (emphasizing the period 2000–2006) from 10 forested sites within the
A
meri
F
lux and
F
luxnet‐
C
anada 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
yr
−1
during the spring transition period and +75 ± 130 g C m
−2
yr
−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
CO
2
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.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
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.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Recent successes in passive remote sensing of far-red solar-induced chlorophyll fluorescence (SIF) have spurred the development and integration of
canopy-level fluorescence models in global ...terrestrial biosphere models (TBMs) for climate and carbon cycle research. The interaction of fluorescence
with photochemistry at the leaf and canopy scales provides opportunities to diagnose and constrain model simulations of photosynthesis and related
processes, through direct comparison to and assimilation of tower, airborne, and satellite data. TBMs describe key processes related to the absorption of
sunlight, leaf-level fluorescence emission, scattering, and reabsorption throughout the canopy. Here, we analyze simulations from an ensemble of
process-based TBM–SIF models (SiB3 – Simple Biosphere Model, SiB4, CLM4.5 – Community Land Model, CLM5.0, BETHY – Biosphere Energy Transfer Hydrology, ORCHIDEE – Organizing Carbon and Hydrology In Dynamic Ecosystems, and BEPS – Boreal Ecosystems Productivity Simulator) and the SCOPE (Soil Canopy Observation Photosynthesis Energy) canopy radiation and vegetation model at a subalpine
evergreen needleleaf forest near Niwot Ridge, Colorado. These models are forced with local meteorology and analyzed against tower-based continuous
far-red SIF and gross-primary-productivity-partitioned (GPP) eddy covariance data at diurnal and synoptic scales during the growing season
(July–August 2017). Our primary objective is to summarize the site-level state of the art in TBM–SIF modeling over a relatively short time period
(summer) when light, canopy structure, and pigments are similar, setting the stage for regional- to global-scale analyses. We find that these models
are generally well constrained in simulating photosynthetic yield but show strongly divergent patterns in the simulation of absorbed photosynthetic
active radiation (PAR), absolute GPP and fluorescence, quantum yields, and light response at the leaf and canopy scales. This study highlights the need for
mechanistic modeling of nonphotochemical quenching in stressed and unstressed environments and improved the representation of light absorption (APAR),
distribution of light across sunlit and shaded leaves, and radiative transfer from the leaf to the canopy scale.
The Western United States is dominated by natural lands that play a critical role for carbon balance, water quality, and timber reserves. This region is also particularly vulnerable to forest ...mortality from drought, insect attack, and wildfires, thus requiring constant monitoring to assess ecosystem health. Carbon monitoring techniques are challenged by the complex mountainous terrain, thus there is an opportunity for data assimilation systems that combine land surface models and satellite‐derived observations to provide improved carbon monitoring. Here, we use the Data Assimilation Research Testbed to adjust the Community Land Model (CLM5.0) with remotely sensed observations of leaf area and above‐ground biomass. The adjusted simulation significantly reduced the above‐ground biomass and leaf area, leading to a reduction in both photosynthesis and respiration fluxes. The reduction in the carbon fluxes mostly offset, thus both the adjusted and free simulation projected a weak carbon sink to the land. This result differed from a separate observation‐constrained model (FLUXCOM) that projected strong carbon uptake to the land. Simulation diagnostics suggested water limitation had an important influence upon the magnitude and spatial pattern of carbon uptake through photosynthesis. We recommend that additional observations important for water cycling (e.g., snow water equivalent, land surface temperature) be included to improve the veracity of the spatial pattern in carbon uptake. Furthermore, the assimilation system should be enhanced to maximize the number of the simulated state variables that are adjusted, especially those related to the recommended observed quantities including water cycling and soil carbon.
Plain Language Summary
The Western United States is dominated by natural lands that play a critical role for carbon balance (e.g., trees, soils), water quality, and timber reserves. This region is also particularly vulnerable to tree death from drought, insect attack, and wildfires, thus requiring constant monitoring to assess its health. Traditional carbon monitoring techniques are usually not possible within mountainous terrain, thus we used satellite observations of leaf area and forest biomass to improve modeled simulations of the Western United States. When we accounted for observations of trees our modeled estimates showed reduced amounts of biomass and relatively small amounts of atmospheric CO2 transfer from the atmosphere to the land (the land absorbs carbon from the atmosphere through photosynthesis). Our best estimate of carbon absorbed by the land was much less than other modeled estimates. This suggests our method better accounted for the current conditions of the trees including death from fire, insect attack, and drought. Our modeled estimate of biomass and carbon balance across the Western United States can be improved further by considering more observations of the land surface related to soil moisture and soil carbon.
Key Points
Assimilating observations of biomass and leaf area reduces simulated biomass and projects a weak land carbon sink across the Western United States
Our estimate of carbon exchange contrasts with an independent FLUXCOM estimate that shows a significant carbon sink in the Western United States
Water cycle observations should be used to complement biomass observations to improve the spatial pattern of modeled carbon fluxes
Full text
Available for:
DOBA, FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK