Short Communication
Short Communication
Aliens in Transylvania: risk maps of invasive alien plant species in Central Romania
expand article infoHeike Zimmermann, Jacqueline Loos, Henrik von Wehrden§, Joern Fischer|
† Leuphana University, Lüneburg, Germany
‡ Institute of Ecology, Faculty of Sustainability, Leuphana University, Lueneburg, Germany
§ Institute of Ecology, Faculty of Sustainability, 21335 Lueneburg, Germany
| Leuphana University, Lueneburg, Germany
Open Access


Using the MAXENT algorithm, we developed risk maps for eight invasive plant species in southern Transylvania, Romania, a region undergoing drastic land-use changes. Our findings show that invasion risk increased with landscape heterogeneity. Roads and agricultural areas were most prone to invasion, whereas forests were least at risk.


Erigeron annuus, MAXENT, Robinia pseudoacacia, Romania, Solidago canadensis, Xanthium strumarium


Species distribution models are a useful tool in biological invasion risk management (Jiménez-Valverde et al. 2011). They allow a rapid assessment of the distribution of alien species, which may serve to identify areas of high invasibility and to understand the mechanisms behind the establishment and spread of alien species (Zimmermann et al. 2011).

We focused our study on southern Transylvania, in Central Romania, where temporary or permanent abandonment of agricultural land is common. Knowledge on the introduction history of invasive plant species and their current distribution in this region is largely missing, although several common alien plant species are among the world’s 100 worst invaders (DAISIE 2013; GISD 2013). To address this knowledge gap, we aimed to identify predictors of the distribution of eight highly invasive species.


Study area

Our study area comprised an area of 7,440 km2 (Fig. 1). It is a heterogeneous, hilly, rural landscape, at elevations from 230 to 1100 m, and characterized by a mosaic of different land cover types supporting high farmland biodiversity (Loos et al. 2014). Small-scale semi-subsistence farming with little use of machinery or agrochemicals, as well as extensively used hay meadows and grasslands are typical for the region. The small-scale mosaic pattern of different land cover types prevails throughout the study area, although the North-West has a relatively higher percentage of arable land, the South has a relatively higher percentage of pastures, and the central part a higher percentage of forests (Hanspach et al. 2014). Like many other parts of Eastern Europe, Romania has experienced drastic socio-economic changes, with the collapse of the communist regime in 1989, and its accession to the European Union in 2007 (Kuemmerle et al. 2009). Among other changes, the communist regime and its collapse triggered mass emigration from some areas, resulting in widespread farmland abandonment (Fig. 2).

Figure 1.

Location of our study area in Romania. Inside the enlarged map of our study area the cities Mediaș and Sighișoara are outlined and black points represent the presence points of all eight study species.

Figure 2.

Population census data from 1990, 2000 and 2005 for a subsample of 22 communes across the study area. One commune typically includes four villages. In Romania the communist regime and its collapse in 1989 triggered mass emigration from these areas.


In summer 2013, we mapped presences of eight prominent alien plant species across the study area using a handheld global positioning system (Table 1, see also Suppl. material 1). Our sampling covered a large variety of landscape elements within 50 km of the town of Sighișoara (Fig. 1). In each local valley, we undertook at least two extended survey walks that lasted between 30 and 180 minutes. We sampled along roads, but also tracked species off-road, by walking towards the top of the hills bordering a given valley (see Suppl. material 2 for details). Furthermore, we took sample points whenever we observed an invasive species during our other field activities. For Robinia pseudoacacia we did not record planted trees, but only naturally dispersed individuals. We did not measure the distance to the next adult tree, but the minimum distance was approx. 100 m and for most recordings adult trees were not in sight. We also included vegetation sampling data from 2012, which was based on a randomized stratified design. The minimum distance between sampling points was 30 m. In combination, our sampling approach covered a wide range of environments across the study area.

Overview of study species and number of sampling points (N Am = North America).

