Timing is Not Everything: Assessing the Efficacy of Pre- Versus Post-Harvest Herbicide Applications in Mitigating the Burgeoning Birch Phenomenon in Regenerating Hardwood Stands

Sweet birch (Betula lenta L.) is aggressively recruiting in temperate forest understories of the eastern United States and often dominates the post-disturbance seedling community, diminishing diversity and hindering sustainable silviculture. The type and timing of silvicultural actions affect birch recruitment via their effects on seedling recruitment, survival, and growth. Here, we examine birch regeneration under two contrasting treatment sequences: preversus post-shelterwood harvest herbicide application (H–S vs. S–H) in combination with white-tailed deer (Odocoileus virginianus Zimmerman) browsing (fenced vs. unfenced) at 22 sites in northwestern Pennsylvania, USA. Additionally, we examine how treatments interact with additional site factors, including potential propagule sources and site productivity (i.e., integrated moisture index). We found the S–H sequence initially reduced birch density by 71% relative to the H–S sequence; however, the magnitude of this reduction waned over five growing seasons. Furthermore, birch proliferated following the H–S sequence only where mature birch were present. Deer browsing reduced birch height by 29% relative to fenced areas protected from browsing; however, by the fifth growing season birch seedlings were over twice as tall as other hardwood species across all treatments. Finally, increasingly mesic sites enhanced birch height growth. In sum, although post-harvest herbicide (S–H) provides short-lived control over birch, land managers should also consider browse pressure, seed source, and site productivity, as these may enhance or diminish the efficacy of post-shelterwood herbicide sequence effects on birch.


Introduction
Within the last two decades, sweet birch (Betula lenta L.; hereafter, birch) has emerged among the few species that consistently and abundantly recruit into mesophytic and northern hardwood forests of the eastern United States. Forest Inventory and Analysis (FIA) data show, for example, sweet birch sapling (2.5-12.5 cm diameter at breast height (DBH)) density increases of 11% from 2004 to 2014 in Pennsylvania and 24% from 1993 to 2007 in New York [1,2]. This upsurge is remarkable given that, unlike other species increasing in forest understories (e.g., beech (Fagus grandifolia L.)), the relative abundance of birch in the canopy is disproportionately low compared to its abundance in forest understories [1]. The proliferation of birch is particularly pronounced following overstory disturbances (e.g., pests, wind, and forest harvests), leading to situations where birch overwhelmingly dominates the regeneration layer [3][4][5][6][7]. Thus, the burgeoning birch phenomenon poses a serious challenge to forest managers who strive to sustain forest diversity and values in the face of seedling demography trends that suggest future forests are becoming increasingly monodominant.
Birch increases in forest understories appear to be the 21st century analogue to the red maple (Acer rubrum L.) surge of the 20th century that has now stabilized or even declined in many regions [1,2,8]. Like red maple, birch may be considered a 'super-generalist' (sensu [9]) in that it possesses wide amplitude in key traits that confer recruitment, survival, and growth advantages under a variety of conditions, including recently disturbed stands. This versatility is, perhaps, most evident in its shade tolerance, which has been categorized as intolerant to intermediate [10,11], yet research demonstrates birch seedlings are able to survive and grow even under extremely low light levels (i.e., <5% of full sun; [12,13]). The species' versatility in ecological attributes is not limited to shade tolerance. For example, birch is not hypothesized to be strongly recruitment limited, as mature trees are prolific and frequent seed producers with broad seed dispersal capabilities [10,14]. Dispersed seeds germinate best on exposed mineral soil, which often increases following harvesting; however, they also germinate readily in undisturbed sites [15,16]. Further, although birch seedlings are browsed by white-tailed deer (Odocoileus virginianus Zimmerman), the species is relatively insensitive to browsing [17] and can dominate forest understories even under moderately high deer populations [18,19]. Finally, it is also possible that birch benefits from changes in soil nutrient availability, including the nitrogen pulse that exists post-overstory disturbance, with enhanced establishment, particularly when light is not limiting [20][21][22].
In actively managed forests, the intensity and timing of forest harvesting, along with the additional silvicultural practices that often accompany cutting, have the potential to either promote or mitigate birch dominance. For example, the degree to which harvesting operations disturb soil, enhance nitrogen mineralization, and increase light may regulate birch recruitment success [23]. The timing of key silvicultural actions may also greatly affect birch establishment. Here, we focus on the regeneration dynamics that occur following two contrasting treatment timing sequences: pre-versus post-shelterwood harvest herbicide applications. Within the Allegheny hardwoods and northern hardwood forest types, the herbicide-shelterwood sequence was the recommended approach during the late 20th century in an effort to chemically control pre-existing and highly recalcitrant herbaceous competing vegetation either before or at the time of the shelterwood harvest [24][25][26]. Over the past two decades, however, managers have shifted to a shelterwood-herbicide sequence in an effort to limit dominance by woody species that establish aggressively post-harvest [27]. The rationale underlying this shift in practice is that fast-growing woody competitors that readily establish from the seed bank following a harvest are subsequently controlled by the broadcast herbicide application and further recruitment is limited relative to the initial post-disturbance pulse given the depletion in the seed bank [28]. Additionally, this approach is pragmatic, as unforeseen harvesting delays (e.g., weather, market volatility) may provide sufficient time for the re-establishment of interfering vegetation where pre-harvest herbicides were applied. Nevertheless, this approach may be of limited utility when the competitor is a species like birch, where seedling recruitment from seeds comes not from an on-site seed bank that may be depleted, but rather from a continuous source of wind-dispersed seeds.
In this study, we used data from nearly two dozen managed stands across multiple ownerships to examine how the application of the predominant even-aged regeneration prescription in the Allegheny hardwood forest type, the shelterwood seed cut and herbicide combination, influences birch establishment. If post-harvest conditions promote birch establishment and growth, then we predict the shelterwood-herbicide sequence will reduce birch seedling densities and height development relative to the herbicide-shelterwood sequence. Moreover, as birch recruitment and growth are known to be influenced by site productivity, browse pressure, and propagule supply, we further predict birch seedling densities and height will be greater in sites with increased soil nitrogen and moisture, greater light availability, low deer browsing, and where mature birch exist in the canopy.

