Past Management Spurs Differential Plant Communities within a Giant Single-Clone Aspen Forest

Sustainable aspen ecosystems hold great promise for global biodiversity conservation. These forests harbor relatively high species diversity, yet are threatened by fire suppression, land development, timber-focused management, extended droughts, and chronic herbivory. “Pando” is a high-profile quaking aspen (Populus tremuloides) forest in Utah, USA which is putatively the ‘largest living organism on earth.’ Pando comprises an estimated 47,000 genetically identical stems, but is threatened by human impacts. Our interest in the present study is whether changes to the giant organism were affecting understorey vegetation and whether discrete zones are displaying divergent community compositions. For instance, recent research has demonstrated strong herbivory impacts that are affecting portions of Pando differentially. This study consists of 20 randomly distributed vegetation survey plots within three de facto management regimes (hereafter, management group or type) along an herbivory protection gradient: No Fence, 2013 Fence (total protection), and 2014 Fence (imperfect protection). The plant survey was supplemented by previously-established forest and herbivore measurements to test for community assemblage explanatory agents. Sixty-eight species were found across the entire study. Analyses indicated strong links between management group orientation, species assemblages, and tree density/canopy openings. We found distinct evidence that within management group species composition was more similar than across groups for two of the three pairings. However, the other pairing, the most successfully protected area and the completely unprotected area, was not statistically distinct; likely a result a deteriorating overstorey in these two areas, whereas the third management type (2014 Fence) exhibited higher canopy cover. Indicator species analysis found that a small group of plant species had statistical allegiances to specific management groups, suggesting resource preference selection within Pando. Ordination analysis searching for causal factors reached two broad conclusions: (1) aspen regeneration, and therefore long-term resilience, is being negatively affected by chronic animal browsing and (2) current understorey species diversity is highest where forest canopy gaps are abundant. Future research at the massive Pando clone will continue informing linkages between understorey communities and overstorey-driven ecological pathways.


Introduction
Increasingly, forest ecology and management is concerned with ripple effects of process-based strategies as they apply to not only primary forest cover, but to dependent species along trophic through understorey plant communities. We approached this is a stepwise fashion, building from group status, to identifying indicators species, to probing casual factors or key ecological gradients. In this manner, the current work used exploratory methods in community ecology to identify key descriptors of Pando's vegetation status. In order to accomplish this task, we developed a complementary study protocol focused on identifying and assessing plant community make-up as it pertains to forest cover, aspen recruitment, and browsing impacts in this unique setting. The following questions, therefore, constituted the central focus of the current study. (1) Are there group differences in understorey composition between management regimes which may indicate potential diverging pathways? (2) Do individual understorey species display preference for de facto management regimes within the Pando clone? (3) What are the key factors (gradients) underlying plant species diversity? (4) Do some plant species exhibit preference for the most prominent explanatory gradients? Depending on results, the first two questions may lead us to ponder whether diverging understorey assemblages at Pando require key ecological linkages which may now be severed within this genetically uniform forest overstorey. Our expectation is that this massive overstorey clone would support closely aligned plant communities throughout its extent, particularly in a setting such as Pando where elevation and topography vary only modestly (Supplementary material, Table S1). Determination of causal factors for plant community divergences (or convergences) will aid our understanding of burgeoning and future restoration related to Pando's declining state [30]. Outcomes here are expected to add to broader discussions pertaining to clonally dominant tree systems and their international contributions to biodiversity [9, 21,22,32].

