Modelling temporal variability in the carbon balance of a spruce/moss boreal forest

A model of the daily carbon balance of a black spruce/feathermoss boreal forest ecosystem was developed and results compared to preliminary data from the 1994 BOREAS field campaign in northem Manitoba, Canada. The model, driven by daily weather conditions, simulated daily soil climate status (temperature and moisture profiles), spruce photosynthesis and respiration, moss photosynthesis and respiration, and litter decomposition. Model agreement with preliminary field data was good for net ecosystem exchange (NEE), capturing both the asymmetrical seasonality and short-term variability. During the growing season simulated daily NEE ranged from -4 g C m^-2 d^-1 (carbon uptake by ecosystem) to + 2 g C m^-2 d^-1 (carbon flux to atmosphere), with fluctuations from day to day. In the early winter simulated NEE values were + 0.5 g C m^-2 d^-1, dropping to + 0.2 g C m^-2 d^-1 in mid-winter. Simulated soil respiration during the growing season (+ 1 to + 5 g C m^-2 d^-1) was dominated by metabolic respiration of the live moss, with litter decomposition usually contributing less than 30% and live spruce root respiration less than 10% of the total. Both spruce and moss net primary productivity (NPP) rates were higher in early summer than late summer. Simulated annual NEE for 1994 was -51 g C m^-2 y^-1, with 83% going into tree growth and 17% into the soil carbon accumulation. Moss NPP (58 g C m^-2 y^-1) was considered to be litter (i.e. soil carbon input; no net increase in live moss biomass). Ecosystem respiration during the snow-covered season (84 g C m^-2) was 58% of the growing season net carbon uptake. A simulation of the same site for 1968–1989 showed = 10–20% year-to-year variability in heterotrophic respiration (mean of + 113 g C m^-2 y^-1). Moss NPP ranged from 19 to 114 g C m^-2 y^-1; spruce NPP from 81 to 150 g C m^-2 y^-1; spruce growth (NPP minus litterfall) from 34 to 103 g C m^-2 y^-1; NEE ranged from +37 to -142 g C m^-2 y^-1. Values for these carbon balance terms in 1994 were slightly smaller than the 1969–89 means. Higher ecosystem productivity years (more negative NEE) generally had early springs and relatively wet summers; lower productivity years had late springs and relatively dry summers.     

Seasonal variations of CO2 and water vapour exchange rates over a temperate deciduous forest

Long-term and direct measurements of CO₂ and water vapour exchange are needed over forested ecosystems to determine their net annual fluxes of carbon dioxide and water. Such measurements are also needed to parameterize and test biogeochemical, ecological and hydrological assessment models. Responding to this need, eddy covariance measurements of CO₂ and water vapour were made ever a deciduous forest growing near Oak Ridge, TN, between April 1993 and April 1994. Periodic measurements were made of leaf area index, stomatal resistance, soil moisture and pre-dawn leaf water potential to characterize the gas exchange capacity of the canopy. Four factors had a disproportionate influence on the seasonal variation of CO₂ flux densities. These factors were photon flux densities (during the growing season), temperature (during the dormant season), leaf area index and the occurrence of drought The drought period occurred during the peak of the growing season and caused a significant decline in daily and hourly CO₂ flux densities, relative to observations over the stand when soil moisture was plentiful. The annual net uptake of carbon was calculated by integrating flux measurements and filling missing and spurious data with the relations obtained between measured CO₂ fluxes and environmental forcing variables. The net flux of carbon for the period between April 1993 and April 1994 was -525 g C m^?2 y^?1. This value represents a net flux of carbon from the atmosphere and into the forest. The net annual carbon exchange of this southern temperate broadleaved forest exceeded values measured over a northern temperate forest (which experiences a shorter growing season and has less leaf area) by 200 g C m^?2 y^?1 (cf. Wofsy et al 1993). The seasonal variation of canopy evaporation (latent heat flux) was controlled mostly by changes in leaf area and net radiation. A strong depression in evaporation rates was not observed during the drought Over a broadleaved forest large vapour pressure deficits promote evaporation and trees in a mixed stand are able to tap a variety of deep and shallow water sources.     

Annual cycles of water vapour and carbon dioxide fluxes in and above a boreal aspen forest

Water vapour and CO₂ fluxes were measured using the eddy correlation method above and below the overstorey of a 21-m tall aspen stand in the boreal forest of central Saskatchewan as part of the Boreal Ecosystem-Atmosphere Study (BOREAS). Measurements were made at the 39.5-m and 4-m heights using 3-dimensional sonic anemometers (Kaijo-Denki and Solent, respectively) and closed-path gas analysers (LI-COR 6262) with 6-m and 4.7-m long heated sampling tubing, respectively. Continuous measurements were made from early October to mid-November 1993 and from early February to late-September 1994. Soil CO₂ flux (respiration) was measured using a LI-COR 6000-09 soil chamber and soil evaporation was measured using Iysimetry. The leaf area index of the aspen and hazelnut understorey reached 1.8 and 3.3, respectively. The maximum daily evapotranspiration (E) rate was 5–6 mm d^?1. Following leaf-out the hazelnut and soil accounted for 22% of the forest E. The estimated total E was 403 mm for 1994. About 88% of the precipitation in 1994 was lost as evapotranspiration. During the growing season, the magnitude of half-hourly eddy fluxes of CO₂ from the atmosphere into the forest reached 1.2 mg CO₂ m^?2 s^?1 (33 μmol C m^?2 s^?1) during the daytime. Downward eddy fluxes at the 4-m height were observed when the hazelnut was growing rapidly in June and July. Under well-ventilated night-time conditions, the eddy fluxes of CO₂ above the aspen and hazelnut, corrected for canopy storage, increased exponentially with soil temperature at the 2-cm depth. Estimates of daytime respiration rates using these relationships agreed well with soil chamber measurements. During the 1994 growing season, the cumulative net ecosystem exchange (NEE) was -3.5 t C ha^?1 y^?1 (a net gain by the system). For 1994, cumulative NEE, ecosystem respiration (R) and gross ecosystem photosynthesis (GEP = R - NEE) were estimated to be -1.3, 8.9 and 10.2 t C ha^?1 y^?1 respectively. Gross photosynthesis of the hazelnut was 32% of GEP.     

Multiple Scales of Temporal Variability in Ecosystem Metabolism Rates: Results from 2 Years of Continuous Monitoring in a Forested Headwater Stream

Headwater streams are key sites of nutrient and organic matter processing and retention, but little is known about temporal variability in gross primary production (GPP) and ecosystem respiration (ER) rates as a result of the short duration of most metabolism measurements in lotic ecosystems. We examined temporal variability and controls on ecosystem metabolism by measuring daily rates continuously for 2 years in Walker Branch, a first-order deciduous forest stream. Four important scales of temporal variability in ecosystem metabolism rates were identified: (1) seasonal, (2) day-to-day, (3) episodic (storm-related), and (4) inter-annual. Seasonal patterns were largely controlled by the leaf phenology and productivity of the deciduous riparian forest. Walker Branch was strongly net heterotrophic throughout the year with the exception of the open-canopy spring when GPP and ER rates were co-equal. Day-to-day variability in weather conditions influenced light reaching the streambed, resulting in high day-to-day variability in GPP particularly during spring (daily light levels explained 84% of the variance in daily GPP in April). Episodic storms depressed GPP for several days in spring, but increased GPP in autumn by removing leaves shading the streambed. Storms depressed ER initially, but then stimulated ER to 2–3 times pre-storm levels for several days. Walker Branch was strongly net heterotrophic in both years of the study, with annual GPP being similar (488 and 519 g O₂ m⁻² y⁻¹ or 183 and 195 g C m⁻² y⁻¹) but annual ER being higher in 2004 than 2005 (−1,645 vs. −1,292 g O₂ m⁻² y⁻¹ or −617 and −485 g C m⁻² y⁻¹). Inter-annual variability in ecosystem metabolism (assessed by comparing 2004 and 2005 rates with previous measurements) was the result of the storm frequency and timing and the size of the spring macroalgal bloom. Changes in local climate can have substantial impacts on stream ecosystem metabolism rates and ultimately influence the carbon source and sink properties of these important ecosystems.     

Seasonal patterns and environmental control of carbon dioxide and water vapour exchange in an ecotonal boreal forest

Carbon dioxide, water vapour, and sensible heat fluxes were measured above and within a spruce dominated forest near the southern ecotone of the boreal forest in Maine, USA. Summer, mid-day carbon dioxide uptake was higher than at other boreal coniferous forests, averaging about – 13 μmol CO₂ m^–2 s^–1. Nocturnal summer ecosystem respiration averaged ≈ 6 μmol CO₂ m^–2 s^–1 at a mean temperature of ≈ 15 °C. Significant ecosystem C uptake began with the thawing of the soil in early April and was abruptly reduced by the first autumn frost in early October. Half-hourly forest CO₂ exchange was regulated mostly by the incident photosynthetically active photon flux density (PPFD). In addition to the threshold effects of freezing temperatures, there were seasonal effects on the inferred photosynthetic parameters of the forest canopy. The functional response of this forest to environmental variation was similar to that of other spruce forests. In contrast to reports of carbon loss from northerly boreal forest sites, in 1996 the Howland forest was a strong carbon sink, storing about 2.1 t C ha^–1.     