Study Area
We conducted our study at 22 northern hardwood forest sites distributed throughout a 6500-km 2 area of northern Pennsylvania, USA. The chosen sites had a mean elevation of 605 m, a humid, temperate climate with average daily temperatures of 9 • C, and average annual precipitation of 1067 mm [29]. Fecal pellet surveys at these sites conducted between 2014-2016 found that deer densities averaged 7.04 ± 0.48 deer/km 2 , levels which are moderate given historic trends in the region [30].
At all sites, managers conducted the initial cut (i.e., the preparatory cut) of a shelterwood sequence to reduce stand relative density (i.e., <75% relative density) and applied broadcast herbicides (tank mix of glyphosate and sulfometuron methyl; [31]) to control interfering plant species and reinitiate seedling recruitment. Herbicide application was accomplished using mist-blowers mounted on skidders.
Where reported, application rates varied from 1.8-3.8 kg/ha of active ingredient for glyphosate and 0.14-0.21 kg/ha of sulfometuron methyl. Twelve sites received a herbicide-shelterwood sequence (H-S) and ten received a shelterwood-herbicide sequence (S-H; see Appendix A). Preparatory cuts in the H-S sequence occurred 0-2 years prior to vegetation monitoring, while herbicide applications in the S-H sequence sites occurred two years prior to vegetation monitoring. Within each site, we established two 0.42-ha (60 × 70 m) plots and randomly assigned one plot a deer exclosure (fence) treatment while the other served as an unfenced control. Exclosure construction was completed by September 2013, one year prior to vegetation monitoring.
Additionally, we calculated various other site variables for each plot. Site productivity was estimated using Iverson et al.'s [32] integrated moisture index (IMI). The IMI combines GIS-derived topographic and soil features of the landscape that govern moisture availability (i.e., direct solar radiation, slope position, curvature, and soil water holding capacity) into a single index of relative moisture and site productivity. We assessed nitrogen availability for each site using soil samples from the A, B, O, and E horizons. Ammonium nitrogen (NH 4 ) and nitrate (NO 3 ) concentration estimates were converted to total nitrogen (mg/kg) using their molar ratios. Finally, we conducted post-harvest variable radius stem surveys using a 10-factor prism to determine whether mature, potentially seed-bearing, birch were present in the overstory of each plot and to calculate overstory relative density. Relative density is an estimate of overstory crowding that integrates both tree size and species identity and serves as our surrogate for light availability [33].

