Microsite Influence on Woody Plant Regeneration in a Pinus palustris Woodland Following Catastrophic Disturbance

Information and material biological legacies that persist after catastrophic forest disturbance collectively constitute the ecological memory of the system and may strongly influence future stand development. Catastrophic disturbances often result in an influx of coarse woody debris (CWD), and this material legacy may provide beneficial microsites that affect successional and structural developmental pathways. We examined how microenvironmental characteristics influence the regeneration of woody plants in a subtropical woodland that experienced a large influx of CWD from a catastrophic wind disturbance. Specifically, we asked (1) what microenvironmental factors best explain woody plant density, richness, and height in the regeneration layer and (2) does woody plant density, richness, and height benefit from the large influx of CWD to a degree that competition dynamics and succession may be modified? Data were collected in a Pinus palustris woodland that had experienced an EF3 tornado and was subjected to a four-year prescribed fire rotation. We documented live woody plants <5 cm diameter at breast height, soil, and site characteristics and tested for differences in seedling and sapling density, species richness, and height in relation to CWD proximity. We used a random forest machine learning algorithm to examine the influence of microenvironmental conditions on the characteristics of woody plants in the regeneration layer. Woody plant density and species richness were not significantly different by proximity to CWD, but plants near CWD were slightly taller than plants away from CWD. The best predictors of woody plant density, richness, and height were abiotic site characteristics including slope gradient and azimuth, organic matter depth and weight, and soil water content. Results indicated that the regeneration of woody plants in this P. palustris woodland was not strongly influenced by the influx of CWD, but by other biological legacies such as existing root networks and soil characteristics. Our study highlights the need to consider ecological memory in forest management decision-making after catastrophic disturbance. Information and material legacies shape recovery patterns, but, depending on the system, some legacies will be more influential on successional and developmental pathways than others.


Introduction
Disturbances are important drivers of forest developmental and successional patterns. Forest canopy disturbances are regularly classified along a gradient based on their spatial extent, magnitude, and frequency. Along this gradient, catastrophic disturbances are among the most infrequent but the most severe, often removing most or all overstory vegetation [1][2][3]. In early forest disturbance ecology theory, catastrophic disturbances were considered to homogenize ecosystems and eliminate nearly all pre-disturbance biotic and abiotic stand conditions [4][5][6][7]. As our understanding of forest disturbance questions: (1) Which microenvironmental factors best explain woody plant density, richness, and height in the regeneration layer? and (2) Does woody plant density, richness, and height benefit from the large influx of CWD to a degree that competition dynamics and succession may be modified? Localized variability of site properties can influence woody plant regeneration [15,36,37]. Therefore, it was anticipated that (1) some soil characteristics, such as bulk density, soil water content, or organic matter depth, would influence woody plant regeneration patterns. CWD may alter microenvironmental conditions and, as a result, facilitate the establishment and growth of a variety of plant species in a range of forest types by providing shade, reducing temperature, and increasing humidity through decreased evapotranspiration [15,38]. Therefore, it was hypothesized that (2) altered microenvironmental conditions as a function of CWD input across the study site would influence patterns in density, richness, and height of the regeneration layer [25,39,40].

Study Site
This study was conducted in the Oakmulgee Ranger District of the Talladega National Forest in central Alabama, USA. The Oakmulgee Ranger District is in the Fall Line Hills [41]. The Fall Line Hills comprise sedimentary rock belts that serve as a transition zone for the Coastal Plain and Appalachian Highlands [42]. Soils are mainly derived from the Maubila series, which is described as deep and moderately well drained with a subangular blocky structure [43]. Specifically, data were collected over an area with two variations of the Maubila series, the Maubila flaggy loam and the Maubila-Smithdale complex [43] (Figure 1). The two soil delineations both have ca. 14 cm of available water storage, and the minimum water table depth for both is 23 cm [43].
Forests 2020, 11, x FOR PEER REVIEW 3 of 19 this study were met by answering the following questions: (1) Which microenvironmental factors best explain woody plant density, richness, and height in the regeneration layer? and (2) Does woody plant density, richness, and height benefit from the large influx of CWD to a degree that competition dynamics and succession may be modified? Localized variability of site properties can influence woody plant regeneration [15,36,37]. Therefore, it was anticipated that (1) some soil characteristics, such as bulk density, soil water content, or organic matter depth, would influence woody plant regeneration patterns. CWD may alter microenvironmental conditions and, as a result, facilitate the establishment and growth of a variety of plant species in a range of forest types by providing shade, reducing temperature, and increasing humidity through decreased evapotranspiration [15,38]. Therefore, it was hypothesized that (2) altered microenvironmental conditions as a function of CWD input across the study site would influence patterns in density, richness, and height of the regeneration layer [25,39,40].