Biomass and carbon storage of the North American deciduous forest

Field measures of tree and shrub dimensions were used with established biomass equations in a stratified, two-stage cluster sampling design to estimate above-ground ovendry woody biomass and carbon storage of the eastern deciduous forest of North America. Biomass averaged 8.1 ± 1.4 (95% C.I.) kg/m² and totaled 18.1 ± 3.1 (95% C.I.) gigatons. Carbon storage averaged 3.6 ± 0.6 (95% C.I.) kg/m² and totaled 8.1 ± 1.4 (95% C.I.) gigatons. These values are lower than previous estimates commonly used in the analysis of the global carbon budget which range from 17.1 to 23.1 kg/m² for biomass and 7.7 to 10.4 kg/m² for carbon storage. These new estimates for the deciduous forest, together with earlier work in the boreal forest begin to reveal a pattern of overestimation of global carbon storage by vegetation in analyses of the global carbon budget. We discuss reasons for the differences between the new and earlier estimates, as well as implications for our understanding of the global carbon cycle.     

The Net Landscape Carbon Balance—Integrating terrestrial and aquatic carbon fluxes in a managed boreal forest landscape in Sweden

The boreal biome exchanges large amounts of carbon (C) and greenhouse gases (GHGs) with the atmosphere and thus significantly affects the global climate. A managed boreal landscape consists of various sinks and sources of carbon dioxide (CO₂), methane (CH₄), and dissolved organic and inorganic carbon (DOC and DIC) across forests, mires, lakes, and streams. Due to the spatial heterogeneity, large uncertainties exist regarding the net landscape carbon balance (NLCB). In this study, we compiled terrestrial and aquatic fluxes of CO₂, CH₄, DOC, DIC, and harvested C obtained from tall‐tower eddy covariance measurements, stream monitoring, and remote sensing of biomass stocks for an entire boreal catchment (~68 km²) in Sweden to estimate the NLCB across the land–water–atmosphere continuum. Our results showed that this managed boreal forest landscape was a net C sink (NLCB = 39 g C m^?2 year^?1) with the landscape–atmosphere CO₂ exchange being the dominant component, followed by the C export via harvest and streams. Accounting for the global warming potential of CH₄, the landscape was a GHG sink of 237 g CO₂‐eq m^?2 year^?1, thus providing a climate‐cooling effect. The CH₄ flux contribution to the annual GHG budget increased from 0.6% during spring to 3.2% during winter. The aquatic C loss was most significant during spring contributing 8% to the annual NLCB. We further found that abiotic controls (e.g., air temperature and incoming radiation) regulated the temporal variability of the NLCB whereas land cover types (e.g., mire vs. forest) and management practices (e.g., clear‐cutting) determined their spatial variability. Our study advocates the need for integrating terrestrial and aquatic fluxes at the landscape scale based on tall‐tower eddy covariance measurements combined with biomass stock and stream monitoring to develop a holistic understanding of the NLCB of managed boreal forest landscapes and to better evaluate their potential for mitigating climate change.     

Seasonal net carbon dioxide exchange of a beech forest with the atmosphere

The seasonal carbon dioxide exchange of a beech forest of Central Italy was studied by means of the eddy covariance technique. Additional measurements of biomass respiration with cuvettes and relationship of carbon dioxide exchanges with temperature and light were used to interpolate missing data during the dormant and part of the growing season. The net ecosystem production of the forest equals 472 g C m^?2 y^?1 while the gross ecosystem production 1016 g C m^?2 y^?1 and respiration 544 g C m^?2 y^?1. These estimates are compared with the net primary production determined by direct biomass sampling which amounts to 802 g C m^?2 y^?1.     

