Vegetative and Edaphic Responses in a Northern Mixed Conifer Forest Three Decades after Harvest and Fire: Implications for Managing Regeneration and Carbon and Nitrogen Pools

: Research Highlights : This experiment compares a range of combinations of harvest, prescribed ﬁre, and wildﬁre. Leveraging a 30-year-old forest management-driven experiment, we explored the recovery of woody species composition, regeneration of the charismatic forest tree species Larix occidentalis Nutt., and vegetation and soil carbon (C) and nitrogen (N) pools. Background and Objectives : Initiated in 1967, this experiment intended to explore combinations of habitat type phases and prescribed ﬁre severity toward supporting regeneration of L. occidentalis. At onset of the experiment, a wildﬁre a ﬀ ected a portion of the 60 research plots, allowing for additional study. Our objective was to better understand silvicultural practices to support L. occidentalis regeneration and to better understand the subsequent impacts of silvicultural practices on C and N pools within the vegetation and soil. Materials and Methods : We categorized disturbance severity based on loss of forest ﬂoor depth; 11 categories were deﬁned, including controls for the two habitat type phases involved. We collected abundance, biomass, and C and N concentrations for the herbaceous layer, shrubs, and trees using nested quadrats and 6 to 10 experimental units per disturbance category plot. Moreover, we systematically sampled woody residue from transects, and forest ﬂoor, soil wood, and mineral soil with a systematic grid of 16 soil cores per disturbance category plot. Results : We found that (1) disturbance severity a ﬀ ected shrub species richness, diversity, and evenness within habitat type phases; (2) L. occidentalis regenerates when ﬁre is part of the disturbance; (3) N-ﬁxing shrub species were more diverse in the hotter, drier plots; (4) recovery levels of C and N pools within the soil had surpassed or were closer to pre-disturbance levels than pools within the vegetation. Conclusions : We conﬁrm that L. occidentalis regeneration and a diverse suite of understory shrub species can be supported by harvest and prescribed ﬁre, particularly in southern and western aspects. We also conclude that these methods can regenerate L. occidentalis in cooler, moister sites, which may be important as this species’ climate niche shifts with climate change. D.S.P.-D. and M.F.J.; investigation: D.S.P.-D., M.F.J., and R.K.D.; data curation, D.S.P.-D., M.F.J., and R.K.D.; formal analysis, R.K.D.; visualizations, R.K.D.; writing—original draft, R.K.D.; writing—review and editing, D.S.P.-D., M.F.J., and R.K.D. All authors read agreed the manuscript.


Introduction
Forest ecosystems are dynamic, driven by disturbance that occurs commonly rather than rarely, and that encompasses a range of spatial, temporal, and severity combinations [1]. The structure, composition, and spatial arrangement of many forested ecosystems, from tropical to boreal, is a  (60 plots total), with the intention to examine the effects of combinations of harvesting and prescribed fire on conifer regeneration, particularly L. occidentalis. The 1967 wildfire burned or re-burned some of the original plots. In 1996, we selected proximate stands of undisturbed forest in the Menziesia ferruginea (MEFE) and Xerophyllum tenax (XETE) phases within the Abies lasiocarpa/Clintonia uniflora habitat type and relocated 9 of the original study plots with various combinations of habitat type phase, anthropogenic activity (undisturbed or harvested stands with or without prescribed fire), and a stochastic event (wildfire) to obtain 11 distinct disturbance categories. See Table 1 for plot identifications and descriptions. Adapted from [46] (p. 6) and [51] (p. 6). Xerophyllum tenax (XETE) phases within the Abies lasiocarpa/Clintonia uniflora habitat type and relocated 9 of the original study plots with various combinations of habitat type phase, anthropogenic activity (undisturbed or harvested stands with or without prescribed fire), and a stochastic event (wildfire) to obtain 11 distinct disturbance categories. See Table 1 for plot identifications and descriptions. Adapted from [46] (p. 6) and [51] (p. 6).  30 48 a Phase within Abies lasiocarpa/Clintonia uniflora habitat type [50]; b [45]; c [51]; d [18]. See methods [45,52]. Original water mass = 3 kg. e DeByle's unpublished data, unless otherwise noted. f Undisturbed forest adjoining East 2, North 8, and North 12. g Undisturbed forest adjoining West 1, South 1, and South 9. h The prescribed fire reduced forest floor by 58% (DeByle's unpublished data) and the wildfire further reduced it to a total of 84% [53]. i Forest floor reductions were "about 100%" [53] and "100%" [19] (p. 67, Table 1, Wildfire). j Stands burned 23 August 1967 [18] and are listed in [45] (p. 45) under "burned by prescription". k Water loss value for the prescribed fire (see footnote d above). We assume that these units had been prepared (i.e., filled water cans installed) for prescribed burning but burned in the wildfire, and thus the wildfire was classified as the prescribed fire [45]. Given the fire danger (buildup index = 239, classed as "extreme" by the National Fire Danger Rating System [54] and [45] (p. 48, Appendix D: Weather Summaries)) and the unfavorable wind conditions (see the Thursday, 24 August 1967, front page account titled "Fire Erupted Tuesday Afternoon" in The Daily Inter Lake (Kalispell, Montana) newspaper), prescribed fire would likely have been avoided. k Water loss value for the prescribed fire (see footnote d above).

Original Study Design
During 1966 and 1967, fifteen 4.05 ha plots were established on each of the four cardinal aspects (60 plots total) to examine the effects of combinations of harvesting and prescribed fire on conifer regeneration, especially for L. occidentalis (Figure 1). The original stands, composed mainly of 200-to 250-year-old L. occidentalis (26%), Pseudotsuga menziesii (Mirb.) Franco var. glauca (Mayr) Franco (31%), and Picea engelmannii (31%) trees in nearly equal volumes, averaged 143 m 3 ha −1 [45]. These stands reflected little to no post-European-settlement disturbance and were either left undisturbed or clearcut. Clearcutting was done by hand; the trees were directionally felled, limbs were removed, and logs were cut to length. Logs were dragged to landings using a stationary tractor-mounted, two-winch cable logging system along designated "skid trails" that were at least 200 m apart. Logging residue was spread evenly across each plot by the hand crew in combination with felling of unmerchantable trees, so that fuel loads would be as uniform as possible [45]. Two levels of prescribed fire (yes or no) were employed during different seasons to achieve various fire severities [45]. The average post-harvest fuel loads were 234 Mg ha −1 , with about 88% being coarse debris >10 cm in diameter [55]. To create a mostly uphill-leading prescribed fire, ignition began in the late afternoon or early evening at the upper edge of a plot, then downslope along the sides, and finally across the bottom [45]. A wildfire on 23 August 1967 burned or re-burned many of the plots, creating a mosaic of unplanned and anthropogenic impacts across the watershed.

