Next Article in Journal
What Factors Shape Spatial Distribution of Biomass in Riparian Forests? Insights from a LiDAR Survey over a Large Area
Next Article in Special Issue
Opportunities and Challenges for Intensive Silviculture in Alberta, Canada
Previous Article in Journal
Is the Seasonal Variation in Frost Resistance and Plant Performance in Four Oak Species Affected by Changing Temperatures?
Previous Article in Special Issue
Nitrogen Recovery from Enhanced Efficiency Fertilizers and Urea in Intensively Managed Black Walnut (Juglans nigra) Plantations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Trade-Offs among Release Treatments in Jack Pine Plantations: Twenty-Five Year Responses

1
Ontario Forest Research Institute, Ontario Ministry of Natural Resources and Forestry, 1235 Queen St. East, Sault Ste. Marie, ON P6A 2E5, Canada
2
Canadian Wood Fibre Centre, Canadian Forest Service, Natural Resources Canada, 1055 du P.E.P.S., P.O. Box 10380, Sainte-Foy Stn., Québec City, QC G1V 4C7, Canada
3
Faculty of Natural Resources Management, Lakehead University, Wood Science Testing Laboratory, 955 Oliver Rd., Thunder Bay, ON P7B 5E1, Canada
*
Author to whom correspondence should be addressed.
Forests 2021, 12(3), 370; https://doi.org/10.3390/f12030370
Submission received: 27 January 2021 / Revised: 8 March 2021 / Accepted: 18 March 2021 / Published: 20 March 2021

Abstract

:
We assessed 27 indicators of plant diversity, stand yield and individual crop tree responses 25 years post-treatment to determine long-term trade-offs among conifer release treatments in boreal and sub-boreal forests. This research addresses the lack of longer-term data needed by forest managers to implement more integrated vegetation management programs, supporting more informed decisions about release treatment choice. Four treatments (untreated control, motor-manual brushsaw, single aerial spray, and complete competition removal) were established at two jack pine (Pinus banksiana Lamb.) sites in Ontario, Canada. Our results suggest that plant diversity and productivity in boreal jack pine forests are significantly influenced by vegetation management treatments. Overall, release treatments did not cause a loss of diversity but benefitted stand-scale yield and individual crop tree growth, with maximum benefits occurring in more intensive release treatments. However, none of the treatments maximized all 27 indicators studied; thus, forest managers are faced with trade-offs when choosing treatments. Research on longer term effects, ideally through at least one rotation, is essential to fully understand outcomes of different vegetation management on forest diversity, stand yield, and individual crop tree responses.

1. Introduction

International agreements such as the Paris Climate Agreement, the UN Sustainable Development Goals and the Montréal Process are providing opportunities and setting challenges for forest managers worldwide. Forests can provide an array of ecosystem services while helping to offset greenhouse gas emissions but are affected by natural and anthropogenic disturbances and climate change. These disturbances can interact with climate change to modify forest composition and productivity, thus affecting the capacity of forests to sustain services such as supplying wood assortments, maintaining biodiversity, and sequestering carbon [1]. Successful stand regeneration can mitigate these effects and is a keystone of sustainable forest management. While natural regeneration plays a substantial role in forest renewal—for example, about 44% of the area harvested in the past 20 years in Canada regenerated naturally [2]—the use of plantation forestry to regenerate specific species compositions offers added advantages for supporting sustainable forest management. Planted forests cover ~291 million ha worldwide, representing 7% of the world forests [3] and are expected to dominate future industrial wood supply [4].
The main reasons for controlling vegetation in boreal forests are to maintain conifer composition and to increase conifer yield. After harvesting, pure jack pine (Pinus banksiana Lamb.) stands may be converted to mixed softwood or mixedwood conditions because of trembling aspen (Populus tremuloides Michx.) ingress [5]. Previous studies from across North America have shown that conifer release post-harvest leads to increases in conifer volumes from 49 to 5478% compared to untreated controls [6], while in Ontario, Canada, conifer release altered individual conifer growth between −49 and 556% compared with trees in untreated controls [7]. According to the Montréal Process sustainable forest management also includes maintaining biodiversity [8]. Trade-offs between biodiversity and productivity have been reported to occur in forest ecosystems. On average, a 10% loss in biodiversity leads to a 3% loss in productivity, although the most productive stands are typically monocultures [9]. Since vegetation management treatments are applied to enhance productivity through removal of competing, or undesirable, species, it seems reasonable to assume biodiversity will be reduced even as productivity increases.
The application of vegetation management treatments such as herbicides or brushsaws involved making trade-offs between criterion, such as stand-scale biodiversity and productivity [6,10,11]. While it is often expected that herbicides reduce biodiversity, e.g., [12], empirical studies have shown that herbicide applications do not typically result in species loss, but rather in shifts in abundances and evenness [13,14,15]. Tenth year post-treatment results have generally supported the hypothesis that forest vegetation release treatments do not produce tree monocultures [15]. However, under some conditions, stand-scale tree monocultures could result from multiple herbicide applications. Longer-term studies are needed to determine if these patterns persist in older plantations.
By properly managing stands, the potential exists to produce healthy trees with better stem form and wood quality to feed current and future industrial needs. Forests generally grow based on resources available to the individual trees. When resources are limited, tree growth slows, and different areas of growth are prioritized than when resources are plentiful. This simple idea also determines, in many cases, the wood properties we see in trees. For example, boreal forests with over 3000 established stems ha−1 at the renewal stage typically produces good wood quality and stem form [16]. In contrast, a forest with less than 3000 established stems ha−1 at the renewal stage produces less desirable wood properties and stem form, with this effect increasing as tree density decreases [17,18]. Stand density strongly influences wood properties such as wood density, number of growth rings per cm, timing of juvenile to mature wood transition [19], stem taper, stem branchiness, total height, diameter at breast height (dbh; 1.3 m), and wood mechanical properties [20]. Typically, as stand density decreases (i.e., below 3000 stems ha−1), stem and wood properties worsen [16,21].
Best forest management practices must address competing demands, such as successful growth of crop trees following aerial herbicide release and increasing public concern about herbicide use in forestry. Thus, the objective of this study was to compare the effects of three release treatments (aerial herbicide spray, motor-manual brushsaw, complete competition removal via herbicides applied by backpack spray) and a control treatment (no vegetation control) on three individual criteria: [1] plant diversity, [2] stand-scale yield, and [3] individual crop tree responses at two jack pine sites. We hypothesized that long-term effects among the three release treatments on individual criteria would be similar, and that all release treatments would have a significantly positive effect on the three criteria compared with the control. We also hypothesized that temporary vegetation control would not lead to overstory monocultures at our study sites. We predicted that, in more intensive treatments (i.e., complete removal), individual crop tree survival and stem quality would increase, but plant diversity and branching would decrease. Furthermore, more intensive treatments would yield more conifers and less hardwood per ha than less intensive treatments. Here, we present data collected from two sites, established in 1992 and 1993, in northern Ontario, Canada, under the Vegetation Management Alternatives Program (VMAP; [22]). This research continues to address the lack of longer-term data needed by forest managers to implement more integrated vegetation management programs in boreal and sub-boreal forests.

2. Materials and Methods

2.1. Experimental Design

A field study was established in jack pine plantations in 1993 and 1994 on two sites to examine the effect of three release treatments and a control on individual crop tree growth, stand-scale growth and yield and plant diversity (Table 1). The first study site (Bending Lake) is located 53 km north of Atikokan, Ontario, Canada (48°57′ N, 92°02′ W) near the southwest shore of Clearwater West Lake. The second study site (E.B. Eddy) is about 120 km northwest of Sudbury, Ontario, in Olinyk (46°47′ N, 82°08′ W; 46°47′ N, 82°10′ W) and Oshell (46°46′ N, 82°04′ W) townships. The study was established as a randomized complete block design with four treatments (including a control) on each of seven blocks (Table S1). Four blocks were at the Bending Lake site and three blocks were at the E.B. Eddy site. The four treatments, in increasing order of intensity, were:
  • No competition control (C), in which, aside from site preparation, post-logging vegetation was left undisturbed.
  • Motor-manual brushsaw cutting (BS), with non-crop vegetation cut at ground level with a motor-manual brushsaw in June 1993 (Bending Lake) or with non-crop vegetation cut at 25 cm above groundline with a Husqvarna 165 clearing saw in October 1994 (E.B. Eddy).
  • Operational single aerial herbicide application (AS), with 4 L ha−1 of Vision® (1.5 hg acid equivalent [a.e.] ha−1) applied in a total spray volume of 35 L ha−1 using a G47T helicopter in late August 1992 (Bending Lake), or with 4 L ha−1 of Vision® (356 g L−1 a.e. glyphosate as a water soluble liquid) applied in a total spray volume of 34 L ha−1 using a Bell 206 helicopter in late August 1993 (E.B. Eddy).
  • Complete competition removal (complete removal; CR), with annual backpack spray applications of a 2% solution (400 mL glyphosate in 20 L water) of Vision® applied in August 1993, June 1995, and June 1996 (Bending Lake) or September 1993 and August 1994 (E.B. Eddy). The amount of Vision® sprayed via backpack application varied among plots depending on how much vegetation survived previous applications; however, comparing the efficacy of different Vision® rates was beyond the scope of this study. Each complete removal treatment plot was randomly located in a single aerial spray plot.
Complete competition removal in not typical of operational practices and is included as an extreme for the purpose of this study. Each treatment plot was a minimum of two ha at Bending Lake (for details, see [23]). At E.B. Eddy, each treatment plot was 50 × 100 m, except for herbicide application (100 × 200 m) and complete removal (40 × 40 m) plots, and berms covered about 5 to 10% of each plot area. In each replicate, treatments were randomly allocated to plots.

2.2. Vegetation Assessment

Before treatments were applied, 20 representative jack pine seedlings were systematically located within each assessment area. Each crop tree was identified with an aluminum tag attached to a numbered assessment pin at its base. At Bending Lake, selection was based on a 10 × 10 m grid of four trees by five trees. At E.B. Eddy, selection was based on an approximated 4 × 4 m grid of four trees by five trees with some trees closer or further apart [23]. Crop trees were recorded as dead, missing, healthy, or unhealthy. For each crop tree, other than the missing ones, total height, height to live crown, diameter at breast height (dbh), crown width, diameter of the five largest branches 50 cm above and below dbh, branching, and stem quality were assessed. After 25 growing seasons, biodiversity plots (10 × 10 m) were established off crop tree pins in a NW, NE, SW, or SE direction. Three of the 20 tree pins were chosen to provide maximum distance among plots. Vegetation was separated into seven layers (dominant tree cover, subdominant tree cover, 2–10 m trees and shrubs, 0.5–2 m trees and shrubs, <0.5 m trees and shrubs, herbs/grasses/sedges/ferns, and bryophytes/lichens) [13,24]. Biodiversity was identified to the species level (or genus if species unknown) and percent cover was estimated. Three replicate growth and yield plots (11.28 m radius) were established off crop tree pins in each treatment plot, except for complete removal plots at E.B. Eddy (which only had two replicate plots due to limited space) to measure forest trees with a dbh ≥ 2.5 cm. To prevent overlap, a minimum 22.56 m radius was required between growth and yield plot centers. Growth and yield surveys were conducted in accordance with the Ontario Growth and Yield Standards as outlined in the Permanent Growth Plot and Permanent Sample Plot Reference Manual [25]. In total, data for 27 growth, yield, and biodiversity variables were collected at the end of the 25th growing season following initial treatment application: 18–22 July 2018 at Bending Lake site and 16–20 July 2019 at E.B. Eddy site. Further explanation of all variables assessed are provided in Table S2.

2.3. Statistical Analyses

To assess differences in species diversity following release treatments, compositional indices (species richness (S), abundance (abun), Shannon entropy (H’), Simpson’s dominance (D) and Heip evenness (EHeip) were calculated using the methods described in [28]. Before calculating the indices, the seven canopy layers assessed in the field were condensed to four: L1 (dominant and subdominant tree cover), L2 (<0.5–10 m shrub layer), L3 (herbs, grasses, sedges and ferns), and L4 (bryophytes and lichens).
Diameter at breast height was measured for all live and dead trees, but height was recorded for only three to five trees per growth and yield plot at both sites. However, only live trees were used to calculate mean total height, quadratic mean diameter (QMD), dbh, height to live crown, branch diameter, rings per cm, inside and outside bark volume, basal area and density (number of trees), and merchantable volume (10 cm top diameter) in each plot. Merchantable volume of individual trees were calculated using Honer’s conversion factors for jack pine [29]. Since height information is needed to estimate stem volume, heights of all live trees were estimated using a height-diameter model developed by [30]. Similarly, heights of hardwood trees were calculated using the models by [31]. Inside bark volumes of conifer trees were calculated using recently published volume equations [32,33] and inside bark volumes of hardwood trees were estimated using Honer’s volume equations [34]. Total inside bark volumes, basal area, and number of trees from each plot were converted to per hectare values for both conifer and hardwood trees.
Statistical analyses were performed with R statistical software (Version 4.0.3, Vienna, Austria) [35] and the lme4 function in the package nmle. Three-way analysis of variance (ANOVA), implemented through generalized linear mixed effects models (GLMM), was used to determine the effects of three release treatments and a control on (1) plant diversity characteristics, (2) stand-scale growth and yield characteristics, and (3) individual crop tree growth characteristics. For each analysis, treatment and site were fixed, and block was the random, independent variable. Residual analysis was performed to test for the assumptions of the three-way ANOVA. Normality was assessed using histograms and the Shapiro-Wilk’s normality test and homogeneity of variances was assessed by Bartlett’s Test. To meet assumptions of normality, the variables rings per cm, conifer density, and hardwood density were log transformed, and the variables gross total volume (GTV), basal area, merchantable volume, and individual jack pine volume were transformed using the cubed root function. When assumptions of normality were not met, we used nonparametric tests when appropriate. Poisson-distributed GLMMs were used to test for the effect of ‘treatment + site + block’ on species richness. Akaike Information Criterion corrected (AICc) score was used to select the most parsimonious model (using the R MuMIn package). p-values for normally distributed GLMMs were calculated using a type “II’’ mixed-model ANOVA, and p-values for Poisson-distributed GLMMs were calculated using the R emmeans package. Statistical significance was accepted at Bonferroni adjusted alpha levels of p < 0.0125 and p < 0.0167 for simple two-way interactions at Bending Lake and E.B. Eddy, respectively. Main effects analyses (pairwise comparisons) were determined using Tukey’s honest significant difference (HSD) test with 95% confidence. Ordinal dependent variables (survival, stem quality, branching) were assessed using ANOVA with planned orthogonal contrasts [36]. Proportion of surviving stems was estimated to calculate individual jack pine survival. Ordinal dependent variables were transformed using the Anscombe transformation to meet assumptions of normality. The orthogonal contrasts compared were within treatments (CR vs. AS, CR vs. BS, AS vs. BS) and between control and treatments (C vs. (CR + AS + BS)/3).
The linear model for the ANOVA was:
Yijk = µ + Bi + S + T + åijk,
where:
Yijk is the calculated response from ith block, jth site, and kth treatment,
µ is the overall mean,
Bi is the random effect of the ith block,
S is the fixed effect of the jth site,
T is the fixed effect of the kth treatment, and
åijk is the pooled interaction effect of the ith block, the jth site, and the kth treatment with error term to test the treatment effect.
Statistical significance of differences among treatments was tested by pooling the interaction term with the error terms in the model. Characteristics for each individual criterion (plant diversity, non-crop tree growth and yield, and individual crop tree growth) were compiled as a percentage deviation from the baseline (ideal/maximum site conditions) (Table S2) and portrayed graphically in the form of radar diagrams. Radar diagrams are recommended as one way to handle multivariate data [37], and in this study the diagrams clearly illustrate differences among the effects of release treatments at the two sites.

3. Results

3.1. Plant Diversity

Plant diversity characteristics are displayed in Figure 1a. The effect of site was significant for softwood richness, tree abundance, herb abundance, and understory evenness (Table S3). Jack pine and trembling aspen as well as five other softwood tree species and four other hardwood tree species (>10 m) were present at both sites: balsam fir (Abies balsamea [L.] Miller.), white spruce (Picea glauca [Moench] Voss), black spruce (Picea mariana [Miller] B.S.P.), red pine (Pinus resinosa Sol. ex Aiton), white pine (Pinus strobus L.), red maple (Acer rubrum L.), white birch (Betula papyrifera Marshall), largetooth aspen (Populus grandidentata Michx.), and pin cherry (Prunus pensylvanica L. f.) (Table S4). However, softwood richness was significantly higher at E.B. Eddy than at Bending Lake in all treatment and control plots, except complete removal. For example, average softwood richness by treatment ranged from 1.8 to 2.6 at Bending Lake compared with 3.0 to 4.2 at E.B. Eddy (Table 2). Hardwood richness was generally higher in all treatment plots than control plots at both sites, and tree abundance (% cover) was significantly higher in treatment plots at E.B. Eddy than in those at Bending Lake. Specifically, complete removal, aerial spray, and motor-manual brushsaw plots at E.B. Eddy had 11.1, 17.4, and 20% higher tree abundance, respectively, than their counterpart plots at Bending Lake.
Understory species richness was not significantly affected by vegetation management treatments and only one exotic species was observed (orange hawkweed, Pilosella aurantiaca (L.) F.W. Schultz & Sch. Bip.) in the herbicide treated plots at E.B. Eddy. However, no exotic species were observed at Bending Lake, and percent cover of the exotic species E.B. Eddy was low (0% in control and motor-manual brushsaw plots, 0–0.1% in aerial spray plots and 0–1% in complete removal plots). Aerial spray and complete removal plots at Bending Lake averaged 14 and 15 herbaceous species per plot, respectively, which was slightly higher than the 12 species found in both motor-manual brushsaw and control plots. Understory diversity and evenness were higher in treated than control plots at E.B. Eddy; however, differences were not significant. At Bending Lake, control plots had the highest understory evenness.

3.2. Stand-Scale Yield

Stand-scale yield characteristics are displayed in Figure 1b. Treatment effect was significant for all stand-scale characteristics (Table S5). Gross total stem density was highest in control plots at both sites, ranging from 3500 to 4000 stems ha−1 (Table 3). Gross totals were higher at E.B. Eddy in all treatments except motor-manual brushsaw, which averaged 2910 stems ha−1 at Bending Lake and 2331 stems ha−1 at E.B. Eddy. All treatment plots had significantly higher conifer density, basal area, and GTVs than control plots at both sites. Motor-manual brushsaw plots had the highest number of conifers ha−1 at both sites, with 2498 and 1842 trees ha−1 at Bending Lake and E.B. Eddy, respectively. Between sites, Bending Lake plots (except complete removal) had more conifers ha−1 and higher conifer basal area ha−1 than those at E.B. Eddy, and were significantly higher for motor-manual brushsaw and aerial spray plots. At Bending Lake, motor-manual brushsaw and aerial spray plots had 656 and 500 more conifers ha−1, respectively, compared with control plots that had 275. Opposing trends were evident for hardwood stand characteristics. Control plots had significantly higher hardwood densities, basal area, and inside bark volumes than treatment plots at both sites, and E.B. Eddy had the highest values for all three hardwood stand characteristics. For example, aerial spray, motor-manual brushsaw, control, and complete removal plots had 4.7×, 2.2×, 1.3×, and 1.2× more hardwoods ha−1 at E.B. Eddy than at Bending Lake.

3.3. Individual Jack Pine Characteristics

Individual jack pine characteristics are displayed in Figure 1c. Controlling for site and block, treatment effect was significant for survival and branching, but not stem quality (Table 4). Effect of site was significant for survival, stem quality, and branching. Planned orthogonal contrasts revealed that treatment plots were associated with significantly higher survival rates compared with control plots. For example, jack pine survival rates in control plots were 47.5% and 36.7% at Bending Lake and E.B. Eddy, respectively, but ranged from 58.8–72.5% and 78.3–84.7% in treatment plots. However, control plots had significantly less branching compared with treatment plots, and higher treatment intensities (i.e., aerial spray, complete removal) were associated with the highest branching. Pairwise comparisons between sites revealed that trees in motor-manual brushsaw and complete removal plots at E.B. Eddy had significantly higher survival rates (26% and 16% higher, respectively) than those at Bending Lake (Table S6). Furthermore, stems were rated lower quality in all plots at Bending Lake by a difference ranging from 0.6 to 1.0 units (on a scale of 1 to 4) and there was significantly more branching in control plots at Bending Lake compared with E.B. Eddy.
Site and treatment significantly affected total height, rings per cm, and merchantable volume, and the interaction effects between site and treatment were significant (Table S7). Treatment also significantly affected quadratic mean diameter (QMD), height to live crown, and branch diameter, but the site x treatment interaction effects were not significant. Aerial spray plots had the tallest trees at both sites (12.7 m at Bending Lake and 12.9 m at E.B. Eddy) and all treatment plots had significantly taller trees (p < 0.05) than control (untreated) plots (Table 5). Between sites, trees at E.B Eddy were significantly taller in control (10.8 m) (p. adj < 0.001) and motor-manual brushsaw (12.4 m) (p. adj < 0.001) plots than those in the Bending Lake control (9.5 m) and motor-manual brushsaw (11.4 m) plots. Rings per cm were significantly higher for trees in control (1.7–2.2) than in treatment (1.3–2.0) plots at both sites, and trees at E.B. Eddy averaged 0.3–0.6 more rings per cm than those at Bending Lake. Merchantable volume was significantly higher in treatment than control plots, ranging from 0.04–0.07 m3 tree−1 in treatment plots compared with 0.03 m3 tree−1 in control plots at both sites. Quadratic mean diameter was significantly larger in treatment than control plots at Bending Lake; however, at E.B. Eddy it was slightly higher in control and motor-manual brushsaw plots than in aerial spray and complete removal plots. Height to live crown was largest in aerial spray plots at both sites and was significantly larger (p < 0.05) in all treatments than in control plots at Bending Lake, and in aerial spray and complete removal plots at E.B. Eddy. Branch diameter was significantly larger (p < 0.05) in all treatments compared with control plots at both sites, and aerial spray plots at Bending Lake had significantly larger branch diameters (p. adj < 0.001) than those in aerial spray plots at E.B. Eddy.

4. Discussion

We analyzed long-term (25 year) individual crop tree growth and stand volume trade-offs between conifers and hardwoods in boreal and sub-boreal Canadian forests following three release treatments applied at two sites. Our results suggest that release treatments did not reduce plant diversity but benefitted stand yield and individual crop tree growth, with maximum benefits occurring in more intensive treatments.

4.1. Effect of Site and Species Ecology

Overall, site significantly influenced many of the diversity and productivity indicators. For example, tree species diversity and understory shrub and herb abundance and stand volumes were higher at E.B. Eddy than at Bending Lake; however, site-related differences must be interpreted cautiously. Although sampling intensity was identical, factors known to influence diversity and productivity differed between sites, for example, experimental plot size, site productivity and heterogeneity, and site preparation and planting timing. Relationships among these factors and species diversity are well known. Larger plots typically show higher species richness than smaller ones [38,39] due to the potential for higher heterogeneity in larger plots. Treatment plots at Bending Lake were minimally two times the area of those at E.B. Eddy; however, the sample area for each treatment was consistent at both sites. The higher diversity indices observed in some plots at E.B Eddy compared with Bending Lake are likely because the former is relatively more fertile based on its geological deposit and soil particle distribution [40]. While high levels of soil N usually decrease species richness (e.g., [41]), plant diversity generally increases from low fertility sites (here represented by the coarse loamy to fine sands at Bending Lake site) to intermediate fertility sites (here represented by very fine sands to silty loams at the E.B. Eddy site) based on a humped-back shaped relationship (e.g., [42]).
Silvicultural treatment differences, such as site preparation method, stand establishment density, and delays in temporary vegetation control, are confounded with site differences. Bending Lake and E.B. Eddy sites were trenched and bladed, respectively. As aspen roots are mainly contained within the top 12 cm of mineral soil [43], site preparation via trenching, which disturbs the mineral soil (not just exposes it, i.e., blading), may stimulate aspen regeneration through reducing the flow of sucker-inhibiting hormones within the roots, encouraging sucker formation [44]. Shade-intolerant hardwoods such as trembling aspen regenerate best in full light conditions [45,46] and moderately disturbed, fertile soils [47]. Both trenching and blading create a berm microtopography, but on our study site blading created 3 to 4 m wide berms that covered 5 to 10% of sample plot area. These wide berms were not planted with jack pine and were not treated during the brushsaw treatment; thus, they created nutrient rich refugia for hardwood and understory species in non-herbicide treated plots. Ref. [24,48] noted that untreated strips act as refugia. Jack pine was planted at 1.8 × 1.8 m at both sites, but due to ingress of naturals establishment densities were higher at Bending Lake. Stand density is known to influence branch habit [49], wood quality, and stand volumes [50,51]. Lastly, temporary vegetation control was delayed by two to four years at E.B. Eddy and by four to five years at Bending Lake (Table 1 and Table S1). Applying release treatments this long after harvesting allows establishment of understory vegetation such as grasses, forbs, and shrubs, which compete with shade-intolerant conifers such as jack pine, subjecting the latter to increased growth loss and mortality [23,52]. Thus, lower crop tree survival in all treatment plots at Bending Lake compared with E.B. Eddy likely resulted from the combination of low fertility and delayed temporary vegetation control at the site.

4.2. Treatment Efficacy

To fully understand the effects of vegetation management treatments on diversity and productivity, knowledge of treatment efficacy (i.e., cost and outcome) is needed. Based on a study of the cost effectiveness of controlling competitive vegetation such as aspen, [15] reported that herbicide treatments applied to white spruce plantations are three times more cost effective than cutting treatments. Ref. [53] included 10th year Bending Lake and E.B. Eddy data in internal rates of return analyses and observed a 4.32% and 2.90% return for aerial herbicide and brushsaw treatments, respectively. Their findings also indicated that both brushsaw and herbicide treatments can be economically viable.
Brushsaw treatments are typically considered non-economical because at least two treatments [54,55] are often needed for success. Our results show that a single treatment application can be highly effective at reducing aspen densities; however, the cutting treatment at Bending Lake was more effective for controlling aspen than that at E.B. Eddy. We presume the difference in effectiveness can be attributed to differences in season, height of cut and site related differences. Research by [55] on aspen suggested that cutting in mid-summer at 50 cm, as applied at Bending Lake, is more effective for controlling aspen than cutting in early fall at 25 cm, as applied at E.B. Eddy.
Efficacy of glyphosate applications is known to be strongly influenced by seasonal timing and rate of application [56] and post-plant timing [52,57]. Although herbicide treatments were applied within the optimal seasonal time period and rates of application, delays in control of competitive vegetation may have significantly reduced the survival and growth of individual jack pine and other conifers. The critical period for jack pine, defined as the time during crop development when interspecific competition from weeds must be controlled to prevent significant yield losses, is two to three years after planting [6]. In our study, the period between planting and herbicide release treatments ranged from two to five years. Delaying temporary vegetation control by as much as two years after the critical period for jack pine may explain the low crop tree survival across all treatment plots. Applying herbicides within the critical period may increase crop tree survival and growth; however, early release may decrease the amount of understory plants such as blueberries [14].

4.3. Treatment Effects on Diversity

Our findings generally support our hypothesis that vegetation management, within the large intensity gradient we sampled, did not lead to overstory monocultures in boreal and sub-boreal forests. In fact, fifth year tree diversity data at E.B. Eddy predicts jack dominated stands would be established following all treatments [23]. However, we found increasing intensity of vegetation management appeared to increase species richness, particularly for hardwood species. This result may be due to the overall reduction in hardwood density in the treatment plots, which might have disrupted the typical dense vegetative reproduction of poplar and allowed other shade intolerant hardwoods time to establish before access to light became limiting. Despite efforts at complete removal in the most intensively treated plots, hardwoods persisted, which supports the hypothesis that tree monocultures do not develop after release treatments [15]. While multiple herbicide applications can reduce hardwood presence post-harvest, as applied in this study they did not remove them completely, probably because of their overall persistence in seedbanks and through sprouting as suggested by the direct regeneration hypothesis [58,59]. The fact that other vegetation components remained in these areas after 25 years suggests that aerial spray and complete removal treatments may offer options for alternative management outcomes. However, verification that these species persist through continued stand development requires further investigation.
Intensification of vegetation management treatments also had positive, or non-significant, effects on shrub and herb layer richness, which was unexpected. This finding may be the result of increases in light availability for the understory through decreases in hardwood proportion in treatment plots [14,60], an increase in substrate diversity from repeated minor disturbance and removal of smaller trees, or the increased overstory diversity, all of which increase vascular plant species richness and shrub species richness in boreal forests [61,62]. Indeed, high stand density, and resulting reduced light, can override the generally positive effects of mixed forests on understory diversity [63]. Non-significant effects of intensification on evenness suggest that while intensification improved species richness of the understory, it did not cause only few species to dominate. Understory vegetation cover in aerial spray and complete removal plots included understory plants such as wild blueberries (Vaccinium angustifolium Aiton) and (Vaccinium myrtilloides Michx.), sweet fern (Comptonia peregrina [L.] J.M. Coult.), beaked hazel (Corylus cornuta Marshall), wild red raspberry (Rubus idaeus L.), bush honeysuckle (Diervilla lonicera Miller), Alnus spp., and Salix spp. These plant species are vital wildlife food sources and support previous research that suggested herbicide treatments maintain wildlife food source availability [23]. Although these results suggest that most understory plant species in mid-successional boreal and sub-boreal forests can withstand intensive silvicultural management regimes, public acceptance of herbicides is low [64]. If herbicide use is not permitted in the management area, or if a chemical-free alternative is preferred, motor-manual brushsaw may offer options for managing mixedwood stands [23]. In either case, vegetation control to maximize planted crop tree growth is not incompatible with maintaining diverse understory communities [65]. However, future studies should consider if or how the understory community composition is altered by release treatments.
The intermediate disturbance hypothesis suggests that diversity should be highest when there are moderate levels of disturbance [66]. However, 25 years after harvesting at our study sites, understory species diversity mostly increased with increasing management intensity. In the first three years where herbicides were ground applied annually in complete removal plots, deciduous tree and tall shrub species abundance was lowest in complete removal plots [23], but quickly rebounded, as glyphosate has no effect on plants germinating post-treatment. If glyphosate application in complete removal plots had continued, these plots would likely have lower species diversity, as only highly tolerant species or good colonizers would persist [67]. Furthermore, multiple studies provide evidence that anthropogenic disturbances such as harvesting and silviculture can facilitate invasions by exotic plant species in boreal forests [13,14,68]. Shade-intolerant exotic species, such as common dandelion (Taraxacum ceratophorum [Ledeb.] DC.), colonize the open spaces created immediately after disturbance [14]. Such species, however, were generally not recorded in the closed canopy conditions observed 25 years post-harvest.

4.4. Treatment Effects on Stand-Scale Yield

Our findings generally support our hypothesis that more intensive treatments would yield higher conifer volume and less hardwood volume than less intensive treatments. Jack pine density, basal area, and stand volume were all significantly higher in treatment plots compared with the control; however, volumes were higher for single applications of glyphosate than the more intensive multiple applications. Why the complete removal treatment did not produce maximum GTV as predicted is unknown; however, the most probable explanation is that jack pine simply allocated excessive amounts of energy to branch growth. Removal of competition at ground level, combined with low stand densities, likely allowed trees in complete removal plots to allocate more wood to branches to support more foliage. This explanation is supported when comparing crown width values among treatment plots. Complete removal plots had the largest crown width (4.6 m) relative to other treatment plots (3.1 m) and control plots (2.6 m) (unpublished data).
Our findings suggest that jack pine volumes could be increased through vegetation management, but there may considerable opportunities to further improve the yield of these jack pine plantations, possibly through a combination of vegetation management and tree improvement. The mean annual volume increment (MAI) for treated jack pine plantations in this study was, on average, 6.3 m3 ha−1 in aerial spray plots and 2.6 m3 ha−1 in control plots. While these numbers clearly show a significant beneficial role of vegetation management, our highest values are still far lower than the MAI of 10.0 m3 ha−1 reported for the best jack pine provenances growing in Ontario [69]. Our findings also suggest that current forest management models, such as Forest Vegetation Simulator (FVS) [70], may underestimate volume yields. Using 10th-year post-treatment data from the Bending Lake and E.B. Eddy sites, [53] applied FVS to predict GTV to an arbitrary age of 70 years. For the period of this study (i.e., 25 years post-treatment), they predicted that GTV for the untreated control and complete removal treatments would yield the lowest and highest volumes, respectively, and that the GTV of these stands would be between 50 and 100 m3 ha−1. These predictions were well under actual values. For example, in the Bending Lake study, the GTVs ranged from 75 m3 ha−1 for the control to 161 m3 ha−1 for the aerial spray. The volume gains associated with the complete removal treatment at year ten post-treatment did not persist as predicted. We presume that underestimation of GTV is attributed to model functions within FVS. This is important, because managed jack pine stands contribute about 600 million m3 GTV [71] to the 720 million m3 GTV of jack pine growing in Ontario [72]. Despite the contribution of jack pine plantations to Ontario’s wood supply, very few studies have reported on the long-term (i.e., greater than 20 year) growth of jack pine plantations [73]; thus, this study helps to fill in an information gap.
In this study, we calculated volumes (GTVs) using real data (i.e., height, diameter, etc., measured in the plots) by applying recently published volume equations developed for plantations. Therefore, GTVs used this analysis are more accurate than those projected using FVS. The volume models used in the FVS calibrated for Ontario incorporated models developed for natural stands and volumes of trees grown in plantations differ from those in natural stands [32,33]. Ref. [71] suggested that the use of tools to estimate the production potential of managed stands may significantly underestimate their productivity. Ref. [73] confirmed this bias for a study of jack pine plantations near Nipigon, Ontario. They reported a mean GTV of 307 m3 ha−1 (range of 2–368 m3 ha−1).
In addition, volumes may be enhanced through complementarity. Recent research on plant species complementarity suggests that productivity of conifer stands may be improved through mixed species plantations [74]. The possibility mixed species forest stand productivity might differ from that of pure stands has long intrigued growth and yield researchers [75]. Ecological theory suggests that when two or more species use resources differently, the resources may be more completely exploited by plants growing in mixture than by those growing as single species, leading to higher productivity in mixed species stands [74,76]. However, recent research has highlighted that complementarity can only occur when species do not overlap in resource acquisition. For example, contrasting shade tolerance levels among constituent species are needed to observe higher productivity with increasing diversity globally [77], and in the boreal forest specifically [78]. This discrepancy may be why stand-scale volumes were higher in near pure stands of jack pine than in intimate mixtures of aspen and jack pine, given that both are intolerant pioneers. Though these stands are relatively young and diversity effects have been frequently observed to increase through time [59,79,80,81] it is unlikely that trees in the control plots (i.e., plots with higher proportions of mixture) would be able to outgrow the near pure jack pine stands in an operational timeframe. Additionally, these stands have reached canopy closure, at which point positive diversity effects become most evident [81]. Future studies should consider mixtures of species that complement jack pine in shade tolerance [78], rooting structure [82], or nitrogen fixation [83]. We found no evidence that mixtures of aspen and pine increase tree growth. Increased productivity in mixtures has not been consistently demonstrated by analyses of the effect of composition on stand growth [75,84].

4.5. Treatment Effects on Individual Jack Pine Characteristics

4.5.1. Survival, Height, and Diameter

Jack pine survival, height, and quadratic mean diameter were all significantly greater in treatment compared to control plots. These findings are generally in line with our hypotheses as well as previous studies that have shown that jack pine survival, height growth, and diameter can be strongly influenced by competition with trembling aspen and other woody species. Typically, jack pine height growth and diameter are more sensitive to competition than survival [7,52,55]. This hierarchical pattern of competition affecting conifers has been well documented in the literature [85]. However, in our study we found a substantial positive effect of treatment on jack pine survival, particularly at E.B. Eddy.
We attribute the relatively low survival rate of jack pine in the controls at the E.B. Eddy site to severe competition from overtopping aspen. While water and nutrients can be limiting for natural jack pine stands [86], after disturbance, light is a limiting resource for individual jack pine [87]. Furthermore, competition for light is asymmetrical, meaning that larger trees obtain disproportionate access to light for their size [88]. After clearcuts without vegetation management, intolerant hardwoods such as trembling aspen can easily recruit at high densities and can outgrow planted jack pine [5]. Once jack pine seedlings are overtopped, it becomes increasingly difficult for them to access light, reducing survival and height and diameter growth. Our results reinforce the assertion that if forest managers aim to regenerate a high proportion of jack pine, vegetation management is essential.
Our findings are comparable to those of [6], who reported gains in wood volume of 49–5478% in northern forests. The results of the single unpublished jack pine study reported in that review showed wood yield gains at 10 years post treatment of 116%, compared to those found in our study, which were wood gains at 25 years post treatment of 288–390%. Ref. [57] have since confirmed that survival, height growth, and diameter of jack pine are increased by herbaceous vegetation control. Continued monitoring of these sites through rotation is essential to determine whether increases in merchantable volume are maintained, or amplified, as stand development occurs.
Lack of significant difference in conifer response to herbicide and brushsaw treatments has been previously reported [7,89]. We attributed the lack of differences in this study to the fact that both the herbicide and mechanical cutting treatments were applied in the proper season and post-treatment jack pine densities were relatively high. Stand density is known to affect individual jack pine survival and growth. When natural pure (or nearly pure) jack pine stands are too dense, survival and height may be reduced due to high intraspecific competition for light and nutrients [90]. Low stand densities can also decrease jack pine height [91], because individuals will allocate more mass to branches when they are not competing as intensely for light [92]. This increased allocation to branches leads to a deterioration in stem form [69]. At optimal stand densities, lower branches on the stem decrease in growth and/or die from shading and competition [21]. Consequently, more mass may be allocated to the stem, increasing diameter. Stem diameter is a major wood quality attribute for the sawmilling industry. With increasing stem diameter, logging and manufacturing costs decrease, whereas lumber volume recovery and grade yield increase significantly.
Despite observed increases in height and quadratic mean diameter in our treatment plots, stand densities were below those recommended (3000 stems ha−1) to maintain jack pine quality wood [17,21], suggesting further room for improvement. Planting density (1.8 × 1.8 m) at our sites was within spacing recommended to maximize jack pine growth and productivity. However, this initial spacing provides the desired stem density of 3000 stems ha−1 when most seedlings survive. Clearly, control plots did not meet this density since jack pine mortality rates were high, possibly due to ingress of aspen and subsequent interspecific competition for resources, reaffirming that site preparation and planting alone cannot regenerate near-pure jack pine stands. Treated plots also did not achieve optimal density, indicating that higher planting densities, tree improvement focused on increasing seedling survival, pre-harvest treatment, or further silvicultural interventions (such as thinning of undesirable species) need to be investigated.

4.5.2. Stem Form, Knot Size, and Branch Diameter

Silviculture activities in jack pine forests have profound effects on how the trees grow and the resulting wood properties. Jack pine is a species that typically has many thousands to tens of thousands of trees on a site after a disturbance, with fire being the main natural disturbance affecting this species. In a study by [93], increasing tree spacing increased jack pine crown length, live crown ratio, and crown width. Additionally, branches were larger in the widest spacing as were stem taper, sweep, and wobbling. This finding provides evidence that jack pine is not a species that benefits from excess space. All these factors lead to less desirable wood for lumber production. In a study by [20] increased black spruce spacing resulted in similar trends for stem features as those found for jack pine, i.e., decreased lumber grades and reduced yield of high-quality lumber. The main downgrades were due to knots and wane, which result from larger knots and more taper in widely spaced trees. They also found that in the lowest stand density, the highest grades of lumber were not found in any of the trees, meaning the value of each tree was significantly lower. In both the [93] and [20] studies, the lowest stand density produced the lowest quality lumber and wood quality was not comparable to that of trees from natural stands.
In many Canadian provinces, jack pine is an important commercial species and represents a substantial component of the total softwood harvest. The wood is used in the lumber and pulp and paper industries, with lower logs going to sawmills and upper logs to pulp mills. The lumber is strong and rated for structural building, while the fiber has enough length for quality paper. Much of the strength properties of jack pine are based on its wood properties and, as discussed above, stand density affects both wood mechanical properties and density [20]. As stand density decreases more wood is produced in each tree, but the mechanical strength of the wood decreases. This finding is significant and requires a more detailed study of wood properties at different stand densities to grasp the full effect of stand density management. Ref. [53] also found that stand density management increased jack pine attributes such as height and diameter, GTV, and projected value of fiber relative to controls.
Although stem quality was unaffected by vegetation management, average stem form was worse at our sites than in similar studies reporting on extrinsic tree quality. For example, sixteen years after herbicide and cutting treatment on a white spruce plantation in northwestern Ontario, 80% of trees in all treatment plots were classified as having good stem form [36] compared with 57, 58, and 62% of living trees in complete removal, motor-manual brushsaw, and aerial spray plots, respectively, after twenty-five years. This result may be due to low stand densities as previously mentioned. Based on a jack pine density management diagram for boreal Ontario [18,94], forest managers should aim for stand densities of at least 3000 stems ha−1 to ensure early crown closure, which is required to minimize branchiness [90]. Therefore, increasing initial planting densities (above the standard 1.8 × 1.8 m spacing that yields 3084 stems ha−1) to account for tree mortality may allow forest managers to optimize stand merchantable volume. Another reason for reduced stem quality and increased branching in treatment plots in our study may be the genetics of our planted crop trees. Jack pine does not usually grow perfectly straight, resulting in irregular shapes, and is usually knotty. However, certain characteristics of jack pine (i.e., growth, branching habit) can be genetically selected for to get better quality wood and increases in overall wood volume [95].
Information about the effects of stand density on wood quality is important to move management towards maximizing growth while maintaining wood quality. Research is underway to look at genetic trials for growth combined with wood quality research to link the growth rates and tree properties to the internal wood properties as a means of selecting families that show increased growth rates while maintaining wood quality.

5. Conclusions and Recommendations

This study offers one of the longest-term evaluations of operational vegetation management options for releasing sub-boreal and boreal jack pine plantations (aerial herbicides, motor-manual cutting), and it includes extremes (untreated control, complete removal) as benchmarks. Moreover, it provides information about treatment effects on plant diversity, stand-scale productivity, and tree-scale growth and quality, in interaction with inherent site characteristics. Our results suggest that plant diversity and productivity in boreal jack pine forests are significantly influenced by vegetation management treatments; however, none of the treatments maximized all 27 indicators studied; therefore, forest managers selecting treatments are faced with trade-offs. With the increase in mass timber production, we must maintain forests that will meet the grades for existing and new products and provide a constant sustainably managed supply to mills. Future studies should investigate the longer term (ideally through at least one rotation) effects of vegetation management treatments on higher jack pine densities (whether through increased initial planting densities, more vigorous thinning regimes, or infill planting).
Research on the criteria and indicators for the conservation and management of boreal forests not discussed here is also needed to make informed policy decisions and progress toward sustainable forest management. These assessments should include forest ecosystem health, soil and water resource conservation, and the maintenance of socio-economic benefits [8]. Given our observed trade-offs among vegetation management treatments for plant diversity, stand-scale growth and yield, and crop tree growth and form, it is highly likely that these treatments have trade-off implications for other ecosystem services such as nutrient cycling and wildlife diversity [96,97,98]. As well, these treatments have socio-economic implications. Investigation of societal and economic values associated with vegetation management alternatives should focus on the maintenance of the “licence to operate” [99,100,101,102]. Coupled with this, recent research suggests that First Nations opposition to herbicide use is complex and multifaceted [103]. Further research is needed to understand these complexities, through collaborative work that meaningfully engages First Nations communities concerned about forest management and that weaves Indigenous ecological knowledge into projects informing vegetation management programs. In addition, research on intrinsic values such as wood density (ring density, earlywood ring density, and latewood ring density), percentage of latewood in the ring, and wood bending properties [104] may also be prudent to support both traditional (i.e., sawmilling, pulp and paper) and emerging (i.e., bioenergy, biochemical extractables) forest sectors.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/12/3/370/s1, Table S1: Summary of jack pine release treatment details applied on two sites in Ontario, Canada, Table S2: Data collected on crop trees, noncrop vegetation growth and yield, and biodiversity during the 25th growing season after aerial herbicide application, manual brushsaw cutting, and complete competition removal at two sites (Bending Lake and E.B. Eddy) in Ontario, Canada [105,106,107], Table S3: p values from three-way analysis of variances (ANOVAs) for biodiversity indicators (n = 85) after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Statistical significance was accepted at Bonferroni adjusted alpha levels of p < 0.0125 and p < 0.0167 for simple two-way interactions at Bending Lake and E.B. Eddy, respectively, Table S4: Species list for Bending Lake and E.B. Eddy sites 25 growing seasons after jack pine release treatments. Bolded names refer to species found only at the Bending Lake site and those with asterisk (*) were found only at the E.B. Eddy site, Table S5: p values from three-way analysis of variances (ANOVAs) for stand growth and yield indicators (n = 82) after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Statistical significance was accepted at Bonferroni adjusted alpha levels of p < 0.0125 and p < 0.0167 for simple two-way interactions at Bending Lake and E.B. Eddy, respectively, Table S6: Adjusted means for individual jack pine survival (n = 560), stem quality (n = 360), and branching (n = 360), grouped by site x treatment variables after tending treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Statistical significance at p < 0.05, calculated only for variables with a significant main effect of site. Stem and branching were rated on a scale of 1 to 4, with 1 having best form and least branching, Table S7: p values from three-way analysis of variances (ANOVAs) for individual jack pine growth indicators (n = 360) after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) located in Ontario, Canada. Statistical significance was accepted at Bonferroni adjusted alpha levels of p < 0.0125 and p < 0.0167 for simple two-way interactions at Bending Lake and E.B. Eddy, respectively.