Vegetation Surveys
We surveyed tree seedling regeneration in the 2014, 2015, and 2017 growing seasons at each plot. Sampling occurred in 1-m 2 circular subplots located at the nodes of a 10 m × 10 m grid established within the central 0.2 ha (40 m × 50 m) area of each plot. We recorded densities of established seedlings (i.e., not new germinants, >5 cm in height) and the height of the tallest individual, by species. Due to limitations in personnel, we sampled the 15 odd-numbered subplots in 2014 and 2017 and all 30 subplots in 2015. Mean density (stems/m 2 ) and average height (cm) were calculated for each site, plot, and year.

Statistical Analyses
We employed an analysis of covariance approach using generalized linear mixed models (PROC GLIMMIX; [34]) to examine the factors that predispose stands to high birch seedling and sapling dominance. To avoid model overfitting when testing multiple covariates, we followed Littell and colleagues' [35] guidelines, which recommend a sequential model-building approach in which the form of each covariate is evaluated independently prior to entering into a final model. This approach first investigated whether the slopes were equal across treatments (i.e., treatment × covariate interaction). If the interaction was deemed non-significant, we then examined whether a common slopes model (i.e., simple covariate) was appropriate. Alternatively, if slopes were unequal then treatment effects were tested at a minimum of three levels of the covariate. For our purposes, we chose the 10th percentile, the median, and the 90th percentile for each chosen covariate, when appropriate. Each potential covariate was tested against our two fixed silvicultural treatment effects: silvicultural sequence (H-S versus S-H) and deer browsing (ambient versus exclosure). Only significant covariates or treatment × covariate effects were included in a final model. The tested co-variables were as follows: birch seed source in the plot (binary: 0 versus 1), site productivity/moisture (continuous: IMI), canopy openness (continuous: relative density), and total nitrogen (continuous). One observation was removed from the input dataset for birch seedling densities, as it was identified as an outlier (i.e., large residual) that exerted high leverage (PROC ROBUSTREG; [34]) Response variables included both established birch seedling density (# individuals/m 2 ) and birch seedling height (cm). Additionally, to explore whether any possible reductions in birch density or height benefited the regeneration of other species, we examined silvicultural sequence and browsing effects on the average stem density and height of all co-occurring species of higher economic regional importance (i.e., desirable). Desirable species in our surveys included Acer rubrum L., A. saccharum Marshall, Fraxinus americana L., Liriodendron tulipifera L., Magnolia acuminata L., Prunus serotina Ehrh., and Quercus rubra L.
Birch density and height were analyzed using a repeated measures analysis of covariance randomized incomplete block factorial design. For both models, sequence, enclosure, year, and the associated interactions were fixed effects, and site was a random effect. Models on desirable species density and height were tested on 2017 data only. Therefore, we used a randomized incomplete factorial design with these two models. Sequence, enclosure, and the associated interaction were fixed effects in these models, and site was the random effect.
Seedling heights were modeled with a lognormal distribution with the identity link function, and average seedling densities were modelled with a gamma distribution using the log link function. Correlations between years were modelled using a spatial power covariance structure. We used the Kenward-Rogers denominator degrees of freedom adjustment method for each model. In addition, normality was statistically tested using the Shapiro-Wilk test, and homogeneity of variance was tested with Levene's test [36]. For significant treatment × year effects, we compared treatments within years using the SLICE option of the LSMEANS statement, and any treatment × covariate tests were accomplished using the LSMESTIMATE function with the AT = option for the three levels of the covariate. We used a critical value of p = 0.10 as significant and employed a Bonferroni correction when examining multiple comparisons.