Field Sampling
In both habitat type phases and proximate the original study plots, we identified stands of undisturbed forest. Using combinations of habitat type phases (MEFE or XETE) and anthropogenic (undisturbed or harvested stands with or without prescribed fire) and stochastic (wildfire) history, we defined 11 distinct disturbance categories. Each category was defined by relative fire severity, which was based on loss of forest floor measured immediately pre-and post-disturbance (1967 or 1968) rather than on water loss using water can analogs (Table 1; [52]) because forest floor data were more complete. Each original 1967 study plot represented a single disturbance category, except for South 12 and South 13, which we combined into a single disturbance category (South 12/13). Within each disturbance category, we systematically placed 6 to 30 experimental units to measure a variety of vegetative metrics and at least 16 soil cores in a systematic grid (Table 1). Between the initial harvest and fire events (1967)(1968)) and our measurements (1996)(1997)(1998), some of the disturbance categories had portions treated with other silvicultural practices (i.e., tree planting, thinning [56]); we avoided these treated portions.

Vegetation
Aboveground Each experimental unit was a nested quadrat [57]. Within the largest quadrat (15 m × 15 m) we situated a 4 m × 4 m quadrat on the downslope corner (northeast for the northern aspect, southeast for the eastern aspect, southwest for the southern aspect, and northwest for the western aspect). Similarly, a 1 m × 1 m quadrat was placed on the downslope corner of the 4 m × 4 m quadrat, with a second 1 m × 1 m quadrat established in the diagonally opposite corner of the 15 m × 15 m quadrat ( Figure 2).
Within each 15 m × 15 m quadrat, we measured trees (basal diameter > 5 cm) by species (Table 2) at diameter breast height (DBH; 137 cm above ground level). Standing dead trees were measured and identified to species. Within each 4 m × 4 m quadrat, we measured basal diameters of tree seedlings (<5 cm basal diameter; hereafter "seedlings"); dead seedlings were not identified to species. We also measured shrubs ≥ 50 cm in height by species (see Table 2 for shrub and tree nomenclature) 2 cm above ground level at seven intervals: ≤0.5, 0.51-1.0, 1.01-1.5, 1.51-2.0, 2.01-3.0, 3.01-5.0, and >5.0 [58]. Dead shrubs were measured but not identified. We removed a branch (≥5 cm in length) from the upper crown of each seedling and shrub (including dead) for further analysis. In each 1 m × 1 m quadrat, all aboveground vegetation was harvested (except lichens and moss, shrubs taller than 50 cm), hereafter named "herbaceous" vegetation. and shrub (including dead) for further analysis. In each 1 m × 1 m quadrat, all aboveground vegetation was harvested (except lichens and moss, shrubs taller than 50 cm), hereafter named "herbaceous" vegetation.    [59]. y Stature classification [58]. For each species of seedling and shrub, branch samples were composited by experimental units; samples from half the units formed one sample and samples from the remaining units forming the second sample. Samples were further composited by the basal diameters (≤2 or >2 cm) of the original plants. Finally, we segregated branch samples into foliage and wood. Therefore, the maximum sample number for any given species within a disturbance category was eight (2 experimental unit composites x 2 basal diameter categories x 2 tissue types). Within each experimental unit, we composited the herbaceous material from both 1 m × 1 m quadrats.

Belowground
We estimated the coarse root biomass using Equation (1): where the constant 0.26 is the root-to-shoot ratio of temperate tree species [60]. Fine root biomass was the average oven-dry weight of material collected in each forest floor and mineral soil core multiplied by the number of cores (area basis) per hectare (see Section 2.3.3).

Woody Residue
To estimate the amount of sound and rotten woody residue ≥ 0.6 cm in diameter (woody residue; WR) on the soil surface, we used wood classification categories [61] and two randomly placed 10.7 m linear transects within each experimental unit. Woody debris < 0.6 cm in diameter (1-h fuel) was sampled as part of the forest floor (see next section). Random samples of WR (across diameter classes) were collected for C and N analyses. We corrected transect biomass using an average of the appropriate specific gravity values and to estimate volume on a per hectare basis [61].

Soil
Within each disturbance plot we collected surface organic horizons and mineral soil from a grid of at least 16 systematically located points (Table 1). At each point, surface organic horizon samples (inclusive of the O i , O e , and O a horizons) were removed from within a 10 cm diameter circle and the depths were recorded. Surface organic horizons included all dead organic matter in the O horizons and any woody debris < 0.6 cm in diameter. Surface mineral soil (0-30 cm) samples were obtained with a large (10 cm diameter × 30 cm depth) corer [62] that allowed us to gather representative soil samples, including large coarse fragments [63]. A total of 160 surface organic horizons and mineral soil samples were taken from the original study plots (n = 10 plots × 16 cores), and 96 were taken from adjacent undisturbed control stands (n = 2 habitat type phases × 3 stands × 16 cores). Soil wood and fine roots found within the mineral soil were collected for processing. Soil wood is highly decomposed woody residue of decay class 5 [64] that is covered by the forest floor material or within the mineral soil and not measured in woody residue transects [35].

Vegetation
Tree and seedling biomass levels per hectare were estimated using the Fire and Fuels Extension (FFE) [65] of the Forest Vegetation Simulator (FVS) [66] with the Northern Idaho/Inland Empire geographic variant [67]. Dead seedlings were analyzed as "other softwoods." We assumed that snags had branches. Aboveground biomass (g) for shrub components (total and foliage) within each experimental unit was approximated with Brown's [58] regression components and Equations (2) and (3): Forests 2020, 11, 1040 9 of 35 where i is the shrub species, D is the basal diameter (cm), and the regression parameters a and b stand for the total (t) and foliage (f ) biomass. We used Brown's [58] diameter class midpoints 0.34, 0.40, and 0.44 for low, medium, and high shrub groups (  (Table 2). We estimated stem biomass (including branches) with Equation (4), and assumed that Biomass stem was the best estimate for dead shrub biomass. For dead shrubs, we used the average of the regression parameters for all shrub species: All tree, shrub, herbaceous, woody residue, and fine root samples were oven-dried at 60 • C to a constant weight, then a sub-sample was finely ground with an 8000D Mixer/Mill (Spex SamplePrep, Metuchen, NJ, USA) and analyzed on a LECO TruSpec CN analyzer (Leco Corp., St. Joseph, MI, USA) to attain total C and total N concentrations. Estimates of biomass ha −1 were multiplied by C and N concentrations to calculate pool sizes.

Species Diversity
We explored species diversity and evenness using Shannon diversity functions [68]. Equation (5) describes species diversity and Equation (6) describes evenness; that is, the equitable distribution of individuals within the plot among species: where H is the diversity index, p i is the proportion of the entire population made up of species i, S is the total number of species within an experimental unit (i.e., richness), and J is the evenness index.

Soil
All forest floor and soil wood samples were removed from the mineral soil and dried at 60 • C. Forest floor and soil wood samples were then processed and analyzed as described above for vegetation. Mineral soil samples were dried to a constant weight at 105 • C, weighed, finely ground to pass a 0.04 mm mesh, and analyzed for C and N concentrations (described above).
The dry mass levels of forest floor, soil wood, and mineral soil in each core were extrapolated to a per hectare basis. Mineral soil C, N, and organic matter contents were corrected for coarse rock fragment contents and re-extrapolated to a per hectare basis using fine fraction bulk density [63,69]. We did not analyze the coarse rock fragments (>2 mm), which have been found to contain appreciable amounts of C and N in some soils [70,71]. The total and fine bulk densities and gravimetric and volumetric rock contents were calculated following Page-Dumroese et al. [63].

Statistical Analysis
For biomass, C and N concentrations, and C and N contents of pool components (e.g., aboveground shrub biomass, coarse roots, mineral soil), pairwise comparisons (α = 0.1) of disturbance categories were computed using the R language [72] and multcomp library [73] after adjustment for family-wise error rate using the Tukey-Kramer method [74]. We prepared a principal coordinates analysis (PCoA) using the Bray-Curtis dissimilarity method to visualize similarities in shrub communities across the disturbance categories [75,76]. Because shrub diversity (H ) and species evenness (J ) within each disturbance category are responses from the same set of variables, we used MANOVA to first test for differences between habitat type phases, and secondly for differences within each habitat type phase among disturbance categories. The resulting p-values for Pillai's trace revealed no significant difference between habitat type phases (p = 0.3568). Thus, we used separate ANOVAs to assess whether more of the variability seen in the preceding MANOVA was due to H or to J . Because data failed the univariate normality test, we separated means using the Kruskal-Wallis method [77]. For tissue N concentrations, we first explored whether any data could be pooled by calculating a mean for each species x tissue type (wood or foliage) x stem basal diameter (<2 or >2 cm) combination and then comparing weighted means among plant types using equal variance t-tests, but found that N concentration was usually highly significant (p < 0.0001) across metrics, so that no pooling combinations of tree, shrub, or N-fixing species could be justified. We then used equal variance t-tests to compare tissue N concentrations across stem basal diameters and ANOVA to compare concentrations within tissues for seedling, shrub, and N-fixing shrub species, separating means with Tukey's Honestly Significant Difference test when appropriate.

MEFE
For this cooler, moister phase, disturbance was only a result of harvesting and prescribed fire. All three disturbance plots were harvested, but East 3 and North 12 also experienced prescribed fire. As shown in Table 1, 42% of the forest floor persisted in East 3 after the prescribed fire in August 1968, while only 9% remained in North 12, as this area was burned with a hotter fire in August 1967.

XETE
Disturbance on the warmer, drier XETE phase was a function of harvesting and prescribed fire, but with the additional impact of wildfire. Immediately after harvest, five of the six XETE plots burned with combinations of prescribed fire and wildfire; the disturbance level increased with fire severity, as measured by the loss of forest floor ( Table 1). The lowest disturbance was to the harvest-only South 11 plot. West 10 and South 8 were harvested and burned with a low severity prescribed fire that removed 16% and 68% of their forest floors, respectively. Fifteen days after the prescribed fire, South 8 was re-burned by wildfire, further reducing its forest floor to 16% of the original level. South 7 and West 2 were harvested, but before they could be treated with prescribed fire, were burned by wildfire that removed 72% and 98% of their forest floors, respectively. South 12 and 13 were not harvested but were burned in the wildfire that resulted in almost total loss of the forest floor.

Herbaceous and Shrubs
Across the experimental landscape, we observed 18 species of shrubs (Table 2); the nine most common "species" (including dead stems as a species) in terms of relative abundance were Alnus viridis ssp. sinuata, "dead" material, and Vaccinium membranaceum (15% each); M. ferruginea (11%); P. myrsinites (9%); S. scouleriana (8%); and Ceanothus velutinus, Shepherdia canadensis, and Spiraea betulifolia (7% each); these accounted for 93% of all stems ha −1 and 96% of the total biomass (Supplementary Materials  Table S1). Three species were uncommon, recorded in low numbers (<125 stems ha −1 ), and only for single disturbance categories: C. canadensis, B. repens, and Juniperus communis. The remaining seven species (Acer glabrum, Amelanchier alnifolia, Lonicera utahensis, R. lacustre, R. gymnocarpa, Rubus parviflorus, and Symphoricarpos albus) occurred more frequently (present on 5 to 7 categories) but contributed < 1% each to the total number of plants observed across all 11 disturbance categories (Table S1). Neither shrub H nor J were significantly affected by habitat type phase (p = 0.3568), with less variability observed in H (p = 0.7316) than in J (p = 0.4819) ( Figure 3). 12 "species" (including dead stems as a species) in terms of relative abundance were Alnus viridis ssp. sinuata, "dead" material, and Vaccinium membranaceum (15% each); M. ferruginea (11%); P. myrsinites (9%); S. scouleriana (8%); and Ceanothus velutinus, Shepherdia canadensis, and Spiraea betulifolia (7% each); these accounted for 93% of all stems ha −1 and 96% of the total biomass (Supplementary Materials Table S1). Three species were uncommon, recorded in low numbers (<125 stems ha −1 ), and only for single disturbance categories: C. canadensis, B. repens, and Juniperus communis. The remaining seven species (Acer glabrum, Amelanchier alnifolia, Lonicera utahensis, R. lacustre, R. gymnocarpa, Rubus parviflorus, and Symphoricarpos albus) occurred more frequently (present on 5 to 7 categories) but contributed < 1% each to the total number of plants observed across all 11 disturbance categories (Table S1). Neither shrub H′ nor J′ were significantly affected by habitat type phase (p = 0.3568), with less variability observed in H′ (p = 0.7316) than in J′ (p = 0.4819) ( Figure 3).  Figure 3. Species richness (S) and box plots for Shannon functions of species diversity (H′) and evenness (J′) for the 11 disturbance categories (see Table 1 Table 1  Although relative abundance and relative biomass across most species usually followed the same trend (i.e., more relative abundance generated more relative biomass), this did not consistently translate for low stature shrubs (listed in Table 2). For example, on the wildfire-only plot (South 12/13), S. betulifolia accounted for 21% of all stems ha −1 (Table S1) (13,646 stems ha −1 ; Table 3) but < 1% of the biomass (Table S1). Often, many low-stature shrubs produced less biomass relative to a few high stature shrubs, e.g., V. membranaceum accounted for 47% of all stems ha −1 and 21% of the biomass on the harvested and prescribed fire plot East 3, whereas A. v. sinuata accounted for 3% of all stems ha −1 but 40% of the biomass. In addition, the high relative abundance of low stature shrubs on the undisturbed plots was often dwarfed in terms of relative biomass by standing dead material, e.g., V. membranaceum accounted for 41% of all stems ha −1 and 11% of the biomass on the MEFE phase, whereas the "dead" category accounted for 7% of all stems ha −1 but 34% of the biomass (Table S1). Table 3. Thirty-year post-treatment measures of the total live stems ha −1 of shrubs, tree seedlings, and trees for the two phases (Menziesia ferruginea, MEFE; Xerophyllum tenax, XETE) of the Abies lasiocarpa/Clintonia uniflora habitat type [50] and the subsequent 11 disturbance categories for the Miller Creek Demonstration Forest, Montana, USA (see Table 1 for descriptions). For treatments: H = harvested; P = prescribed fire; W = wildfire; C = control. N-fixing plants shown in bold.  (Figure 3) and varied species composition among MEFE plots (Figure 4). On the harvest-only plot East 2, being mainly V. membranaceum, M. ferruginea, and A. v. sinuata ( Figure 5 and Table S1), species evenness (J ) was lower than on the harvested East 3 that had prescribed fire. The occurrence of fire on East 3 (and the 57% loss of forest floor) shifted the shrub community composition ( Figure 4) and increased S to 10 ( Figure 3), although abundance and biomass were reduced for V. membranaceum and M. ferruginea and increased for N-fixing A. v. sinuata ( Figure 5 and Table S1). The combination of harvest with a more severe prescribed fire on North 12 (91% loss of forest floor) further increased species richness (S = 13), but because V. membranaceum and M. ferruginea were nearly eliminated and A. v. sinuata and S. scouleriana predominated (nearly 90% of the biomass; Figure 5 and Table S1), J was significantly lower than when harvested and burned by prescription in East 3 and the undisturbed plots (UM). The total shrub biomass values significantly increased with fire severity, with the lowest value observed in the harvest-only plot East 2 (no loss of forest floor) with 5.  (Table S4). The relative biomass of the "dead" material was low across the disturbed plots (≤ 6%) compared to the undisturbed plots (UM; 34%) The amount of herbaceous material was low across all MEFE plots (Table 4).   Table 1 for descriptions).   Table 1 for descriptions).
Polygons delineate the range of data for the 11 disturbance categories in the Miller Creek Demonstration Forest (see Table 1 for descriptions).    Table 1 for disturbance category descriptions.

XETE
Without fire, the herbaceous biomass was 2-to 5-fold greater than when fire occurred (Table 4). For shrubs, harvesting alone stimulated shrub development compared to the control; the amount of shrub stems ha −1 increased from 1584 in the undisturbed plots (UX) to 28,750 in the harvest-only South 11 plot (Table S1), while aboveground shrub biomass significantly increased from 0.1 (0.3) to 2.6 (2.4) Mg ha −1 (Table S4). Moreover, H increased by almost 2-fold and J increased by 50% (Figure 3). Fire further stimulated shrub occurrence; relative abundance and relative biomass each increased from about 30% to 2.5-fold compared to the harvest-only South 11 (Table S1), while aboveground shrub biomass in non-harvested South 12/13, the most severely disturbed plot with 100% of the forest floor consumed by wildfire, was significantly different from harvest-only South 11 (7.2 (6.0) vs. 2.6 (2.4) Mg ha −1 ) (Table S4). Salix scouleriana, absent on undisturbed plots (UX) and harvest-only South 11, became abundant on southern aspects after fire, but its total abundance decreased with increasing fire severity ( Figure 6 and Table 3, and Table S1). Although the principle component analysis indicated that XEFE plots (with the exception of West 10, which was harvested and had prescribed fire) were quite similar in terms of composition (Figure 4), only four species were detected (A. alnifolia, S. canadensis, S. betulifolia, and V. membranaceum) after harvest alone (South 11) compared to 8 to 13 species on the burned plots. The H and J functions both increased with the addition of fire, with maximum H occurring in plots with the most severe disturbance (i.e., 100% loss of forest floor by wildfire in South 12/13) (Figure 3).  x-axis values progress from undisturbed forest (left) to the highest severity disturbance (right); see Table 1 for disturbance category descriptions.  Table 1 for disturbance category descriptions.
Eight species were present as trees ( Table 3). The two Populus species, tremuloides and trichocarpa, were each observed in four plots, while T. brevifolia was observed in just one. The other five species (A. lasiocarpa, L. occidentalis, P. engelmannii, Pinus contorta, P. menziesii) were seen in nearly all plots, had relative abundance values of 13% to 25%, and accounted for 94% of all trees and 99% of the tree biomass observed (Table S3).

MEFE
Increasing levels of fire severity promoted the occurrence of more tree seedlings (Table 3 and  Table S2). On the harvest-only East 2 plot, A. lasiocarpa, P. engelmannii, and P. menziesii dominated in terms of relative abundance (Figure 7 and Table S2) and stems ha −1 (Table 3), but A. lasiocarpa and P. menziesii accounted for almost all seedling biomass (Figure 7 and Table S2). The addition of a low-severity prescribed fire increased the relative abundance and biomass of P. engelmannii and decreased that of P. menziesii in East 3 (harvest and prescribed fire, 58% loss of forest floor), whereas the most severe fire (North 12; harvest and prescribed fire, 91% loss) resulted in a plot dominated by A. lasiocarpa (Figure 7).  Table 1 for disturbance category descriptions.  Note: x-axis values progress from undisturbed forest (left) to the highest severity disturbance (right); see Table 1 for disturbance category descriptions.
In contrast to shrub biomass, tree biomass decreased with increasing disturbance severity (Table 4). In the low-disturbance, harvest-only plot East 2, the same species that dominated as seedlings in terms of composition also dominated as trees, but most of the tree biomass was associated with A. lasiocarpa ( Figure 8 and Table S3). The increasing disturbance severity from East 2 to East 3 to North 12 had mixed effects on L. occidentalis, while its relative abundance increased from 2% to 28% and its relative biomass grew from 4% to 75% (Figure 8 and Table S3); the most trees ha −1 were observed in East 3, where the forest floor was reduced about 57% by prescribed fire ( Table 3). The fine root biomass was highest in East 3, which experienced a moderate disturbance level.  Table 1 for disturbance category descriptions.  Table 1 for disturbance category descriptions. Note: x-axis values progress from undisturbed forest (left) to the highest severity disturbance (right); see Table 1 for disturbance category descriptions.

XETE
Absence of fire eliminated the occurrence of any seedlings 30 years after harvest in South 11 ( Table 3). As the fire severity increased, the dominant seedling species in terms of relative abundance ranged from A. lasiocarpa (West 10; harvest and prescribed fire, 53% loss of forest floor) to P. contorta (South 7; harvest and wildfire, 72% loss) to P. tremuloides (South 8, harvest, prescribed fire, and wildfire, 84% loss), with L. occidentalis being most abundant in the sites with the highest fire severity (West 2, harvest and wildfire, 98% loss; South 12/13, wildfire, 100% loss) (Figure 9 and Table S3).
Regarding the aboveground relative biomass of trees, A. lasiocarpa was dominant in the harvestonly plot South 11, whereas any occurrence of fire stimulated regeneration of L. occidentalis ( Figure  10, Tables 3 and S3). As was the case for MEFE, the highest numbers of L. occidentalis seedlings and trees ha −1 were observed in the plot with the lowest fire severity (i.e., West 10, harvest and prescribed fire, 53% loss of forest floor), where the number was 3 times that of any other disturbance category ( Figure 10 and Table 3).  Table 1 for disturbance category descriptions.  Table 1 for disturbance category descriptions.
Regarding the aboveground relative biomass of trees, A. lasiocarpa was dominant in the harvest-only plot South 11, whereas any occurrence of fire stimulated regeneration of L. occidentalis ( Figure 10, Table 3 and Table S3). As was the case for MEFE, the highest numbers of L. occidentalis seedlings and trees ha −1 were observed in the plot with the lowest fire severity (i.e., West 10, harvest and prescribed fire, 53% loss of forest floor), where the number was 3 times that of any other disturbance category ( Figure 10 and Table 3).  Table 1 for disturbance category descriptions.

