An uphill battle? The elevational distribution of alien plant species along rivers and roads in the Austrian Alps

Ever-increasing international trade and anthropogenic activity has led to the relocation of thousands of plant species worldwide. So far, the harsh climate of the European Alps historically has restricted the establishment of alien plants. However, new opportunities created by rising temperatures and increasing human activity might allow alien plants to spread further upwards. Here, the distribution of alien plants along an altitudinal gradient in two Austrian valleys is analyzed. Specifically, the distribution along two contrasting corridors (roads, rivers) and the spread of alien plants into adjacent habitats is examined. Following the MIREN sampling protocol, 20 transects composed of three plots along each river and main road, were established in each study region. Plant species cover and a range of site-specific factors were collected. In total, 641 plant species were recorded, of which 20 were alien. Alien species richness along roads was slightly higher compared to rivers, and the composition of the alien flora differed markedly between roads and rivers. Further, alien plant species richness decreases with distance to roads and rivers (indicating that adjacent habitats are less invaded), as well as with increasing elevation. Mowing along roadsides resulted in lower alien plant species cover, but higher alien plant species richness. Finally, compositional dissimilarity between sites showed that elevation, proximity of a plot to a river or road, and alien plant cover are important factors for higher dissimilarity. This study demonstrates that both natural (rivers) and manmade (roads) corridors play an essential role in the upward spread of different alien plants in mountains.


introduction
Rising intensity of international trade and anthropogenic activity increases the humanmediated transport of plant species to regions beyond their native range. The introduction of alien plant species accelerated over the past centuries with the expansion of the global trade network (Hulme 2009;Seebens et al. 2015Seebens et al. , 2017Amano et al. 2016;Kueffer 2017;van Kleunen et al. 2018). Many plants that are transported by anthropogenic means face difficulties in reproducing and finally establishing themselves successfully in their new location due to climatic or other environmental constraints (Caley et al. 2007;Stotz et al. 2016). Nevertheless, globally at least 13,168 alien plant species managed to establish (i.e. form viable populations) in their new environments , and there is no sign of a saturation of alien plant species accumulation (Seebens et al. 2017).
Mountains are generally subject to lower levels of plant invasion in comparison to lowland areas, and invasions into high elevations are rare (Lembrechts et al. 2016a;Haider et al. 2018). However, there is strong evidence that mountain systems are becoming more susceptible to plant invasions (Pauchard et al. 2009;Kueffer et al. 2013;Pauchard et al. 2016) particularly due to climate change (Walther et al. 2009;Bellard et al. 2018). Increasing temperatures may translate into greater habitat availability for alien species (Lenoir et al. 2008;Petitpierre et al. 2015;Carboni et al. 2018). At the same time, climate change alters the distribution of native biota , thereby making mountains more vulnerable to invasions (Diez et al. 2012). For instance, it has been shown that several alien plant species in the Swiss Alps are occurring at higher elevations compared to a few decades ago (Becker et al. 2005). In addition, many newly introduced alien species undergo a 'lag phase', i.e. range expansion, that often follows with substantial delay after initial introduction (Aikio et al. 2010;Rouget et al. 2016;Alexander et al. 2018). This implies that many alien species are still expanding, and that the equilibrium distribution of alien plant species will be considerably larger than their current one.
Corridors connecting different habitats (e.g. rivers, roads) play a decisive role in the spread of alien plants (Pattison et al. 2017;Rauschert et al. 2017;Follak et al. 2018). So far, studies in mountains have mainly focused on road networks as facilitators for the dispersal of alien plant species . Vehicular traffic plays a critical role in the rapid dispersal of seeds (von der Lippe and Kowarik 2007; Lemke et al. 2019). Frequent disturbances along roadsides create opportunities for the establishment of alien species (Hansen and Clevenger 2005), and these roadside habitats allow the spread of alien plants into higher elevations Lembrechts et al. 2016a;Seipel et al. 2016). From there, species might subsequently spread into adjacent less-disturbed habitats . However, while it is well-understood that European riverine habitats show high levels of plant invasion (Hejda et al. 2009;Kalusová et al. 2013), comparatively less work has been done on studying the role of rivers rather than of roads for the spread of alien plants in mountains (but see Siniscalco et al. 2011 andBarni et al. 2012).
Riverbanks typically harbor many native plant species of conservation interest. These experience increased competition with the arrival of alien plant species in these habitats (Naiman and Décamps 1997). Studies in lowland areas have shown that the constant movement of rivers provides opportunities for rapid seed dispersal of alien plant species into previously uninvaded areas (Pyšek and Prach 1994;Aronson et al. 2017). River corridors do not solely facilitate the spread of alien plants downstream but have also been observed to do so upstream (Wang et al. 2011;Osawa et al. 2013). Changes in riverside vegetation by alien plants can lead to erosion and may ultimately disturb ecosystem services (Richardson et al. 2007;Greenwood et al. 2018).
In the face of increasing numbers of alien species, a series of control programs have been enacted in Austria to prevent the establishment, hinder further spread, or to eradicate (invasive) alien species (Schiffleithner and Essl 2016). One of these methods is mowing, which is a common practice at roadsides of the European Alps for both conservation efforts and road safety (ASFINAG 2019). However, mowing can have diverging effects on alien and native plants depending on the individual species and the timing of the mowing event, indicating that this method should be considered on a case-by-case basis (Vitalos and Karrer 2009;Kettenring and Adams 2011;Song et al. 2018).
The aim of this research was to analyze and compare the role of natural (rivers) and man-made (roads) corridors for the spread of alien plant species in mountain valleys. Specifically, the following research questions were addressed: 1) What is the elevational distribution of alien plant species along rivers and roads? 2) What is the effect of roadside mowing on the distribution of alien plant species richness and cover? 3) Do plant communities along rivers and roads, and with different proximities to these corridors, show a dissimilarity in composition? 4) Which factors, including the occurrence of alien plant species, affect plant community composition?

Study regions
Climate. In Austria, where 60% of the country has an elevation above 500 m, alien species are mostly confined to the lowlands and large river valleys in the Alps (Walter et al. 2005). For this study, we selected the main rivers and roads of two mid-size Austrian valleys as study regions: the Lech valley located in the Northern Alps, and the Isel valley located in the Central Alps (Fig. 1A, B, D). Both study regions represent rural, somewhat agricultural, landscapes dominated by forests, mountain pastures and grasslands, while the major rivers are still highly dynamic and in near-natural condition. The Lech valley has a predominantly oceanic climate and calcareous bedrock. The mean annual temperature is 7.3 °C and it has a mean annual precipitation of 1397 mm, near the town of Reutte at 842 m a.s.l. Further up-valley in Warth at 1478 m a.s.l., the mean annual temperature drops to 4.7 °C, with a mean annual precipitation of 1867 mm (ZAMG 2019). The Isel valley is dominated by siliceous bedrock. This valley sees a mean annual temperature of 8.5 °C and a mean annual precipitation of 944 mm, at an elevation of 661 m a.s.l. In the village of Virgen at 1212 m a.s.l., the mean annual temperature is 6.7 °C, with a mean annual precipitation of 850 mm (ZAMG 2019). Topography. The Lech river originates to the southwest of the village Lech at 1840 m a.s.l. and runs northeast, stretching 256 km in total before it flows into the Danube. The Austrian-German border (824 m a.s.l.) marks the lowest point included in this study, and the distance between this point and the source is approximately 90 km. The end of the riverbed marks the highest point included (1816 m a.s.l). The Lech river is characterized by the way it alternately flows through gorges and sections of broad river valleys with a braided riverbed (Müller and Bürger 1990). The main road mostly runs along the river at the bottom of the valley but ascends to slopes along the gorges.
The Isel river and its accompanying road have a length of approximately 57 km. Its main source is the Umbalkees glacier (2400 m a.s.l.) in the Hohe Tauern mountains. The Isel river flows into the Drau river near the city of Lienz (673 m a.s.l.). This point therefore marked the lowest elevation in the study. As high alpine elevations are unlikely to contain alien plant species (Becker et al. 2005;Pauchard et al. 2009;Alexander et al. 2016), the upmost transect was established near the Johannis mountain hut at 2121 m a.s.l. The rise in elevation along the course of the river is roughly exponential as is characteristic for glacier valleys (Montgomery 2002).

