Plant Species Over-Occupancy Indicates River Valleys Are Natural Corridors for Migration

River valleys are considered natural corridors for migration of plant species, however, there is a lack of studies conrming higher colonisation rates of plant species in these areas. We compare plant species-richness between ancient and recent forests, and those located in a river valley with those located in the surrounding landscape. We hypothesise that, close to a river, higher plant species-richness will be associated with recent forests thus indicating a higher colonisation rate. The study area includes part of the Elbe River Valley and its surrounding landscape in the Czech Republic. We sampled an equal number of recent and ancient forests but lying at different distances from the river. We used generalised linear models to test the effect of distance from the river in dependence upon forest continuity (recent/ancient forest) on two plant species-richness categories, i.e. richness of forest species and overall species richness. In the surrounding landscape, higher richness of forest species was associated with ancient forests, whereas overall species richness was comparable. In the river valley, richness of forest species as well as overall species richness was higher in the recent forests. Recent forests in the river valley were more saturated by plant species than those in the surrounding landscape, indicating that in the river valley the colonisation rate of plant species is higher. These results conrm the importance of river valleys as natural corridors for migration of plant species.


Introduction
It is well known that the plant species diversity is markedly different in a river valley compared to in the surrounding at landscape (Blažková 1964;Naiman et al. 1993;Renöfält et al. 2005;Zelený and Chytrý 2007;Nobis et al. 2016). Hence, river valleys are frequently identi ed as regional hotspots of biodiversity and are also regarded as natural corridors for plant species migration (Jeník and Slavíková 1964;Tockner and Ward 1999;Ward et al. 2001;Saarela et al. 2017). However, we lack studies which estimate whether the colonisation rate of plant species is actually higher in a river valley, than in the surrounding landscape.
The linear connectivity of a river valley allows plant species to migrate in both downstream and upstream directions and also along interconnecting gullies and so colonise suitable habitats (Gallé et al. 1995;Ward et al. 2002;Czochański and Wiśniewski 2018). Here, the colonisation rate decreases as the target habitats become more distant from the riverbank which serves as a source of diaspores. Correspondingly, the decrease in habitat connectivity outside of a river valley has a negative effect on colonisation rate (Yao et al. 1999;Renöfält et al. 2005;Vellend 2016;Nobis et al. 2016;Nobis et al. 2017). For example, the oak-hornbeam forests of central Europe offer highly connected habitats along the slopes of river valleys, whereas outside of these a patchy agricultural landscape dominates (Chytrý 2013). Therefore, outside of a river valley, compared to colonisation, both the distance from the riverbank and also habitat fragmentation increase the importance of local community processes including competition and ecological drift. The ecological drift and competition at any particular site increase the chances of species extinction. Thus, the balance between the rate of colonisation and that of extinction determines the composition and species richness of a local plant community (Vellend 2016).
Whether or not a plant species reaches a remote habitat, depends on its colonisation ability. For example, both anemochorous and zoochorous species produce many seeds which can move over large distances, however, colonisation success is also affected by the seedling establishment conditions. Since prehistory in Europe, forest cover has gradually decreased accompanied by a corresponding increase in cropland cover (Bellemare et al. 2002;Pongratz et al. 2008). Since the intensi cation of farming practices in the 20th century, many croplands formerly established on forest soils have been abandoned and the land has returned to forest (Mather et al. 1999). Therefore, two types of forests can be found in European landscapes that differ in respect to their times of origin. Ancient forests are associated with land has enjoyed continuous forest land use, while recent forests are associated with land that has suffered temporarily use as pasture, or for arable cropping or for horticulture (Koerner et al. 1997;Webb 1998). River valleys have been less exposed to these human-induced changes due to their rugged topography, although evidence of some traditional farming practices can still be seen in the more accessible locations (Douda 2010). A good number of studies have shown that higher plant species richness is associated with ancient forests due to their long and continuous history (Peterken and Game 1984;Dzwonko and Gawronsky 1992;Hermy et al. 1999;Wulf 2003;Sciama et al. 2009;Matuszkiewicz et al. 2013).
Furthermore, particular groups of species are related to ancient forests, namely ancient forest species that have been found to serve well as indicators of forest communities and these are characterised by high diversity Douda 2010;Stefańska-Krzaczek et al. 2016). The recent forests are dependent on diaspore input from ancient forests, i.e. colonisation success increases in proximity to an ancient forest, and recent forests continuous with ancient forests may also be colonised by herbaceous species (Dzwonko 1993;Wulf and Heinken 2008;Orczewska 2010;Orczewska and Fernes 2011).
Disentangling the importance of the processes that drive the plant species richness of riverine communities, requires information on plant species colonisation. In previous studies, a number of authors have described the patterns of plant species richness along a gradient starting in a river valley and continuing into the surrounding landscape (i.e., a riverine gradient, Blažková 1964;Ložek 1974;Chytrý and Tichý 1998;Renöfält et al. 2005;Zelený and Chytrý 2007;Nobis et al. 2016). Here we seek to provide an additional step in understanding how well plant species richness copes with the riverine phenomenon by recording the species saturation of forest communities along a riverine gradient. We compare plant species richness between ancient and recent forests where ancient forests are considered to be areas of continuously forested land, at least since the second half of the 19th century, whereas recent forests are those that have originated after that time. This method allows us to estimate how successfully recent forests are colonised at different distances from the river and to associate the riverine gradient with the importance of the colonisation processes. The ancient forests provide us information on plant species richness determined by the long-term balance between colonisation, extinction and environmental conditions. Recent forests inform us as to whether the plant species richness here has reached a value close to the natural state in an ancient forest due to the rapid spread of plant species or whether the dispersal of plant species is more limited, so that the recent forests are less saturated than the ancient forests. Use of Ellenberg indicator values (EIVs) helps us to understand patterns of plant species composition in respect to changing environmental conditions along the riverine gradient.
Assessing the distribution of plant dispersal vectors in forests located in, and outside, the river valley further contributes to a better explanation of the colonisation processes. Speci cally, we ask: How do plant species-richness and the prevailing dispersal mode differ between ancient and recent forests, with respect to the distance from a river? We would expect that plant species-richness will be more similar between ancient and recent forests if these are located close to a river, due to their connectivity with the river corridor that will enhance the colonisation rate of plant species. Due to higher connectivity between recent and ancient forests, we also expect higher proportions of species with limited dispersal, such as myrmecochorous forest species, in recent forests in a river valley. Conversely, we would expect that a recent forest far from a river, will be less saturated with plant species, compared with an ancient forest, due to more limited dispersal. Therefore, if present, forest species will be represented dominantly by anemochorous species.