Total Vegetation
The three disturbed MEFE plots averaged just 20% of the total vegetation biomass found in the undisturbed, 225-year-old stands (UM; Table 4). In contrast, the six disturbed XETE plots averaged about 34% of the total vegetation biomass of the average found for the undisturbed plots (UX). The biomass in the least disturbed harvest-only South 11 plot was just 19% of the UX value, whereas the two plots with the most disturbance (West 2, harvest and wildfire; South 12/13, wildfire) had almost 50% of the biomass value, 28% of the total vegetation C pool value, and 43% of the total vegetation N pool value found for UX (Table 4).
Although we segregated aboveground plant tissues by species and into four categories (wood or foliage from plants with basal diameters < 2 cm or > 2 cm), we found that overall basal diameter was not significant when sampling foliage or wood (branches) from tree seedlings. Basal diameter  Table 1 for disturbance category descriptions.

Total Vegetation
The three disturbed MEFE plots averaged just 20% of the total vegetation biomass found in the undisturbed, 225-year-old stands (UM; Table 4). In contrast, the six disturbed XETE plots averaged about 34% of the total vegetation biomass of the average found for the undisturbed plots (UX). The biomass in the least disturbed harvest-only South 11 plot was just 19% of the UX value, whereas the two plots with the most disturbance (West 2, harvest and wildfire; South 12/13, wildfire) had almost 50% of the biomass value, 28% of the total vegetation C pool value, and 43% of the total vegetation N pool value found for UX (Table 4).
Although we segregated aboveground plant tissues by species and into four categories (wood or foliage from plants with basal diameters < 2 cm or > 2 cm), we found that overall basal diameter was not significant when sampling foliage or wood (branches) from tree seedlings. Basal diameter was, however, significant for shrubs (excluding N-fixing species; Table 5); N concentrations for wood and foliage ranged from 0.38 to 0.80 and 1.16 to 2.45, respectively, and we noted significant differences among species (Table S5). As with tree species, basal diameter was not significant for overall N concentrations for the wood or foliage of N-fixing shrub species (p = 0.9161 and 0.6327, respectively), but we observed differences among species (Table 6). In particular, the wood of S. canadensis had a significantly higher N concentration, whereas C. velutinus had the lowest foliage N concentration.

Woody Residue and Soil
Regardless of habitat type phase and disturbance, the amount of woody residue was similar and low (1-9 Mg ha −1 ), at about 25% of the undisturbed MEFE (UM) and XETE (UX) plots (Table 4 and Table S4). Across the landscape, the forest floor levels decreased with increasing amounts of disturbance at the onset of the study. Despite this, in the moister MEFE sites the amount of forest floor exceeded the levels found in undisturbed plots (UM), whereas in XETE sites the levels were about one-third of the undisturbed plots (UX); overall, forest floor levels had recovered to about 80% of the undisturbed plots (UM and UX). Soil wood provided the smallest contribution to soil biomass (Table 4), the levels of which were about 60% of those in undisturbed plots (UM and UX). Mineral soil biomass values were about 25% higher in MEFE plots than XETE plots, but regardless of disturbance or habitat type phase, the values usually increased (up to 40%) compared to the undisturbed plots. Overall, the soil biomass levels were about the same in disturbed plots compared to undisturbed plots.

MEFE
Surface coarse woody residues and mineral soil biomass values were less affected by disturbance than the forest floor and soil wood categories were. Woody residue biomass values were similar for all disturbance treatments and were 21% of the values in undisturbed plots (UM) ( Table 4). Forest floor levels were significantly lower in the harvested North 12 plot (116 (2) Mg ha −1 ), which had 91% of its forest floor removed by prescribed fire, compared with either of the other two plots that were harvested and burned by prescription (East 2 and East 3), which averaged 185 Mg ha −1 and showed significantly higher values than undisturbed plots (139 (60) Mg ha −1 ) after 30 years (Table S4). Soil wood values ranged from 23 to 69 Mg ha −1 (Table S4). Overall, total soil biomass levels (inclusive of forest floor, soil wood, and mineral soil) in the least disturbed plots (East 2 and East 3, harvested and prescribed fire with > 42% of the forest floor remaining) exceeded those measured on the undisturbed plots (UM). North 12, which experienced the most disturbance (harvest and prescribed fire, 91% loss of forest floor), had 88% of the total profile biomass of UM (Table 4 and Table S4).

XETE
Woody residue biomass in South 12/13, the uncut plot consumed by wildfire, was 2-to 6-fold higher than the other disturbed plots, but the level was still lower than in the undisturbed plots (UX) ( Table 4). Despite soil wood in West 10 (harvest and prescribed fire) being 3 times higher than in UX, and despite some XETE plots having no soil wood, we found no significant differences among plots (Table S4). Unlike in the MEFE plots, the mineral soil biomass values were similar among the disturbed and undisturbed plots (Table 4).

C and N Pools
The total vegetation biomass across the disturbed plots was about 25% of that for the undisturbed plots, whereas the soil biomass level was similar to that for undisturbed plots ( Table 4). The overall biomass (vegetation + soil) for disturbed plots was about 55% of that for the undisturbed plots. Carbon generally followed this same pattern. The total N in the vegetation was low, at about 37% of that in undisturbed forest plots, regardless of the habitat type phase. The total soil N exceeded the values in the undisturbed plots by 22% (MEFE) and 8-fold (XETE).

MEFE
Aboveground biomass (plus roots) had 83% less C and 32% less N on average in the disturbed plots than in undisturbed plots (Table 4). For total soil, disturbed plots averaged nearly the same C pool as undisturbed plots (UM). Soil samples in plots that experienced harvest and prescribed fire where > 42% of the forest floor remained (East 2 and East 3) had about 120% of the total N pool values as those found in UM plots. Soil total C and N pools in the most disturbed plot (North 12; harvested and prescribed fire, 91% loss of forest floor), however, were only about two-thirds those of UM after 30 years, with less than half of the N as in the other two MEFE plots (Table 4). For total pools (vegetation plus soil), the two plots with the lowest relative disturbance levels (East 2 and East 3) averaged about 58% and 111% of the C and N pools, respectively, of the undisturbed plots (UM), whereas when most (91%) of the forest floor was removed by prescribed fire (North 12), the total C and N pools were slower to recover and reached only 35% and 60% of the levels of the undisturbed plots ( Table 4). As expected, the greatest coarse woody residue C and N pools were located in undisturbed plots (UM). In all MEFE phase plots, the dominant C and N pools were found in the soil.