Sampling design
Transect layout. The study design along rivers and roads followed a similar approach to ensure that results were comparable. The underlying method was structured following the method designed by the Mountain Invasion Research Network (MIREN) (Seipel et al. 2012). Each transect had a T-shape and consisted of three plots, all sized 50 m x 2 m. The first plot (river/road plot) was placed parallel to the road or river. The second plot (intermediate plot) was positioned at a 90° angle adjacent to the center of the river/road plot. The third plot (interior plot) was set up next to the intermediate plot in the same orientation (Fig 2). Roadside plots were placed next to the unvegetated edge on the road verge. The road verge, in contrast to the road edge, was vegetated but showed signs of mechanically induced disturbance (e.g. road construction/ maintenance, overrunning). In practice, the road edge did not exceed a few cm in the Isel valley and not more than 60 cm in the Lech valley, while the road verge was up to 101 cm wide in the Isel valley and up to 158 cm wide in the Lech valley.
Mowing protocols along roadsides were designed and executed by the ASFINAG, the Austrian publicly-owned corporation responsible for road maintenance. A total of 14 transect locations along roads in the Isel valley and 15 in the Lech valley were subject to mowing. The use of similar mowing machines in both areas led to comparable and systematic mowing regimes (ASFINAG 2019). In all cases, it was clearly visible up to where the vegetation was mown ('mowing line', see picture in Suppl. material 2: Appendix 1). We placed the fourth plot (mowing line plot) next to the mowing line in the first strip of unmown vegetation, parallel to the road plot (Fig. 2B). On average, the distance between the road edge and the mowing line was 451 cm in the Isel valley, and 495 cm in the Lech valley.
Transect positioning. Transects locations were established using ArcGIS 10.2.2 (ESRI 2014), by dividing the total length of both the rivers and roads into four equal sections. These four sections were subsequently divided into five equal-elevation parts, each covering the same elevational range. Transects were placed in the center of each of the resulting 20 elevational bands, except for the first and last transect of each road and river that were placed at the lowest and highest points of the study region (Fig. 1C). The side of the river or road at which the transect was established was chosen at random using a binary random number generator, unless one side was inaccessible (e.g. Figure 2. Study design illustrating the layout of the river (A) and road (B) transects. All plots had a size of 50 m × 2 m. The measurement points for the coordinates of each plot are visualized by a red rhombus. A the river plot was positioned at the line where 50% of the surface of the riverbank was covered by herbaceous species. The riverbank width was the area which was covered in gravel and had very little or no vegetation B the road plot was positioned adjacent to the non-vegetated road edge at the road verge. The road verge was characterized by the first occurrence of vegetation. To quantify the effect of mowing, a plot was established directly beyond the mowing line (i.e. in the first unmown vegetation). The mowing line plot followed the same orientation as the road plot. due to steep cliffs). Transects were relocated when the river or road was surrounded by agricultural land or built-up areas. Transects were preferably placed at least 100 m away from secondary roads or streams. If these requirements were not met, the transect was moved to the closest location where they were fulfilled.
Plot coordinates were recorded in the field with a Garmin eTrex 10 GPS device and verified in ArcGIS. Elevation and location for river and road plots were measured in their center. For intermediate and interior plots, this was recorded at the end of the plot most distant from the river or road (Fig. 2). From these values, the average elevation of the intermediate and interior plots was subsequently calculated for the statistical analysis. The width of the river, riverbank, road, road verge and road edge were measured in meters by taking the average of six measurements at every 10 m mark of the river/road plots. ArcGIS was used to check or perform the measurement of the river and road width, in case in-field measuring proved to be infeasible (see Suppl. material 1 for all plot information, including environmental data).

Species data collection
The cover of all vascular plant species was recorded in every plot. Plants were determined at species level. Taxonomy and nomenclature followed Fischer et al. (2008). Species occurrences were grouped by height and life form into the: (1) herb layer (herbaceous species of any height and woody species with or lower than 30 cm in height); (2) shrub layer (woody species between 30 and 300 cm in height); and (3) tree layer (woody species with or taller than 300 cm in height). Species were categorized as native (including archeophytes, i.e. plant species introduced before 1492) and alien (= neophytes), i.e. plant species introduced after 1492 (Essl and Rabitsch 2002). For each species, cover values (in %) were estimated visually (Suppl. material 2: Appendix 2). The ratio of woody species to herbaceous species was subsequently calculated by subtracting the cover of the herb layer from the relative cover of trees together with that of shrubs.

