PLoS ONE
Home Combining botanical collections and ecological data to better describe plant community diversity
Combining botanical collections and ecological data to better describe plant community diversity
Combining botanical collections and ecological data to better describe plant community diversity

Competing Interests: The authors have declared no competing interest.

Article Type: Research Article Article History
  • Altmetric
Abstract

In this age of rapid biodiversity loss, we must continue to refine our approaches to describing variation in life on Earth. Combining knowledge and research tools from multiple disciplines is one way to better describe complex natural systems. Understanding plant community diversity requires documenting both pattern and process. We must first know which species exist, and where (i.e., taxonomic and biogeographic patterns), before we can determine why they exist there (i.e., ecological and evolutionary processes). Floristic botanists often use collections-based approaches to elucidate biodiversity patterns, while plant ecologists use hypothesis-driven statistical approaches to describe underlying processes. Because of these different disciplinary histories and research goals, floristic botanists and plant ecologists often remain siloed in their work. Here, using a case study from an urban greenway in Colorado, USA, we illustrate that the collections-based, opportunistic sampling of floristic botanists is highly complementary to the transect- or plot-based sampling of plant ecologists. We found that floristic sampling captured a community species pool four times larger than that captured using ecological transects, with rarefaction and non-parametric species estimation indicating that it would be prohibitive to capture the “true” community species pool if constrained to sampling within transects. We further illustrate that the discrepancy in species pool size between approaches led to a different interpretation of the greenway’s ecological condition in some cases (e.g., transects missed uncommon cultivated species escaping from nearby gardens) but not others (e.g., plant species distributions among functional groups were similar between species pools). Finally, we show that while using transects to estimate plant relative abundances necessarily trades off with a fuller assessment of the species pool, it is an indispensable indicator of ecosystem health, as evidenced by three non-native grasses contributing to 50% of plant cover along the highly modified urban greenway. We suggest that actively fostering collaborations between floristic botanists and ecologists can create new insights into the maintenance of species diversity at the community scale.

Alba,Levy,Hufft,and Wan: Combining botanical collections and ecological data to better describe plant community diversity

Introduction

What is the minimum sampling effort needed to adequately document plant community richness and composition? This question forms a fulcrum upon which multiple branches of plant science have revolved for decades [14]. We still do not have a universal solution, because the answer depends both on the characteristics of the sampled community and the goals of the researcher [5]. Further, there are several dimensions of biodiversity including richness, abundance, and evenness, each of which can manifest differently across temporal and spatial scales, as well as study systems [6,7]. Given the complexity, and increasing urgency, of describing Earth’s biodiversity, it is necessary to continue refining our sampling approaches.

Combining knowledge and research tools from multiple disciplines is one way to better describe complex systems [810]. One instance where collaboration remains elusive is between floristic botanists and plant ecologists. While these groups flank each other on the spectrum of biodiversity scientists, they are often siloed, in part because of their different disciplinary histories and research goals [11]. Here we explore the unique research lenses and sampling approaches that floristic botanists and plant ecologists use to describe plant community diversity. We then illustrate how these different approaches are complementary in describing diversity using research conducted along an urban greenway. We close by discussing the circumstances under which collaboration is likely to be most beneficial in this time of rapid biodiversity loss.

Understanding plant community diversity requires documenting both pattern and process. We must first know which species exist, and where (i.e., taxonomic and biogeographic patterns), before we can determine why they exist there (i.e., ecological and evolutionary processes). Botanical specimens, which are routinely collected by floristic botanists, form the backbone of what we know about plant taxonomy and biogeography, from species discovery to the generation of exhaustive and meticulously vouchered local and regional floras [12] (Fig 1). While such data can in aggregate be used to test hypotheses about what shapes biodiversity, the primary data (i.e., physical specimens and species lists) are not designed to uncover site-level ecological processes. In contrast, ecological data are collected for the express purpose of answering a question or testing a hypothesis [13]. This often requires statistical design and analysis, including the use of replicated transects or plots, which are distributed in a manner that reduces biases in the resulting data [5,14,15]. Unavoidably, this requirement drastically shrinks the area that plant ecologists can sample [16], which has led to a multitude of papers concerned with determining “how much sampling is enough?” (with “enough” often approximated by the asymptote on a species accumulation or rarefaction curve) [1720].

Conceptual diagram illustrating the different research foci, sampling approaches, and research outputs of floristic botanists and plant ecologists, placed within the context of the species-rank abundance curve.
Fig 1

Conceptual diagram illustrating the different research foci, sampling approaches, and research outputs of floristic botanists and plant ecologists, placed within the context of the species-rank abundance curve.

This curve typifies many plant communities in which a few species are common, while many species are rare. Plant ecology research often focuses on common species that drive ecological processes. Describing these processes requires transects or plots that allow for statistical hypothesis-testing to answer specific questions. This approach necessarily shrinks the area that can be sampled, causing species far out on the tail of rarity to be missed. Conversely, floristic botany focuses on species discovery, taxonomy, systematics, and the building of floras, which do not rely on statistically designed field sampling. This frees floristic botanists to search exhaustively for unique species in a study area, moving them farther out on the tail of rarity for a better estimate of a site’s “true” species pool. Each approach yields unique outputs. For example, the ecological approach may quantify contemporary relationships between plant relative abundances and their environment (top left of figure), while the floristic approach generates specimens, phylogenies, or floras with species distribution maps (top right of figure). The research outputs of each field could be more mutually informative (cycling arrows), e.g., if a collected specimen was linked to a plot-based estimate of that species relative abundance at the collection location.