Study Site
The Pando aspen clone is found on the Fishlake National Forest in south-central Utah (UTM 434701 E, 4264266 N). The average elevation is 2707 m and slopes are gentle (0-10 degrees) with a generally southeast aspect. Soils are mixed gravelly and cobbly in both loamy A-horizons and clayey B-horizons (Draft Survey, Fishlake National Forest, National Cooperative Soil Survey UT651) originating from tertiary volcanic materials. The forest floor is modestly vegetated, with ample exposed volcanic boulders and bare ground patches. Common juniper (Juniperus communis) and mountain big sagebrush (Artemisia tridentata) are common shrubs across the area. Annual precipitation is attributed mostly to winter snow, with a prominent pulse of summer rainfall (466 mm year −1 precipitation, SNOTEL 1149, 6 years. average).
Pando is located near Fish Lake, a popular resort area. An access road bisects Pando and there are camp sites and private summer homes within the clone. These tourist activities provide indirect protection for browsing ungulates, as hunting of game species is prohibited in such popular areas. Small clearfell-coppice cuts were performed near the edge of the clone in 1987 and 1988, but were unfenced. An additional cut within Pando took place in 1992, though managers at that time elected to fence this treatment due to heavy browsing (and consequently complete forest loss) following the earlier cuts. Currently, domestic cattle (Bos spp.) forage at Pando under a U.S. Forest Service grazing allotment for approximately two weeks annually. Mule deer (Odocoileus hemionus) and Rocky Mountain elk (Cervus elaphus) access this area freely during the typically seven-month snow-free season. Elk appear to be absent in the immediate area based on previous scat surveys within the study area, though we know they are common further from areas of elevated human use. Visitors to the area are generally unaware of Pando's reputation as the 'world's largest organism' [28] due to a lack of signage or other interpretive materials. Thus, while Pando is easily accessible from the access roadway, few people intentionally visit the clone and therefore direct human impacts are minimal.

Field Methods
Data collection for the current study took place over a four year period in two basic sampling efforts: a previously-established network of sample plots to assess base forest conditions and the current understorey vegetation assessment. The original plot selection was based on a 50 m 2 grid of potential sample points with random sub-selection from that matrix. This process resulted in 65 total forest measurement plots across three distinct management regime areas. These areas are referred to as management groups here, although they were not originally designated as experimental zones, but they remain de facto groupings due to the resulting differences in their reproductive success [31]. Those management areas are referred to as: No Fence, 2013 Fence, and 2014 Fence (Box 1; Figure 1). These distinct zones, having previously shown disparate regeneration based on ungulate browsing patterns [31] may also be developing distinct understorey plant communities (the purpose of this study). While the designation No Fence is self-evident, the 2013 Fence (experimental manipulations) and 2014 Fence (no recent manipulations) areas are distinguished not only by recent human disturbance, but by their relative success in protection from browsing. The 2013 Fence was very successful at facilitating aspen regeneration whether plots were disturbed or undisturbed [30], while the 2014 Fence was much less so [31].

