Next Article in Journal
Modeling Crown Biomass for Four Pine Species in China
Previous Article in Journal
Building Resilience into Sitka Spruce (Picea sitchensis (Bong.) Carr.) Forests in Scotland in Response to the Threat of Climate Change
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stand Dynamics, Humus Type and Water Balance Explain Aspen Long Term Productivity across Canada

1
Chaire Industrielle CRSNG-UQAT-UQAM en Aménagement Forestier Durable, Université du Québec en Abitibi-Témiscamingue (UQAT), 445 boul. de l'Université, Rouyn-Noranda, QC J9X-5E4, Canada
2
Institut des Sciences de la Forêt Tempérée, Université du Québec en Outaouais, 58 rue Principale, Ripon, QC J0V-1V0, Canada
3
Centre d'étude de la Forêt, Université du Québec à Montréal, Case postale 8888, Succursale Centre-Ville, Montréal QC H3C-3P8, Canada
*
Author to whom correspondence should be addressed.
Forests 2015, 6(2), 416-432; https://doi.org/10.3390/f6020416
Submission received: 6 September 2014 / Accepted: 2 February 2015 / Published: 9 February 2015

Abstract

:
This study examined the relative importance of soil, stand development and climate hypotheses in driving productivity for a species that is widely distributed in North America. Inventory plots, 3548 of such, either dominated by aspen or made up of species mixture of which aspen occurs in dominant canopy position were sampled along a longitudinal gradient from Quebec to British Columbia. Site index (SI), was used as a measure of productivity, and soil, climate and stand attributes were correlated with site index in order to determine their effects on productivity. Results show a decline in productivity with high moisture deficit. Soil humus correlates significantly with SI but does not sufficiently capture differential rates of litter deposition and decomposition effects over the long-term. Consequently, aspen composition, stand ageing, and stand structural changes dominate variability in productivity. Within the context where deciduous cover has being increasing, there are implications for forest productivity.

1. Introduction

The inherent nutritive potential of a soil plays a major role in the nutrition and growth of trees [1] hereafter described as the soil hypothesis. Tree growth is affected by biotic factors such as the activity of soil microbes [2], which is influenced by litter quality and therefore changes in stand structure and composition [3]. Stand structure and compositional dynamics represent developmental changes that are driven by disturbances, tolerance, and competition [4,5,6]. Consequently, areas with long fire cycles tend to develop complex stand structures and are characterized by higher species richness [7], thus stand development hypothesis. Also, climate directly impacts the rate at which a tree acquires and assimilates nutrients [8], and therefore climate hypothesis. The contribution of each factor to the overall health and growth of a tree vary both spatially and temporally.
The period 1983–2012 was the warmest 30-year period in the last 1400 years for the northern hemisphere. This is because atmospheric CO2 has increased by 40% since pre-industrial times [9]. Warming is expected to continue beyond 2100 and northern latitudes are predicted to warm more rapidly and also the contrast between drier and wetter regions is expected to be more pronounced [9]. These changes have consequences for forest productivity as lengthening of growing seasons leads to increased CO2 uptake. Conversely, higher soil temperature and increased biomass return to the forest floor would lead to an increase in decomposition activities and therefore higher respiratory losses [8]. Effective resource monitoring is essential in better understanding how forest ecosystems are responding to global changes. Developing tools capable of tracking these changes is therefore paramount to achieving this goal.
Recent studies report that climate warming will result in an expansion of deciduous forests and reduction in conifer cover [10]. A climate-driven increase in boreal fire activity is also expected to promote the expansion of aspen stands [11] in the future. Aspen (Populus tremuloides Michaux) is the most widely distributed broadleaf species in North America [12] and therefore ecologically important. The wood is light and soft and used for lumber, oriented strand-board, and pulp [13,14]. Across Canada, the species may occur as single species stand, or in mixed species stands, commonly in combination with Balsam fir (Abies balsamea (L.) Mill.), Black spruce (Picea mariana Mill. BSP), White spruce (P. glauca Moench Voss), Yellow or White birch (Betula alleghaniensis Britt. and Betula papyrifera Marsh., respectively). In most of interior and western Canada moisture balance was reported to be the major factor influencing aspen growth and mortality [15,16,17], while warming was identified as the major driver of aspen growth in aspen dominated stands in eastern Canada [18]. However, within mixed species stands, direct effects of warming on productivity were minimal [19]. Aspen growth could also be significantly impacted by forest tent caterpillar (Malacosoma disstria Hbn.) defoliation to the extent that the direct effects of climate could become difficult to detect at shorter time scales [20]. In better tracking the response of aspen forests to environmental changes over the longer term, models need to reflect a wide range of environmental conditions that are relevant for the species growth and productivity.
Site index (SI), defined as mean height of trees at a reference age (normally 50 years), is a measure of long-term productivity of forest stands. It is widely applied within even aged single species stands [21], but also applicable to mixed species forest stands [19]. The direct relationship between stand height growth and volume growth, and the non-dependence of height growth on stocking [21], makes SI a useful and widely applied measure of forest productivity. The goals of this work were to apply the concept of SI in measuring productivity within aspen forests and identify the major factors driving productivity along climatic gradients and also along gradients of varying stand structure, composition, age, and soil attributes. This study sought to (1) determine the relative importance of the soil, developmental and climate hypotheses in driving aspen productivity across Canada; (2) parameterize a productivity model that could be used in monitoring the productivity of aspen forests across Canada; based on identified factors and (3) identify potential areas of emphasis owing to peculiar sensitivity to the identified drivers.
Aspen is a fast growing pioneer, intolerant hardwood species with preference for moist-textured soils [12]. It is therefore hypothesized that both warming and moisture balance would be a significant factor controlling aspen growth. Being widely distributed, aspen is found on sites where substrate nutritional attributes vary geographically [22]. It is therefore expected that aspen growth would vary geographically with soil nutritional status explaining the majority of this spatial variability. Given that aspen mixed species forest stands are widespread across Canada [6] and that productivity is influenced by competition and disturbances [19], it is hypothesized that developmental changes would dominate the dynamics of site productivity within mixed species stands. However, the relative contributions of each of these three factors to the overall health and productivity of aspen forests, is unknown.

