Introduction

Climate change, caused by increasing atmospheric CO2 concentration resulting from fossil fuel use, is a pressing concern of international scope (IPCC 2014). Managing terrestrial carbon sinks to ameliorate the problem of increasing CO2 concentration is of prime interest (IPCC 2014). The world’s forests have been estimated to contain about 80% of all aboveground carbon and 40% of all belowground carbon (Dixon et al. 1994), and forests are recognized as important global carbon (C) sinks and sources (Landsberg and Gower 1997) that play a key role in the atmospheric carbon dynamics and climate change mitigation.

Around the globe, boreal forests cover 1220 million ha and are the world’s second largest forest biome (Landsberg and Gower 1997). It is therefore important to consider their management in relation to the goal of increasing carbon stocks and sequestration in forests. Forested land covers 348 million ha in Canada (Natural Resources Canada 2012), including 223 million ha of boreal forest (Brandt 2009), and of this total, mixedwood forest occupies some 152 million ha (Baldwin et al. 2012). Canada’s managed boreal forest occupies 54% of the nation’s total boreal forested area, with 28 Pg of carbon in biomass, organic matter and soil pools (Kurz et al. 2013). Boreal mixedwood stands in the province of Ontario, which are defined as those in which no single species exceeds 80% of basal area (MacDonald 1995), occupy 15.3 million ha or 53% of the forested land base (Towill 1996).

The commonest tree species in mature Canadian boreal mixedwood stands are trembling aspen (Populus tremuloides Michx.), white birch (Betula papyrifera Marsh.), black spruce [Picea mariana (Mill.) B.S.P.], white spruce [Picea glauca (Moench) Voss], and balsam fir [Abies balsamea (L.) Mill.]; these stands are disturbance-related and represent a mid-successional stage in forest development. Early stages are generally aspen-dominated and late stages are spruce-fir-birch dominated. The two natural disturbances that most commonly occur in boreal mixedwood stands are wildfire and spruce budworm (Choristoneura fumiferana Clem.) infestations, while manmade disturbances include harvesting and planting, which in some cases may lead to species conversion.

Boreal mixedwoods are complex forests in which different species accumulate biomass at different rates and utilize different resources at a site (Paré and Bergeron 1995) and have been considered the most productive forest ecosystems in North American boreal forest (Chen and Popadiouk 2002), although more recent research suggests that aspen-dominated stands may be more productive (Laganière et al. 2015). The diversity of species supported on a boreal mixedwood site suggests efficient use of the land available for tree growth, with each species exploiting its niche within the ecosystem such that total production on the site reflects the best growth that each tree species can offer (Légaré et al. 2004). Boreal mixedwoods often occupy the more fertile sites, and, because of their extensive areal coverage, play an important role in carbon sequestration in Canadian forests. Regionally, Liu et al. (2002) noted that Ontario’s boreal forests were relatively young, with an average age of 36.2 year in 1990, a further indication of their potential for C sequestration and storage. Knowledge of the total ecosystem carbon stocks in the extensive Canadian boreal mixedwood biome is therefore needed to design forest management strategies to enhance sequestration and storage at different stages of development and contribute to climate change mitigation. The need for further quantification of the climate change mitigation potential of Canada’s boreal zone and the benefits of implementing forestry management practices has been noted by Lemprière et al. (2013). Quantification of the carbon stocks and net primary productivity at different stages of boreal mixedwood stand development are also of value for carbon budget model development and validation (e.g., Li et al. 2002). Modelling Canada’s forest carbon content and projecting its’ overall effectiveness in sequestering carbon from the atmosphere is foundational for policy discussions around possible climate change mitigation measures.

In this study, we quantified ecosystem carbon stocks and aboveground net primary productivity in typical Ontario boreal mixedwood stands at different stages of development, including a mature stand (> 65 year) before and after harvesting, a juvenile stand replanted with softwood seedlings, and an intact mature stand. Specifically, our goals were to assess the effect of clear-cut harvesting on ecosystem carbon stocks and aboveground net primary productivity, to quantify ecosystem carbon stocks and aboveground net primary productivity in a juvenile stand replanted with softwood seedlings and in mature Ontario boreal mixedwood forest, and to assess stocks in detritus and mineral soil in a replanted juvenile stand and a mature mixedwood stand.

Materials and methods

Study sites

Study site details are provided in Table 1. The sites were located within about 40 km of each other, near Timmins in northeastern Ontario. The mature stand at the McKeown Lake site (48.166°N, 81.550°W) was sampled at age 69 years in 2006, and was then clear-cut harvested in late 2008 and early 2009, and re-assessed at age 2 in 2011. The juvenile stand, which was planted with spruce seedlings after clear-cut harvesting to achieve stand conversion (Childerhose Township site 48.133°N, 81.628°W), was quantified at ages 20 (2006) and 25 years (2011). Harvested boreal mixedwood stands are sometimes planted with coniferous seedlings to add commercially desirable species and enhance stand productivity. These two sites were approximately 2 km away from each other and shared many climatic and soil descriptive features. Both are ground moraines of silty to sandy tills and predominantly humo-ferric podzol (Bélisle 1980; Soil Classification Working Group 1998). Soil moisture is fresh, the topography is simple, and site elevation lies between 365 m and 380 m a.s.l. (Chapman and Thomas 1968). The other mature stand, at the Groundhog River site (48.217° N, 82.156° W), was measured at stand age 74 (2003) and 81 years (2010). The soil at this site is a sandy loam and characterised as a dystric brunisol with pockets of orthic gleysol. Soil moisture is fresh to very fresh, the topography is simple and flat, and site elevation is 340 m a.s.l.

Table 1 Study design

Site climate

The weather pattern at the sites is a modified continental climate. These sites typically experience 160 growing days, i.e., the number of days when the mean temperature is above 5.6 °C (42 °F). Growing days typically begin around May 3 and end around October 12 (Chapman and Thomas 1968). Normals for the period 1971–2000 are 831 mm of annual precipitation, with 558 mm as rainfall and 3.1 m as snowfall (Environment Canada 2013). These data were collected at Timmins Airport, about 40 km from the Childerhose and McKeown Lake sites and 70 km from the Groundhog River site.

Sampling methodology

Sampling plots were established to assess vegetation, woody debris and soil properties using the Canadian National Forest Inventory plot establishment protocol (National Forest Inventory 2008). Twenty-five sampling plots were distributed over the three sites, with 8 plots installed on a 300 m × 300 m grid at both the juvenile site and the harvested mature site, and 9 plots at the other mature site. Each sampling plot comprised a circular plot with 11.28 m radius (area 0.04 ha), in which all trees (i.e., diameter breast height [dbh] ≥ 9 cm) and stumps were measured, and a nested inner circular plot with 3.99 m radius (area 0.005 ha) that was used for measuring saplings (dbh < 9 cm). Within each sampling plot, four circular mini-plots with radius 0.56 m (area 0.0001 ha) were used for measuring shrubs and trees < 1.3 m in height, and destructive biomass sampling of shrubs, herbs and mosses. Woody debris was assessed along two 30-m transects intersecting at right angles at plot centre, following procedures described by Van Wagner (1968) and McRae et al. (1979). Size classes for fine, small and coarse woody debris refer to diameter ranges of < 1 cm, 1–7 cm and > 7 cm, respectively. In addition, eight litterfall traps (sampling area 0.25 m2) were sampled annually at each sampling plot over a 5-year period. A soil pit was dug at the installation of each sampling plot to sample the forest floor and mineral soil layers to a depth of 60 cm. The soil pit assessments were made at age 20 in the juvenile stand at the Childerhose site and at age 69 and 74 in the mature stands at McKeown Lake and Groundhog River sites. Stand characteristics at each site are provided in Tables 2 and 3. The mature sites were about 85% fully stocked, as determined by comparison of tree numbers, basal area and age for individual species.

Table 2 Mean (± SE) stand density and basal area
Table 3 Mean (± SD) tree diameter (dbh ≥ 9 cm)

Above- and belowground biomass were calculated using allometric equations of the form M = aDb, where M is dry mass of biomass component (kg), D is dbh (cm), and a and b are constants. This equation provided accurate biomass predictions with lower data requirements by precluding the need to measure plant height. Equations based on in-house measurements of trees and saplings, sampled in northern Ontario by Great Lakes Forestry Centre researchers, were available for most species. The equations used to calculate above- and belowground tree biomass for trembling aspen, white birch, white spruce, black spruce and balsam fir, with dbh ranges, are given in Table 4. Where regional biomass equations were not available, those reported in the literature were employed (Ohmann et al. 1981; Crow and Erdmann 1983; Smith and Brand 1983; Perala and Alban 1994; Roussopoulos and Loomis 1979; Ter-Mikaelian and Korzukhin 1997; Brassard et al. 2011). Oven-dried samples of litterfall, forest floor organic layer, shrubs, herbs and mosses and soil samples were weighed and underwent chemical analysis to quantify carbon content at the Great Lakes Forestry Centre Soils Laboratory. Carbon content of trees, saplings and woody debris was taken as 50% of biomass (Amiro et al. 2001). Annual aboveground net primary productivity (AGNPP) was estimated from measurements made at each site. The repeated measurements of trees and sapling dbh were employed with allometric equations to estimate the annual aboveground biomass increment, and combined with measurements of understorey biomass and annual litterfall, to estimate total annual AGNPP. Statistical analysis of the data provided mean values with standard errors or deviations for each variable at the site of interest, for the stand age at quantification.

Table 4 In-house allometric equation coefficients and coefficient of determination for commonest trees (AG: aboveground biomass, BG: belowground biomass)

Results

Carbon stocks

Aggregate stocks

Total carbon stocks (± SE) in the recently harvested stand were 134 ± 13 Mg C ha−1, with 27, 35 and 38% in vegetation, detritus and mineral soil, respectively (Table 5, Fig. 1a, b). Carbon stocks in the preharvest mature stand were reduced by 90 Mg C ha−1 or 40% by harvesting. In the juvenile stand of age 20, total stocks were 119 ± 13 Mg C ha−1, with 26, 34 and 40% in vegetation, detritus and mineral soil. During juvenile stand development the regrowth of the trees and saplings compensated for carbon loss due the decomposition of woody debris and stumps, and total carbon stocks increased by 2% p.a. to 130 Mg C ha−1 by age 25, with vegetation comprising 34%. Mature stands had an average of 208 Mg C ha−1 with 59, 16 and 25% in vegetation, detritus and mineral soil. Between age 74 and 81, total stocks increased by 13.3 Mg C ha−1 or 1% p.a. Combined stocks in the detritus and mineral soil layers were 89, 94 and 82 Mg C ha−1 at ages 20, 69 and 74, equivalent to 74, 42 and 43% of carbon in the stand. At all ages, these two components made up a large fraction of total stocks, to which the forest canopy added differing amounts according to stand age.

Table 5 Mean (± SE) carbon stocks in stands
Fig. 1
figure 1

Mixedwood stand carbon stocks versus age by substrate, and proportion by substrate: mixedwood stand species composition versus age, by stem count and by basal area

Vegetation components

Live and dead trees and saplings comprised 7, 43 and 126 Mg C ha−1, equivalent to 5.3, 42 and 61% of carbon stocks at 2, 25 and 81 year, respectively (Table 5, Fig. 1a, b). Of these, stocks in live trees and saplings contained 1.6, 41 and 108 Mg C ha−1, or 1.2, 33 and 52% of carbon stocks. The clear-cut harvest reduced stocks in live trees and saplings by 98%, from 117 to 1.6 Mg C ha−1 (Fig. 1a). In the juvenile stand, stocks in live trees and saplings increased by 71% from 24 to 41 Mg C ha−1, equivalent to 3.5 Mg C ha−1 p.a. (15% p.a.) between age 20 and 25, whereas in the mature stand, stocks in live trees and saplings increased from 101 to 108 Mg C ha−1, or about 1% p.a. between age 74 and 81.