Box 1. Treatment Groups.
No Fence-area is unrestricted from browsing herbivores.
2013 Fence-area is restricted from browsing herbivores AND implemented fire, shrub removal, and 50% overstory tree cutting. Approximately half this area was undisturbed as a control.
2014 Fence-area is restricted from browsing herbivores and no recent disturbance has taken place. A previous study, however, determined that mule deer were entering and browsing aspen suckers (cattle had no access).
Forests 2018, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/forests but they remain de facto groupings due to the resulting differences in their reproductive success [31]. Those management areas are referred to as: No Fence, 2013 Fence, and 2014 Fence (Box 1; Figure 1). These distinct zones, having previously shown disparate regeneration based on ungulate browsing patterns [31] may also be developing distinct understorey plant communities (the purpose of this study). While the designation No Fence is self-evident, the 2013 Fence (experimental manipulations) and 2014 Fence (no recent manipulations) areas are distinguished not only by recent human disturbance, but by their relative success in protection from browsing. The 2013 Fence was very successful at facilitating aspen regeneration whether plots were disturbed or undisturbed [30], while the 2014 Fence was much less so [31].
No Fence-area is unrestricted from browsing herbivores.
2013 Fence-area is restricted from browsing herbivores AND implemented fire, shrub removal, and 50% overstory tree cutting. Approximately half this area was undisturbed as a control.
2014 Fence-area is restricted from browsing herbivores and no recent disturbance has taken place. A previous study, however, determined that mule deer were entering and browsing aspen suckers (cattle had no access). Each forest mensuration plot consisted of two fixed area transects (2 × 30 m) aligned perpendicular to each other to capture terrain deviations which may affect site/growth performance. From 2014-2017 field technicians conducted basic site description, mensuration, and reproduction, browse, and animal scat surveys described in detail by Rogers and Gale [30]. In brief, survey crews measured mature tree diameters at breast height (dbh), tree status (live/dead), number of aspen Each forest mensuration plot consisted of two fixed area transects (2 × 30 m) aligned perpendicular to each other to capture terrain deviations which may affect site/growth performance. From 2014-2017 field technicians conducted basic site description, mensuration, and reproduction, browse, and animal scat surveys described in detail by Rogers and Gale [30]. In brief, survey crews measured mature tree diameters at breast height (dbh), tree status (live/dead), number of aspen suckers within height groups, tree and shrub cover, and percent of browsed suckers. Additionally, we counted ungulate pellet groups within fixed area transects. All site, tree, reproduction, browse, and scat data were summarized as per ha values for analysis.
For the current study, a randomly-selected subset of the complete Pando plot network was used to gather new detailed understorey vegetation. With our current sampling regime, we were unable to sample the full mensuration grid of 65 plots due to time and expertise limitations. On 20 plots distributed across the three management zones we identified all plant species and estimated their cover using a dual collection method following procedures of the U.S. Forest Service, Forest Health Monitoring program [33]. Species identification and abundance/cover were obtained on three 1 m 2 quadrats set at 5 m from permanent plot center markers at compass readings of 0 • , 120 • , and 240 • . An additional timed search of the entire vegetation plot (7 m radius) was conducted for species not previously noted in sample quadrats. This combined small-area (quadrats) and larger timed-survey area (plot) two-step process is thought to yield the most thorough species capture [33,34]. Visual estimates of abundance in small-and large-areas were classified in percent cover groups favoring lower coverage accuracy as follows: <1, 1-5, 6-10, 11-25, 26-50, and >50. Abundance data were pooled by species for total sample plot (i.e., quadrats and large-area) using a standardized algorithm to group total abundance ratings into classes specified above, as per the Forest Health Monitoring technique [33]. All vegetation plots were sampled during June 2018 to take advantage of maximum plant expression prior to mid-to-late summer drying and desiccation. Plants not positively identified in the field were confirmed at the Utah State University Herbarium at a later date. Species data were summarized using standard diversity indices: species richness, species evenness, Shannon's diversity index, and Simpson's diversity index [35]. Nomenclature throughout this study followed Welsh et al. [36].

Analytical Methods
Our analytical objective was to comprehensively understand the various management regimes at Pando through the lens of understorey plant communities. We took three distinct analytical approaches: management group differences, individuals species affinity distinctions, and (exploratory) causal factor analysis, which may explain putative differences wrought by management groups. The first two methods employed only plant communities to understand group and species preferences, while the third (causal factor) utilized numerous site, tree, herbivory, and plant community measures. The mechanism for accomplishing this was the PC-ORD (v.7.0, [37]) software package, specifically techniques designed to understand group preferences and explanatory variable "data landscapes" (gradients) using ordination methods [35,37]. In this regard, being a first-time measurement of plant species, data analysis should be considered exploratory in nature (i.e., the study design was not a conventional, hypothesis-driven, framework).
Our initial question was to discern whether species presence and abundance was more similar on plots within management groups as compared to across groups. Multi-response permutation procedures (MRPP) is a nonparametric test for describing within-group agreement of variables in contrast to adjoining data groups. We selected MRPP using the Sørensen distance measure because it is less inclined to exaggeration based on outliers and zero values for cover estimates [38]. All categorical variables were removed from the data matrix except for our group designator, thus this analysis was performed on a matrix of 20 plots by 68 species. MRPP produces a T score indicating the degree of difference between group pairs, an A-value which is the chance-corrected within group agreement (effect size), as well as a p-value establishing level of test significance [35].
The next analytical question involved individual plant species habitat preferences; do certain species favor one management trajectory over another within this genetically uniform forest? Indicator species analysis (ISA) is a multivariate approach testing for no difference between a priori groups in terms of species affinity, or "faithfulness," based on presence and abundance of taxa within groups [39]. We chose ISA to gain a better understanding of species preferences (if any) for current management regimes at Pando. Perfect affinity is defined as always being present and exclusive in the identified group [35]. The statistical significance of the maximum indicator value for each species was tested by 5000 runs of a Monte Carlo randomization procedure. The resulting p-value represents the probability that the calculated indicator value for any species is greater than chance occurrence. Output includes the group for which the maximum indicator value is found, the indicator score for that group, and the associated p-value for each species.
Our final goal (questions 3-4) was to seek explanatory gradients within a comprehensive data set encompassing site description, mensuration, regrowth, browsing, animal presence, and understorey species diversity indices. The incorporation of diversity indices, based on all plant species tallied, differentiated this present analysis from past Pando studies addressing primarily regeneration and animal interactions [30,31]. Diversity indices further provided potential plant community-based indicators at the management group level to distinguish broader diversity patterns within this unique landscape. For exploratory analyses of data gradients addressing causal factors we selected nonmetric multi-dimensional scaling (NMS). NMS is an ordination technique suited for ecological questions where multiple data types are used, normal distributions and variances may be absent, and data sets may contain large numbers of zeros [40]. The aim of NMS is to reduce complex multidimensional data sets to a few key explanatory gradients (axes). The primary matrix for this analysis consisted of 20 plots (rows) pitted against 18 sampling variables (columns). A broader suite of 26 field variables was trimmed down by two categories: raw tally scores (e.g., stem counts) that were redundant to per ha values and location variables (location coordinates and elevation) that detracted from ordination analysis due to their small variation (i.e., having little explanatory value). Two secondary matrices of environmental variables and individual species were plotted over main NMS results in a joint plot graph for display purposes. The lowest stress solution was derived from 250 runs with real plot data. "Stress" is defined as a quantitative assessment of final ordination monotonicity, a measure of how well real data fit the ordination [35,38]. The lowest stress solution was tested with a Monte Carlo run of an additional 250 randomized iterations to evaluate the probability of the final NMS solution being greater than chance (i.e., yielding a p-value). Orthogonal rotation of the final ordination was used to maximize alignment between the strongest environmental variables (i.e., Pearson r values) and primary axes. The lowest number of dimensions (axes) was selected when adding another dimension would have decreased the final stress by <5 [35]. Results were for all tests were considered significant where p < 0.05.

Is Plant Community Composition Different between Groups?
We recorded a total plant diversity (γ) of 68 species, five of which could only be identified to genus, on a total of 20 field plots within the three management regimes of the Pando aspen clone. The average species richness per plot (α) was 23.15 (SD = 3.89). Beta diversity (β = α/γ) at 0.34 represents an estimate of "community turnover" or (in this case) a relatively high species overlap among all sample locations at Pando. Supplementary Tables S1 present mean and variance (SD) values for analysis variables, as well as cover and frequency values for common plant species.
We used multi-response permutation procedures (MRPP) to test for within-group similarity in understorey species composition. Results of this test are shown as an overall score, as well individual pairings of three group combinations (Table 1). Larger negative "T" scores indicate greater separation (dissimilarity) between groups. "A" scores signify a chance-corrected within-group agreement score: 1.0 = perfect agreement, 0 = no better than expected by chance, and an A value of >0.3 is considered a strong indicator of within group species homogeneity [35]. An overall MRPP test result indicates strong within group homogeneity, while comparisons between the 2014 fenced plots and the other two groups also yielded significant within group agreement (Table 1). Understorey vegetation assemblages, when compared between the No Fence and 2013 Fence groups, displayed no significant difference.