Study Site
This study was conducted in the Oakmulgee Ranger District of the Talladega National Forest in central Alabama, USA. The Oakmulgee Ranger District is in the Fall Line Hills [41]. The Fall Line Hills comprise sedimentary rock belts that serve as a transition zone for the Coastal Plain and Appalachian Highlands [42]. Soils are mainly derived from the Maubila series, which is described as deep and moderately well drained with a subangular blocky structure [43]. Specifically, data were collected over an area with two variations of the Maubila series, the Maubila flaggy loam and the Maubila-Smithdale complex [43] (Figure 1). The two soil delineations both have ca. 14 cm of available water storage, and the minimum water table depth for both is 23 cm [43]. The regional climate is humid mesothermal, characterized by year-round rainfall; long, hot summers; and brief, mild winters [44]. The 30 year normal annual precipitation and temperature are The regional climate is humid mesothermal, characterized by year-round rainfall; long, hot summers; and brief, mild winters [44]. The 30 year normal annual precipitation and temperature are 1376 mm and 17 • C, respectively [45]. February has the highest precipitation at 139 mm, and October has the lowest at 87 mm [45]. The coldest month is January, with a mean temperature of 7 • C, and the warmest month is July, with a mean temperature of 27 • C [45]. Since 2000, severe droughts have occurred in the study area every 5-7 years, and the last severe drought was in 2017 [46].
The study was conducted in the central Pinus palustris belt of Alabama [47], and the plant communities liken to the Quercus-Pinus forest region [48]. Pinus palustris historically dominated the forest canopy of the study area; however, Pinus taeda, Pinus echinata, and various hardwood species can reach subcanopy and canopy positions [49,50]. The ecosystem supports a diverse herbaceous layer and an open mid-story [51,52]. Fire is an important disturbance in the advancement of Pinus palustris systems, maintaining the characteristically diverse understory and preventing succession to hardwood dominance [53]. Additionally, Pinus palustris woodlands require frequent low-intensity surface fires to decrease interspecific competition and expose bare mineral soil for seed germination [54].
Tornadoes occur frequently in this region [47,55]. The National Climatic Data Center has recorded 27 tornado events in the county since 1950, and three of these events were of EF3 or greater magnitude [56]. On 27 April 2011 an EF3 long-tracked wedge tornado, one of the 362 confirmed tornadoes of the 2011 super outbreak, moved through the Oakmulgee Ranger District. In fact, this was one of three tornadoes to impact the forest during the 2011 super outbreak. The EF3 tornado had maximum wind speeds of 233 kph and a maximum path width of 1609 m [57]. The tornado created a distinct edge between the undisturbed forest and the swath of catastrophically disturbed forest, and reduced basal area from 22 m 2 ha −1 in the undisturbed forest to 1 m 2 ha −1 in the tornado swath [58].
To promote Pinus palustris recovery, low-intensity prescribed fires are set every 2-4 years in P. palustris-dominated stands across the Oakmulgee Ranger District. Since the tornado in 2011, two operational, prescribed fires have been conducted, in April 2014 and April 2018. We sampled in the second growing season after the 2018 prescribed fire.

Field Methods
In the summer of 2019, 350 quadrats (300 quadrats contiguous to CWD and 50 quadrats located away from CWD) were established in stands that had experienced catastrophic wind disturbance and were not salvage-harvested on the Oakmulgee Ranger District. CWD sampling only occurred in a single compartment on the district to ensure consistency in prescribed fire, land-use, and management history. Transects, beginning from random starting points, were used to locate and sample CWD that met the following criteria: >3 m in length and ≥10 cm in diameter.
Our sampling scheme, modified from Chećko et al. [59], consisted of a group of three 1 m 2 quadrats placed equidistant on the ground on both sides of the CWD to analyze microsite characteristics. Thus, six quadrats were established around one CWD fragment with three quadrats on one side of the log and the other three quadrats on the other side of the CWD (Figure 2). The quadrats were uniformly established by placing Quadrat 1 on the downslope side of the CWD, and then moving clockwise to assign the five remaining quadrats. A seventh quadrat was established at least 2 m away from the focal log and any other CWD. Forests 2020, 11, x FOR PEER REVIEW 5 of 19 For every quadrat contiguous to CWD, multiple attributes of the focal CWD were recorded: tree species to the lowest taxon possible, diameter, length, presence of rot, decay class, and azimuth. Azimuth of CWD was always recorded from upslope to downslope. The area where data were collected was under a uniform prescribed fire plan; therefore, all of the CWD analyzed was subjected to prescribed fire. However, if the CWD had visual signs of combustion where a quadrat was established, it was scored as charred [60,61]. CWD was also categorized as being either on the ground (i.e., resting on the ground) or elevated (i.e., the quadrat was placed adjacent to a section of the log that was suspended). Decay class was determined using a seven decay class system adopted from Ulyshen et al. [62] that accounted for the relatively rapid decomposition of sapwood, but common decomposition of heartwood that occurs in pines of the southeastern USA [63].
For every 1 m 2 quadrat, including the seventh quadrat, the presence, species, and height of all live woody plants <5 cm diameter at breast height (dbh) were noted. Individual plants were classified as either seedlings (<1 m in height) or saplings (≥1 m in height, <5 cm dbh). Organic matter depth was documented to 0.25 cm using a ruler. We took four organic matter depth measurements equidistant from a transect perpendicular to the log through each quadrat (i.e., four readings quadrat −1 ). This method allowed for the calculation of an organic matter depth average for each quadrat and the understanding of how organic matter depth changed with distance from CWD. Organic matter measurements were not taken in portions of a quadrat where pieces of bark had recently sloughed off during CWD decay. The slope aspect and topographic position (linear, concave, or convex) for each quadrat were recorded. To examine how the amount of CWD influenced plant and soil traits, quadrats located alongside CWD were classified as ungrouped, grouped, or extremely grouped by adapting the methods of Pinno & Das Gupta [37] and Chećko et al. [59]. A quadrat was considered grouped if CWD or a stump was less than 1 m away from the quadrat edge, and classified as extremely grouped if the additional CWD was within the quadrat boundaries.
Soil samples were taken from the top of the soil profile using a 4 cm by 7.5 cm steel pipe with a volume of 94.25 cm 3 , following the methods of LaFevor [64]. Samples were collected at the center of each of the seven quadrats. Rocks occasionally impeded soil sample collection and forced the collection to occur closer to a quadrat edge than the true center. Samples were collected by removing all organic matter from the top of the soil, hammering down the steel pipe, removing the sample with a trowel, and then placing the sample in a labeled, double-sealed plastic bag. After collection, the samples were transferred to a separate bag for protection and moved to the laboratory to be weighed, dried at 105 °C for 24 h, and then reweighed [65,66]. The samples were then dry-sieved (2 mm) to remove all organic matter and rocks, which were subsequently weighed. For every quadrat contiguous to CWD, multiple attributes of the focal CWD were recorded: tree species to the lowest taxon possible, diameter, length, presence of rot, decay class, and azimuth. Azimuth of CWD was always recorded from upslope to downslope. The area where data were collected was under a uniform prescribed fire plan; therefore, all of the CWD analyzed was subjected to prescribed fire. However, if the CWD had visual signs of combustion where a quadrat was established, it was scored as charred [60,61]. CWD was also categorized as being either on the ground (i.e., resting on the ground) or elevated (i.e., the quadrat was placed adjacent to a section of the log that was suspended). Decay class was determined using a seven decay class system adopted from Ulyshen et al. [62] that accounted for the relatively rapid decomposition of sapwood, but common decomposition of heartwood that occurs in pines of the southeastern USA [63].
For every 1 m 2 quadrat, including the seventh quadrat, the presence, species, and height of all live woody plants <5 cm diameter at breast height (dbh) were noted. Individual plants were classified as either seedlings (<1 m in height) or saplings (≥1 m in height, <5 cm dbh). Organic matter depth was documented to 0.25 cm using a ruler. We took four organic matter depth measurements equidistant from a transect perpendicular to the log through each quadrat (i.e., four readings quadrat −1 ). This method allowed for the calculation of an organic matter depth average for each quadrat and the understanding of how organic matter depth changed with distance from CWD. Organic matter measurements were not taken in portions of a quadrat where pieces of bark had recently sloughed off during CWD decay. The slope aspect and topographic position (linear, concave, or convex) for each quadrat were recorded. To examine how the amount of CWD influenced plant and soil traits, quadrats located alongside CWD were classified as ungrouped, grouped, or extremely grouped by adapting the methods of Pinno & Das Gupta [37] and Chećko et al. [59]. A quadrat was considered grouped if CWD or a stump was less than 1 m away from the quadrat edge, and classified as extremely grouped if the additional CWD was within the quadrat boundaries.
Soil samples were taken from the top of the soil profile using a 4 cm by 7.5 cm steel pipe with a volume of 94.25 cm 3 , following the methods of LaFevor [64]. Samples were collected at the center of each of the seven quadrats. Rocks occasionally impeded soil sample collection and forced the collection to occur closer to a quadrat edge than the true center. Samples were collected by removing all organic matter from the top of the soil, hammering down the steel pipe, removing the sample with a trowel, and then placing the sample in a labeled, double-sealed plastic bag. After collection, the samples were transferred to a separate bag for protection and moved to the laboratory to be weighed, dried at 105 • C for 24 h, and then reweighed [65,66]. The samples were then dry-sieved (2 mm) to remove all organic matter and rocks, which were subsequently weighed.

