We upscaled FLUXNET observations of carbon dioxide, water, and energy fluxes to the global scale using the machine learning technique, model tree ensembles (MTE). We trained MTE to predict site‐level ...gross primary productivity (GPP), terrestrial ecosystem respiration (TER), net ecosystem exchange (NEE), latent energy (LE), and sensible heat (H) based on remote sensing indices, climate and meteorological data, and information on land use. We applied the trained MTEs to generate global flux fields at a 0.5° × 0.5° spatial resolution and a monthly temporal resolution from 1982 to 2008. Cross‐validation analyses revealed good performance of MTE in predicting among‐site flux variability with modeling efficiencies (MEf) between 0.64 and 0.84, except for NEE (MEf = 0.32). Performance was also good for predicting seasonal patterns (MEf between 0.84 and 0.89, except for NEE (0.64)). By comparison, predictions of monthly anomalies were not as strong (MEf between 0.29 and 0.52). Improved accounting of disturbance and lagged environmental effects, along with improved characterization of errors in the training data set, would contribute most to further reducing uncertainties. Our global estimates of LE (158 ± 7 J × 1018 yr−1), H (164 ± 15 J × 1018 yr−1), and GPP (119 ± 6 Pg C yr−1) were similar to independent estimates. Our global TER estimate (96 ± 6 Pg C yr−1) was likely underestimated by 5–10%. Hot spot regions of interannual variability in carbon fluxes occurred in semiarid to semihumid regions and were controlled by moisture supply. Overall, GPP was more important to interannual variability in NEE than TER. Our empirically derived fluxes may be used for calibration and evaluation of land surface process models and for exploratory and diagnostic assessments of the biosphere.
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
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.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, INZLJ, KILJ, NLZOH, NMLJ, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZRSKP
We report estimates of the amount, distribution, and uncertainty of aboveground biomass (AGB) of the different ecoregions and forest land cover classes within the North American boreal forest, ...analyze the factors driving the error estimates, and compare our estimates with other reported values. A three-phase sampling strategy was used (i) to tie ground plot AGB to airborne profiling lidar metrics and (ii) to link the airborne estimates of AGB to ICESat-GLAS lidar measurements such that (iii) GLAS could be used as a regional sampling tool. We estimated the AGB of the North American boreal forest at 21.8 Pg, with relative error of 1.9% based on 256 GLAS orbits (229 086 pulses). The distribution of AGB was 46.6% for western Canada, 43.7% for eastern Canada, and 9.7% for Alaska. With a single exception, relative errors were under 4% for the three regions and for the major cover types and under 10% at the ecoregion level. The uncertainties of the estimates were calculated using a variance estimator that accounted for only sampling error, i.e., the variability among GLAS orbital estimates, and airborne to spaceborne regression error, i.e., the uncertainty of the model coefficients. Work is ongoing to develop robust statistical techniques for integrating other sources of error such as ground to air regression error and allometric error. Small ecoregions with limited east-west extents tended to have fewer GLAS orbits and a greater percent sampling error. AGB densities derived from GLAS agreed closely with the estimates derived from both forest inventories (<17%) and a MODIS-based interpolation technique (<26%) for more southern, well-inventoried ecoregions, whereas differences were much greater for unmanaged northern and (or) mountainous ecoregions.
Full text
Available for:
BF, DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Aboveground dry biomass was estimated for the 1.3 M km
2 forested area south of the treeline in the eastern Canadian province of Québec by combining data from an airborne and spaceborne LiDAR, a ...Landsat ETM+ land cover map, a Shuttle Radar Topographic Mission (SRTM) digital elevation model, ground inventory plots, and vegetation zone maps. Plot-level biomass was calculated using allometric relationships between tree attributes and biomass. A small footprint portable laser profiler then flew over these inventory plots to develop a generic airborne LiDAR-based biomass equation (
R
2
=
0.65,
n
=
207). The same airborne LiDAR system flew along four portions of orbits of the ICESat Geoscience Laser Altimeter System (GLAS). A square-root transformed equation was developed to predict airborne profiling LiDAR estimates of aboveground dry biomass from GLAS waveform parameters combined with an SRTM slope index (
R
2
=
0.59,
n
=
1325).
Using the 104,044 quality-filtered GLAS pulses obtained during autumn 2003 from 97 orbits over the study area, we then predicted aboveground dry biomass for the main vegetation areas of Québec as well as for the entire Province south of the treeline. Including cover type covariances both within and between GLAS orbits increased standard errors of the estimates by two to five times at the vegetation zone level and as much as threefold at the provincial level. Aboveground biomass for the whole study area averaged 39.0
±
2.2 (standard error) Mg ha
−
1
and totalled 4.9
±
0.3 Pg. Biomass distributions were 12.6% northern hardwoods, 12.6% northern mixedwood, 38.4% commercial boreal, 13% non-commercial boreal, 14.2% taiga, and 9.2% treed tundra. Non-commercial forests represented 36% of the estimated aboveground biomass, thus highlighting the importance of remote northern forests to C sequestration. This study has shown that space-based forest inventories of northern forests could be an efficient way of estimating the amount, distribution, and uncertainty of aboveground biomass and carbon stocks at large spatial scales.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
► 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
Assisted population migration (APM) is the intentional movement of populations within a species range to sites where future environmental conditions are projected to be more conducive to growth. APM ...has been proposed as a proactive adaptation strategy to maintain forest productivity and to reduce the vulnerability of forest ecosystems to projected climate change. The validity of such a strategy will depend on the adaptation capacity of populations, which can partially be evaluated by the ecophysiological response of different genetic sources along a climatic gradient. This adaptation capacity results from the compromise between (i) the degree of genetic adaptation of seed sources to their environment of origin and (ii) the phenotypic plasticity of functional trait which can make it possible for transferred seed sources to positively respond to new growing conditions. We examined phenotypic variation in morphophysiological traits of six seed sources of white spruce (Picea glauca Moench Voss) along a regional climatic gradient in Québec, Canada. Seedlings from the seed sources were planted at three forest sites representing a mean annual temperature (MAT) gradient of 2.2°C. During the second growing season, we measured height growth (H2014) and traits related to resources use efficiency and photosynthetic rate (A max). All functional traits showed an adaptive response to the climatic gradient. Traits such as H2014, A max, stomatal conductance (g s ), the ratio of mesophyll to stomatal conductance, water use efficiency, and photosynthetic nitrogen-use efficiency showed significant variation in both physiological plasticity due to the planting site and seed source variation related to local genetic adaptation. However, the amplitude of seed source variation was much less than that related to plantation sites in the area investigated. The six seed sources showed a similar level of physiological plasticity. H2014, A max and g s , but not carboxylation capacity (V cmax), were correlated and decreased with a reduction of the average temperature of the growing season at seed origin. The clinal variation in H2014 and A max appeared to be driven by CO2 conductance. The presence of locally adapted functional traits suggests that the use of APM may have advantages for optimizing seed source productivity in future local climates.
Although mature black spruce forests are a dominant cover type in the boreal forest of North America, it is not clear how their carbon (C) budgets vary across the continent. The installation of an ...eddy covariance flux tower on an Old Black Spruce (OBS) site in eastern Canada (EOBS, Québec) provided a first opportunity to compare and contrast its annual (2004) and seasonal C exchange with two other pre-existing OBS flux sites from different climatic regions located in Saskatchewan Southern OBS (SOBS) and Manitoba Northern OBS (NOBS). Although there was a relatively uniform seasonal pattern of net ecosystem productivity (NEP) among sites, EOBS had a lower total annual NEP than the other two sites. This was primarily because warmer soil under a thicker snowpack at EOBS appeared to increase winter C losses and low light suppressed both NEP and gross ecosystem productivity (GEP) in June. Across sites, greater total annual GEP and ecosystem respiration (R) were associated with greater mean annual air temperatures and an earlier beginning of the growing season. Also, GEP at all three sites showed a stronger relationship with air temperature in spring and early summer compared with later in the growing season, highlighting the importance of springtime conditions to the C budget of these boreal ecosystems. The three sites had different parameter estimates describing the responses of R and GEP at the half hour time scale to near surface temperature and light, respectively. On the other hand, the responses of both R and GEP to temperature at the monthly scale did not differ among sites. These results suggest that a general parameterization could be sufficient at coarse time resolutions to model the response of C exchange to environmental factors of mature black spruce forests from different climatic regions.
Full text
Available for:
BFBNIB, FZAB, GIS, IJS, KILJ, NLZOH, NUK, OILJ, SBCE, SBMB, UL, UM, UPUK
Knowledge of thermal acclimation of physiological processes of boreal tree species is necessary to determine their ability to adapt to predicted global warming and reduce the uncertainty around the ...anticipated feedbacks of forest ecosystems and global carbon cycle to climate change. The objective of this work was to examine the extent of thermal acclimation of net photosynthesis (An) and dark respiration (Rd) of two distant white spruce (Picea glauca) seed sources (from south and north of the commerial forest zone in Québec) in response to latitudinal and seasonal variations in growing conditions.
The temperature responses of An, its biochemical and biophysical limitations, and Rd were measured in 1-year-old needles of seedlings from the seed sources growing in eight forest plantations along a regional thermal gradient of 5.5 °C in Québec, Canada.
The average optimum temperature (Topt) for An was 19 ± 1.2 °C and was similar among seed sources and plantation sites along the thermal gradient. Net photosynthesis at Topt (Aopt) varied significantly among plantation sites and was quadratically related to the mean July temperature (MJT) of plantation sites. Topt for mesophyll conductance, maximum electron transport rate and maximum rate of carboxylation were 28, 22 and 30 °C, respectively. Basal respiration rate (Rd at 10 °C) was linearly and negatively associated with MJT. Q10 of Rd (the rate of change in Rd with a 10 °C increase in temperature) did not show any significant relationship with MJT and averaged 1.5 ± 0.1. The two seed sources were similar in their thermal responses to latitudinal and seasonal variations in growing conditions.
The results showed moderate thermal acclimation of respiration and no evidence for thermal acclimation of photosynthesis or local genetic adaptation for traits related to thermal acclimation. Therefore, growth of local white spruces may decline in future climates.
Full text
Available for:
BFBNIB, IZUM, KILJ, NMLJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
The boreal forest accounts for one-third of global forests, but remains largely inaccessible to ground-based measurements and monitoring. It contains large quantities of carbon in its vegetation and ...soils, and research suggests that it will be subject to increasingly severe climate-driven disturbance. We employ a suite of ground-, airborne- and space-based measurement techniques to derive the first satellite LiDAR-based estimates of aboveground carbon for the entire circumboreal forest biome. Incorporating these inventory techniques with uncertainty analysis, we estimate total aboveground carbon of 38±3.1Pg. This boreal forest carbon is mostly concentrated from 50 to 55° N in eastern Canada and from 55 to 60° N in eastern Eurasia. Both of these regions are expected to warm >3°C by 2100, and monitoring the effects of warming on these stocks is important to understanding its future carbon balance. Our maps establish a baseline for future quantification of circumboreal carbon and the described technique should provide a robust method for future monitoring of the spatial and temporal changes of the aboveground carbon content.
•We used ground measurements with LiDAR to map circumboreal forest carbon with errors.•We estimate total aboveground carbon of 38 ± 3.1 petagrams.•Carbon is concentrated from 50–55° N in eastern Canada and from 55–60° N in eastern Eurasia.•We establish a baseline and provide a robust method for future monitoring of aboveground carbon.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK