Regeneration patterns of key pine species in a mixed-pine forest indicate a positive effect of variable retention harvesting and an increase in recruitment with time

Many fire-dependent forests have experienced significant declines in species, structural, and functional diversity. These changes are attributed in part to traditional management approaches that were dominated by even-aged regeneration methods such as clearcutting. Variable retention harvesting (VRH) is an ecologically based forestry practice that involves retention of some mature overstory trees and other biological structures in the postharvest stand to emulate the effects of natural disturbance events. In this study, we examined the effect of a VRH treatment on recruitment of historically dominant pine species and understory vegetation two- and six years after its implementation in a naturally regenerated mixed pine forest in the Upper Peninsula of Michigan, USA. We found greater regeneration of red pine (Pinus resinosa L.) and eastern white pine (Pinus strobus L.) in VRH stands compared to the unharvested controls. Although red pine recruitment was very low during the first 2 years, both species increased significantly from year two to year six. Recruitment responses for the two species did not differ significantly between the aggregated and dispersed retention treatments. We also found high recruitment of deciduous tree species, consisting primarily of sprouts, and red pine exhibited a negative correlation with these sprouts. Treatment, time, and the interaction of treatment and time were significant factors in red pine recruitment. Time had a strong effect on the understory vegetation, with significant decreases in shrubs and ferns over time. The results suggest that VRH enhanced recruitment of red pine and eastern white pine and that this recruitment increased with time. We identify competition from sprouts and unfavorable seedbed conditions as the main factors that limited recruitment of red pine during the first 2 years. To reduce competing vegetation and improve seedbed conditions, we suggest following the VRH with a prescribed fire, mechanically removing sprouts, and implementing herbicide treatments. At a time when changes in global climate exacerbates the effects of traditional drivers of forest degradation, and with the need to maintain biodiversity, new ecologically based forest practices such VRH have the potential to facilitate regeneration of native trees and enhance the resilience of many fire-dependent forests.


Background
Changes in fire regimes, land use practices that replace forested areas with human structures, and forest management techniques that simplify forest ecosystems have adversely affected many pine-dominated forests in the Great Lakes region of the United States (Whitney 1987;Cleland et al. 2004;Drobyshev et al. 2008a). The adverse effects include shifts in plant species composition, inadequate regeneration of historically dominant pine species, structurally simplified stands, declines in wildlife habitat, and accumulation of fuels (Frelich 1995;Drobyshev et al. 2008b;Corace III et al. 2014;Nyamai et al. 2014). Traditional management practices in these forests have been based on models whose effects are not analogous to those of the natural disturbance processes that shaped their development (Seymour et al. 2002;. Clearcutting as one of these practices has been associated with impacts that simplify the structural and compositional complexity of a stand, resulting in reduced biodiversity, habitat functions, and ecological resilience of these forests among other things (Lindenmayer and Franklin 2002;Puettmann et al. 2009;Lindenmayer et al. 2012). As a result, new forest management practices that impact stands in ways that more closely emulate the effects of natural disturbance events have been adopted in many places (Franklin et al. 2007;Gustafsson et al. 2012;Lindenmayer et al. 2012;Fedrowitz et al. 2014;Puettmann et al. 2015;. Among these practices is "Retention forestry", a management approach that emphasizes creation of multiaged stands and retention of mature overstory trees in varying structural patterns to achieve various objectives (Franklin et al. 2007;Gustafsson et al. 2012;Lindenmayer et al. 2012;Fedrowitz et al. 2014;Mori and Kitagawa 2014).
The more widely adopted retention forestry practice is variable retention harvesting (VRH). VRH is defined in Helms (1998) as "an approach to harvesting based on the retention of structural elements or biological legacies from the harvested stand for integration into the new stand to achieve various ecological objectives". Implementation of VRH involves retention of some mature overstory trees and other biological structures in various structural patterns during harvesting operations to achieve specific management objectives (Franklin et al. 2007;. The three key characteristics of VRH are 1) what is retained, 2) how much is retained, and 3) the spatial pattern of retention. Pattern of retention refers to the spatial distribution of retained structures, which can be either aggregated or dispersed. Aggregated retention comprises a mosaic of open patches where nearly all trees are removed, and patches where trees are retained in clusters. In dispersed retention, the retained trees are uniformly distributed across the stand. This variation in pattern of retained structures is thought to provide opportunities to emulate various structural features created in stands by different types and intensities of natural disturbance events (Franklin et al. 2007).
VRH has been applied in and effects on various forest components studied across the Great Lakes region (Bebber et al. 2005;Peck et al. 2012;Roberts et al. 2016Roberts et al. , 2017, in other parts of the United States (Battaglia et al. 2002;Aubrey et al. 2009;Franklin and Donato 2020), and across the globe (Mori and Kitagawa 2014;Beese et al. 2019;Scott et al. 2019). Within the northern Lake States region, majority of these studies have been conducted in forests that in the past were managed using even-aged regeneration methods such as clearcutting, followed by planting pine seedlings . To expand our knowledge and inform management, there is need to examine the effectiveness of VRH in restoring key species in naturally regenerated forests that have not been clearcut or planted since historical logging by EuroAmerican settlers. Additionally, given its documented benefits Fedrowitz et al. 2014), including potential to enhance the resilience of forest ecosystems to the impacts of climate change (Nagel et al. 2017), information on the viability of VRH in diverse ecosystem types would be beneficial to its expansion to other geographical and management contexts. To address these needs, we conducted an experimental study that examined the effect of VRH on recruitment of key species in a naturally regenerated mixed pine forest in the eastern Upper Peninsula of Michigan, USA.
The VRH in our study area was implemented as a restoration application to promote recruitment of red pine (Pinus resinosa L.) and eastern white pine (Pinus strobus L.). The two pine species were targeted because they are the historically dominant species in these mixed-pine forests under unaltered conditions. A secondary objective was to reduce jack pine (Pinus banksiana Lamb.) and broadleaved tree species that had increased in abundance in these stands, contributing flammable material and creating hazardous fire conditions (Drobyshev et al. 2008b). The VRH created both aggregated and dispersed retention patterns in the harvested stands. We sought to answer the following questions: 1) what is the effect of VRH on recruitment of red pine, eastern white pine, and other vegetation two-and six years following application, and does this effect differ by pattern of retention? 2) how does recruitment change between year two and year six following the harvest? We examined responses using the two time periods because they represent a gradation of site conditions and competitive interactions following disturbance, which may impact tree recruitment differently. Examining recruitment over time also allows us to examine potential trends in recruitment of the two target species, particularly given that red pine is often difficult to regenerate in many areas (Ahlgren 1976;Godman and Mattson 1976). We hypothesized that there would be greater recruitment of red pine and eastern white pine in the VRH stands compared to the unharvested control stands due to increased availability of resources including greater light penetration into the understory and decreased competition since the stand density would be reduced. We also hypothesized that recruitment would be greater in the initial years (e.g., during the first year or two) because the newly altered overstory and forest floor conditions are likely to be more favorable for recruitment of the two pine species during this initial phase of recolonization. Finally, we hypothesized that the effects of VRH on recruitment would differ between the two retention patterns because of the variation in stand structure created by the two retention treatments.

Study area
The study was conducted at the Seney National Wildlife Refuge (SNWR), located in the Upper Peninsula of Michigan, USA (46.271594°N, 86.057078°W). SNWR lies within the Seney Lake Plain ecoregion and is characterized by soils and physiographic features shaped by postglacial erosion and soil formation processes (Albert 1995). All the stands included in this study were on soils ranging from well to excessively drained coarse sands. Temperatures vary between an average minimum of − 13°C in the coldest month and 27°C in the warmest month, average annual precipitation is 78 cm, and growing season averages about 119 days (MRCC 2020). The stands were also within the eastern white pine/blueberry (Vaccinium)/trailing arbutus (Epigaea repens) habitat type as described by Burger and Kotar (2003). Presettlement mixed-pine forest ecosystems were dominated by red pine and eastern white pine, with minor components of jack pine and a variety of deciduous species (Whitney 1987). Current stands have been significantly altered, having a large component of jack pine and deciduous species in many areas. These altered naturally regenerated stands, coupled with available information on stand structure, composition and disturbance history from prior studies within SNWR (Drobyshev et al. 2008a, b;Corace III et al. 2013;Nyamai et al. 2014) makes the study area a good model system to test emerging restoration and management techniques such as VRH.