Laboratory and Analytical Methods
We studied a natural wind disturbance to better understand how catastrophic disturbance and a large influx of CWD might influence woody plant regeneration patterns; therefore, experimental replicability was not possible. The ecological knowledge to be gained from studying natural disturbance events often outweighs potential issues of pseudoreplication; however, we sought to mitigate possible issues of pseudoreplication through the statistical treatment of our data [67]. Seedling and sapling density and richness were determined for each quadrat. Density was calculated as the sum of woody seedlings and saplings in each quadrat, and richness was calculated as the number of woody seedling and sapling species. Organic matter depth was averaged by quadrat. Relative density was calculated for woody plant seedlings and saplings. Soil water content was calculated by subtracting the soil wet weight by the dry weight and dividing the difference by the dry weight.
For this study, significance was determined at p < 0.05. To assess significant differences in microsite characteristics, the six quadrats adjacent to a piece of CWD were classified as "near" and the seventh quadrat at least 2 m away from CWD was classified as "away." We calculated the average density, richness, and height of the six near quadrats to compare to the paired away quadrat of each log (n = 50). Analysis was further stratified by seedlings and saplings. Shapiro-Wilks tests for normality and Levene's tests for equality of variance were used to ensure assumptions of parametric tests were met. Data were not normally distributed and variances were unequal. Therefore, a paired-samples Wilcoxon signed-rank test was used to statistically analyze the differences in woody plant density, richness, and height near CWD versus away [68].
To determine the influence of microsite characteristics on woody plants, we applied a random forest (RF) machine learning algorithm. RF is a complex, non-parametric tool derived from classification and regression trees and is capable of modeling complex interactions between variables with high accuracy [69,70]. The algorithm combines many trees, each generated by bootstrap samples, and retains samples for internal cross-validation (out of bag estimates, OOB). The final model output is the mean result of all trees. Model outputs include a variable importance measure, which ranks each predictor variable by mean increase in error (% IncMSE) when observed values of predictor variables are permuted through the model [69]. These values can be used to infer which predictor variable has the greatest influence on the independent variable within the model. RF was selected because of the large number of recorded CWD and soil characteristics (n = 18), and the hypothesized complex interactions between microsite characteristics and woody plant communities. We fit three independent RF models to predict woody plant density, richness, and height, incorporating 500 trees into each model. Density and richness were modeled for each individual quadrat (n = 300) and height was modeled for each individual woody plant (n = 2609). We did not include the 50 quadrats that were located away from CWD in the RF models because they were not associated with any CWD characteristics. For each model, we removed predictor variables that had negative % IncMSE values, as these variables had no influence on the woody plant characteristics, and thus reduced the predictive ability and increased the error of the models. To account for potential multicollinearity of predictor variables, we tested for correlation among all predictor variables and removed variables that were strongly correlated to one or more predictor variables. All models were cross-validated with a randomly selected training and testing subset of the data. RF was conducted using the package "randomForest" in R version 1.2.5001 [71].

Relative Density of Woody Plant Species
Vaccinium arboreum and Rhus copallinum were the most common woody plant seedling species (Table 1). Rhus copallinum, Vaccinium arboreum, and Pinus palustris seedlings exhibited the greatest relative frequencies, with Rhus copallinum seedlings appearing in 54% of quadrats, Vaccinium arboreum seedlings in 38%, and Pinus palustris seedlings in 33% of quadrats. A total of 15 woody plant species were reported only as seedlings, including Pinus palustris and five Quercus species. The three woody plant species with the highest densities in the sapling stage were Rhus copallinum, Oxydendrum arboreum, and Vaccinium arboreum ( Table 2). Rhus copallinum saplings appeared in 25% of quadrats and Vaccinium arboreum saplings were found in 14% of quadrats. About 36% of woody plant species in the sapling stage were found in only one quadrat, including Asimina parviflora, Hamamelis virginiana, Liquidambar styraciflua, Pinus taeda, Quercus marilandica, Quercus stellata, Quercus velutina, and Rhus glabra.

CWD Proximity and Woody Plant Density
Seedling density was not significantly different by proximity to CWD (p = 0.23, Figure 3a). The mean density of seedlings near CWD was 7.98 stems m −2 ± 0.59 (SE), and the mean density of woody plant seedlings away from CWD was 10.10 stems m −2 ± 1.38. The greatest seedling density was found in a plot located away from CWD, where there were 54 Vaccinium arboreum seedlings. There were 21 quadrats that had a seedling density ≥ 20 stems m −2 , and there were 317 quadrats with a seedling density <20 stems m −2 . There were 29 quadrats that were reported to have only one seedling, which included quadrats that were both near CWD and away from CWD. Sapling density was not significantly different by proximity to CWD (p = 0.17, Figure 3b). The mean density of saplings near CWD was 0.74 stems m −2 ± 0.14, and the mean density of saplings away from CWD was 0.5 stems m −2 ± 0.13. The greatest sapling density recorded in a quadrat was attributed to 19 Oxydendrum arboreum stems m −2 , and the quadrat with the second highest density had 10 woody plant stems m −2 . All other quadrats had woody plant sapling densities that were less than 10 stems m −2 . Almost 45% of quadrats had a sapling density of 1 stem m −2 .
(a) (b) Figure 3. Mean density of (a) seedlings (live woody stems <1 m height) and (b) saplings (live woody stems ≥1 m height, <5 cm dbh) of quadrats located near and away from coarse woody debris (CWD) in a Pinus palustris woodland following catastrophic disturbance. No significant differences (p < 0.05) were found.

Proximity to CWD and Woody Plant Richness
Seedling richness did not differ between quadrats near CWD and quadrats located away from CWD (p = 0.78, Figure 4a). The mean richness near CWD was 0.30 species m −2 ± 0.05 (SE), and the mean richness away from CWD was 0.34 species m −2 ± 0.07. One quadrat had a richness of 7 species m −2 , and it was near CWD. Every other quadrat had a richness value of <7. Quadrats with only one species present comprised 25% of all quadrats with woody plant seedlings. Woody plant sapling richness was also not significantly different between quadrats near CWD and quadrats away (p = 0.80, Figure 4b). Sapling richness near CWD was 0.31 species m −2 ± 0.04 and sapling richness away from CWD was 0.32 species m −2 ± 0.08. Sapling richness was documented to be only either 2 species m −2 or 1 species m −2 . Quadrats with a sapling richness of 2 species m −2 comprised 13% of quadrats, and quadrats with a richness of 1 species m −2 comprised 87% of quadrats with saplings present.  . Mean species richness of (a) seedlings (live woody stems <1 m height) and (b) saplings (live woody stems ≥1 m height, <5 cm dbh) by quadrats located near and away from coarse woody debris (CWD) in a Pinus palustris woodland following catastrophic disturbance. No significant differences (p < 0.05) were found.