Floristic botanists, quite powerfully, are not tied to plots. They are free to make informed decisions about how much area must be covered to adequately sample a plant community of interest, which can often be completely traversed using an opportunistic sampling approach. Thus, botanical collectors can cover a considerably larger area than ecologists, thereby finding species further out on the tail of rarity (Fig 1) and generating a better estimate of the species pool in a given study area. Further, because botanists collect high-veracity specimens that are in flower or fruit, they usually have the structures needed for a species-level identification. Conversely, as part of their effort to reduce sampling biases, ecologists must identify everything captured in a plot regardless of phenological stage, potentially reducing the veracity of the identification or constraining it to the genus level (especially for graminoids). Still, ecologists typically capture the most abundant species in a community, with such species often driving ecological processes including nutrient cycling and response to disturbances like fire and grazing [21]. Importantly, only by using plot- or transect-based approaches can researchers move beyond documenting plant presence (richness) to quantitatively relating plant relative abundances to environmental variation.

Here we illustrate that the sampling approaches of botanists and ecologists are complementary, with the former providing critical context about the community species pool, and the latter layering on quantitative data for a subset of abundant species. We demonstrate the utility of implementing both approaches with a case study from an urban greenway in Colorado, USA. First, we compare species richness estimates generated using opportunistic sampling to those generated using ecological (line-point intercept) transects, assuming that the intensive floristic search will best approximate the “true” community species pool (sensu [28]). We then use rarefaction and species richness estimation to determine how closely the transect-based species-pool estimate reflects the larger species pool. This opportunistic versus bounded comparison of species pools is unusual (although not unique; see [5,22] for examples with mosses and lichens), with most vascular plant studies instead using plots of increasing size to estimate species pools [23]. Second, we explore whether the composition of the species pools captured by each approach leads to different ecological interpretations of the greenway’s flora. To this end, we compare the taxonomic coverage of each approach, as well as how species are distributed in relationship to functional group, biogeographic origin, floristic quality (based on Coefficients of Conservatism), and Wetland Indicator Status. Third, we integrate abundance estimates from the ecological transects into our interpretation of the greenway’s ecological condition to show that abundance is a desirable, if not necessary, complement to richness in any setting with applied conservation or management goals. Finally, we highlight areas of opportunity for collaboration between floristic botanists and plant ecologists.

Methods

Study system and field sampling

We sampled plant communities along the High Line Canal greenway, a 66-mile recreational trail that passes through 11 municipalities in the Denver-Metro area of Colorado. The trail runs alongside a 71-mile earthen canal (owned by Denver Water), which was excavated in the late 1800s to support agriculture and human settlement in what was historically native plains and foothills shrubland vegetation (see S1 Fig for a map of the greenway in relationship to EPA Ecoregions). The greenway thus represents a human-created waterway that is highly managed, yet supports a species pool that contains native flora (see [24] for habitat descriptions).

The greenway’s length (with the Canal’s inception located at 39.48362, -105.11293) is demarcated by mile markers that we used to generate a random subset of 45 locations at which the botanical and ecological field crews could synchronize their data collection. As is typical for collections-based floristic surveys, the botanical crew sampled exhaustively from early spring (7 May) through late summer (28 September) to capture early-, middle-, and late-season species. Starting at each of the 45 mile markers, the crew walked in the Canal’s downstream direction, searching the greenway for newly encountered species to collect and accession to the Kathryn Kalmbach Herbarium (KHD) at Denver Botanic Gardens [25]. Permission to collect plant specimens was provided by Denver Water. Most mile marker locations were sampled once during the inventory, but a few were revisited if they occurred in a vegetation type that would not be re-encountered later in the season at the other mile markers (e.g., mile markers zero and one at the inception of the Canal were the only locations in the foothills shrubland Ecoregion; S1 Fig).

We used a staggered sampling design in which 5 mile markers spanning the southwestern to northeastern extent of the Canal were sampled every other week from May to September. The floristic survey was carried out over 57 days, comprising 850 search-hours and an estimated distance covered of 42 miles (calculated from our daily starting and stopping waypoints logged with a GPS unit; S2 Fig). The botanical crew consisted of two botanists trained in the local flora and one to two additional non-botanists who assisted with specimen collection. All members of the crew searched for species within an ~50 to 75-foot-wide viewshed moving from the bed of the Canal, up the Canal bank, across the greenway trail, and over to the property line that marked the end of Denver Water’s ownership (S3 Fig). High-veracity (with identifying structures) herbarium specimens were accessioned for every species encountered during the floristic survey (numbering 1570 specimens, including duplicates; collections data available; 26). Identifications were made using [2729].

The ecological sampling was carried out over 10 days, from May 22, 2018 to June 6, 2018, to capture a snapshot of plant communities around peak biomass. This method of deploying a concerted sampling effort over a short time period is common in ecological sampling, because it is often of interest to detect treatment differences that could be obscured by confounding time lags between sample dates (as opposed to the floristic botany goal of exhaustively delineating a species pool over time). At each of the 45 miles markers, we laid a 12 m × 2 m transect, the length of which captured habitat variation across the greenway corridor (S3 Fig). We used the line-point intercept method [30] to make field observations of plant species presence every 0.25 m along the 12 m transect (as well as bare ground, plant litter, and rocks, which we do not report herein) [31]. In the associated data set [32], the “first hit” was used to generate the reported percent cover estimates (number of hits per species per total number of hits), while the “second hit” was used to add species to our presence list. We also searched each of the two, one-meter-wide belt transects for additional species that were not encountered along the line-point transect. Voucher specimens were collected for the species encountered during the observational ecological sampling (collected outside the transects so as not to influence long-term sampling). However, given the short time period of the ecological sampling, not all specimens had flowers or fruits, and therefore were not of sufficient quality to be curated. All specimens were kept during the field season and subsequent analyses to facilitate identifications, but only higher quality specimens were accessioned to the herbarium [collections data available; 32]. Please note that one example specimen exists for potentially hundreds of field observations (i.e., each time a species was encountered along the transects).

We chose the line-point method as the most appropriate for our system, with its narrow and steep canal bank that could not accommodate other plot designs. Additionally, our aim with the ecological transects was to estimate not only species presence, but also composition. For questions about composition, the line-point method is highly repeatable across individuals and rapidly deployed, thereby maximizing sampling replicates across many locations in a single season [30]. Such transects will capture fewer species than other methods (e.g., Modified-Whitaker plots); however, any bounded sampling approach will cover considerably less area than can be achieved with opportunistic sampling based in the floristic tradition of using the habitat itself as the sampling unit [5]. Importantly, we note that it was not our goal to equalize the temporal or spatial scales of the two sampling approaches (which in our experience is not often done in practice), but rather to sample in a manner that is broadly consistent with collections-based versus ecological disciplines. The particulars of our comparison, such as the sampling window and the use of transects rather than any number of plot types, contextualize the results.

Community species pools and ecological metrics

We first estimated the species pools captured by the collections-based floristic and quantitative ecological sampling approaches, and then compared the pools in ecologically meaningful ways. Species pools are hierarchical and scale dependent, having been variously defined, but they are typically partitioned from larger scale (regional pools), to mid-scale (local pools), to smaller scale (actual or community pools). Here, we define a community species pool according to [33] as “the set of species present in a target community,” with our target community being the urban greenway. We consider the species found during the intensive floristic sampling as a best estimate of the greenway’s “true” community species pool and expect that the ecological transects will capture a subset of this larger pool. (We acknowledge that even the intensive floristic sampling will not capture the true pool, but the goal is to employ realistic sampling schemes used in botanical floristics and plant ecology to see how they compare).

After delineating the species pools, we chose several ecologically informative metrics to compare them: species distributions among families; plant functional group based on longevity and growth form; floristic quality (based on Coefficients of Conservatism or C values) [34]; Wetland Indicator Status [35]; and biogeographic origin (assigned as native or introduced to Colorado according to [27], and if introduced, whether it is cultivated). Plant functional groups are extensively used to aggregate large numbers of species into a few classes that are expected to respond similarly to changes in their environment, or to similarly affect their environment [3638]. Floristic Quality Analysis is often used in the conservation realm to assess an area’s conservation value [34]. Sites with high floristic quality are relatively pristine, having departed little from the disturbance regime that existed prior to European settlement. Related to this, species in a community can be ranked on a scale of zero to 10 according to their “conservatism,” or their fidelity to habitats that are more (or less) degraded by human use. Species that can only persist in undegraded, native habitats are assigned high C values, while ruderal species, which can withstand substantial degradation, are assigned low scores (rankings below adapted from [39]).

    0–3: Introduced species (always = 0), plus native species that occur in moderately to highly degraded sites (1–3)

    4–6: Native species that show some affinity to natural areas and are often dominant or are present across a wide range of habitats and environments

    7–8: Native species associated mostly with natural areas, but that can sometimes persist in degraded habitat

    9–10: Native species that tolerate very little or no habitat degradation

Wetland Indicator Status ranks species according to their dependence on saturated soils, or wetland conditions. This metric is meaningful in our study system because the Canal represents a novel waterway in an otherwise semi-arid landscape. The status rankings are as follows. 1) Obligate Wetland: almost always a hydrophyte, rarely found in uplands; 2) Facultative Wetland: usually a hydrophyte but occasionally found in uplands; 3) Facultative: commonly occurs as either a hydrophyte or non-hydrophyte; 4) Facultative Upland: occasionally a hydrophyte, but usually occurs in uplands; 5) Upland: rarely a hydrophyte, almost always found in uplands [40]. Our study area occurs on the interface of the Western Mountain Valleys and Coasts and Great Plains regions, which can occasionally have different wetland indicators assigned for the same species. If there was a discrepancy between regions, we chose the more hydrophilic option to generate a conservative list in terms of species reliance upon water.

Data analysis

We calculated the expected number of species in our pooled samples using species rarefaction. The rarefaction curve was produced from 1000 random resamples drawn without replacement from the pool of species in the transects. We extrapolated out to 1.5X the original sample size (68 transects), as extrapolation past doubling or tripling of the reference sample size is not recommended due to increased uncertainty [41]. Asymptotic species richness was estimated using the Chao2 estimator in EstimateS version 9.1.0 [42] using the incidence of each species within each sampling transect.

We used Pearson’s Chi-square tests of independence ([43]; chisq.test in R v. 3.6.2) to explore whether the floristic and ecological sampling approaches generated different distributions of species in relationship to taxonomic coverage, functional groups, C values, Wetland Indicator Status, and biogeographic origin. For each data set (floristic and ecological) we summed species frequencies from the 45 mile markers to alleviate low cell counts within mile markers. When the sampling approaches generated significantly different (P ≤ 0.05) distributions of species among groups, we used the adjusted standardized residuals to assess which groups contributed to the disparity (with residuals exceeding an absolute value of approximately 2 considered important; [44]). We could not assign functional groups to nine species that were only identified to the level of genus. We also could not assign C values to a subset of species that did not have them available (n = 18 species or 4% of all collections and six species or 5% from transects). The same was true for Wetland Indicator Status, in particular for cultivated species, which are not assigned this type of indictor (n = 76 species, or 17% of all collections, and 10 species, or 8% from transects).

For the purpose of assessing species composition along the greenway, we calculated abundance using percent cover from the ecological transects. Abundance estimates were calculated by taking the number of hits per species divided by the total number of hits sampled over the entire Canal [raw data available; 31].

Results

We found 452 species using the opportunistic sampling approach used in floristic botany [collections data available; 45] and 126 species using ecological transects (see S1 Table for full species list). Species richness modeled from the transect data underestimated the floristic estimate by 41% (mean = 184; 95% CI lower bound = 151; 95% CI upper bound = 253; Fig 2). This marked underestimate manifested despite using the Chao2 estimator, which statistically accounts for the fact that uncommon species will likely be missed. The species rarefaction curve showed that a 51% increase in our sampling effort, from 45 to 68 transects, would only capture 31% of the species pool observed during the opportunistic sampling (mean = 141; 95% CI lower bound = 125; 95% CI upper bound = 157; Fig 2). While neither the species rarefaction curve nor the Chao2 estimator reached a definitive asymptote (although they were distinctly leveling off), our transect sampling effort was based on what we could reasonably achieve with the available financial and personnel resources. We expect that other researchers are similarly constrained in most applied situations.