Carbon stocks in shrubs and herbaceous vegetation contained 67, 3 and 0.8% of the total stocks in living biomass in the 2-, 25- and 81-year-old stands. Removal of the tree canopy during harvest resulted in adding biomass to this stand component, and at age 2, stocks in this component totalled 3.3 Mg C ha−1, or 2.5% of stocks, compared to 0.2% before harvest. In the juvenile stand, stocks in this component were 1.8 and 1.3 Mg C ha−1 at ages 20 and 25, or 1.5 and 1% of total carbon. Shrubs and ground vegetation contained 0.8 Mg C ha−1, or 0.4% of total stocks in the mature stand at age 74 and 81. As expected, this stand component only contributed a significant fraction of carbon stocks during the early stages of stand development.

Stumps contributed appreciably to carbon stocks only in the early stages of stand development after harvest; at age 2 stumps contained 26.1 Mg C ha−1 or 20% of total carbon (Table 5). Stumps were still present in the juvenile stand, with 3 and 0.5% of total stocks at ages 20 and 25. At age 69 and later, decomposition incorporated the stump biomass into the forest floor organic layer and mineral soil layer.

Detritus and soil components

Woody debris carbon was greatest after harvesting, containing 14.9 Mg C ha−1 or 11% of stocks. At age 2, coarse and small debris accounted for 93% of debris mass. In the juvenile stand, woody debris contained 4.0 and 0.9 Mg C ha−1 at age 20 and 25, or 6.4 and 1.2% of stocks. Woody debris carbon in mature stands was greater than in juvenile stands with average totals of 8 Mg C ha−1, or 4% of stocks. In the mature stands, coarse and small debris comprised an 87% of total woody debris mass. For all stand ages, fine woody debris ranged between 6.1 and 21% of total woody debris carbon stocks. Woody debris at the McKeown lake site increased from 11.5 Mg C ha−1 at age 69 to 14.9 Mg C ha−1 at age 2 due to the addition of downed branches and debris resulting from harvesting, although the difference was not statistically significant (t test, p = 0.05).

The forest floor organic layer contained 36.4, 30.9 and 22.1 Mg C ha−1 at ages 20, 69 and 74, or 31, 14 and 11% of stocks; mean layer depths (± SD) were 10.4 ± 3.3, 7.5 ± 3.0 and 9.3 ± 6.0 cm. The mineral soil layer to a depth of 60 cm averaged 51 Mg C ha−1, and comprised 41, 23 and 27% of stocks at age 20, 69 and 74 years. Mineral soil carbon in the 0–30 cm layer averaged 35 Mg C ha−1, between 15 and 32% of carbon stocks in the stand. Mean combined stocks (± SE) in the forest floor organic and mineral soil layers were 85 ± 12, 82 ± 12 and 75 ± 7 Mg C ha−1 at age 20, 69 and 74 year, equivalent to 71, 37 and 38% of total stocks.

Aboveground net primary productivity

In the harvested stand mean annual aboveground NPP (AGNPP) was 2.1 Mg C ha−1, compared to 3.7 and 2.6 in the juvenile and mature stand (Table 6). Thus, the juvenile stand was most productive, while the harvested and mature stands were adding to aboveground carbon at 59 and 73% of the rate in the juvenile stand. Tree and sapling growth comprised 7, 73 and 31% of annual AGNPP in the harvested, juvenile and mature stands. Annual litterfall in the harvested stand averaged 0.3 Mg C ha−1, compared to 1.1 in the juvenile stand and 1.8 in mature stands contributing 15, 30 and 68%, respectively. The largest contributors to annual AGNPP in the harvested, juvenile and mature stands were the understorey, tree growth and litterfall, respectively. As expected, litterfall in mature canopies was greater than in the harvested and juvenile stands (Table 6).

Table 6 Mean (± SE) annual aboveground net primary production

Forest stand characteristics