Data analysis
Four separate analyses were conducted to evaluate alien species distribution and community composition at the plot level. We briefly introduce the different modelling approaches and expand on the details in the subsequent sections for the respective analyses. Bayesian modeling was used to analyze 1) patterns in alien species richness in relation to elevation and proximity to rivers and roads, and 2) the effect of mowing on both alien species richness and cover. Non-metric dimensional scaling (NMDS) was applied to 3) visualize dissimilarities in plot composition among corresponding plots and those from contrasting groups (i.e. corridor type, plot type). Lastly, permutational multivariate analysis of variance (PERMANOVA) was used to 4) examine whether the presence of alien species, among other factors, is associated with a change in plot composition. All analyses were performed with the software R, version 3.6.1 (R Development Core Team 2019).

Alien species distribution
The probability of alien species occurrence in each plot type in relationship to elevation was estimated through Bayesian inference implemented in the 'brms' R package (Bürkner 2020). Given the high number of plots with zero alien plant species (72%) we used a zero-inflated Poisson model (ZIP; Lambert 2012; Bürkner 2018). The ZIP model included elevation and plot type (road/river, intermediate, and interior) as explanatory variables, and alien species richness as response variable. To account for spatial autocorrelation, we included a nested random effect with the corridor (i.e. river or road) nested in the region (i.e. Isel or Lech). Bayesian approaches are able to account for a low number of groups (< 5) and low numbers of observations per group, respectively four and sixty in this study (Harrison et al. 2018).
Weakly informative priors were set for the random and fixed effects. A student-t distribution with ten degrees of freedom, a mean of zero and a standard error of three was incorporated as the prior for the explanatory variables. For the grouping ('random') variables, a half-Cauchy distribution with a shape parameter of one was specified. To control for over-dispersion and unbalanced grouping of zeroes at higher elevations, a smoothing spline with five dimensions was fitted for elevation (Opitz et al. 2013).
To assess model convergence, each model ran four times with 5000 iterations (burn-in 2500 iterations), after which chain convergence was assessed visually (Suppl. material 2: Appendix 3). Post-run diagnostics were performed to determine whether the posterior predictions reflect the observed zero-inflation, maximum value, and the mean well. Additionally, the amount of autocorrelation was visually examined to determine the robustness of the model.
Mowing. Bayesian inference was used to assess the effect of mowing on alien species cover (%) and alien species richness. A zero-inflated beta (ZINB) model and zeroinflated Poisson (ZIP) model were applied, respectively (Ospina and Ferrari 2012;Bürkner 2018). The mowing line plot and road plot were compared for all transects subject to mowing (Fig. 2B). Again, ZINB/ZIP models with the same weakly informative priors were used. As grouping variable we specified transect.
Model convergence was evaluated by running each model four times with 5000 iterations (burn-in 2500 iterations) and a consecutive visual assessment of chain convergence (Suppl. material 2: Appendices 4, 5). The zero inflation, maximum value, and mean of the posterior predictions were compared to their observed values to determine their fit. Finally, the model was checked visually for excessive autocorrelation.
Community composition analysis. Community dissimilarity across plots was assessed via Bray-Curtis dissimilarity using non-metric multidimensional scaling (NMDS) with Wisconsin double standardization and square root transformation via the metaMDS function from R package 'vegan' (Oksanen et al. 2019). The optimal number of dimensions for the NMDS-model was selected according to its corresponding stress value using the function dimcheckMDS of the R package 'goeveg' (Goral and Schellenberg 2018). The plot type (river/road, intermediate, interior) and corridor type (river, road) were used as grouping factors and for each region a separate NMDS-model was run. Alien species were included in the analysis, thus making up part of the plot community composition.