Species richness estimated using transect-based rarefaction (black line shown with 95% CI in red shading) based on 45 sampled transects and extrapolated up to 68 transects (values after the vertical black line).
Fig 2

Species richness estimated using transect-based rarefaction (black line shown with 95% CI in red shading) based on 45 sampled transects and extrapolated up to 68 transects (values after the vertical black line).

An asymptotic richness estimate was calculated using the Chao2 estimator (grey line with 95% CI in green shading).

The taxonomic coverage of the ecological versus botanical sampling became less complete moving hierarchically from family, to genus, to species (with transects capturing 50% of families [39 versus 78], 36% of genera [103 versus 289], and 28% of species [126 versus 452]; see S4 Fig for species distributions among families). The distribution of species among functional groups was statistically similar for the two sampling approaches (Fig 3A; χ2 = 5.6, df = 7, P = 0.59), although perennial grasses were weakly over-represented (by 6%) along transects, reflecting the pattern found for family distributions (see Poaceae, S4 Fig). The floristic quality of the greenway, based on the distribution of species among C values (Fig 3B; χ2 = 9.9, df = 10, P = 0.45), did not differ significantly between the two sampling approaches. Wetland Indicator Status significantly differed (Fig 3C; χ2 = 10.6, df = 4, P = 0.03), with transects underestimating obligate wetland species by 5% (adjusted residuals = -2.05) and upland species by 10% (adjusted residuals = 2.0) relative to the larger species pool. In terms of biogeographic origin, the transects over-estimated (although non-significantly) the proportion of introduced species (55%) relative to the opportunistic sampling (46%; χ2 = 2.75, df = 1, P = 0.10). The composition of introduced species captured by the two approaches differed in an important respect, with the floristic approach capturing a substantial number of uncommon, non-native garden cultivars (79 cultivars out of 208 introduced species, or 38%) that the transects missed (four cultivars out of 68 introduced species, or 6%).

A-C. Distribution of species captured by the floristic botany (collections) and ecological (transects) sampling approaches in relationship to functional groups, floristic quality (based on Coefficients of Conservatism or C values), and Wetland Indicator Status. The distribution of species among functional groups and C values did not significantly differ, while the distribution among wetland indicator status differed. Functional groups: PF = perennial forb; W = woody (and perennial); AF = annual forb; PG = perennial grass; VFG = variable forbs and grasses (annual to short-lived perennials; includes only two species of grasses); AG = annual grass; BF = biennial forb. Wetland Indicator Status: OBL = obligate wetland; FACW = facultative wetland; FAC = facultative; FACU = facultative upland; UPL = upland. See text for statistics and definitions of C values and wetland indicators.
Fig 3

A-C. Distribution of species captured by the floristic botany (collections) and ecological (transects) sampling approaches in relationship to functional groups, floristic quality (based on Coefficients of Conservatism or C values), and Wetland Indicator Status. The distribution of species among functional groups and C values did not significantly differ, while the distribution among wetland indicator status differed. Functional groups: PF = perennial forb; W = woody (and perennial); AF = annual forb; PG = perennial grass; VFG = variable forbs and grasses (annual to short-lived perennials; includes only two species of grasses); AG = annual grass; BF = biennial forb. Wetland Indicator Status: OBL = obligate wetland; FACW = facultative wetland; FAC = facultative; FACU = facultative upland; UPL = upland. See text for statistics and definitions of C values and wetland indicators.

In terms of abundance, the greenway revealed a typical species-rank abundance curve in which a few species dominated, while a long tail of uncommon species contributed to species richness (Fig 4). The three most abundant species, which comprised ~50% of plant cover during our sampling window, were three non-native grasses: Bromus inermis L. (smooth brome), Bromus tectorum L. (cheatgrass), and Agropyron cristatum (L.) Gaertn. (crested wheatgrass).

Species relative abundances (percent cover) estimated using the ecological transects showing that only three species (the introduced grasses Bromus inermis L., Bromus tectorum L., and Agropyron cristatum (L.) Gaertn.) made up nearly half of the greenway’s percent cover during the 2018 sampling window.
Fig 4

Species relative abundances (percent cover) estimated using the ecological transects showing that only three species (the introduced grasses Bromus inermis L., Bromus tectorum L., and Agropyron cristatum (L.) Gaertn.) made up nearly half of the greenway’s percent cover during the 2018 sampling window.

Compare with the conceptual diagram of the species-rank abundance curve in Fig 1.

Discussion

Here we illustrate that sampling approaches from floristic botany and plant ecology capture complementary dimensions of biodiversity. As hypothesized, opportunistic sampling generated a markedly more robust (nearly four times larger) empirical estimate of the greenway’s community species pool than did the transects. Yet, the transects revealed that only three species of introduced grasses comprised 50% of plant cover at the time of sampling. Generating this critical indicator of the greenway’s ecological condition necessarily traded off with a fuller assessment of the species pool. This long-acknowledged trade-off is typically addressed by using rarefaction and accumulation curves and species richness estimators based on transect or plot data [46]. Our findings suggest that for sampling designs with limited areal coverage, such tools may not adequately capture the long tail of uncommon species that contribute disproportionately to species richness [13,47,48]. This issue was raised by Heilmann-Clausen and Læssøe [49], who clarified that species accumulation curves and species richness estimators address “how many species will be recorded if [a particular] sampling regime is followed in perpetuity or extended to cover all available habitat,” rather than telling the size of the species pool in the system. Related, Newmaster et al. (2005) illustrated that rarefaction and Chao estimates of common forest moss species reached an asymptote at small sample sizes (25 plots), while estimates of richness for rare species never leveled off. Thus, the very transect- or plot-based studies that often rely on estimator tools to determine the completeness of sampling likely produce underestimates of the true species pool [50,51]. Of course, many ecological questions can be rigorously answered without exhaustive documentation the species pool. Still, it is worthwhile to address these interpretive nuances when presenting results based on species estimators.