Species Preference for Management Groups
Of the 68 distinct species represented in our understorey plant survey at Pando only seven showed a solid affinity for particular management groups ( Table 2). Indicator species analysis results found that four species met a higher significance threshold (<0.05 p-value) for group allegiance. Populus tremuloides (understorey only) and Potentilla hippiana favored the 2013 fenced area and Calamagrostis spp. and Lupinus argenteus were statistically associated with the 2014 fenced group. Additionally, three other species fell just below the significance threshold and were also only statistically aligned with the fenced area survey groups (Table 2). No species showed ISA preference for the No Fence area, suggesting a generalist plant assemblage or perhaps areas without species-specific habitat.

Exploring Factors Affecting Species Make-Up and Habitat Preferences
Ordination analysis resulted in a two-dimensional solution on a matrix of 20 plots by 18 variables. We projected the results of the NMS in ordination "plot-data space" as a joint plot highlighting environmental variables as overlay vectors (Figure 2a). The final NMS solution produced a stress value of 8.84 with an instability of 0.00. A Monte Carlo test of 250 random data runs versus the real data set verified a significant NMS outcome (p = 0.03). A two-axis solution described about 87% of ordination variance (axis 1: r 2 = 0.489; axis 2: r 2 = 0.377; and orthogonality = 99.8). Length and direction of vectors corresponds to environmental variable strength and relationship to the two-dimensional plot-data space. Only environmental variables making the strongest contributions to data ordinations are shown (i.e., r 2 ≥ 0.3) in the joint plot (Figure 2a). Table 3 presents NMS results by axes for all environmental variables. Most measures tested in this ordination contributed very little to environmental gradients (NMS axes) at Pando. Strong positive and negative responses to axis 1 (Figure 2a) represent factors working in opposition in terms of their influence on the totality of plot variable values at Pando. Axis 1 describes a gradient of regeneration abundance with a negative correlation to deer and cattle presence. This result is nearly identical to that reported previously [31]. Axis 2 illustrates an inverse correlation between aspen cover (collectively) and understorey species diversity. A moderately strong correlation of aspen recruitment to axis 2 may present a spurious result given the low overall recruitment occurrence among our sample sites (with the exception of a single anomalous plot in a 1992 historic clearfell-coppice cut area [30]). Figure 2b projects species vectors as an overlay on the main ordination to gauge how particular plants align with explanatory variables. Only those species with r 2 > 0.3 are shown. While most plants showed no particular statistical relation to prime NMS axes, P. tremuloides, Cirsium clavatum, Festuca thurberi, and Penstemon watsonii displayed robust correlations (r > 0.5 or < −0.5) to these environmental gradients (Table 4). single anomalous plot in a 1992 historic clearfell-coppice cut area [30]). Figure 2b projects species vectors as an overlay on the main ordination to gauge how particular plants align with explanatory variables. Only those species with r 2 > 0.3 are shown. While most plants showed no particular statistical relation to prime NMS axes, P. tremuloides, Cirsium clavatum, Festuca thurberi, and Penstemon watsonii displayed robust correlations (r > 0.5 or < −0.5) to these environmental gradients (Table 4).

