The Impact of Variable Retention Harvesting on Growth and Carbon Sequestration of a Red Pine (Pinus resinosa Ait.) Plantation Forest in Southern Ontario, Canada

Abstract

As atmospheric carbon dioxide concentrations continue to rise and global temperatures increase, there is growing concern about the sustainability, health, and carbon sequestration potential of forest ecosystems. Variable retention harvesting (VRH) has been suggested to be a potential method to increase forest biodiversity, growth, and carbon (C) sequestration. A field trial was established in an 88-year-old red pine (Pinus resinosa Ait.) plantation in southern Ontario, Canada, using a completely randomized design to examine the response of tree productivity and other forest values to five harvesting treatments: 33% aggregate retention (33A), 55% aggregate retention (55A), 33% dispersed retention (33D), and 55% dispersed retention (55D) in comparison to an unharvested control (CN). In this study, we explored the impacts of VRH on aboveground stem radial growth and annual C increment. Standard dendrochronological methods and allometric equations were used to quantify tree- and stand-level treatment effects during a five-year pre-harvest (2009–2013) and post-harvest (2014–2018) period. Tree-level growth and C increment were increased by the dispersed retention pattern regardless of retention level. At the stand level, the total C increment was highest at greater retention levels and did not vary with retention pattern. These results suggest that the choice of retention level and pattern can have a large influence on management objectives as they relate to timber production, climate change adaptation, and/or climate change mitigation.

Introduction

Atmospheric carbon dioxide (CO₂) concentrations have increased from approximately 280 ppm in the mid-1700s to >411 ppm in 2020 due to fossil fuel combustion and other anthropogenic activities (Dlugokencky and Tans, 2018; Lindsey, 2018). Depending on global efforts to reduce emissions in the coming decades, atmospheric CO₂ concentrations are predicted to further increase to as high as 1000 ppm by 2100 (IPCC, 2014; Lindsey, 2018). Forests play an important role in the global carbon (C) cycle and can be managed to offset a proportion of these emissions by sequestering atmospheric C via photosynthesis into biomolecular C sinks (Unwin and Kriedemann, 2000; Pan et al., 2011). However, climate change and the associated projected increase in frequency, extent, and severity of natural disturbances represent a serious threat to the ecological functioning of forest ecosystems and their capacity to act as significant C sinks (Gauthier et al., 2014; Peterson et al., 2014).

A variety of silvicultural approaches can be used to manage forests to minimize the adverse effects of climate change stresses on productivity (Puettmann, 2011; Peterson et al., 2014; Williamson et al., 2019) and maintain or enhance the rate of C sequestration (Colombo et al., 2005; Peterson et al., 2014). Reduction of stand density by thinning and use of partial harvesting systems can reduce resource competition and increase the growth and C sequestration rates of residual trees and stands (Davis et al., 2009; Dwyer et al., 2010; Nunery and Keeton, 2010), as well as improve soil moisture availability, tree water status, and resilience to periodic drought (D’Amato et al., 2013; Sohn et al., 2016; Bradford and Bell, 2017). Variable retention harvesting (VRH) is a relatively new approach to partial harvesting used to balance economic with ecological objectives by creating a post-harvest forest structure that emulates that resulting from the natural disturbance regime (Franklin et al., 2007; Palik and D’Amato, 2019). Specifically, this approach retains a proportion of canopy trees and other biological structures at various densities and spatial patterns to enhance the structural complexity, biodiversity, and resilience of the resulting stand, while maintaining acceptable levels of productivity (Franklin et al., 2007; Palik and D’Amato, 2019). The primary considerations in application of VRH systems are the trees and structures to be retained, and the level, spatial pattern, and the period of retention best suited to the forest type being managed (Palik and D’Amato, 2019). Generally, VRH is used to create stands containing areas of both uniformly dispersed trees of different densities and aggregated patches of residual trees differing in size and shape (Palik and D’Amato, 2019; Franklin and Donato, 2020).

The effect of VRH on productivity and other forest values has been the focus of several, operational scale experiments established in North America over the past two decades (Palik et al., 2002, 2014; Aubry et al., 2009; Puettmann et al., 2016; Xing et al., 2018; Beese et al., 2019). The impact of VRH on productivity varies with forest type, tree species, the level and spatial pattern of tree retention, time after treatment, and other factors (Roberts and Harrington, 2008; Xing et al., 2018). In general, VRH has been shown to increase stem growth rates of individual trees as compared with unharvested areas (Roberts and Harrington, 2008; Powers et al., 2010; Xing et al., 2018). Tree growth rates and mortality are negatively related to retention level, but this relationship varies with species (Maguire et al., 2006; Xing et al., 2018). As well, trees located near the perimeter of retained patches have shown higher (Powers et al., 2009) or no difference in growth rate (Roberts and Harrington, 2008) compared to those in the interior. At the stand level, unharvested control plots have greater total biomass or volume growth compared to VRH treatments regardless of the level or pattern of retention (Maguire et al., 2006; Palik et al., 2014). Given the initial emphasis on VRH effects on tree regeneration and biodiversity, it remains unclear which level and pattern of tree retention best balances timber production with climate change adaptation and mitigation objectives for different forest types and site conditions in natural stands and plantations (Xing et al., 2018; Franklin and Donato, 2020).

Red pine (Pinus resinosa Ait.) plantations were established in many areas in the Great Lakes region of North America beginning in the early 1900s to control soil erosion on abandoned agricultural land and to rehabilitate large cutover and burned forests (Cayford and Bickerstaff, 1968; Johnson, 1995; Buckman et al., 2006). Red pine was favored for these efforts due to its commercial value, genetic uniformity, suitability to even-aged management, and ability to grow well in drought-prone, nutrient-poor soils (Stiell, 1959; Buckman et al., 2006). Many of the red pine plantations established on degraded agricultural land are located south of its natural range in highly fragmented, urbanized landscapes where the area of forested land has been dramatically reduced (Borczon, 1982; Elliott, 1998). These plantations are typically being managed to restore them to native forest types through traditional row or selection thinning methods coupled with regeneration of desired species by planting or release of advance reproduction (Parker et al., 2008; Abella, 2010).

In this study, we examined the impact of VRH on residual tree- and stand-level growth and C increment in an 88-year-old red plantation established beyond its natural range limit in southern Ontario, Canada. We used dendrochronological methods to examine VRH effects because tree rings provide fine spatial and temporal resolution in radial stem growth (Davis et al., 2009; Dye et al., 2016; Pompa-García et al., 2018) and are reliable proxy measures of tree- and stand-level C sequestration rates (Babst et al., 2014). The specific objectives of this study were to determine the effect of: (1) the level and spatial pattern of retention on post-harvest tree- and stand-level stem growth and C increment during the 5-year post-harvest period; and (2) tree location within aggregate patches of residual trees on tree stem growth and C increment.

Materials and Methods

Site Description

The study site is located within the Turkey Point tract of the St. Williams Conservation Reserve (SWCR) located in Norfolk County, in southern Ontario, Canada (42°42′16″ N, 80°21′29″ W) (Figure 1A). The SWCR consists of approximately 1034 ha of Crown Land managed in collaboration with the Ontario Ministry of Northern Development, Mines, Natural Resources and Forestry (NDMNRF). The primary management objective of SWCR is to increase biodiversity, protect species-at-risk, and restore these lands to their native ecological communities (SWCRMP (St. Williams Conservation Reserve Management Plan), 2007). This study was conducted in a 20.2 ha red pine plantation established in 1931. The area was most likely occupied by mixed oak and pine forests prior to settlement of the region. This plantation is part of the Turkey Point Flux Station and is associated with the Global Fluxnet and Global Water Futures research network (Peichl et al., 2010). The plantation was thinned around 1959–1960 when every fourth row of trees was removed. The plantation was thinned a second time in February 2014 to demonstrate and test the VRH approach with the long-term objective of restoring the area to a native forest type while enhancing biodiversity through the creation of a structurally complex residual tree canopy. A typical operational approach to the management of this plantation would be to reduce stand density according to guidelines for uniform shelterwood with the objective of initiating regeneration (Ontario Ministry of Natural Resources and Forestry (OMNRF), 2015).

FIGURE 1

A pre-harvest survey of the area conducted in 2011 using a series (n = 60) of 12.62 m radius plots indicated the canopy was almost entirely composed of dominant and codominant red pine, with a few scattered white pine (Pinus strobus L.) and black oak (Quercus velutina Lam.). The plantation was overstocked with an average total basal area (BA) of 38.4 ± 5.5 m²/ha and a density of 641.3 ± 103.2 trees/ha. Red pine averaged 28.6 ± 4.3 cm in diameter at 1.3 m (DBH) and ranged from 17.2 to 47.2 cm. Height of red pine averaged 23.8 ± 2.8 m. The subcanopy consisted of naturally established Prunus, Acer, and Quercus spp. of intermediate crown class and averaged about 70 trees/ha and 13.2 cm in DBH across the study area. The VRH study uses a completely randomized experimental design, with four replicates of five treatments (Figure 1B). Treatments were applied to 20 plots, each averaging 0.89 ha in area and separated by fire access roads (Figure 1B). The treatments include: 33% dispersed retention (33D), 55% dispersed retention (55D), 33% aggregate retention (33A), 55% aggregate retention (55A), and an unharvested control (CN). The dispersed treatments were created according to provincial guidelines for using the shelterwood silvicultural system for regeneration of red and white pine forests (Ontario Ministry of Natural Resources and Forestry (OMNRF), 2015). Tree marking of the dispersed treatments was aimed at uniform spacing of residual trees with overstory crown closure of 33 or 55%. The aggregate treatments were harvested to retain four (two large and two small) uncut, roughly circular patches over 33 or 55% of the plot area. Over the entire study area, 4,737 red pine with DBH of 16–48 cm were harvested, with ca. 74% trees removed having DBH of 24 to 32 cm. Pre- and post-harvest structural features and photographs of the VRH treatments are presented in Table 1 and Figures 2A–F, respectively.

FIGURE 2

Dendrochronological Methods

Cores were collected from at least 30 red pine trees within each of the 20 treatment plots at ∼1.3 m above ground (DBH) during the spring and summer of 2019 using a 5 mm increment borer (Haglöf, Sweden). For each tree, two cores were collected at approximately 90° from each other to account for variability within individual trees. In the 33A and 55A treatments, cores were collected from an equal number of trees from the interior and exterior of the patch to account for potential spatial variation in radial growth with location. As well, exterior trees sampled were evenly distributed around the patch perimeter to account for possible effects of aspect on radial growth (Roberts and Harrington, 2008; Powers et al., 2009). The DBH was measured for each tree sampled for dendrochronological assessment of annual C increment.

Tree cores were prepared for ring width measurement according to standardized dendrochronological methods (Stokes and Smiley, 1968; Speer, 2010). Cores were sanded with increasingly finer grades of sandpaper until the tree rings were clearly visible (Stokes and Smiley, 1968) and visually crossdated using the list method (Yamaguchi, 1991). Ring widths were measured to the nearest 0.001 mm using CooRecorder and CDendro (Larsson, 2018). Visual crossdating was verified using the COFECHA software program (Holmes, 1983, 1994) before assigning calendar dates to each tree ring (Speer, 2010).

Climate

The climate in southern Ontario is classified as Warm-Summer Humid Continental Climate (Kottek et al., 2006) and is characterized by cool winters and hot and humid summers (Crins et al., 2009). Climate data for the study area for the period of the tree-ring chronologies were obtained from Environment Canada climate stations located in Delhi, ON and Tillsonburg, ON located within 34 km from TP31 (Supplementary Table 1 and Supplementary Figure 1)¹. From 1937 to 2018, the mean annual temperature was 8.0°C, the mean annual maximum temperature was 13.2°C, and the mean annual minimum temperature was 2.9°C. Precipitation was evenly distributed throughout the year and the total annual precipitation averaged 962 mm.

Percent Growth Change

Percent growth change (%GC) over two discrete 5-year periods for each sampled tree was calculated using ring widths and the radial growth average technique (Nowacki and Abrams, 1997) as:

where M₁ is the mean pre-harvest radial growth (2009–2013) and M₂ is the mean post-harvest radial growth (2014–2018).

Annual Stem Biomass and C Increment

Tree ring derived DBH values and allometric equations were used to estimate annual stemwood biomass and C increment for each sample tree. It is important to note that estimating aboveground biomass from annual growth rings is a source of uncertainty (Martin-Benito et al., 2021). The allometric equations developed by Lambert et al. (2005) were selected for use because they were constructed from a large number of red pine trees (n = 371), primarily from Ontario, with an average DBH and height representative of the study plantation. Allometric equations using tree height and DBH measurements provide more accurate estimates of aboveground biomass, but annual height measurements of our sample trees were not available (Lambert et al., 2005; Cienciala et al., 2008). The allometric equation used to estimate the biomass of the individual tree stemwood is represented by the general form:

where y_wood is the dry biomass of the stemwood (kg), β_wood1 and β_wood2 are model parameters with coefficient estimates, D is DBH (cm), and e_wood is the error term (Lambert et al., 2005). For a mean prediction of biomass, e_wood was assumed to be 0, as it is assumed e_wood follows a normal distribution with a mean of 0. Based on the β_wood1 and β_wood2 estimates for red pine, the allometric equation used for this study is as follows (Lambert et al., 2005):

Calculating Tree-Level Stemwood C Increment

The DBH of each tree in each study year (2009–2018) was calculated as:

where DBH_k is the value for any k year (2009–2017), j is the final year of our assessment period (2018) and the upper limit of the summation, and RWi is the sum of the ring widths from year i to year j (Dai et al., 2013). The estimated DBH for each year and the allometric equation were used to estimate the total stemwood biomass increment for each year as:

where k represents the year of interest. Although C concentration may vary among species and plant parts (Lamlon and Savidge, 2005; King et al., 2007), C increment was estimated by multiplying, annual stemwood biomass by 0.5, assuming that approximately 50% of red pine stem dry mass is C (King et al., 2007; Dai et al., 2013). Using the generic value of 50% may not accurately estimate the absolute amount of C sequestered, but it is suitable for comparison of VRH treatment effects on C increment.

Stand-Level C Estimates

Stem density (trees/ha) of living red pine canopy trees for all treatment plots during the pre-harvest period was estimated using 2011 pre-harvest inventory data. The total number of red pine trees removed from each treatment plot was subtracted from its pre-harvest stand density to obtain initial (2014) post-harvest density. Annual changes in living red pine stand density during the post-harvest period were estimated using a subsample (n = 10–38) of canopy trees in each plot that were monitored for mortality (Supplementary Table 2). Because of the relatively small number of trees monitored, stand density for each plot was estimated using the mean percent mortality for a given treatment (Supplementary Table 2). Mean annual tree-level C increment (kg C/tree) for each treatment plot was multiplied by stem density to calculate pre- and post-harvest stand-level estimates of mean annual C increment (kg C/ha).

Statistical Analysis