Results
Initial (2014) birch seedling densities were 71% lower in areas treated with the S-H sequence relative to the H-S sequence; however, seedling densities were statistically equivalent between treatments by 2015 and grew even more equitable by 2017 (sequence × year effect; Table 1; Figure 1). Within the H-S sequence areas, birch established at high densities only where residual birch existed in the overstory (Table 1; Figure 2). In fact, areas with the H-S sequence that lacked mature birch in the canopy had birch seedling densities comparable to S-H areas ( Figure 2). Overall, deer browsing had no effect on birch seedling densities. Excluding browsers in the H-S treatment did result in seedling densities that were two-to three-fold greater than other treatments (exclosure × sequence effect; Table 1; Figure 3); however, these initially significant differences (p < 0.10) were non-significant following the adjustment for multiple comparisons. Site productivity (IMI) and relative density were identified as potential covariates in the exploratory model-building phase for birch density; however, neither was significant when entered into the full model.    Birch seedlings grew over time, attaining an average height of 80.1 ± 11.9 cm by 2017, regardless of silvicultural treatment sequences or browsing pressure (Table 2). Deer browsing reduced birch seedling heights by 29% when averaged across all sample periods (Table 2, Figure 4). By 2017, birch seedlings were 61.7 ± 11.3 and 103.92 ± 19.5 cm in control and fenced areas, respectively. Birch seedlings attained larger heights as site productivity (IMI) increased, irrespective of treatments, but the relationship had high variability (F 1,105 = 13.27; R 2 = 0.11; Table 2; Figure 5). No other potential covariates were identified in the model-building phase as important factors driving the birch height response.   Although combined seedling densities of desirable woody regeneration were, on average, nearly an order of magnitude greater than the birches in 2017 (9.5 ± 1.0 stems/m 2 versus 1.2 ± 0.3 stems/m 2 , respectively), these were unaffected by either silvicultural sequence or browsing treatments (Appendix B). Similarly, desirable seedling heights were unaffected by treatments (Appendix B). The seedling heights of desirable species were uniformly short (Mean height in 2017: 29.3 ± 3.2 cm) among treatments and averaged less than half the height attained by birch stems.