XETE
After 30 years, total C and N pools (vegetation plus soil) of the disturbed plots were about 44% and 69%, respectively, of the undisturbed plots (UX) ( Table 4). Total soil N increased with the increasing relative biomass of N-fixing plants ( Figure 11). Overall, the mineral and organic portions of the soil in these six disturbed plots averaged about 70% of the total soil C pool and about 80% of the total soil N pool values as those found in UX (Table 4). Compared to UX, the total C pool in the most disturbed plot (South S12/13, 100% of forest floor consumed by wildfire) was about 59% of the UX value, whereas for the plots with the least disturbance (South 11, harvest-only; West 10, harvest and burned by prescription; >53% of forest floor remained) the values were about 84% of the UX value.
In addition, the total N in the vegetation plus soil in the UX was 32% greater than the disturbed plots, with most located in the soil pool. The total C and N pools were generally lower in XETE plots as compared to MEFE plots.
Forests 2020, 11, x FOR PEER REVIEW 25 of 37 most disturbed plot (South S12/13, 100% of forest floor consumed by wildfire) was about 59% of the UX value, whereas for the plots with the least disturbance (South 11, harvest-only; West 10, harvest and burned by prescription; >53% of forest floor remained) the values were about 84% of the UX value. In addition, the total N in the vegetation plus soil in the UX was 32% greater than the disturbed plots, with most located in the soil pool. The total C and N pools were generally lower in XETE plots as compared to MEFE plots.
Soil nitrogen pool (kg ha -1 )  Figure 11. The relationship between the relative biomass of nitrogen-fixing shrubs (Alnus viridis ssp. sinuata, Ceanothus velutinus, and Shepherdia canadensis) and the soil nitrogen pool for 30-year-old disturbances in the Xerophyllum tenax (XETE) habitat type phase within the Abies lasiocarpa/Clintonia uniflora habitat type. See Table 1 for disturbance category descriptions.

Shrub Regeneration
In forest sites, shrubs have a key role in supporting forest function by providing wildlife habitats, cycling nutrients, and adding to overall forest biodiversity [78][79][80]. After disturbances that remove the tree overstory, the shrub cover can expand rapidly due to increased access to light, water, and nutrients, and because shrubs regenerate promptly by sprouting from the root collar, from rhizomes, or from the germination of seeds from the seed bank [81][82][83][84]. After disturbances in the northern Rocky Mountains, shrubs reach their maximum potential in about 10-30 years [85][86][87]. In a nearby study site in northwestern Montana with similar species to those in our study, the shrub community after biomass removal was persisting at pre-harvest levels 38 years later [88]. Our measurements from five out of six disturbed XETE plots 30 years after disturbance revealed that the relative abundance of dead shrubs was high (>18%), and that in plots with the greatest disturbance levels (i.e., West 2, harvest and wildfire; South 12/13, wildfire; >98% loss of forest floor) the relative abundance of dead shrubs surpassed the relative abundance of any individual live species. This decline of shrubs, which can occur rapidly [87], is likely the result of canopy closure [85,89], marking entry into the stem exclusion phase of forest stand development [90]. We did not, however, see a similar amount of mortality in the more mesic MEFE sites, which also had comparatively higher biomass production. These results contrast with Alldredge et al. [87]; their productivity models, which were developed using data obtained from a range of silvicultural methods, suggest higher rates of productivity for southern aspects than for northern aspects and a faster decline of shrubs on northern aspects.
In the undisturbed MEFE plots (UM), the overall number of shrub species was high (11); the most prevalent shrub species were M. ferruginea and V. membranaceum. While both species can Figure 11. The relationship between the relative biomass of nitrogen-fixing shrubs (Alnus viridis ssp. sinuata, Ceanothus velutinus, and Shepherdia canadensis) and the soil nitrogen pool for 30-year-old disturbances in the Xerophyllum tenax (XETE) habitat type phase within the Abies lasiocarpa/Clintonia uniflora habitat type. See Table 1 for disturbance category descriptions.

Shrub Regeneration
In forest sites, shrubs have a key role in supporting forest function by providing wildlife habitats, cycling nutrients, and adding to overall forest biodiversity [78][79][80]. After disturbances that remove the tree overstory, the shrub cover can expand rapidly due to increased access to light, water, and nutrients, and because shrubs regenerate promptly by sprouting from the root collar, from rhizomes, or from the germination of seeds from the seed bank [81][82][83][84]. After disturbances in the northern Rocky Mountains, shrubs reach their maximum potential in about 10-30 years [85][86][87]. In a nearby study site in northwestern Montana with similar species to those in our study, the shrub community after biomass removal was persisting at pre-harvest levels 38 years later [88]. Our measurements from five out of six disturbed XETE plots 30 years after disturbance revealed that the relative abundance of dead shrubs was high (>18%), and that in plots with the greatest disturbance levels (i.e., West 2, harvest and wildfire; South 12/13, wildfire; >98% loss of forest floor) the relative abundance of dead shrubs surpassed the relative abundance of any individual live species. This decline of shrubs, which can occur rapidly [87], is likely the result of canopy closure [85,89], marking entry into the stem exclusion phase of forest stand development [90]. We did not, however, see a similar amount of mortality in the more mesic MEFE sites, which also had comparatively higher biomass production. These results contrast with Alldredge et al. [87]; their productivity models, which were developed using data obtained from a range of silvicultural methods, suggest higher rates of productivity for southern aspects than for northern aspects and a faster decline of shrubs on northern aspects.
In the undisturbed MEFE plots (UM), the overall number of shrub species was high (11); the most prevalent shrub species were M. ferruginea and V. membranaceum. While both species can regenerate from root crowns, M. ferruginea also reproduces from seeds [91], whereas V. membranaceum also reproduces from rhizomes [92]. Although harvesting had little influence on their relative abundance (East 2) or overall species biodiversity, their relative abundance and relative biomass levels were reduced by the increasing levels of fire severity that removed the forest floor, likely destroying the seed bank and rhizomes. In contrast, A. v. sinuata, an actinorhizal N-fixing species found in moist sites, increased in relative abundance as fire severity increased. This species is an aggressive colonizer after fire because of its ability to sprout from the root crown and its prolific production of wind-borne seeds [93]. Reproduction through wind-borne seeds is characteristic of large, woody, actinorhizal N-fixing plants, which leads to them being more evenly distributed across the landscape and favoring earlier development in succession [94].
In the drier XETE sites, any level of disturbance reduced V. membranaceum, the most abundant shrub in the undisturbed plots. Removal of the overstory in South 11 and harvest and fire in the other XETE plots nearly eliminated this shade-tolerant, rhizomatous, late seral species. Disturbance did, however, result in four species (S. scouleriana, A. v. sinuata, S. canadensis, and C. velutinus) that were not observed in undisturbed plots (UX) becoming abundant, three of which are actinorhizal N-fixers. We did not observe S. scouleriana in any of the undisturbed plots (UX or UM), in agreement with other observations [87]. Thus, even though S. scouleriana can sprout from the root crown [95], its high relative abundance in the disturbed XETE plots suggests that its small, fluffy, airborne seeds are its prominent method of early successional colonization [81]. As in some MEFE plots, A. v. sinuata was a robust colonizer, but only on the western aspect; the southern aspect was likely too dry [93]. Shepherdia canadensis can regenerate by sprouting from the root crown and via its seeds [96]. Its high relative abundance in South 11, where harvest was the sole disturbance, could be a function of existing plants able to grow more vigorously without the overstory competition for resources [84], likely because this species is more prevalent in later successional stages, and as we observed can be present in old-growth stands [97]. Its presence may also be a result of the off-site transfer of seeds onto the disturbed forest floor [96]. Because seed scarification is often used to stimulate germination [98], it is likely that seeds persist in the seed bank. However, the abundance of S. canadensis in sites with frequent fire (<50 year intervals [99]) suggests a shorter longevity in the seed bank than C. velutinus, for which some evidence suggests can persist for 200 + years [100].
Ceanothus velutinus presents its greatest potential on southern aspects [101], where in the first three years after severe fire it can account for 40% of the relative abundance [102], and within seven years can bring site N concentrations (i.e., soil, litter, and biomass) to pre-burn levels [103]. While these four species were absent in the undisturbed XETE plots (UX), they became prevalent post-disturbance, leading to an overall increase in shrub diversity compared to the UX plots. In the harvest-only plot South 11, the resurgence of S. canadensis with 71% relative abundance yielded the lowest within-plot species diversity (H ) and evenness (J ) in the study. It is noteworthy that any disturbance on the southern aspects resulted in robust development of an N-fixing shrub species, and unlike a previous study [104], we observed mixtures of N-fixing shrub species in some plots.
We observed a strong relationship between the relative abundance of the actinorhizal N-fixing shrubs (A. v. sinuata, C. velutinus, and S. canadensis) and the amount of N in the total soil pool. While Newland and Deluca [104] found that potentially mineralizable N (PMN) and the ratio of PMN to total N were greater in non-burned stands than burned stands when sampling herbaceous and woody N-fixing species in western Montana, they found no correlation between PMN and N-fixer abundance (i.e., cover or frequency). In developing floodplain soils in Alaska, S. canadensis contributed about 50% of the total N accretion obtained from inception through to 300 years later, with maximum accretion occurring when this shrub dominated the community [105]. Reported N-fixation rates for S. canadensis and C. velutinus vary widely, likely the result of differing shrub densities, community dynamics, measuring techniques, and other variables [106][107][108][109]. Nonetheless, recent studies have demonstrated that where productivity is limited by N availability, non-N-fixing plants respond positively to increased soil N generated by N-fixing plants [110,111], although changes to soil biota [111], increased soil organic matter levels [110], subsequent slow release of chemically bound organic N [112], and service as nurse plants for other species' establishment and growth [97,113] may also be important.