The question then becomes whether and when the transect-based sampling constraint, and any attendant underestimate of the species pool, affects how the ecology of the sampled area is interpreted (at least in terms of species richness). To assess this, we grouped species into ecologically meaningful categories to reduce complexity and uncover general patterns that are independent of species identity per se [52,53]. We found that based on functional groups and floristic quality (i.e., C-values), the ecological condition of the greenway appears similar between the two species pools. This comparable delineation of functional groups is especially desirable given that functional group identity and diversity are routinely used to gauge community response to disturbances such as fire and biological invasions [54,55]. In our system, functional traits such as longevity and woodiness (and attendant traits like rooting depth) likely shape soil and hydrologic conditions along the Canal banks. In terms of floristic quality, both sampling approaches captured the bell-shaped distribution of native species and similarly indicate that few conservative (C-value of 7 or more) species remain along this highly modified urban corridor. These findings corroborate previous work showing that floristic quality performs well when using transects or plots, because mean C-values are less dependent on the area sampled than is species richness [34].

While transects robustly captured patterns associated with functional groups and floristic quality, differences between the sampling approaches arose for Wetland Indicator Status and the presence of uncommon cultivated species. In particular, patchily distributed wetland areas and the obligate wetland species they harbor went largely undetected by the transects. This is a non-trivial miss in a Canal system where episodic drought and various management strategies strongly affect the hydrologic regime and thus persistence of sensitive wetland areas. Transects also failed to detect the nascent incursion of garden plants from adjacent private properties onto the Canal banks, potentially hindering Early Detection and Rapid Response (EDRR) management interventions [56,57]. These examples illustrate that without knowing the more complete species pool, the ability to thoroughly uncover sensitive ecological conditions along the greenway would be hampered. Still, if only species richness from the floristic inventory were used to assess the greenway, the high abundance of non-native grasses would go unreported. Without these abundance data, it would be impossible to assess costs associated with potential control or revegetation efforts, as well as to relate these dominant grasses to ecological processes of interest. For example, we are currently asking how the implementation of green stormwater infrastructure will affect vegetation along the Canal banks, and in turn, how in situ vegetation will affect stormwater infiltration, retention, and removal of pollutants. We are now positioned to take a two-pronged approach to this question by exploring how the most abundant species might shape stormwater dynamics, while also integrating information about the identity and location of uncommon species likely to be particularly responsive to hydrologic changes (e.g., obligate wetland species) and disturbance from infrastructure installations (e.g., establishing individuals of ruderal non-natives).

Broadening out from our greenway example, when should floristic botanists and ecologists develop on-the-ground collaborations to better describe contemporary biodiversity? We suggest that any question about the maintenance of community-level species diversity would benefit from a paired approach, as it is ultimately the interplay of local (competition, predation, microenvironmental variation) and regional (immigration and extinction) processes that shape biodiversity [5860]. To better integrate across spatial scales, it would be powerful to link quantitative data from bounded sampling to floristically based best estimates of the community species pool that functions as the backdrop for species immigration into embedded transects or plots. (While many terrestrial plant studies sample hierarchically across plot sizes to infer species pool sizes, they fall short of breaking free of plots to approach a comprehensive site- or habitat-level survey; [23]). Moreover, pairing floristic and ecological approaches addresses the call for increased metric complementarity in assessing biodiversity [61]. In particular, the historically heavy reliance upon species richness as a sole indicator of ecosystem health or biodiversity change has proven insufficient, as it fails to capture changes in other key phenomena such as species turnover and changes in species relative abundances [61,62]. Re-imagining floristic surveys of species richness as integral components of hypothesis-driven ecological work that uses other metrics can lead to new insights.

For example, to achieve metric complementarity in a restoration context, a floristic inventory could be used to assess the feasibility of passive restoration (which depends on the community species pool; [63,64]), while paired ecological sampling could quantify the degree of habitat degradation and monitor effectiveness of restoration efforts. Invasive species management would also benefit from a combined sampling approach, where initial arrivals of rare introduced species to an area are captured during unbounded floristic sampling bouts, while the spread and population biology of already established populations are monitored using transects or plots. Indeed, integrating researchers versed in alpha taxonomy is critical in invasion biology, as the resolution of taxonomically challenging groups, including those that hybridize, is critical to proper ecological interpretation [10]. A further boon of integrating floristic botany into settings where the species pool is of interest is excellent temporal sampling of early- to late-season bloomers, which can be missed using the “peak biomass” approach typical of ecological sampling [65]. Moreover, validation of transect- or plot-based data with vouchered specimens is tantamount to institutional knowledge that can be readily accessed by all researchers who carve out projects from a particular locale. The long-standing specimen and data curation practices used in natural history collections have achieved a level of standardization and data-sharing not yet realized in the ecological realm. However, improvements are being made in this arena, such as the application of the event-based Darwin Core data standard when publishing ecological data, as we have done herein [26,45].

To most informatively link plot- and site-scale diversity data, best practices would be to pair the replicable plot data with equally replicable floristic sampling of the community species pool [e.g., 22,62,63]. Collections-based floristic botany has not fully adopted standardized field sampling practices across individual collectors [66], which can limit the use of collections data for analyzing hypothesis-driven questions about community-level change in diversity [67]. The value of specimens and species lists generated using the opportunistic sampling approach could be increased by reporting the spatial extent of the surveyed area [68], as well as the intended goal of the collection event (e.g., an exhaustive inventory [implying species absence] versus targeted sampling based on an investigator’s taxa of interest). Such reporting practices can be achieved within the Ecological Metadata Language standard contained within a data package [69] and would provide additional context for downstream uses of aggregated data (e.g., species pool estimates derived from geo-referenced specimen databases) [70]. However, it must be considered that standardizing collections data (or at least reporting accurate areal coverage of survey sites) requires some degree of bounding that is both time-consuming and at odds with maximizing the number of species encountered [68]. Thus, it may only be worth bounding collections-based sampling when the specific question calls for it (e.g., quantifying immigration into experimental plots from the surrounding species pool).

