INTRODUCTION
The extent of urban areas is increasing worldwide and most of the human population now lives in towns and cities (Johnson and Munshi-South 2017). Urbanization has a strong impact on natural habitats and resources (Murgui and Hedblom 2017), and such impact can have effects on not only species distribution and abundance, but also on the composition of biological communities (Lancaster and Rees 1979, Dowd 1992). For instance, the introduction of tree brown snake (Boiga irregularis) was responsible for forest bird extinction on Guam Island (Savidge 1987); and the release of pig, black rat and other species on islands have had a negative effect on endemic species (Blackburn et al. 2009, Donald et al. 2010).
Urban populations are expected to rapidly grow in the future (United Nations Organization c2012), and this trend suggests that their activities will substantially transform natural habitats to meet a growing population's demands with a higher levels of resources consumption (Hutton and Leader-Williams 2003). Therefore, the emergence of new small urban areas, near big towns and cities, would be an ideal environment to facilitate the expansion of the Eurasian collared dove (Streptopelia decaocto). This species came from South-Western Asia (Baptista et al. 1997) and has expanded through Europe due to anthropic corridors that facilitate invasions according to processes that have been known for more than 40 years (Udvardy 1969, Blackburn et al. 2009). Climate and biogeographic changes have also been proved to play a role in the expansion of this species (Ziska and Dukes 2014). For example, at the beginning of the twentieth century in Europe, increased rainfall weakened the North Balkan barrier and opened the way for the Eurasian collared-dove to expand from Turkey (Bernis et al. 1985) and its spread since the 1930s (Hengeveld and Van den Bosch 1991).
This species is one of the most successful invasive birds for its ability to travel long distances, and also due to its multi-brooding behavior (Bernis et al. 1985). It is believed to have arrived in the Iberian Peninsula from France (Bernis et al. 1985). In Spain, the first individuals were observed in Asturias in 1960, and rapidly spread along the Cantabrian coast by crossing the central part of the Spain and eventually colonizing the Mediterranean coast (Rocha-Camarero and de Trucios 2002). It has recently colonized North America (Hengeveld 1993, Fujisaki et al. 2010), and it has been postulated that this species expansion will continue in this continent, including Central and South America (Romagosa and Labisky 2000, Fujisaki et al. 2010).
The link between the Eurasian collared-dove and urban areas has been known for centuries (Hengeveld 1993, Romagosa and Labisky 2000). However, studies about its preferences in urban landscapes are scarce. In Great Britain, the Eurasian collared-dove is linked to populated zones, especially parks with trees of the genus Pinus, which provide food resources (Coombs et al. 1981). This species avoids areas where intensive agriculture takes place (Hudson 1972) but, in Spain, it inhabits other environments, such as the dehesas. Spanish dehesas are wooded open savannas where the main tree species are the holm oak evergreen tree (Quercus ilex) and cork oak tree (Quercus suber). This habitat is used by domestic herds (Rocha-Camarero and de Trucios 2002, Gaspar et al. 2009, Oviedo et al. 2013).
Analyzing the characteristics of the urban landscape used by this species would be useful to determine which factors influence the presence and abundance of the Eurasian collared dove in colonized towns and cities of Europe. Given the prediction that this species will keep expanding, this study aims to determine the most important factors that could favor Eurasian collared-dove expansion in its non-native environment. At the same time, we discuss the Eurasian collared-dove's microhabitat preferences in anthropic areas that, to our knowledge, have been poorly studied.
MATERIAL AND METHODS
The study was carried out in 46 localities within the geographical framework of the Valencia Region in eastern Spain (37°51' - 40°47' North and 02°09' West - 04°12' East; Table 1; Fig. 1), which occupies the central sector of eastern Iberian Peninsula and covers an area of 23,255 km2 (Murgui 2001, Aguilella et al. 2010).
The abundance and presence of the dove were determined using point counts in urban and periurban areas of the 46 pre-selected localities. Then, for five minutes after a one-minute settlement period, Eurasian collared-doves were counted once in 220-point counts without regarding distance. This number of point counts was randomly distributed in the set of localities according to their area and human population (Table 1). The distance between count points was at least 200 m. Counting periods started at sunrise and lasted for 3 h. Birds were recorded by both sight and sound (Hutto et al. 1986, Kasprzykowski and Golawski 2009). These point counts were carried out from autumn 2015 to late winter 2016 (before the 2016 breeding season) by one observer per point count. Point counts were carried out by four independent observers who followed the same census protocol (Table 1).
The geographical coordinates and altitude of each point count were obtained from SigPac (SigPac c2016), while mean minimum temperature and rainfall were obtained from World Climate Data (WCD c2016). The area and human population of each locality were obtained from the Spanish National Statistics Office (INE c2016).
Two types of areas in towns and cities were sampled based on type of buildings, traffic and pedestrian traffic. Urban areas were characterized by the presence of commercial and residential buildings with more than five floors, public green areas and continuous traffic (vehicles and pedestrians), while periurban areas were composed of detached and semidetached residential houses with private gardens, less traffic where vehicles and pedestrians are more occasional (Supplementary Material Photography 1).
Eighteen explanatory variables were estimated within a 50-meter radius from the point counts to detect those factors with a positive or negative association with the presence and abundance of the Eurasian collared-dove. Explanatory variables were grouped into three categories: resources, climate and urban structure, and human presence (Table 2).
Statistical analysis
Hierarchical partitioning (HP) analysis (Walsh and Mac Nally c2015) was used to assess the variables that best explained the presence and abundance of Eurasian collared doves. HP is used in ecological analysis to find environmental key factors for the presence or abundance of species (López-Iborra et al. 2011; Bernat-Ponce et al. 2018). HP fits all the possible regression models for a group of predictor variables, and estimates the independent effect (I) of each predictor on the response variable and percentage (I %) of the total I of the variable group (Chevan and Sutherland 1991, Mac Nally 2000, Quinn and Keough 2002). This work applied logistic regression and log-likelihood as goodness-of-fit measures for presence HP. For abundance data HP, Poisson regressions were used and R2 as the goodness-of-fit measure (Jongman et al. 1995). In both cases, the significance of effects was determined with 999 randomizations (Mac Nally 2002). For a significant contribution of the variable, the observed I had to be over the 95 percentile with Z-scores ≥ 1.64 (Hallstan et al. 2010). HP analyses were carried out using the "hier part" package (Walsh and Mac Nally c2015) in RStudio 3.5.1 (RStudio Team c2018).
A separate HP was performed for each group of variables. Then an HP analysis was performed which included the variables that were significant in the analyses per group. Since all groups included less than nine variables, HP analysis were not affected by the errors produced when more variables were analyzed (Walsh and Mac Nally c2015). The abundance HP analysis only included those localities (n = 33 localities: 119-point counts) with the presence of Eurasian collared-dove. A spatial term was included in all the models as a cubic function of the coordinates to control the effect of spatial autocorrelation (Legendre 1993). The directionality of the effect of each variable was obtained from univariate logistic or Poisson regression for the presence or abundance models, respectively.
To detect multicollinearity, the variance inflation factor (VIF) in each variable group was calculated using the 'vif function of the package 'car' (Fox and Weisberg 2019). Multicollinearity was considered high when VIF > 5 (Zuur et al. 2010). Altitude and mean minimum annual temperature correlated highly (VIF > 5) in the presence and abundance analyses and, consequently, altitude was not included in the HP analyses. Considering that mean minimum annual temperature was not the same at the same altitude at different latitudes, we believed that it had more biological sense to include temperature in the analyses.
RESULTS
The species was detected in 33 out of 46 of the studied localities and 54.1 % of point counts (n = 220). The number of Eurasian collared doves counted in points with presence ranged between 0 and 36.
The HP analysis showed that at least one variable in each group had a significant positive association with the presence of the Eurasian collared dove (Table 3). The presence HP analysis showed in the resources group that the positive and significant variables were parks, schools, and exotic vegetation. In the climate and urban structure group, only mean minimum annual temperature had a positive and significant relation. In the human factors group, town/city population was the only significant and positive variable. Of all the analyses, only restaurant presence (resources), and pedestrian traffic (human presence) had a negative and significant relation (Table 3). The combined HP analysis revealed that seven variables and the spatial term were significant. Apart from the spatial term, population in towns/cities and park presence explained approximately the 30 % of the deviance on the presence of Eurasian collared dove (Table 4). Exotic vegetation and restaurants had relatively small but significant contributions according to their results in final model.
The same analysis in relation to abundance identified three significant variables with positive associations with the Eurasian collared dove in the resources group (water, native vegetation, exotic vegetation) (Table 3). The urban area variable (climate and urban structure) was negative (Table 3) while the variables related with the human presence group had not effect on the abundance of Eurasian collared-dove. The final combined HP showed that exotic and native vegetation gave the highest independent contributions and together explained 35 % in deviance (Table 4).
DISCUSSION
In our study area Eurasian collared-doves inhabit urban and periurban landscapes, which agrees with previous studies in other countries (Coombs et al. 1981, Cramp 1985, Hengeveld 1993, Blair 1996, Romagosa and Labisky 2000). Further, it has a clear preference for green areas, such as parks and urban areas with native and exotic vegetation. This preference can be explained by the species nesting in trees and feeding especially on the ground of green areas (Coombs et al. 1981). Urban parks are usually associated with tree presence. Furthermore, exotic vegetation revealed a positive association with both dove presence and abundance and it is often composed by trees which are suitable for nesting and roosting. This kind of vegetation is always linked to artificial habitats managed by humans, especially in private gardens that are part of the suburban environment.
The coastal zones of our study area have a climate with warm temperatures during summer and in winter these temperatures are mild (Rivas-Martínez 1987, Pérez-Cueva 1994).
However, the mean minimum temperature decreases when altitude increases and both variables, altitude and-mean minimum annual temperature were highly negatively correlated. The expansion of the Eurasian collared-dove is associated with specific environmental gradients, such as temperature and precipitation (Bled et al. 2011). Our study showed that the presence of the Eurasian collared dove is restricted by the mean minimum annual temperature in the studied localities, and is less likely to establish in low temperatures (Altwegg et al. 2008). Meteorological varia-bles, such as precipitation, is a negative predictor for its presence in USA (Fujisaki et al. 2010, Bled et al. 2011, Scheidt and Hurlbert 2014). However, this variable was not important in our study, which was characterized by a Mediterranean climate, where precipitations are scarce and irregular, and temperatures are warmer (Pérez-Cueva 1994). Thus, in warmer and drier areas, precipitation has no effect and, thus, temperature could have a stronger influence on presence.
Birds respond differently to landscape characteristics (Villegas and Garitano-Zavala 2010). In urbanized areas, green areas are associated with parks and avenues (Fernández-Juricic 2000, González-Oreja et al. 2012, Aram et al. 2019), while in suburban areas there are not only parks, but also many private gardens. Therefore, the surface covered by vegetation increases in these suburban areas, where vegetation is a mixture of both native and exotic vegetation (Bernat-Ponce et al. 2018) which would explain the relevance of these areas for the species abundance. Further, the presence of trees provides the necessary places for Eurasian collared-doves to build nests (Coombs et al. 1981).
Cities and towns with a bigger human population increase the probability of finding Eurasian collared doves. A larger human population is associated with larger cities and towns, which also means more periurban areas and number of schools (Giménez García et al. 2018), which are variables that positively associated with the presence of doves, respectively. Schools are potential resources suppliers in urban areas because of the food scraps left by children in playgrounds. Moreover, big cities contain more parks, avenues and other types of green areas (Russo and Cirella 2018). However, the probability of finding Eurasian-collared doves in areas with heavier pedestrian traffic and restaurants is reduced. Areas with a heavier pedestrian traffic coincide with those areas where buildings, restaurants and commercial centers are more concentrated, which diminish available area, increase disturbances, and limit access to natural food and shelter resources. The negative link of restaurants with Eurasian collared-doves presence could be a consequence of increased pedestrian traffic.
The success of invasive species, especially in arid areas like the Mediterranean Region, is limited by lack of water (Macmillen 1962, George 1975, Taylor et al. 2006). Thus, an additional water resource, especially during warm seasons, must supply their daily water requirements. Our results showed that water presence had a positive association with Eurasian collared-dove abundance. Usually most detached and semidetached houses in periurban areas have swimming pools, and these constructions constitute a new source of water. Water presence must be favorable during warm seasons because it facilitates access to water. Therefore, the positive selection of periurban areas in relation to the negative association with urban areas might be related to presence of water resources, the increase in green areas and the lesser impact of pedestrian traffic. Although Eurasian collared-doves also occupy urban areas, our results agree with previous studies (Hengeveld 1993, Romagosa and Labisky 2000) as suburban areas provide the necessary resources for species to avoid urban areas (Blair 1996, Chapman and Reich 2007). Furthermore, the resources offered by this environment, habitat and food resources would facilitate its use as new corridors for invasive species (Ziska and Dukes 2014).
In conclusion, our study indicates the key factors, such as parks and number of inhabitants, for the presence and, exotic vegetation and water availability, for the abundance of a species that is expanding in Europe and is invasive in America. Our results can be useful to predict which non colonized areas will be more plausibly colonized by the Eurasian collared dove. Furthermore, our results could be used to predict the parts of invaded areas that will harbor more individuals, which would be key for managing the Eurasian collared dove in its native and its non-native distribution area.