Clear-cut harvest reduced the stem density and basal area (BA) by over 90%, and dead trees accounted for 84% of postharvest BA (Table 2, Fig. 1c, d). Live trees constituted only 19% of the 100 ha−1 remaining in the harvested stand. In the replanted juvenile stand, the maturing saplings caused stem density and BA to more than triple between ages 20 and 25, from 250 to 910 ha−1, and 2.7 to 9.2 m2 ha−1. Mean spruce tree diameters were similar at year 20 and 25 as more saplings reached the tree dbh threshold (Table 3). In mature stands, the average stem density was 1167 (ha−1), with 92% live trees. Trembling aspen and white birch made up 51% of live stems, and white and black spruce and balsam fir accounted for a further 44% of stems. Basal area in the mature stands averaged 34 m2 ha−1, and total BA only increased by 10% between stand age 74 and 81, while mean live tree diameter increased by 5% between year 74 and 81; trembling aspen showed the greatest increase, of 14%, whereas spruce dbh only increased by 3% (Table 3).

Discussion

Boreal forest carbon stocks

Measured carbon stocks in boreal forest stands at various ages were generally higher in mixedwood than in coniferous stands, and mature trembling aspen stands also held higher stocks than coniferous stands (Fig. 2a). In addition, for a given BA, mixed wood stands generally had higher stocks than coniferous stands (Fig. 2b). The mature mixedwood stands studied by Martin et al. (2005), Morrison et al. (2001), Hazlett et al. (2005), Laganière et al. (2015) and those quantified in the present study contained carbon stocks (± SD) of 213 ± 52 Mg C ha−1, and had an average stand age (± SD) of 78 ± 14 year. Mature trembling aspen stands studied by Gower et al. (1997) and Laganiere et al. (2015) contained mean carbon stocks (± SD) of 239 ± 89 Mg C ha−1. Estimated overall carbon stocks in various managed Canadian boreal ecosystems ranged from 140 to 240 Mg C ha−1, with mean carbon density of 193 (Kurz et al. 2013). Thus, mature mixedwood and trembling aspen stands are typically at or above the average density for Canadian boreal ecosystems.

Fig. 2
figure 2

Total carbon stocks in various boreal forest stand types, by age and by basal area (MW mixedwood, TA trembling aspen, JP jack pine, BS black spruce: plus sign Payne, times sign Martin, asterisk Howard, open circle Gower, closed circle Hazlett, open square Laganiere Ontario, closed square Laganiere Quebec, open diamond Morrison (1993), closed diamond Morrison (2001), open inverted triangle Seedre et al. (2014) cut, closed inverted triangle Seedre et al. (2014) fire)

Vegetation carbon

Carbon stocks in vegetation for various boreal forest ecosystems typically are higher in this component of mixedwood and trembling aspen stands, as compared to jack pine and black spruce, for a wide range of ages and BA (Fig. 3a, b). In mature mixedwood stands over 60 years of age (mean age ± SD: 77 ± 9), the average carbon stocks (± SD) in vegetation were 114 ± 51 Mg C ha−1, and on average this comprised 53% of carbon stocks present. Mature jack pine stand vegetation held stocks of 87 ± 67 Mg C ha−1, representing 59% of total. As might be expected, measurements show a tighter correlation with BA as opposed to age, with roughly 40 Mg C ha−1 of vegetation carbon added for every 10 m2 ha−1 of BA (Fig. 3b).

Fig. 3
figure 3

Vegetation carbon stocks in various boreal forest stand types, by age and by basal area; and basal area by age [MW mixedwood, TA trembling aspen, JP jack pine, BS black spruce: plus sign Payne, times sign Martin, asterisk Howard, open circle Gower, closed circle Hazlett, open square Laganiere Ontario, closed square Laganiere Quebec, open diamond Morrison (1993), closed diamond Morrison (2001), open inverted triangle Seedre et al. (2014) cut, closed inverted triangle Seedre et al. (2014) fire, closed triangle Strukelj et al. (2015)]

Published basal area measurements in various boreal stand types at different ages show that mixedwood and trembling aspen stands typically have higher BA values than jack pine and black spruce stands of similar age (Fig. 3c). Mean basal area (± SD) for mature mixedwood stands was 36 ± 9 m2 ha−1 in studies by Hazlett et al. (2005), Laganiere et al. (2015), Martin et al. (2005), Morrison et al. (2001), and the present investigation, whereas in mature coniferous stands studied by Gower et al. (1997), Howard et al. (2004), Laganiere et al. (2015) and Morrison et al. (1993) mean BA was 31 ± 14 m2 ha−1.

Detritus and mineral soil carbon

Forest floor carbon stocks averaged 32 ± 13 Mg C ha−1 (± SD) in mixedwood stands as compared to 16 ± 6 in jack pine stands (Fig. 4a, b). In addition, these forest floor organic layer stocks were stable across the observed range of stand age and BA. This characteristic can simplify estimation of carbon stocks held in boreal stands of different ages when planning forest management for optimum carbon sequestration. Carbon stocks in detritus (woody debris plus the forest floor organic layer, Fig. 4c) were generally higher in mixedwood stands (42 ± 13 Mg C ha−1; ± SD), as compared to jack pine stands (19 ± 9 Mg C ha−1), and in mixedwood stands tended to decrease slightly with age. The overall mean (± SD) proportion of stocks found in detritus was 29 ± 12%, compared to 20 ± 9% in the forest floor organic layer alone.

Fig. 4
figure 4

Forest floor organic layer in various boreal forest stand types by age and by basal area, and detritus and mineral soil carbon stocks by age (MW mixedwood, TA trembling aspen, JP jack pine, BS black spruce: plus sign Payne, times sign Martin, asterisk Howard, open circle Gower, closed circle Hazlett, open square Laganiere Ontario, closed square Laganiere Quebec, circled vertical line Lee, open diamond Morrison (1993), closed diamond Morrison (2001), open inverted triangle Seedre et al. (2014) cut, closed inverted triangle Seedre et al. (2014) fire)

Litter degradation rate is found to vary with tree species, with litter from conifers being slower to degrade than that from deciduous species (Paré et al. 2006). Thus the forest floor layer at the juvenile site restocked with conifers would be expected to degrade more slowly than that at the mature stands comprising of similar proportions of softwood and hardwood trees. Woody debris input is found to increase at harvesting, as in the present study (Brais et al. 2004; Moroni 2006). Brais et al. found decomposition rates of woody debris of about 30% per year in both a harvested and a mature boreal mixedwood stand in Quebec, whereas in other studies decomposition took up to 30 years, depending on climatic conditions and species (Carlton and Pickford 1982; Prescott et al. 1989; Morrison 1991; Lee et al. 2002). In a survey of coarse and fine woody debris across the United States, Woodall and Liknes (2008) found that production and decay rates are correlated positively with availability of moisture, and negatively with increases in maximum temperature.

Measurements of carbon stocks in mineral soil in boreal stands averaged 62 ± 15 Mg C ha−1 (± SD) and contained 32 ± 8% of total carbon stocks (Fig. 4d). As would be expected due to slow accumulation and removal processes, mineral soil stocks did not show obvious age dependence. Because of the magnitude of these mineral soil carbon stocks, it has been suggested that carbon sequestration in boreal forest soils may be a useful strategy to mitigate climate change (Lal 2005). However, improved forest management strategies leading to enhanced carbon sequestration in soil require a better understanding of the dynamics of belowground allocation, heterotrophic respiration and carbon stabilization in the soil (Noormets et al. 2015).

Aboveground NPP

Annual aboveground NPP measured in various boreal stand types across a range of age and BA was generally higher in mixedwood and trembling aspen stands than in jack pine and black spruce stands (Fig. 5a, b). These measurements generally show lower values of AGNPP when the canopy is young (< 10 years) or old (> 90 years), with higher values in intermediate years, as exemplified in the present study where carbon stocks in live trees and saplings increased 15% p.a. in the juvenile stand, compared to 1% p.a. in the mature stand (Fig. 1a, Table 5). And for a given BA, mixedwood and aspen canopies typically had higher AGNPP than coniferous canopies. In addition, with increasing BA, AGNPP in mixedwood and aspen canopies were reduced, but in coniferous canopies AGNPP tended to be more constant across the observed BA range. Boreal forest stand management strategies may be adapted to improve carbon sequestration by taking into account these trends.