Our experience is that physical collections are not often considered by ecological principal investigators as essential to their field protocols (despite substantial movement in this direction by, e.g., the National Ecological Observatory Network). This is partly because it is no small task to integrate the disparate training, project planning, data curation, and analyses implemented in floristic botany and ecology [7173]. Thus, cross-disciplinary partnerships across non-profit, governmental, and academic institutions are key. We suggest that ecologists reach out to their campus or regional herbaria, connect with curators and collections managers, and dedicate a line item in their research budgets for vouchered floristic (or faunistic) surveys of their study sites (see [74] for discussion of under-funding in collections-based research). We similarly suggest that curators and collections managers build relationships with ecology principal investigators and members of their labs, sharing their skills as integral assets to be included a priori into the proper design of biodiversity-focused ecological fieldwork (i.e., it is imperative to quash the “end-user mode” attitude that views botanists as simply providing identification services to those who use their keys and field guides [75,76]). For example, gaps in biodiversity data, including both species discovery and ecological monitoring, are high in tropical relative to temperate ecosystems [77,78], and would thus be best addressed by teams of floristic botanists and ecologists. As large-scale digitization of collections data has revealed, there are myriad, previously unimaged ways that natural history collections can inform ecological questions [79]. We believe our understanding of biodiversity 100 years from now can only benefit from thoughtful co-exploration of today’s ecosystems by floristic botanists and ecologists.

Acknowledgements

Dr. Melissa Islam was instrumental in providing intellectual and logistical guidance for this project. We thank a cadre of Denver Botanic Gardens volunteers for their many hours spent in the field helping make collections, as well as mounting and digitizing specimens for accession into the Kathryn Kalmbach Herbarium. We also thank several Gardens field seasonal staff including Audrey Spencer, Margo Yousse, Meghan McGill, and Katherine Fu, and intern Jodie Kirschner, for contributions to data collection, plant identifications, and map-making. Dr. Janet Wingate, Loraine Yeatts, and Nancy Downs provided expert guidance on difficult identifications of graminoids and cultivated plants.

References

BHopkins. The concept of minimal area. J Ecol. 1957;45: 441449.

SHoffmann, LSteiner, AHSchweiger, AChiarucci, CBeierkuhnlein. Optimizing sampling effort and information content of biodiversity surveys: A case study of alpine grassland. Ecol Inform. 2019;51: 112120. 10.1016/j.ecoinf.2019.03.003

DMStoms. Scale dependence of species richness maps. Prof Geogr. 1994;46: 346358. 10.1111/j.0033-0124.1994.00346.x

SACain. The species-area curve. Am Midl Nat. 1938;19: 573581.

SGNewmaster, RJBelland, AArsenault, DHVitt, TRStephens. The ones we left behind: Comparing plot sampling and floristic habitat sampling for estimating bryophyte diversity. Divers Distrib. 2005;11: 5772. 10.1111/j.1366-9516.2005.00123.x

JMChase, BJMcGill, DJMcGlinn, FMay, SABlowes, XXiao, et al Embracing scale-dependence to achieve a deeper understanding of biodiversity and its change across communities. Ecol Lett. 2018;21: 17371751. 10.1111/ele.13151

MCrawley, JHarral. Scale dependence in plant biodiversity. Science. (80-) 2001;291: 864868. 10.1126/science.291.5505.864

ORestif, DTSHayman, JRCPulliam, RKPlowright, DBGeorge, ADLuis, et al Model-guided fieldwork: Practical guidelines for multidisciplinary research on wildlife ecological and epidemiological dynamics. Ecol Lett. 2012;15: 10831094. 10.1111/j.1461-0248.2012.01836.x

AJMSchoot Uiterkamp, CVlek. Practice and outcomes of multidisciplinary research for environmental sustainability. J Soc Issues. 2007;63: 175197. 10.1111/j.1540-4560.2007.00502.x

10 

PPyšek, PEHulme, LAMeyerson, GFSmith, JSBoatwright, NRCrouch, et al Hitting the right target: Taxonomic challenges for, and of, plant invasions. AoB Plants. 2013;5: 125. 10.1093/aobpla/plt042

11 

JMajor. Plant ecology as a branch of botany. Ecology. 1958;39: 352363.

12 

JMHeberling, BLIsaac. Herbarium specimens as exaptations: New uses for old collections. Am J Bot. 2017;104: 963965. 10.3732/ajb.1700125

13 

AChiarucci, GBacaro, SMScheiner. Old and new challenges in using species diversity for assessing biodiversity. Philos Trans R Soc B Biol Sci. 2011;366: 24262437. 10.1098/rstb.2011.0065

14 

SHurlbert. Pseudoreplication and the design of ecological field experiments. Ecol Monogr. 1984;54: 187211.

15 

RHeffner, MButler, CReilly. Pseudoreplication Revisited. Ecology. 1996;77: 25582562.

16 

CGolodets, JKigel, YSapir, MSternberg. Quantitative vs qualitative vegetation sampling methods: A lesson from a grazing experiment in a Mediterranean grassland. Appl Veg Sci. 2013;16: 502508. 10.1111/avsc.12005

17 

KHeck, Gvan Belle, DSimberloff. Explicit calculation of the rarefaction diversity measurement and the determination of sufficient sample size. Ecology. 1975;56: 14591461.

18 

