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
A large area of boreal jack pine (Pinus banksiana Lamb.) forest in Canada is recovering from clear-cut harvesting, and the carbon (C) balance of these regenerating forests remains uncertain. Net ...ecosystem CO₂ exchange was measured using the eddy-covariance technique at four jack pine sites representing different stages of stand development: three postharvest sites (HJP02, HJP94, and HJP75) and one preharvest site (OJP). The four sites, located in the southern Canadian boreal forest, Saskatchewan, Canada, are typical of low productivity jack pine stands and were 2, 10, 29, and 90 years old in 2004, respectively. Mean annual net ecosystem production (NEP) for 2004 and 2005 was -137±11, 19±16, 73±28, and 22±30 g C m⁻² yr⁻¹ at HJP02, HJP94, HJP75 and OJP, respectively, showing the postharvest jack pine stands to be moderate C sources immediately after harvesting, weak sinks at 10 years, moderate C sinks at 30 years, then weak C sinks at 90 years. Mean annual gross ecosystem photosynthesis (GEP) for the 2 years was 96±10, 347±20, 576±34, and 583±35 g C m⁻² yr⁻¹ at HJP02, HJP94, HJP75, and OJP, respectively. The ratio of annual ecosystem respiration (R) to annual GEP was 2.51±0.15, 0.95±0.04, 0.87±0.03, and 0.96±0.03. Seasonally, NEP peaked in May or June at all four sites but GEP and R were highest in July. R at a reference soil temperature of 10 °C, ecosystem quantum yield and photosynthetic capacity were lowest for the 2-year-old stand. R was most sensitive to soil temperature for the 90-year-old stand. The primary source of variability in NEP over the course of succession of the jack pine ecosystem following harvesting was stand age due to the changes in leaf area index. Intersite variability in GEP and R was an order of magnitude greater than interannual variability at OJP. For both young and old stands, GEP had greater interannual variability than R and played a more important role than R in interannual variation in NEP. Based on year-round flux measurements from 2000 to 2005, the 10-year stand had larger interannual variability in GEP and R than the 90-year stand. Interannual variability in NEP was driven primarily by early-growing-season temperature and growing-season length. Photosynthesis played a dominant role in the rapid rise in NEP early in stand development. Late in stand development, however, the subtle decrease in NEP resulted primarily from increasing respiration.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Accurate estimation of latent heat flux (LE) based on remote sensing data is critical in characterizing terrestrial ecosystems and modeling land surface processes. Many LE products were released ...during the past few decades, but their quality might not meet the requirements in terms of data consistency and estimation accuracy. Merging multiple algorithms could be an effective way to improve the quality of existing LE products. In this paper, we present a data integration method based on modified empirical orthogonal function (EOF) analysis to integrate the Moderate Resolution Imaging Spectroradiometer (MODIS) LE product (MOD16) and the Priestley-Taylor LE algorithm of Jet Propulsion Laboratory (PT-JPL) estimate. Twenty-two eddy covariance (EC) sites with LE observation were chosen to evaluate our algorithm, showing that the proposed EOF fusion method was capable of integrating the two satellite data sets with improved consistency and reduced uncertainties. Further efforts were needed to evaluate and improve the proposed algorithm at larger spatial scales and time periods, and over different land cover types.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
More accurate estimation of the carbon dioxide flux depends on the improved scientific understanding of the terrestrial carbon cycle. Remote-sensing-based approaches to continental-scale estimation ...of net ecosystem exchange (NEE) have been developed but coarse spatial resolution is a source of errors. Here we demonstrate a satellite-based method of estimating NEE using Landsat TM/ETM+data and an upscaling framework. The upscaling framework contains flux-footprint climatology modeling, modified regression tree (MRT) analysis and image fusion. By scaling NEE measured at flux towers to landscape and regional scales, this satellite-based method can improve NEE estimation at high spatial-temporal resolution at the landscape scale relative to methods based on MODIS data with coarser spatial–temporal resolution. This method was applied to sixteen flux sites from the Canadian Carbon Program and AmeriFlux networks located in North America, covering forest, grass, and cropland biomes. Compared to a similar method using MODIS data, our estimation is more effective for diagnosing landscape NEE with the same temporal resolution and higher spatial resolution (30m versus 1km) (r2=0.7548 vs. 0.5868, RMSE=1.3979 vs. 1.7497gCm−2day−1, average error=0.8950 vs. 1.0178gCm−2day−1, relative error=0.47 vs. 0.54 for fused Landsat and MODIS imagery, respectively). We also compared the regional NEE estimations using Carbon Tracker, our method and eddy-covariance observations. This study demonstrates that the data-driven satellite-based NEE diagnosed model can be used to upscale eddy-flux observations to landscape scales with high spatial–temporal resolutions.
•NEE at landscape scale with high spatial-temporal resolution.•Generated NEE prediction model only depended on remote-sensing data.•Upscaling framework used footprint climatology modeling, MRT analysis & image fusion.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
▶ Climate models predict rising temperatures and more frequent and prolonged droughts, particularly in the northern hemisphere and in the Canadian Prairies. However, few studies have examined the ...interannual variation in evapotranspiration (
E) of northern terrestrial ecosystems in relation to drought. This study analyses multi-year (1998-2006), eddy-covariance measurements to quantify the seasonal and interannual variability in
E from grassland and mature aspen, black spruce, and jack pine forest ecosystems in the dry interior plains of western Canada. It also investigates the response of
E to the historic 2001-2003 drought in this region. Leaf area index (LAI) was a primary factor controlling the difference in
E among different ecosystems in the same ecozone. Annual
E was higher and more variable for the aspen forest (405 ± 84
mm, mean ± s.d.) and grassland (395 ± 90 mm) than the black spruce (374 ± 34 mm) and jack pine (300 ± 20 mm) forests. Interannual variation of
E was controlled by early spring soil temperature at all four sites, with warm springs enhancing annual
E. Only the aspen forest and grassland showed a significant suppression of
E by summer drought, related to reduced surface conductance at the aspen site, and to reduced surface conductance and early leaf senescence at the grassland. We conclude that the potential for drought impacts on annual
E of northern ecosystems is greatest for grasslands, moderate for deciduous broadleaf aspen forests, and smallest for coniferous black spruce and jack pine forests. The forests of the Boreal Plains, adjacent to the prairie region, may ameliorate the onset of drought through the recycling of moisture to the atmosphere, whereas the prairie grasslands have only limited capacity to counteract drought through moisture recycling.
Climate models predict rising temperatures and more frequent and prolonged droughts, particularly in the northern hemisphere and in the Canadian Prairies. However, few studies have examined the interannual variation in evapotranspiration (E) of northern terrestrial ecosystems in relation to drought. This study analyses multi-year (1998-2006), eddy-covariance measurements to quantify the seasonal and interannual variability in E from grassland and mature aspen, black spruce, and jack pine forest ecosystems in the dry interior plains of western Canada. It also investigates the response of E to the historic 2001-2003 drought in this region. Leaf area index (LAI) was a primary factor controlling the difference in E among different ecosystems in the same ecozone. Annual E was higher and more variable for the aspen forest (405
±
84
mm, mean
±
s.d.) and grassland (395
±
90
mm) than the black spruce (374
±
34
mm) and jack pine (300
±
20
mm) forests. Interannual variation of E was controlled by early spring soil temperature at all four sites, with warm springs enhancing annual E. Only the aspen forest and grassland showed a significant suppression of E by summer drought, related to reduced surface conductance at the aspen site, and to reduced surface conductance and early leaf senescence at the grassland. We conclude that the potential for drought impacts on annual E of northern ecosystems is greatest for grasslands, moderate for deciduous broadleaf aspen forests, and smallest for coniferous black spruce and jack pine forests. The forests of the Boreal Plains, adjacent to the prairie region, may ameliorate the onset of drought through the recycling of moisture to the atmosphere, whereas the prairie grasslands have only limited capacity to counteract drought through moisture recycling.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
► We analyzed the effects of climatic factors and vegetation type on evapotranspiration and water use efficiency for mature forest, grassland and peatland sites across a continental-scale transect in ...Canada. ► Despite highly variable precipitation among sites (250−1450mm), annual evapotranspiration was limited to 400−500mm. ► The seasonal course of net radiation does not account for the seasonal variation in evapotranspiration of all ecosystems. ► Analysis of daytime dry-foliage Priestley–Taylor α and canopy conductance indicated mainly stomatal limitation to transpiration at most sites. ► Water use efficiency was relatively constant with values in the range of 2.6–3.6g C kg−1 H2O, while climate, differing LAI and biomass abundance associated with different plant functional types were responsible for values outside this range.
The effects of climatic factors and vegetation type on evapotranspiration (E) and water use efficiency (WUE) were analyzed using tower-based eddy-covariance (EC) data for nine mature forest sites, two peatland sites and one grassland site across an east–west continental-scale transect in Canada during the period 2003–2006. The seasonal pattern of E was closely linked to growing-season length and rainfall distribution. Although annual precipitation (P) during the observation period was highly variable among sites (250−1450mm), minimum annual E was not less than 200mm and was limited to 400−500mm where annual P exceeded 700mm. Site-specific interannual variability in E could be explained by either changes in total P or variations in solar irradiance. A highly positive linear correlation was found between monthly mean values of E and net radiation (Rn) at the grassland site (AB-GRL), the two peatland sites (AB-WPL and ON-EPL), and only one of the forest sites (coastal Douglas-fir, BC-DF49) whereas a hysteretic relationship at the other forest sites indicated that E lagged behind the typical seasonal progression of Rn. Results of a cross-correlation analysis between daily (24-h) E and Rn revealed that site-specific lag times were between 10 and 40 days depending on the lag of vapour pressure deficit (D) behind Rn and the decoupling coefficient, Ω. There was significant seasonal variation in daytime mean dry-foliage Priestley–Taylor α with maxima occurring in the growing season at all sites except BC-DF49 where it was relatively constant (∼0.55) throughout all years. Annual means of daytime dry-foliage α mostly ranging between 0.5 and 0.7 implied stomatal limitation to transpiration. Increasing D significantly decreased canopy conductance (gc) at the forest sites but had little effect at the peatland and grassland sites, while variation in soil water content caused only minor changes in gc. At all sites, a strong linear correlation between monthly mean values of gross primary production (GPP) and E resulted in water use efficiency being relatively constant. While at most sites, WUE was in the range of 2.6–3.6gCkg−1 H2O, the BC-DF49 site had the highest WUE of the twelve sites with values near 6.0gCkg−1 H2O. Of the two peatland sites, AB-WPL, a western treed fen, had a significantly higher WUE (∼3.0gCkg−1 H2O) than ON-EPL, an eastern ombrotrophic bog (∼1.8gCkg−1 H2O), which was related to peatland productivity and plant functional type.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK, ZRSKP
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⁻² yr⁻¹; while temperate forests were either a medium or strong C sink, ranging from 150 to 500 g C m⁻² yr⁻¹, 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.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
The leaf area index (LAI) and the clumping index (CI) provide valuable insight into the spatial patterns of forest canopies, the canopy light regime and forest productivity. This study examines the ...spatial patterns of LAI and CI in a boreal mixed-wood forest, using extensive field measurements and remote sensing analysis. The objectives of this study are to: (1) examine the utility of airborne lidar (light detection and ranging) and hyperspectral data to model LAI and clumping indices; (2) compare these results to those found from commonly used Landsat vegetation indices (i.e. the normalized difference vegetation index (NDVI) and the simple ratio (SR)); (3) determine whether the fusion of lidar data with Landsat and/or hyperspectral data will improve the ability to model clumping and LAI; and (4) assess the relationships between clumping, LAI and canopy biochemistry.Regression models to predict CI were much stronger than those for LAI at the site. Lidar was the single best predictor of CI (r² > 0.8). Landsat NDVI and SR also had a moderately strong predictive performance for CI (r² > 0.68 with simple linear and non-linear regression forms), suggesting that canopy clumping can be predicted operationally from satellite platforms, at least in boreal mixed-wood environments. Foliar biochemistry, specifically canopy chlorophyll, carotenoids, magnesium, phosphorus and nitrogen, was strongly related to the clumping index. Combined, these results suggest that Landsat models of clumping could provide insight into the spatial distribution of foliar biochemistry, and thereby photosynthetic capacity, for boreal mixed-wood canopies. LAI models were weak (r² < 0.4) unless separate models were used for deciduous and coniferous plots. Coniferous LAI was easier to model than deciduous LAI (r² > 0.8 for several indices). Deciduous models of LAI were weaker for all remote sensing indices (r² < 0.67). There was a strong, linear relationship between foliar biochemistry and LAI for the deciduous plots. Overall, our results suggest that broadband satellite indices have strong predictive performance for clumping, but that airborne hyperspectral or lidar data are required to develop strong models of LAI at this boreal mixed-wood site.
Full text
Available for:
BFBNIB, GIS, IJS, KISLJ, NUK, PNG, UL, UM, UPUK
Net ecosystem productivity (NEP) during August 2003 was measured by using eddy covariance above 17 forest and 3 peatland sites along an east-west continental-scale transect in Canada. Measured sites ...included recently disturbed stands, young forest stands, intermediate-aged conifer stands, mature deciduous stands, mature conifer stands, fens, and an open shrub bog. Diurnal courses of NEP showed strong coherence within the different ecosystem categories. Recently disturbed sites showed the weakest diurnal cycle; and intermediate-aged conifers, the strongest. The western treed fen had a more pronounced diurnal pattern than the eastern shrub bog or the Saskatchewan patterned fen. All but three sites were clearly afternoon C sinks. Ecosystem respiration was highest for the young fire sites. The intermediate-aged conifer sites had the highest maximum NEP (NEPmax) and gross ecosystem productivity (GEPmax), attaining rates that would be consistent with the presence of a strong terrestrial C sink in regions where these types of forest are common. These results support the idea that large-scale C cycle modeling activities would benefit from information on the age-class distribution and disturbance types within larger grid cells. Light use efficiency followed a pattern similar to that of NEPmax and GEPmax. Four of the five recently disturbed sites and all three of the peatland sites had low water use efficiencies.
Full text
Available for:
BF, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Components of the surface energy balance of a mature boreal jack pine forest and a jack pine clearcut were analysed to determine the causes of the imbalance that is commonly observed in ...micrometeorological measurements. At the clearcut site (HJP02), a significant portion of the imbalance was caused by: (i) the overestimation of net radiation (R n ) due to the inclusion of the tower in the field of view of the downward facing radiometers, and (ii) the underestimation of the latent heat flux (λE) due to the damping of high frequency fluctuations in the water vapour mixing ratio by the sample tube of the closed-path infrared gas analyzer. Loss of low-frequency covariance induced by insufficient averaging time as well as systematic advection of fluxes away from the eddy-covariance (EC) tower were discounted as significant issues. Spatial and temporal distributions of the total surface-layer heat flux (T), i.e. the sum of sensible heat flux (H) and λE, were well behaved and differences between the relative magnitudes of the turbulent fluxes for several investigated energy balance closure (C) classes were observed. Therefore, it can be assumed that micrometeorological processes that affected all turbulent fluxes similarly did not cause the variation in C. Turbulent fluxes measured at the clearcut site should not be forced to close the energy balance. However, at the mature forest site (OJP), loss of low-frequency covariance contributed significantly to the systematic imbalance when a 30-min averaging time was used, but the application of averaging times that were long enough to capture all of the low-frequency covariance was inadequate to resolve all of the high-frequency covariance. Although we found qualitative similarity between T and the net ecosystem exchange (NEE) of carbon dioxide (CO₂), forcing T to closure while retaining the Bowen ratio and applying the same factor to CO₂ fluxes (F C ) cannot be generally recommended since it remains uncertain to what extent long wavelength contributions affect the relationship between T, F C and C.
Full text
Available for:
DOBA, EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, IZUM, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UILJ, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