Group and Species Preferences in a Single-Genotype Forest
The work presented here constitutes the first complete, or baseline, survey of plant species at the iconic Pando aspen clone in central Utah, USA. The massive Pando aspen clone has likely survived for millennia [41], though as a result of recent human actions there are clearly at least three distinct plant communities supported under a genetically uniform forest canopy (Table 1). Our results also indicate distinct plant communities, interestingly, within separate fenced exclosures at Pando. Through the lens of understorey community "loyalty" (MRPP), we note that the 2014 fenced plots were statistically more similar to each other than they were to either No Fence or 2013 Fence locations (Table 1). We did not find a group distinction when the No Fence and 2013 Fence plant communities were compared to each other. Given the high rate of aspen reproduction now occurring within the 2013 Fenced exclosure [30], after only five years we must conclude that early-stage regeneration is not yet affecting plant community composition. However, light exposure to the forest floor appears to be a stronger causal factor in plant make-up at this point (Figure 2a). Our results may be describing a secondary gradient of pioneer plants (high diversity and light) favoring current canopy gaps and younger forests, while a shade-tolerant understorey community is tracking greater tree cover, although such conclusions seem speculative given the preliminary nature of this study.
Four species showed clear group preference (fidelity) within management areas at Pando (Table 2) and an additional three species also seem to favor specific conditions although they fell just beyond our declared statistical threshold (p < 0.05). Based on previous work at Pando [30,31], it is no surprise that P. tremuloides in the understorey is thriving in the 2013 Fence management area. This finding, however, is likely due more to the success of the early exclusion of browsers via fencing, then to any strictly ecological factor. Welsh et al. [36] suggest that Potentilla hippiana, Bromus ciliatus, and Artemisia tridentata favor meadows, open slopes, and partially forested (high light) habitats; thus, this fits the relatively open canopy of the 2013 Fence regime. However, similar adjacent habitat found in the No Fence zone did not show similar species fidelity leading us to wonder if perhaps low browsing level or rapid aspen regeneration are facilitating improved habitat for these species? It will be interesting to track how species composition changes over time with increased vertical complexity in recruitment facilitated by fencing. While two species favoring the 2014 Fence area have broad tolerance of light (Calamagrostis spp. and Lupinus argenteus), Silene menziesii is known to be a more shade-preferring species in montane settings [36]. Higher aspen canopy cover [31] in the 2014 Fence, alongside moderate correlations between Silene menziesii and Axis 2 ( Table 4) at least partially corroborate this ISA group fidelity for shaded habitat.
Both MRPP and ISA are analyses driven solely by plant community make-up. They do not attempt to address covariation or more plainly, causation. Thus, we turned to NMS as a powerful ordination tool which combines multiple factors, alongside plant community composition, to explore causality. However, seemingly divergent results using two different analytical approaches here requires further explanation. Specifically, different species are highlighted in results of indicator species analysis (ISA; Table 2) vs. nonmetric multidimensional scaling (NMS; Figure 2b, Table 4). ISA is more pointedly directed toward discriminating between management groups based on plant communities only. While NMS (Figure 2b) may suggest alignment of certain groups in corroboration of ISA results (Table 2; i.e., P. tremuloides and Bromus ciliatus), the NMS ordination is a more comprehensive analysis taking into account an array of potential explanatory factors. In Figure 2b, the 2013 Fence plots group nicely, but the other two management group plots do not. This explains the poor performance of the variable "Fence" in the NMS final result (Table 3). Other individual species with significant correlations to the NMS (Figure 2b) align, to varying degrees with key gradients (axes), but those alignments are independent, for the most part, of group preferences exhibited in ISA results. Thus, Festuca thurberi seems to show preference for areas with higher herbivore frequency, Aquilegia caerulea and Bromus ciliatus thrive in high aspen regeneration locations, and Penstemon watsonii performs best where diversity indices are elevated (i.e., high light environments). Overall, then, we must view ISA and NMS as presenting two independent looks at habitat preferences within Pando's disparate management regimes.