Species Family Common name Life strategy Origin Reproduction/dispersal Presence points
Amaranthus retroflexus L. Amaranthaceae Redroot amaranth annual herb
~1 m
N Am monoecious, wind pollinated, dispersed by wind, water and animals 45
Asclepias syriaca L. Apocynaceae Common milkweed perennial herb
~1–2 m
N Am insect pollinated, seeds wind dispersed, and vegetative reproduction (rhizomes) 65
Conyza canadensis (L.) Cronquist Asteraceae Canadian horseweed annual herb
~1 m
N Am insect pollinated, self- and cross-fertilization, seeds wind dispersed 35
Erigeron annuus
(L.) Pers.
Asteraceae Annual fleabane annual herb
~1 m
N Am insect pollinated, self- and cross- fertilization, winged achenes dispersed by wind and animals 475
Fallopia japonica (Houtt.) Ronse Decr. Polygonaceae Japanese knotweed perennial herb
~3 m
Asia insect pollinated, dioecious, winged achenes dispersed by wind, water, animals, and reproduces vigorously by rhizomes 69
Robinia pseudoacacia L. Fabaceae Black locust deciduous tree
~30 m
N Am insect pollinated, seeds wind dispersed, reproduces vigorously by root suckering and stump sprouting 264
Solidago canadensis L. Asteraceae Canadian goldenrod perennial herb
~2.5 m
N Am insect pollinated out-crossing, wind dispersed achene with pappus, and vegetative reproduction (rhizomes) 298
Xanthium strumarium L. Asteraceae Common cocklebur annual herb
~1 m
N Am wind-pollinated, monoecious, self- and cross- fertilization, apomixis, seeds dispersed by animals and water 236


We derived invasion risk maps for each species individually. To this end, we applied the Maximum Entropy algorithm (MAXENT), which is based on presence only data to map the likely current distribution for each species in our study area (Phillips et al. 2006; Elith et al. 2011; Merow et al. 2013). Our approach to model the likely current distribution based on presence only data is in line with our assumption that the expansion of our species has not yet reached equilibrium, making the use of observed absences misleading (Jiménez-Valverde et al. 2011). Clamping was activated as well as random seed, and we made an internal validation with 50 replicates and 20 % random test percentage. Predictors encompassed information on the topography, land use and potential distribution corridors (Table 2). Slope and terrain ruggedness (the variation in altitude) did not improve any of the models. Our final models therefore included four predictors (road and village distance, heterogeneity, and land cover classes). Each of these predictors was re-sampled to a cell size of 30 x 30 m within ArcGIS, which thus equals the resolution of our risk maps. To calculate the potential areas suitable for each species a threshold rule was applied to convert continuous suitability surfaces into binary outputs. We selected the threshold “maximum training sensitivity plus specificity logistic threshold”, which focuses on the correct classification of presences and background points (Jiménez-Valverde and Lobo 2007).

For each grid cell, we then calculated the mean probability of occurrence over all eight species. This resulted in a map of general invasion risk for the study area, referred to “the invasibility map” hereafter.

Predictors for the MAXENT model. All predictors have a 30 x 30 m resolution. (h.s. = habitat suitability).

Predictor Description Relative importance in the MAXENT model
Road distance minimum distance to the closest road 18 to < 40 % (A. retroflexus, A. syriaca, C. canadensis)
45–48 % (E. annuus, R. pseudoacacia, S. canadensis)
>50%–56% (F. japonica, X. strumarium)
for all species h.s. high with decreasing distance
Village distance minimum distance to the closest village 15 % (S. canadensis) high h.s. at 3-4 km
for remaining species values <10%
Heterogeneity (CNES 2007, Distribution Spot Image SA) variation in the panchromatic channel of SPOT 5 satellite imagery >20 %–35% (A. syriaca, C. canadensis, E. annuus, F. japonica, S. canadensis, X. strumarium)
>40%–46% (A. retroflexus, R. pseudoacacia)
for all species h.s. high with increasing heterogeneity
Corine land cover classes
(Corine 2006 Land Cover Map, EEA (2006) Corine land cover 2000 -- A seamless vector database (European Environment Agency, Copenhagen)
(1) broad leaved forest
(2) coniferous and mixed forest
(3) water
(4) inland marshes
(5) natural areas (sparsely vegetated, bare rocks, natural grasslands)
(6) transitional woodland-shrub habitats
(7) artificial surfaces
(8) agriculture
(9) pasture
(10) Land principally occupied by agriculture with significant areas of natural vegetation
15 % (X. strumarium high h.s. for class 8)
22% (A. retroflexus high h.s. for class 4)
23 % (A. syriaca high h.s. for classes 8, 9)
35% (C. canadensis high h.s. for classes 8, 10)
for remaining species values <15 %


Single distribution models of the eight study species all had high discrimination performances with AUC values ranging from 0.8 to 0.9 (Elith 2002). Standard deviation of the output grids showed no signs of autocorrelation (see also Suppl. material 2). Road distance and heterogeneity were the most important predictors (Table 2). Corine land cover classes substantially improved the MAXENT model of three species, with agricultural areas, pastures and land principally occupied by agriculture with significant areas of natural vegetation, having the highest risk of invasion. Forest was at least risk of being invaded. With 31 % potential habitat, Erigeron annuus has the largest potential distribution and with 8 %, Asclepias syriaca had the lowest (Fig. 3 a, b). Risk maps of all species show a high invasibility of distribution corridors such as roads and rivers (see Suppl. material 3). The overall invasibility map therefore shows a network structure with the highest invasibility close to distribution corridors stretching across the whole study area (Fig. 3c). However, the invasibility was highest in heterogeneous areas, which is where it also reached furthest away from roads.

Figure 3.

The likely current distribution for a Erigeron annuus and b Asclepias syriaca in our study area. c Shows the invasibility, calculated as the mean probability of occurrence over all eight study species.


Our risk maps show that the eight invaders considered have great potential to further expand their distributions. All except for one study species are wind dispersed, which is an effective long distance dispersal method (Cain et al. 2000). Roads serve as invasion corridors, enabling the species establishment inside disturbed road margins (Birdsall et al. 2011). Environments at greatest risk of invasion away from roads appeared to be heterogeneous agricultural areas. Landscape heterogeneity is beneficial for native species in farming landscapes, because it offers a greater variety of habitats (Fischer and Lindenmayer 2007). However, invasive species also appear to benefit from this variety. A positive correlation of native and invasive plant species richness in relation to landscape heterogeneity was also shown for the Rocky Mountain National Park in the USA (120-960 m radius scale; Kumar et al. 2006), for rural and urban areas in Germany (32 km2 scale; Deutschewitz et al. 2003), as well as for rural areas and woodlands in Catalonia (100 km2 scale; Bartomeus et al. 2012). Furthermore, heterogeneous landscapes offer more edge environments which are subject to higher levels of propagule pressure and disturbance, and therefore a higher abundance of invasive plant species (Vilà and Ibáñez 2011). Like many invasive species, our study species are adapted to disturbance through their short life cycles, rapid germination, persistent seed banks, or vegetative reproduction (Rejmánek and Richardson 1996; Kolar and Lodge 2001).

We observed that areas with a high heterogeneity often coincides with areas that experienced the most widespread emigration following the collapse of communism. Socio-economic effects at regional or local scales are rarely considered in invasion science (Vilà and Pujadas 2001, Guo et al. 2012). This is a potentially important oversight. For example, farmland abandonment and land-use change can facilitate high abundances of alien plant species (Hobbs 2000; Cramer et al. 2008). Decades of cultivation may deprive native seed banks and altered soil conditions may no longer be suitable for native species (Cramer et al. 2008). In comparison to native species, species new to this region may possess traits that make them better colonizers of abandoned fields. They benefit from increased human mediated propagule pressure (propagule bias), are often stronger competitors (enemy release, novel weapons), and once established may further transform their environment, making it even less suitable for native species – thereby creating self-perpetuating alien monocultures (Catford et al. 2012). Our study area continues to undergo socio-economic changes. With our main project (“Fostering sustainable development in Eastern Europe” we aim to provide leverage points for enabling sustainable land use practices. Based on the results of this study we recommend to preserve the remaining intact forests, which seem to be least in risk of invasions and to further investigate the role of land abandonment and land use change in our study area.


We thank all the people in the field who helped us to track down aliens in Transylvania. We are also thankful to P. Brandt and P. Fust for their technical advice. This study was funded by a Leuphana small research grant (HZ) and through a Sofja Kovalevskaja Award by the Alexander von Humboldt Foundation (JL, JF).


  • Bartomeus I, Sol D, Pino J, Vicente P, Font X (2012) Deconstructing the native–exotic richness relationship in plants. Global Ecology and Biogeography 21: 524–533. doi: 10.1111/j.1466-8238.2011.00708.x
  • Birdsall JL, McCaughey W, Runyon JB (2011) Roads impact the distribution of noxious weeds more than restoration treatments in a lodgepole pine forest in Montana, USA. Restoration Ecology 20: 517–523. doi: 10.1111/j.1526-100X.2011.00781.x
  • Cain ML, Milligan BG, Strand AE (2000) Long-distance seed dispersal in plant populations. American Journal of Botany 87: 1217–1227. doi: 10.2307/2656714
  • Catford JA, Daehler CC, Murphy HT, Sheppard AW, Hardesty BD, Westcott DA, Rejmánek M, Bellingham PJ, Pergl J, Horvitz CC, Hulme PE (2012) The intermediate disturbance hypothesis and plant invasions: Implications for species richness and management. Perspectives in Plant Ecology, Evolution and Systematics 14: 231–241. doi: 10.1016/j.ppees.2011.12.002
  • Cramer VA, Hobbs RJ, Standish RJ (2008) What´s new about old fields? Land abandonment and ecosystem assembly. Trends in Ecology & Evolution 23: 104–112. doi: 10.1016/j.tree.2007.10.005
  • Deutschewitz K, Lausch A, Kühn I, Klotz S (2003) Native and alien plant species richness in relation to spatial heterogeneity on a regional scale in Germany. Global Ecology and Biogeography 12: 299–311. doi: 10.1046/j.1466-822X.2003.00025.x
  • Elith J (2002) Quantitative methods for modeling species habitat: comparative performance and an application to Australian plants. In: Ferson S, Burgman M (Eds) Quantitative methods for conservation biology. Springer, New York, 39–58.
  • Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ (2011) A statistical explanation of MaxEnt for ecologists. Diversity and Distributions 17: 43–57. doi: 10.1111/j.1472-4642.2010.00725.x
  • Fischer J, Lindenmayer DB (2007) Landscape modification and habitat fragmentation: a synthesis. Global Ecology and Biogeography 16: 265–280. doi: 10.1111/j.1466-8238.2007.00287.x
  • Guo Q, Rejmanek M, Wen J (2012) Geographical, socioeconomic, and ecological determinants of exotic plant naturalization in the United States: insights and updates from improved data. NeoBiota 12: 41–55. doi: 10.3897/neobiota.12.2419
  • Hanspach J, Hartel T, Milcu A, Mikulcak F, Dorresteijn I, Loos J, von Wehrden H, Kuemmerle T, Abson DJ, Kovács-Hostyánszki A, Báldi A, Fischer J (2014) A holistic approach to studying social-ecological systems and its application to Southern Transylvania. Ecology and Society 19(4): 32. doi: 10.5751/ES-06915-190432
  • Hobbs RJ (2000) Land-use changes and invasions. In: Mooney HA, Hobbs RJ (Eds) Invasive Species in a Changing World. Island Press, Washington DC, 55–64.
  • Jiménez-Valverde A, Lobo JM (2007) Threshold criteria for conversion of probability of species presence to either–or presence–absence. Acta Oecologia 31: 361–369. doi: 10.1016/j.actao.2007.02.001
  • Jimenez-Valverde A, Peterson AT, Soberon J, Overton JM, Aragon P, Lobo JM (2011) Use of niche models in invasive species risk assessments. Biological Invasions 13: 2785–2797. doi: 10.1007/s10530-011-9963-4.
  • Kolar CS, Lodge DM (2001) Progress in invasion biology: predicting invaders. Trends in Ecology & Evolution 16: 199–204. doi: 10.1016/S0169-5347(01)02101-2
  • Kuemmerle T, Müller D, Griffiths P, Rusu M (2009) Land use change in Southern Romania after the collapse of socialism. Regional Environmental Change 9: 1–12. doi: 10.1007/s10113-008-0050-z
  • Kumar S, Stohlgren TJ, Chong GW (2006) Spatial heterogeneity influences native and nonnative plant species richness. Ecology 87: 3186–3199. doi: 10.1890/0012-9658(2006)87
  • Loos J, Dorresteijn I, Hanspach J, Fust P, Rakosy L, Fischer J (2014) Low-intensity agricultural landscapes in Transylvania support high butterfly diversity: implications for conservation. PloS One 9: e103256. doi: 10.1371/journal.pone.0103256
  • Merow C, Smith MJ, Silander JA (2013) A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography 36: 1058–1069. doi: 10.1111/j.1600-0587.2013.07872.x
  • Phillips SJ, Anderson RP, Schapire RE (2006) Maximum entropy modeling of species geographic distributions. Ecological Modelling 190: 231–259. doi: 10.1016/j.ecolmodel.2005.03.026
  • Rejmánek M, Richardson DM (1996) What attributes make some plant species more invasive? Ecology 77: 1655–1661. doi: 10.2307/2265768
  • Vilà M, Ibáñez I (2011) Plant invasions in the landscape. Landscape Ecology 26: 461–472. doi: 10.1007/s10980-011-9585-3
  • Vilà M, Pujadas J (2001) Land-use and socio-economic correlates of plant invasions in European and North African countries. Biological Conservation 100: 397–401. doi: 10.1016/S0006-3207(01)00047-7
  • Zimmermann H, von Wehrden H, Damascos MA, Bran D, Welk E, Renison D, Hensen I (2011) Habitat invasion risk assessment based on Landsat 5 data, exemplified by the shrub Rosa rubiginosa in southern Argentina. Austral Ecology 36: 870–880. doi: 10.1111/j.1442-9993.2010.02230.x