Larix occidentalis Regeneration
Natural regeneration of L. occidentalis is negatively correlated with the amount and diversity of understory species [114]. A primary objective of the original 1967 treatments was to elucidate silvicultural practices to encourage natural regeneration of L. occidentalis, with the general finding that although germination occurred in a variety of seedbed types, seedling survival was limited to mineral soil seedbeds [18,115]. This is not surprising, given its high level of shade and drought intolerance [116]. On L. occidentalis sites, fire can reduce forest floor and microsite heterogeneity (i.e., increases levels of bare mineral soil) and competing vegetation [117], which would reduce competition for light and water. In East 2 and South 11, where the only disturbance was harvest, little regeneration of L. occidentalis was apparent, likely because of the lack of mineral soil substrate and the high abundance of shrubs. The addition of any level of fire stimulated the regeneration of L. occidentalis. On the harvested East 3 plot, which experienced low-severity fire, L. occidentalis regenerated. The relative abundance of biomass for L. occidentalis compared to other tree species was similar to the levels found in undisturbed MEFE plots. The likely rapid recovery of other vegetation in East 3 probably decreased the number of suitable microsites for L. occidentalis seed germination and persistence, as revealed by a subsequent lack of seedlings 30 years later. Harvested plot North 12 experienced a hotter fire, as evidenced by the greater loss of forest floor; although colonization by the high-stature shrubs A. v. sinuata and S. scouleriana greatly reduced the number of L. occidentalis trees ha −1 , the rapid juvenile growth of this species [118] did allow some trees to out compete the shrubs, leading to its high relative biomass.
In the hotter, drier XETE plots, where fire-induced changes to the forest floor may have persisted longer, more L. occidentalis seedlings were observed compared to MEFE plots, especially in South 12/13, which experienced a complete loss of forest floor in the 1967 stand-replacing wildfire. While all levels of disturbance in XETE plots fostered L. occidentalis regeneration, we noted a trend of decreasing L. occidentalis relative abundance and relative biomass with increasing fire severity, however the actual numbers of trees ha −1 were quite variable. Undoubtedly other factors may be involved, including disturbance plot hysteresis. For example, in South 13, P. contorta was a component of the original stand (7%) [51], and its serotinous cones likely contributed to the dominance of this species in this disturbance category.
The ability of L. occidentalis to readily regenerate when presented with early seral stand conditions as the result of fire has allowed this species to respond favorably to short-term changes in climate [119]. The continuing shift toward a hotter and drier climate, which facilitates more frequent and intense fires, has benefitted seedling recruitment, indicating that climatic factors are apparently less important to regeneration than those associated with overstory competition [119].

Soil Responses and C and N Pools
Low-temperature surface fires alter forest ecosystems; surface combustion reduces the levels of WR, forest floor, soil wood, and understory vegetation, and thereby releases nutrients, which in turn increase nutrient availability. As fires increase in severity, the result is more oxidation of OM within the mineral soil, volatilization of nutrients, and reduced nutrient availability, which may persist for a few years or have a longer term impact [120]. This is evident in our moister MEFE plots; the harvested East 3 plot, which experienced lower severity prescribed fire, had a soil N pool similar to the undisturbed plots (UM); whereas harvested North 12, which experienced higher severity prescribed fire, did not. Similarly, while the overall vegetation C pool in disturbed MEFE plots averaged only 17% of that in UM, the soil C pools were equal to or exceeded those of UM plots, except under the highest disturbance level in North 12. In the drier XETE plots, the soil C values were greater in the undisturbed plots, but N values were generally lower.
In the MEFE plots, the proportion of biomass in the WR, surface organic horizons, and soil wood were greater in disturbed plots than in the undisturbed plots (UM), but the distribution of these substrates was different. For example, WR levels were almost 5 times greater in UM plots and 3 times greater in undisturbed XETE plots than in the disturbed plots. Across both MEFE and XETE plots, the total WR biomass was less than half that of undisturbed plots, with the greatest residue levels occurring in the uncut plot burned by wildfire (South 12/13). Woody residue amounts are a function of forest type, site type, time since disturbance, type of disturbance, and decomposition [121][122][123]; therefore, in the same climate, plots with the greatest severity of disturbance, whether from harvest or fire, will have less WR. The woody residue, forest floor, and soil wood account for 44-84% of C in several mid-successional stands across the western USA [29]. Similarly, the C pools in disturbed plots represented 88% and 48% of the C pools in MEFE and XETE plots, respectively, as compared to the undisturbed plots for either habitat type phase. The forest floor and soil wood levels in the moister MEFE plots had nearly recovered from burning.
The recovery of the forest floor and soil wood is important for overall site productivity [35], because of the physical (i.e., water-holding capacity), chemical (i.e., nutrient cycling), and biological (i.e., non-symbiotic N-fixation [124]) functions they support. Therefore, even though we saw reductions in WR, forest floor, and soil wood levels relative to undisturbed MEFE plots, the C and N pools we observed 30 years after harvesting and prescribed fire in the MEFE plots suggest that these activities were not detrimental to stand productivity. Our results are similar to those in a moist, east-facing L. occidentalis forest in western Montana [21]. After three decades, however, the XETE disturbed plots showed less recovery than the moister MEFE plots. In particular, most of the soil wood was removed in all XETE plots (except West 10) after harvesting and fire. Soil wood and forest floor material are critical for supporting ectomycorrhizal activity, maintaining soil moisture through the summer drought, and facilitating seedling survival and understory development [125].
In the western USA, N is a limiting nutrient in soils [126], with N pools typically being much larger in the mineral soil than surface organic layers [29]. We did not, however, observe this pattern in our MEFE or XETE plots, where more N was located in the forest floor than in the mineral soil samples to a 30 cm depth. This means that soil N pools are at-risk if another fire occurs, the result of which could compromise stand recovery.
It is important to recognize how stands accumulate C and N after disturbance, as this has implications for climate change. Increased fire frequency in Canada has converted boreal forests to a net source for C [127]. Unless high-severity fires predominate, woody residues are resistant to burning, meaning this pool would be a resilient source of ecosystem C. Further, in wetter ecosystems, such as our MEFE plots, the forest floor would likely recover before the next disturbance. In the drier XETE habitat type phase, the C and N soil pools in harvest-only South 11 were close to initial levels within the forest floor, however if wildfire occurs this pool will be vulnerable.