One-way analysis of variance was used to quantify VRH treatment effects on %GC and tree- and stand-level annual C increment. Where significant treatment effects were exhibited, a post hoc Tukey test was performed to quantify differences among treatment means (computed with https://astatsa.com/OneWay_Anova_with_TukeyHSD/). Differences among treatments in annual C increment between pre- and post-harvest periods, and between interior and exterior trees within and among the aggregated treatments, were quantified using t-tests.

Results

Treatment Effects on Tree-Ring Growth Responses

The tree-ring chronologies for all 20 treatment plots spanned the period from 1936 to 2018 and exhibited comparable descriptive statistics (Table 2). The mean series intercorrelation was 0.646 and ranged from 0.551 to 0.736, indicating that ring growth of all sampled red pine trees was influenced by a common growth factor. The mean annual ring width growth increment was 1.64 mm and ranged from 1.47 to 1.89 mm. Autocorrelation was high (<0.871) for all treatment plots, with a mean value of 0.903 (Table 2), indicating that tree growth in a given year is influenced by growth in the previous year. The DBH of the 600 trees sampled averaged 31.1 ± 4.4 cm and ranged from 21.1 to 45.8 cm (Table 2).

Prior to harvest, ring-width index (RWI) was similar among treatments (Figure 3). Following harvest, RWI of all VRH treatments exhibited similar year-to-year high-frequency variability as the CN treatment. Beginning in 2016, the third growing season after harvest, the 33D and 55D treatments exhibited a trend toward higher RWI as compared with the 33A, 55A, and CN treatments (Figure 3). Treatment type had a significant (p < 0.05) effect on mean annual %GC over the post-harvest period, with values being significantly higher in the 33D than the 55A treatment (Figure 4A). In the aggregate treatments, mean annual %GC of the post-harvest period of exterior trees was significantly (p < 0.05) higher than interior trees in the 33A treatment but did not differ (p = 0.139) with tree position in the 55A treatment (Figure 4B). For both aggregate treatments combined, mean annual %GC was significantly (p ≤ 0.001) higher in the exterior trees.

FIGURE 3

FIGURE 4

Treatment Effects on C Increment

Treatment effects on pre- and post-harvest rates of mean annual C increment depended on the spatial scale at which this variable was expressed (Figure 5). At the individual tree level, mean annual C increment in the post-harvest period was significantly lower in the CN (p ≤ 0.001), 33A (p = 0.026), and 55A (p ≤ 0.001) treatments compared to the pre-harvest period (Figure 5A). Pre- and post-harvest mean annual C increment per tree did not differ in the 33D (p = 0.724) and 55D (p = 0.528) treatments. Mean annual C increment in the exterior trees in the post-harvest period was significantly higher than the interior trees in the 33A [t(38) = 4.623, p < 0.001], 55A [t(38) = 4.793, p ≤ 0.001] and both treatments combined [t(38) = 6.677, p < 0.001]. At the stand scale, mean annual C increment was significantly (p ≤ 0.001) lower in the post-harvest period for all five treatments (Figure 5B).

FIGURE 5

Treatments also differed in tree- and stand-level mean annual C increment within the post-harvest period (Figure 6). At the individual tree scale, the two dispersed retention treatments had significantly (p < 0.05) higher mean annual C increment than the CN and the two aggregate treatments (Figure 6A). For the aggregate treatments, tree location had a significant effect on growth, with exterior trees having significantly (p < 0.001) higher mean annual C increment than the interior trees (Figure 6B). By comparison, at the stand scale, mean annual C increment was significantly (p < 0.01) higher in the CN treatment than all VRH treatments except the 55D treatment (Figure 6C). Mean annual C increment did not differ between dispersed and aggregate treatments of the same retention level, but mean annual C increment was significantly (p < 0.05) higher in the 55D than the 33D and 33A treatments (p < 0.01).

FIGURE 6

Patterns of inter-annual variability in mean annual C increment rates at the individual tree- and stand-level during the 10-year pre- and post-harvest period suggest that the two dispersed retention treatments began to exhibit greater C uptake compared to the three other treatments beginning in 2016, the third growing season after harvest (Figure 7).

FIGURE 7

Discussion

Impact of Retention Level and Pattern on Stem Growth

Stand density reduction through thinning or partial harvest systems is employed in part to release residual canopy trees from neighborhood competition and increase tree growth rates (Gilmore and Palik, 2006). The pattern and intensity of tree removal determine the degree of release from competition. Annual variation in RWI from residual chronologies indicates that red pine radial growth exhibited a strong response to both interannual climatic variability and stand density changes over our 10-year study period. The shared annual variability in RWI is associated with climate-growth relationships established for red pine. Cooler, wetter growing season conditions favor higher growth in red pine growing in natural stands (Graumlich, 1993; Girardin et al., 2006) and plantations (Magruder et al., 2013; Ashiq and Anand, 2016). The gradual decline in RWI exhibited prior to harvest for all treatments is likely due to high inter-tree competition and growth suppression when stand density and BA were comparatively high. Post-harvest tree RWI varied with VRH treatments and indicated that stem radial growth rate of residual canopy trees depended primarily on the pattern of retention as it influenced the degree of reduction in neighborhood competition. The trend toward higher RWI in the dispersed treatments is most likely due to near-complete release from competition by adjacent trees in the uniformly spaced trees. In contrast, only the perimeter trees in the aggregate treatments experience a partial release from competition and more limited increase in resources (Chen et al., 1992; Sherich et al., 2007; Powers et al., 2009). The positive mean %GC values for the dispersed treatments and negative mean %GC for the aggregate treatments are consistent with trends in RWI. Further, although exterior trees had higher mean %GC than interior trees, mean values of %GC for exterior trees were negative and largely similar to that of the CN treatment. Interestingly, %GC of exterior trees was higher than interior trees in the 33A treatment and both aggregate treatments combined, but not the 55A treatment, suggesting that exterior trees in smaller patches with higher perimeter to area ratios may experience greater resource availability than in larger patches.

Impact of Retention Level and Pattern on Tree- and Stand-Level C Increment

Not surprisingly, the pattern of annual changes in C increment over our 10-year study period was consistent with that for RWI, with only trees in dispersed treatments exhibiting a gradual increase in C increment after harvest. The significantly lower post-harvest C increment for the control and aggregate treatments indicate that trees in residual patches respond similarly to the unharvested treatment plots, despite the comparatively higher radial growth observed for exterior trees that experienced a partial reduction in competition. These results are consistent with the findings for natural red pine stands (Powers et al., 2009, 2010) and other tree species (Deal and Tappeiner, 2002; McDonald and Urban, 2004; Sherich et al., 2007; Roberts and Harrington, 2008). In contrast to our findings, reduced growth of Douglas-fir in aggregate patches was attributed to poor growth response of exterior trees (Maguire et al., 2006).

The similarity in post-harvest mean C increment for trees in the dispersed treatments is likely a reflection of the delayed, gradual increase in radial growth following harvest that was not expressed until the third growing season after treatment. This response is probably due to so-called “thinning shock” that results from the abrupt exposure of trees to higher radiation, air temperature, wind speed, and a more evaporative atmospheric environment in thinned stands (Youngblood, 1991; Latham and Tappeiner, 2002; Bebber et al., 2004). This gradual, growth response is thought to be related to a brief period of physiological and morphological acclimation of trees to the new post-harvest microclimate when carbohydrate allocation to root growth may be favored over aboveground growth (Bladon et al., 2006; Powers et al., 2010). The comparatively greater variation observed in mean C increment in the dispersed treatments is probably associated with both heterogeneity in growth resources and differences in acclimation capacity among residual trees of different size and crown class (Roberts and Harrington, 2008). Trees of dominant and co-dominant crown classes with higher live crown ratio and leaf area show more immediate, positive growth response to thinning than subordinate, smaller diameter, less vigorous canopy trees of intermediate crown class (Horton and Bedell, 1960; Gilmore and Palik, 2006).

The two dispersed treatments exhibited significantly higher mean C increment in the post-harvest period than the control and aggregate treatments, while the C increment of the aggregate treatments did not differ from the unharvested control treatment, consistent with the findings of Roberts and Harrington (2008). This indicates that in our study, mean C increment is dependent on the pattern of retention, and not the level of retention. Maguire et al. (2006) also found higher stem growth in the dispersed compared to aggregate treatments due to increased access to resources and minimal crown overlap.

Expressed at the stand level, mean C increment was higher during the pre-harvest than the post-harvest period for all five treatments. Decreasing mean annual C increment in the control treatment reflects increasing stand basal area and tree growth suppression in unharvested plots, which is the rationale behind the use of frequent, periodic thinning to maximize growth and reduce mortality in managed red pine stands (Gilmore and Palik, 2006). Despite increased growth rates due to complete or partial release of residual canopy trees from competition, the greatly reduced stem density due to the thinning process reduced the total stand-level C increment in most VRH treatments (Maguire et al., 2006). This is not surprising given that thinning increases individual tree growth rates but rarely increases stand-level rates (Gilmore and Palik, 2006; Davis et al., 2009), since it takes time for these faster individual rates to compensate for lower densities.

Mean stand-level C increment in the post-harvest period was significantly higher in the unharvested control treatment than all other treatments, regardless of level or pattern of retention. Again, this is an expression of the much higher stand density in the control treatment despite lower tree C increment. In contrast to findings with tree-level C increment, stand-level C increment is driven by retention level. Stand-level C increment was significantly lower in the 33% compared with the 55% retention level regardless of pattern. In natural red pine stands, stand-level growth was higher in unharvested versus VRH treatments but was only marginally higher in dispersed treatments compared to aggregate treatments of the same retention level (Palik et al., 2014).

Conclusion

In this study of VRH effects on tree- and stand-level growth and C increment of red pine, we used tree rings and allometric equations to estimate stem radial growth and C increment, respectively. While using dendrochronological methods to estimate annual productivity has potential limitations (Martin-Benito et al., 2021), our comparatively large sample of randomly selected residual trees and our focus on a short, 10-year growth period, coupled with pre- and post-harvest tree and stand measurements, eliminated or reduced these sources of error (Babst et al., 2014; Dye et al., 2016). The results of this study indicate that VRH has a significant effect on tree- and stand-level stem growth and C increment of an 88-year-old red pine plantation that was significantly related to the level and spatial pattern of retention. Dispersed VRH treatments at both levels of retention were more effective in promoting post-harvest tree-level growth and C increment than aggregate treatments. In contrast, stand-level growth and C increment were highest for greater levels of retention, regardless of pattern. We recognize that stand C sequestration depends on above and below ground growth response of all forms of vegetation as well as treatment effects on C contained in dead biomass and mineral soil. Our focus on living red pine canopy trees clearly does not quantify the potential effects of VRH on total C accumulation.

The results of this study are applicable to management of red pine plantations in southern Ontario and other parts of the Great Lakes region where balancing timber production with climate change adaptation and/or mitigation is an objective (D’Amato et al., 2011). Where improved tree growth rates are desired, dispersed retention was the superior treatment. Where mitigation through enhanced stand-level C sequestration is a primary objective, a higher level of retention is desirable. In stands where adaptation and mitigation are of high priority, a mixture of dispersed and aggregate treatments of different retention levels may best meet both objectives.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1

  AbellaS. R. (2010). Thinning Pine Plantations to Reestablish Oak Openings Species in northwestern Ohio.Environmental Management46391–403. 10.1007/s00267-010-9538-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AshiqM. W.AnandM. (2016). Spatial and temporal variability in dendroclimatic growth response of red pine (Pinus resinosa Ait.) to climate in northern Ontario, Canada.Forest Ecology and Management372109–119. 10.1016/j.foreco.2016.04.008

  • CrossRef
  • Google Scholar

• 3

  AubryK. B.HalpernC. B.PetersonC. E. (2009). Variable-retention harvests in the Pacific Northwest: a review of short-term findings from the DEMO study.Forest Ecology and Management258398–408. 10.1016/j.foreco.2009.03.013

  • CrossRef
  • Google Scholar

• 4

  BabstF.AlexanderM. R.SzejnerP.BouriaudO.KlesseS.RodenJ.et al (2014). A tree-ring perspective on the terrestrial carbon cycle.Oecologia176307–322. 10.1007/s00442-014-3031-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  BebberD. P.ThomasS. C.ColeW. G.BalsillieD. (2004). Diameter increment in mature eastern white pine Pinus strobus L. following partial harvest of old-growth stands in Ontario.Canada. Trees1829–34. 10.1007/s00468-003-0274-y

  • CrossRef
  • Google Scholar

• 6

  BeeseW. J.DealJ.DunsworthB. G.MitchellS. J.PhilpottT. J. (2019). Two decades of variable retention in British Columbia: a review of its implementation and effectiveness for biodiversity conservation.Ecological Processes81–22. 10.1186/S13717-019-018-9

  • CrossRef
  • Google Scholar

• 7

  BladonK. D.SilinsU.LandhäusserS. M.LieffersV. J. (2006). Differential transpiration by three boreal tree species in response to increased evaporative demand and variable retention harvesting.Agricultural and Forest Meteorology138104–119. 10.1016/j.agformet.2006.03.015

  • CrossRef
  • Google Scholar

• 8

  BorczonE. L. (1982). Evergreen challenge: The Agreement Forest Story.Ontario, OT: Queen’s Printer.

  • Google Scholar

• 9

  BradfordJ. B.BellD. M. (2017). A window of opportunity for climate-change adaptation: easing tree mortality by reducing basal area.Frontiers in Ecology and Environment15:11–17. 10.1002/fee.1445

  • CrossRef
  • Google Scholar

• 10

  BuckmanR. E.BishawB.HansonT. J.BenfordF. A. (2006). Growth and yield of red pine in the Lake States. USDA For. Serv., Gen. Tech. Rep. NC-271.St. Paul, MN: U.S. Department of Agriculture, 114.

  • Google Scholar

• 11

  CayfordJ. H.BickerstaffA. (1968). Man-made forests in Canada. Dept. Fisheries and Forestry, Forestry Branch Pub. No. 1240.Ottawa, ON: Canadian Forest Service Publications, 68.

  • Google Scholar

• 12

  ChenJ.FranklinJ. F.SpiesT. A. (1992). Vegetation responses to edge environments in old-growth Douglas-fir forests.Ecological Applications2387–396. 10.2307/1941873

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  CiencialaE.ApltauerJ.ExnerováZ.TatarinovF. (2008). Biomass functions applicable to oak trees grown in Central-European forestry.Journal of Forest Science54109–120. 10.17221/2906-jfs

  • CrossRef
  • Google Scholar

• 14

  ColomboS. J.ParkerW. C.LuckaiN.DangQ.CaoT. (2005). The Effects of Forest Management on Carbon Storage in Ontario’s Forests. Ontario Ministry of Natural Resources, Ont. Min. Nat. Resour., Appl. Res. Develop. Br., Sault Ste. Marie, Ontario. Clim. Change Res. Rep. CCRR-03.Peterborough, ON: Ontario Ministry of Natural Resources, 114.

  • Google Scholar

• 15

  CrinsW. J.GrayP. A.UhligP. W. C.WesterM. C. (2009). Ecosystems of Ontario. Part 1: Ecozones and Ecoregions. Ontario, Ministry of Natural Resources, Inventory, Monitoring and Assessment Section. Technical report SIB TER IMA TR-01.Peterborough, ON: Ontario Ministry of Natural Resources.

  • Google Scholar

• 16

  DaiL.JiaJ.YuD.LewisB. J.ZhouL.ZhouW.et al (2013). Effects of climate change on biomass carbon sequestration in old-growth forest ecosystems on Changbai Mountain in Northeast China.Forest Ecology and Management300106–116. 10.1016/j.foreco.2012.06.046

  • CrossRef
  • Google Scholar

• 17

  D’AmatoA. W.BradfordJ. B.FraverS.PalikB. J. (2011). Forest management for mitigation and adaptation to climate change: Insights from long-term silviculture experiments.Forest Ecology and Management262803–816. 10.1016/j.foreco.2011.05.014

  • CrossRef
  • Google Scholar

• 18

  D’AmatoA. W.BradfordJ. B.FraverS.PalikB. J. (2013). Effects of thinning on drought vulnerability and climate response in northern temperate forest ecosystems.Ecological Applications231735–1742. 10.1890/13-0677.1

  • CrossRef
  • Google Scholar

• 19

  DavisS. C.HesslA. E.ScottC. J.AdamsM. B.ThomasR. B. (2009). Forest carbon sequestration changes in response to timber harvest.Forest Ecology and Management2582101–2109. 10.1016/j.foreco.2009.08.009

  • CrossRef
  • Google Scholar

• 20

  DealR. L.TappeinerJ. C. (2002). The effects of partial cutting on stand structure and growth of western hemlock–Sitka spruce stands in southeast Alaska.Forest Ecology and Management159173–186. 10.1016/S0378-1127(00)00727-1

  • CrossRef
  • Google Scholar

• 21

  DlugokenckyE.TansP. (2018). Trends in Atmospheric Carbon Dioxide. National Oceanic & Atmospheric Administration, Earth System Research Laboratory (NOAA/ESRL).Washington DC: NOAA.

  • Google Scholar

• 22

  DwyerJ. M.FenshamR.BuckleyY. M. (2010). Restoration thinning accelerates structural development and carbon sequestration in an endangered Australian ecosystem.Journal of Applied Ecology47681–691. 10.1111/j.1365-2664.2010.01775.x

  • CrossRef
  • Google Scholar

• 23

  DyeA.Barker PlotkinA.BishopD.PedersonN.PoulterB.HesslA. (2016). Comparing tree-ring and permanent plot estimates of aboveground net primary production in three eastern US forests.Ecosphere7e01454.

  • Google Scholar

• 24

  ElliottK. A. (1998). The forests of southern Ontario.The Forestry Chronicle74850–854. 10.5558/tfc74850-6

  • CrossRef
  • Google Scholar

• 25

  FranklinJ. F.DonatoD. C. (2020). Variable retention harvesting in the Douglas-fir region.Ecological Processes98. 10.1186/s13717-019-0205-5

  • CrossRef
  • Google Scholar

• 26

  FranklinJ. F.MitchellR. J.PalikB. J. (2007). Natural disturbance and stand development principles for ecological forestry. US Department of Agriculture, Forest Service, Northern Research Station, General Technical Report NRS-19.Newton Square, PA: US Department of Agriculture.

  • Google Scholar

• 27

  GauthierS.BernierP.BurtonP. J.EdwardsJ.IsaacK.IsabelN.et al (2014). Climate change vulnerability and adaptation in the managed Canadian boreal forest.Environmental Reviews22256–285. 10.1139/er-2013-0064

  • CrossRef
  • Google Scholar

• 28

  GilmoreD. W.PalikB. J. (2006). A Revised Manager’s Handbook for Red Pine in the North Central Region. US Department of Agriculture, Forest Service, General Technical Report NC-264.Washington, DC: US Department of Agriculture, 55.

  • Google Scholar

• 29

  GirardinM.-P.TardiffJ. C.FlanniganM. D.BergeronY. (2006). Synoptic-scale atmospheric circulation and boreal Canada summer drought variability of the past three centuries.J. Climate191922–1947. 10.1175/jcli3716.1

  • CrossRef
  • Google Scholar

• 30

  GraumlichL. J. (1993). Response of tree growth to climatic variation in the mixed conifer and deciduous forests of the upper Great Lakes region.Canadian Journal of Forest Research23133–143. 10.1139/x93-020

  • CrossRef
  • Google Scholar

• 31

  HolmesR. L. (1983). Computer-assisted quality control in tree-ring dating and measurement.Tree-Ring Bulletin & Tree-Ring Research4369–78. 10.1016/j.dib.2018.08.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HolmesR. L. (1994). Dendrochronology Program Library User’s Manual. Laboratory of Tree-Ring Research.Arizona City, AZ: University of Arizona.

  • Google Scholar

• 33

  HortonK. W.BedellG. H. D. (1960). White and Red Pine. Ecology, Silviculture and Management. Can. Dept. North. Affairs and Nat. Resour., For. Br., Bull.Ottawa, ON: Canadian Forest Service Publications, 124.

  • Google Scholar

• 34

  IPCC (2014). “Climate change 2014: synthesis report,” in Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edsPachauriR. K.MeyerL. A. (Geneva: IPCC), 151.

  • Google Scholar

• 35

  JohnsonJ. E. (1995). “The Lake States Region,” in Regional Silviculture of the United States, 3rd Edn, ed.BarrettJ. W. (New York, NY: John Wiley & Sons, Inc), 81–127.

  • Google Scholar

• 36

  KingJ. S.GiardinaC. P.PregitzerK. S.FriendA. L. (2007). Biomass partitioning in red pine (Pinus resinosa) along a chronosequence in the Upper Peninsula of Michigan.Canadian Journal of Forest Research.3793–102. 10.1139/x06-217

  • CrossRef
  • Google Scholar

• 37

  KottekM.GrieserJ.BeckC.RudolfB.RubelF. (2006). World Map of the Köppen-Geiger climate classification updated.Meteorologische Zeitschrift.15259–263. 10.1127/0941-2948/2006/0130

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  LambertM. C.UngC. H.RaulierF. (2005). Canadian national tree aboveground biomass equations.Canadian Journal of Forest Research351996–2018. 10.1139/x05-112

  • CrossRef
  • Google Scholar

• 39

  LamlonS. H.SavidgeR. A. (2005). A reassessment of carbon content in wood: variation within and between 41 North American species.Biomass and Bioenergy25381–388. 10.1016/s0961-9534(03)00033-3

  • CrossRef
  • Google Scholar

• 40

  LarssonL. A. (2018). Cybis CooRecorder v9.3.1 and Cybis CDendro v9.3.1.Saltsjobaden: Cybis Elektronik and Data AB.

  • Google Scholar

• 41

  LathamP.TappeinerJ. (2002). Response of old-growth conifers to reduction in stand density in western Oregon forests.Tree Physiology22137–146. 10.1093/treephys/22.2-3.137

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  LindseyR. (2018). Climate change: Atmospheric carbon dioxide. National Oceanic and Atmospheric Administration (NOAA).Washington, DC: NOAA.

  • Google Scholar

• 43

  MagruderM.ChhinS.PalikB.BradfordJ. B. (2013). Thinning increases climatic resilience of red pine.Canadian Journal of Forest Research43878–889. 10.1139/cjfr-2013-0088

  • CrossRef
  • Google Scholar

• 44

  MaguireD. A.MainwaringD. B.HalpernC. B. (2006). Stand dynamics after variable-retention harvesting in mature Douglas-fir forests of Western North America.Allg. Forst-u. J.-Ztg.177120–131.

  • Google Scholar

• 45

  Martin-BenitoD.PedersonN.FérrizM.Gea-IzquierdoG. (2021). Old forests and old carbon: A case study on the stand dynamics and longevity of aboveground carbon.Science of The Total Environment765142737. 10.1016/j.scitotenv.2020.142737

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  McDonaldR. I.UrbanD. L. (2004). Forest edges and tree growth rates in the North Carolina Piedmont.Ecology852258–2266. 10.1890/03-0313

  • CrossRef
  • Google Scholar

• 47

  NowackiG. J.AbramsM. D. (1997). Radial-growth averaging criteria for reconstructing disturbance histories from presettlement-origin oaks.Ecological Monographs67225–249. 10.1890/0012-9615(1997)067[0225:rgacfr]2.0.co;2

  • CrossRef
  • Google Scholar

• 48

  NuneryJ. S.KeetonW. S. (2010). Forest carbon storage in the northeastern United States: Net effects of harvesting frequency, post-harvest retention, and wood products.Forest Ecology and Management2591363–1375. 10.1016/j.foreco.2009.12.029

  • CrossRef
  • Google Scholar

• 49

  Ontario Ministry of Natural Resources and Forestry (OMNRF) (2015). Forest Management Guide to Silviculture in the Great Lakes-St. Lawrence and Boreal Forests of Ontario.Peterborough, ON: OMNRF.

  • Google Scholar

• 50

  PalikB. J.D’AmatoA. W. (2019). Variable retention harvesting in Great Lakes mixed-pine forests: emulating a natural model in managed ecosystems.Ecological Processes81–15.

  • Google Scholar

• 51

  PalikB. J.MitchellR. J.HiersK. J. (2002). Modelling silviculture after natural disturbance to sustain biodiversity in the longleaf pine (Pinus palustris) ecosystem: balancing complexity and implementation.Forest Ecology and Management155347–356. 10.1016/s0378-1127(01)00571-0

  • CrossRef
  • Google Scholar

• 52

  PalikB. J.MontgomeryR. A.ReichP. B.BoydenS. B. (2014). Biomass growth response to spatial pattern of variable-retention harvesting in a northern Minnesota pine ecosystem.Ecological Applications242078–2088. 10.1890/13-1173.1

  • CrossRef
  • Google Scholar

• 53

  PanY.BirdseyR. A.FangJ.HoughtonR.KauppiP. E.KurzW. A.et al (2011). A Large and Persistent Carbon Sink in the World’s Forests.Science333988–993. 10.1126/science.1201609

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  ParkerW. C.ElliottK. A.DeyD. C.BoysenE. (2008). Restoring southern Ontario forests by managing succession in conifer plantations.The Forestry Chronicle8483–94. 10.5558/tfc84083-1

  • CrossRef
  • Google Scholar

• 55

  PeichlM.ArainM. A.BrodeurJ. J. (2010). Age effects on carbon fluxes in temperate pine forests.Agricultural and Forest Meteorology1501090–1101. 10.1016/j.agrformet.2010.04.008

  • CrossRef
  • Google Scholar

• 56

  PetersonD. L.VoseJ. M.Patel-WeynandT (Eds) (2014). Climate Change and United States Forests.New York, NY: Springer, 261.

  • Google Scholar

• 57

  Pompa-GarcíaM.Venegas-GonzáLezA.JúniorA. A.Sigala-RodríguezJ. A. (2018). Dendroecological approach to assessing carbon accumulation dynamics in two Pinus species from northern Mexico.Tree-ring research74196–209. 10.3959/1536-1098-74.2.196

  • CrossRef
  • Google Scholar

• 58

  PowersM. D.PregitzerK. S.PalikB. J.WebsterC. R. (2010). Wood δ13C, δ18O and radial growth responses of residual red pine to variable retention harvesting.Tree Physiology30326–334. 10.1093/treephys/tpp119

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  PowersM. D.WebsterC. R.PregitzerK. S.PalikB. J. (2009). Spatial dynamics of radial growth and growth efficiency in residual Pinus resinosa following aggregated retention harvesting.Canadian journal of forest research39109–117. 10.1139/X08-162

  • CrossRef
  • Google Scholar

• 60

  PuettmannK. J. (2011). Silvicultural challenges and options in the context of global change: “Simple” fixes and opportunities for new management approaches.J. Forestry109321–331.

  • Google Scholar

• 61

  PuettmannK. J.AresA.BurtonJ. I.DodsonE. K. (2016). Forest restoration using variable density thinning: lessons from Douglas-fir stands in western Oregon.Forests7:310. 10.3390/f7120310

  • CrossRef
  • Google Scholar

• 62

  RobertsS. D.HarringtonC. A. (2008). Individual tree growth response to variable-density thinning in coastal Pacific Northwest forests.Forest Ecology and Management2552771–2781. 10.1016/j.foreco.2008.01.043

  • CrossRef
  • Google Scholar

• 63

  RudolfP. O. (1990). Red Pine. In, Burns, Russell M., and Barbara H. Honkala, tech. coords. 1990. Silvics of North America: 1. Conifers; 2. Hardwoods. Agriculture Handbook 654, Vol. 2. Washington, DC: U.S. Department of Agriculture, Forest Service, 877.

  • Google Scholar

• 64

  SherichK.PocewiczA.MorganP. (2007). Canopy characteristics and growth rates of ponderosa pine and Douglas-fir at long-established forest edges.Canadian Journal of Forest Research372096–2105. 10.1139/X07-105

  • CrossRef
  • Google Scholar

• 65

  SohnJ. A.HartigF.KohlerM.HussJ.BauhusJ. (2016). Potential of forest thinning to mitigate drought stress: a meta-analysis.Forest Ecology and Management380261–273. 10.1016/j.foreco.2016.07.046

  • CrossRef
  • Google Scholar

• 66

  SpeerJ. H. (2010). Fundamentals of Tree-Ring Research.Arizona City: AZ: University of Arizona Press.

  • Google Scholar

• 67

  StiellW. M. (1959). Seeding and Planting red and White pine. Dept. Northern Affairs and National Resources, Forestry Branch, Forest Research Division, Tech. Note No. 80.Ottawa, ON: Canadian Forest Service Publications, 18.

  • Google Scholar

• 68

  StokesM. A.SmileyT. L. (1968). An Introduction to Tree-Ring Dating.Chicago, IL: University of Chicago Press.

  • Google Scholar

• 69

  SWCRMP (St. Williams Conservation Reserve Management Plan) (2007). Executive Summary, Ontario Ministry of Natural Resources.Ontario, OT: Queen’s Printer.

  • Google Scholar

• 70

  UnwinG. L.KriedemannP. E. (2000). Principles and Processes of Carbon Sequestration by Trees (State Forests of New South Wales Technical Paper No. 64).West Pennant Hills, NSW: State Forests of New South Wales.

  • Google Scholar

• 71

  WilliamsonT. B.JohnstonM. H.NelsonH. W.EdwardsJ. E. (2019). Adapting to climate change in Canadian forest management: Past, present and future.The Forestry Chronicle9576–89. 10.5558/tfc2019-015

  • CrossRef
  • Google Scholar

• 72

  XingD.NielsenS. E.MacdonaldS. E.SpenceJ. R.HeF. (2018). Survival and growth of residual trees in a variable retention harvest experiment in a boreal mixedwood forest.Forest Ecology and Management411187–194. 10.1016/j.foreco.2018.01.026

  • CrossRef
  • Google Scholar

• 73

  YamaguchiD. K. (1991). A simple method for cross-dating increment cores from living trees.Canadian Journal of Forest Research21414–416. 10.1139/x91-053

  • CrossRef
  • Google Scholar

• 74

  YoungbloodA. (1991). Radial growth after a shelterwood seed cut in a mature stand of white spruce in interior Alaska.Canadian Journal of Forest Research21410–413. 10.1139/x91-052

  • CrossRef
  • Google Scholar