Proximity to CWD and Woody Plant Height
Seedling height was significantly greater in quadrats that were near CWD than in quadrats that were away from CWD (p < 0.001, Figure 5a). The mean height for seedlings near CWD was 33.57 cm ± 0.51 (SE), and seedlings away from CWD had a mean height of 25.20 cm ± 0.93. Seedling heights ranged from 1 cm to 99 cm. Similarly, sapling height was significantly greater near CWD than away (p < 0.001, Figure 5b). Saplings near CWD had a mean height of 131.77 cm ± 2.13 and saplings away from CWD had a mean height of 120.92 cm ± 4.65.

Proximity to CWD and Woody Plant Richness
Seedling richness did not differ between quadrats near CWD and quadrats located away from CWD (p = 0.78, Figure 4a). The mean richness near CWD was 0.30 species m −2 ± 0.05 (SE), and the mean richness away from CWD was 0.34 species m −2 ± 0.07. One quadrat had a richness of 7 species m −2 , and it was near CWD. Every other quadrat had a richness value of <7. Quadrats with only one species present comprised 25% of all quadrats with woody plant seedlings. Woody plant sapling richness was also not significantly different between quadrats near CWD and quadrats away (p = 0.80, Figure 4b). Sapling richness near CWD was 0.31 species m −2 ± 0.04 and sapling richness away from CWD was 0.32 species m −2 ± 0.08. Sapling richness was documented to be only either 2 species m −2 or 1 species m −2 . Quadrats with a sapling richness of 2 species m −2 comprised 13% of quadrats, and quadrats with a richness of 1 species m −2 comprised 87% of quadrats with saplings present.

Proximity to CWD and Woody Plant Richness
Seedling richness did not differ between quadrats near CWD and quadrats located away from CWD (p = 0.78, Figure 4a). The mean richness near CWD was 0.30 species m −2 ± 0.05 (SE), and the mean richness away from CWD was 0.34 species m −2 ± 0.07. One quadrat had a richness of 7 species m −2 , and it was near CWD. Every other quadrat had a richness value of <7. Quadrats with only one species present comprised 25% of all quadrats with woody plant seedlings. Woody plant sapling richness was also not significantly different between quadrats near CWD and quadrats away (p = 0.80, Figure 4b). Sapling richness near CWD was 0.31 species m −2 ± 0.04 and sapling richness away from CWD was 0.32 species m −2 ± 0.08. Sapling richness was documented to be only either 2 species m −2 or 1 species m −2 . Quadrats with a sapling richness of 2 species m −2 comprised 13% of quadrats, and quadrats with a richness of 1 species m −2 comprised 87% of quadrats with saplings present.  . Mean species richness of (a) seedlings (live woody stems <1 m height) and (b) saplings (live woody stems ≥1 m height, <5 cm dbh) by quadrats located near and away from coarse woody debris (CWD) in a Pinus palustris woodland following catastrophic disturbance. No significant differences (p < 0.05) were found.

Proximity to CWD and Woody Plant Height
Seedling height was significantly greater in quadrats that were near CWD than in quadrats that were away from CWD (p < 0.001, Figure 5a). The mean height for seedlings near CWD was 33.57 cm ± 0.51 (SE), and seedlings away from CWD had a mean height of 25.20 cm ± 0.93. Seedling heights ranged from 1 cm to 99 cm. Similarly, sapling height was significantly greater near CWD than away (p < 0.001, Figure 5b). Saplings near CWD had a mean height of 131.77 cm ± 2.13 and saplings away from CWD had a mean height of 120.92 cm ± 4.65.

Proximity to CWD and Woody Plant Height
Seedling height was significantly greater in quadrats that were near CWD than in quadrats that were away from CWD (p < 0.001, Figure 5a). The mean height for seedlings near CWD was 33.57 cm ± 0.51 (SE), and seedlings away from CWD had a mean height of 25.20 cm ± 0.93. Seedling heights ranged from 1 cm to 99 cm. Similarly, sapling height was significantly greater near CWD than away (p < 0.001, Figure 5b). Saplings near CWD had a mean height of 131.77 cm ± 2.13 and saplings away from CWD had a mean height of 120.92 cm ± 4.65.

Influence of Microenvironmental Conditions on Woody Plant Characteristics
The RF models explained 17.2% of the variance in woody plant density, 18.4% of the variance in woody plant richness, and 21.7% of the variance in woody plant height. The woody plant density model included 16 of the 18 predictor variables. Topographic position and CWD species had negative %IncMSE values, and were thus removed from the final density model. Slope, decay class, and slope azimuth had the greatest impact on woody plant density ( Figure 6). Of the soil characteristics, wet bulk density and soil water content were the most important predictors of density ( Figure 6). The woody plant richness model included 12 of the 18 predictor variables. Topographic position, CWD species, presence of char, presence of rot, the quadrat slope position relative to CWD, and proximity to additional CWD were excluded from the final model. Wet bulk density, organic matter weight, and CWD azimuth were the most important predictors of woody plant richness (Figure 7). The woody plant height model had the highest predictive capability and explained the most variance of the three models (21.7%). Each of the 17 predictor variables were included in the height model. The three most influential variables on woody plant height were soil-related. These included organic matter depth, soil water content, and organic matter weight (Figure 8). Of the three models, the microsite characteristics that consistently had the greatest impact on woody plants (based on %IncMSE) were slope, organic matter depth and weight, and wet bulk density.

Influence of Microenvironmental Conditions on Woody Plant Characteristics
The RF models explained 17.2% of the variance in woody plant density, 18.4% of the variance in woody plant richness, and 21.7% of the variance in woody plant height. The woody plant density model included 16 of the 18 predictor variables. Topographic position and CWD species had negative %IncMSE values, and were thus removed from the final density model. Slope, decay class, and slope azimuth had the greatest impact on woody plant density ( Figure 6). Of the soil characteristics, wet bulk density and soil water content were the most important predictors of density ( Figure 6). The woody plant richness model included 12 of the 18 predictor variables. Topographic position, CWD species, presence of char, presence of rot, the quadrat slope position relative to CWD, and proximity to additional CWD were excluded from the final model. Wet bulk density, organic matter weight, and CWD azimuth were the most important predictors of woody plant richness (Figure 7). The woody plant height model had the highest predictive capability and explained the most variance of the three models (21.7%). Each of the 17 predictor variables were included in the height model. The three most influential variables on woody plant height were soil-related. These included organic matter depth, soil water content, and organic matter weight (Figure 8). Of the three models, the microsite characteristics that consistently had the greatest impact on woody plants (based on %IncMSE) were slope, organic matter depth and weight, and wet bulk density.

Influence of Microenvironmental Conditions on Woody Plant Characteristics
The RF models explained 17.2% of the variance in woody plant density, 18.4% of the variance in woody plant richness, and 21.7% of the variance in woody plant height. The woody plant density model included 16 of the 18 predictor variables. Topographic position and CWD species had negative %IncMSE values, and were thus removed from the final density model. Slope, decay class, and slope azimuth had the greatest impact on woody plant density ( Figure 6). Of the soil characteristics, wet bulk density and soil water content were the most important predictors of density ( Figure 6). The woody plant richness model included 12 of the 18 predictor variables. Topographic position, CWD species, presence of char, presence of rot, the quadrat slope position relative to CWD, and proximity to additional CWD were excluded from the final model. Wet bulk density, organic matter weight, and CWD azimuth were the most important predictors of woody plant richness (Figure 7). The woody plant height model had the highest predictive capability and explained the most variance of the three models (21.7%). Each of the 17 predictor variables were included in the height model. The three most influential variables on woody plant height were soil-related. These included organic matter depth, soil water content, and organic matter weight (Figure 8). Of the three models, the microsite characteristics that consistently had the greatest impact on woody plants (based on %IncMSE) were slope, organic matter depth and weight, and wet bulk density.

Relative Density of Woody Plant Species
Our findings indicate that Pinus palustris was primarily in the grass stage during the second growing season after the prescribed fire (2350 P. palustris seedlings ha −1 [72]), and the Pinus taeda individuals were mostly new seeds that had germinated since the April 2018 burn (the second prescribed fire since the catastrophic wind disturbance). Although not recorded in our study, the lack of Pinus palustris saplings was also attributed in part to the recruitment of individuals to the tree size class (>5 cm dbh) eight years post-disturbance. Working in the same stands, Kleinman et al. [72] noted increased post-fire Pinus palustris sapling density in catastrophically disturbed sites, but found that Pinus palustris individuals remained outnumbered by hardwood species (56 P. palustris saplings ha −1 vs. 776 hardwood saplings with canopy potential ha −1 ). However, these sites did exhibit fire-mediated sapling mortality, which indicated that hardwood and herbaceous species produced surface fuels that aided in the spread of fire in addition to the limited Pinus litter [72]. We suggest that continued use of frequent prescribed fire will allow Pinus palustris individuals to recruit into larger size classes while hardwood individuals will continue to be top-killed and thus restricted to smaller size classes. Therefore, we suggest that fire is more important than microenvironmental conditions for the

Relative Density of Woody Plant Species
Our findings indicate that Pinus palustris was primarily in the grass stage during the second growing season after the prescribed fire (2350 P. palustris seedlings ha −1 [72]), and the Pinus taeda individuals were mostly new seeds that had germinated since the April 2018 burn (the second prescribed fire since the catastrophic wind disturbance). Although not recorded in our study, the lack of Pinus palustris saplings was also attributed in part to the recruitment of individuals to the tree size class (>5 cm dbh) eight years post-disturbance. Working in the same stands, Kleinman et al. [72] noted increased post-fire Pinus palustris sapling density in catastrophically disturbed sites, but found that Pinus palustris individuals remained outnumbered by hardwood species (56 P. palustris saplings ha −1 vs. 776 hardwood saplings with canopy potential ha −1 ). However, these sites did exhibit fire-mediated sapling mortality, which indicated that hardwood and herbaceous species produced surface fuels that aided in the spread of fire in addition to the limited Pinus litter [72]. We suggest that continued use of frequent prescribed fire will allow Pinus palustris individuals to recruit into larger size classes while hardwood individuals will continue to be top-killed and thus restricted to smaller size classes. Therefore, we suggest that fire is more important than microenvironmental conditions for the

Relative Density of Woody Plant Species
Our findings indicate that Pinus palustris was primarily in the grass stage during the second growing season after the prescribed fire (2350 P. palustris seedlings ha −1 [72]), and the Pinus taeda individuals were mostly new seeds that had germinated since the April 2018 burn (the second prescribed fire since the catastrophic wind disturbance). Although not recorded in our study, the lack of Pinus palustris saplings was also attributed in part to the recruitment of individuals to the tree size class (>5 cm dbh) eight years post-disturbance. Working in the same stands, Kleinman et al. [72] noted increased post-fire Pinus palustris sapling density in catastrophically disturbed sites, but found that Pinus palustris individuals remained outnumbered by hardwood species (56 P. palustris saplings ha −1 vs. 776 hardwood saplings with canopy potential ha −1 ). However, these sites did exhibit fire-mediated sapling mortality, which indicated that hardwood and herbaceous species produced surface fuels that aided in the spread of fire in addition to the limited Pinus litter [72]. We suggest that continued use of frequent prescribed fire will allow Pinus palustris individuals to recruit into larger size classes while hardwood individuals will continue to be top-killed and thus restricted to smaller size classes. Therefore, we suggest that fire is more important than microenvironmental conditions for the recovery of these Pinus palustris stands. Many individuals on this study site were hardwood sprouts (Tables 1  and 2). Hardwood regeneration can be prolific after prescribed fire because of sprouting [73][74][75]. The high frequency of hardwood sprouts after fire is typically attributed to hardwood rootstocks that can survive fire and facilitate hardwood sprouting [76,77]. High densities of hardwoods have been documented in the second growing season after fire in mixed Quercus-Pinus systems [73,78] and indeed have been recorded in these Pinus palustris woodlands [79].
For both seedlings and saplings, Vaccinium arboreum and Rhus copallinum were among the woody plant species with the greatest density. Vaccinium arboreum seedlings can form large colonies from root sprouts [80]. Similarly, the underground rhizomes of Rhus copallinum allow it to resprout readily after fire, and past research has illustrated that the species can dominate extremely disturbed sites [81]. The prevalence of Vaccinium arboreum and Rhus copallinum on this site was also shown in past studies [82], and is likely a result of the recent catastrophic disturbance that eliminated most all overstory vegetation and subsequent prescribed fire. The high abundance of sprouters recorded in this study may also indicate that microenvironmental variability influences non-sprouters more drastically than sprouters, although this was not directly measured.

Proximity to CWD and Woody Plant Density, Richness, and Height
Seedling and sapling density were not significantly different near CWD and away. The difference in seedling density was about two quadrat −1 , and the difference in sapling density was less than one sapling quadrat −1 . Only one Pinus individual was classified as a sapling; therefore, the majority of saplings were hardwoods and likely originated from rootstocks that produced shoots that had grown to the sapling layer since the prescribed fire of the prior growing season. We hypothesized that sapling density would be greater near CWD because CWD might provide some protection to root and bud systems found belowground from prescribed fire char (60% CWD was uncharred and 40% was charred). CWD has been found to protect mycorrhizal fungi and fungal hyphae from fire, and to aid in the re-establishment of these species after fire [83]. The protection of root and bud systems could have allowed individuals to resprout more readily after fire, and protect them from future prescribed fire by helping individuals to reach greater height classes more quickly.
No significant difference was found in seedling and sapling richness by proximity to CWD. This pattern may have been because the majority of individuals in our study site originated from root stocks and not from seed. CWD can be beneficial for many small-seeded species, because the seeds can be trapped in the ridges of CWD [59,84]. However, the majority of woody plants in this study originated from an underground root system; therefore, the CWD had little influence on hardwood species density. The catastrophic wind disturbance removed nearly all mature overstory trees, and thus reduced seed availability for natural regeneration of the overstory. The residual trees that survived the tornado were primarily Pinus palustris, but residual canopy basal area was low (1 m 2 ha −1 ) [72]. Therefore, it was likely that most of the hardwood individuals recorded were biological legacies from the pre-disturbed stand, and not a function of intrastand microenvironmental variability. Species richness can potentially be increased by CWD providing more microtopographic diversity [85][86][87]. However, past research has found that CWD may not create sufficient microtopographic diversity to influence species richness patterns [88]. Other studies found that CWD was not a significant influence on vascular plant richness [89,90], and our results corroborated these findings.
Seedling and sapling height were found to be significantly greater in quadrats near CWD than in quadrats located away from CWD. However, the pattern may not be ecologically meaningful because the observed difference in seedling height was only 8 cm, and for sapling height the difference was 11 cm (i.e., seedlings near CWD were ca. 32% taller than seedlings away from CWD and saplings near CWD were ca. 9% taller than saplings away from CWD). Past CWD studies conducted in various forest types have found that CWD presence can improve individual growth rates [84,91,92]. CWD has been shown to provide a windbreak for seedlings located one to seven barrier heights on the leeward side of the log and one barrier height on the windward side, and this protection from wind can contribute to greater plant growth [93]. Additionally, a sheltering object can alter temperature and humidity values over a relatively short distance when compared to wind velocity [94]. The presence of CWD can result in a more favorable microclimate for seedling growth by reducing temperature changes and increasing humidity through decreased evapotranspiration [15]. Regardless of the height differences by CWD proximity, we posit that taller, fire-sensitive woody plants near CWD are as likely to be top-killed by a subsequent prescribed fire as woody plants away from CWD.