Carbon isotope composition of boreal plants: functional grouping of life forms

 We tested the hypothesis that life forms (trees, shrubs, forbs, and mosses; deciduous or evergreen) can be used to group plants with similar physiological characteristics. Carbon isotope ratios (δ¹³C) and carbon isotope discrimination (Δ) were used as functional characteristics because δ¹³C and Δ integrate information about CO₂ and water fluxes, and so are useful in global change and scaling studies. We examined δ¹³C values of the dominant species in three boreal forest ecosystems: wet Picea mariana stands, mesic Populus tremuloides stands, and dry Pinus banksiana stands. Life form groups explained a significant fraction of the variation in leaf carbon isotope composition; seven life-form categories explained 50% of the variation in δ¹³C and 42% of the variation in Δ and 52% of the variance not due to intraspecific genetic differences (n=335). The life forms were ranked in the following order based on their values: evergreen trees<deciduous trees=evergreen and deciduous shrubs=evergreen forbs<deciduous forbs=mosses. This ranking of the life forms differed between deciduous (Populus) and evergreen (Pinus and Picea) ecosystems. Furthermore, life forms in the Populus ecosystem had higher discrimination values than life forms in the dry Pinus ecosystem; the Picea ecosystem had intermediate Δ values. These correlations between Δ and life form were related to differences in plant stature and leaf longevity. Shorter plants had lower Δ values than taller plants, resulting from reduced light intensity at lower levels in the forest. After height differences were accounted for, deciduous leaves had higher discrimination values than evergreen leaves, indicating that deciduous leaves maintained higher ratios of intracellular to ambient CO₂ (c _i/c _a) than did evergreen leaves in a similar environment within these boreal ecosystems. We found the same pattern of carbon isotope discrimination in a year with above-average precipitation as in a year with below-average precipitation, indicating that environmental fluctuations did not affect the ranking of life forms. Furthermore, plants from sites near the northern and southern boundaries of the boreal forest had similar patterns of discrimination. We concluded that life forms are robust indicators of functional groups that are related to carbon and water fluxes within boreal ecosystems. Received: 15 April 1996 / Accepted: 16 November 1996     

Carbon, Water, and Energy Exchanges of a Hybrid Poplar Plantation During the First Five Years Following Planting

Eddy covariance was used to measure above-canopy exchanges of CO₂ and water vapor at an operational plantation of hybrid poplar (variety ??Walker??) established on marginal agricultural land in east central Alberta, Canada. Winter ecosystem respiration (R _e) rates were inferred from seasonal changes in the normalized respiration rate at 10°C (R ₁₀) for the growing season and observations of soil CO₂ concentration measured with solid-state probes. Over five consecutive growing seasons following planting, gross ecosystem production (GEP) increased each year, ranging from 21?g?C?m^?2?y^?1 in year 1 to 469?g?C?m^?2?y^?1 in year 5. During this period, the annual carbon balance shifted from a net source of greater than 330?g?C?m^?2 in year 1 to approximately C-neutral in year 5. Total carbon (C) release over 5?years likely exceeded 630?g?C?m^?2. Intra- and interannual variations in temperature and soil water availability greatly affected annual C balance each year. GEP and R _e were particularly sensitive to temperature during spring and to soil water availability in summer: year 5 was notable because a cold spring and accumulating drought caused growth and carbon uptake to fall well below their potential. Annual evapotranspiration (ET) increased slightly with leaf area, from 281?mm in year 1 to 323?mm in year 4, but in year 5 it declined, while exceeding total precipitation (P). This trend of increasing annual ET/P suggests that annual GEP could become increasingly water-limited in years with below normal precipitation, as the plantation achieves maximum leaf area. Measured canopy albedos did not change appreciably over three winters, suggesting that estimates of increased radiative forcing resulting from afforestation in high latitudes could be exaggerated in regions where fast-growing deciduous plantations are managed on short (~20-year) rotations.     

Carbon Balance in Forest Ecosystems of Southern Part of Moscow Region under a Rising Aridity of Climate