Discussion
Our study demonstrates that a post-shelterwood herbicide application (S-H) initially reduced birch recruitment into stands. Mechanistically, post-shelterwood harvest applications coincide with the period during which birch recruitment is greatest, and therefore exert the strongest effect. Moreover, all sites in this study had sulfometuron methyl combined with glyphosate. Although glyphosate does not have residual activity, sulfometuron methyl does possess residual soil activity, which lasts for several months [37,38] and thus may further limit birch recruitment. Although mist-blown herbicide effectiveness on interfering vegetation has been tested in pre-(e.g., [39]) and post-harvest contexts (e.g., [40]), to our knowledge, this is the first study to explicitly examine how variation in the timing of the harvest and herbicide sequence affects woody regeneration dynamics. In an experiment designed to test herbicide effects on plant diversity, Ristau and colleagues [41] had ten sites evenly split into either the H-S and S-H sequences. Although they found sequence did not affect overall plant diversity or composition, closer examination of the data suggests the H-S sequence promotes birch (H-S: 5180 ± 1986 stems/ha versus S-H: 3070 ± 934 stems/ha, respectively; Ristau unpublished data). Our work also complements and expands findings of Kelty and Nyland [42] and Ray and colleagues [43], who found pre-shelterwood harvest herbicide applications promoted yellow birch (B. allegheniensis) in northern hardwood stands. Thus, our research supports management recommendations that control of competing woody vegetation may be most effective when the herbicide follows the harvest [28,40,44,45].
Our results provide strong evidence that browsing reduces average birch seedling height development and potentially limits the increase in birch seedling density observed in the H-S sequence. These findings concur with prior work that found browsing limits yellow and sweet birch [18,[46][47][48]. Nevertheless, even in areas exposed to browsing birch seedlings were, on average, taller than the desirable regeneration by 2017. These results suggest birch seedlings may be either less preferred by deer, relatively resilient to browsing, or both. Indeed, Bressette and colleagues [17] rated both yellow and sweet birch as minimally sensitive to browsing, whereas some of the species we categorized as desirable regeneration (e.g., maples, northern red oak, and white ash (Fraxinus americana L.)) are often much more sensitive to browsing (see also [49]). This may explain why birches often dominate post-disturbance regeneration even under average deer densities much higher than those in our study areas (11-18 deer/km 2 ; [18,19]). These findings demonstrate the challenge of sustaining diverse tree regeneration in Allegheny Plateau forests, as high browse impacts severely affect the regeneration of most tree species [18], and in present-day forests the fast-growing and relatively browse-insensitive birches dominate at moderate or even low (i.e., fences) browse impacts.
Our results also demonstrate that birch recruitment success is substantially enhanced by the presence of residual mature birch trees in the post-shelterwood canopy, specifically following the H-S sequence. The likely explanation is that these residual trees shed seeds that successfully germinated under the favorable seedbed conditions created by the harvest [50] and were not subsequently controlled by herbicide. These results are surprising given that a leading forest dynamics model (SORTIE; [51]) finds that birch is not strongly recruitment limited as a consequence of being a prolific seeder with extremely broad dispersal capabilities [14,52]. Our findings suggest that despite their broad dispersal ability, seed rain can still be positively influenced by the proximity of mature, seed-bearing trees (see [53]). Culling any mature birch trees from stands during the initial shelterwood cut will greatly limit, although not eliminate, birch recruitment.
Although birch grew taller in more productive sites, this relationship exhibited considerable variation. In our work, productivity was defined by a GIS-derived integrated moisture index [54]. This index is heavily weighted towards topographic and soil features of the landscape that govern moisture availability (i.e., direct solar radiation, slope position, curvature, and soil water-holding capacity). Thus, our findings suggest birch does have an affinity to sites with greater moisture but can nonetheless thrive across a wide variety of sites. Our results broadly concur with the existing literature that finds sweet and yellow birch thrive in cooler, moist sites at higher elevations with well-drained soils [10,55]. This versatility with regard to soil moisture relationships supports existing findings documenting the encroachment and survival of birch on drier sites, including ridgetops, even following droughts [56,57].
Contrary to expectations, neither overstory relative density nor total nitrogen were identified as factors influencing birch establishment dynamics. The lack of a relative density effect is likely attributable to birch exhibiting a fairly labile response to variation in light conditions [20,58] coupled with post-shelterwood residual relative densities in our study that are intended to increase understory light to stimulate regeneration. Hence, in these sites, light is not limiting, particularly for taxa that perform well under a gamut of light levels. The lack of a nitrogen effect may be attributable to at least two causes. First, given that harvesting in the S-H stands occurred several years (Mean: 4.1 years) prior to the initial (2014) survey, whereas harvesting in the H-S stands occurred much closer to the sampling period (Mean: 1.4 years), this may mean that nitrogen is temporally confounded with treatment sequence, as overstory disturbances typically create an ephemeral pulse in soil ammonium and nitrate [22,59]. This explanation is supported by the fact that nitrogen concentrations in the H-S sequence stands were approximately 33% greater than those with the S-H sequence (20.6 versus 15.1 ppm or mg/kg, respectively). Moreover, soil nitrogen may not be a key limiting resource for birch, particularly under higher light conditions [20,60]. While it is important to note that these caveats constrain our ability to test the importance of these variables across broader resource gradients, our goal was to assess their effects under real-world, operational management scenarios in this region. Hence, we caution readers not to extrapolate our findings to infer that, in general, canopy openness and nitrogen have no effect on birch establishment and growth.

Conclusions
Our results provide empirical confirmation that the shelterwood-herbicide sequence initially limits the establishment of birch seedlings relative to the herbicide-shelterwood sequence. This approach provides only a temporary reprieve, as within two years birch densities recovered to levels comparable to the herbicide-shelterwood. More importantly, given that in contemporary Allegheny and northern hardwood forests, birch grows faster than most of its competitors following overstory disturbance [4,46], birch rapidly overtop most other species. In fact, five growing seasons post-herbicide, birch seedlings were already over twice as tall as the desirable regeneration. Birch recruitment and growth were also influenced by propagule supply and deer browsing. Hence, shelterwood harvests are at greater risk of birch dominance when there are residual birch in the stands as possible seed sources and where browse pressure is low. Nevertheless, we caution that even in sites where birch seed source is lacking, where modest browse pressure limits birch growth, and when managers use a post-shelterwood herbicide (H-S) sequence to control birch, there is no guarantee that birch will not become established and grow to be the tallest competitors in the regenerating community. Thus, follow-up treatments such as broadcast herbicide or targeted cleanings (e.g., stem injection, basal spray) that remove undesirable, competing woody vegetation may be necessary to promote desirable regeneration [61].  Table A1. Site characteristics of the 22 experimental sites. Sequence denotes whether sites received a herbicide preceding the shelterwood harvest (H-S) or if the herbicide was applied post-harvest (S-H). Browsing refers to the 0.42-ha fenced and unfenced areas. IMI is the integrated moisture index, a metric of soil moisture availability and productivity based on GIS-derived metrics [32]. Basal area is in m 2 per hectare. Relative density is a measure of overstory crowding [33]. Birch seed source (Absent/Present) denotes whether there were any residual canopy birch (B. allegheniensis or B. lenta) in the canopy post-shelterwood harvest. We also present elevation (meters), slope (degrees), and aspect for each site.