Key Factors Influence Pando's Plant Make-Up
Human interventions, whether intended or not, have negatively affected Pando's ability to reproduce and thrive [30,31]. Stable aspen forests sustain themselves via relatively continuous (i.e., not episodic) recruitment; in fact, regular recruitment and a multi-story stand structure are thought to be key components of healthy stable aspen types [1,4,6,42]. Past research has described a rapidly deteriorating aspen clone resulting from increased deer and cattle herbivory, compounded by other human developments, such as added home sites, a recreational campground, and increased road traffic [35]. These features, collectively, have increased human presence and decreased mule deer hunting in a localized area to the point that the Pando grove has become a refuge from human predation, or safe-zone, for deer. Past work has shown that higher presence of deer scat related closely to higher browse level and reduced recruitment [31]. Existing monitoring data suggests that this situation is a prime cause of concentrated, unfettered, herbivory leading to grossly unbalanced aspen stem demographics: while mature/dying aspen are abundant, there are few recruitment stems to replace the decaying overstorey. This basic formula is clearly illustrated by key indicators opposing each other within our NMS ordination; regeneration was shown to increase on all plots where deer and cattle scat were less abundant (Figure 2a, Table 3). While this is the strongest explanatory factor presented here, it was also the central issue revealed using this data set in a previous study [35]. Browsing herbivore effects on aspen recruitment, however, are broadly implicated in the distinct understorey plant assemblages.
The addition of a detailed understorey plant survey has allowed us to not only establish group and species habitat preferences (Tables 1 and 2), but to incorporate common diversity metrics into a broader assessment of the Pando ecosystem (Tables 1 and 4; Figure 2a,b). We conclude that rapid die-off of the forest overstorey appears to be creating broad canopy openings that support relatively high understorey plant diversity. In contrast, a lack of regeneration alongside relatively high canopy cover, basal area, and total tree counts is reflected in lower diversity scores, in particular Simpson's diversity index (axis 2, Table 3; Figure 2a and Table S1). Other diversity indices corroborate a positive relationship to greater canopy openness, but our results correlate evenness negatively to Axis 1 in association with less regeneration and more herbivore presence. This evenness result is somewhat confounding to other diversity indices, but may suggest greater disparity in abundance (not richness) of species with a recent surge in aspen regeneration, most pointedly within the 2013 Fence regime group.
In terms of the explanatory value of our results, there is still much to be explored. The period since recent fencing has been erected, while inadvertently initiating divergent pathways, is likely too short for establishment of stable understorey communities (i.e., a composition in transition). Many species present, and perhaps some not yet established, may be displaying ambiguous habitat preferences at this stage within rapidly changing (e.g., 2013 Fence) environments. We also must consider the perspective of microhabitat variability-in soil types, soil depth, or terrain-inherent within even a genetically homogenous aspen clone. How would a map of native understorey variability be projected over the current human-induced larger scale compositional groupings now evident at Pando? While it may be difficult to thoroughly detect microhabitat variability with the current monitoring approach, we are still left with important questions unanswered regarding species diversity in this unique forest. Beyond plant habitat considerations, it is likely that deer (and cattle outside fences) are making preferential food choices regarding plant consumption which are affecting current understorey communities. For instance, it has been suggested that indicator plants favoring the 2014 Fence zone may be "deer resistant" and therefore are defining that (now) unique community. This assertion, if true, adds complexity to our simpler model of aspen tree canopy and structure being the key driver in defining Pando's understorey composition. While such knowledge gaps are not trivial to researchers, they also present a solid argument for management actions respective of the entire Pando system, inclusive of currently unknown explanatory features of plant diversity.

