3.1. Species Richness
A total of 291 plant species were identified within the 54 sites surveyed. For the dry basins, plant species with a frequency of occurrence of more than 70% were Festuca rubra, Lathyrus pratensis, Vicia cracca, Plantago major, Taraxacum officinale, Trifolium pratense and Trifolium repens,which are mostlygrass and melliferous species (Supplementary Information S3). For wet basins, the wetland species Carex scoparia, Galium palustre, Juncus vaseyi, Scirpus atrovirens, Typha sp., Juncus effusus, Lythrum salicaria and Equisetum arvensis were the most frequent (more than 70%). For channelized wet basins, plant species with a frequency of occurrence of more than 70% were Solidago canadensis, Vicia cracca, Juncus vaseyi and Symphyotrichum lanceolatum, as well as Equisetum arvensis, Lythrum salicaria and Typha sp., that were found in all these sites. Few species were frequent in ponds, with only Lythrum salicaria and Typha sp. showing a frequency of occurrence of more than 70% of the sites.
We identified 84 families within the subset of 28 sites selected for the macroinvertebrate survey. For dry basins, a total of 20 families were found at a frequency of occurrence of more than 70%, among them Apidae as well as 14 families found on all sites, such as Brachycera, Hymenoptera and Nematocera (see Supplementary Information S4). For wet basins, a total of 11 families were captured on more than 70% of the sites, including Aphididae, Brachycera, Cicadellidae and Heteroptera, which were found on all sites. For channelized wet basins, 25 families were sampled at a frequency of occurrence of more than 70%, including 12 families found on all sites, such as Brachycera, Hymenoptera, Nematocera and Lepidoptera. Ponds had only seven families frequently encountered on more than 70% of the sites. Only Aphididae, Brachycera and Nematocera were captured at all ponds.
Differences in specific richness between infrastructure types were assessed using species richness accumulation curves with a significant increase in richness as the number of quadrats sampled increases (Fig.2). Total plant richness was the lowest for dry basins (Fig.2a). Wet basins had greater diversity than ponds, with channelized wet basins showing intermediate values (Fig 2a). The higher specific richness in wet basins was associated with a greater FACW species richness (Fig.2c). When the banks were excluded, ponds had significantly fewer species than all other types of basins (Fig.2b), with the lowest FACW species richness as well (Fig.2d). Wet basins and channelized wet basins were similar in terms of plant richness whether the banks of the sites were included or not (Fig.2a, 2b). The same analyses of OBL species richness showed no significant differences between wet basins, channelized wet basins and ponds, whether including or excluding the banks (results not shown). A comparison between ponds, wet basins and channelized wet basins focusing on OBL and FACW species proportion of the total pool of sampled species, excluding the banks of the infrastructure, showed a greater proportion of FACW species in wet basins than in ponds (Fig.2g) and inversely, a greater proportion of OBL species in ponds than in wet basins (Fig.2h), with channelized wet basins showing intermediate proportion values for both. Overall, proportion of OBL species decreased with the number of quadrats for the three infrastructure types.
Macroinvertebrate richness showed no difference among the four types of stormwater infrastructures (Fig.2e). However, when the banks were excluded, dry basins showed greater family richness than ponds (Fig.2f), indicating that the banks attracted more macroinvertebrate families compared to the bottom of the ponds.
3.2. Beta Diversity and Community Composition
Beta diversity was used to assess the variability in community composition among the studied sites. Dry basins had significantly lower beta diversity than the other types of infrastructure when data from the banks were included (Table 1; Fig.3a). When the banks were excluded, channelized wet basins had the greatest beta diversity, followed by wet basins, whereas dry basins and ponds had the lowest (Table 1; Fig.3b). Ponds had a high variability in beta diversity. Indeed, four ponds were markedly different in community composition from the others and caused greater variability, even though average beta diversity was low. Plant communities differed between types, except for wet basins and channelized wet basins, which had a similar composition (Table 1; Fig.3a, 3b; centroid position) whether the data from the banks were considered or not.
Beta diversity analyses for macroinvertebrates yielded patterns like those observed for the vegetation. Dry basins had significantly lower beta diversity than the other infrastructure types when data from the banks were included (Table 1; Fig.3c). When excluding the banks, channelized wet basins had the highest beta diversity, followed by wet basins and dry basins, while ponds showed an intermediate beta diversity between that of dry and wet basins (Table 1; Fig.3d). Dry basins, channelized wet basins and ponds each hosted different communities, whether considering banks data or not (Table 1; Fig.3c, 3d, centroid position). Nevertheless, wet basins showed comparable communities with both ponds and channelized wet basins. The relationship between plant and macroinvertebrate communities was supported by a significant Mantel correlation between the two taxa groups (observation: 0.58, p-value < 0.001).
The differences in vegetation community composition between infrastructure types are shown in a PCA (Fig.4a). The first two axes explain 29.4 % of the variation. The results revealed the importance of grasses as well as melliferous plants such as Trifolium repens, Prunella vulgaris and Taraxacum officinale in dry basins. OBL species characterizing ponds included Typha sp., Schoenoplectus tabernaemontani and the invasive FACW species Phragmites australis. Typha sp. and Phragmites australis tended to cover higher portions of the communities in channelized wet basins and ponds compared to the other types of infrastructures (Fig.5). Frequency of occurrence of Typha sp. was 0%, 70%, 100%, 100% for dry basins, wet basins, channelized wet basins and ponds, respectively, and for Phragmites australis, 0%, 20%, 53%, 50%. Similar communities in wet basins and channelized wet basins were represented by FACW species such as Juncus vaseyi, Scirpus rubrocinctus, S. atrovirens, Juncus effusus, Phalaris arundinacea, Carex scoparia, Symphyotrichum lanceolatum and Galium palustris. Generalist species such as Equisetum arvense, Trifolium hybridum and Lathyrus pratensis were also commonly present in these types of infrastructures (Fig.4a).
Macroinvertebrate community composition was also represented using a PCA, with differences in families plotted according to the type of infrastructure (Fig.4b). The first two axes explain 40.5 % of the variation. Hymenoptera and Apidae tended to be associated with dry basins, whereas Nematocera, Brachycera and Coenagriidae, which represent mosquitoes, flies and dragonflies respectively, were more commonly found in ponds. Heteroptera, Cicadellidae, Membracidae and Miridae were well represented in wet basins and channelized wet basins. Pollinators were supplemented by higher abundances of Apidae and Hymenoptera (including several species of pollinators) in dry basins compared to ponds (Fig.5). However, the high variability in abundance among dry basins must be considered. Frequency of occurrence of Apidae and Hymenoptera was 71%, 28%, 28%, 14% and 100%, 86%, 100%, 86% for dry basins, wet basins, channelized wet basins and ponds, respectively.
3.3. Landscape Variable Analyses
Out of the 18 variables tested (Supplementary Information S2), only the presence of natural areas within a 2000 m radius (Nat2000) was found to be associated with plant composition in stormwater infrastructures and accounted for 15.2% of the variation (Fig.6; first axis). Dry basins were excluded from this analysis because of the particularly homogenous and poorly diversified plant composition in these highly disturbed sites (frequent maintenance). A higher proportion of natural areas favoured FACW species such as Scirpus rubrotinctus, S. atrocinctus, Carex echinacea, C. garberi, C. crinita, whereas generalist species such as Lathyrus pratensis, Phleum pratense, Agrostis stolonifera, Vicia sepium, Solidago canadensis, and Lotus corniculatus, were noticeably present with less natural areas around. The second axis seemed to reflect the water availability gradient (8.7% of the variability explained) with the OBL species Schoenoplectus tabernaemontani and Typha sp. associated with ponds and a higher number of FACW species such as Carex scoparia, C. vulpinoidea and Phalaris arundinacea occurring in wet basins, whether in the presence of a channel or not. The same analysis was conducted including only wet basins and channelized wet basins sites, and Nat2000 was again the only significant variable (15.5% of variation explained, plot not shown). For ordination analyses on ponds alone, no variable was found to be significantly associated with plant composition. The same analyses were conducted on macroinvertebrate data and showed a significant influence of natural areas as well as vegetation cover mapped from NDVI index on a finer scale within a 500 m radius. The two variables explained 41% of the variation (plot not shown). Neither local variables such as site surface area or distance to main roads, nor the distance between one site and other sites were significantly associated with plant and macroinvertebrate community structures.