2. Methodology

2.1. Study Area

The study area is a longitudinal transect extending from Québec (from longitude 64° W) to British Columbia (to longitude 125° W) and from latitude 42° N in Ontario to 60° N in British Columbia (BC), Figure 1.
Figure 1. Map showing the distribution of plots used in this study. Plot densities are higher in the east and the area of individual plots tends to be higher in the west. The deeper the shade of blue, the higher the amount of precipitation (1981–2010 climate normal).
Figure 1. Map showing the distribution of plots used in this study. Plot densities are higher in the east and the area of individual plots tends to be higher in the west. The deeper the shade of blue, the higher the amount of precipitation (1981–2010 climate normal).
Forests 06 00416 g001
The transect is located within three boreal ecological zones; plots in Quebec, Ontario and eastern Manitoba were within the boreal shield ecozone, those in mid-western Manitoba, Saskatchewan and Alberta were within the boreal plains while plots in BC were within the boreal cordillera ecological zone. The Boreal cordillera is dominated by mountains and valleys with wide lowlands and relatively sparse populations, mainly the Yukon population. The Boreal plains is the low lying interior of Canada, agriculture is an important activity here however, forestry remains the main industry with about 80% of the area forested. The boreal shield derives its name from the fact that geographically, at this point the boreal zone meets the Canadian shield rock.
The study area is comprised of 32-ecological regions (areas characterized by a particular distribution of vegetation types) and 126-ecological districts (land units characterized by particular relief, soil deposits, and drainage regimes). Drainage is largely imperfect and the landscape is generally gently sloped (<10%). Total annual precipitation varies along the study transect from about 600 mm in the west to about 1300 mm in the eastern part of the transect (Figure 1). Annual snow fall is about 160 mm in the boreal cordillera, 90 mm in the plains with an eastward increasing trend within the boreal shield to about 400 mm. Annual mean temperature varies from about −1 °C in the western part of the transect to 2.5 °C in the east.

2.2. Tree Measurement Data

In this study, inventory data was obtained from the provincial jurisdictions in all six provinces along the study transect. Quebec’s Ministère des Ressources Naturelles (MRN) started a decadal province-wide measurement program in 1970, where circular plots (400 m2) were positioned in forested lands through a stratified sampling scheme. Data used in this study comes from its third inventory program (1992–2002). In Ontario, the Ontario Forest Research Institute has a forest Growth and Yield Program where 6400 m2 main plots are sampled from managed and natural forests within which are nested three 400 m2 growth plots and nine 25 m2 shrub/regeneration plots. Data from the growth plots were used in this study. In Manitoba, data was collected from 500 plots by Louisiana Pacific Canada Ltd., Swan Valley Forest Resources Division. In Saskatchewan, trees were sampled within plots of varying sizes from 100–1000 m2 depending on the stand density and height. Historically, data was collected by the Saskatchewan Ministry of Environment but responsibility was transferred to licensed holders. Data from different sources were assembled by Weyehaeuser into one database for Saskatchewan. In Alberta, 200–3400 m2 plots were sampled and re-measured every five to ten years from 1960 by the Ministry of Environment and Sustainable Resource Development. Finally, data from British Columbia was obtained by the Ministry of Forests, Lands and Natural Resource Operations from plots monitored every five to fifteen years from 1920. The area of plots varied from 400–1000 m2, depending on stand density.
From a total of 8560 inventory plots within the transect area, 3548 plots had aspen in dominant or co-dominant canopy positions. Plots were identified as either pure aspen stands (19% of such) or mixture species stands (81% of such). A pure aspen plot was defined as having at least 75% of total plot basal area contributed by aspen. Plot density varied along the transect (Figure S1). Consequently, the number of trees selected for top height estimation varied by plot. By definition, top height refers to the 100 largest trees per hectare, the four largest trees per 400 m2, or the 30 largest trees per 3000 m2. For each plot, data was collected on trees with a diameter at breast height (DBH) > 9 cm. For age determination, trees were cored at breast height (1.3 m) or stump height (0.30 m) corresponding to breast height age or total age, respectively. Since the study focused on aspen, in stands where aspen proportion was relatively low, it was expected that plot mean of aspen age would be different from mean age estimated from all trees.

2.3. Factors Potentially Linked to Productivity

2.3.1. Climate Data