Fig. 5
figure 5

Aboveground NPP in various boreal forest stand types, by age and by basal area (MW mixedwood, TA trembling aspen, JP jack pine, BS black spruce: plus sign Payne, times sign Martin, asterisk Howard, open circle Gower, open square Laganiere Ontario, closed square Laganiere Quebec, open diamond Morrison (1993), closed diamond Morrison (2001), open triangle Bond-Lamberty)

The age-related decline in annual AGNPP in the present study is consistent with other studies that have shown that while older forests have the higher C stocks, younger forests have higher C absorption capacities (Bond-Lamberty et al. 2004; He et al. 2012). Because annual NPP is the sum of above- and belowground values, the estimates given (Table 6) are lower than the total carbon sequestration in the stand. However, published measurements show that aboveground NPP is typically about two-thirds of total annual NPP (e.g., Bond-Lamberty et al. 2004; Li et al. 2002). For example, in post-fire regenerating black spruce forest, which was mixedwood in composition up to age 50, the aboveground NPP exceeded 78% of total NPP for stands aged less than 37 years, and averaged 60% for older stands (Bond-Lamberty et al. 2004). Other field measurements showed that belowground NPP is less than half of total NPP in jack pine, black spruce and trembling aspen stands (Li et al. 2002). So aboveground annual NPP estimates provide a useful basis for estimating total annual NPP for the stand.

Implications for boreal forest stand management

Carbon stocks in mature boreal mixedwood and trembling aspen stands (Fig. 2a; mean 213 and 239 Mg C ha−1) were found to be at or above the average value for those found in various managed boreal forest ecosystems (193 Mg C ha−1, Kurz et al. 2013). And carbon stocks in vegetation in mixedwood and trembling aspen stands were typically greater than those found in jack pine and black spruce stands (Fig. 3a). Similarly, carbon stocks in the forest floor organic layer and detritus were typically greater in boreal mixedwood stands, compared to those found in jack pine and black spruce stands (Fig. 4a, c). Annual AGNPP was also greater in mixedwood and trembling aspen stands than in jack pine and black spruce at all stages of stand development (Fig. 5a). Overall, measurements from this and other studies demonstrate the potential value of boreal mixedwood and trembling aspen stands in helping achieve climate change mitigation through enhanced carbon sequestration and provide data needed for quantifying climate change mitigation potential with a view to modifying forest management practices (Lemprière et al. 2013). Investigation of stand thinning as a way to increase carbon sequestration, by taking advantage of the greater NPP in juvenile canopies without large reductions in carbon stocks, may be a worthwhile opportunity for future study.

Conclusion

Published measurements and those from the present investigation show that total carbon stocks, and those in the vegetation and forest floor organic layer, were typically higher in boreal mixedwood and trembling aspen stands than in coniferous canopies. Of these stocks, the average fraction present in the forest floor organic layer and mineral soil was 20 and 32%, and these components provide relatively stable carbon storage over the stand lifetime, including harvesting. In addition, measurements of annual AGNPP were typically higher in boreal mixedwood and trembling aspen canopies as compared to boreal coniferous canopies. These measurements of AGNPP also reveal higher values in intermediate age stands, as compared those that were recently harvested or mature. The value of boreal mixedwood forest as an above-average ecosystem for carbon sequestration in climate change mitigation is well demonstrated by results from this and other Canadian boreal forest ecosystem studies. Carbon stocks present in boreal mixedwood stands are determined by the particular tree species present and their proportions, as well as the local soil conditions and climate; thus boreal mixedwood carbon stocks over a widespread geographical area can only be estimated from individual stand measurements. However, measurements from the present study, together with other Canadian boreal ecosystem carbon sequestration measurements, can aid decision-making for management practices in Canada’s boreal mixedwood forests, e.g., in relation to restocking and timing of harvesting to optimise carbon sequestration and climate-change mitigation.