MFMoro, DJLde Sousa, LQMatias. Rarefaction, richness estimation and extrapolation methods in the evaluation of unseen plant diversity in aquatic ecosystems. Aquat Bot. 2014;117: 4855. 10.1016/j.aquabot.2014.04.006

19 

AChiarucci. To sample or not to sample? That is the question…for the vegetation scientist. Folia Geobot. 2007;42: 209216.

20 

NGotelli, RColwell. Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecol Lett. 2001; 379391.

21 

KJGaston. Common Ecology. Bioscience. 2011;61: 354362. 10.1525/bio.2011.61.5.4

22 

MESCáceres, RLücking, GRambold. Efficiency of sampling methods for accurate estimation of species richness of corticolous microlichens in the Atlantic rainforest of northeastern Brazil. Biodivers Conserv. 2008;17: 12851301. 10.1007/s10531-008-9342-3

23 

H V.Cornell, SPHarrison. What are species pools and when are they important? Annu Rev Ecol Evol Syst. 2014;45: 4567. 10.1146/annurev-ecolsys-120213-091759

24 

NatureServe. International ecological classification standard: Terrestrial ecological classifications. NatureServe Central Databases. Arlington, VA, U.S.A.; 2009.

25 

Alba C, Islam M. Field collecting protocol for vascular plants: Denver Botanic Gardens Kathryn Kalmbach Herbarium of Vascular Plants. 2019. 10.17504/protocols.io.4f4gtqw

26 

GBIF Occurrence at GBIF.org. 2020.

27 

JAckerfield. Flora of Colorado. First Edit BLipscomb, editor. Fort Worth: BRIT Press; 2015.

28 

Wingate J. Illustrated Keys to the Grasses of Colorado. Denver; 1994.

29 

Wingate J. Sedges of Colorado. Denver; 2017.

30 

JHerrick, JVan Zee, SMcCord, ECourtright, JKarl, LBurkett. Monitoring manual for grassland, shrubland, and savanna ecosystems. Second Edi Las Cruces: USDA-ARS Jornada Experimental Range; 2017.

31 

Hufft R, Alba C, Sahud A. Vegetation monitoring protocol for measuring and collecting ecological data. 2019;version 2. 10.17504/protocols.io.4jmguk6

32 

CAlba, RLevy, BHufft. Ecological data from: Combining botanical collections and ecological data to better describe plant community diversity. Dryad. 2020 10.5061/dryad.15dv41nw6

33 

MZobel, Evan der Maarel, CDupré. Species pool: The concept, its determination and significance for community restoration. Appl Veg Sci. 1998;1: 5566. 10.2307/1479085

34 

GSpyreas. Floristic Quality Assessment: a critique, a defense, and a primer. Ecosphere. 2019;10 10.1002/ecs2.2825

35 

US Army Corps of Engineers. National Wetland Plant List. 2018. http://wetland-plants.usace.army.mil/nwpl_static/v34/species/species.html?DET=001100

36 

LMVerheijen, RAerts, GBönisch, JKattge, PMVan Bodegom. Variation in trait trade-offs allows differentiation among predefined plant functional types: Implications for predictive ecology. New Phytol. 2016;209: 563575. 10.1111/nph.13623

37 

JHCCornelissen, SLavorel, EGarnier, SDíaz, NBuchmann, DEGurvich, et al A handbook of protocols for standardised and easy measurement of plant functional traits worldwide. Aust J Bot. 2003;51: 335380. 10.1071/BT02124

38 

NPérez-Harguindeguy, SDíaz, EGarnier, SLavorel, HPoorter, PJaureguiberry, et al New handbook for standardised measurement of plant functional traits worldwide. Aust J Bot. 2013;61: 167234. 10.1071/BT12225

39 

JRocchio, DAnderson, DBuckner, KCarsey, DClark, JColes, et al Floristic quality assessment indices for Colorado plant communities. Color Nat Herit Program, Color State Univ. 2007;1: 1245.

40 

RLichvar, DBanks, WKirchner, NMelvin. The National Wetland Plant List: 2016 wetland ratings. Phytoneuron. 2016;30: 117.

41 

RKColwell, AChao, NJGotelli, SYLin, CXMao, RLChazdon, et al Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. J Plant Ecol. 2012;5: 321. 10.1093/jpe/rtr044

42 

Colwell R. EstimateS: Statistical estimation of species richness and shared species from samples. 2016;Version 9. http://viceroy.eeb.uconn.edu/estimates/index.html

43 

MLMchugh. The Chi-square test of independence. Biochem Medica. 2013;23: 1439. 10.11613/bm.2013.018

44 

DSharpe. Chi-square test is statistically significant: Now what? Pract Assessment, Res Eval. 2015;20: 110.

45 

GBIF Occurrrence at GBIF.org. 2020.

46 

AChao, RKColwell, CLin, NJGotelli. Sufficient sampling for asymptotic minimum species richness estimators. Ecology. 2009;90: 11251133. 10.1890/07-2147.1

47 

WUlrich, BKusumoto, SFattorini, YKubota. Factors influencing the precision of species richness estimation in Japanese vascular plants. Divers Distrib. 2020;26: 769778. 10.1111/ddi.13049

48 

AChiarucci, NJEnright, GLWPerry, BPMiller, BBLamont. Performance of nonparametric species richness estimators in a high diversity plant community. Divers Distrib. 2003;9: 283295. 10.1046/j.1472-4642.2003.00027.x

49 

JHeilmann-Clausen, TLæssøe. On species richness estimates, climate change and host shifts in wood-inhabiting fungi. Fungal Ecol. 2012;5: 641646. 10.1016/j.funeco.2011.10.003

50 

MPalmer. The estimation of species richness by extrapolation. Ecology. 1990;71: 11951198.

51 

AChao, C-HChiu. Nonparametric estimation and comparison of species richness eLS. Chichester: John Wiley & Sons, Ltd; 2016 pp. 111.