Author Contributions

Conceptualization, H.D.D. and F.W.B.; methodology, H.D.D., F.W.B., E.B.S., M.L., M.S. and J.D.; software, H.D.D.; validation, N.T., E.B.S. and M.L.; formal analyses, H.D.D. and M.S.; investigation, F.W.B. and J.D.; resources, F.W.B. and J.D.; data curation, J.D.; writing—original draft preparation, H.D.D., F.W.B., N.T., E.B.S., M.L. and M.S.; writing—review and editing, J.D.; visualization, H.D.D. and F.W.B.; supervision, F.W.B.; project administration, F.W.B. and J.D.; funding acquisition, F.W.B. and J.D. All authors have read and agreed to the published version of the manuscript.

Funding

Data collections for this study were funded by the Ontario Ministry of Natural Resource and Forestry.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors acknowledge Bob Wagner for initiating and administering the Vegetation Management Alternatives Program (VMAP), without which the Bending Lake and E.B. study sites would not have been conceived. Trial establishments and treatment installations were administered by Chris Holstedt and Andrée Morneault at Bending Lake and E.B. Eddy, respectively. The authors also thank Justin Viljakainen, Cole Barban, Amy Bolduc, Madison McCaig, Elaine Mallory, Erika Mihell, Christine St. Jules, Ana Kostic, Kayla Hayden, Riley Stobie, and Sierra Rutland for their help in collecting the twenty-fifth-year measurements. Thanks also to Lisa Buse (OFRI) and Mike Brienesse (OMNRF) for their thoughtful reviews of the draft manuscript.