This study estimates carbon balance in a mixed mature forest on sod-podzolic sandy-loamy soil (Albeluvisols sandy, the Prioksko-Terrasny Nature Reserve) and in a secondary deciduous forest at the Experimental Field Station of the Institute of Physicochemical and Biological Problems of Soil Sciences, Russian Academy of Sciences (gray forest loamy soil, Luvisols loamy). CO₂ emissions from soils have been continuously measured every week since 1998. Net primary production was estimated in 2000–2014 by remote sensing using the normalized difference vegetation index. Long-term weather monitoring has revealed a distinctive trend toward increasing aridity of climate in the southern part of Moscow region in the observation period (1998–2014). Based on long-term ground-based and satellite monitoring data, this study shows that in the growing season, mixed and deciduous forests of the southern part of Moscow region are the sink of carbon with a mean flux of 41–112 g C m^–2, depending on the contribution of root respiration. Taking into account the CO₂ emissions from soils during the cold season, the forests are very likely to function as sources of atmospheric carbon at an amount of 30–100 g C m^–2 yr^–1, sometimes reaching very significant values of C flux (170–300 g C m^–2 yr^–1). In mature forest ecosystems, a significant influence on the quantitative estimate of the C balance is hampered by the CO₂ emission activity from coarse woody debris, which can reach up to 14% of the total losses of C during the decomposition of soil organic matter in the mixed forest, which turns it into a persistent source of CO₂ to the atmosphere. It is shown that the sink function of the forest ecosystems was more pronounced in dry years, whereas the excessive wetness diminishes their sink potential, turning the forests into sources of carbon dioxide.     

Seasonal variation in leaf properties and ecosystem carbon budget in a cool-temperate deciduous broad-leaved forest: simulation analysis at Takayama site, Japan

Seasonal changes in gross primary production (GPP) and net ecosystem production (NEP) in temperate deciduous forests are mostly driven by environmental conditions and the phenology of leaf demography. This study addresses another factor, temporal changes in leaf properties, i.e., leaf aging from emergence to senescence. A process-based model was used to link the ecosystem-scale carbon budget with leaf-level properties on the basis of field observation and scaling procedures; temporal variations in leaf thickness (leaf mass per area, LMA), photosynthetic rubisco (V_cmax) and electron-transport (J_max) capacity, and dark respiration (R_d) were empirically parameterized. The model was applied to a cool-temperate deciduous broad-leaved forest at Takayama, in central Japan, and validated with data of net ecosystem CO₂ exchange (NEE=–NEP) measured using the eddy-covariance method. NEP of the Takayama site varied seasonally from 3 g C m^–2 day^–1 net source in late winter to 5 g C m^–2 day^–1 net sink in early to mid-summer. A sensitivity experiment showed that removing the leaf-aging effect changed the seasonal CO₂ exchange pattern, and led to overestimation of annual GPP by 6% and annual NEP by 38%. We found that seasonal variation in V_cmax affected the seasonal pattern and annual budget of CO₂ exchange most strongly; LMA and R_d had moderate influences. The rapid change in V_cmax and R_d during leaf emergence and senescence was important in evaluating GPP and NEP of the temperate deciduous forest.     

An initial intercomparison of micrometeorological and ecological inventory estimates of carbon exchange in a mid-latitude deciduous forest

The role of mid‐latitude forests in the sequestration of carbon (C) is of interest to an increasing number of scientists and policy‐makers alike. Net CO₂ exchange can be estimated on an annual basis, using eddy‐covariance techniques or from ecological inventories of C fluxes to and from a forest. Here we present an intercomparison of annual estimates of C exchange in a mixed hardwood forest in the Morgan‐Monroe State Forest, Indiana, USA for two years, 1998 and 1999. Based on eddy‐covariance measurements made at 1.8 times canopy height from a tower, C uptake by the forest was 237 and 287 g C m^?2 y^?1 for 1998 and 1999, respectively. For the same time period, biometric and ecophysiological measures and modelled estimates of all significant carbon fluxes within deciduous forests were made, including: change in living biomass, aboveground and belowground detritus production, foliage consumption, and forest floor and soil respiration. Using this ecological inventory method for these same two time periods, C uptake was estimated to be 271 and 377 g C m^?2 y^?1, which are 14.3% and 31.4% larger, respectively, than the tower‐based values. The relative change between this method's annual estimates is consistent with that of the eddy‐covariance based values. Our results indicate that the difference in annual C exchange rates was due to reduced heterotrophic soil respiration in 1999.     

Exceptional carbon uptake in European forests during the warm spring of 2007: a data–model analysis