Study area
The study was conducted about halfway along the course of Elbe River located in the northern part of the Czech Republic (14°09'25''-14°09'47'' E, 50°41'53''-50°42'29'' N). The Elbe River Valley is considered one of the most important natural corridors for species migration in Europe. The study area covers ca. 10 km 2 and includes the western slope of the river valley and its surrounding landscape (Fig. 1). The altitude ranges from 153 m to 340 m above mean sea level. The climate is subcontinental with annual mean temperature ca. 9°C. The mean annual precipitation ranges from 450 to 800 mm (Bína and Demek 2012). The whole area is a part of a protected landscape České středohoří (IUCN category V) and comprises many primary but also secondary ecosystems including ancient and recent forests. It therefore represents a suitable model system for studying patterns of plant species richness in forests along a riverine gradient.

Sampling
To select the vegetation plots for random-strati ed sampling, the study area was overlaid with a square grid with cells of 200 × 200 m. For each cell, forest continuity (ancient/recent forest) was determined based on a comparison of current with historical maps, speci cally the third military mapping of Austria-Hungary from 1877 − 1880 on a scale of 1:75 000 (http://oldmaps.geolab.cz/). The grid cells were grouped into ve 600-m-wide strips, corresponding to their different distances from the river. In each strip, the cells of ancient/recent forest were numbered. Four of these cells were then selected at random and vegetation plots were established at their centres. In cases where ancient/recent forest was absent in a cell centre, the vegetation plot was established in the nearest location in which the forest occurred ( Fig. 1). In this way, eight vegetation plots (four ancient forest plots and four recent forest plots) were established in each 600 m strip. Data was collected in July 2016. A list of vascular plants (shrubs and trees) with species cover was recorded for each plot using the nine-step ordinal Braun-Blanquet scale (Braun-Blanquet 1964;Podani 1984). The standard plot size used for sampling was 400 m 2 .

Delimitation of variables
Forest species were selected from the overall dataset based on a study by Chytrý and Tichý (2003), which categorised the a nity of plant species to forest habitats. The number of species belonging to the forest species category (richness of forest species) and the number of all species (overall species richness) were calculated for each vegetation plot. The full list of species recorded, with forest species highlighted, is provided in Supplement Table 1. The nomenclature follows Kaplan et al. (2019). Table 1 The effects of forest continuity, distance from the river and their interaction on plant species richness, dispersal and Ellenberg indicator values is shown by generalised linear models. Light ns ns ns *** = P < 0.001, ** = P < 0.01, * = P < 0.05, ns = non-signi cant The shortest distance from each plot centre to the river was calculated. This distance ranged from a minimum of 56 m (for a plot lying adjacent to the river) to a maximum of 2840 m (for the most distant plot).
For each species, dispersal mode was assigned according to the LEDA database (Kleyer et al. 2008) and the number of species within each mode (autochory, anemochory, endozoochory, epizoochory, myrmecochory) was then calculated for each vegetation plot. For each vegetation plot, the mean EIVs for nutrients, soil reaction, moisture and light were also calculated and weighted on species cover (Ellenberg et al. 1992).
It is generally well known that the plant species-richness of a recent forest is affected by proximity to ancient forest. Speci cally, the closer a recent forest is to an ancient forest, the higher is the number of plant species present (Dzwonko 1993). Therefore, it is likely that the plant species-richness of a recent forest will be affected not only by distance from a river but also by distance from the nearest ancient forest. To control this effect, the distance from each recent forest plot to the nearest ancient forest plot was calculated (distances ranged from a minimum of 179 m to a maximum of 968 m) and its effect on plant species-richness in the recent forest was also assessed.

Statistical analyses
To estimate how successfully recent forests are colonised at different distances from the river, the effect of distance from the river in dependence upon forest continuity (recent/ancient forest) on richness of forest species/overall species richness was tested using generalised linear models with Poisson response distributions.
Several additional analyses were carried out to explain the observed patterns: First, to relate the colonisation rate to the dispersal traits of recorded species, the effect of distance from the river in dependence upon forest continuity and on richness of autochorous, anemochorous, endozoochorous, epizoochorous and myrmecochorous species was tested using generalised linear models with Poisson response distributions.
Second, to assess the change in environmental conditions along the riverine gradient, the effect of distance from the river in dependence upon forest continuity on distributions of the EIVs for nutrients, soil reaction, moisture and light were tested using generalised liner models. In this case, the Gaussian response distribution was used, as the original EIVs were converted on weighted means.
Third, to assess whether the proximity to an ancient forest contributes to higher plant species-richness of recent forests, the effect of the interaction between the distance from the river and the distance from an ancient forest (and the separate effects) on richness of forest species/overall species richness in recent forests were tested using generalised linear models with Poisson response distributions. Similar analyses were also used to test the effects of the putative interactions of the richness of autochorous, anemochorous, endozoochorous, epizoochorous and myrmecochorous species in recent forests.
Fourth, detrended correspondence analysis (DCA) was used to assess the overall pattern of variation in species composition. The analyses were carried out using the decorana function in the R package Vegan (Oksanen et al. 2020). To interpret the DCA results in terms of environmental gradients, the signi cant Ellenberg indicator values were then plotted onto a DCA ordination diagram as supplementary variables using the env t function from the same R package. Env t provides a vector for each variable that maximises correlation between that variable and the projections of ordination points onto that vector. Signi cance was assessed using permutation of variables. All statistical analyses were carried out in the statistical program R (R Development Core Team 2019).

Results
In general, signi cant positive correlations were found between distance from the river and each speciesrichness category (i.e., overall and forest-species richness). Higher numbers of species were associated with plots located far from the river (outside the river valley), whereas lower number of species persisted on plots located closer to the river (in the river valley) (Table 1; Fig. 2). Both species-richness categories differed signi cantly for ancient and recent forests (i.e., forest continuity) ( Table 1; Fig. 2).
Higher richness of forest species was associated with recent forests in the river valley, whereas outside of the river valley higher richness of forest species was associated with ancient forests. Speci cally, the number of species in recent forests located closer to the river (0.056-1.06 km) was 1.6-times higher compared with ancient forests at the same distances. On the other hand, the number of species in recent forests located far from the river (1.07-2.84 km) was 1.14-times lower (Fig. 2a).
In contrast, higher overall species-richness was associated with recent forests in the river valley as well as those outside the river valley, however, the difference in the numbers of species between ancient and recent forest decreased gradually as distances increased. Speci cally, the numbers of species in recent forests located close to the river (0.056-1.06 km) were 1.8-times higher compared with in ancient forests at the same distance, whereas the numbers of species in recent forests located far from the river were comparable with the numbers of species in ancient forests (Fig. 2b).
A signi cant positive correlation was found between distance from the river and the dispersal mode. However, we recorded signi cantly higher numbers of autochorous, endozoochorous and myrmecochorous species associated with recent forests along the entire length of the riverine gradient. In cases of endozoochory and myrmecochory, the numbers of species were higher in plots located close to the river (Table 1; Fig. 3).
Signi cant positive correlations with distance from the river were found only for the EIVs for nutrients. However, these EIVs did not differ signi cantly between ancient and recent forests (Table 1; Fig. 4). This result was also supported by the DCA ordination, indicating that the riverine gradient is correlated with the EIV for nutrients, which explains the most important part of the variability in the species-composition data ( Table 2; Fig. 5). In recent forests, the effect of interaction between the distance from the river and the distance from the ancient forest, was not signi cant, either for richness of forest species or for overall species richness. Similar results were also found for the richness of autochorous, anemochorous, endozoochorous, epizoochorous and myrmecochorous species in recent forests (see Supplement Table 2).

Discussion
Our results support the view that colonisation patterns of plant species strongly re ect the riverine gradient. Our results show a higher saturation of plant species in recent forests and this nding is consistent with our predictions that recent forests in river valleys are exposed to higher migration. Conversely, the proportion of plant species in the recent forests decreased with increasing distance from the river valley. The results indicate that migration processes are more important in river valleys than away from them, con rming the view that river valleys are important migration corridors compared to the surrounding at landscapes.
The lower richness of forest species found at some distances from the river valley is consistent with other studies which found evidence for low saturation in recent forests (Peterken and Game 1984;Dzwonko and Gawronsky 1992;Hermy et al. 1999;Wulf 2003;Sciama et al. 2009;Matuszkiewicz et al. 2013). In most of these studies, the dearth of forest species is interpreted as being a consequence of limited dispersal or of unsuitable recruitment conditions (Ehrlén and Eriksson 2000;Ehrlén et al. 2006).
According to the dispersal limitation hypothesis (1), forest species that are unable to rapidly colonise recently developed forests, are primarily associated with ancient forests Stefańska-Krzaczek et al. 2016). The reason for this is that these species are adapted to small scale disturbance regimes in temperate forests, where species usually spread only over quite short distances as they colonise canopy gaps Bellemare et al. 2002). For example, many forest species are dispersed locally by ants, some species of which are very sensitive to habitat fragmentation and rarely travel through treeless areas , Mitchell et al. 2002, Douda 2010). The observed patterns in richness of myrmecochorous species support this view as the values between ancient and recent forest are more similar in the river valley (Fig. 3e).
Meanwhile, the recruitment limitation hypothesis (2) states that in a recent forest site, the altered soil characteristics resulting from the former agricultural use are retained. These include, increased pH, base saturation and a greater availability of phosphate and nitrogen (Dzwonko and Loster 1994;Richter et al. 1994;Koerner et al. 1997;Honnay et al. 1999). These altered soil conditions restrict the recruitment and population growth of forest species, even though they are able to colonise the site (Hurtt and Pacala 1995).
The second hypothesis does not seem able to account satisfactorily for the patterns of species saturation we observed. We would also expect a poor saturation of recent forests in the river valley if the recruitment limitation due to altered soil conditions. Instead, the contrast between the levels of forest species in the river valley and in the surrounding landscape suggests that dispersal limitation occurs outside of the river valley, with easier dispersal within the river valley. Contrary to the recruitment limitation hypothesis, we did not nd differences in the EIVs, which might indicate the different levels of nutrients, light availability or moisture found between recent and ancient forests.
The patterns of overall species richness also supports easier dispersal within the river valley, as higher numbers of species are associated with recent forests in the river valley. However, outside the river valley, the effects of dispersal limitation are not so clear, as the values between ancient and recent forests are highly comparable. It should be recognised that in this study we include all plant species recorded in the forests, not just the species belonging to the forest-species category. In other words, plant species could still persist in recent forests that occurred in a previously open habitat.
The reverse species richness patterns were observed in the river valley where more plant species occurred in recent forests than in ancient forests. The over-occupancy of recent forests could be explained by a combination of a fast colonisation rate and a slow extinction rate in the river valley. It is well known that competition among species can be relatively weak early in forest succession during the initial establishment of plant populations, whereas the importance of interspeci c competition for light and nutrients grows as a population becomes more established (Slatkin 1974;Goldberg and Barton 1992). At the start of a succession (i.e., in recent forests), a higher colonisation rate leads to an accumulation of species, whereas in later stages of succession (i.e., in ancient forests) the increasing role of competition limits the successful recruitment of new colonisers and leads to elimination of some species from the community. The number of species in ancient forests is determined by the long-term balance between species carrying capacity, the colonisation rate and the extinction rate (Connell and Slatyer 1977;Tilman 1990;Groenendael 1998, Kneitel andChase 2004). We expect low levels of interspeci c competition at the start of a succession leading to a rapid recruitment of plant species in recent forests that outnumber those in ancient forests. That this occurred only in the river valley, indicating that colonisation rate is higher there than it is in the surrounding landscape. Moreover, a higher colonisation rate in the river valley should also be re ected in a higher proportion of plant species with various dispersal modes in recent forests. This was con rmed at least in the cases of the autochorous, endozoochorous and myrmecochorous species (Table 1; Fig. 3).
In the future, we would expect a gradual decrease in plant species-richness in the recent forests as it declines to the level observed in ancient forests. However, population decline in perennial species, and especially in forest species, is a very slow process, so the gradual extinction of many forest species is much delayed (Tilman et al. 1994;Vellend et al. 2006).
River valleys may favour species migration for a number of reasons. Among these are the increased nescale habitat heterogeneity in such a variable topography, both the mesoclimate and the geological conditions. The topography of a river valley is characterised by steep slopes with exposed rocky outcrops, at alluvial areas and variable landform shapes (Jeník and Slavíková 1964). Also, frequent temperature inversions result from the valley shape, and these increase the contrast between a cold and wet valley bottom and a dry, more continental upper valley edge (Chytrý and Tichý 1998). Similarly, north facing and south facing slopes represent suitable habitats for shade-demanding and light-demanding species so enabling the close coexistence of these species within a river valley (Blažková 1964;Chytrý and Tichý 1998;Zelený and Chytrý 2007). All these abiotic variables create a set of suitable environmental conditions, which may host a large spectrum of species with quite different habitat requirements. This nes-scale heterogeneity may enhance the connectivity of a river valley for multiple species, rather than presenting a barrier to their migration.
Contrary to the expectation based on numerous studies (Blažková 1964;Ložek 1974;Chytrý and Tichý 1998;Renöfält et al. 2005;Zelený and Chytrý 2007;Nobis et al. 2016), we found the distance from the river was positively associated with plant species-richness. The likely explanation for this is that environmental conditions are more favourable here, so more plant species can establish successfully. This explanation is supported by the increasing EIVs for nutrients with distance from the river as shown in Fig. 4. Moreover, higher resource availability in the surrounding landscape is also suggested by the DCA ordination in which the nutrient and water availabilities are also highly correlated with distance from the river (Fig. 5).
Higher productivity in the surrounding landscape appears to be associated mainly with geological conditions. The geological base for our whole study area comprised two main geological formations: volcanic rocks (predominant in the river valley) and clay sediments (predominant in the surrounding landscape). These have fundamentally different in water holding properties (see Supplement Fig. 1). Clay sediments retain water in the upper levels of the soil horizon and thus create the appropriate moisture conditions for plant species during the larger part of the growing season. Conversely, highly waterpermeable volcanic rocks can lose water earlier in the growing season, and this can lead to severe water stress and to a decline in species richness (Cannon 1971;Lite et al. 2005). In general, forest species are adapted to deep, moist soils rather than to the more extreme habitats of exposed rocky outcrops, where dry-grassland species increase their proportion (Chytrý and Tichý 2003;Douda et al. 2017).
Many studies mention river valleys as important migration corridors for plant species (Jeník and Slavíková 1964;Tockner and Ward 1999;Ward et al. 2001;Saarela et al. 2017). However, none of these compare the colonisation rates in river valleys with those in the surrounding landscape. This also includes studies of invasive plant species for which dispersal along rivers is fundamental to their rapid spread outside of their native range (Hood and Naiman 2000;Richardson et al. 2007;Osawa et al. 2013;Aronson et al. 2017). In this study, we used an original approach to compare community saturation by plant species along a riverine gradient. Ancient forests served as a control where species richness is determined by the balance between species carrying capacity, colonisation rate and extinction rate. Under these conditions, the proportion of plant species in recent forests compared with that in ancient forests becomes a measure of colonisation rate. As predicted, the recent forests in the river valley were saturated with plant species as they had been exposed to a high colonisation rate. This result con rms the strong importance of river valleys as natural corridors for species migration. River valleys may also serve as refuges for endangered plant species, with rapid dispersal among river sites helping prevent their extinction. Unfortunately, recent large constructions for irrigation, for hydropower and reservoirs for urban water supply have cause major fragmentation of these river-valley habitats . We hope that our results will serve to underpin conservation efforts and for the prioritisation of river valleys in conservation strategies.

Declarations
Funding This study was supported by the project -The understanding of long-term sustainability of ecosystem services of European oodplain forests (LTC18045) of the Ministry of Education, Youth and Sports of the Czech Republic, and Vegetation changes and biodiversity conservation of oodplain forests situated along the regulated rivers (2020B0035) of the Faculty of Environmental Sciences, Czech University of Life Sciences Prague.  Study area with 40 vegetation plots distinguished as being either recent or ancient forest. Note: The designations employed and the presentation of the material on this map do not imply the expression of any opinion whatsoever on the part of Research Square concerning the legal status of any country, territory, city or area or of its authorities, or concerning the delimitation of its frontiers or boundaries. This map has been provided by the authors. The effect of distance from the river on a) richness of forest species, b) overall species richness distinguished as being either recent or ancient forest.

Figure 4
The effect of distance from the river on Ellenberg indicator values for nutrients.