Broader Implications for Sustained Aspen Diversity
Much has been written about trophic implications of overabundant herbivores (e.g., [24,43]), and there is recognized disagreement on explanatory factors leading to system collapses or rebounds [44][45][46] that may affect aspen ecosystem diversity. At landscapes scales, diversity in aspen canopy conditions (i.e., successional stage, canopy cover, drought response, and human alterations) are expected to have concomitant implications for understorey species [9,18]. Here we examined a limited forest landscape, ostensibly one of unique tree species and genotypic uniformity that would be expected (if ever were the case) to have minimal variation in tree cover and thus, dependent plant communities. It may be that variation within even this limited expanse of forest encourages similar diversity, although perhaps not on the basis of spatially distinct understory plant communities. In contrast, we would anticipate that a resilient stable aspen forest (i.e., one reproducing continuously) would exhibit small-patch dynamics rather than large, human-induced, distinct vegetative communities, as appears to be the current situation. In fact, the present inadvertent management zones-different conditions arising from human actions affecting growth and survival-also appear to depict varying degrees of community separation based on plant composition and abundance (Tables 1 and 2). Other researchers have similarly found that herbivore effects on aspen are consequential to plant [9,16], lichen [47], insect [11], and avian diversity [24], but none have examined such factors in a genetically uniform aspen forest.
Lessons gleaned from this small-scale, partially controlled (i.e., a single genotype), experiment have wider ramifications for aspen ecology and management broadly. Still, some caution is warranted in projecting our findings too far afield, given the unique nature of this large single-genotype forest. With this disclaimer broader findings from this study are as follows. First, aspen forests have traditionally been treated uniformly in terms of management prescriptions. A 1992 clearfell-coppice cutting, while appropriately protected with fencing, has resulted in a single-storied forest today. Future management actions within Pando, using the wisdom of hindsight alongside results here, would do well to manage stable aspen within functionally appropriate aspen prescriptions [4]; meaning selection and group tree felling would promote more resilient uneven-aged forests. Second, herbivory is now a common threat to survival of aspen forests, not only in North America, but across Europe as well [10,32,48,49]. Human desire for large ungulates for subsistence, aesthetics, or recreational purposes, while simultaneously eliminating key predators, in some locales, has resulted in runaway populations and concomitant vegetation impacts [43], such as the changes in understorey vegetation implicated here. In areas, undergoing rapid climate warming such as the broad region of the current study, combined effects of extended drought and herbivory are particularly detrimental to sustained aspen communities [19,24,49]. Third, we have seen here that understorey plant assemblages (and presumably dependent fauna) will be altered by aspen's response to human decisions-in cutting, recreational use, herbivore management, and other developments (i.e., roads, houses, campgrounds, and trails)-in the form of canopy openings and varying degrees of recruitment success. Human activities, overall, may also act as vectors for exotic species introductions that can have lasting effects on biodiversity in aspen communities [9].
A key resource value of aspen forests is their well-known ability to harbor a large number of plant and animal species [21,23]. In the current study area, as well as in aspen communities at scales from local to continental, we should anticipate that inappropriate management actions will have cascading effects on a rich and obligate biota. This study, focused on the single-genotype Pando forest, suggests that aspen ecosystems at-large will benefit from functionally appropriate management, prioritization of monitoring (particularly for herbivore impacts), and an understanding of critical linkages between sustainable overstorey conditions and understorey composition.

Conclusions
Previous works have demonstrated marked changes in the famed Pando aspen clone in recent decades [30,31]. We set out here to understand whether alterations in aspen overstorey and reproduction were affecting dependent plant communities. Our results indicate strong links between this unique, genetically uniform, forest and understorey species responses to recent de facto management regimes reflecting levels of canopy coverage and herbivory. We found distinct evidence that within management group species composition was more similar than across groups for two of the three pairings. Interestingly, there was greater similarity between the most successfully protected area and the completely unprotected area; this is likely a result of more open canopy covers in these two areas, whereas the third management type (2014 Fence) displayed the greatest overstorey shading. Indicator species analysis found that a small group of plant species showed statistical affinities for specific management groups, suggesting species resource preferences currently exhibited in the respective areas within Pando. Finally, an exploratory ordination analysis searching for causal factors resulted in two broad conclusions: (1) aspen regeneration, and therefore long-term resilience, is being negatively affected by chronic animal browsing and (2) current understorey species diversity is highest where forest canopy gaps are more abundant (i.e., where there is greater tree canopy closure, basal area, and total number of trees, there is less plant cover and richness).
We are interested in how plant species will be influenced by changing conditions at Pando over longer periods. Initial results reported here suggest Pando is moving toward a future of divergent ecological pathways despite its genetically uniform overstorey. This study constitutes a first approximation using exploratory analytical methods; repeat measures may show continued community divergence, or perhaps a convergence as structural complexity is restored at Pando. The current trajectory of increased mortality and absent successful (overall) recruitment bodes poorly for Pando's future. While greater light availability may temporarily increase plant diversity, the long-term impacts of a collapsing aspen cover are expected to decrease overall biodiversity. As a baseline study, we have documented here certain patterns that are expected to change as this forest declines or flourishes in the coming decades.
Author Contributions: P.C.R. conceived of the study, collected field data, obtained funding, performed analysis, and co-wrote the manuscript. J.Š. collected field data, identified plant vouchers, performed analysis, and co-wrote the manuscript.