Effect of alien species on plot composition
A permutational multivariate analysis of variance (PERMANOVA) was performed with the adonis function of the R package 'vegan' to test whether elevation, plot type, woody-to-herbaceous species ratio, and alien species cover are related to a difference in plot composition (Anderson 2008). Alien plant species were included in the dissimilarity matrices. A total of 9999 permutations were run for each test.
The corridor was treated as a stratifying term, thereby limiting the permutations to specific roads and rivers. The betadisper function of the R package 'vegan' based on 9999 permutations allowed to check for homogeneity of dispersions across the different plot types and between rivers and roads. A pairwise comparison of the levels of plot type and corridor type was made using the pairwise.adonis2 function of the R package 'pairwiseAdonis' (Arbizu 2019), with post-hoc Holm-Bonferroni adjustment of p-values (Holm 1979).

Alien and native plant species richness and composition
In total, 641 plant species were recorded in the study plots, of which 20 (= 3.1%) were alien plant species (Table 1). The total number of species (n = 500) in the Lech valley was almost identical to that of the Isel valley (n = 499). The number of species per plot in the Lech valley was significantly higher than the number of species per plot in the Isel valley (t(358) = -5.69, p < 0.001). Plant species numbers along roads (n = 511) were slightly higher than along rivers (n = 499). The road in the Isel region had the largest alien species richness (n = 15), followed by the Isel river (n = 10), the road in the Lech region (n = 5), and the Lech river (n = 4).
The most common alien species found along rivers were Solidago canadensis, Impatiens parviflora, and Impatiens glandulifera. Along roads, the most common alien species were Matricaria discoidea, Erigeron canadensis, and Solidago canadensis. In general, the upper distribution limits of alien plants were lower along rivers (833 m) than along roads (985 m). Table 1. Alien plant species recorded in the study sites. Given are the species names of all species introduced after 1942 (= neophytes), its family, the observed elevational range (lower -upper limit in m a.s.l.) in which it occurs, the total frequency of occurrence in all plots, and the occurrence for each river and road in both regions.

Species name
Family Elevational range (m)

The effect of elevation and plot type on alien species distribution
Alien species richness declined with increasing elevation and with distance from the corridor (Fig. 3, Table 2). Furthermore, alien species richness was higher in the Isel region (intercept = 1.14) compared to the Lech region (intercept = 0.84), as well as along roads (intercept Isel = 1.49; intercept Lech = 1.01) compared to rivers (intercept Isel = 0.83; intercept Lech = 0.74; Suppl. material 2: Appendix 6, Fig. 4).

Compositional dissimilarity
The NMDS for plot type in the Isel region shows a partial separation of river/road plots from the other two plot types (Fig. 6A) and no clear distinction between intermediate and interior plots. The floristic composition of plots along roads and rivers (Fig. 6B) largely overlap, although there is a subset of road plots that is clearly separated.  In the Lech region, again, river/road plots appear more distinct from the other plot types (Fig. 6C), and intermediate and interior plots have largely overlapping compositions. Our results suggest that communities along the river and the accompanying road in the Lech region are slightly more similar than those in the Isel region ( Fig. 6D and paragraph below).

Factors changing plot species composition
The PERMANOVA test results (Table 3) for the Isel region show that the occurrence of alien species has a significant effect on plot species composition. The same applies to elevation and plot type. The betadisper test implies that there are heterogenous variances among plot types (F = 6.76, p = 0.002) and the river and road (F = 4.63, p = 0.033). This means that the significant effect of plot type can be explained by a difference in species composition across plot types and possibly the difference of in-group dispersion between the plot types. The Holm-Bonferroni corrected pairwise test results (Table 3) Figure 6. NMDS based on Bray-Curtis dissimilarities for the Isel and the Lech region. Illustrated are the species compositions in the Isel region predicted through the three different plot types (A) and two corridor types (B), plus the species compositions according to plot type (C) and corridor type (D) in the Lech region. Table 4. PERMANOVA test results (n = 9999) based on Bray-Curtis dissimilarities for the Lech region. Shown are the effects of elevation, plot type, woody-herbaceous ratio, and alien species cover on plot species community. In addition, a pairwise comparison between the levels of plot type and corridor type is displayed. The unadjusted p-value is listed next to the p-value after Holm-Bonferroni correction.  Table 3. PERMANOVA test results (n = 9999) based on Bray-Curtis dissimilarities for the Isel region. Shown are the effects of elevation, plot type, woody-herbaceous ratio, and alien species cover on plot species community. In addition, a pairwise comparison between the levels of plot type and corridor type is displayed. The unadjusted p-value is listed next to the p-value after Holm-Bonferroni correction. show that directly adjacent river/road plots differ from both intermediate and interior plots in species composition, while there is no difference between the latter two. Species composition of transects along the river and the road are significantly different. The cover of alien species in the Lech region also has a significant effect on plot species composition, as does elevation and plot type ( Table 4). The exception between the two regions is that the interaction between elevation and alien species cover for the Lech region is significant as well. The variances of the different plot types display heterogeneity (F = 6.85, p = 0.002). Contrastingly, those between the river and road do not (F = 1.81, p = 0.182). The Bonferroni-corrected pairwise test results show that river/road plots differ significantly from intermediate and interior plots, while the latter two do not differ from one another (Table 4). Comparing river and road transects shows a significant difference.