Harvest treatments and field methods
The study was conducted within nine stands ranging in size from 3.2 to 6.4 ha. Selected stands were randomly assigned to treatments, with three receiving the aggregated retention, three receiving the dispersed retention, and three unharvested controls. For the aggregated retention, overstory basal area was reduced by creating within the stand a mosaic of completely open patches (~0.3 ha) and patches with trees retained in clusters ( Fig. 1). For the dispersed retention, overstory basal area was reduced by leaving retained trees uniformly dispersed across the stand (Fig. 1). Both treatments targeted 30% retention of the initial overstory basal area in each stand (which is~22 m 2 per ha). We used unharvested stands for controls instead of clearcuts because the focus of this study was to understand what happens when VRH is implemented in current degraded stands compared to when nothing is done and stands left as they are. Additionally, clearcutting is not a practice conducted in our study system since it is a national wildlife refuge where more ecologically based management approaches are required. Other studies have also used unharvested stands as controls (e.g., Halpern et al. 2005;Roberts et al. 2016Roberts et al. , 2017Lindenmayer et al. 2019).
Deciduous trees and jack pine were primarily targeted for removal, although in very few cases some red pine trees were removed to achieve the retention target. Additionally, all subcanopy stems including saplings of species targeted for removal were taken out. As a result, retained overstory trees in the VRH stands primarily consisted of red pine, eastern white, and any snags of all tree species. Logging was conducted during the growing season and a "timber processor" with a mechanized head on an extendable boom was used. Trees were delimbed to approximately 2.6 m and removed from the stands with a forwarder. All equipment used had tires and not tracks. The cut trees were de-limbed on site, merchantable boles taken out, and the smaller-size woody debris and slash were left on site and spread across the ground surface.
Within each stand, we established a network of six 25m 2 (5 m by 5 m) plots that we used to collect data on tree seedlings and cover of herbs, shrubs, and ferns. Since the aggregated retention consists of both open patches and patches with tree clusters, plots in these stands were established by assigning three plots in the open patches and three within the tree clusters. Although edge effects can be an important factor in the aggregated retention (Powers et al. 2008;Franklin et al. 2018), our objective was to capture responses that represent the average conditions of both the open and retained tree patches represented in the aggregated treatment. To establish plots, we established a grid system over each stand in ArcGIS 10.1 (2012) and generated several grid intersections. We then randomly selected six grid intersections that served as centers of six replicate study plots in each stand. The result was a total of 54 plots in this study (n = 18 in aggregated, n = 18 in dispersed, and n = 18 in control). Within each plot, we recorded count of all tree seedlings (stems < 2.5 cm dbh). For deciduous species that exhibited sprouting, each individual stem was counted separately. To discriminate between seedlings recruited pre-and posttreatment, we estimated the age of tree seedlings by counting number of whorls for pine species, and number of terminal bud scars for deciduous species. The stands did not appear to have much advance regeneration based on visual assessments, and nearly all the seedlings recorded were younger than 2 years at the time of initial sampling, based on our age estimation. We also estimated the percent cover of herbs, shrubs, and ferns. Taxonomic authority used in plant identification followed the PLANTS database (USDA NRCS 2014).

Statistical analyses
To prepare data for analysis, we summarized seedling and cover data by plot. We then examined differences between treatments using Analysis of Variance (ANOVA) in SAS (SAS v. 9.4 SAS Institute, Cary, NC) and conducted pairwise comparisons of these treatments using Tukey's HSD post-hoc test. ANOVA was conducted on the data for each time periods individually, then on the combined data of both time periods to examine treatment and time interactions. Understory vegetation particularly shrubs in the Ericaceae family have been found to influence recruitment of pine trees (Roberts et al. 2016(Roberts et al. , 2017. We were thus interested in how these understory species might influence red pine and eastern white pine recruitment in the VRH stands. To examine that, we conducted a Canonical Correspondence Analysis (CCA) using CANOCO v. 5.12 (Braak and Smilauer, Ithaca, NY, USA). CCA is a direct gradient ordination analysis used to examine how measured predictor variables drive patterns of species occurrence and distribution in space. We used the data for herbs, shrubs, ferns, red maple (Acer rubrum L.), northern red oak (Quercus rubra L.) and aspen (Populus spp.) as predictor variables for the occurrence of red pine and eastern white pine. Red maple, northern red oak, and aspen were included as predictors because of their abundance based on the seedling density data. The effect of predictor variables on the response variables was tested using an unrestricted random Monte Carlo permutation test with 999 permutations and the default linear combination (LC) site scores.

