Connectivity and seasonality cause rapid taxonomic and functional trait succession within an invertebrate community after stream restoration

PLOS ONE, May 2018

General colonization concepts consent that a slow process of microhabitat formation and subsequent niche realization occurs during early stages after new habitat is released. Subsequently, only few species are able to colonize new habitat in the early onset of succession, while species richness increases steadily over time. Although most colonization studies have been performed in terrestrial ecosystems, running water ecosystems are equally or even more prone to colonization after disturbance due to their dynamic nature. We question how invertebrate succession patterns reconcile with general colonization concepts. With this study we provide insight into the colonization process in newly created lowland stream trajectories and answer how within-stream bio- and functional diversity develops over time. Our results show a rapid influx of species, with a wide range of functional traits, during the first season after water flow commenced. During more than two years of regular monitoring, immigration rates were highest in autumn, marking the effects of seasonality on invertebrate dispersal. Biodiversity increased while abundance peaks of species alternated between seasons. Moreover, also days since start of the experiment explains a considerable part of the variability for taxa as well as traits. However, the relative trait composition remained similar throughout the entire monitoring period and only few specific traits had significantly higher proportions during specific seasons. This indicates that first phase colonization in freshwater streams can be a very rapid process that results in a high biodiversity and a large variety of species functional characteristics from the early onset of succession, contradicting general terrestrial colonization theory.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0197182&type=printable

Connectivity and seasonality cause rapid taxonomic and functional trait succession within an invertebrate community after stream restoration