Discussion
We found a strong decrease in alien species richness with elevation along rivers and roads, and in both study regions. Our findings along roads are in line with previous studies on the functionality of roads as corridors for mountain invasions Seipel et al. 2016). The distribution pattern of alien plant species along the rivers also resembles findings in other regions (Barni et al. 2012). Further, our results for rivers show that invasion patterns are similar to those encountered along mountain roads in the study area, and that habitats at high elevations are less invaded than lowand mid-elevation habitats. This decrease of alien species richness and elevation can be attributed to a decrease in temperature, as well as a decrease in human influence (Marini et al. 2009;Marini et al. 2012;Dainese et al. 2014) and is supported when comparing the lower lying Isel valley and the Lech valley in our study. However, given that we only investigated two valleys, more research is needed to establish that the observed pattern along rivers, in relation to their accompanying roads, can be generalized.

The role of rivers and roads in the spread of alien species
The higher alien plant species richness observed in plots directly located parallel to roadsides and riverbanks, as opposed to plots located further away from these networks, strongly highlights their eminent role for alien species spread in mountain valleys (Seipel et al. 2012). While it has been known that roads and rivers are important linear structures for the spread of alien species in lowland areas of central and northern Europe (Pyšek and Prach 1994;Follak et al. 2018), we show here that they also produce similar patterns for the distribution of alien plants in the Austrian Alps. We found that alien plants spread further upward along roads compared to rivers, likely because passive dispersal of species along roads is bidirectional (e.g. upwards and downwards by cars) but predominantly unidirectional along rivers (e.g. downstream drift; Wang et al. 2011;Alexander et al. 2016). Microclimatic conditions of roadsides might be more favorable with higher temperatures along roads compared to the surrounding landscape (Forman et al. 2003).
Remarkably, we found that the contribution of roads and rivers to the spread of alien plants is complementary, as different alien plants are spreading along each of these corridors. This reflects differences in disturbance and habitat characteristics along rivers and roads. Riverine habitats are subject to substantial, mainly natural, levels of disturbance caused by hydro-morphological dynamics that create ample opportunities for the spread of alien plant species (Pyšek and Prach 1994;Aronson et al. 2017).
Roadside habitats on the other hand are shaped by extensive anthropogenic habitat modifications and disturbances associated with the maintenance of road infrastructures. They are characterized by warm microclimate, frequent mowing, and often saline soil conditions resulting from the application of de-icing salt in winter that allows for the spread of salt-tolerant species (Pavol et al. 2013;Lázaro-Lobo and Ervin 2019). A process that is clearly reflected by the abundance of the facultative halophyte Puccinellia distans, which was observed in 60% of the road plots in the Lech valley and 65% of the road plots in the Isel valley.
Matricaria discoidea was commonly found along roads but absent on riverbanks, and Erigeron canadensis and Galinsoga ciliata were almost solely present next to roads. These species share a strong resistance to trampling caused by traffic, as frequently seen on roadsides (Carni and Mucina 1998). While this feature is lacking on riverbanks, these zones do have their own unique characteristics that promote the settlement of another subset of alien plant species. For example, their wet and open habitat, without mowing, gives species such as Cornus sericea, Fallopia japonica, and Impatiens glandulifera the ability to form thick stances that can remain unnoticed by conservationists (Pyšek and Prach 1994). Thus, while the specific types of disturbance regimes differ strongly between habitats associated with rivers and roads, they both facilitate the spread of alien plant species.