Tree seedling densities
Two years after treatment, recruitment of red pine was low overall, with seedlings observed only in the dispersed retention (104 seedlings·ha − 1 ) and no treatment effect (F 2,51 = 1.89, P = 0.162; Fig. 2, Table 1). Compared to red pine, we found higher density of eastern white pine seedlings and significant treatment effects (F 2,51 = 3.51, P = 0.037) in year two. There were 659 seedlings·ha − 1 in the aggregated retention, 763 seedlings·ha − 1 in the dispersed retention, and 138 seedlings·ha − 1 in the control (Table 1). These densities only differed significantly (P = 0.043) between the dispersed retention and the control in year two, with no significant difference between the two retention treatments (P = 0.911; Fig. 2). After 6 years, both red pine and eastern white pine had increased significantly in the VRH stands ( Fig. 2, Table 1). Unlike year two, there was significant treatment effects for red pine in year six (F 2,51 = 4.76, P = 0.013), the highest seedling density was found in the aggregated retention (312 seedlings·ha − 1 ), the density in the dispersed retention doubled (208 seedlings·ha − 1 ), and there was still no recruitment in the control stands (Fig. 2, Table 1). As was in year two, recruitment did not differ significantly between the two retention treatments in year six (P = 0.57). We also found significant treatment effects for eastern white pine in year six (F 2,51 = 5.20, P = 0.008), with recruitment only differing significantly between the aggregated retention and the control (P = 0.007) but not between the two retention treatments (P = 0.636; Fig. 2). Similar to red pine in year six, the highest seedling density for eastern white pine was found in the aggregated retention (1180 seedlings·ha − 1 ) with the least recruitment in control stands (104 seedlings·ha − 1 ; recruitment actually declined in year six in these stands (Table 1). Among the other tree species, red maple exhibited the most recruitment in VRH stands in both time periods, with highest density in the aggregated retention (Table 1) in both years, and significant treatment effects in year two (F 2,51 = 5.39, P = 0.007) and year six (F 2,51 = 4.55, P = 0.015). Recruitment for red maple differed significantly between the two retention treatments in year two (P = 0.02) but not in year six (P = 0.086) and increased between the 2 years (Table  1). Northern red oak and aspen showed similar trends of high recruitment in VRH stands and increases in recruitment with time since harvest (Table 1). When the two time periods were considered together, both treatment and time were significant factors in red pine recruitment with time having a stronger effect, and a significant time by treatment interaction (Table 2). Conversely, treatment was the only significant factor in eastern white pine recruitment (Table 2). Red maple, northern red oak, and aspen were all strongly affected by both treatment and time, although the interaction between these two factors was not significant for these species (Table 2).

Cover of shrubs, herbs, and ferns
For the understory vegetation, lowbush blueberry (Vaccinium angustifolium Aiton) and velvet-leaf blueberry Fig. 2 Comparison of red pine and eastern white pine seedling density (seedlings·ha −1 ) among treatments 2 and 6 years since VRH treatment in a mixed-pine forest ecosystem in the eastern Upper Peninsula of Michigan, USA. Letters a and b are used to indicate significant differences among treatments for the 2 year data, while letters x and y are used to indicate significant differences among treatments for the 6 year data (P < 0.05) (Vaccinium myrtilloides Michx.) dominated the shrub layer, while bracken fern (Pteridium aquilinum (L.) Kuhn) was most common in the herbaceous layer (Table 3). VRH did not appear to have a significant effect on shrub recruitment in the initial periods following harvest (F 2,51 = 0.37, P = 0.69), but there was a treatment effect after 6 years (F 2,51 = 7.67, P = 0.001; Fig. 3). Shrub cover decreased an average of 78% in VRH stands between the two time periods. Herbaceous vegetation response was nearly opposite that of shrubs. There were significant treatment effects on herbaceous vegetation after 2 years (F 2,51 = 3.76, P = 0.03), but the difference was only between the aggregated retention and the control (P = 0.04; Fig. 3).
Although there was an average of 65% increase in herbaceous vegetation from year two to year six in VRH stands, there was no treatment effect for year six (F 2,51 = 1.02, P = 0.36; Fig. 3). Bracken fern exhibited similar patterns to those of the shrubs; there were no treatment effects after 2 years (F 2,51 = 2.73, P = 0.07), but there were treatment effects in year six (F 2,51 = 5.42, P = 0.007), with the cover differing only between the dispersed retention and the control (P = 0.006; Fig. 3). There was an average of 86% decrease in the cover of ferns in VRH stands between the two periods. When the two time periods were considered together, treatment had a significant influence on ferns only, time had a very strong influence on all three vegetation groups, and the interaction between treatment and time was a significant factor only for herbaceous vegetation ( Table 2).
Influence of other vegetation on red pine and eastern white pine recruitment CCA analyses revealed significant patterns of species occurrences along gradients associated with the predictor variables for both the two-year (F = 5.6, P = 0.001) and the sixyear (F = 3.8, P = 0.001) time periods. The four canonical axes for the two-year data accounted for 42% of the variation in the species composition data, while that of the six-year data accounted for 33%. Two years after VRH, red pine seedlings exhibited a strong positive correlation with the first gradient (score = 0.42; Fig. 4). Along this gradient, red pine seedlings were high in areas with high shrub cover but low cover of herbs and ferns, as well as low densities of maple, oak, and aspen seedlings (Fig. 4, Table 4). Eastern white pine on the other hand had a strong positive

Bracken fern
Pteridium aquilinum (L.) Kuhn association with the second gradient (score = 0.45), with high seedling density in areas with aspen but with low density and cover of all the other species (Fig. 4, Table 4). During the first 2 years, shrubs, ferns, oak, and maple exerted significant influence on recruitment of red pine and eastern white pine, while aspen and herbs did not have significant effects (Table 4). In contrast to responses in the first 2 years, red pine had a strong negative correlation with the first gradient after 6 years (score = − 0.56; Fig. 4), with high density in areas of high herbaceous cover but low densities of maple, oak, and aspen seedlings (Fig. 4, Table  4). Eastern white pine exhibited a marginal positive association (score = 0.12) with the first gradient and a moderate positive association (score = 0.27) with the second gradient ( Fig. 4) in this time period. Most of the eastern white pine seedlings were positively associated with the other cooccurring tree species including maple, oak and aspen (Fig. 4, Table 4). Only maple, oak, and ferns exerted  (Table 4).

VRH effects on recruitment
When used for regeneration purposes, the application of VRH is typically based on the premise that by altering the structure of the overstory, the understory resource environment is likely to be affected in ways that favor recruitment and establishment of desired tree species. Our findings suggest that VRH can promote recruitment of restoration target species, in our case red pine and eastern white pine. Further, our findings indicate a significant time and interactive treatment by time effects for red pine, with greater recruitment over time seemingly driven by opportunistic resource availability. Given the silvical characteristics of red pine and eastern white pine, we attribute the responses observed in this study to several factors as discussed below.

Dynamics of resource competition
Following VRH operations, greater survival and growth of planted pine seedlings have been attributed to increases in resource availability (Boyden et al. 2012;Peck et al. 2012;Montgomery et al. 2013). We, therefore, anticipated that the structural changes brought about by the VRH at our study site would create conditions that would promote recruitment of both red pine and eastern white pine (Ahlgren 1976;Carey 1993). Our findings, however, show that red pine exhibited very low regeneration in VRH stands during the first two years but increased significantly by the sixth year, contrary to our hypothesis. In contrast, eastern white pine showed better recruitment in the VRH stands in both years, also increasing substantially from year two to year six. We suggest that the recruitment patterns observed for the two species was primarily driven by competition and believe that the contrasts in responses observed between red pine and eastern white pine are likely consistent with their differences in sensitivity to competition and substrate requirements for germination (Ahlgren 1976;Carey 1993). The primary source of this competition for resources is likely recruited seedlings of deciduous species such as red maple, northern red oak, and aspen (Abrams 1998;Fei and Steiner 2009). We attribute the massive recruitment of these deciduous species to resprouting, and believe that it was disturbance mediated, likely by the VRH. It is likely that the sprouts exerted significant pressure on resources during the initial years, placing red pine at a greater competitive disadvantage (Dietrich et al. 2017). Ordination analysis also revealed significant negative influence of oak and maple on red pine, which supports the suggested impact of these species on red pine recruitment. Although eastern white pine increased significantly in VRH stands between year two and year six, its density was higher in VRH stands in both time periods relative to those of red pine, suggesting that competition was not as significant a factor in its recruitment as it was for red pine (Ahlgren 1976). The poor recruitment of the two pine species in the control stands coupled with increases in recruitment of competitor deciduous species over time may also be informative from a management perspective. They suggest the potential for significant and possibly irreversible alterations in species composition and stand structure in these forests if no restoration interventions are implemented. We found a strong effect of time and a significant interaction of treatment and time on red pine recruitment, while treatment was the only significant factor in eastern white pine recruitment. Both pine species also increased substantially between year two and year six. These findings suggest that although VRH enhanced recruitment of both pine species, these effects varied over the six-year duration, particularly for red pine. The new cohort of competing species discussed not previously as it is now may have exerted such intense pressure for resources during the initial years, that only a few individuals of red pine were able to survive in the competitive microsite conditions. Although this resource pressure appears to have persisted given the concurrent increases in recruitment of competitor species, we think that its effects were moderated over time by availability of recruitment opportunities. These opportunities may have been created for example by decreases in recruitment of understory vegetation, and/or as a result of decay of slash and other debris on the forest floor. The newly available space and nutrients likely supported additional red pine germination and establishment after the initial years. The potential for red pine to be more opportunistic for available resources long after the immediate period following disturbance as our findings suggests is an interesting finding given its restrictive site requirements (Hauser 2008). Besides competition by other tree species, it could be argued that advance regeneration in these stands can be an additional factor that could have had some effect on regeneration of the two pine species. Advance regeneration of deciduous species has been shown to respond to canopy gaps, making them more competitive for resources during the initial years (Reuling et al. 2019). We do not believe, however, that this was a significant factor in this study given the targeted reduction during harvest of unwanted advance regeneration.
We suggest some management recommendations for this, and similar mixed-pine ecosystems based on our discussion above regarding the dynamics of space and resource use between the pine and the deciduous species. When increasing natural regeneration of red pine is an objective, we recommend taking steps during or after the logging operation to limit resprouting of deciduous tree species. Prescribed fire that targets fire-intolerant competitors such as maple can be implemented following VRH treatments. Managers can also explore the possibility of mechanical removal of sprouts and/or herbicide treatment among other options. Since eastern white pine appears to do just fine with the VRH treatment based on our findings, the additional interventions we suggest here may not be as critical for eastern white pine as they are for red pine in cases where the primary objective is to regenerate eastern white pine.

Effect of understory vegetation
Understory vegetation in the herb and shrub layers can also be an important source of competition for naturally regenerating pine species. Retention harvests have been shown to influence recruitment and growth of understory vegetation (Halpern et al. 2012). By preempting resources, the understory vegetation can then place seedlings of woody species at a competitive disadvantage and limit treatment effects (Weyenberg et al. 2004;Montgomery et al. 2013;Roberts et al. 2017). We found no significant VRH effects on shrubs and ferns in year two, but the effects were significant in year six. Both shrubs and ferns decreased significantly in VRH stands between the two time periods, and only the herbaceous component increased. Ordination analysis revealed a strong positive association between red pine and shrubs in the initial years, but this spatial relationship changed by year six where red pine begins to associate more with the herbaceous vegetation. Eastern white pine on the other hand did well in areas with less shrubs in the initial years but the association became more positive by year six. The positive association between red pine and shrubs was contrary to our expectation that recruitment of both pine species would be limited in areas with high shrub cover (which was mostly blueberry). We had thought that the shrub component would be a source of resource competition, in addition to other potential inhibition mechanisms (Mallik 2003). We are not certain about the exact mechanism behind the positive relationship in the initial years for red pine. It could be that the shrubs provided some sort of refuge to newly germinated seedlings in the initial years, possibly just the right amount of shade from direct sunlight. Such a relationship would not be common or typical for an intolerant species like red pine though.
The effect of time was stronger for all three understory vegetation groups relative to the effect of treatment. Shrubs and ferns decreased by 78% and 86% respectively, while herbaceous vegetation increased by 65%. We argue that decreases in shrubs and ferns likely contributed to the recruitment dynamics observed for red pine specifically. As described earlier, we think that these decreases potentially freed some resources that red pine then took advantage of. The strong effect of time also suggests that recruitment of shrubs and ferns may increase in the initial years after VRH but decreases over time. While we highlight this temporal aspect of VRH effects, we caution that any facilitative effects may not translate to significant recruitment impacts for the pine species if the sprouting of competing tree species is not addressed. We also note that there were no shrub reduction treatments before the VRH application in our study, as has been done in other studies within the region (Boyden et al. 2012;Peck et al. 2012;Montgomery et al. 2013;Roberts et al. 2017).

Seedbed requirements
Compared to eastern white pine, red pine requires specific seedbed conditions, with more successful natural germination on exposed mineral soil, low duff and litter cover on the forest floor, and little competition from other vegetation (Ahlgren 1976;Ahlgren and Ahlgren 1981;Mallik and Roberts 1994). Low intensity surface fires typically create these conditions and facilitate natural recruitment of both red pine and eastern white pine in our study area (Drobyshev et al. 2008a;Nyamai et al. 2014). Mechanical scarification of the soil surface has also been used in other pine systems to promote successful seedling emergence and establishment (Landhäusser 2009). In this study, nearly all small-sized woody debris were left on site following the VRH. The reasoning at the time was that these materials would serve to provide habitat for other organisms and would break down over time. Although we did not examine the influence of litter on the forest floor, it is reasonable to think that this debris covering the ground surface likely contributed to the unfavorable seedbed conditions, possibly by inhibiting seed germination (Halpern and McKenzie 2001).
It has also been suggested that harvesting operations that do not significantly impact the forest floor to expose a mineral soil seedbed can provide a competitive advantage to other species that regenerate from root suckers and stump sprouts (Wiensczyk et al. 2011;D'Amato et al. 2012). We suspect that the VRH implemented in this case did not impact the seedbed enough to create adequate conditions favorable for red pine germination, and instead may have provided an advantage to competitor tree species. The fact that no post-harvest scarification or other forest-floor treatments were conducted on the treated stands may also have contributed to the unfavorable seedbed conditions. Scarification of the forest floor after a logging operation can be critical in removing surface litter and exposing mineral soil, both of which can enhance successful germination of lightseeded species such as the pines we studied (Reuling et al. 2019). Techniques such as mastication that involve shredding, chopping, or chipping of the slash have also been used in conifer systems in other places (Reiner et al. 2009), and can be explored in our study area.

Seed source and overstory effects
Although resources at the microsite level may be strongly impacted by factors such as competition and seedbed characteristics as discussed above, availability of adequate seed sources is also important. Residual trees on the harvested stands are likely to be the primary source of seeds for new colonization in mixed-pine stands, although some seeds may be dispersed from other sources (Carey 1993;Hauser 2008). We believe that the residual red pine and eastern white pine trees in our VRH stands can provide an adequate seed source for new recruitment based on visual assessment and the number of retained mature trees (basal area of 22 m 2 per ha). Potentially supporting this assessment, a study in our region found that seed trees left after logging activities were associated with quicker regeneration and a more diverse growth structure in the post-disturbance regenerating population (Palik and Pregitzer 1994). Although seed supply would be a function of seed crop productivity (Ahlgren 1976;Hauser 2008), we do not think that lack of a seed source was a significant contributor to the low recruitment of red pine. Nevertheless, additional studies on seed crop dynamics of the two pines species in this study would likely shed more light on the quantity and viability of seeds produced by the residual trees.
The nature of the overstory composition and structure can also influence light, temperature, and moisture conditions in the understory (Oliver and Larson 1996). This VRH treatment targeted a 30% retention of the initial stand basal area in each of the two retention treatments. This retention level was selected as a middle ground between higher retention levels that may not adequately modify understory environmental conditions (Aubry et al. 2009), and lower retention levels that may create extreme changes, thereby giving competitor species an advantage (Peck et al. 2012;Reuling et al. 2019). Given an adequate seed source, coupled with its versatile characteristics (Abrams 2001), eastern white pine has been shown to have better recruitment even under a high overstory cover (Dovčiak et al. 2003). Comparatively, red pine is more shade intolerant with greater recruitment success in more open stands where there is greater light penetration reaching the forest floor (Hauser 2008). We cannot conclusively infer the impact of the 30% retention level on the observed recruitment patterns in this study as we did not have comparative tests of other retention levels. Given the silvical characteristics of the two pine species however, we think that this retention level provided adequate overstory conditions for their recruitment. We also believe that the previously described factors of competition and seedbed condition were greater drivers of recruitment in these stands compared to the impact of the level of retention.

Spatial pattern of retention
The two retention treatments (aggregated and dispersed) have been thought to have different effects on resource availability and plant competitive interactions (Boyden et al. 2012;Montgomery et al. 2013). By allowing for spatial variation in overstory structure and composition, the two retention patterns can provide contrasts in postharvest stand structure. This diversity of structures can provide opportunities for some ecological objectives to be met by dispersing retained structures and other objectives by aggregating them (Franklin et al. 2007;Halpern et al. 2012). Variations in spatial pattern of post-treatment structures can also provide heterogeneity in competitive environments across the stand, which can be targeted to favor recruitment of specific species (Boyden et al. 2012). Further, the retention patterns may provide opportunities to emulate patterns created on the landscape by different disturbance events or by varied intensities of such events (Franklin et al. 2007).
Literature suggests that red pine regeneration would be more successful in the aggregated retention while eastern white pine is likely to do well in both retention patterns . This difference has been attributed to the potential for greater changes in competitive environments in gaps created by the aggregated retention treatment, leading to favorable conditions for regeneration of pine species such as red pine that have more specific germination requirements (Battaglia et al. 2002;. While we largely expected that the differences described above would hold for our study, we did not find consistent differences in the effects of VRH between the two retention treatments. Recruitment of red pine and eastern white pine did not differ significantly between the two retention treatments during both time periods as we had hypothesized. These findings compare with those of other studies that have found no significant differences among retention patterns when survival, growth, and performance of planted red pine, eastern white pine, and jack pine seedlings were examined (Powers et al. 2011;Montgomery et al. 2013). Other studies have found differences among retention patterns with regards to growth of planted pine seedlings and resource availability Boyden et al. 2012;Peck et al. 2012). We think that the relative impact of the two retention treatments in our study may have been moderated by the strong effect of competitor vegetation and the unfavorable seedbed conditions. These factors may have had a strong confounding effect, potentially masking differences in treatment effects between the retention patterns. It is also possible that there are just no marked differences in the responses of these two pine species between the aggregated and dispersed retention patterns when VRH is implemented at the 30% retention level. If the latter is the case, and depending on the scale and management objectives, such similarities in species responses between retention patterns may provide managers with the flexibility to accommodate multiple management objectives (Franklin et al. 2018), especially where such objectives are complimentary.
We would like to conclude our discussion by noting that this study examined natural regeneration of red pine and eastern white pine in response to VRH, with no pre-harvest site preparation treatments or seeding of the desired species. We also examined responses only up to six years since harvest, which is a relatively short successional timeframe for long-lived species such as red pine and eastern white pine. A longer-term study, with additional components that examines for example changes in the down woody debris or seed production by the retained mature trees could provide additional insights into other drivers of recruitment in these stands. We also note that VRH as a forest management practice has certain limitations that are often a function of the condition of the forest, the size of the management unit, the landscape context, and management objectives, as highlighted by Palik et al. (2020). They advise that before embarking on a VRH treatment, an assessment needs to be made as to whether the range of structures created by the retention treatments represents characteristics typical of the ecosystem following a natural disturbance. Also, because you are creating a matrix of varied structures within a given forest area (i.e., some trees removed and others retained in various patterns), a large management unit is likely to be more supportive of the practical and operational aspects of VRH compared to a smallersized unit. Finally, if the forest is managed for commercial timber production, the potential economic cost and priority tradeoffs needed to diversify structural and compositional features in the stand should be considered carefully when adopting VRH as a management practice.

Conclusions
Our findings indicate greater recruitment of both red pine and eastern white pine in VRH stands and significant increases in recruitment with time. Although both species increased significantly from year two to year six, red pine exhibited very low recruitment during the first two years while eastern white pine regeneration was higher in both years. We attribute the low red pine recruitment in the initial years to competition from sprouts of deciduous tree species, and unfavorable seedbed conditions. Eastern white pine on the other hand likely did better due to its less strict germination and establishment requirements, including seedbed conditions, competition, and shade tolerance. The lack of recruitment of red pine, decline in recruitment of eastern white pine with time, and significant increases in recruitment of deciduous trees with time in the unharvested stands all suggest a potential for significant alterations in species composition and stand structure beyond thresholds of recovery, in the absence of restoration interventions.
We hope that our findings contribute to expanding the application of VRH in similar forest types in the region and beyond. We believe that adoption and expansion of ecologically based silvicultural approaches is critical as we seek to enhance resilience and adaptive capacity of forest ecosystems across the globe so they can better survive the impacts of climate change and minimize their vulnerability. Finally, while the general concepts and principles of VRH are typically similar and hence replicable in similar ecosystem types, treatments are likely to be most effective when tailored to suit the specific ecosystem type and management objectives.