Niche Models Di ﬀ erentiate Potential Impacts of Two Aquatic Invasive Plant Species on Native Macrophytes

: Potamogeton crispus (curlyleaf pondweed) and Myriophyllum spicatum (Eurasian watermilfoil) are widely thought to competitively displace native macrophytes in North America. However, their perceived competitive superiority has not been comprehensively evaluated. Coexistence theory suggests that invader displacement of native species through competitive exclusion is most likely where high niche overlap results in competition for limiting resources. Thus, evaluation of niche similarity can serve as a starting point for predicting the likelihood of invaders having direct competitive impacts on resident species. Across two environmental gradients structuring macrophyte communities—water depth and light availability—both P. crispus and M. spicatum are thought to occupy broad niches. For a third dimension, phenology, the annual growth cycle of M. spicatum is typical of other species, whereas the winter-ephemeral phenology of P. crispus may impart greater niche di ﬀ erentiation and thus lower risk of native species being competitively excluded. Using an unprecedented dataset comprising 3404 plant surveys from Minnesota collected using a common protocol, we modeled niches of 34 species using a probabilistic niche framework. Across each niche dimension, P. crispus had lower overlap with native species than did M. spicatum ; this was driven in particular by its distinct phenology. These results suggest that patterns of dominance seen in P. crispus and M. spicatum have likely arisen through di ﬀ erent mechanisms, and that direct competition with native species is less likely for P. crispus than M. spicatum. This research highlights the utility of ﬁne-scale, abundance-based niche models for predicting invader impacts.


Introduction
Invasive species are among the greatest threats to biodiversity worldwide [1,2]. Across taxonomic groups and trophic levels, invaders have consistent negative effects on native communities' species richness, evenness, and other measures of diversity [3]. Despite these clear overall patterns, important gaps remain in our ability to predict how particular invaders will interact with particular native species. The strength and directionality of individual invader-native interactions, and the attendant impacts of invasions on communities as a whole, can be difficult to assess and are often poorly understood [1,4]. Variation in impacts of invasive species arises through multiple mechanisms by which invasive species interact with native species and alter invaded communities [4][5][6]. These mechanisms range from altering disturbance regimes [7] to better tolerating diseases [8] to directly competing for limiting resources [9]. Direct competition for limiting resources can cause competitive displacement of native species-one of multiple mechanisms that can be detrimental to native communities and give rise to a pattern of invader dominance [10,11]. occurrences at the local scale (i.e., within-lake distributions). To date, these data have been disparate in their storage and analysis-we retrieved and united these disconnected data to develop an extensive dataset covering a wide geographic range and time period. We incorporated these data into the probabilistic niche modeling framework, which allows robust evaluation of niche differentiation among species. This approach generated fine-scale niche models for Minnesota macrophytes incorporating both species occupancy along each niche axis and abundance within occupied segments of the niche (Figure 1). The use of a novel dataset that unites disparate but consistently collected local-scale plant data enables niche models to be built based on neighborhood-level species occurrences and fine-scale environmental variation, i.e., within-lake differences in habitat conditions. The probabilistic nature of this approach is made possible by inclusion of abundance data (how often was a species observed in particular conditions relative to all other conditions) and results in niche models that predict likelihood of occurrence along each modeled niche axis, as opposed to only defining niche boundaries [55]. Finally, the probabilistic niche modeling framework enabled attribution of niche differences to either differences in the niche space occupied ("non-shared niche differences") or differences in the abundances of two species in co-occupied niche space ("shared niche differences") ( Figure 1). This extensive dataset and use of probabilistic niche models allowed (1) development of within-lake habitat models (as opposed to coarse regional models based on climate or other large-scale variables), (2) incorporation of abundance along niche axes to account for differential distribution within fundamental niches, and (3) discrimination of niche differentiation as both differences in niche boundaries and differences in relative abundances within shared niche space.
Diversity 2020, 12, x FOR PEER REVIEW 3 of 18 managers across Minnesota, USA, all implementing a unified survey methodology to collect data on plant occurrences at the local scale (i.e., within-lake distributions). To date, these data have been disparate in their storage and analysis-we retrieved and united these disconnected data to develop an extensive dataset covering a wide geographic range and time period. We incorporated these data into the probabilistic niche modeling framework, which allows robust evaluation of niche differentiation among species. This approach generated fine-scale niche models for Minnesota macrophytes incorporating both species occupancy along each niche axis and abundance within occupied segments of the niche (Figure 1). The use of a novel dataset that unites disparate but consistently collected local-scale plant data enables niche models to be built based on neighborhoodlevel species occurrences and fine-scale environmental variation, i.e., within-lake differences in habitat conditions. The probabilistic nature of this approach is made possible by inclusion of abundance data (how often was a species observed in particular conditions relative to all other conditions) and results in niche models that predict likelihood of occurrence along each modeled niche axis, as opposed to only defining niche boundaries [55]. Finally, the probabilistic niche modeling framework enabled attribution of niche differences to either differences in the niche space occupied ("non-shared niche differences") or differences in the abundances of two species in cooccupied niche space ("shared niche differences") ( Figure 1). This extensive dataset and use of probabilistic niche models allowed 1) development of within-lake habitat models (as opposed to coarse regional models based on climate or other large-scale variables), 2) incorporation of abundance along niche axes to account for differential distribution within fundamental niches, and 3) discrimination of niche differentiation as both differences in niche boundaries and differences in relative abundances within shared niche space.

Figure 1.
Graphical depiction of niches modeled using probability density functions (PDFs) [53] for two hypothetical species along one niche axis. Each PDF is built from occurrence (histogram of occurrences shown for Sp 1 only) and associated environmental data and PDFs integrate to 1. Dissimilarity is estimated as follows: total dissimilarity = 1 -shared area (a), non-shared niche dissimilarity is the ratio of dissimilar PDF area in niche space occupied by only one of the two species (d, e) to all dissimilar PDF area (d + b + c + e), and shared dissimilarity is the ratio of dissimilar PDF area in niche space occupied by both species (b, c) to all dissimilar PDF area (d + b + c + e).
We used the probabilistic niche modeling framework to assess the potential for competition between the focal invasive species and native macrophyte species. We developed niche models for P. crispus, M. spicatum, and 32 native macrophyte species that were sufficiently represented in our dataset for modeling purposes. For each invader, we calculated its niche overlap with each native species across all three niche dimensions as a whole and for each individual dimension. We further distinguished niche differences between species that arose from differences in shared vs. non-shared Figure 1. Graphical depiction of niches modeled using probability density functions (PDFs) [53] for two hypothetical species along one niche axis. Each PDF is built from occurrence (histogram of occurrences shown for Sp 1 only) and associated environmental data and PDFs integrate to 1. Dissimilarity is estimated as follows: total dissimilarity = 1 − shared area (a), non-shared niche dissimilarity is the ratio of dissimilar PDF area in niche space occupied by only one of the two species (d, e) to all dissimilar PDF area (d + b + c + e), and shared dissimilarity is the ratio of dissimilar PDF area in niche space occupied by both species (b, c) to all dissimilar PDF area (d + b + c + e).
We used the probabilistic niche modeling framework to assess the potential for competition between the focal invasive species and native macrophyte species. We developed niche models for P. crispus, M. spicatum, and 32 native macrophyte species that were sufficiently represented in our dataset for modeling purposes. For each invader, we calculated its niche overlap with each native species across all three niche dimensions as a whole and for each individual dimension. We further distinguished niche differences between species that arose from differences in shared vs. non-shared niche space. Using these results, we tested the hypothesis that M. spicatum would have greater niche overlap with native species than P. crispus, suggesting greater potential for direct competitive displacement by M. spicatum.

Macrophyte Occurrence Data
In Minnesota, many organizations (state agencies, local units of government, consultants, and others) conduct point-intercept surveys of aquatic plant communities in lakes. These surveys are used for a wide variety of purposes, including baseline monitoring [56], long-term ecological trend monitoring (SLICE: https://www.dnr.state.mn.us/fisheries/slice/index.html), university research projects [57,58], and monitoring of management outcomes [59]. Recent work has illustrated the research opportunities these datasets provide, as they are cohesive surveys conducted using a common method over long time scales and wide geographic ranges [59][60][61][62]. We compiled plant occurrence and depth data from point-intercept surveys performed by agencies, researchers, and lake managers in the period 2000-2018. These data are representative of the geographic distribution of lakes in Minnesota ( Figure 2). Diversity 2020, 12, x FOR PEER REVIEW 4 of 18 niche space. Using these results, we tested the hypothesis that M. spicatum would have greater niche overlap with native species than P. crispus, suggesting greater potential for direct competitive displacement by M. spicatum.

Macrophyte Occurrence Data
In Minnesota, many organizations (state agencies, local units of government, consultants, and others) conduct point-intercept surveys of aquatic plant communities in lakes. These surveys are used for a wide variety of purposes, including baseline monitoring [56], long-term ecological trend monitoring (SLICE: https://www.dnr.state.mn.us/fisheries/slice/ index.html), university research projects [57,58], and monitoring of management outcomes [59]. Recent work has illustrated the research opportunities these datasets provide, as they are cohesive surveys conducted using a common method over long time scales and wide geographic ranges [59][60][61][62]. We compiled plant occurrence and depth data from point-intercept surveys performed by agencies, researchers, and lake managers in the period 2000-2018. These data are representative of the geographic distribution of lakes in Minnesota ( Figure 2).  All data preparation, analysis, and visualization of results were conducted in R version 3.6.1 [64]. All code used for analyses and visualizations are included in Supplementary Material (Supplementary Code). Data used for analysis were collected via the point-intercept method, in which surveyors navigate to a predetermined grid of points in a lake and use a metal rake to sample macrophytes growing in that location. Use of point-intercept data ensured unbiased sample selection (point locations are pre-determined prior to surveys) and consistent methodology across surveys [65]. We analyzed data at the point scale, i.e., each record comprised one sampling point with associated survey information (date, surveyor, lake identifier), and the water depth and plant species recorded at the point. The average number of points sampled in point-intercept surveys was 104 (± 108, 1 SD; distribution shown in Supplementary Code). For our analysis, we treated each plant observation at each survey point as a data point (mean taxa observed per sample point 2.05 ± 1.62, 1 SD; Supplementary Code) and compiled associated data for each of these observations. Plant occurrences from each survey point were translated to scientific names from codes or common names using keys provided by the data contributors and referencing common names in regional identification guides [66][67][68][69]. We standardized species names based on the Taxonomic Name Resolution Service [70] using the "taxize" package in R 3.6.1 [64,71]. In total, we compiled 3404 surveys comprising 2514 lake-years of data (each lake × year combination comprises a lake-year) from 1526 Minnesota lakes surveyed at least once between 2000 and 2018, representing 353,148 total survey points, with a combined 564,038 macrophyte observations across 218 taxa (163 identified to species, 56 to higher taxa when surveyors were unable to resolve to species).

Environmental Data
Water depth data were collected by surveyors at each sampling point. We omitted sample points with no depth information, a depth of 0, or depths > 5 standard deviations above the mean (depth > 7 m, 0.4% of all occurrences). If units were not specified for depth (i.e., whether in feet or meters), we plotted the distribution of depths and plant occurrences observed in the survey and compared these to lake bathymetry and reasonable colonization depths of plants. For example, if a survey with unlabeled depths occurred in a lake with a maximum depth of 3 m, and unlabeled depths ranged from 0 to 8.5, we concluded that depth must have been recorded in feet. Secchi depth data for each lake-year were collected from several sources and cleaned to ensure consistency. For most surveys, we used an empirical Secchi dataset compiled to develop predictive models of lake clarity for Minnesota (Kelsey Vitense, University of Minnesota, unpublished data). If a lake-year in our dataset did not have a corresponding value in this dataset, we used Secchi data collected by the surveyor if available. We only used Secchi measurements taken in summer (July, August, or September) to avoid measurements sensitive to large seasonal changes in water clarity [72]. Finally, summer Secchi data for a three-year window (the survey year, the year prior, and the year after) were used to calculate a mean Secchi value for each survey (lake-year) to minimize the effects of sampling error and variability in Secchi readings [73].
Three-year average Secchi depths for each survey and point-specific water depth were used to estimate light availability for each survey point [73][74][75][76], calculated as the proportion of surface irradiance remaining at the substrate (hereafter, light availability). Light at Secchi depth was assumed to be 10% of surface irradiance, which enabled us to use Secchi depth measurements to estimate a lake-level light decay constant, K d [75]: We then used the Lambert-Beer relationship to calculate the proportion of surface irradiance (light availability) remaining at the substrate for each plant observation's associated depth to substrate [76] using the lake-wide K d and assuming water clarity to be constant throughout each waterbody: Secchi data were sufficient in coverage to calculate light availability for 82.9% of plant observations in our dataset. We examined phenology using growing degree days (GDD), a measure of the heat energy received by a plant over a given time period (e.g., within a year) [77]. Growing degree days are a more robust way to assess phenology compared to time of year or day of the year, which do not account for geographic location or interannual variability in temperature [78,79]. Daily minimum and maximum temperatures were required to calculate GDD. For each lake, we identified the nearest weather station with comprehensive daily weather data using the "get_weather" function in the chillR package in R [80]. We used the same function to download daily weather data for each weather station from the Integrated Surface Dataset [81]. A total of 16 weather stations were used with a mean distance of 64 km (± 31 km, 1 SD) from associated lakes. Some stations were located outside of Minnesota, but these political boundaries are not associated with abrupt disjuncts in climate or environmental conditions, as for example if the borders were associated with abrupt changes in elevation ( Figure 2). Using associated station data, we calculated GDD as, GDD = (maximum daily temperature + minimum daily temperature)/2 − base temperature, (3) where base temperature is temperature below which no plant growth is expected. We used 4 • C as the base temperature, following Boissezon et al. [82], and 1 March as the "biofix date", i.e., the date on which GDD began accumulating each year. For some daily weather records, minimum and/or maximum temperatures were missing (<0.2% of all records). When this occurred, we substituted the previous day's minimum and/or maximum temperature. We used the cumulative yearly GDD for each vegetation survey date as the GDD value associated with each survey point.
Following the collection, cleaning, and merging of all datasets, we retained a total of 462,118 plant occurrences from 292,824 sample points for which we had complete depth, light, and phenology data.

Data Analysis
We constructed probabilistic niche models by fitting probability density functions (PDFs) to the occurrence of each species and axis of interest using a kernel density estimation (KDE) method [83][84][85]. These PDFs represent the realized niche occupied by that species based on the three niche dimensions we evaluated (depth, light, and phenology). The result is a model wherein any location within the realized niche of a species is represented by a probability of its occurrence at that location in three-dimensional niche space. We then compared niches between species to reveal where species are likely to compete for the same resources, i.e., where the species have overlapping niche requirements. We conducted niche modeling and dissimilarity evaluations among species using the "TPD" package in R [83].
To characterize the availability of potential niche space represented in our dataset, we first modeled the probabilistic niche of all survey points, i.e., generated a representation of the total "niche" occupied by all survey points. Probability distributions of sampling effort indicated sampling effort was non-uniform (Supplementary Code). Specifically, sampling was biased toward shallower sites and those with lower light availability, and survey efforts peaked mid-summer at 1500 GDD. Nonetheless, sampling effort was sufficient to provide good coverage of the multi-dimensional niche space across all three dimensions (Figure 3, Figures S1-S3). In addition, correlations among niche axes indicated no issues with collinearity for use in niche models (Light-GDD r = −0.01, Light-Depth r = −0.41, GDD-Depth r = −0.05; visualizations presented in Supplementary Code).
We standardized niche data (centered data on the variable mean and scaled it to one standard deviation) to ensure equal weights among variables measured on different scales [54,85]. To ensure sufficient fit of PDFs, we limited our analysis to species with a minimum number of occurrences equal to e n niche dimensions , or 1000 occurrences [85,86]. This restricted our dataset to 371,053 occurrences of 34 species used in the analysis.
For each species × niche dimension, we fit a PDF to the abundance data across that dimension using KDE, such that the final probability function for each dimension integrated to 1 (i.e., the area under the curve = 1; Figure 1). Differences between the probabilistic niche models produced for P. crispus and M. spicatum for each niche dimension were evaluated using a Kolmogorov-Smirnov test. We also constructed PDFs for each species across all three niche dimensions simultaneously, then estimated the dissimilarity of species' niche models using the proportion of overlap in the volume (three niche dimensions) of the two species' niche volumes [53]. In three-dimensional niche space, the total dissimilarity between two species is equivalent to 1 − the niche model PDF volume shared by the two species. We standardized niche data (centered data on the variable mean and scaled it to one standard deviation) to ensure equal weights among variables measured on different scales [54,85]. To ensure sufficient fit of PDFs, we limited our analysis to species with a minimum number of occurrences equal  Dissimilarity among species was further decomposed into two elements: the dissimilarity of abundances within niche space shared between the two species and the dissimilarity measured as non-shared niche space occupied by only one of the two species (Figure 1) [53,83]. In this way, we were able to parse the relative importance of two different ways that niche space can be partitioned between species. Finally, we calculated the dissimilarity of each species from all other species (the total community) as the mean of all pairwise PDF dissimilarity values for that species. A two-proportion z-test was used to evaluate the hypothesis that the dissimilarity of P. crispus was more strongly driven by contributions of non-shared dissimilarity than that of M. spicatum. Species dissimilarity was further subdivided as the dissimilarity contributed individually by each of the three niche dimensions.

Results
For the 34 species evaluated, we produced PDFs for the three niche dimensions considered (Figure 3, Figures S1-S3). These distributions were constructed for multi-dimensional niche space ( Figure 3) and each individual niche dimension (Figure 4). In both cases, there was substantial phenological niche separation for P. crispus, for which most occurrences were under low-GDD values (~750 GDD). Potamogeton crispus also had peaks in its abundance under low-light conditions and water depth of~2.5 m (Figure 3). In comparison, M. spicatum showed broad phenology with two peaks, one near 900 GDD and the other at approximately 2000 GDD (Figure 4). Relative to P. crispus, there was lower occupancy under low-light conditions and broader occupancy across depths (Figures 3 and 4).
On average, P. crispus was the most dissimilar to all species of any species evaluated, whereas M. spicatum was much more similar to other plant species ( Figure 5). The dissimilarity of P. crispus was driven by a lower use of niche space shared with native species than M. spicatum. This difference in niche use was reflected in the strong relative contribution of non-shared niche space ( Figure 5), indicating that a difference in absolute niche boundaries, rather than relative abundance within shared niche space, drove the dissimilarity of P. crispus. Whereas the total dissimilarity of P. crispus to native species was strongly influenced by non-shared niche space, the dissimilarity of M. spicatum was primarily driven by differences in shared niche space. In comparisons with 25 of 32 native species, non-shared niche dissimilarity was greater for P. crispus than it was for M. spicatum. In addition, when compared to each native species, the relative contribution of non-shared niche dissimilarity for P. crispus exceeded that of shared niche dissimilarity 72% of the time (24 of 33 cases) compared to 0% for M. spicatum (0 of 33 cases); this difference was significant (two proportion z-test with continuity correction, p < 0.001).
The distribution of P. crispus and M. spicatum differed across all three niche dimensions evaluated. For P. crispus, niche differentiation was greatest in phenology (Figures 4 and 6); P. crispus exhibited earlier annual growth (mean = 747 GDD) than M. spicatum (mean = 1351 GDD), and native species as a whole (mean = 1254 GDD). Potamogeton crispus also differed more from other species than did M. spicatum in terms of depth and light (Figures 4 and 6). Overall, P. crispus exhibited earlier phenology and occupied locations that were deeper and had lower light availability than M. spicatum (Figures 3 and 4, Kolmogorov-Smirnov tests: p < 0.001 for all three dimensions, D values = 0.41, 0.16, and 0.31 for GDD, depth, and light availability, respectively).   The distribution of P. crispus and M. spicatum differed across all three niche dimensions evaluated. For P. crispus, niche differentiation was greatest in phenology (Figures 4 and 6); P. crispus exhibited earlier annual growth (mean = 747 GDD) than M. spicatum (mean = 1351 GDD), and native species as a whole (mean = 1254 GDD). Potamogeton crispus also differed more from other species than did M. spicatum in terms of depth and light (Figures 4 and 6). Overall, P. crispus exhibited earlier phenology and occupied locations that were deeper and had lower light availability than M. spicatum (Figures 3 and 4, Kolmogorov-Smirnov tests: p < 0.001 for all three dimensions, D values = 0.41, 0.16, and 0.31 for GDD, depth, and light availability, respectively).

Discussion
Niche differentiation of P. crispus was greater than that of M. spicatum, suggesting that P. crispus may be less likely to directly displace native species through competition for limiting resourcesthus, previously observed detrimental impacts on native assemblages for each of these species have likely arisen through differing mechanisms. Niche models exhibited three main patterns in support of this conclusion: 1) a lower overall similarity of the P. crispus niche to those of native macrophyte species, 2) a greater contribution of non-shared niche space (vs. shared niche space) to the overall dissimilarity of P. crispus than was observed for M. spicatum, and 3) less similarity between P. crispus

Discussion
Niche differentiation of P. crispus was greater than that of M. spicatum, suggesting that P. crispus may be less likely to directly displace native species through competition for limiting resources-thus, previously observed detrimental impacts on native assemblages for each of these species have likely arisen through differing mechanisms. Niche models exhibited three main patterns in support of this conclusion: (1) a lower overall similarity of the P. crispus niche to those of native macrophyte species, (2) a greater contribution of non-shared niche space (vs. shared niche space) to the overall dissimilarity of P. crispus than was observed for M. spicatum, and (3) less similarity between P. crispus and other macrophytes within each individual niche axis, with the phenological component making the greatest contribution to niche dissimilarity.
The expectation that invasive species' competitive superiority will lead to the exclusion of co-occurring native species and give rise to invader dominance is often the basis of concerns regarding non-native plant invasions [7, 11,87]. By comparing habitat-use strategies of invasive species to resident species of invaded communities, we assessed the potential risk of displacement of native species through competitive interactions. This study provides evidence for large differences in niche overlap with native species for two well-established invasive species that are the subject of intensive management efforts, and indicates that their impacts on native plant communities are likely to differ in magnitude with respect to direct competition and/or arise through different mechanisms [9,88].
Competition for limiting resources can lead to loss of native species following the addition of a competitively superior invader to a community [16], but thoroughly evaluating the extent to which resource competition drives impacts of invasive species to native communities requires experimental manipulations [11,89,90] rather than the observational data used in this study. Further, there are numerous indirect, non-competitive mechanisms by which invasive species can affect native organisms (e.g., [91]). For example, community impacts of P. crispus could arise through increased turbidity of lakes driven by its mid-summer senescence and release of biomass-contained nutrients that stimulate algal productivity [9,92,93]. With the niche-based approach we used, it would not be possible to predict such-potentially substantial-invader impacts. Our approach using observational data cannot address alternative, non-competitive mechanisms of interactions between invasive and native species, rather it provides an estimation of whether direct competition for limited niche space per se is likely to drive displacement of particular native species by particular invasive species. Thus, our findings should not be interpreted as meaning that P. crispus is not detrimental to native macrophytes or is necessarily less detrimental than M. spicatum.
Nonetheless, the distinctiveness of the niche occupied by P. crispus suggests that its potential for direct competitive displacement of native species is lower than that of M. spicatum [14]. Despite this, it is critical to note that both species have succeeded in establishing in North America and both can exhibit dominance of lake plant communities [48,94]. Our results suggest that P. crispus may achieve high abundance through occupancy of a relatively vacant niche [12,86], whereas M. spicatum reaches high abundances under conditions favored by common native species, which may indicate competitive superiority [87] and greater likelihood of driving native species declines via competitive exclusion.
We also found that the dissimilarity of P. crispus to the native community was driven much more by occupancy of non-shared (novel) niche space than in the case of M. spicatum, for which there was greater occupancy of niche space well used by native species. These differences are important to consider in the context of the influence these invaders may have on native species [88,95] and suggest that the competitive threat to native species posed by P. crispus may be lower. However, because the phenological niche we evaluated (GDD) has a strong association with annual growth patterns, P. crispus may be dominating a phenological niche that is important for critical early life stages of native species that are underrepresented through rake-based surveys. Specifically, rake-based methods under sample low-growing or recently sprouted plants [96]. Yet, it is precisely during short-statured, early life stages when macrophyte germinants may be most susceptible to shading by P. crispus individuals that are further into their annual life cycle. Such a phenomenon would be missed by the sampling methods and phenological measures applied in this study. However, in the case of P. crispus, we found that the niche occupied by its early phenology was also highly dissimilar from native species in terms of depth. Because of this, we suspect that P. crispus is capitalizing on seasonal water clarity variations to grow at deeper locations than would be possible for plants with later phenology that encounter lower light availability [72]. In other words, P. crispus may be exploiting a deeper, earlier niche to grow in locations that have insufficient light availability by midsummer (due to high algal biomass in productive, nutrient-rich lakes in temperate regions) [97]. This suggests that P. crispus may have a minimal competitive effect on an early or critical phase of native growth because it is not only temporally separated from native species but also spatially separated.
Interestingly, for both P. crispus and M. spicatum, phenological differentiation was the greatest contributor to niche dissimilarity, though it was exhibited in different ways. In P. crispus, the strong effect of an early-season growth strategy was evident [51], whereas the phenological niche we observed for M. spicatum was highly bimodal. This bimodality is consistent with previous descriptions of M. spicatum phenology [98,99]. However, other studies have described high biomass throughout the growing season, accentuated by two summer peaks [46,[98][99][100], whereas our analysis identified a larger apparent gap in its presence in surveys, despite representative sampling effort across GDD. This pattern of bimodality has previously been attributed to sloughing of plant parts between two flowering peaks [100], but the pattern observed in our results shows a lack of presence of M. spicatum in mid-season surveys (our dataset did not include biomass at occurrences). Future work could follow multiple populations through the growing season to evaluate whether observed bimodality arises from of a two-peaked phenology, a variation in the phenological pattern among lakes (populations tending to have either early-or late-season phenology), or a product of both.
A variety of mechanisms may allow each species to compensate for stressors along each of these three axes. For example, M. spicatum exhibits no significant photosynthetic advantage for low-light conditions, but instead tolerates low-light conditions through the ability to elongate and form a canopy at the water surface where light is sufficient for growth [101,102]. For P. crispus, an inverted reproductive cycle and the ability to grow under low temperatures enable it to persist in deeper and more eutrophic waters at times of year when algal biomass is lower and light availability is sufficient for photosynthesis [51,52,72,103]. In short, P. crispus can capitalize on early phenology to grow in locations that have insufficient light later in the growing season. In ecosystems with strong seasonality, as in northern temperate lakes, the influence of temperature on community structure can occur through differences in occupancy of different phenological niches [41,42]. The role of phenology in invasions remains an area of active research, especially as it relates to competition avoidance in early stages of invasion [104].
Even with this dataset, which is unusually large in comparison to other macrophyte studies, there are limitations to the retrospective niche modeling approach using observational data. For example, data demands of the PDF niche model used here meant that only 34 of 163 recorded species had sufficient data to develop three-dimensional niche models. Had we considered a fourth niche dimension (e.g., soluble nutrient concentrations), modeling would have required 10,000 occurrences per species and limited our analysis to just 12 species. In addition, the scale of the dataset meant that measurements of some environmental data were not available. For example, water chemistry data have been collected in a subset of the lakes analyzed, but the need for complete data records for our analysis meant that these data were not useable for our study. We believe that our use of weather station data to calculate annual GDD is reasonable, but lake-level water temperatures would be more accurate-these data, however, are not available. As a consequence, the models we have derived may be shifted with respect to the actual thermal regime of a lake at any given time (i.e., average air temperature will generally increase faster than lagging water temperatures throughout the year), but because of the comparative nature of this study, all modeled niches are subjected to the same direction and magnitude of this potential bias. It is important to note that we have not captured every factor that influences macrophyte niches in this study. For some cases, such as free-floating plants (which have no physiological reason to be constrained by water depth, yet show a strong depth pattern in their niche occupancy; Figures S1 and S3), observed patterns likely reflect the effects of an unobserved covarying factor, such as wave-induced surface movement. In addition, nutrient loading has been associated with dominance of invasive species [8]. Light is often, but not always driven by nutrient concentrations in Minnesota lakes [105], thus future work should explicitly explore the mechanism by which nutrient enrichment promotes invader dominance in these systems. We acknowledge these limitations and suggest that each be addressed or expanded upon in the future.
It is rarely feasible to evaluate the full niche of an invasive species in its invaded range, where it is unlikely to have yet reached the full extent of its potential spread. Moreover, it is difficult to infer whether co-occurrences with native species translate to antagonistic interactions [106]. However, we can make preliminary inferences regarding interactions between native and invasive species when there is suitable information to characterize the niche of an invading species, as enabled by the large state-wide dataset we used here. By evaluating multiple niches across many species, we were able to draw inferences about the extent to which species may be competing or avoiding competition. Experiments to directly quantify the strength of competitive interactions between these invaders and resident plant communities should subsequently test our finding that P. crispus exhibits less competitive interaction with native species than M. spicatum.