Conflicts of Interest

The following authors, H.D.D., N.T., M.L., M.S., J.D., and E.S. declare no conflict of interest. F.W.B., an employee of the funder, had a role in the design of the study; in the collection, analyses, and interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Gauthier, S.; Bernier, P.; Kuuluvainen, T.; Shvidenko, A.Z.; Schepaschenko, D.G. Boreal forest health and global change. Science 2015, 349, 819–822. [Google Scholar] [CrossRef]
  2. Natural Resources Canada. The State of Canada’s Forests: Annual Report 2019; Natural Resources Canada, Canadian Forest Service: Ottawa, ON, Canada, 2020; p. 84. Available online: https://cfs.nrcan.gc.ca/publications?id=40084 (accessed on 1 December 2020).
  3. FAO; UNEP. The State of the World’s Forests: Forests, Biodiversity and People; Food and Agriculture Organization of the United Nations: Rome, Italy, 2020; p. 214. [Google Scholar]
  4. Carle, J.; Holmgren, P. Wood from planted forest: A global outlook 2005–2030. For. Prod. J. 2008, 58, 6–18. [Google Scholar]
  5. Hearnden, K.W.; Millson, S.V.; Wilson, W.C. A Report on the Status of Forest Regeneration: Ontario Independent Forest Audit; Ontario Ministry of Natural Resources, Ontario Independent Forest Audit Committee: Sault Ste. Marie, ON, Canada, 1992; p. 117. [Google Scholar]
  6. Wagner, R.G.; Little, K.M.; Richardson, B.; McNabb, K. The role of vegetation management for enhancing productivity of the world’s forests. Forestry 2006, 79, 57–79. [Google Scholar] [CrossRef]
  7. Bell, F.W.; Dacosta, J.; Penner, M.; Morneault, A.; Stinson, A.; Towill, B.; Luckai, N.J.; Winters, J. Longer-term volume trade-offs in spruce and jack pine plantations following various conifer release treatments. For. Chron. 2011, 87, 235–250. [Google Scholar] [CrossRef] [Green Version]
  8. Canadian Forest Service. Criteria and Indicators for the Conservation and Sustainable Management of Temperate and Boreal Forests: The Montreal Process; Natural Resources Canada, Canadian Forest Service, Headquarters, Policy, Planning and International Affairs Branch: Ottawa, ON, Canada, 1995; p. 28. [Google Scholar]
  9. Liang, J.; Crowther, T.W.; Picard, N.; Wiser, S.; Zhou, M.; Alberti, G.; Schulze, E.-D.; McGuire, A.D.; Bozzato, F.; Pretzsch, H.; et al. Positive biodiversity-productivity relationship predominant in global forests. Science 2016, 6309, aaf8957. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Lieffers, V.J.; Armstrong, G.W.; Stadt, K.J.; Marenholtz, E. Forest regeneration standards: Are they limiting management options for Alberta’s boreal mixedwoods? For. Chron. 2008, 84, 78–82. [Google Scholar] [CrossRef] [Green Version]
  11. Willig, M.R. Biodiversity and productivity. Science 2011, 333, 1709–1710. [Google Scholar] [CrossRef] [PubMed]
  12. May, E.E. At the Cutting Edge—The Crisis in Canada’s Forests; Key Porter Books: Toronto, ON, Canada, 2005; p. 432. [Google Scholar]
  13. Bell, F.W.; Newmaster, S.G. The effects of silviculture disturbance on the diversity of seed-producing plants in the boreal mixedwood forests. Can. J. For. Res. 2002, 32, 1180–1191. [Google Scholar] [CrossRef]
  14. Mallik, A.U.; Bell, F.W.; Gong, Y. Effectiveness of delayed brush cutting and herbicide treatments for vegetation control in a seven-year-old jack pine plantation in north-western Ontario, Canada. Silva Fenn. 2002, 35, 505–519. [Google Scholar] [CrossRef] [Green Version]
  15. Dampier, J.E.E.; Luckai, N.; Bell, F.W.; Towill, W.D. Do tree-level monocultures develop following Canadian boreal silviculture? Tree-level diversity tested using a new method. Biodivers. Conserv. 2007, 16, 2933–2948. [Google Scholar] [CrossRef]
  16. Zheng, J.; Blanco, J.A.; Wei, X.; Liu, C. Sustainable management of Metasequoia glyptostroboides plantation forests in Shanghai. Forests 2018, 9, 64. [Google Scholar] [CrossRef] [Green Version]
  17. Fahlvik, N.; Ekö, P.-M.; Pettersson, N. Influence of precommercial thinning grade on branch diameter and crown ratio in Pinus sylvestris in southern Sweden. Scand. J. For. Res. 2005, 20, 243–251. [Google Scholar] [CrossRef]
  18. Sharma, M.; Zhang, S.Y. Stand density management diagram for jack pine stands in Eastern Canada. North J. Appl. For. 2007, 24, 22–29. [Google Scholar] [CrossRef] [Green Version]
  19. Briggs, D.G.; Smith, W.R. Effects of silvicultural practices on wood properties in conifers: A review. In Douglas-Fir: Stand Management for the Future; Hanley, D.P., Johnson, J.A., Eds.; College of Forest Resources, University of Washington: Seattle, WA, USA, 1985; No. 55; pp. 108–117. [Google Scholar]
  20. Zhang, S.Y.; Chauret, G.; Ren, H.Q.; Desjardins, R. Impact of initial spacing on plantation black spruce lumber grade yield, bending properties, and MSR yield. Wood Fiber Sci. 2002, 34, 460–475. [Google Scholar]
  21. Beniušiene, L.; Šilinskas, B.; Beniušis, R.; Aleinikovas, M.; Petrauskas, E.; Varnagiryte-Kabašinkiene, I. Scots pine stem parameters in sites with different stand densities in Lithuania. Forests 2020, 11, 716. [Google Scholar] [CrossRef]
  22. Wagner, R.G.; Buse, L.J.; Lautenschlager, R.A.; Bell, F.W.; Hollstedt, C.; Strobl, S.; Morneault, A.; Lewis, W.; Ter-Mikaelian, M.T. Vegetation Management Alternatives Program: 1994–1995 Annual Report; Ontario Ministry of Natural Resources, Ontario Forest Research Institute: Sault Ste. Marie, ON, Canada, 1995; p. 99. [Google Scholar]
  23. Pitt, D.G.; Morneault, A.E.; Bunce, P.; Bell, F.W. Five years of vegetation succession following vegetation management treatments in a jack pine ecosystem. North J. Appl. For. 2000, 17, 100–109. [Google Scholar] [CrossRef] [Green Version]
  24. Newmaster, S.G.; Bell, F.W. The effects of silvicultural disturbances on cryptogram diversity in the boreal-mixedwood forest. Can. J. For. Res. 2002, 32, 38–51. [Google Scholar] [CrossRef]
  25. Forest Co-op Growth and Yield Science Unit. Annual Report; The Forest Ecosystem Science Co-operative Inc.: Thunder Bay, ON, Canada, 2009. [Google Scholar]
  26. Sims, R.A.; Towill, W.D.; Baldwin, K.A.; Wickware, G.M. Field Guide to the Forest Ecosystem Classification for Northwestern Ontario; Canadian Forest Service, Ontario Ministry of Natural Resources: Sault Ste. Marie, ON, Canada, 1989; p. 191. ISBN 0777859904. [Google Scholar]
  27. Chambers, B.; Naylor, B.; Nieppola, J.; Merchant, B.; Uhlig, P. Field Guide to Forest Ecosystems of Central Ontario; SCSS FG-01; Ontario Ministry of Natural Resources, Southcentral Science Centre: North Bay, ON, Canada, 1997; p. 22. [Google Scholar]
  28. Bell, F.W.; Hunt, S.; Dacosta, J.; Sharma, M.; Larocque, G.R.; Winters, J.A.; Newmaster, S.G. Effects of silviculture intensity on plant diversity response patterns in young managed northern temperate and boreal forests. Écoscience 2014, 21, 327–339. [Google Scholar] [CrossRef]
  29. Honer, T.G.; Ker, M.F.; Alemdag, L.S. Metric Timber Tables for Commercial Tree Species of Central and Eastern Canada; Info. Rep. M-X-140; Environment Canada, Canadian Forest Service, Maritimes Forest Research Centre: Fredericton, NB, Canada, 1983; p. 139. [Google Scholar]
  30. Sharma, M. Comparing height-diameter relationships of boreal tree species grown in plantations and natural stands. For. Sci. 2016, 62, 70–77. [Google Scholar] [CrossRef]
  31. Sharma, M.; Parton, J. Height-diameter equations for boreal tree species in Ontario using a mixed-effects modeling approach. For. Ecol. Manag. 2007, 249, 187–198. [Google Scholar] [CrossRef]
  32. Sharma, M. Inside and outside bark volume models for jack pine (Pinus banksiana) and black spruce (Picea mariana) plantations in Ontario, Canada. For. Chron. 2019, 95, 50–57. [Google Scholar] [CrossRef]
  33. Sharma, M. Increasing volumetric prediction accuracy—An essential prerequisite for end-product forecasting in red pine. Forests 2020, 11, 1050. [Google Scholar] [CrossRef]
  34. Honer, T.G. Standard Volume Tables and Merchantable Conversion Factors for the Commercial Tree Species of Central and Eastern Canada; Info. Rep. FMR-X-5; Canadian Department of Forestry and Rural Development, Forest Management Research and Services Institute: Ottawa, ON, Canada, 1967; p. 21. [Google Scholar]
  35. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. Available online: https://www.R-project.org/ (accessed on 9 January 2021).
  36. Homagain, K.; Shahi, C.K.; Leitch, M.; Luckai, N.J.; Bell, F.W. Differences in extrinsic tree quality and value of fibre production following alternative vegetation management treatments in northwestern Ontario. For. Chron. 2011, 87, 251–259. [Google Scholar] [CrossRef] [Green Version]
  37. Thompson, I.D.; Guariguata, M.B.; Okabe, K.; Bahamondez, C.; Nasi, R.; Heymell, V.; Sabogal, C. An operational framework for defining and monitoring forest degradation. Ecol. Soc. 2013, 18, 20. [Google Scholar] [CrossRef]
  38. Condit, R.; Hubbell, S.P.; Lafrankie, J.V.; Sukumar, R.; Manokaran, N.; Foster, R.B.; Ashton, P.S. Species-area and species-individual relationships for tropical trees: A comparison of three 50-ha plots. J. Ecol. 1996, 84, 549–562. [Google Scholar] [CrossRef] [Green Version]
  39. Yang, T.-R.; Lam, T.Y.; Su, S.-H. A simulation study on the effects of plot size and shape on sampling plant species composition for biodiversity management. J. Sustain. For. 2018, 38, 116–129. [Google Scholar] [CrossRef]
  40. Pausas, J.G.; Austin, M.P. Patterns of plant species richness in relation to different environments: An appraisal. J. Veg. Sci. 2001, 12, 153–166. [Google Scholar] [CrossRef]
  41. Yang, Z.; Hautier, Y.; Borer, E.T.; Zhang, C.; Du, G. Abundance- and functional-based mechanisms of plant diversity loss with fertilization in the presence and absence of herbivores. Oecologia 2015, 179, 261–270. [Google Scholar] [CrossRef]
  42. Lichter, J. Primary succession and forest development on coastal lake Michigan sand dunes. Ecol. Monogr. 1998, 68, 487–510. [Google Scholar] [CrossRef]
  43. Kemperman, J. Sucker-Root Relationships in Aspen; Forest Research Note; Paper 4773; no. 12; Ministry of Natural Resources: ON, Canada, 1978. [Google Scholar]
  44. Farmer, R.J. Aspen root sucker formation and apical dominance. J. For. Sci. 1962, 8, 403–410. [Google Scholar]
  45. Frey, B.R.; Lieffers, V.J.; Landhäusser, S.M.; Comeau, P.G.; Greenway, K.J. An analysis of sucker regeneration of trembling aspen. Can. J. For. Res. 2003, 33, 1169–1179. [Google Scholar] [CrossRef]
  46. Perala, D.A. Populus tremuloides Michx.-Quaking Aspen. In Silvics of North America: Volume 2, Hardwoods; Agriculture Handbook 654; Burns, R.M., Honkala, B.H., Eds.; United States Department of Agriculture, Forest Service: Washington, DC, USA, 1990; pp. 555–569. [Google Scholar]
  47. Wang, G.G. Early regeneration and growth dynamics of Populus tremuloides suckers in relation to fire severity. Can. J. For. Res. 2003, 33, 1998–2006. [Google Scholar] [CrossRef]
  48. McCavour, M.J.; Paré, D.; Messier, C.; Thiffault, N.; Thiffault, E. The role of aggregated forest harvest residue in soil fertility, plant growth, and pollination services. Soil Sci. Soc. Am. J. 2014, 78, S196–S207. [Google Scholar] [CrossRef]
  49. Pickett, S.T.A.; Kempf, J.S. Branching patterns in forest shrubs and understory trees in relation to habitat. New Phytol. 1980, 86, 219–228. [Google Scholar] [CrossRef]
  50. Gregg, B.M.; Dougherty, P.M.; Hennessey, T.C. Growth and wood quality factors in relation to stand density and climatic factors. Can. J. For. Res. 2011, 18, 851–858. [Google Scholar] [CrossRef]
  51. Man, R.; Rice, J.A.; MacDonald, G.B. Early effects of pre- and post-harvest herbicide application and partial cutting in regenerating aspen—Jack pine mixtures in northeastern Ontario. Can. J. For. Res. 2011, 41, 1082–1090. [Google Scholar] [CrossRef]
  52. Fu, S.; Chen, H.Y.H.; Bell, F.W.; Sharma, M.; Delaney, J.R.; Peterson, G. Effects of timing of glyphosate application on jack pine, black spruce, and white spruce plantations in northern Manitoba. For. Chron. 2008, 84, 37–45. [Google Scholar] [CrossRef] [Green Version]
  53. Homagain, K.; Shahi, C.K.; Luckai, N.J.; Leitch, M.; Bell, F.W. Benefit-cost analysis of vegetation management alternatives: An Ontario case study. For. Chron. 2011, 87, 260–273. [Google Scholar] [CrossRef] [Green Version]
  54. Jobidon, R. Stump height effects on sprouting of mountain maple, paper birch and pin cherry—10 year results. For Chron. 1997, 73, 590–595. [Google Scholar] [CrossRef]
  55. Bell, F.W.; Morneault, A.E.; Pickering, S.M. Response of immature trembling aspen to season and height of cut. North J. Appl. For. 1999, 16, 108–114. [Google Scholar] [CrossRef] [Green Version]
  56. Bell, F.W.; Pitt, D.G. Seasonal susceptibility of boreal plants: Red raspberry phenology as a bioindicator of optimum within-season timing of glyphosate applications. For. Chron. 2007, 83, 733–741. [Google Scholar] [CrossRef] [Green Version]
  57. Wagner, R.G.; Robinson, A.P. Critical period of interspecific competition for four northern conifers: 10-year growth response and associated vegetation dynamics. Can. J. For. Res. 2006, 36, 2474–2485. [Google Scholar] [CrossRef] [Green Version]
  58. Ilisson, T.; Chen, H.Y.H. The direct regeneration hypothesis in northern forests. J. Veg. Sci. 2009, 20, 735–744. [Google Scholar] [CrossRef]
  59. Taylor, A.R.; Gao, B.; Chen, H.Y.H. The effect of species diversity on tree growth varies during forest succession in the boreal forest of central Canada. For. Ecol. Manag. 2020, 455, 117641. [Google Scholar] [CrossRef]
  60. Man, R.; Bell, F.W. Temporal changes of understory plant community in response to pre- and post-harvesting herbicide treatments and partial cutting in aspen-dominated boreal mixedwood stands. Eur. J. For. Res. 2018, 137, 337–348. [Google Scholar] [CrossRef]
  61. Kumar, P.; Chen, H.Y.H.; Thomas, S.C.; Shahi, C. Linking resource availability and heterogeneity to understorey species diversity through succession in boreal forest of Canada. J. Ecol. 2018, 106, 1266–1276. [Google Scholar] [CrossRef]
  62. Zhang, Y.; Chen, H.Y.H.; Taylor, A.R. Positive species diversity and above-ground biomass relationships are ubiquitous across forest strata despite interference from overstorey trees. Funct. Ecol. 2017, 31, 419–426. [Google Scholar] [CrossRef]
  63. Hedwall, P.; Holmström, E.; Lindbladh, M.; Felton, A. Concealed by darkness: How stand density can override the biodiversity benefits of mixed forests. Ecosphere 2019, 10, e02835. [Google Scholar] [CrossRef] [Green Version]
  64. Wagner, R.; Flynn, J.; Mertz, C.; Slovic, P.; Gregory, R. Acceptable practices in Ontario’s forests: Differences between the public and forestry professionals. New For. (Dordr.) 1998, 16, 139–154. [Google Scholar]
  65. Jobidon, R.; Cyr, G.; Thiffault, N. Plant species diversity and composition along an experimental gradient of northern hardwood abundance in Picea mariana plantations. For. Ecol. Manag. 2004, 198, 209–221. [Google Scholar] [CrossRef]
  66. Connell, J. Diversity in tropical rainforests and coral reefs. Science 1978, 199, 1302–1310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Collins, S.; Glenn, S.; Gibson, D. Experimental analysis of intermediate disturbance and initial floristic composition: Decoupling cause effect. Ecology 1995, 76, 486–492. [Google Scholar] [CrossRef]
  68. Haeussler, S.; Bedford, L.; Leduc, A.; Bergeron, Y.; Kranabetter, M. Silvicultural disturbance severity and plant communities of the southern Canadian boreal forest. Silva Fenn. 2002, 76, 307–327. [Google Scholar] [CrossRef] [Green Version]
  69. Magnussen, S.; Smith, V.G.; Yeatman, C.W. Tree size, biomass, and volume growth of twelve 34-year-old Ontario jack pine provenances. Can. J. For. Res. 1985, 15, 1129–1136. [Google Scholar] [CrossRef]
  70. ESSA Technologies Ltd. FVSOntario User Guide: Forest Stand Growth Projection Modelling System (Version 3.08); ESSA Technologies Ltd.: Vancouver, BC, Canada, 2008. [Google Scholar]
  71. Guo, J.; Wang, J.R. Comparison of height growth and growth intercept models of jack pine plantations and natural stands in northern Ontario. Can. J. For. Res. 2006, 36, 2179–2188. [Google Scholar] [CrossRef]
  72. Zhang, L.; Peng, C.; Dang, Q. Individual-tree basal area growth models for jack pine and black spruce in northern Ontario. For. Chron. 2004, 80, 366–374. [Google Scholar] [CrossRef] [Green Version]
  73. Morris, D.; Reid, D.; Kwiaton, M.; Hunt, S.; Gordon, A. Comparing growth patterns of jack pine and black spruce in mixed natural stands and plantations. Écoscience 2014, 21, 1–10. [Google Scholar] [CrossRef]
  74. Barry, K.E.; Mommer, L.; van Ruijven, J.; Wirth, C.; Wright, A.J.; Bai, Y.; Connolly, J.; De Deyn, G.B.; de Kroon, H.; Isbell, F.; et al. The future of complementarity: Disentangling causes from consequences. Trends Ecol. Evol. 2019, 34, 167–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Pretzch, H.; Schütze, G. Crown allometry and growing space efficiency of Norway spruce (Picea abies [L.] Karst.) and European beech (Fagus sylvatica L.) in pure and mixed stands. Plant Biol. 2005, 7, 628–639. [Google Scholar] [CrossRef] [Green Version]
  76. Kelty, M. Comparative productivity of monocultures and mixed-species stands. In The Ecology and Silviculture of Mixed-Species Forests; Kelty, M.J., Larson, B.C., Oliver, C.D., Eds.; Springer: Dordrecht, The Netherlands, 1992; Volume 40, pp. 125–141. ISBN 978-90-481-4135-7. [Google Scholar]
  77. Zhang, Y.; Chen, H.Y.H.; Reich, P.B. Forest productivity increases with evenness, species richness and trait variation: A global meta-analysis. J. Ecol. 2012, 100, 742–749. [Google Scholar] [CrossRef]
  78. Searle, E.B.; Chen, H.Y.Y. Complementarity effects are strengthened by competition intensity and global environmental change in the central boreal forests of Canada. Ecol. Lett. 2020, 23, 79–87. [Google Scholar] [CrossRef] [PubMed]
  79. Reich, P.B.; Tilman, D.; Isbell, F.; Mueller, K.; Hobby, S.E.; Flynn, D.F.B.; Eisenhauer, N. Impacts of biodiversity loss escalate through time as redundancy fades. Science 2012, 336, 589–592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Chen, X.; Chen, H.Y.H.; Searle, E.B.; Chen, C.; Reich, P.B. Negative to positive shifts in diversity effects on soil nitrogen over time. Nat. Sustain. 2020, 1–8. [Google Scholar] [CrossRef]
  81. Jucker, T.; Koricheva, J.; Finér, L.; Bouriaud, O.; Iacopetti, G.; Coomes, D. Good things take time—Diversity effects on tree growth shift from negative to positive during stand development in boreal forests. J. Ecol. 2020, 108, 2198–2211. [Google Scholar] [CrossRef]
  82. Brassard, B.W.; Chen, H.Y.H.; Cavard, X.; Laganière, J.; Reich, P.B.; Bergeron, Y.; Pare, D.; Yuan, Z. Tree species diversity increases fine root productivity through increased soil volume filling. J. Ecol. 2013, 101, 210–219. [Google Scholar] [CrossRef]
  83. Urli, M.; Thiffault, N.; Houle, D.; Gauthier, S.; Bergeron, Y. Role of green alder in boreal conifer growth: Competitor or facilitator? FACETS 2020, 5, 166–181. [Google Scholar] [CrossRef]
  84. Griess, V.; Knoke, T. Growth performance, windthrow, and insects: Meta-analyses of parameters influencing performance of mixed-species stands in boreal and northern temperate biomes. Can. J. For. Res. 2011, 41, 1141–1158. [Google Scholar] [CrossRef]
  85. Wagner, R.G. Competition and critical-period thresholds for vegetation management decisions in young conifer stands. For. Chron. 2000, 76, 961–968. [Google Scholar] [CrossRef] [Green Version]
  86. Kenkel, N.; Walker, D.; Watson, P.; Caners, R.; Lastra, R. Vegetation dynamics in boreal forest ecosystems. Coenoses 1997, 12, 97–108. [Google Scholar]
  87. Logan, K.T. Growth of Tree Seedlings as Affected by Light Intensity. II. Red Pine, White Pine, Jack Pine, and Eastern Larch; Publication 1160; Canadian Department of Forestry: Ottawa, ON, Canada, 1966; p. 19. [Google Scholar]
  88. Groot, A.; Adhikary, S.; Sharma, M.; Luckai, N.; Bell, F.; Larocque, G. Effect of species composition on the production rate and efficiency of young Picea glauca–Populus tremuloides forests. For. Ecol. Manag. 2014, 315, 1–11. [Google Scholar] [CrossRef]
  89. Jobidon, R.; Trottier, F.; Charette, L. Dégagement chimique ou manuel de plantations d’épinette noire? Étude de cas dans le domaine de la sapinière à bouleau blanc au Québec. For. Chron. 1999, 76, 973–979. [Google Scholar] [CrossRef] [Green Version]
  90. Bella, I.E.; DeFranceschi, J.P. Spacing Effects 15 Years after Planting Three Conifers in Manitoba; Information Report NOR-X-223; Canadian Forestry Service, Northern Forest Research Centre: Edmonton, AB, Canada, 1980; p. 10. [Google Scholar]
  91. Bella, I.E. Spacing Effects 20 Years after Planting Three Conifers in Manitoba; Forest Management Note 39; Canadian Forestry Service, Northern Forestry Centre: Edmonton, AB, Canada, 1987; p. 12. [Google Scholar]
  92. Cremer, K.W.; Borough, C.J.; McKinnell, F.H.; Carter, P.R. Effects of stocking and thinning on wind damage in plantations. N. Z. J. For. 1982, 12, 244–268. [Google Scholar]
  93. Tong, Q.J.; Zhang, S.Y. Impact of initial spacing and precommercial thinning on jack pine tree growth and stem quality. For. Chron. 2005, 81, 418–427. [Google Scholar] [CrossRef] [Green Version]
  94. Archibald, D.J.; Bowling, C. Jack Pine Density Management Diagram for Boreal Ontario; TN-005 (NWST TN-31); Ontario Ministry of Natural Resources, Northeast and Northwest Science and Technology: Thunder Bay, ON, Canada, 1995; p. 20. [Google Scholar]
  95. Lu, P.; Bell, F.; Charrette, P.; Thompson, M. Performance of jack pine (Pinus banksiana) rooted cuttings from proliferated dwarf shoots versus seedlings 8 years after planting. Can. J. For. Res. 2012, 42, 1–6. [Google Scholar] [CrossRef]
  96. Lautenschlager, R.; Sullivan, T. Effects of herbicide treatments on biotic components in regenerating northern forests. For. Chron. 2002, 78, 695–731. [Google Scholar] [CrossRef] [Green Version]
  97. Swift, K.I.; Bell, F.W. What are the environmental consequences of using silviculturally effective vegetation management alternatives? For. Chron. 2011, 87, 201–216. [Google Scholar] [CrossRef]
  98. Venier, L.A.; Thompson, I.D.; Fleming, R.; Malcolm, J.; Aubin, I.; Trofymow, J.A.; Langor, D.; Sturrock, R.; Patry, C.; Outerbridge, R.O.; et al. Effects of natural resource development on the terrestrial biodiversity of Canadian boreal forests. Environ. Rev. 2014, 22, 457–490. [Google Scholar] [CrossRef]
  99. Bell, F.W.; Thiffault, N.; Szuba, K.; Luckai, N.J.; Stinson, A. Synthesis of silviculture options, costs, and consequences of alternative vegetation management practices relevant to boreal and temperate conifer forests: Introduction. For. Chron. 2011, 87, 155–160. [Google Scholar] [CrossRef] [Green Version]
  100. Thiffault, N.; Roy, V. Living without herbicides in Québec (Canada): Historical context, current strategy, research and challenges in forest vegetation management. Eur. J. For. Res. 2011, 130, 117–133. [Google Scholar] [CrossRef]
  101. Wyatt, S.; Rousseau, M.H.; Nadeau, S.; Thiffault, N.; Guay, L. Social concerns, risk and the acceptability of forest vegetation management alternatives: Insights for managers. For. Chron. 2011, 87, 274–289. [Google Scholar] [CrossRef]
  102. Moffat, K.; Lacey, J.; Zhang, A.; Leipold, S. The social licence to operate: A critical review. Forestry 2015, 89, 477–488. [Google Scholar] [CrossRef] [Green Version]
  103. Kayahara, G.J.; Armstrong, C.L. Understanding First Nations rights and perspectives on the use of herbicides in forestry: A case study from northeastern Ontario. For. Chron. 2015, 91, 126–140. [Google Scholar] [CrossRef] [Green Version]
  104. Schneider, R.; Zhang, S.Y.; Swift, D.E.; Bgin, J.; Lussier, J.-M. Predicting selected wood properties of jack pine following commercial thinning. Can. J. For. Res. 2008, 38, 2030–2043. [Google Scholar] [CrossRef]
  105. Pielou, E.C. Ecological Diversity; John Wiley and Sons: New York, NY, USA, 1975; p. 165. [Google Scholar] [CrossRef]
  106. Simpson, E.H. Measurement of diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  107. Heip, C. A new index measuring evenness. J. Mar. Biolog. Assoc. UK 1974, 54, 555–557. [Google Scholar] [CrossRef]