Climate data was generated from the BioSIM software [23] using the “climatic annual” model. This model computes annual climate variables from a specified number of weather stations nearest to the point of interest and corrects for differences in elevation and aspect of the weather stations to the point of interest. It also accounts for the distance of the weather stations to the point of interest by assigning greater interpolation weights to nearer stations. Two climate normals, 1961–1990 and 1981–2010 were used in generating climate data.
Climate variables included annual means of lowest daily minimum temperature, daily minimum temperature, daily mean temperature, daily maximum temperature, highest daily maximum temperature, annual total precipitation (mm), annual snowfall (mm), annual mean of daily mean dew point temperature (°C), annual mean of daily mean relative humidity (mm), mean wind speed, annual growing degree-day above 5 °C, frost days (number of days in the year when the daily minimum temperature is <0), frost free days (longest uninterrupted period without frost in the year), growing season days (the period between the last three consecutive days with frost in the spring and the first three consecutive days with frost in the fall), annual total potential evapotranspiration (mm), and aridity index, which is accumulated monthly water deficit (Thornwaite potential evapotranspiration-monthly precipitation, in mm; zero if negative). Other variables included annual mean of daily vapor pressure deficit (kpa), annual total radiation (MJ/m2), number of days with precipitation >1 mm, consecutive days without precipitation (precipitation <1 mm) and climate moisture index (total annual precipitation–potential annual evapotranspiration) [15].

2.3.2. Soil Data

As part of the inventory program, data on soil texture, drainage class, slope class, surface deposit type, humus type, and depth of organic layer were collected. In dealing with differences in classification systems across the provinces, a particular system was selected and others re-classified along the selected system. For the purposes of this study, the Saskatchewan drainage classification system was adopted; classes 1–7 denote respectively; very rapid, rapid, well, moderately well, imperfect, poor, and very poor drainage classes. Slope percentages were grouped into 10 classes where 1 is level (<0.6%) and 10 is very steep slope (>100%). Humus types were also regrouped into five classes; Mull (MU), Moder (MD), Mor (MO), Anmoor (AN), and peatland humus (OT). The peatland humus type is a combination of the organic humus type and the peaty humus types. Both soil humus and texture are indicative of soil nutritional status [3,24].

2.3.3. Stand Attributes

Sampled stands vary in stand age, diameter structure and species composition. Diameter structure was measured with the Shannon diversity index computed in Equation (1) [25]:
S h = i = 1 s p i In ( p i )
where p i is the basal area proportion of the diameter class i relative to the total stand basal area, and s is the number of diameter classes.

2.4. Modeling Productivity

2.4.1. Model Selection

Across the provinces covered by the study transect, there are several site index models potentially applicable in this study. A literature search allowed consideration of five models. A model is considered because it is used in growth and yield estimation within a particular province or because of its potentially wide applicability. Considered were i) the model by Pothier and Savard [26] currently applied in growth and yield estimation in Quebec, ii) the model by Plonski [27,28] widely applied in Ontario, iii) the model by Huang et al. [29] used in growth and yield modelling (GYPSY model) in Alberta and validated for parts of central and western Canada [30], iv) the model by Nigh et al. [31], otherwise called the site tools software used widely in British Columbia and recently validated for central and western Canada [30], and finally, 5th model) the model by Garcia [30], which was calibrated using data spanning from Manitoba to British Columbia and thus valid for central and western Canada.
Models were first applied to data from the particular region for which they were calibrated. The site index values produced in these cases were called reference SIs. For a given region with a reference SI, the other models were applied to obtain four alternative SIs per region. The reference and alternative values were used in computing root mean square error for each alternative model and by region. The model selected was the one which overall, produced the minimum root mean square error. As the model of Nigh et al. [31] produced the minimum root mean square error, it was selected for site index estimation.

2.4.2. Estimating Site Index

The model of Nigh et al. [31] is as follows:
H = 1.3 + ( S I 1.3 ) [ 1 + e a 0 + a 1 * I n ( 49.5 ) + a 2 * I n ( S I 1.3 ) 1 + e a 0 + a 1 * I n ( A B H 0.5 ) + a 2 * I n ( S I 1.3 ) ]
where H is plot dominant height, SI is site index at reference age of 50-years, ABH is breast height age and a0, a1, a2 are species specific model parameters. The height of aspen trees was available in the database from which plot dominant height was computed.
Existing site index models, especially those that produce anamorphic SIs, are not particularly suitable for juvenile forest stands [32]. Previous work addressed this problem either by excluding juvenile stands from site index modelling or using growth intercept approach to estimate site index [33]. Growth intercept approach computes SI from the early (<50 years) average height growth. Here, we applied the growth intercept method in obtaining SI values for juvenile (<30 years, Supplementary information section 2.0) forest stands. Stem analysis data from Anyomi et al. [34] was used for this purpose, and the resultant SI was denoted as the periodic growth index ( S i H ). Growth intercepts were computed for each breast height age and correlated with S i H . Results showed strong correlation between growth intercept and S i H (R2 = 0.81), Equation (3).
S i H = 1.3 + e b 1 G I b 2
where S i H is periodic growth index, G I is growth intercept (cm/year) and b1, b2 are parameter values. The parameter values of Equation (3) can be found in Table S1. Equation (3) was therefore used in estimating SI for juvenile stands.

2.4.3. Modeling Approach