May Connectivity and seasonality cause rapid taxonomic and functional trait succession within an invertebrate community after stream restoration Judith J. Westveer 0 1 Harm G. van der Geest 0 1 E. Emiel van Loon 0 1 Piet F. M. Verdonschot 0 1 0 Institute for Biodiversity and Ecosystem Dynamics, University of Amsterdam , Amsterdam , The Netherlands , 2 Wageningen Environmental Research , Wageningen UR, Wageningen , The Netherlands 1 Editor: Judi Hewitt, University of Waikato , NEW ZEALAND General colonization concepts consent that a slow process of microhabitat formation and subsequent niche realization occurs during early stages after new habitat is released. Subsequently, only few species are able to colonize new habitat in the early onset of succession, while species richness increases steadily over time. Although most colonization studies have been performed in terrestrial ecosystems, running water ecosystems are equally or even more prone to colonization after disturbance due to their dynamic nature. We question how invertebrate succession patterns reconcile with general colonization concepts. With this study we provide insight into the colonization process in newly created lowland stream trajectories and answer how within-stream bio- and functional diversity develops over time. Our results show a rapid influx of species, with a wide range of functional traits, during the first season after water flow commenced. During more than two years of regular monitoring, immigration rates were highest in autumn, marking the effects of seasonality on invertebrate dispersal. Biodiversity increased while abundance peaks of species alternated between seasons. Moreover, also days since start of the experiment explains a considerable part of the variability for taxa as well as traits. However, the relative trait composition remained similar throughout the entire monitoring period and only few specific traits had significantly higher proportions during specific seasons. This indicates that first phase colonization in freshwater streams can be a very rapid process that results in a high biodiversity and a large variety of species functional characteristics from the early onset of succession, contradicting general terrestrial colonization theory. - Funding: This study was funded by the foundation for applied water research (http://stowa.nl/ foundation_for_applied_water_research_stowa/). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Introduction Colonization is a key concept in community ecology and its study has revealed many mechanisms by which ecological communities commence and develop during ongoing succession [1±8]. Most theories describe either facilitating, tolerating or inhibiting interactions between species during succession. The main resemblance in each theory, however, is that the early successional stage is characterized by a slow sequence of microhabitat formation and houses a species poor community. Successive to the colonization of vegetation is the appearance of invertebrates, as plants provide nutrients and refuge [ 9, 10 ]. While substantial redundancy exists among theories on succession, there are still many knowledge gaps. Especially when it comes to identifying key parameters that determine the recolonization of fauna [11]. Remarkably, most colonization studies have been performed in terrestrial ecosystems, even though disturbance and subsequent succession is common in lotic ecosystems due to their dynamic nature [ 12, 13 ]. Aquatic invertebrates have traits to quickly disperse to recently released habitat [14±18], but previous studies on stream colonization have shown that succession does not merely depend on species characteristics or time itself [19]. Invertebrate succession is rather shaped by a combination of 1) the distance between the regional species pool and new habitat [20±22], 2) the presence of a suitable habitat for settling and colonization [ 23, 24 ], 3) dispersal capacity and life history traits of present species [17, 25±28], and 4) the timing of dispersal [ 29, 30 ]. Unconnected terrestrial habitat patches rely on wind-mediated dispersal, along with animal and sometimes even human vectors [31±34]. However, freshwater streams are often connected throughout a catchment and water flow is regarded as the number one facilitator of dispersal in streams [ 27 ]. As general colonization concepts seem to match terrestrial colonization more than flowing water colonization, we aim to unravel invertebrate successional patterns in newly created trajectories in temperate lowland streams. In order to do this, we monitored three restored stream trajectories that have been connected to an existing stream channel. As previous studies have shown, the recolonizing community is shaped by the distance and composition of the surrounding species pool [ 21, 22 ] and habitat formation is affected by the inflow of allochtonous material providing a fast formation of refugia [ 35, 36 ]. This leaves us to address 1) how does taxonomic and functional trait diversity develop over time in a new stream trajectory, 2) which traits regarding dispersal capacity and life history are important during early succession, and 3) how dispersal success is affected by timing. Our expectation is a rapid increase in taxonomic as well as functional diversity right after water flow commences. New species, with different suits of functional traits, will arrive throughout the first years while the habitat is developing and more niches become available. Metacommunity studies often neglect temporal patterns (such as ecological colonization and succession) in the interpretation of community composition, since it requires long-term and extensive monitoring. While the outcome of many metacommunity studies have led to an increased understanding of how spatial dynamics and local interactions structure communities, more empirical evidence is needed in order to fully comprehend the composition patterns in newly created habitats [37]. Therefore, we focus on temporal shifts in community composition in a restored stream. The observed temporal shift cannot be viewed independently of larger scale spatial context, but the outcome will support a mechanistic explanation for the speed at which stream communities develop after disturbance and lead to a clearer understanding of early succession patterns in lotic systems. Material and methods Study area The area studied is located in the Leuvenumse stream (52Ê18'55.17ºN; 5Ê42'33.63ºE), in the province of Gelderland, the Netherlands. This is a slow-flowing, meandering lowland stream 2 / 17 (flow velocities ranged from 8 cm/s in winter up to 46 cm/s during spates in summer), fed by both precipitation and ground water. On average, the studied stream is about 4 meters wide. The deepest part of the channel (thalweg) varies between 15±53 cm, depending on season and location. Parts of the catchment are used by agriculture. The stream itself alternates through open pasture and deciduous and coniferous woodland. The streambed substrate consists of sand and clay with patches of gravel and coarse particulate organic matter. Throughout the monitoring period, the water temperature shifted seasonally from a minimum of 3.9 ÊC in winter up to 17.6 ÊC in summer. Besides water temperature, dissolved oxygen, pH and conductivity were measured at each sampling moment (Table 1). In the autumn of 2014, three former stream trajectories were re-connected to the existing stream by blocking parts of the main channel with sand (Fig 1). This caused all stream water to flow through the re-connected trajectories that had been dry forest floor without any pre-restoration communities up until then. All data regarding invertebrate community and morphological, physical and chemical characteristics derive from these three trajectories. Permission to enter the national park and sample biota in the re-connected trajectories was provided by local authority NGO Natuurmonumenten with permit number 035-2014-430. Sample collection and processing Within the first 6 months after water flow commenced (SeptemberÐMarch 2014), the invertebrate community was sampled every two weeks in all three re-connected stream trajectories. After 6 months, samples were taken once per season as we assumed that community saturation slowed down the colonization rate. To sample the entire invertebrate community, three multiplate samples and two handnet samples were taken at each trajectory and each sampling moment. Hester-Dendy multiplate samplers (total sampling area is 0.10 m2 per sampler) were used to exclude habitat variation during colonization and placed on the stream bottom of the three new stream trajectories as soon as water started flowing (t = 0). Net samples (mesh size 250 μm, width of net is 30 cm) consisted of a multi-habitat sample composed of 1 m of mineral and 1 m of organic habitat, and were included to cover all habitat variation present at each sampling moment. Samples were transported to the laboratory, where they were rinsed and the invertebrates picked out and stored in ethanol (70%) within 12 hours after being collected from the stream. In the laboratory, each individual was identified up to species-level if possible (Crustacea, Gastropoda, Hirudinae, Insecta: Trichoptera, Ephemeroptera, Plecoptera, Coleoptera, Simuliidae). Juveniles and some invertebrate-groups were identified up to genus-level (Bivalvia, Insecta: Chironomidae). Oligochaetes and Arachnida were excluded from the analyses due to limited identification at species-level. Care was taken to ensure that taxonomic resolution was sufficient according to Haase et al. [ 38 ]. Identification was performed with the use of a dissecting microscope (120x magnification) and a light microscope (300x magnification). 3 / 17 Fig 1. Simplified map of the Leuvenumse stream catchment with the three reconnected former stream trajectories (indicated with red circles). The arrows show the direction of the water flow. The inset shows the most upstream location; the existing channel was blocked with sand (red dots), after which the water took its natural course through the former streambed (green dots). 4 / 17 Functional trait data Data on functional traits was provided by freshwaterecology.info, a database that contains taxonomic and ecological information on freshwater macroinvertebrates [ 39 ], combined with trait-database EKOO [ 40 ] and Tachet's et al. database [ 25 ]. All trait data is number coded using a 10-point system or was transformed to fit a 10-point system. We used the following six trait categories for our analyses, based on the level of importance for colonization and presence of existing knowledge for our community: feeding type, locomotion type, dispersal mode, rarity-score, number of reproductive life cycles per year and r/K-strategy. Selected functional trait data was available for respectively 94%, 92%, 68%, 65%, 41% and 23% of the 96 taxa that were found in the new trajectories during the sampling period. Statistical analyses The multi-habitat and multiplate samples were pooled since the combination of taxa found by both methods represented the local community at each sampling moment. It appeared that habitat variation did not affect the data composition. No aquatic plants were present in any of the sampling sites, and streambed substrate coverage remained largely similar over time (data not included). The abiotic parameters (Table 1) did not shift over time, apart from some specific and well-known relationships between season and parameter (temperature, water depth) and were therefore not included in the analysis as an explanatory factor. Samples of each of the three trajectories were pooled due to high community similarities (Jaccard similarity coefficient = 0.80). This enhanced the robustness of the analyses, and averaged out any local and small-scale effects. The development of diversity and abundance over time is described with a species accumulation curve, while explicitly keeping track of new colonists. Functional diversity is investigated with the same linear model (involving season and time as predictors) on a range of trait characteristics. The development of diversity at the species level is explored by clustering omnipresent species abundance patterns over time, with a non-metric multidimensional scaling (NMDS) ordination plot [ 41 ] and a distance based redundancy analysis [ 42 ]. To unravel community shifts in functional trait composition over time, community weighted trait means included taxa abundance per sampling moment. The explained variance in the NMDS by seasonality as well as days since the start of the experiment were evaluated. Next, the partial contribution by season and days since start were evaluated by a distance based RDA and tested by a permutationbased ANOVA at a 0.05 significance level. To test if specific functional traits were of increasing or decreasing importance during succession, linear regressions of functional traits over time were calculated. Only taxa presence/ absence data was considered in this case, multiplied with specific trait values and divided by the total number of species present per sampling moment. All analyses were conducted in R, using functions from the package vegan for NMDS and distance-based redundancy analysis. Results Species composition In total, 96 different taxa were found in the new trajectories during the entire monitoring period of 735 days. Species diversity ranged between 23 and 45 taxa per sampling moment (Fig 2). 5 / 17 Fig 2. Number of taxa over time, measured 18 times from the moment water started flowing through the reconnected trajectories of the Leuvenumse stream. All three trajectories were pooled due to community similarities (Jaccard Jij = 0.8). Number of taxa is plotted per sampling moment (absolute), # new taxa arriving (immigration) and # of total taxa (cumulative) over time (days since start). The colonization process started with a rapid increase in taxonomic diversity, with 23 taxa being present within the first 14 days after water flow commenced (t = 1). After 60 days (t = 4), already more than 50% of all taxa collected was present. The highest numbers of newly identified immigrating taxa is at the start in autumn (t = 1, 2 and 3) and subsequent autumn sampling moment (t = 14), exactly one year later. Towards the end of the experiment, a decreasing amount of new colonizers was observed. A linear model was used to predict differences in immigration rate based on season and time since start. A significant regression was found (F4,13 = 5.12, p = 0.01, R2 = 0.61), with more taxa arriving in autumn than winter. Insects were most abundant throughout the study, ranging from 141 to 949 individuals per sample moment (average 640 ± sd 211, Fig 3). Bivalves showed abundance peaks with high numbers of individuals at the start of the colonisation (max. 642 individuals), but low abundances towards later sampling moments (average 105 ± sd 158 individual bivalves per moment). Crustaceans steadily became more abundant over time, ranging from 17 at t = 1 to 407 individuals at t = 13 (average of 44 ± sd 98 individual crustaceans per moment). The 6 / 17 Fig 3. Changes in abundance (log (# individuals)) per taxonomical macroinvertebrate class over time (# days since start). Insects form the largest class of colonist from the early onset of succession, maintaining equal levels of large abundance over time. depression in bivalve and crustacean abundance at t = 6 is probably due to the occurrence of a heavy spate at the moment of monitoring. The four most abundant orders of insects showed annual patterns with a large variation in abundance over time (Fig 4). Some taxa showed abundance peaks midwinter (Panel B. Plecoptera; N. cinerea), while other taxa were most abundant in summer (Panel C. Trichoptera; H. radiatus) or spring (Panel D. Coleoptera; O. villosus). Peaks in larval abundance are often found right after adults have oviposited their eggs. Univoltinistic patterns are found for most taxa, while some (e.g. Micropsectra sp.) show semivoltine or flexivoltine life cycles with more than one peak in abundance per year. Annual patterns were similar in both years. Functional trait composition Of the 29 traits considered in this study, divided into six different trait categories, 25 traits were present from the first sample moment onward. The overall majority of these traits (20 out of 29 traits) remained present without significant decreases and increases in relative occurrence over time (Fig 5). We do see significant increases in number of shredders, active filter 7 / 17 Fig 4. Changes in abundance over time (days since start) of four omnipresent taxa of insects, divided by order. A: Diptera, B: Ephemeroptera/ Plecoptera, C: Trichoptera, D: Coleoptera. feeders and burrowers (Table 2) over time. Significant increases were also found for taxa with univoltine and semivoltine reproductive cycles over time. Semi-rare taxa increased over time while rare and very rare taxa did not show significant changes. Passive filter feeders and aerial passive dispersers decreased over time. Predator taxa were most abundant (average of 34% of the total community) throughout the whole study period compared to gatherers, grazers, miners, filter feeders and shredders. Active walking taxa were most abundant (48%) throughout the whole period compared to active swimming, drifting, sessile and burrowing taxa. Taxa with a r-strategy were most abundant (87%) compared to K-strategists. Univoltinistic taxa were most abundant (62%) compared to bivoltine, trivoltine, semivoltine and multivoltine taxa. Taxa with different rarity-scores were all present in equal ratios (18%-22%), apart from the very rare taxa, which are not abundant (6%) at all times (Fig 5). Equal ratios were also present for aquatic passive and aerial active taxa, which both account for 32% of the total community. A NMDS of the taxa (Fig 6A) and trait composition (Fig 6B, Table 3) show that seasonality as well as days since start of experiment are important factors shaping the community composition. For the taxa-ordination, the stress value at the optimal solution is 0.159. The R-squared 8 / 17 Fig 5. Changes in functional trait composition (relative trait presence) over time of six different trait categories. Linear regressions of functional traits over time were calculated to test if specific functional traits were of increasing or decreasing importance during succession. Only taxa presence/absence data was considered in this case, multiplied with specific trait values and divided by the total number of species present per moment. Regression lines and coefficient are only shown for significant regressions. A: Feeding group, significant regressions: shredder, active filter feeder, passive filter feeder, B: Locomotion type, sig. regr.: burrowing, C: Dispersal mode, sig. regr.: aerial passive, aerial active, D: Rarity, sig. regr: semi-rare, E: Reproductive cycle, sig. regr: semivoltine, univoltine, F: r-K selection strategy, no sig. regr. 9 / 17 Trait class Shredder active filter feeder passive filter feeder Burrowing semi-rare aerial passive aerial active Semivoltine Univoltine for the non-metric fit between observed dissimilarity and ordination distance is 0.975. Season explains the variation of the taxa in the reduced space with an r2 of 0.598, while days since start explains it with an r2 of 0.710. When analysing the contribution by season and days since start in a distance-based RDA, both factors also turn out to be significant (permutation-based pvalue of 0.001 for season and 0.002 for days), with r2 values of 0.16, 0.40 and 0.51 for model including respectively days since start, season and days since start+season. For the trait-ordination, the stress value at the optimal solution is 0.11. The R-squared for the non-metric fit between observed dissimilarity and ordination distance is 0.987. Season explains the variation of the traits in the reduced space with an r2 of 0.418 and days since start with an r2 of 0.574. A distance based-RDA on the traits shows that both factors are significant (permutation-based p-value of 0.002 for season and 0.004 for days), with r2 values of 0.27, 0.41 and 0.58 for model including respectively days since start, season and days since start+season. So, the differences between the observation times in terms of taxa as well as trait diversity can be represented well in a reduced (two-dimensional) space both season and days since start of the experiment explain a considerable part of the variability for taxa as well as traits. Discussion This study gives insight into the shifts in community composition on a taxonomic as well as a functional trait-level during the early stages of colonization of a temperate lowland stream. Our results show that new trajectories of temperate streams are characterized by a rapid colonization by multiple species with different traits in large abundances. This observational study was performed in a highly connective system with constant supply of allochtonous organic material. These parameters could be the agent of success that caused a high rate of succession. Circumstances that are thought to increase succession rate are small distances to the regional source pool [ 22, 34 ], species with active dispersal mechanisms [ 43, 44 ], specific seasonal effects [ 30 ] and suitable new habitat. We will discuss how these four elements have played a role in our study and how they affect the process of succession in temperate lowland streams. Distance and dispersal mechanisms Parkyn and Smith [ 45 ] consider dispersal constraints low if a new habitat has an upstream supply of invertebrates and additional aerial recolonists are no further than 2 km away. Under such conditions the new habitat could resemble a reference community within 10±50 years after disturbance. Sundermann et al. [ 22 ] and Winking et al. [ 17 ] consider dispersal constraints high when there is no direct upstream community and recolonists are more than 5 km 10 / 17 Fig 6. Non-metric multidimensional scaling plot (NMDS; a rank-based approach to ordination, representing the pairwise dissimilarity between all cases). For both the species (A) and the trait-based (B) analysis, the data is fitted to two axes, using the Bray-Curtis distance metric. 11 / 17 away. Those communities are very unlikely to resemble a reference community ever again. Recent studies have shown that succession can be a slow process if the connectivity between streams is poor [ 45, 46, 47 ]. In our case, the re-connected trajectories were directly, and imperatively, connected to the upstream and downstream existing stream channel and therefore dispersal constraints were considered low. Furthermore, hydrological connectivity and riparian vegetation was present at all times, which enabled rapid colonization of aerial and aquatic dispersers from the start. Even though there is still a debate on which pathway makes up the majority of dispersal, it has been found that short-distance dispersal mainly occurs due to aquatic drift [ 22, 45, 48 ] while long-distance dispersal by adult insect flight is often the primary mechanism for recolonization of restored trajectories in separate catchments [ 49, 50 ]. However, many species are incapable of actively dispersing themselves and rely on animal vectors, wind, or water flow to provide passive transport between sites [19]. In terrestrial ecosystems, wind-mediated dispersal is the prime mechanism that enables early stage succession [ 51 ]. In running waters, flow facilitates colonization from the first moment on [ 52, 53 ] and connectivity to source populations will further speed up arrival. All four dispersal mechanisms (aquatic/aerial active/passive) were present among taxa from the onset of succession. Aerial passive dispersers decreased over time, all other dispersal mechanisms remained present at equal ratios over time. This finding contradicts recent studies where aerial active generalists were found to colonize most rapidly and weakly dispersing generalists immigrated much later [ 17, 28 ]. This contradiction and our observation leads us to hypothesize that categorizing species by dispersal mechanisms alone does not give sufficient information on actual dispersal rate. A clearer view of how the community developed during the process of succession is obtained by looking at the functional trait composition. We found that while some traits were dominant over others, overall trait composition did not change over time. For terrestrial systems, new habitat is colonized by species that have rapid life cycles and are r-strategists. In later stages, species with competitive traits, such as predatory behavior, are able to settle in the new habitat. Analogous, Gore [ 20 ] showed that collector-gatherers and collector-filterers were the initial colonizers of stream habitat instead of predators. Surprisingly, the traits found to be dominant in this study contradict the colonizer-competitor trade-off concept [ 7, 54, 55 ] and the findings of Gore [ 20 ]. From the onset of succession, predatory species were present in the new habitat. The diversity of predators was even higher than the diversity of filter feeders. In 12 / 17 accordance with other studies [ 20, 28, 48, 56 ] our study shows that the trait composition of the newly arriving species community probably depends on the regional species pool in the proximity, rather than develop according to previously theorized sequence of functional traits. Succession and seasonal effects Our results indicate that there is a temporal trend for the duration of the project (viz the results by distance-based RDA), which might be successional. On top of that there are clear seasonal effects (also shown in distance-based RDA). The results show that most species either immigrate in the first season or exactly one year later in the same season. Both immigration events took place in autumn. This indicates the importance of seasonality during the colonization process in temperate running waters. Seasonality can effect dispersal processes directly, in terms of optimal hydrological connectivity [ 57 ] and species availability due to life-cycle timing [58±60], and indirectly, in terms of increases and decreases of specific environmental parameters that trigger organisms to disperse [ 30 ]. Peak flows after heavy rainfall are common in Dutch lowland streams during the end of summer and autumn. The sudden increase in current velocity can cause species to become dislodged from the substrate and enter a state of drift [ 61 ]. Additionally, an increase in the volume of water due to melting or rainfall can expand potential habitat surface and increase hydrological connectivity. This improves instream dispersal and colonization patterns [ 62 ]. In the late summer and early autumn temperate lowland streams are populated by many juvenile and first larval stages specimens of many species [63± 65]. These young animals are often susceptible to drift, either passively or actively as dispersal trait, and prone to transportation downstream where they can settle. Throughout the year, also some new species immigrate into the denuded habitat. This can be explained by both ongoing habitat development and the ability of species to inhabit these new habitats due to their functional traits [66, 67]. Interestingly, species from different insect orders show annual peaks in abundance, but alternate in the timing of these peaks. By partitioning resources and interspecific competition, more species are able to coexist while total abundance levels are stable and succession continues [ 12 ]. The result that species and trait composition is correlated to seasonality indicates that seasonal life history characteristics provide the opportunity to colonize, survive and reproduce and not just time since existence or dispersal capacity alone. This finding is supported by previous studies on stream restoration effects [ 56 ], and can be explained by the different life history strategies needed to exploit (a)biotic resources [ 67, 68 ]. Surprisingly, little attention has been given to seasonal patterns in macroinvertebrate colonization in temperate streams. In our study, the combination of strong hydrological connectivity, high availability of juveniles, and availability of habitats and resources of re-connected trajectories set the requisites for enhanced autumn colonization. Habitat suitability Once the macroinvertebrates arrive at the newly restored habitat, they have to be able to settle and reproduce before recruitment can be considered successful. However, the presence of favorable in-stream and riparian habitat alone is not enough for community recovery [69]. Distances between leaf patches and sizes of these organic refuges have a profound effect on species survival [70, 71]. Species might not be able to recolonize the stream before multiple new microhabitats have established in the stream [ 35, 36 ]. The rapid increase in biodiversity and functional diversity right after water flow commenced indicates that microhabitat formation was not a limiting factor for the first species to arrive. High hydrological connectivity between the upstream source pool of organic matter and the reconnected trajectory together 13 / 17 with the input of allochtonous material from the autumn-shed leaves along this temperate lowland stream could have provided a fast formation of suited habitats. To conclude, this study shows that temperate lowland stream colonization is marked by a rapid increase in species richness and abundance right after water flow commences, when hydrological connectivity is not a constraint. The community is heterogeneous in terms of functional trait diversity from the early onset, with no clear sequence in colonizer to competitor trait characteristics. The immigration rate of new species is affected by seasonality. To strengthen the interpretation of these observations, data from nearby established streams help to standardize the community composition by season and determine at what time point the assembled communities of new streams become statistically indistinguishable from established streams. This could unravel the proportion of both spatial and temporal effects on successional patterns. Nonetheless, the observed patterns show that colonization and subsequent succession rate can be high in restored stream trajectories, yet effects of the regional species pool might be a limiting factor in community restoration and worth to explore in future studies. Acknowledgments This research was funded by the Dutch foundation for applied water research (STOWA). We would like to acknowledge the NGO Natuurmonumenten and Waterboard Vallei and Veluwe for field support, Dorine Dekkers for her help with the identification of aquatic macroinvertebrates. Furthermore, we would like to thank Paula dos Reis Oliveira, Stella Blok, Mariska Beekman, Remy Bos, Luiza Dolcetti Domingos and Daniela Cueva Arias for their assistance in the field. All data is stored at Zenodo data repository and openly accessible via DOI code: 10.5281/ zenodo.1149125. Author Contributions Conceptualization: Judith J. Westveer, Harm G. van der Geest, Piet F. M. Verdonschot. Data curation: Judith J. Westveer, E. Emiel van Loon. Formal analysis: Judith J. Westveer, E. Emiel van Loon. Supervision: Harm G. van der Geest, Piet F. M. Verdonschot. Writing ± original draft: Judith J. Westveer. Writing ± review & editing: Harm G. van der Geest, E. Emiel van Loon, Piet F. M. Verdonschot. 14 / 17 15 / 17 37. 16 / 17 Higler LWG. Verspreidingsatlas Nederlandse Kokerjuffers (Trichoptera). EIS Nederland press; 2008. 1. Clements FE. Plant Succession : An Analysis of the Development of VegetationÐCarnegie Institution of Washington; 1916 . 2. Egler FE . Vegetation science concepts I. Initial floristic composition, a factor in old-field vegetation development . Plant Ecology . 1954 ; 4 : 412 ± 417 3. MacArthur RH , Wilson EO . The theory of island biogeography . Princeton University Press; 1967 . 4. Odum EP . Strategy of ecosystem development . Science . 1969 ; 164 : 262 ± 270 . PMID: 5776636 5. Connell J , Slatyer R . Mechanisms of Succession in Natural Communities and Their Role in Community Stability and Organization . Am Nat . 1977 ; 111 ( 982 ): 1119 ± 44 . 6. Tilman D. The resource-ratio hypothesis of plant succession . Am. Naturalist . 1985 ; 125 : 827 ± 852 7. Huston MA , Smith T. Plant succession: life history and competition . Am. Naturalist . 1987 ; 130 : 168 ± 198 . 8. Cadotte MW . Competition-Colonization Trade-offs and Disturbance Effects at Multiple Scales . Ecology . 2007 ; 88 : 823 ± 829 . PMID: 17536699 9. Southwood TRE , Brown VK , Reader PM . The relationships of plant and insect diversities in succession . Biol J Linn Soc . 1979 ; 12 ( 4 ): 327 ± 48 . 10. Schowalter TD . Insect ecology, an ecosystem approach . 4th ed. Academic Press; 2016 11. Pulsford SA , Lindenmayer DB , Driscoll DA . A succession of theories: purging redundancy from disturbance theory . Biol Rev . 2016 ; 91 : 148 ± 67 . https://doi.org/10.1111/brv.12163 PMID: 25428521 12. Townsend CR . The patch dynamics concept of stream community ecology . J North Am Benthol Soc . 1989 ; 8 ( 1 ): 36 ± 50 . 13. Lake PS . Disturbance, patchiness, and diversity in streams . Am Benthol Soc . 2000 ; 19 ( 4 ): 573 ± 92 . 14. Waters TF . Recolonization of denuded stream bottom areas by drift . Trans. Am. Fish. Soc . 1964 ; 93 : 311 ± 15 . 15. Gore JA , Milner AM . Island biogeographical theory: can it be used to predict lotic recovery rates? Environ . Management. 1990 ; 14 : 737 ± 753 . 16. Gustafsson S , OÈ sterling M , Skurdal J , Schneider LD , Calles O . Macroinvertebrate colonization of a nature-like fishway: The effects of adding habitat heterogeneity . Ecol Eng . 2013 ; 61 : 345 ± 53 . 17. Winking C , Lorenz AW , Sures B , Hering D. Recolonisation patterns of benthic invertebrates: A field investigation of restored former sewage channels . Freshw Biol . 2014 ; 59 ( 9 ): 1932 ± 44 . 18. Baumgartner SD , Robinson CT . Short-term colonization dynamics of macroinvertebrates in restored channelized streams . Hydrobiologia . 2017 ; 784 ( 1 ): 321 ± 35 . 19. Bilton DT , Freeland JR , Okamura B. Dispersal in Freshwater Invertebrates . Annu Rev Ecol Syst . 2001 ; 32 : 159 ± 81 . 20. Gore JA . Benthic invertebrate colonization: source distance effects on community composition . Hydrobiologia . 1982 ; 94 : 183 ± 93 . 21. Lake PS , Bond N , Reich P. Linking ecological theory with stream restoration . Freshw Biol . 2007 ; 52 ( 4 ): 597 ± 615 . 22. Sundermann A , Stoll S , Haase P . River restoration success depends on the species pool of the immediate surroundings . Ecol Appl . 2011 ; 21 ( 6 ): 1962 ± 71 . PMID: 21939037 23. JaÈhnig SC , Lorenz AW , Hering D. Restoration effort, habitat mosaics, and macroinvertebratesÐdoes channel form determine community composition? Aquat Conserv Mar Freshw Ecosyst . 2009 ; 19 : 157 ± 69 . 24. Heino J . Environmental heterogeneity, dispersal mode, and co-occurrence in stream macroinvertebrates . Ecol Evol . 2013 ; 3 ( 2 ): 344 ± 55 . https://doi.org/10.1002/ece3.470 PMID: 23467653 25. Tachet , H. Invertebres d'Eau Douce . Systematique, Biologie, Ecologie. CNRS Editions; 2010 . 26. Van Leeuwen CA , Huig N , Van der Velde G , Van Alen T , Wagemaker C , Sherman CDH et al. How did this snail get here? Several dispersal vectors inferred for an aquatic invasive species . Freshw. Biol . 2013 ; 58 : 88 ± 99 . 27. Heino J , Melo AS , Siqueira T , Soininen J , Valanko S , Bini LM . Metacommunity organisation, spatial extent and dispersal in aquatic systems: Patterns, processes and prospects . Freshw Biol . 2015 ; 60 ( 5 ): 845 ± 69 . 28. Li F , Sundermann A , Stoll S , Haase P. A newly developed dispersal metric indicates the succession of benthic invertebrates in restored rivers . Sci Total Environ . 2016 ; 569 ± 570 : 1570 ±8. https://doi.org/10. 1016/j.scitotenv. 2016 . 06 .251 PMID: 27443455 29. Peckarsky BL . Colonization of natural substrata by stream benthos . Can. J. Fish. Aquat. Sci . 1986 ; 43 : 700 ± 709 . 30. MigueÂlez D , Valladares LF . Seasonal dispersal of water beetles (Coleoptera) in an agricultural landscape: a study using Moericke traps in northwest Spain . Ann Entomol Fr. 2008 ; 44 ( 3 ): 317 ± 26 . 31. Maguire B. The passive dispersal of small aquatic organisms and their colonization of isolated bodies of water . Ecol. Monogr . 1963 ; 33 : 161 ± 185 32. Bullock JM , Clarke RT. Long distance seed dispersal by wind: measuring and modelling the tail of the curve . Oecologia . 2000 ; 124 : 506 ±521 https://doi.org/10.1007/PL00008876 PMID: 28308389 33. Havel JE , Shurin JB . Mechanisms, effects, and scales of dispersal in freshwater plankton . Limnol. Oceanogr . 2004 ; 49 : 1229 ± 1238 . 34. HorvaÂth Z , Vad CF , Ptacnik R . Wind dispersal results in a gradient of dispersal limitation and environmental match among discrete aquatic habitats . Ecography . 2016 ; 39 : 726 ± 32 . https://doi.org/10.1111/ ecog.01685 PMID: 28529408 35. Lancaster J , Hildrew AG . Characterizing in-stream Flow Refugia . Can. J. Fish. Aquat. Sci . 1993 ; 50 : 1663 ± 1675 . 36. Lancaster J , Hildrew AG . Flow refugia and the microdistribution of lotic macroinvertebrates . J North Am Benthol Soc . 1993 ; 12 ( 4 ): 385 ± 93 . Winegardner AK , Jones BK , Ng ISY , Siqueira T , & Cottenie K. The terminology of metacommunity ecology . Trends in Ecology and Evolution . 2012 ; 27 ( 5 ), 253 ± 254 . https://doi.org/10.1016/j.tree. 2012 . 01 . 007 PMID: 22325446 38. Haase P , Murray-Bligh J , Lohse S , Pauls S , Sundermann A , Gunn R , et al. Assessing the impact of errors in sorting and identifying macroinvertebrate samples . Hydrobiologia . 2006 ; 566 : 505 ± 21 . 39. Schmidt-Kloiber A , Hering D. www . freshwaterecology.infoÐversion 7 .0, accessed in 2016/ 2017 , an online tool that unifies, standardizes and codifies more than 20,000 European freshwater organisms and their ecological preferences . Ecological Indicators 2015 ; 53 : 271 ± 282 . 40. Verdonschot PFM . Ecological characterization of surface waters in the province of Overijssel (The Netherlands) . University of Wageningen Press; 1990 . 41. Kenkel NC , OrloÂci L. Applying metric and nonmetric multidimensional scaling to ecological studies: some new results . Ecology . 1986 ; 67 : 919 ± 928 . 42. McArdle BH , Anderson MJ . Fitting multivariate models to community data: a comment on distancebased redundancy analysis . Ecology . 2001 ; 82 : 290 ± 297 . 43. Miller MF , Labandeira CC . Slow crawl across the salinity divide: Delayed colonization of freshwater ecosystems by invertebrates . GSA Today . 2002 ; 12 ( 12 ): 4 ± 10 . 44. Kappes H , Haase P . Slow, but steady: dispersal of freshwater molluscs . Aquat Sci . 2012 ; 74 :1± 14 . 45. Parkyn SM , Smith BJ . Dispersal Constraints for Stream Invertebrates: Setting Realistic Timescales for Biodiversity Restoration . Environ Manage . 2011 ; 3 : 602 ± 614 46. Milner AM , Robertson AL , Monaghan KA , Veal AJ , Flory EA . Colonization and development of an Alaskan stream community over 28 years . Front Ecol Environ . 2008 ; 6 ( 8 ): 413 ± 9 . 47. Haase P , Hering D , JaÈhnig SC , Lorenz AW , Sundermann A . The impact of hydromorphological restoration on river ecological status: a comparison of fish, benthic invertebrates, and macrophytes . Hydrobiologia . 2013 ; 704 : 475 ± 488 . 48. Tonkin JD , Stoll S , Sundermann A , Haase P . Dispersal distance and the pool of taxa, but not barriers, determine the colonization of restored river reaches by benthic invertebrates . Freshw. Biol . 2014 ; 59 : 1843 ± 1855 . 49. Schultheis AS , Weigt LA , Hendricks AC . Gene flow, dispersal, and nested clade analysis among populations of the stonefly Peltoperla tarteri in the Southern Appalachians . Mol. Ecol . 2002 ; 11 : 317 ± 327 . PMID: 11928706 50. Hughes JM. Constraints on recovery: using molecular methods to study connectivity of aquatic biota in rivers and streams . Freshw. Biol . 2007 ; 52 : 616 ± 631 51. Taylor BW . An example of long distance dispersal . Ecology . 1954 ; 35 : 569 ± 572 , 52. Mackay RJ . Colonization by Lotic Macroinvertebrates: A Review of Processes and Patterns . Can J Fish Aquat Sci . 1992 ; 49 ( 3 ): 617 ± 28 . 53. Allan JD . Stream Ecology: Structure and Function of Running Waters. Chapman and Hall , 1995 54. Tilman D. Constraints and tradeoffs: Toward a Predictive Theory of Competition and Succession . Oikos. 1990 ; 58 ( 1 ):3± 15 . 55. Tilman D. Competition and biodiversity in spatially structured habitats . Ecology . 1994 ; 75 ( 1 ):2± 16 . 56. Leps M , Sundermann A , Tonkin JD , Lorenz AW , Haase P . Time is no healer: Increasing restoration age does not lead to improved benthic invertebrate communities in restored river reaches . Sci Total Environ . 2016 ; 557 ± 558 : 722 ± 32 . https://doi.org/10.1016/j.scitotenv. 2016 . 03 .120 PMID: 27046138 57. Gallardo B , DoleÂdec S , Paillex A , Arscott DB , Sheldon F , Zilli F , et al. Response of benthic macroinvertebrates to gradients in hydrological connectivity: A comparison of temperate, subtropical, Mediterranean and semiarid river floodplains . Freshw Biol . 2014 ; 59 ( 3 ): 630 ± 48 . 58. Hynes HBN . Biology of Plecoptera. An. Rev. of Entomology . 1976 ; 15 : 25 ± 42 . 59. Jacobi G , Cary S. Winter stoneflies (Plecoptera) in seasonal habitats in New Mexico, USA . J. North Am. Benthol. Soc . 1996 ; 15 : 690 ± 699 . 60. LoÂpez-RodrÂõguez MJ , Tierno de Figueroa JM , Fenoglio S , Alba-Tercedor J . Life strategies of 3 Perlodidae species (Plecoptera) in a Mediterranean seasonal stream in southern Europe . J. North Am. Benthol. Soc . 2009 ; 28 : 611 ± 625 61. Verdonschot PFM , Besse-Lototskaya AA , Dekkers DBM , Verdonschot RCM . Mobility of lowland stream Trichoptera under experimental habitat and flow conditions . Limnologica . 2012 ; 42 ( 3 ): 227 ± 234 . 62. Knispel S , Castella E. Disruption of a longitudinal pattern in environmental factors and benthic fauna by a glacial tributary . Freshw. Biol . 2003 ; 48 : 604 ± 618 . Buffagni A. Distribution and Ecological Preferences of European Freshwater Organisms. Volume 3Ð Ephemeroptera . Edited by Schmidt-Kloiber, A . and D. HeringÐPensoft Publishers; 2009 . Graf W. Distribution and Ecological Preferences of European Freshwater Organisms . Volume 2ÐPlecoptera. Edited by Schmidt-Kloiber A . and D. HeringÐPensoft Publishers; 2008 Townsend CR , Scarsbrook MR , Dole dec S. Quantifying Disturbance in Streams: Alternative Measures of Disturbance in Relation to Macroinvertebrate Species Traits and Species Richness . J North Am Benthol Soc . 1997 ; 16 ( 3 ): 531 ± 44 . Verberk WCEP , Siepel H , Esselink H . Life-history strategies in freshwater macroinvertebrates . Freshw Biol . 2008 ; 53 ( 9 ): 1722 ± 38 . 68. Jacobsen D , Cressa C , Mathooko JM , Dudgeon D. Macroinvertebrates : composition, life histories and production . In: Dudgeon D , Tropical stream ecology . Academic Press. 2008 ; pp 65±105 Palmer MA, Menninger HL , Bernhardt E. River restoration, habitat heterogeneity and biodiversity: a failure of theory or practice? Freshw Biol . 2010 ; 55 : 205 ± 22 . Westveer JJ , Verdonschot PFM , Verdonschot RCM . Substrate homogenization affects survival and fitness in the lowland stream caddisflies Micropterna sequax and Potamophylax rotundipennis: a mesocosm experiment . Freshw Sci . 2017 ; 36 : 585 ± 94 . J Anim Ecol . 2000 ; 69 : 442 ± 57 .


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0197182&type=printable

Judith J. Westveer, Harm G. van der Geest, E. Emiel van Loon, Piet F. M. Verdonschot. Connectivity and seasonality cause rapid taxonomic and functional trait succession within an invertebrate community after stream restoration, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0197182