Influence of Microenvironmental Conditions on Woody Plant Characteristics
The RF results indicated that seedling and sapling patterns were influenced by CWD, along with other site variables. Slope, slope azimuth, and decay class were the best predictors of woody plant density. Wet bulk density, organic matter weight, and CWD azimuth were the best predictors of woody plant richness. Organic matter depth, soil water content, and organic matter weight were the best predictors of woody plant height. Previous CWD studies have found that site properties can be more dominant factors in determining plant community characteristics, and that CWD may be a minor contributing factor [36,37]. However, in each of the three models, at least two of the top five most influential predictor variables were CWD attributes. Therefore, it could be argued that CWD is nearly as influential as other microsite properties on woody plant density, richness, and height in the regeneration layer.
The Pinus palustris ecosystem in the southeastern USA typically supports some of the lowest volumes of CWD when compared to other forest types and regions [62,95]. Little has been reported on the role of CWD in Pinus palustris woodlands specifically; however, one herpetofaunal study was completed on upland Pinus taeda stands in the southeastern US coastal plain. Davis et al. [96] found that the amphibians and reptiles of the study site were not reliant on CWD for habitat. They concluded that the herpetofauna must have adapted to the characteristically low CWD volumes that are typical of the Pinus palustris ecosystem that historically dominated the study area [96,97]. The results from our study may indicate a similar relationship between CWD and woody plants. It could be that ecosystems in the southeastern US, which do not have historically high volumes of CWD, support woody plant species and animal populations that have largely adapted to survive without CWD. Additionally, our results indicate that CWD inputs may not differentially benefit woody plant species that occur in this system. The variance explained by the RF models for woody plant density, richness, and height were relatively low-not higher than 22%. This low value may indicate that we did not measure particular variables that have a more significant influence on woody plant regeneration. Undocumented soil characteristics may have a greater influence on woody plant assemblages than soil water content or wet bulk density. For example, soil hardness and water repellency can play important roles in woody plant regeneration [15]. Additionally, CWD may have modified soil nutrients that influence woody plant density, richness, or height. Carbon and nitrogen soil content, soil base saturation, and pH may increase with an influx of CWD [25,98,99].

Conclusions
We analyzed woody plant regeneration patterns in a Pinus palustris woodland that experienced a catastrophic wind disturbance. The stands were under a four-year prescribed fire rotation and had been burned twice since the catastrophic event. In this analysis, the best predictors of woody plant density, richness, and height were site characteristics-slope, wet bulk density, and soil water content-and CWD attributes-decay class, CWD azimuth and CWD diameter. Results also indicated that proximity to CWD alone may not drive woody plant regeneration patterns. We suggest that woodland ecosystems located in the humid subtropics, that experience frequent prescribed fire, and that have historically low volumes of CWD support woody plant species that may not be dependent on and do not differentially benefit from CWD. A large influx of CWD has been found to alter carbon and nitrogen storage [100], fungi [101], animals [25], soil fertility [36], and soil erosion [102], and these variables can all influence competitive dynamics and favor some woody plant species over others. However, in the Pinus palustris woodland studied here, pre-disturbance biotic and abiotic conditions seem to be more impactful on woody plant communities than the CWD pulse, such as existing root networks in hardwood species and soil characteristics.
Salvage logging is a common forest operation that occurs after catastrophic disturbances to mitigate the loss of profits from downed or damaged trees and to reduce the risk of insect outbreak and fire [103]. Salvage logging is a contested practice because the ecological consequences are not yet fully understood [104][105][106]. Past research at this study site found that salvage logging did not negatively influence woody plant diversity; however, it did negatively impact ground flora diversity and macrofungal species richness and fruiting abundance [79,101]. Additionally, research conducted in other Pinus spp. stands found that salvage logging can benefit open-habitat avian species [107,108]. The results from our study indicate that retention of CWD after catastrophic disturbance may not alter the regeneration potential of woody plants in Pinus palustris woodlands. The decision to salvage harvest should be based on possible impacts on organisms of management concern. Here, we only documented seedlings and saplings, and we acknowledge that findings are contingent upon the response variables studied [109] and that salvage harvesting influences biophysical conditions other than just CWD loading.
Climate change and other stressors are modifying forest disturbance regimes, and the need to understand disturbance processes and patterns and their influence on forest composition, structure, and function is essential for forest management [110,111]. Information and material legacies from pre-disturbance conditions are critical in driving successional and developmental pathways in developing stands [9,112]. CWD is considered a material biological legacy that can strongly impact forest ecosystem functioning [11,14]. However, in our study we did not find that an influx of CWD significantly influenced woody plant seedling and sapling density or richness in a Pinus palustris woodland. This finding indicates that CWD as a component of ecological memory did not strengthen or weaken resilience in our study system. With a predicted increase in tornado frequency in the southeastern US [35], it is essential to understand how catastrophic disturbances and their legacies will influence resilience of these systems.