Soil variables were first accounted for in the model (soil hypothesis) and then other variables (climate and stand developmental hypothesis) were introduced stepwisely. An approach in which the mean value of site index ( S I ¯ ) across the study transect is modified by retained explanatory variables ( x i ) was applied [34] in modeling productivity, Equations (4) and (5).
S I = S I ¯ i = 1 n f i ( x i )
f i ( x i ) = 1 + β l . x i ( x i x ¯ i x ¯ i ) + β q . x i ( x i x ¯ i x ¯ i ) 2
where β l . X i and β q . X i represent the linear and quadratic effects of the explanatory variable x i on SI. Categorical variables were converted to dummy variables so that they could be used in the NLMIXED procedure in SAS (see Supplementary information section 3.0). The product of n modifiers has a value close to unity when the variable x i is equal to the average x ¯ i observed value, and increases or decreases when moving away from the average.
A variable was retained in the model if it significantly (p < 0.05) related with SI and contributed substantially (>1%) to the total explained variance (or in reducing RMSE). Inclusion or exclusion of a modifier was determined with a likelihood-ratio test [35].
In verifying if some regional-level factors contribute to stand-level productivity, which were not already included in the model, a random effect variable ( s ) was introduced into Equation (4) and tested at defined landscape scales, Equation (6).
S I = ( S I ¯ + s ) i = 1 n f i ( x i )
The parameters of Equation (6) were estimated with the NLMIXED procedure to, offer the potential to capture both fixed and random effects that may occur at different spatial scales, in this case stand and regional-levels respectively. Ecological district and ecological region were the regional-level categories considered. These models were rated against the model that included only fixed effects of local variables with a likelihood-ratio test.
The sensitivity of aspen forests to climatic and non-climatic (stand attributes) growth driving factors was measured using the coefficient of determination (R2) and the Akaike information criterion (AIC). Annual mean temperature was grouped into five classes; class 1 (<−1 °C), class 2 (−1 to 0 °C), class 3 (0 to 1 °C), class 4 (1 to 2 °C), class 5 (>2 °C) and similarly total annual precipitation was grouped into classes 1 (500–700 mm), 2 (700–900 mm), 3 (900–1100 mm), 4 (1100–1300 mm) and 5 (1300–1500 mm).

3. Results

3.1. Variability in SI

SI significantly correlated with stand variables (Figure 2A–D). SI decreased with stand age (AIC = 19,048.5; r = −0.51) as well as with an increase in proportion of spruce basal area(AIC = 19,275.8; r = −0.28). Aspen dominated stands exhibited higher SI which decreased with a decrease in proportion of aspen basal area (AIC = 19,355.1; r = 0.24). There was also a positive effect of stand horizontal structure measured with the Shannon index (AIC = 19,535.1; r = 0.05) such that SI slightly increased with structural complexity.
Figure 2. Variability in productivity with stand attributes (A) stand age, basal area proportions of (B) spruce and (C) aspen, (D) Shannon index and spatial elements; (E) latitude and (F) longitude and (G) humus type (H) aridity index (mm). Red line denotes linear or quadratic fits and grey lines refer to the number of plots under each x-category. Error bars denote standard error.
Figure 2. Variability in productivity with stand attributes (A) stand age, basal area proportions of (B) spruce and (C) aspen, (D) Shannon index and spatial elements; (E) latitude and (F) longitude and (G) humus type (H) aridity index (mm). Red line denotes linear or quadratic fits and grey lines refer to the number of plots under each x-category. Error bars denote standard error.
Forests 06 00416 g002aForests 06 00416 g002b
SI also slightly varied by latitude (ns) and longitude (Figure 2E–F) (AIC = 20107.1; r = 0.05). Aspen stands located from longitude −110° W to −120° W were the most productive with mean SIs of over 25 m. SI slightly increased from longitude −108° W eastwards and also slightly increased with latitude with the most productive sites located above latitude 55° N. A significant positive correlation with mull humus type and a negative correlation with moder were also observed (AIC = 15,996; R2 = 6%). There was also a negative effect of moisture balance on SI (AIC = 19,934.3; r = −0.22), particularly when moisture deficit was beyond 60 mm (Figure 2G–H).

3.2. Major Drivers of SI

Of the climate variables considered in this study, only relative humidity and potential evapotranspiration were not significantly correlated with SI and aspect was the only site variable not significantly related to SI. Soil texture was a major explanatory variable (R2 = 8%), but there were 15-textural classes and in ensuring parsimony with respect to model parameters, soil humus type was rather retained. Soil effects were first accounted for in the model, humus was therefore the first explanatory variable retained, explaining 6% of variability in SI. Indeed, results suggest that aspen growth is better on mull humus type: mean SI was higher on mull compared to the other humus types. Higher variability in SI was observed on Anmoor humus type (Figure 2G). Stand age was the second most correlated variable after accounting for soil nutritional status and it explained 29% of the residual variance. Stand diameter structure measured with Shannon index and aspen basal area proportion subsequently explained 11% and 4% respectively of residual variance. Finally, aridity index explained 3% of the remaining variability in SI after controlling for soil and stand features. Even though other climatic variables were additionally significant, their contributions to total explained variance were marginal (<1%), hence they were not retained in the model. The selected variables, Equation (7), explained 53% of variability in aspen SI.
S I = ( S I ¯ + i = 1 n β i z H u m , i )    f A g e f S h a n f G a s p e n f A r i d i t y
where S I is predicted site index, S ¯ is the mean site index, i = 1 n β i z H u m , i accounts for humus effect on SI, f A g e , f S h a n , f G a s p e n , f A r i d i t y are respectively modifiers for stand age, Shannon index, aspen basal area proportion and aridity index.
Ecological district ( s D ) as a random parameter, contributed 12% to the explained variability in stand-level SI, raising the total explained variability in SI to 65%. However, much of the variability captured by this parameter is related to plot density at that scale (Figure S2).
S I = ( S I ¯ + i = 1 n β i z H u m , i + s D ) f A g e f S h a n f G a s p e n f A r i d i t y
where s D is random effect parameter that varies by ecological district (Table S2), parameter values of Equation (8) can be found in Table 1.
Table 1. Parameter values of model [8].
Table 1. Parameter values of model [8].
VariableParameterEstimateStandard ErrorP-value
Soil humus type β Moder −2.2890.484<0.0001
β Mor −2.2630.473<0.0001
β Mull −1.7360.4910.000
β Anmoor −2.3951.5810.132
β Peantland −2.2280.6400.000
Aspen basal area proportion β l.Gaspen 0.0790.006<0.0001
Stand age β l.Age −0.5120.011<0.0001
β q.Age 0.2940.012<0.0001
Shannon index β l.Shan 0.4370.019<0.0001
β q.Shan 0.2750.036<0.0001
Aridity index β l.Aridity −0.0310.0090.012
β q.Aridity 0.0140.0070.002
Mean Site Index S I ¯ 18.6280.468<0.0001