Temperate and boreal forests undergo drastic functional changes in the springtime, shifting within a few weeks from net carbon (C) sources to net C sinks. Most of these changes are mediated by temperature. The autumn 2006–winter 2007 record warm period was followed by an exceptionally warm spring in Europe, making spring 2007 a good candidate for advances in the onset of the photosynthetically active period. An analysis of a decade of eddy covariance data from six European forests stands, which encompass a wide range of functional types (broadleaf evergreen, broadleaf deciduous, needleleaf evergreen) and a wide latitudinal band (from 44° to 62°N), revealed exceptional fluxes during spring 2007. Gross primary productivity (GPP) of spring 2007 was the maximum recorded in the decade examined for all sites but a Mediterranean evergreen forest (with a +40 to +130 gC m^?2 anomaly compared with the decadal mean over the January–May period). Total ecosystem respiration (TER) was also promoted during spring 2007, though less anomalous than GPP (with a +17 to +93 gC m^?2 anomaly over 5 months), leading to higher net uptake than the long‐term mean at all sites (+12 to +79 gC m^?2 anomaly over 5 months). A correlative analysis relating springtime C fluxes to simple phenological indices suggested spring C uptake and temperatures to be related. The CASTANEA process‐based model was used to disentangle the seasonality of climatic drivers (incoming radiation, air and soil temperatures) and biological drivers (canopy dynamics, thermal acclimation of photosynthesis to low temperatures) on spring C fluxes along the latitudinal gradient. A sensitivity analysis of model simulations evidenced the roles of (i) an exceptional early budburst combined with elevated air temperature in deciduous sites, and (ii) an early relief of winter thermal acclimation in coniferous sites for the promotion of 2007 spring assimilation.     

Net primary production and net ecosystem production of a boreal black spruce wildfire chronosequence

Net primary production (NPP) was measured in seven black spruce (Picea mariana (Mill.) BSP)‐dominated sites comprising a boreal forest chronosequence near Thompson, Man., Canada. The sites burned between 1998 and 1850, and each contained separate well‐ and poorly drained stands. All components of NPP were measured, most for 3 consecutive years. Total NPP was low (50–100 g C m^?2 yr^?1) immediately after fire, highest 12–20 years after fire (332 and 521 g C m^?2 yr^?1 in the dry and wet stands, respectively) but 50% lower than this in the oldest stands. Tree NPP was highest 37 years after fire but 16–39% lower in older stands, and was dominated by deciduous seedlings in the young stands and by black spruce trees (>85%) in the older stands. The chronosequence was unreplicated but these results were consistent with 14 secondary sites sampled across the landscape. Bryophytes comprised a large percentage of aboveground NPP in the poorly drained stands, while belowground NPP was 0–40% of total NPP. Interannual NPP variability was greater in the youngest stands, the poorly drained stands, and for understory and detritus production. Net ecosystem production (NEP), calculated using heterotrophic soil and woody debris respiration data from previous studies in this chronosequence, implied that the youngest stands were moderate C sources (roughly, 100 g C m^?2 yr^?1), the middle‐aged stands relatively strong sinks (100–300 g C m^?2 yr^?1), and the oldest stands about neutral with respect to the atmosphere. The ecosystem approach employed in this study provided realistic estimates of chronosequence NPP and NEP, demonstrated the profound impact of wildfire on forest–atmosphere C exchange, and emphasized the need to account for soil drainage, bryophyte production, and species succession when modeling boreal forest C fluxes.     

Paired-tower measurements of carbon and energy fluxes following disturbance in the boreal forest

Disturbances by fire and harvesting are thought to regulate the carbon balance of the Canadian boreal forest over scales of several decades. However, there are few direct measurements of carbon fluxes following disturbances to provide data needed to refine mathematical models. The eddy covariance technique was used with paired towers to measure fluxes simultaneously at disturbed and undisturbed sites over periods of about one week during the growing season in 1998 and 1999. Comparisons were conducted at three sites: a 1‐y‐old burned jackpine stand subjected to an intense crown fire at the International Crown Fire Modelling Experiment site near Fort Providence, North‐west Territories; a 1‐y‐old clearcut aspen area at the EMEND project near Peace River, Alberta; and a 10‐y‐old burned, mixed forest near Prince Albert National Park, Saskatchewan. Nearby mature forest stands of the same types were also measured as controls. The harvested site had lower net radiation (R_n), sensible (H) and latent (LE) heat fluxes, and greater ground heat fluxes (G) than the mature forest. Daytime CO₂ fluxes were much reduced, but night‐time CO₂ fluxes were identical to that of the mature aspen forest. It is hypothesized that the aspen roots remained alive following harvesting, and dominated soil respiration. The overall effect was that the harvested site was a carbon source of about 1.6 gC m^?2 day^?1, while the mature site was a sink of about ?3.8 gC m^?2 day^?1. The one‐year‐old burn had lower R_n, H and LE than the mature jackpine forest, and had a continuous CO₂ efflux of about 0.8 gC m^–2 day^?1 compared to the mature forest sink of ? 0.5 g C m^?2 day^?1. The carbon source was likely caused by decomposition of fire‐killed vegetation. The 10‐y‐old burned site had similar H, LE, and G to the mature mixed forest site. Although the diurnal amplitude of the CO₂ fluxes were slightly lower at the 10‐y‐old site, there was no significant difference between the daily integrals (? 1.3 gC m^?2 day^?1 at both sites). It appears that most of the change in carbon flux occurs within the first 10 years following disturbance, but more data are needed on other forest and disturbance types for the first 20 years following the disturbance event.     

