Modeling Analysis of Primary Controls on Net Ecosystem Productivity of Seven Boreal and Temperate Coniferous Forests Across a Continental Transect

Process-based models are effective tools to synthesize and/or extrapolate measured carbon (C) exchanges from individual sites to large scales. In this study, we used a C- and nitrogen (N)-cycle coupled ecosystem model named CN-CLASS (Carbon Nitrogen-Canadian Land Surface Scheme) to study the role of primary climatic controls and site-specific C stocks on the net ecosystem productivity (NEP) of seven intermediate-aged to mature coniferous forest sites across an east–west continental transect in Canada. The model was parameterized using a common set of parameters, except for two used in empirical canopy conductance–assimilation, and leaf area–sapwood relationships, and then validated using observed eddy covariance flux data. Leaf Rubisco-N dynamics that are associated with soil–plant N cycling, and depend on canopy temperature, enabled the model to simulate site-specific gross ecosystem productivity (GEP) reasonably well for all seven sites. Overall GEP simulations had relatively smaller differences compared with observations vs. ecosystem respiration (RE), which was the sum of many plant and soil components with larger variability and/or uncertainty associated with them. Both observed and simulated data showed that, on an annual basis, boreal forest sites were either carbon-neutral or a weak C sink, ranging from 30 to 180 g C m−2 yr−1; while temperate forests were either a medium or strong C sink, ranging from 150 to 500 g C m−2 yr−1, depending on forest age and climatic regime. Model sensitivity tests illustrated that air temperature, among climate variables, and aboveground biomass, among major C stocks, were dominant factors impacting annual NEP. Vegetation biomass effects on annual GEP, RE and NEP showed similar patterns of variability at four boreal and three temperate forests. Air temperature showed different impacts on GEP and RE, and the response varied considerably from site to site. Higher solar radiation enhanced GEP, while precipitation differences had a minor effect. Magnitude of forest litter content and soil organic matter (SOM) affected RE. SOM also affected GEP, but only at low levels of SOM, because of low N mineralization that limited soil nutrient (N) availability. The results of this study will help to evaluate the impact of future climatic changes and/or forest C stock variations on C uptake and loss in forest ecosystems growing in diverse environments.Earth and Planetary Science

Terrestrial Biosphere Model Performance for Inter-Annual Variability of Land-Atmosphere CO2 Exchange

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.Organismic and Evolutionary Biolog

Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis

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.Engineering and Applied SciencesOrganismic and Evolutionary Biolog

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

The FLUXNET2015 dataset provides ecosystem-scale data on CO2, 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.Peer reviewe

A model-data comparison of gross primary productivity: Results from the North American Carbon Program site synthesis

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 ε[subscript max] (LUE), V[subscript cmax] (unstressed Rubisco catalytic capacity) or J[subscript max] (the maximum electron transport rate)

Climate control of terrestrial carbon exchange across biomes and continents

Peer reviewe

Prediction of Macronutrients at the Canopy Level Using Spaceborne Imaging Spectroscopy and LiDAR Data in a Mixedwood Boreal Forest

Information on foliar macronutrients is required in order to understand plant physiological and ecosystem processes such as photosynthesis, nutrient cycling, respiration and cell wall formation. The ability to measure, model and map foliar macronutrients (nitrogen (N), phosphorus (P), potassium (K), calcium (Ca) and magnesium (Mg)) at the forest canopy level provides information on the spatial patterns of ecosystem processes (e.g., carbon exchange) and provides insight on forest condition and stress. Imaging spectroscopy (IS) has been used particularly for modeling N, using airborne and satellite imagery mostly in temperate and tropical forests. However, there has been very little research conducted at these scales to model P, K, Ca, and Mg and few studies have focused on boreal forests. We report results of a study of macronutrient modeling using spaceborne IS and airborne light detection and ranging (LiDAR) data for a mixedwood boreal forest canopy in northern Ontario, Canada. Models incorporating Hyperion data explained approximately 90% of the variation in canopy concentrations of N, P, and Mg; whereas the inclusion of LiDAR data significantly improved the prediction of canopy concentration of Ca (R2 = 0.80). The combined used of IS and LiDAR data significantly improved the prediction accuracy of canopy Ca and K concentration but decreased the prediction accuracy of canopy P concentration. The results indicate that the variability of macronutrient concentration due to interspecific and functional type differences at the site provides the basis for the relationship observed between the remote sensing measurements (i.e., IS and LiDAR) and macronutrient concentration. Crown closure and canopy height are the structural metrics that establish the connection between macronutrient concentration and IS and LiDAR data, respectively. The spatial distribution of macronutrient concentration at the canopy scale mimics functional type distribution at the site. The ability to predict canopy N, P, K, Ca and Mg in this study using only IS, only LiDAR or their combination demonstrates the excellent potential for mapping these macronutrients at canopy scales across larger geographic areas into the next decade with the launch of new IS satellite missions and by using spaceborne LiDAR data