Figure 1. Radar diagrams proposed as a visual illustration of the effects of four release treatments on (a) plant biodiversity indicators, (b) stand-scale growth and yield indicators, and (c) individual crop tree growth indicators for two jack pine sites located in Ontario, Canada: Bending Lake and E.B. Eddy. Interior colored lines represent scoring for each criterion separated by treatment (aerial spray, motor-manual brushsaw, complete removal, and no competition control) 25 growing seasons after treatments applied. Values were calculated as percentage deviation from the baseline (ideal or maximum site conditions).
Figure 1. Radar diagrams proposed as a visual illustration of the effects of four release treatments on (a) plant biodiversity indicators, (b) stand-scale growth and yield indicators, and (c) individual crop tree growth indicators for two jack pine sites located in Ontario, Canada: Bending Lake and E.B. Eddy. Interior colored lines represent scoring for each criterion separated by treatment (aerial spray, motor-manual brushsaw, complete removal, and no competition control) 25 growing seasons after treatments applied. Values were calculated as percentage deviation from the baseline (ideal or maximum site conditions).
Forests 12 00370 g001
Table 1. Description of jack pine release study areas in northwestern Ontario, Canada.
Table 1. Description of jack pine release study areas in northwestern Ontario, Canada.
SiteBlockGeologyPreharvest Stand ConditionsHarvestSite PreparationPlanting
Bending Lake1–4 (48°57′ N, 92°02′ W)Flat to gently rolling; 1–10% bedrock outcrops; rapidly drained, coarse loamy to fine sandV-17 : jack pine mixedwood/shrub rich forestClearcut; tree length, conventional cut and skid, 1986–1987Mechanical; heavy drags (barrels and chains), 1987Block 1: Jack pine; Cerkon shelter cone seeding, 3000 seedlings ha−1 (1.8 × 1.8 m spacing), 1988
Blocks 2–4: Jack pine; overwintered Spencer-Lemaire 6 container stock; 3000 seedlings ha−1 (1.8 × 1.8 m spacing), 1988
E.B. Eddy1 (46°47′ N, 82°08′ W)Glacial outwash plain; fine sands to silty loams (ES13.1) *FRI 1974 ±: Pw5Pj3Bw2
Site Class 2; age 150; height 24 m; stocking 0.7
Clearcut; tree-length, conventional cut and skid, 1988 *Young’s teeth, 1989 *Jack pine; current 408 paper pot, planted with Potti Putki, June 27–29, 1989
2 (46°47′ N, 82°10′ W)Glacial outwash plain; silty loam (ES15.1) FRI 1974±: Pj7Sb2Po1
Site Class 2; age 85; height 20 m; stocking 0.8
Clearcut; tree-length, conventional cut and skid, 1989 Young’s teeth, 1989 Jack pine; overwintered FH408 paper pot, planted with Potti Putki, May 13–16, 1991
3 (46°46′ N, 82°04′ W)Glacial outwash plain; silty, very fine sands to silty loams (ES15.1) FRI 1974±: Pj6Sb2Po2
Site Class 2; age 90;
height 21 m; stocking 0.6
Clearcut; tree-length, conventional cut and skid, 1989 Young’s teeth, 1989 Jack pine; current 165 Jiffy pot, planted with Potti Putki, June 14, 1991
[26]. * [27]. ± Forest Resource Inventory stand description: Bw = white birch, Sb = black spruce, Pj = jack pine, Po = poplar, Pw = white pine. Superscripts refer to the proportion of total stocking (out of 10) occupied by a species. [23].
Table 2. Estimated marginal means (± SE) for biodiversity indicators (n = 85) grouped by site x treatment variables after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) in each site. Asterisk (*) indicates statistically significant results from pairwise comparisons between sites.
Table 2. Estimated marginal means (± SE) for biodiversity indicators (n = 85) grouped by site x treatment variables after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) in each site. Asterisk (*) indicates statistically significant results from pairwise comparisons between sites.
SiteTreatmentSoftwood Richness (Number of Species)Hardwood Richness (Number of Species)Shrub Layer Richness (Number of Species)Herb Layer Richness (Number of Species)Tree Abundance (Percent Cover)Shrub Layer Abundance (Percent Cover)Herb layer Abundance (Percent Cover)Understory Diversity (H’)Understory Dominance (D)Understory Evenness (EHeip)
Bending LakeControl1.8 * (0.2)2.4 (0.3)9.0 (0.5)11.5a (0.6)52.0 (3.7)41.9 (5.0)47.2 (4.2)1.48 (0.18)0.34 (0.06)0.41 * (0.05)
Brushsaw1.8 * (0.2)3.0 (0.3)10.3 (0.3)11.8a (0.6)35.4 (3.7) *48.8 (6.8)34.2 (4.8)1.29 (0.23)0.36 (0.08)0.31 (0.04)
Aerial spray2.0 * (0.2)3.3 (0.2)9.8 (0.4)14.2b (0.8)39.0 (4.0) *52.2 * (7.1)49.5 (6.5)1.52 (0.22)0.38 (0.08)0.35 (0.07)
Complete removal2.6 (0.3)3.0 (0.4)8.7 (0.5)15.1b (0.7)43.2 (4.0) *54.1 * (7.8)61.8 * (10.5)1.59 (0.22)0.40 (0.07)0.39 (0.05)
E.B. EddyControl3.7 (0.3)2.0 (0.4)8.2 (0.4)13.7 (0.5)48.3 (2.9)43.2 (4.3)51.1 (10.0)1.23 (0.28)0.48 (0.12)0.22 (0.04)
Brushsaw4.2 (0.3)3.0 (0.2)8.9 (0.9)12.7 (1.3)54.4 (3.8)41.2 (5.6)33.0 (7.0)1.45 (0.30)0.41 (0.11)0.31 (0.05)
Aerial spray3.1 (0.5)3.3 (0.4)9.2 (0.6)14.0 (0.8)56.4 (2.6)34.5 (6.4)44.2 (7.7)1.44 (0.29)0.40 (0.10)0.32 (0.05)
Complete removal3.0 (0.4)4.0 (0.2)8.2 (0.9)13.7 (0.7)54.3 (3.8)32.1 (5.7)28.5 (3.8)1.56 (0.23)0.30 (0.08)0.39 (0.05)
* p values for fixed effects provided in Table S3.
Table 3. Estimated marginal means (± SE) for stand growth and yield indicators grouped by site x treatment combinations after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) within each site. Asterisk (*) indicates statistically significant differences from pairwise comparisons between sites.
Table 3. Estimated marginal means (± SE) for stand growth and yield indicators grouped by site x treatment combinations after vegetation management treatments applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) within each site. Asterisk (*) indicates statistically significant differences from pairwise comparisons between sites.
SiteTreatmentGross Total Stem Density (Stems ha−1) Conifer Density (Stems ha−1)Hardwood Density (Stems ha−1)Conifer Basal Area (m2 ha−1)Conifer Inside Bark Volume (m3 ha−1)Hardwood Basal Area (m2 ha−1)Hardwood Inside Bark Volume (m3 ha−1)
(n = 151)(n = 81)(n = 70)(n = 81)(n = 81)(n = 70)(n = 70)
Bending LakeControl3518c (228)1233a (121)2285b * (287)11.8a (1.5)63.6a (8.1)10.6b * (1.3)57.0b * (6.9)
Brushsaw2910bc * (183)2498c * (157)452a * (73)27.7b (0.7)148b (5.2)0.69a (0.2)3.53a (0.9)
Aerial spray2402ab (177)2236bc * (145)180a * (53)28.4b * (0.8)160b * (4.7)0.50a (0.2)2.90a (1.2)
Complete removal2045a * (102)1794b (71)430a (166)28.2b (0.8)154b (5.5)0.60a (0.2)3.27a (1.0)
E.B. EddyControl4011a (344)958a (163)3053b (353)10.8a (2.1)63.0a (10.9)13.8b (2.2)89.1b (19.0)
Brushsaw2331b (224)1842b (114)989a (219)26.9b (0.9)154b (7.1)2.38a (0.7)13.6a (4.1)
Aerial spray2561b (149)1736b (108)850a (143)22.7b (1.1)128b (4.9)1.65a (0.3)8.24a (1.5)
Complete removal2725b (177)1925b (89)550a (173)27.3b (0.9)152b (8.8)1.31a (0.4) 6.90a (2.3)
* p values for fixed effects provided in Table S5.
Table 4. Analysis of variance results with orthogonal contrasts for individual jack pine tree survival (n = 560), stem quality (n = 360), and branching (n = 360) after vegetation management treatments (C = control, BS = motor-manual brushsaw, AS = aerial spray, CR = complete removal) applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada.
Table 4. Analysis of variance results with orthogonal contrasts for individual jack pine tree survival (n = 560), stem quality (n = 360), and branching (n = 360) after vegetation management treatments (C = control, BS = motor-manual brushsaw, AS = aerial spray, CR = complete removal) applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada.
SourcedfSurvivalStem QualityBranching
F p-ValueFp-ValueFp-Value
Site *15.80.01630.41<0.0014.120.048
Treatment *317.71<0.0011.140.3348.67<0.001
CR vs. AS11.6510.099NANA2.390.017
CR vs. BS11.0680.286NANA3.310.001
AS vs. BS11.790.074NANA0.780.435
C vs. treatments17.06<0.001NANA3.8<0.001
* Means (±SE), separated by site and treatment, provided in Table S6.
Table 5. Estimated marginal means (±SE) for individual jack pine growth characteristics (n = 360), grouped by site x treatment variables after vegetation management applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) within each site. Asterisk (*) indicates statistically significant differences from pairwise comparisons between sites.
Table 5. Estimated marginal means (±SE) for individual jack pine growth characteristics (n = 360), grouped by site x treatment variables after vegetation management applied on two sites (Bending Lake and E.B. Eddy) in Ontario, Canada. Different lowercase letters indicate statistically significant differences from Tukey’s honest significant difference (HSD) test (pairwise comparisons between each treatment and control) within each site. Asterisk (*) indicates statistically significant differences from pairwise comparisons between sites.
SiteTreatmentTotal Height (m) Quadratic Mean Diameter (cm)Height to Live Crown (m)Branch Diameter (mm)Rings (Per cm)Merchantable Volume (m3 ha−1)
Bending LakeControl9.5a * (0.4)11.0a (0.4)6.1a (1.1)11.5a (2.0)2.2b * (0.1)33.3a * (3.3)
Brushsaw11.4b * (0.2)11.9b (0.3)7.2b (0.8) 16.6b * (2.4)2.0ab * (0.1)97.4b * (6.1)
Aerial spray12.7c (0.1)12.7b (0.3)7.8b (1.2)38.9d * (2.5)1.7a * (0.1)114.1c * (7.4)
Complete removal11.9bc (0.2)14.2c (0.2)7.3b (1.1)22.8c (2.8)1.6a * (0.1)129.2d (5.1)
E.B. EddyControl10.8a (0.4)13.5ab (0.3)7.0a (0.2)13.3a (1.3)1.7b (0.1)32.6a (5.5)
Brushsaw12.4b (0.2)13.4ab (0.5)7.0a (0.2)19.8b (0.7)1.4a (0.1)93.9b (5.8)
Aerial spray12.9b (0.2)12.0a (0.3)7.6b (0.2) 23.5c (0.9) 1.2a (0.1)118.1c (7.4)
Complete removal12.3b (0.3)12.9a (0.2)7.2ab (0.2)25.4c (1.0)1.3a (0.1)127.1c (4.8)
* p values for fixed effects provided in Table S7.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Deighton, H.D.; Bell, F.W.; Thiffault, N.; Searle, E.B.; Leitch, M.; Sharma, M.; Dacosta, J. Trade-Offs among Release Treatments in Jack Pine Plantations: Twenty-Five Year Responses. Forests 2021, 12, 370. https://doi.org/10.3390/f12030370

AMA Style

Deighton HD, Bell FW, Thiffault N, Searle EB, Leitch M, Sharma M, Dacosta J. Trade-Offs among Release Treatments in Jack Pine Plantations: Twenty-Five Year Responses. Forests. 2021; 12(3):370. https://doi.org/10.3390/f12030370

Chicago/Turabian Style

Deighton, Holly D., Frederick Wayne Bell, Nelson Thiffault, Eric B. Searle, Mathew Leitch, Mahadev Sharma, and Jennifer Dacosta. 2021. "Trade-Offs among Release Treatments in Jack Pine Plantations: Twenty-Five Year Responses" Forests 12, no. 3: 370. https://doi.org/10.3390/f12030370

APA Style

Deighton, H. D., Bell, F. W., Thiffault, N., Searle, E. B., Leitch, M., Sharma, M., & Dacosta, J. (2021). Trade-Offs among Release Treatments in Jack Pine Plantations: Twenty-Five Year Responses. Forests, 12(3), 370. https://doi.org/10.3390/f12030370

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop