Urban sprawl facilitates invasions of exotic plants across multiple spatial scales

Exotic plant invasions are considered one of the major threats to biodiversity causing important impacts at the population, community, and ecosystem levels. Understanding the drivers of plant invasions across multiple spatial and temporal scales often requires a landscape approach. The effect of landscape composition on biological invasion has been extensively studied, whereas landscape configuration effects were seldom considered or the analyses were limited to single species. Here, we aimed to analyze how the expansion of urban and agricultural areas can affect exotic species richness (both neophytes and archaeophytes) at three spatial scales, namely regional (scale: 37.5 km2), landscape (scale: 7.1 km2) and local (scale: 100 m2). We considered the possible contribution of urban and agricultural areas both in terms of composition (i.e. habitat cover) and configuration (i.e. shape complexity of patches). First, we found that increasing urbanization coupled with high shape complexity of urban elements were major drivers of both neophyte and archaeophyte invasions across heterogeneous landscapes. In particular, shape complexity seemed to be a key driver of plant invasions at large spatial scale, whereas the type of recipient habitat and urban cover determined the exotic success at the patch level. Second, archaeophytes were also affected by agriculture land use, i.e. agricultural patches shape complexity increased their spread at both regional and landscape scales. High shape complexity of highly disturbed habitats is expected to increase the exchange surface that exotic plant use to spread their propagules across the landscape mosaics. Our findings suggest that urban planning aimed at curbing urban fragmentation by both reducing shape complexity and diffuse urban sprawl might greatly improve the resistance of landscapes to biological invasions.


Introduction
Exotic plant invasions are considered one of the major threats to biodiversity (Vilà et al. 2010;Simberloff et al. 2013;Vilà and Hulme 2017;Pyšek et al. 2020) causing important impacts at the population, community and ecosystem levels (Vilà et al. 2011;Boscutti et al. 2020;Pellegrini et al. 2021a). Understanding the drivers of plant invasions across multiple spatial and temporal scales often requires a landscape approach (Vilà and Ibáñez 2011). First, the role of landscape composition, i.e. proportion of landscape covered by different habitats, has received increasing attention, as human-impacted landscapes are expected to harbor a larger pool of exotic species and increase the likelihood of propagule dissemination to natural areas (Jauni and Hyvönen 2010;Basnou et al. 2015;Andersen et al. 2015;Pellegrini et al. 2021b). Anthropogenic land-uses such as agricultural and urban habitats represent primary sources of alien plants by providing empty niches released from competition with natives (Lonsdale 1999;Davis et al. 2000;Vilà et al. 2007;Polce et al. 2011;Boscutti et al. 2018). In particular, urban land use encompasses a large array of disturbed habitats (e.g. roadsides, gardens) that are usually highly invaded (Chytrý et al. 2008;Lambdon et al. 2008;Aikio et al. 2012) and that can act as source of exotic plants that can spread into surrounding habitats (Lippe and Kowarik 2008; Botham et al. 2009;Säumel and Kowarik 2010;Kühn et al. 2017;Guo et al. 2019). While most of these studies have focused on landscape composition, the effects of landscape configuration of urban and agriculture areas have been seldom investigated empirically (but see Chen et al. 2017;O'Reilly-Nugent et al. 2016).
Landscape configuration effects have been mostly investigated to study plant invasions within natural fragmented habitats embedded in intensive matrices (Hastings et al. 2005). This body of research showed that small fragments usually have a higher density of exotic species than larger fragments (Guirado et al. 2006;Ohlemüller et al. 2006), and that an increase in exotic plant density, richness and cover is often observed at the edge between natural habitats and more intensive land-uses such as urban or agricultural areas (Cilliers et al. 2008;Boscutti et al. 2018). These results suggest that increased patch shape complexity is expected to promote the spread of exotic plants across human-impacted landscape mosaics (Yates et al. 2004;Song et al. 2005;Malavasi et al. 2014). In particular, landscapes with large contagion, high shape complexity and high patch density are expected to promote exotic invasions by increasing the dispersal opportunities of exotics between urban areas and the surrounding habitats and the availability of disturbed edges (Aikio et al. 2012;Bar-Massada et al. 2014). Upscaling the effect of spatial configuration from the habitat level to the whole landscape is an urgent task to predict the effect of land-use changes on plant invasions across multiple spatial scales.
When evaluating the drivers of exotic plant invasions, it is important to consider how these effects vary across spatial scales Pauchard and Shea 2006;Foxcroft et al. 2009). Most of previous studies adopting a multi-scale approach has focused on the relationship between native and alien species richness (Campos et al. 2013), whereas only a few have tested for scale-dependence in the effect of abiotic drives of plant invasions (but see Cabra-Rivas et al. 2016). At broad spatial scales, climate, soil and habitat types are expected to be more influential in determining species distribution (Araújo et al. 2005;Guisan et al. 2006;Bailey et al. 2017) while at finer scales, disturbance and biotic interactions should become more important (Bartuszevige et al. 2006;Milbau et al. 2009).
In the European floras, the status of exotic plants is commonly further classified into two categories: archaeophytes (i.e. plants introduced before 1500) and neophytes (i.e. plants introduced after 1500). Archaeophytes and neophytes often exhibit consistent distribution patterns, suggesting that similar mechanisms can be responsible for their spread (Chytrý et al. 2008). Nonetheless, archaeophytes are often ruderal or segetal species mainly linked to agricultural land uses (Chytrý et al. 2008;Sheppard and Brendel 2021;Zając et al. 2009), while neophytes are expected to be more related to urban areas. In this light, studying the response of these two groups to landscape composition and configuration might provide further insight into the invasion mechanisms.
In this study, we aim to analyze how the expansion of urban and agricultural areas can affect exotic species richness (i.e. neophytes and archaeophytes) at three spatial scales, namely regional (scale: 37.5 km 2 ), landscape (scale: 7.1 km 2 ) and local (scale: 100 m 2 ). We considered the possible contribution of urban and agriculture areas both in terms of composition (i.e. habitat cover) and configuration (i.e. patch shape complexity) (Fig. 1). We hypothesize that increasing cover of urban areas and increasing shape complexity will increase invasion of exotic species across the landscape by increasing colonization opportunities linked to disturbance and by facilitating dispersal of exotic plants across the interface where human settlements adjoin or intermix with other ecosystems. We also expect that these effects should be stronger for neophytes than for the more agriculturedependent archaeophytes.
The diversity patterns of exotics (i.e. neophyte = post-1500 immigrants; archaeophytes = pre-1500 immigrants) in the study region (ca. 3,000 km 2 ) was studied at three spatial scales: (i) regional (ii) landscape (iii) local habitat. First, the whole area was split into a regular grid of 83 cells, each of 3′ × 5′ (arcmin), corresponding to ca. 35.7 km 2 (i.e. regional scale) (Fig. 2). Within each grid cell, exotic species distribution was extracted from an ongoing floristic inventory of the Friuli Venezia Giulia, which contained over 290,000 records collected in the period 1980-2017. Second, 15 circular landscapes (ca. Fig. 1 Representation of tested gradients of landscape composition (e.g. cover of urban areas) and configuration (e.g. shape complexity of urban elements) at regional scale, land use considered were: urban (dark grey), agricultural (green) and seminatural (white) urban patch shape complexity urban land use cover For each landscape, a detailed land use map (scale 1:2000) was realized distinguishing agricultural (i.e. annual and perennial crops), urban and semi-natural habitat (i.e. forests hedgerows and meadows). In each landscape, 20 sampling points (10 × 10 m) were selected for a total of 300 (20 × 15) plots for flora survey at local and landscape scale. Within each landscape, the number of points per habitat type reflected the relative abundance of the habitat. The surveyed habitat included both agricultural (i.e. annual crops, N = 122; perennial crops N = 55) and semi-natural habitats (forests, N = 71; hedgerows, N = 13; meadows, N = 28). The complete list of sampled habitats and number of plots per landscape can be found in supplementary materials (Online Resource 1). Plots were randomly distributed within each landscape excluding landscape edges. At this scale, we did not sample urban habitats (e.g. parks, gardens, yards), as we wanted to test the habitat and landscape invasibility of non-urban habitats.
In summer 2017, all the plots were surveyed in the field recording all occurring vascular species. During the surveys, 11 plot were discarded due to accessibility issues. Number of exotic species were calculated at the plot, landscape and regional scale by pooling species list at different spatial scales. Nomenclature followed the Italian check-list (Bartolucci et al. 2018;Galasso et al. 2018). Exotic status (i.e. neophyte and archaeophyte) was classified according to the regional checklists (Poldini 2002;Buccheri et al. 2019). Landscape composition and configuration were assessed at all investigated scales in terms of land-use cover and patch shape complexity using (i) the official regional land use map, for the regional scale (scale 1:25,000, IRDAT FVG, realized in 2017), and, (ii) the detailed land use maps realized for each of the 15 surveyed landscapes, for the landscape and plot scales. Landscape composition was quantified as cover of urban land use in each cell grid (35.7 km 2 ) and landscape (7.1 km 2 ). At habitat level, urban cover was calculated in a buffer with a radius of 250 m around each plot, which has been already proved to be a sensitive scale when analyzing exotic plant spread in similar landscapes (Kumar et al. 2006;Boscutti et al. 2018). Urban areas ranged between 1.9-57.3%, 0.8-24.5% and 0.0-51.4% at the regional, landscape and plot scale, respectively. Arable land ranged between 26.8-92.4%, 46.4-97.9% and 0.3-100.0%; at regional, landscape and plot scale, respectively. The cover of urban areas and agricultural land was not highly correlated at the three scales of analysis (regional, r = − 0.45; landscape, r = − 0.21; plot scale, r = − 0.17).
Landscape configuration was assessed in terms of patch shape complexity by using the Area Weighted Mean Patch Fractal Dimension (AWMPFD) index (McGarigal and Marks 1995) for both urban and agriculture patches. The patch fractal dimension represents a comparison of the perimeter and area of each shape, similarly to the shape index, but testing the exponent of that relationship (Milne 1988). Fractal dimension is a measure of the complexity of shapes that gives more emphasis on the complexity of form while the common shape index stresses the elongation of patches. A fractal dimension greater than 1 for a 2-dimensional landscape mosaic indicates a departure from a Euclidean geometry (that is, an increase in patch shape complexity). AWMPFD approaches 1 for shapes with very simple perimeters, such as circles or squares, and approaches 2 for shapes with highly convoluted, plane-filling perimeters, this index is weighted by patch area so that larger patches weigh more than smaller patches. In the study, AWMPFD will be referred as shape complexity. All the analyses were performed using the Patch Analyst exten-sion© for ArcGIS 10.0 (ESRI). Preliminarily, we also assessed other landscape configuration metrics, namely the total edge (TE), as the sum of the lengths (m) of all edge segments involving the corresponding patch type, and the Area-Weighted Mean Shape Index (AWMSI), an area-weighted index of shape complexity (McGarigal and Marks 1995). All indices showed a high correlation with AWMPFD (r > 0.5) at all scales and were therefore excluded from the analyses (Online Resource 2).

Data analysis
To explore the effects of landscape composition (proportion of urban and agriculture areas) and shape complexity at different spatial scales on exotic plant richness (i.e. neophytes and archaeophytes), Linear Models (LM) and Generalized Linear Models were run in R statistical environment (R Core Team 2020). At the regional and landscape scale, we fitted a LM testing the effect of proportion of urban and agriculture areas and shape complexity on the species number of both neophytes and archaeophytes.
At the habitat scale, we fitted a generalized mixedeffects model with a Poisson distribution testing proportion of urban and agricultural areas and habitat type and their two-way interactions on the number of exotics (i.e. neophytes and archaeophytes) in the 10 × 10 m plot. Then, we simplified the models by removing one-by-one the non-significant interaction terms (p > 0.05) using a manual backward selection procedure. We included landscape ID as random factor to account for the spatial dependence in the design. The interaction between shape complexity and habitat type was not tested at the plot scale, where the low number of patches did not allow a solid analysis of landscape configuration. Full models are shown in Online Resource 3.The model assumptions were checked using diagnostic plots of model residuals using the car and DhARMA packages in R. Due to the better residual diagnostics outcome, LMs were applied at regional and landscape scale even if species richness measures were counts.
Our findings highlighted scale-dependent effects of landscape composition and patch shape complexity on exotic species richness (Table 1). At the regional Table 1 Results of the models testing the effects of urban and agriculture cover (%) and patch shape complexity (Area Weighted Mean Patch Fractal Dimension-AWMPFD) on neophyte and archaeophyte species richness at regional (a), landscape (b), and local (c) scale At local scale, the interaction between patch shape complexity and habitat type was not considered scale (n = 83 grid cells), we found both urban cover and urban shape complexity to have a positive effect on the number of neophytes ( Fig. 3a-b). The model explained 37.3% of the total variation. Archaeophyte species richness increased in areas with high urban cover and agriculture shape complexity (Fig. 3c-d).
The model accounted for the 44.1% of the total variation. At the landscape scale (n = 15 landscapes), only shape complexity of urban and agriculture patches affected neophyte and archaeophyte species richness, showing a linear positive relationship ( Fig. 4; neophytes: R 2 = 37.8%, archaeophytes: R 2 = 19.7%). At the local scale (n = 289 plots), semi-natural habitats were less invaded than agricultural ones for both neophyte and archaeophyte species (Fig. 5a-c). Neophyte species number increased in plots surrounded by high percentage of urban areas (Fig. 5b), irrespective of habitat type. Archaeophyte species number increased in plots surrounded by high percentage of both agriculture and urban areas (Fig. 5d-e). In particular, (c) (d) Fig. 3 Relationships between neophyte species number, urban cover (a) and urban shape complexity (b), and archaeophyte species number, urban cover (c) and agriculture shape com-plexity (c), at regional scale. Shape complexity is expressed as Area Weighted Mean Patch Fractal Dimension (AWMPFD). Shaded area is the 95% confidence interval urban cover showed a stronger effects on archaeophytes invading semi-natural areas, when compared to agricultural ones (Fig. 5e).

Discussion
The contribution of landscape composition on biological invasion has been extensively studied, whereas landscape configuration has been seldom considered (Vilà and Ibáñez 2011) or limited to single species (Song et al. 2005;Montti et al. 2017). We found that increasing urbanization coupled with high shape complexity of urban elements were major drivers of neophyte invasions. In particular, shape complexity seemed to be the most important abiotic driver of plant invasions at large spatial scale whereas the type of recipient habitat and urban cover determined the exotic success at scale of habitat patches. Moreover, we found that increasing shape complexity further increased the spread of archaeophytes at both the regional and landscape scales. At the regional scale, we found that both increasing urban cover and patch shape complexity contributed to increase the number of exotic species. Urban areas are often over-represented among the sites of first introduction due to high propagule pressure (Pyšek 1998;Chytrý et al. 2008;Aikio et al. 2012). Urban land use creates a mosaic of different humanmade habitats (e.g. roadsides, gardens) that are often highly invaded and that are expected to promote the spread of exotic species also to the surrounding habitats (Lippe and Kowarik 2008; Lambdon et al. 2008;Botham et al. 2009;Säumel and Kowarik 2010). A key role is played by roads that are often considered effective dispersal corridors (Arévalo et al. 2005;Follak et al. 2018;McDougall et al. 2018; Lázaro-Lobo and Ervin 2019). A high patch shape complexity in the landscape is generally associated with a low overall disturbance intensity but more diffused across the landscape matrix (Krummel et al. 1987;O'Neill et al. 1988), which in rural landscapes has been proved to favor overall plant species richness (Moser et al. 2002).
This pattern was even stronger at the landscape scale where we only found that shape complexity of urban and agriculture areas affected the spread of neophyte and archaeophyte plants species, respectively. This suggests a possible relationship between the success of spread and dispersal of exotic species and the configuration of the invasion foci (i.e. urban settlements for neophytes and arable land for archaeophytes). So far, a few studies have related plant invasion with landscape shape complexity and patch edge density (Bartuszevige et al. 2006;Kumar et al. 2006). Indeed, the shape complexity of urban areas might have amplified the exchange surface between habitats that neophytes use to spread their propagule across the landscape mosaics. A low-density of urban land-cover is often characterized by fragmented and disturbed natural areas within a loose matrix of Relationships between exotic species number (i.e. neophytes or archaeophytes) and shape complexity of urban (a) and agriculture (b) patches at landscape scale. Shape complexity is expressed as Area Weighted Mean Patch Fractal Dimension (AWMPFD). Shaded area is the 95% confidence interval gardening and housing areas (Gavier-Pizarro et al. 2010). This might produce a greater probability of exotic ornamental species to spread from gardens and urban areas into adjacent natural areas (Gavier-Pizarro et al. 2010;González-Moreno et al. 2013 (Marco et al. 2010;Guo et al. 2019).
In addition, urban shape complexity is also linked to a higher development of road network, which often determines the invasion success of exotic species within semi-natural habitats (Hansen and Clevenger 2005;Dainese et al. 2017;Riitters et al. 2018). In contrast to neophytes, the richness of archaeophytes at the landscape scale increased with high patch complexity of agricultural areas. Other studies showed that archaeophytes are more linked to agricultural than to urban land use (Chytrý et al. 2008). Our findings suggest that, at large spatial scale, increasing the exchange surface between arable land and seminatural habitats provide a greater chance for archaeophytes to spread across the territory. Previous studies have showed archaeophytes to be particularly common in landscapes with small fields rather than in landscapes dominated by large fields (Baessler and Klotz 2006). At the local scale, we found that agricultural habitats (i.e. annual and perennial crops) were more invaded than semi-natural habitats by both neophyte and archaeophyte species. In general, highly humanimpacted habitats are widely recognized for their high level of plant invasion (e.g. Aikio et al. 2012;Chytrý et al. 2008;Lambdon et al. 2008;Maskell, Bullock, Smart, Thompson, and Hulme 2006). At this scale also patch features might influence invasion (Milbau et al. 2009) and local variations in soil quality and vegetation structure might become pivotal in the invasion success (Ohlemüller et al. 2006;Vitti et al. 2020). Our findings suggest that urban areas may be the main source of exotic plant invasion toward other habitats. Invasion by exotic species often begins in human-disturbed environments, where their propagules are unintentionally introduced or where cultivated species frequently escape (Pyšek 1998;Hulme 2005;Čeplová et al. 2017;Geppert et al. 2021). Other studies showed that semi-natural habitats are more prone to be invaded when included in a urban-rather than agricultural-dominated matrix (Cilliers et al. 2008) and that at the local scale the proximity to roads might be pivotal in the invasion process (Zeeman et al. 2018). A consistent pattern was found for the two exotic groups when considering urban and habitat effects. In addition, archaeophytes were also favored by high agriculture cover in the surrounding of the plot (Preston et al. 2004). Our findings would suggest that arable land might be concurrently a source and a sink of archaeophyte plants (Hulme 2005), as it also emerged in studies on the effects of anthropic disturbance in agricultural landscapes (Boscutti et al. 2018).

Conclusions
Clarifying the mechanisms through which exotic species spread across the landscape might inform an effective management/planning aimed at curbing biological invasions. Our study elucidated how plant invasion was affected by alterations of landscape composition and configuration. We showed that urbanization was a major driver of plant invasions showing that its shape complexity of urban areas was fundamental in the process. In particular, landscape shape complexity seemed to be a facilitator of plant dispersal and spread at the landscape scale whereas the habitat characteristics and landscape composition determine the exotic success at the local scale. Such evidence suggests that urban planning should consider reducing urban fragmentation to reduce the exchange surface between anthropogenic and natural habitats. These results allow to gain a better understanding of the way landscape composition and configuration contribute in modifying plant invasions, providing far reaching consequences in invasion ecology and landscape planning.

Fig. 5
Relationships between neophyte species number, habitat type (a) and urban cover (b), and archaeophyte species number, habitat type (c), agriculture cover (d), and interaction between urban cover and habitat type (e), at local scale ◂