Management Implications
When this experiment was initiated in the late 1960s, the management focus for federal forests in the USA was extraction of timber, and L. occidentalis, with its desirable growth characteristics, was seen as an important species [128]. As noted, however [129], traditional forest management techniques often target a specific outcome that can lead to unforeseen negative consequences. Often, the result is more homogeneity among stands across the landscape [129], associated with loss of biodiversity, and subsequently resilience, which is necessary to retain as a vital pillar of climate change adaptation in forests [130]. Despite the use in western Montana of prescribed fire following clearcut harvesting, which are silvicultural techniques that emulate stand-replacing wildfire, researchers suggest that silvicultural practices may have caused a fundamental loss of N-fixing species across the landscape, with possible repercussions in terms of long-term nutrient cycling and forest productivity [104].
As concluded in [129], silvicultural techniques that retain patterns and processes of natural disturbances can be successful, however the fundamental differences between practiced silviculture techniques and natural disturbances mean that full ecosystem benefits may not be provided. Of course, the natural disturbance being emulated and the realization of benefits are both social constructs often lacking consensus [8,131].
Nonetheless, the novelty of our experiment was the opportunity to compare a range of silviculture techniques (e.g., harvest and prescribed fire) in combination with a naturally occurring, stand-replacing wildfire. Our results suggest that harvest plus prescribed fire is a viable method for maintaining L. occidentalis as part of a mixed conifer forest, and can also support diverse understory shrub species, especially N-fixing species. While we observed composition shifts with varied disturbance levels, the homogenous overall increase in shrub biodiversity, abundance, and biomass after fire suggests attainment of a heterogenetic landscape [83]. In this experiment, many of the soil metrics we measured were approaching, or had already surpassed, those of the undisturbed forest. In contrast, the vegetation total biomass and C levels were a fraction of those observed in the undisturbed forest. Given the wood volume of the undisturbed forest this is not surprising, and gives support to others who concluded that achieving pre-disturbance C pools requires mature trees [20,24,26,132].
Although the recovery trajectory of disturbed forests can vary in the northern Rocky Mountains, "the eventual outcome is relatively consistent", which is a "dominant conifer overstory with an understory of shade-tolerant shrubs and herbs" [81] (p. 357). Our results appear to concur with that conclusion, and confirm that L. occidentalis regeneration can be supported by harvest and prescribed fire in mixed conifer forests, particularly on southern and western aspects. As the climate niche for L. occidentalis continues to retract within its current range and regeneration shifts to the cool, dry subset of the sites now occupied by this species [119], land managers may need to place additional focus on eastern and northern aspects. Our results show that L. occidentalis can be regenerated with harvest and prescribed fire in these cooler, moister sites as well, although additional treatment(s) may be necessary to control competition that would otherwise prevent this seral, shade-intolerant species from establishing its renowned place in western USA landscapes and from fulfilling desired ecosystem functions.

Conclusions
Our study evaluated the 30-year responses of vegetation and soil to 11 disturbances of varying severity, which occurred within the context of a large-area (240 ha) study initiated in 1967. Although the original objective of this management-driven study was to evaluate the effects of prescribed burning at various times of the year on the establishment and growth of L. occidentalis following clearcut harvesting, a 1967 wildfire within the study site allowed the effects of wildfire to be incorporated. Thus, we leveraged this experiment to measure changes in plant communities (species composition and biomass) and soils (C and N pools) in different habitat type phases in combination with harvesting and fire.
Our first hypothesis and its corollaries were that increasing levels of disturbance severity would increase species richness, diversity, and evenness more in mesic sites than in xeric sites, as identified by habitat type phase, and that disturbance severity would influence the presence and abundance of N-fixing plant species and their influence on soil N pools. Our results provided mixed support. While disturbance severity in the cooler, moister habitat type phase had no effect on species diversity, it did increase diversity in the hotter, drier habitat type phase, whereas species evenness showed the opposite pattern. In the cooler, moister sites, extremes of disturbance severity skewed the individual abundance of shrub species, whereas abundance was more even in the hotter, drier sites. Overall, species diversity and evenness across all disturbance severities and habitat type phases were similar. Regardless of disturbance type (i.e., harvest, prescribed fire, or wildfire), the abundance of N-fixing shrubs increased and a positive effect of their relative abundance on the soil N pool was noted.
We hypothesized that L. occidentalis would regenerate prevalently in more severely disturbed sites and that these conducive site conditions would further support longer term recruitment. Indeed, regardless of the habitat type phase, harvesting without fire resulted in few L. occidentalis trees and no seedlings 30 years hence, whereas any severity of fire (58 to 100% losses of forest floor) promoted the density, relative abundance, and relative biomass of L. occidentalis. In the hotter, drier sites, seedling recruitment continued through the first 30 years, whereas seedling recruitment in the cooler, moister sites was only observed with the greatest fire severity (91% loss of the forest floor). Thus, the trend of greater L. occidentalis regeneration in cooler, drier sites in response to climate change may require land managers to consider higher severity prescribed burns on northern and eastern aspects to promote L. occidentalis regeneration.
We accept our final hypothesis and its corollaries, which stated that long-term C and N pool recovery depends on disturbance severity, but in different ways. Less severe disturbance favors aboveand belowground C pool development, whereas N pools increase with increasing disturbance severity. The recovery of soil C is linked to vegetation regeneration, and additional fires in this forest type may not allow for complete re-accumulation of C in above-and belowground pools. Nitrogen pools after harvesting resulted in greater N for all horizons, except soil wood. We note that although prescribed fire can increase aboveground N pools, high burn severity reduced it in the moister plots and increased it in the drier plots. Because conifer forest soil N limits productivity, the size of the N pool after disturbance and the speed of recovery will also affect future aboveground vegetative re-growth. Finally, we note that in many plots, mineral soil pools were near or above the undisturbed levels after 30 years, but these pools are located in forest floor horizons vulnerable to fire.