52 

MWCadotte, KCarscadden, NMirotchnick. Beyond species: Functional diversity and the maintenance of ecological processes and services. J Appl Ecol. 2011;48: 10791087. 10.1111/j.1365-2664.2011.02048.x

53 

DHooper, MSolan, SSymstad, SDiaz, MGessner, NBuchmann, et al Species diversity, functional diversity and ecosystem functioning. Biodivers Ecosyst Funct Synth Perspect. 2002;17: 195208. 10.1126/science.312.5775.846a

54 

ASymstad. A test of the effects of functional group richness and composition on grassland invasibility. Ecology. 2000;81: 99109.

55 

NJEnright, JBFontaine, BBLamont, BPMiller, VCWestcott. Resistance and resilience to changing climate and fire regime depend on plant functional traits. J Ecol. 2014;102: 15721581. 10.1111/1365-2745.12306

56 

RWestbrooks. New approaches for early detection and rapid response to invasive plants in the United States. Weed Te. 2004;18: 14681471.

57 

JKReaser, SWBurgiel, JKirkey, KABrantley, SDVeatch, JBurgos-Rodríguez. The early detection of and rapid response (EDRR) to invasive species: a conceptual framework and federal capacities assessment. Biol Invasions. 2020;22: 119. 10.1007/s10530-019-02156-w

58 

MLoreau, NMouquet. Immigration and the maintenance of local species diversity. Am Nat. 1999;154: 427440. 10.1086/303252

59 

SLCollins, SMGlenn, JMBriggs. Effect of local and regional processes on plant species richness in tallgrass prairie. Oikos. 2002;99: 571579.

60 

MJCaley, DSchluter. The relationship between local and regional diversity. Ecology. 1997;78: 7080.

61 

HHillebrand, BBlasius, ETBorer, JMChase, JADowning, BKEriksson, et al Biodiversity change is uncoupled from species richness trends: Consequences for conservation and monitoring. J Appl Ecol. 2018;55: 169184. 10.1111/1365-2664.12959

62 

SLLHill, MHarfoot, APurvis, DWPurves, BCollen, TNewbold, et al Reconciling biodiversity indicators to guide understanding and action. Conserv Lett. 2016;9: 405412. 10.1111/conl.12291

63 

KDHoll, TMAide. When and where to actively restore ecosystems? For Ecol Manage. 2011;261: 15581563. 10.1016/j.foreco.2010.07.004

64 

ASundermann, SStoll, PHaase. River restoration success depends on the species pool of the immediate surroundings. Ecol Appl. 2011;21: 19621971. 10.1890/10-0607.1

65 

JDoležal, VLanta, OMudrák, JLepš. Seasonality promotes grassland diversity: Interactions with mowing, fertilization and removal of dominant species. J Ecol. 2019;107: 203215. 10.1111/1365-2745.13007

66 

PSastre, JMLobo. Taxonomist survey biases and the unveiling of biodiversity patterns. Biol Conserv. 2009;142: 462467. 10.1016/j.biocon.2008.11.002

67 

AChiarucci, IBonini. Quantitative floristics as a tool for the assessment of plant diversity in Tuscan forests. For Ecol Manage. 2005;212: 160170. 10.1016/j.foreco.2005.03.041

68 

MWPalmer, PGEarls, BWHoagland, PSWhite, TWohlgemuth. Quantitative tools for perfecting species lists. Environmetrics. 2002;13: 121137. 10.1002/env.516

69 

Jones M, O’Brien M, Mecum B, Boettiger C, Schildhauer M, Maier M, et al. Ecological Metadata Language version 2.2.0. In: KNB Data Repository. 2019.

70 

HBruelheide, BJiménez-Alfaro, UJandt, FMSabatini. Deriving site-specific species pools from large databases. Ecography (Cop). 2020;43:12151228. 10.1111/ecog.05172

71 

MD’Antraccoli, GBacaro, ETordoni, GBedini, LPeruzzi. More species, less effort: Designing and comparing sampling strategies to draft optimised floristic inventories. Perspect Plant Ecol Evol Syst. 2020;45: 125547 10.1016/j.ppees.2020.125547

72 

JCTouchon, MWMcCoy. The mismatch between current statistical practice and doctoral training in ecology. Ecosphere. 2016;7: 111. 10.1002/ecs2.1394

73 

RBritz, AHundsdorfer, UFritz. Funding, training, permits—the three big challenges of taxonomy. Megataxa. 2020;1: 4952. 10.11646/megataxa.1.1.10

74 

SSaunders. Taxonomy at a crossroads: Communicating value, building capability, and seizing opportunities for the future. Megataxa. 2020;1: 6366. 10.11646/megataxa.1.1.13

75 

RSluys. The unappreciated, fundamentally analytical nature of taxonomy and the implications for the inventory of biodiversity. Biodivers Conserv. 2013;22: 10951105. 10.1007/s10531-013-0472-x

76 

MCEbach, AGValdecasas, QDWheeler. Impediments to taxonomy and users of taxonomy: Accessibility and impact evaluation. Cladisitics. 2011;27: 550557. 10.1111/j.1096-0031.2011.00348.x

77 

OLPhillips, RVasquez Marti­nez, PNunez Vargas, ALMonteagudo, MChuspe Zans, WGaliano Sanchez, et al Efficient plot-based floristic assessment of tropical forests. J Trop Ecol. 2003;19: 629645. 10.1017/S0266467403006035

78 

BCollen, MRam, TZamin, LMcRae. The tropical biodiversity data gap: Addressing disparity in global monitoring. Trop Conserv Sci. 2008;1: 7588. 10.1177/194008290800100202

79 

SAJames, PSSoltis, LBelbin, ADChapman, GNelson, DLPaul, et al Herbarium data: Global biodiversity and societal botanical needs for novel research: Global. Appl Plant Sci. 2018;6: 18. 10.1002/aps3.1024