3.3. Sensitivity to Climate

Sensitivity to aridity index (climatic variable retained in the model, Equations (7) and (8) as illustrated in Figure 3A,C reveal a pattern similar to the sensitivity to stand dynamics (effects of the changes in stand age, structure and composition; Figure 3B,D). Sites with annual mean temperature of less than −1 °C (Classes 1 & 2) were the most sensitive along the temperature gradient while areas with annual total precipitation ranging from 900 mm to 1300 mm (Classes 3 & 4) were relatively less sensitive to the effects of both climate and stand dynamics (Figure 3C,D).
Figure 3. Changes in the strength of productivity response to aridity index (A,C) and stand dynamics (the combined effects of stand age, structure and composition) (B,D); measured with co-efficient of determination (A,B) and Akaike information criterion (C,D). The solid red line depicts the variability in the strength of the relationships along the temperature gradient, while the broken dark blue line reveals changes in the strength of the relationship along the precipitation gradient.
Figure 3. Changes in the strength of productivity response to aridity index (A,C) and stand dynamics (the combined effects of stand age, structure and composition) (B,D); measured with co-efficient of determination (A,B) and Akaike information criterion (C,D). The solid red line depicts the variability in the strength of the relationships along the temperature gradient, while the broken dark blue line reveals changes in the strength of the relationship along the precipitation gradient.
Forests 06 00416 g003

4. Discussion

4.1. Stand Development Explains Long-Term Productivity

Consistent with expectations, results suggest that humus type is an important driver of aspen productivity. Humus forms vary spatially, mainly by geology and topography, and temporally, by climate and litter type. The nature of sub-soils influences litter quality and also the accessibility to earthworms, such that more acidic soils tend to promote mor humus formation while mull humus develops in more alkaline environments [24]. Mull humus type is transformed organic material, mixed with mineral soil through activities of soil fauna and since it leads to more complete decomposition, nutrient availability is high and thereby supports greater productivity [3]. Mors on the other hand are primarily products of fungal decomposition and because of incomplete decomposition, nutrients are fairly immobilized. Higher average SI observed on mull humus type is therefore due to its nutritional potential. Understanding the spatial repartitioning of humus forms and their changes through time within the boreal forest, would greatly assist in the long-term management of aspen productivity. That said, analysis also revealed that humus forms do not entirely reflect litter decomposition regime along the studied transect.
Results show that stand age, changes in species composition, and stand structure, correlate significantly with SI. Stand ageing, composition, and structural changes are dynamics that characterize stand development, and ultimately stand succession [36]. In this study, SI declined over time; we observed mean SI of 21 m (standard error = 0.46) at mean stand age of 30 years, a mean of 17 m (0.38) at a stand age of 60 years and a mean SI of 14 m (1.40) at a stand age of 90 years. Three processes can explain the age-related decline in productivity, (1) physiological constraint to CO2 sequestration [37]; (2) accumulation of organic material that results in soil cooling and nutrient immobilization [38]; and (3) stand age as a driver of leaf litter quality decline [39].
Because a significant proportion of the studied plots (81%) has aspen trees in a mixture with other tree species, in line with our hypothesis, SI also vary by the composition of aspen. A mean site index of 20 m (0.78) was observed in pure aspen stands (aspen basal area proportion >75%) compared to a mean of 18 m (0.90) in mixed species stands. Compositional effects reflect the differential rates of litter deposition vis a vis rates of decomposition because hardwood leaves decompose faster than conifer leaves [3]. Changes in the species composition of the dominant cohort [40] could also directly impact SI since tree species inherently vary in their growth potential. However, given that the effect of stand structure on SI is marginal (Figure 2D), it implies that compositional changes reflect more accurately the decomposition regime driven by changes in litter quality than by the direct impact of dominant cohort displacement on SI.
Finally and also consistent with research hypothesis, results suggest moisture deficit constraints on aspen SI, especially with deficits higher than 60 mm. By reference to an SI of 19 m, aspen height dropped by 5 m with a doubling in moisture deficit from 60–120 mm (Table 1; Figure 2). This effect is equivalent to the loss of 13 years in height growth. This observation adds to a growing body of evidence on the role of moisture in aspen growth and productivity [15,16,17], especially considering that aspen does not tolerate imperfectly drained soils and have only been observed on moderate to well drained soils. Yet, within the ensemble of environmental factors influencing aspen long-term productivity, the contribution of this variable was only 3%.

4.2. Aspen Forests in Cooler, Drier Environments are Most Sensitive to Productivity Drivers

Results suggest that, the productivity of aspen forests located in areas with mean annual temperature less than −1 °C was the most sensitive to growth drivers while relatively warmer regions of aspen stands (>2 °C) are less responsive to changes in growth drivers (Figure 3). Consequently, most of northern aspen stands (>latitude 54° N) of BC, Alberta and Saskatchewan were relatively more responsive to growth drivers when compared to those of Manitoba, Ontario and Quebec (Figure 4A). Along the transect from east to west, two major areas were most responsive to growth drivers, north east of Alberta and western Ontario, because of the relatively low precipitation amounts these areas receive annually (<900 mm). These observations were true, when only moisture balance was the explanatory variable in the model as well as when both moisture balance and stand attributes were considered. Overall, stands in western Quebec and eastern Ontario were most resistant (less responsive) to the impacts of moisture and stand dynamics.
Figure 4. (A) Northern lying (deeper shades) aspen forests are more sensitive to productivity drivers relative to more southern stands (lighter shades) along a north south temperature gradient; (B) Areas that are more responsive to growth drivers along an east-west precipitation gradient. The darker the shade, the higher the sensitivity of the stand.
Figure 4. (A) Northern lying (deeper shades) aspen forests are more sensitive to productivity drivers relative to more southern stands (lighter shades) along a north south temperature gradient; (B) Areas that are more responsive to growth drivers along an east-west precipitation gradient. The darker the shade, the higher the sensitivity of the stand.
Forests 06 00416 g004

4.3. Implications within the Context of Global Changes

Results from this study show that aspen composition has a positive influence on productivity and that stand dynamics have greater impact on forest productivity than the direct effects of climate. While stand dynamics may be directly influenced by climate [10], disturbance events are increasingly common within the circumboreal region [41,42] and substantially affect stand attributes [36,43]; this study shows that such dynamics control productivity of aspen forests across Canada. A recent study [11] suggests that the fire regime in the last millennia exceeds the limit observed within the last 10,000 years. Such intense fires promote shifts to an increasingly deciduous-dominated landscape [11,36]. Aspen abundance is also favored by logging [44], and considering that the demand for aspen wood is expanding, primarily because it is a relatively low-cost, easy-to-use wood, an increase in logging would promote aspen expansion. With an increase in aspen abundance, this work reports that aspen productivity would substantially increase. Considering that plots cover a wide spatial extent and data was sampled by different provincial officials, differences in plot selection and measurement protocols could potentially affect our results.

5. Conclusions

Analyses of inventory data along an east-west transect show that variability in aspen long-term growth is mainly driven by species compositional changes and stand ageing effects. While soil humus plays a major role in long-term growth, it does not entirely reflect litter deposition and decomposition regime effects on productivity. Extreme dry events have major negative consequences for aspen growth but, within the ensemble of growth driving factors, its effects are relatively weak. It is concluded that stand development drives aspen productivity along an east-west precipitation gradient. An empirical model has been calibrated which could aid in monitoring forest site productivity of aspen forests in Canada.

Acknowledgments

The authors are grateful to all provincial offices for making available data used in this study and Pierre Nlungu Kweta for assistance in procuring these data. We are also grateful to Pierre Mérian for the climate data. This project was financed by MITACS accelerate programme IT02897.

Author Contributions

Kenneth A. Anyomi & Yves Bergeron worked on project proposal also reviewed by François Lorenzetti. Alain Leduc supervised acquisition of inventory data, François Lorenzetti supervised acquisition of climate data. All authors were involved in data analysis, and the discussion of results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kimmins, J.P. Importance of soil and role of ecosystem disturbance for sustained productivity of cool temperate and boreal forests. Soil Sci. Soc. Am. J. 1996, 60, 1643–1654. [Google Scholar] [CrossRef]
  2. Van der Heijden, M.G.A.; Bardgett, R.D.; van Straalen, N.M. The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 2008, 11, 296–310. [Google Scholar]
  3. Prescott, C.E.; Zabek, L.M.; Staley, C.L.; Kabzems, R. Decomposition of broadleaf and needle litter in forests of British Columbia: Influences of litter type, forest type, and litter mixtures. Can. J. For. Res. 2000, 30, 1742–1750. [Google Scholar] [CrossRef]
  4. Hart, S.A.; Chen, H.Y.H. Fire, logging and overstory affect understory abundance, diversity and composition in boreal forest. Ecol. Monogr. 2008, 78, 123–140. [Google Scholar] [CrossRef]
  5. Rohne, B.; Bigler, C.; Wunder, J.; Brang, P.; Bugmann, H. Fifty years of natural succession in Swiss forest reserves: Changes in stand structure and mortality rates of oak and beech. J. Veg. Sci. 2012, 23, 892–905. [Google Scholar] [CrossRef]
  6. St.-Clair, S.B.; Cavard, X.; Bergeron, Y. The role of facilitation and competition in the development and resilience of aspen forests. For. Ecol. Manag. 2013, 299, 91–99. [Google Scholar] [CrossRef]
  7. Bouchard, M.; Pothier, D.; Gauthier, S. Fire return intervals and tree species development in the North Shore region of eastern Quebec. Can. J. For. Res. 2008, 38, 1621–1633. [Google Scholar] [CrossRef]
  8. Grant, R.F.; Margolis, H.A.; Barr, A.G.; Black, T.A.; Dunn, A.L.; Bernier, P.Y.; Bergeron, O. Changes in net ecosystem productivity of boreal black spruce stands in response to changes in temperature at diurnal and season time scales. Tree Physiol. 2008, 29, 1–17. [Google Scholar] [CrossRef] [PubMed]
  9. IPCC. Summary for Policymakers. In Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  10. Tang, G.; Beckage, B.; Smith, B. The potential transient dynamics of forests in New England under historical and projected future climate change. Clim. Chang. 2012, 114, 357–377. [Google Scholar] [CrossRef]
  11. Kelly, R.; Chipman, M.L.; Higuera, P.E.; Stefanova, I.; Brubaker, L.B.; Hu, F.S. Recent burning of boreal forests exceeds fire regime limits of the past 10,000 years. PNAS 2013. [Google Scholar] [CrossRef]
  12. Peterson, E.B.; Peterson, N.M. Ecology, Management and Use of Aspen and Balsam Poplar in the Prairie Provinces, Canada; Special Report; Forestry Canada, Northwest Region Northern Forestry Centre: Edmonton, Alberta, Canada, 1992; Volume 1, p. 252. [Google Scholar]
  13. Ondro, W.J. Utilization and Market Potential of Poplar in Alberta; Information Report; Forestry Canada, Northern Forestry Centre: Edmonton, Alberta, Canada, 1989; p. 72. [Google Scholar]
  14. Christersson, L. Wood production potential in poplar plantations in Sweden. Biomass Bioenergy 2010, 34, 1289–1299. [Google Scholar] [CrossRef]
  15. Hogg, E.H. Climate and the southern limit of the western Canadian boreal forest. Can. J. For. Res. 1994, 24, 1835–1845. [Google Scholar] [CrossRef]
  16. Hogg, E.H.; Barr, A.G.; Black, T.A. A simple soil moisture index for representing multi-year drought impacts on aspen productivity in the western Canadian interior. Agric. For. Meteorol. 2013, 178–179, 173–182. [Google Scholar] [CrossRef]
  17. Worrall, J.J.; Gerald, E.R.; Hamann, A.; Hogg, E.H.; Marchetti, S.B.; Michaelian, M.; Gray, L.K. Recent declines of Populus tremuloides in North America linked to climate. For. Ecol. Manag. 2013, 299, 35–51. [Google Scholar] [CrossRef]
  18. Lapointe-Garant, M.-P.; Huang, J.-G.; Gea-izquierdo, G.; Raulier, F.; Bernier, P.Y.; Berninger, F. Use of tree rings to study the effect of climate change on trembling aspen in Québec. Glob. Chang. Biol. 2010, 16, 2039–2051. [Google Scholar] [CrossRef]
  19. Anyomi, K.A.; Raulier, F.; Bergeron, Y.; Mailly, D. The predominance of stand composition and structure over direct climatic and site effects in explaining aspen (Populus tremuloides Michaux) site index within boreal and temperate forests of western Quebec, Canada. For. Ecol. Manag. 2013, 302, 390–403. [Google Scholar] [CrossRef]
  20. Lorenzetti, F.; Université du Québec en Outaouais, Ripon, Canada. Unpublished work. 2015.
  21. Skovsgard, J.P.; Vanclay, J.K. Forest site productivity: A review of the evolution of dendrometric concepts for even-aged stands. Institute of chartered foresters. Forestry 2007, 81. [Google Scholar] [CrossRef]
  22. Paré, D.; Bergeron, Y.; Longpré, M.-H. Potential productivity of aspen cohorts originating from fire, harvesting and tree-fall gaps on two deposit types in northwestern Quebec. Can. J. For. Res. 2001, 31, 1067–1073. [Google Scholar] [CrossRef]
  23. Régnière, J.; Saint-Amant, R. BioSIM 9-Manuel de l’Utilisateur. Ressources Naturelles Canada, Service Canadien des Forêts; Centre de Foresterie des Laurentides: Quebec, Canada, 2008; p. 81. [Google Scholar]
  24. Ponge, J.-F.; Jabiol, B.; Gégout, J.-C. Geology and climate conditions affect more humus forms than forest canopies at large scale in temperate forests. Geoderma 2011, 162, 187–195. [Google Scholar] [CrossRef]
  25. Boucher, D.; Gauthier, S.; de Grandpré, L. Structural changes in coniferous stands along a chronosequence and a productivity gradient in the northeastern boreal forest of Québec. Ecoscience 2006, 13, 172–180. [Google Scholar] [CrossRef]
  26. Pothier, D.; Savard, F. Actualisation des Tables de Production Pour les Principales Espèces Forestières du Québec; Direction de la Recherche Forestière; Ministère des Ressources Naturelles du Québec: Quéebec, QC, Canada, 1998; p. 183. [Google Scholar]
  27. Plonski, W.L. Introduction to Plonski’s Normal Yield Tables. Available online: http://flash.lakeheadu.ca/~fluckai/pow.html (accessed on 3 February 2015).
  28. Payandeh, B. Metric site index formulae for major Canadian timber species. Dept. Fish. Environ. Can. For. Serv. 1977, 33, 37–39. [Google Scholar]
  29. Huang, S.; Meng, S.X.; Yang, Y. A Growth and Yield Projection System (GYPSY) for Natural and Post-Harvest Stands in Alberta; Forestry Division, Alberta Sustainable Resource Development; Technical Report Publications: Edmonton, Canada, 2009; p. 25. [Google Scholar]
  30. Garcia, O. Building a dynamic growth model for trembling aspen in western Canada without age data. Can. J. For. Res. 2013, 43, 256–265. [Google Scholar]
  31. Nigh, G.D.; Krestov, P.V.; Klinka, K. Trembling aspen height-age models for British Columbia. Northwest Sci. 2002, 76, 202–212. [Google Scholar]
  32. Nigh, G.D.; Love, B.W. A model for estimating juvenile height of lodgepole pine. For. Ecol. Manag. 1999, 123, 157–166. [Google Scholar] [CrossRef]
  33. Mailly, D.; Gaudreault, M. Modèles de croissance internodale pour l’épinette noire, le pin gris et le sapin baumier au Québec. For. Chron. 2005, 81, 114–124. [Google Scholar] [CrossRef]
  34. Anyomi, K.A.; Raulier, F.; Mailly, D.; Girardin, M.; Bergeron, Y. Using height growth to model local and regional response of trembling aspen (Populus tremuloides Michx.) to climate within the boreal forest of western Québec. Ecol. Model. 2012, 243, 123–132. [Google Scholar] [CrossRef]
  35. Bates, D.M.; Watts, D.G. Non-Linear Regression Analysis and Its Applications; John Wiley and Sons: New York, NY, USA, 1988; p. 359. [Google Scholar]
  36. Johnstone, J.F.; Hollingsworth, T.N.; Chapin, F.S.; Mack, M.C. Changes in fire regime break the legacy lock on developmental trajectories in Alaskan boreal forest. Global Chang. Biol. 2010, 16, 1281–1295. [Google Scholar] [CrossRef]
  37. Bond, B.J. Age-related changes in photosynthesis of woody plants. Trends Plant Sci. 2000, 5, 349–353. [Google Scholar] [CrossRef] [PubMed]
  38. Perrette, G.; Lorenzetti, F.; Moulinier, J.; Bergeron, Y. Site factors contribute to aspen decline and stand vulnerability following a forest tent caterpillar outbreak in the Canadian Clay Belt. For. Ecol. Manag. 2014, 323, 126–137. [Google Scholar] [CrossRef]
  39. Trap, J.; Hättenschwiler, S.; Gattin, I.; Aubert, M. Forest ageing: An unexpected driver of beech leaf litter quality variability in European forests with strong consequences on soil processes. For. Ecol. Manag. 2013, 302, 338–345. [Google Scholar] [CrossRef]
  40. Raulier, F.; Lambert, M.-C.; Pothier, P.; Ung, C.-H. Impact of dominant tree dynamics on site index curves. For. Ecol. Manag. 2003, 184, 65–78. [Google Scholar] [CrossRef]
  41. Stocks, B.J.; Mason, J.A.; Todd, J.B.; Bosch, E.M.; Wotton, B.M.; Amiro, B.D.; Flannigan, M.D.; Hirsch, K.G.; Logan, K.A.; Martell, D.L.; et al. Large forest fires in Canada, 1959–1997. J. Geophys. Res. 2003, 107, FFR5.1–FFR5.12. [Google Scholar]
  42. Bond-Lamberty, B.; Peckham, S.D.; Ah, D.E.; Gower, S.T. Fire as the dominant driver of central Canadian boreal forest carbon balance. Nature 2007, 450, 89–92. [Google Scholar] [CrossRef] [PubMed]
  43. Pan, Y.; Chen, J.M.; Birdsey, R.; McCullough, K.; He, L.; Deng, F. Age structure and disturbance legacy of North American forests. Biogeosciences 2011, 8, 715–732. [Google Scholar] [CrossRef]
  44. Laquerre, S.; Harvey, B.D.; Leduc, A. Spatial analysis of response of trembling aspen patches to clear-cutting in black spruce-dominated stands. For. Chron. 2011, 87, 77–85. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Anyomi, K.A.; Lorenzetti, F.; Bergeron, Y.; Leduc, A. Stand Dynamics, Humus Type and Water Balance Explain Aspen Long Term Productivity across Canada. Forests 2015, 6, 416-432. https://doi.org/10.3390/f6020416

AMA Style

Anyomi KA, Lorenzetti F, Bergeron Y, Leduc A. Stand Dynamics, Humus Type and Water Balance Explain Aspen Long Term Productivity across Canada. Forests. 2015; 6(2):416-432. https://doi.org/10.3390/f6020416

Chicago/Turabian Style

Anyomi, Kenneth A., François Lorenzetti, Yves Bergeron, and Alain Leduc. 2015. "Stand Dynamics, Humus Type and Water Balance Explain Aspen Long Term Productivity across Canada" Forests 6, no. 2: 416-432. https://doi.org/10.3390/f6020416

APA Style

Anyomi, K. A., Lorenzetti, F., Bergeron, Y., & Leduc, A. (2015). Stand Dynamics, Humus Type and Water Balance Explain Aspen Long Term Productivity across Canada. Forests, 6(2), 416-432. https://doi.org/10.3390/f6020416

Article Metrics

Back to TopTop