Productivity of forests in the Eurosiberian boreal region and their potential to act as a carbon sink –- a synthesis

Based on review and original data, this synthesis investigates carbon pools and fluxes of Siberian and European forests (600 and 300 million ha, respectively). We examine the productivity of ecosystems, expressed as positive rate when the amount of carbon in the ecosystem increases, while (following micrometeorological convention) downward fluxes from the atmosphere to the vegetation (NEE = Net Ecosystem Exchange) are expressed as negative numbers. Productivity parameters are Net Primary Productivity (NPP=whole plant growth), Net Ecosystem Productivity (NEP = CO₂ assimilation minus ecosystem respiration), and Net Biome Productivity (NBP = NEP minus carbon losses through disturbances bypassing respiration, e.g. by fire and logging). Based on chronosequence studies and national forestry statistics we estimate a low average NPP for boreal forests in Siberia: 123 gC m^–2 y^–1. This contrasts with a similar calculation for Europe which suggests a much higher average NPP of 460 gC m^–2 y^–1 for the forests there. Despite a smaller area, European forests have a higher total NPP than Siberia (1.2–1.6 vs. 0.6–0.9 × 10¹⁵ gC region^–1 y^–1). This arises as a consequence of differences in growing season length, climate and nutrition. For a chronosequence of Pinus sylvestris stands studied in central Siberia during summer, NEE was most negative in a 67-y old stand regenerating after fire (– 192 mmol m^–2 d^–1) which is close to NEE in a cultivated forest of Germany (– 210 mmol m^–2 d^–1). Considerable net ecosystem CO₂-uptake was also measured in Siberia in 200- and 215-y old stands (NEE:174 and – 63 mmol m^–2 d^–1) while NEP of 7- and 13-y old logging areas were close to the ecosystem compensation point. Two Siberian bogs and a bog in European Russia were also significant carbon sinks (– 102 to – 104 mmol m^–2 d^–1). Integrated over a growing season (June to September) we measured a total growing season NEE of – 14 mol m^–2 summer^–1 (– 168 gC m^–2 summer^–1) in a 200-y Siberian pine stand and – 5 mol m^–2 summer^–1 (– 60 gC m^–2 summer^–1) in Siberian and European Russian bogs. By contrast, over the same period, a spruce forest in European Russia was a carbon source to the atmosphere of (NEE: + 7 mol m^–2 summer^–1 = + 84 gC m^–2 summer^–1). Two years after a windthrow in European Russia, with all trees being uplifted and few successional species, lost 16 mol C m^–2 to the atmosphere over a 3-month in summer, compared to the cumulative NEE over a growing season in a German forest of – 15.5 mol m^–2 summer^–1 (– 186 gC m^–2 summer^–1; European flux network annual averaged – 205 gC m^–2 y^–1). Differences in CO₂-exchange rates coincided with differences in the Bowen ratio, with logging areas partitioning most incoming radiation into sensible heat whereas bogs partitioned most into evaporation (latent heat). Effects of these different surface energy exchanges on local climate (convective storms and fires) and comparisons with the Canadian BOREAS experiment are discussed. Following a classification of disturbances and their effects on ecosystem carbon balances, fire and logging are discussed as the main processes causing carbon losses that bypass heterotrophic respiration in Siberia. Following two approaches, NBP was estimated to be only about 13–16 mmol m^–2 y^–1 for Siberia. It may reach 67 mmol m^–2 y^–1 in North America, and about 140–400 mmol m^–2 y^–1 in Scandinavia. We conclude that fire speeds up the carbon cycle, but that it results also in long-term carbon sequestration by charcoal formation. For at least 14 years after logging, regrowth forests remain net sources of CO₂ to the atmosphere. This has important implications regarding the effects of Siberian forest management on atmospheric concentrations. For many years after logging has taken place, regrowth forests remain weaker sinks for atmospheric CO₂ than are nearby old-growth forests.     

A long-term record of carbon exchange in a boreal black spruce forest: means, responses to interannual variability, and decadal trends