Effects of roadside management
Roadside management (i.e. mowing) increased alien species richness and decreased alien species cover. This might be explained by the selective pressure that mowing imposes on plants. Especially late-flowering annual species such as Impatiens glandulifera and Erigeron annuus fail to propagate when mown before seed development (Regan et al. 2006). While in our case roadside mowing was likely intended for road safety and not alien species management, in other regions it has proven to be an effective management method for limiting the spread of alien plants when the timing is correct (Milakovic et al. 2014). This way alien plant species dispersal along roads by passing vehicular traffic can be substantially reduced (von der Lippe et al. 2013). At the same time, the mowing machinery itself can act as a vector for seed dispersal across regions, a process that needs to be accounted for in management strategies (Vitalos and Karrer 2009). In addition, many short-lived alien species occurring at roadsides (e.g. Erigeron canadensis, Matricaria discoidea, Galinsoga spp.) respond adaptively to frequent mowing by developing dwarf individuals that are able to set fruits at branches growing below cutting height of conventional roadside mowing machines (Milakovic et al. 2014). Future management options in the study systems should include a rigorous monitoring of alien plant species along roadsides and their response to mowing activity. Based on these monitoring efforts, targeted measures for a mowing regime that accounts for both road safety and alien species management should be developed, including conservation scientists in the process. These measures should be consistently re-evaluated and if necessary adapted to changing conditions and potential new alien species.

Effects of alien species on resident vegetation
The decreasing effect of alien species on community composition with increasing elevation likely results from lower temperatures in high-elevation plots along roads and rivers that limits alien species performance and thus spread. For example, in Switzerland, Erigeron annuus is found to have a higher winter seedling mortality at elevations above 1000 m (Trtikova et al. 2010). Additionally, less propagule pressure and roadside disturbance due to lower traffic intensity likely reduces alien species effects (Chytrý et al. 2008;Dullinger et al. 2009;Lembrechts et al. 2016b). In the future, however, it is possible that the effect of alien plants at higher elevations will increase as more species are projected to establish in these habitats (Pyšek et al. 2012), and more alien plants will spread into surrounding mountain vegetation . The expenditure of agricultural lands in the Austrian Alps is one of the factors that can advance the introduction and distribution of alien plant species (Krausmann et al. 2003). Domesticated animals can function as a vector for alien plant seeds, by both endozoochory (ingestion) and epizoochory (carrying in fur) (Vavra et al. 2007). Similarly, tourism is an important driver of alien species invasions, as tourists can spread seeds by means of clothes, vehicles, etc. (Anderson et al. 2015). Further, the horticulture embodied by mountain villages also forms a pathway for high elevation introductions of alien plants (Kueffer et al. 2013). In our study regions, alien species appear to spread more strongly into habitats when there was some form of disturbance (i.e. logging) present, but more research and data is necessary to analyze such a relationship.

Conclusions
This study is, to our knowledge, the only one to date comparing alien species distribution across two highly different potential invasion corridors -a natural one (rivers) and a man-made one (roads) -in mountain regions. While both contribute to alien species' spread in mountain valleys, they do so for different alien species, and most species are confined to plots directly adjacent to roadsides and riverbanks. The alien plant species occurring in plots along rivers are distinct from the ones along roads and had lower upper limits than the respective sites along roads. Increased human activity and rising temperatures in alpine valleys are likely to reduce environmental constraints currently in place for alien plant species. The geographical isolation of mountain valleys and their early stage of alien plant invasion provide opportunities for conservation efforts. Nonetheless, further investigation into the spreading patterns of alien plants along alpine rivers is key to effective control.
Dr. Leni Duistermaat provided valuable suggestions for the fieldwork and comments on earlier versions of the manuscript. We particularly thank Dr. Sylvia Haider for input on the MIREN sampling protocol. We would also like to thank Dr. Johannes Wessely for fruitful discussions on the statistical methods. We appreciate the constructive comments of the Editor in Chief, Ingolf Kühn, and of two anonymous reviewers. The authors declare that there are no competing interests.
supplementary material 1