We present a decadal (1994–2004) record of carbon dioxide flux in a 160‐year‐old black spruce forest/veneer bog complex in central Manitoba, Canada. The ecosystem shifted from a source (+41 g C m⁻², 1995) to a sink (−21 g C m⁻², 2004) of CO₂ over the decade, with an average net carbon balance near zero. Annual mean temperatures increased 1–2° during the period, consistent with the decadal trend across the North American boreal biome. We found that ecosystem carbon exchange responded strongly to air temperature, moisture status, potential evapotranspiration, and summertime solar radiation. The seasonal cycle of ecosystem respiration significantly lagged that of photosynthesis, limited by the rate of soil thaw and the slow drainage of the soil column. Factors acting over long time scales, especially water table depth, strongly influenced the carbon budget on annual time scales. Net uptake was enhanced and respiration inhibited by multiple years of rainfall in excess of evaporative demand. Contrary to expectations, we observed no correlation between longer growing seasons and net uptake, possibly because of offsetting increases in ecosystem respiration. The results indicate that the interactions between soil thaw and water table depth provide critical controls on carbon exchange in boreal forests underlain by peat, on seasonal to decadal time scales, and these factors must be simulated in terrestrial biosphere models to predict response of these regions to future climate.     

Net ecosystem productivity of boreal jack pine stands regenerating from clearcutting under current and future climates

Life cycle analysis of climate and disturbance effects on forest net ecosystem productivity (NEP) is necessary to assess changes in forest carbon (C) stocks under current or future climates. Ecosystem models used in such assessments need to undergo well-constrained tests of their hypotheses for climate and disturbance effects on the processes that determine CO₂ exchange between forests and the atmosphere. We tested the ability of the model ecosys to simulate diurnal changes in CO₂ fluxes under changing air temperatures (T_a) and soil water contents during forest regeneration with eddy covariance measurements over boreal jack pine (Pinus banksiana) stands along a postclearcut chronosequence. Model hypotheses for hydraulic and nutrient constraints on CO₂ fixation allowed ecosys to simulate the recovery of C cycling during the transition of boreal jack pine stands from C sources following clearcutting (NEP from −150 to −200 g C m⁻² yr⁻¹) to C sinks at maturity (NEP from 20 to 80 g C m⁻² yr⁻¹) with large interannual variability. Over a 126-year logging cycle, annualized NEP, C harvest, and net biome productivity (NBP=NEP–harvest removals) of boreal jack pine averaged 47, 33 and 14 g C m⁻² yr⁻¹. Under an IPCC SRES climate change scenario, rising T_a exacerbated hydraulic constraints that adversely affected NEP of boreal jack pine after 75 years. These adverse effects were avoided in the model by replacing the boreal jack pine ecotype with one adapted to warmer T_a. This replacement raised annualized NEP, C harvest, and NBP to 81, 56 and 25 g C m⁻² yr⁻¹ during a 126-year logging cycle under the same climate change scenario.     

A new seasonal‐deciduous spring phenology submodel in the Community Land Model 4.5: impacts on carbon and water cycling under future climate scenarios

A spring phenology model that combines photoperiod with accumulated heating and chilling to predict spring leaf‐out dates is optimized using PhenoCam observations and coupled into the Community Land Model (CLM) 4.5. In head‐to‐head comparison (using satellite data from 2003 to 2013 for validation) for model grid cells over the Northern Hemisphere deciduous broadleaf forests (5.5 million km²), we found that the revised model substantially outperformed the standard CLM seasonal‐deciduous spring phenology submodel at both coarse (0.9 × 1.25°) and fine (1 km) scales. The revised model also does a better job of representing recent (decadal) phenological trends observed globally by MODIS, as well as long‐term trends (1950–2014) in the PEP725 European phenology dataset. Moreover, forward model runs suggested a stronger advancement (up to 11 days) of spring leaf‐out by the end of the 21st century for the revised model. Trends toward earlier advancement are predicted for deciduous forests across the whole Northern Hemisphere boreal and temperate deciduous forest region for the revised model, whereas the standard model predicts earlier leaf‐out in colder regions, but later leaf‐out in warmer regions, and no trend globally. The earlier spring leaf‐out predicted by the revised model resulted in enhanced gross primary production (up to 0.6 Pg C yr^?1) and evapotranspiration (up to 24 mm yr^?1) when results were integrated across the study region. These results suggest that the standard seasonal‐deciduous submodel in CLM should be reconsidered, otherwise substantial errors in predictions of key land–atmosphere interactions and feedbacks may result.