Space invaders : House Sparrow densities along three urban-agricultural landscapes

One of the most successful avian urban invaders is the House Sparrow (Passer domesticus). Here, we examine House Sparrow numbers in the urban-agricultural landscape of three cities, one where the sparrow is native (Barcelona) and two where it is nonnative (Los Angeles, Mexico City), comparing between five different land classes based on land-use, location relative to urban center, and urbanization intensity. We estimated sparrow densities using distance-sampling corrected procedures and compared sparrow numbers between the three cities and between the land classes within each city. This two-year survey showed the highest density of sparrows occurring in Mexico City and the lowest in Los Angeles. In both cities where the sparrow is nonnative, we recorded the greatest abundances in highly urbanized areas. In Barcelona, where the sparrow is native, abundances did not differ among land classes. We also related sparrow numbers to specific environmental variables, describing vegetation structure, urban infrastructure, and human activity, using classification and regression trees. We found several variables related to differences in House Sparrow abundances among the three cities. For Los Angeles and Mexico City, the most important variables explaining House Sparrow numbers were urban infrastructure (i.e., maximum building height, high urbanization) and/or human activity (i.e., pedestrians), all of which related positively with sparrow abundance. For Barcelona, sparrow numbers showed both positive and negative relationships with vegetation (i.e., tree abundance, maximum tree height, maximum herbaceous plant height) and human-related variables (i.e., cars, pedestrians). Although the scope of this study is geographically limited to three cities, our results suggest that House Sparrow numbers between cities and land classes areas vary considerably, and are likely determined by a complex array of variables. Our study provides some insight to the growth and decline of House Sparrow populations in their nonnative and native ranges, respectively. Alerte aux envahisseurs: Densités de moineaux domestiques au long de trois paysages agrico-urbains RÉSUMÉ. Un des envahisseurs aviaires ayant le plus de succès est le Moineau Domestique (Passer domesticus). Ici, nous examinons les densités dans les zones agrico-urbaines de trois villes, une où le moineau est natif (Barcelone) et deux où l'espèce est considérée envahissante (Los Angeles, Mexico), en comparant entre cinq classifications terrestres établit selon l'utilisation du territoire, localisation relative au centre urbain, et l'intensité d'urbanisation. Nous avons estimé les densités de moineaux en utilisant la méthode d'échantillonnage selon la distance corrigée et avons comparé le nombre de moineau entre les trois villes, ainsi qu'entre les classifications terrestres dans chaque ville. Cet échantillonnage de deux ans a démontré la plus grande densité à Mexico et la plus basse à Los Angeles. Chez les deux villes où le moineau est considéré comme envahissant, les densités les plus élevées furent enregistrées dans les zones très urbanisées. À Barcelone, où le moineau est natif, les abondances ne démontrent pas de différences entre les classifications terrestres. Nous avons aussi relié le nombre de moineau à des variables environnementales spécifiques, décrivant la structure de la végétation, les infrastructures urbaines, les activités humaines, en utilisant les méthodes de classification et régression arborescentes. Nous avons trouvé plusieurs variables reliées à des différences chez le Moineau Domestique parmi les trois villes. À Los Angeles et Mexico, les variables les plus importantes expliquant les densités du Moineau Domestique sont les infrastructures urbaines (c.-à-d., hauteur maximum de l'édifice, haut degré d'urbanisation) et/ou l'activité humaine (c.-à-d., piétons), tous démontrent un lien positif avec l'abondance de moineaux. En ce qui concerne Barcelone, des liens positives et négatives sont associés avec la végétation (c.-à-d., abondance d'arbre, hauteur maximum d'arbre, hauteur maximum des plantes herbacées) ainsi que les variables anthropogéniques (c.-à-d., voitures, piétons). Bien que la portée de cette étude soit limitée géographiquement à trois villes, nos résultats suggèrent une considérable variation du nombre de Moineaux Domestiques entre les villes, ainsi que les classifications terrestres, et est probablement déterminée par un éventail complexe de variables. Notre étude apporte un aperçu de la croissance et du déclin des populations de Moineaux Domestiques dans leurs aires de répartition respectives, qu'elle soient invasives ou naturelles. .


INTRODUCTION
Urbanization is a process that drastically alters pre-existing systems and concentrates significant amounts of resources to satisfy the needs of modern life (McDonnell andMacGregor-Fors 2016, Seto andRamankutty 2016).This generates enormous environmental impacts at various geographic scales, ranging from the local to the global (Seto and Ramankutty 2016).In fact, urbanization is related to a broad range of long-term cascading ecological impacts, especially for native flora and fauna (Czech et al. 2000, Berkowitz et al. 2003, Jetz et al. 2007, Grimm et al. 2008, Pickett et al. 2008, Aronson et al. 2014, Seto and Ramankutty 2016).This can importantly affect the ecology of an entire geographic region, and, as a result, threaten native species (Emlen 1974, Beissinger and Osborne 1982, Czech et al. 2000, Van Rensburg et al. 2009, Maxwell et al. 2016).As urbanization poses novel pressures to biodiversity, some species are able to adapt, and even thrive, within urban centers, while those unable to do so decline in number (Chace and Walsh 2006, Evans et al. 2009, Fischer et al. 2015).In particular, several bird species have been shown to successfully exploit areas of human development, to the degree that they have established populations beyond their native geographic ranges (Blair 1996, Shochat 2004, Kark et al. 2007, Blackburn et al. 2009).
The House Sparrow (Passer domesticus) has been introduced anthropogenically worldwide, making it one of the most prolific and successful introduced bird species in the world (Anderson 2006, Lowther and Cink 2006, Aronson et al. 2014).Native to Eurasia, the House Sparrow has been dominant and successful at exploiting urban areas of North America, often repelling potential competitor avifauna (Gavareski 1976, Gowaty 1984, Clergeau et al. 1998, White et al. 2005, MacGregor-Fors et al. 2010).The success of the House Sparrow as an urban exploiter species in North America has been mainly attributed to four natural history traits: (1) it is a dietary generalist that feeds on grains, insects, fruits, and even human litter (Gavett andWakeley 1986, Kimball 1997); (2) it has colonial-communal nesting strategies that allow for rapid proliferation upon invading new areas (Kalinoski 1975, McGillivray 1980, Gowaty 1984); (3) it can effectively expand its ranges throughout human-altered landscapes (Kark et al. 2007); and (4) it is a territorial species that actively destroys the nests of other species (McGillivray 1980, Gowaty 1984, Kimball 1997).
Since its introduction in the northeastern United States in the 1800s, the House Sparrow has steadily expanded its range across human-disturbed landscapes, establishing large populations in urban centers (Johnston and Selander 1973, Baughman 2003, Lowther and Cink 2006).Although records of its arrival to Mexico are sparse, it presumably arrived in Mexico City by the 1930s (Wagner 1959).We have no information about its presence in Central America until the 1970s (Davis 1972).Although potential habitat for the species in Mexico only excludes the Southern Baja Peninsula, northwestern Mexico, and Yucatan Peninsula (Navarro and Peterson 2007), its presence is closely associated with that of human development, especially from within the heavily developed urban centers out to isolated farmlands.They are rarely found in undisturbed areas (Summers-Smith et al. 2017).In its native Eurasian range, House Sparrow populations have experienced considerable declines as a result of a suite of complex anthropogenic factors (Robinson et al. 2005, Brichetti et al. 2007, Ferrer et al. 2013, PECBMS 2013, 2016).
Previous studies have investigated House Sparrow success in relation to specific habitat features.Recent studies performed in Mexican cities, i.e., Mexico City, Morelia, that examined the effects of built cover and building height showed that avian communities in areas invaded by the sparrows had lower species richness in comparison to those of noninvaded areas (Ortega-Álvarez and MacGregor-Fors 2009).This indicates that House Sparrow success in invading nonnative areas, mainly through its synergistic interactions with human activity, is closely related to changes in the diversity, composition, and structure of native bird communities (MacGregor-Fors et al. 2010).More generally, another study found that urban residential bird communities can experience lower species richness but greater abundances, with those in commercial areas favoring generalist and urban exploiter species (Ortega-Álvarez and MacGregor-Fors 2009).Yet, the House Sparrow is comparably less abundant in other Mexican cities, e.g., Xalapa (Escobar-Ibáñez and MacGregor-Fors 2016).Also, House Sparrow abundances have been found to vary between the four urban land-uses of Mexico City, with the highest numbers recorded in residential-commercial areas and the lowest in green areas (Ortega-Álvarez and MacGregor-Fors 2011).This study showed that the variables determining House Sparrow numbers vary among urban land-uses, with some variables consistently positively related, e.g., number of passing pedestrians or building height, and others consistently negatively related, e.g., tree density and cover.Similar patterns of habitat relationships have also been reported in this species' native European distribution (Mason 2006, Chamberlain et al. 2007, Tratalos et al. 2007, Murgui 2009, Murgui and Macías 2010, García-Rodríguez 2011, Šálek et al. 2015), where variables associated with vegetation, but not specifically with urbanization, have been shown to drive their numbers.The declines of native House Sparrow populations have been recorded in both urban and agricultural areas, although the processes are seemingly different.For instance, its populations have experienced considerable declines since the replacement of horses by cars in the 1920s, followed by nonurban population declines associated with the reduced availability of both grain and invertebrates due to the industrialization of farmland practices, e.g., heavier use of pesticides and herbicides, reduced grain spillage, improved grain storage, from the 1980s onward (Summers-Smith 2003).The renewed decline of its urban populations since the 1990s has been much more complex and highly variable with the reduced number of shrubs and nesting sites, increased pesticides use in greenspaces, and increased concentration of lower socioeconomic status areas all being related to lower sparrow numbers (Summers-Smith 2003, Balmori and Hallberg 2007, Shaw et al. 2008, Murgui and Macías 2010, García-Rodríguez 2011).However, it is unclear which of these reasons are directly causing the decline (De Laet and Summers-Smith 2007).
The success of the House Sparrow in North America has been extensively studied in terms of its multiple introductions and subsequent differentiation among populations (Johnston and Selander 1964, Brown and Wilson 1975, Baughman 2003, Lowther and Cink 2006).Yet, how the numbers of this sparrow might differ along urban-agricultural landscapes where it is native and nonnative is not well understood.By setting up our study in a comparative framework that includes different urbanagricultural land classes, based on land-use, location relative to urban centers, and urbanization intensity, we provide insight to the paradox of the growth and decline of House Sparrows in their nonnative and native ranges, respectively.Moreover, there is an important lack of knowledge regarding the environmental variables related to their numbers along urban-agricultural landscapes.Thus, the aim of our study is to document how the House Sparrow numbers differ in five different land classes, i.e., intra-urban high urbanization, intra-urban low urbanization, peri-urban high urbanization, peri-urban low urbanization, and agricultural, in three cities, one located within its native distribution, Barcelona; two in its North American invasion, Los Angeles, Mexico City.Specifically, we ask the following: (1) how do House Sparrow densities differ between three cities, (2) how do their abundances differ among land classes, defined by geographic location and urbanization intensity, and (3) which variables, e.g., vegetation structure, urban infrastructure, and/or human activity, are related to House Sparrow numbers across the studied urban-agricultural landscapes?

Study areas
This study was conducted in three urban-agricultural landscapes: Los Angeles (California, United States), Mexico City (Mexico), and Barcelona (Catalonia, Spain; Fig. 1, Table 1).We worked in urban-agricultural landscapes, including urban agglomerations rather than governmentally delimited cities, to assess House Sparrow abundances.Focusing on governmentally delimited urban boundaries did not allow us to consider House Sparrow abundance variations when comparing urban cores and periurban areas.

Focal land classes within urban-agricultural landscapes
We established our study sites in urban, e.g., residential or commercial, and agricultural land-uses, because the House Sparrow's native and nonnative distributions are generally limited to these types of areas (Cramp 1998, Anderson 2006, Lowther and Cink 2006, Murgui 2009).Specifically, for each urban area, we assessed House Sparrow numbers between five land classes: (1) intra-urban high urbanization; (2) intra-urban low urbanization; (3) peri-urban high urbanization; (4) peri-urban low urbanization; and (5) agricultural.We categorized these land classes based on three dichotomous environmental conditions that have been shown to drive bird distributions in urban areas: land-use, i.e., urban, agricultural; location relative to urban centers, i.e., intra-urban, peri-urban (López 2010, Puga-Caballero et al. 2014); and urbanization intensity, i.e., low, high urbanization (Chace andWalsh 2006, Evans et al. 2009).
We selected sites within each land class using high quality satellite imagery based on land-use, location relative to urban centers, and urbanization intensity variables.We established a ~50% built cover threshold to differentiate between sites with low and high urbanization, i.e., 0-50% = low urbanization; 51-100% = high urbanization.To distinguish between intra-and peri-urban areas, we used methods to delineate the limit of the peri-urban area of an urban agglomeration (MacGregor-Fors 2010).Briefly, this method calculates the representative area in which the urban core intermingles with adjacent nonurban systems, representing an ecologically meaningful ecotone for birds (Puga-Caballero et al. 2014).Regarding agricultural sites, we selected them in croplands with few built elements (0-30% built cover) along the city outskirts.We verified that our sites matched the correct parameters of the land class by manually quantifying built cover, i.e., all hard and relatively impervious surfaces, including buildings and paved streets, on the satellite imagery.

Point-count surveys
We surveyed House Sparrows over a two-year period, for the winter seasons (December-January) of 2012-2013 and 2013-2014, and the breeding seasons (May-June) of 2013 and 2014.We used limited radius point-counts, at which we recorded all House Sparrows seen or heard within a 50-m radius in a five-min period (following Ralph et al. 1993), from 07:00 to 11:00.In order to calculate distance-sampling corrected densities, we measured the radial distances from each recorded sparrow to the observer using laser rangefinders (Bushnell Yardage Pro Sport 450, Table 1.Geographical and population traits for the three studied urban-agricultural landscapes.Overland Park, KS).We established our point-count locations at a minimum separation distance of 200 m to avoid issues of pseudoreplication (Ralph et al. 1993(Ralph et al. , 1995)).We performed 10 point-counts within one-or two-morning periods, depending on weather conditions and logistical reasons, for each land class for the three surveyed cities.This resulted in a total of 150 pointcounts per survey season in all three cities.Each site was repeated once per season, for a total of 600 point-count repetitions.

Environmental variables
In order to assess which environmental variables were related to House Sparrow numbers in the three surveyed cities, we reviewed previous studies relating the presence and abundances of urbanexploiter species to environmental variables (e.g., Kark et al. 2007, Evans et al. 2009, Murgui and Macías 2010, Ortega-Álvarez and MacGregor-Fors 2009, 2011).We identified 12 specific variables that describe vegetation structure, urban infrastructure, and human activity, which we sampled within the 50-m radius area at which we surveyed the sparrows: (1) tree cover, (2) tree abundance, (3) maximum tree diameter at breast height (DBH), (4) maximum tree height, (5) shrub cover, (6) maximum herbaceous plant height, (7) number of lamp poles, (8) number of telephone and electric power poles, (9) number of telephone and electric power cables, (10) maximum building height, (11) number of passing pedestrians in five min, and (12) number of passing cars in five min.

Data analysis
We estimated House Sparrow densities using the distancesampling corrected procedures in Distance 6.0 (Thomas et al. 2010).This software estimates density by calculating the detection probability of individuals (ind/ha) at increasing distances from the observer, and then by standardizing the detection rates along the concentric surveyed area (Buckland et al. 2001).Although detectability can differ among land classes, our data were not sufficient for performing distance-corrected analyses by "land class" but were sufficient for "city" and "season."With its capabilities of correcting for detection probability as well as standardizing by the surveyed area, we believe that the use of Distance is crucial for avoiding important biases in our results, assuming pooled "land class" detection probabilities.We used the key function/series expansion combination that best fitted our data for each city/season based on the Akaike Information Criterion (AIC) and coefficient of variance.To determine statistical differences between the density estimations between the three cities, we compared their 84% confidence intervals, following MacGregor-Fors and Payton (2013).Given that overlapping 95% confidence intervals fail in detecting statistical differences with an α ≤ 0.05, we used 84% confidence intervals, which have been shown to approximate overlapping confidence intervals with an α ≈ 0.05 for both symmetrical and asymmetrical intervals; for House Sparrow density estimations, confidence intervals are asymmetrical, from normal and log-normal distributions (Payton et al. 2003, MacGregor-Fors andPayton 2013).Hence, when the 84% confidence intervals for the density estimations of two groups did not overlap, we considered them to be statistically different from one another; when their confidence intervals did overlap, we considered them to not be statistically different from one another.
For comparisons between land classes within cities, we were unable to use distance-sampling corrected density estimations because of the low numbers of recorded House Sparrows in some of the land classes; in agricultural areas, sparrows were found in more sparsely distributed and tighter groups.Thus, we used the observed House Sparrows individuals/point-count, rather than density estimations, for making comparisons between land classes within cities.We performed general linearized mixed models (GLMM) of House Sparrow numbers using Poisson family and a log-link function (package lme4; Bates et al. 2016) in R (R Core Team 2017).We used House Sparrow individuals/point-count as the response variable and the land class, season, and year as predictor variables.We selected the most parsimonious model based on the Akaike information criterion (AIC).We considered point-counts as a random variable and used a dummy data variable to control for model overdispersion.We used the best model to perform pairwise Tukey comparisons for establishing statistical differences in the abundance of House Sparrows per point-count between land classes within the urban-agricultural landscapes.
Finally, to assess the relationships between House Sparrow abundances per point-count and environmental variables at each city, we performed two analyses.Given that exploratory analyses allowed us to identify potential nonlinear relationships between our independent variables and House Sparrow numbers in the three studied cities, we first performed a generalized additive mixed model (GAMM) considering "city" as a random factor.GAMMs are similar to GLMs in having different error structures and link functions, but the shape of relationships is not specified by explicit functional forms; instead, nonparametric "smoothers" to describe the specified nonlinear relationship, as well as random (nonfixed) factors/effects (Crawley 2013).Based on the set of significantly related variables identified through the GAMM, we performed classification and regression tree (CART) analyses using R (R Core Team 2017) for the three studied cities, considering "city," "land class," and "season" as independent variables, besides the other eight independent variables identified by the GAMM.This exploratory general CART showed "land class" to be, hierarchically, the most important variable explaining variance in House Sparrow numbers in the three cities, followed by "city," passing pedestrians, maximum building height, passing cars, maximum tree height, maximum herbaceous plant height, tree abundance, and poles.Thus, we performed three additional CARTs, one per city, considering the eight independent variables identified by the GAMM, as well as "land class."We used CARTs because they allow the interpretation of datasets where complex nonlinear relationships occur between the set of response and predictor variables (De'ath and Fabricius 2000).These trees are analogous to multiple regression models, specifically to those using the forward selection of predictor variables (Crawley 2007(Crawley , 2013)).By using binary recursive partitioning to identify the successive critical threshold values for a set of predictor variables, CARTs split the dataset of the response variable in a dichotomous and hierarchical manner; as a result, the rank order of only one numeric predictor variable is identified at each split of the tree, showing not only relationships, but also scenarios under which independent variables can relate differently to the dependent one (Palomino and Carrascal 2007).We used the function "rpart" for R (Therneau et al. 2015, R Core Team 2017) because it carries out analysis of variance with the two-level variables associated with each split, which makes it more suitable for anticipating the results of simplified models compared to functions like "tree" (Crawley 2013).This procedure is well suited for assessing relationships between a dependent variable (in this case, House Sparrow abundance per point-count) and multiple independent variables, which can be a mix of continuous and categorical variables (in this case, the 12 environmental variables sampled within the 50-m radius area of where sparrows were surveyed) even if they are highly correlated (Andersen et al. 2000, Jackson andBartolome 2002).

RESULTS
After quantifying built cover using the satellite imagery for all survey sites, we confirmed that our defined land classes within the urban-agricultural landscapes are unique, based on land-uses, location relative to urban centers, and urbanization intensity (Table 2).Built cover in the surveyed sites was very low in agricultural fields (average 5.7% ± SE 1.6), with all sites showing values < 30%.For low urbanization sites, average built cover was 13.43% (± SE 2.1); all values, except for one site in Los Angeles, i.e., intra-urban low urbanization, 55.1% built cover, had < 50% built cover.For high urbanization sites, average built cover was 77.3% (± SE 2.1), with only two peri-urban sites having < 50% built cover in Barcelona, i.e., 41.9% and 45.9% built cover.
House Sparrow densities differed significantly between the three surveyed cities (Fig. 2).Within each city, we found no difference between seasons, with the exception of the 2014 breeding season in Mexico City.The number of House Sparrows per hectare in Mexico City (84% CI range for both seasons in both years: 9.5-33.3sparrows/ha) was significantly higher when compared to the numbers estimated for Los Angeles (84% CI range for both seasons in both years: 0.4-3.1 sparrows/ha) and Barcelona (84% CI range for both seasons in both years: 3.2-7.1 sparrows/ha).We found differences in House Sparrow abundances per pointcount when comparing between land classes within both Los Angeles and Mexico City; however, no differences were found when comparing between land classes within Barcelona (Tables 3 and 4).We recorded the highest House Sparrow abundances in the intra-urban high urbanization land class within both Mexico City and Los Angeles (Fig. 3).In Mexico City, we recorded a high number of House Sparrows in the peri-urban high urbanization land class, with only a few or no individuals in the remaining land classes.Furthermore, in Mexico City we recorded no House Sparrows in agricultural areas.To compare sparrow densities between cities, we considered two groups statistically different from one another when their 84% confidence intervals did not overlap, while we considered them not statistically different from one another when their confidence intervals did overlap.
Results of the GAMM show that 8 out of the 12 assessed variables were significantly related with House Sparrow numbers, showing the importance of both "city" and "land class," as well as no effect of "season" (Table 5).Results of the CARTs for each city using the eight variables identified by the GAMM as well as "land class" show that the most important variables for explaining sparrow number differences, being the earliest splits in the trees, were maximum tree height in Barcelona, maximum building height in Los Angeles, and "land class" (separating agriculture and low urbanization from high urbanization sites) in Mexico City (Fig. 4).CARTs show that both vegetation and human-related variables explain shifts in House Sparrow numbers in Barcelona, while high urbanization scenarios, as indicated by results of the GLMM (Tables 3 and 4), show a positive relationship with their     City all showed positive relationships with House Sparrow abundances.Meanwhile, the only vegetation variable selected by the CARTs for these cities (i.e., tree abundance) was negatively related.In contrast, the CART for Barcelona shows a complex array of positive and negative relations with both vegetation and human-related variables.

DISCUSSION
Our study, focused on House Sparrow abundances in different land classes along three urban-agricultural landscapes, shows the highest densities in Mexico City and the lowest densities in Los Angeles.Thus, depending on specific habitats and scenarios, House Sparrow numbers can resemble those of both urban exploiters/dwellers (Blair 1996, Fischer et al. 2015) and urban adapters/users (Blair 1996, Fischer et al. 2015), varying between cities located outside of its native range and even between urbanagricultural land classes within those cities.This agrees with a number of other studies that have shown that sparrow numbers can vary in both their native and nonnative distribution despite their reputation for high adaptability (Johnston and Selander 1964, Johnston and Selander 1973, Kendeigh 1976, Martin et al. 2004) and exploitation of urban-and agriculture-related resources (Blair 1996, Herrando et al. 2012, Fischer et al. 2015).
We recorded a high variation in House Sparrow densities in the two cities located outside of the sparrow's native distribution.Regarding House Sparrow numbers in the United States, data from the North American Breeding Bird Survey (Sauer et al. 2012) show that House Sparrows in the Los Angeles area are (1) scarce in the northern parts of the urban continuum (~1-3 sparrows/2.5 hr of point-count surveys, at which every bird seen or heard within a 0.25-mi radius in a 3-min period was recorded) and ( 2) moderately abundant in the southern parts (~10-30 sparrows/2.5 hr).This pattern in sparrow density shows that its numbers and invasive success are not homogeneous throughout urban territories and are not solely due to latitude or productivity but possibly more to specific landscape-and local-scale environmental variables (e.g., Kark et al. 2007).In such a case, this would resemble some of the patterns observed in its native areas (Murgui 2009(Murgui , Šálek et al 2015)).
Regarding the House Sparrow's native range, our density results for Barcelona are similar to those found by Murgui (2006) in another Spanish city (Valencia), who reported densities of 0.9-8.5 sparrows/ha.Although we did not find differences in sparrow abundance between the studied land classes, other ecological patterns, such as urban metapopulation dynamics, could be driving this result (Chávez-Zichinelli et al. 2010).
When analyzing House Sparrow numbers at a finer scale, within the land classes of each urban-agricultural landscape, the highest abundance values for both Mexico City and Los Angeles were from the high urbanization land classes, for both intra-and periurban locations.This result agrees with previous studies suggesting that this sparrow can take advantage of a wide array of resources in areas where other species are not present because of the lack of suitable habitat components and the presence of urban-related hazards (Blair 1996, Kark et al. 2007).Contrasting with our observations in the agricultural sites of Barcelona, we only recorded a few individuals in the agricultural land class for Los Angeles and no individuals in the agricultural land class for Mexico City.The latter agrees with our nonsystematic observations in both areas from over the past four years: that the sparrows are often found in lower densities in these agricultural areas (I.M-F. and J. G-H.L., personal observation).
Given that 8 of the 12 environmental variables we assessed showed significant relationships with House Sparrow numbers, it is clear that vegetation structure, urban infrastructure, and human activity are related to changes in their abundance, consistent with previous studies that have reported some of these variables related with House Sparrow numbers in urban areas (Murgui 2009, MacGregor-Fors et al. 2010).Yet, the types of variables and scenarios related to variations in House Sparrow numbers differed importantly among cities.The most relevant variables in explaining variation in House Sparrow abundances in Los Angeles and Mexico City were those related to urban infrastructure, i.e., maximum building height and high urbanization, and/or human activity, i.e., passing pedestrians, all of them showing a positive relationship with House Sparrow numbers in the two cities.Furthermore, in Mexico City, tree abundance, the only vegetation variable explaining some variance of the recorded House Sparrow abundance, showed a negative relationship.In contrast, our findings for Barcelona show a complex array of scenarios under which "land class" did not play a crucial role, with no urban infrastructure variable showing an important relationship with House Sparrow numbers.Vegetation variables played the most important role, explaining variance in the sparrow's abundances, both positively and negatively.For instance, the scenario in which we recorded lower average House Sparrow numbers per point-count (0.3) involved > 96, > 7 m trees.
On the other end, the scenario in which we recorded highest average House Sparrows per point-count (9.0) has > 56, 8-10 m trees.This further demonstrates the species' close association with highly developed and human-frequented sites, especially in nonnative areas where such features represent important foraging and nesting resources (Anderson 2006).
Although this study only considers three urban-agricultural landscapes, our results suggest that the densities and abundances of House Sparrows located outside of their native ranges are determined by environmental variables that are unique to different cities and to different land classes, locations relative to urban centers, and urbanization intensities.There are a number of further questions this study raises related to the potential effects of House Sparrow density on the native biota, as well as to its role as a driver of biodiversity in urban-agricultural landscapes.
Comparative studies focused on the individual (e.g., body condition, personality), population (e.g., trends in distribution and growth, intra-specific interactions, genetics), and community approaches (e.g., diversity, interspecific interactions) related to the House Sparrow in its native and nonnative areas would provide further insight for the development of population management tools in addressing two growing concerns: (1) their recent declining native populations in Europe (Murgui 2006, De Laet andSummers-Smith 2007) and ( 2) the problems they pose to the local biodiversity in invaded areas (Grussing 1980).

Fig. 1 .
Fig. 1.Map of the three studied urban-agricultural landscapes.Gray polygons of insets correspond to the urban agglomerations.

Fig. 2 .
Fig. 2. House Sparrow (Passer domesticus) distance-sampling corrected density estimations for Los Angeles, Mexico City, and Barcelona.To compare sparrow densities between cities, we considered two groups statistically different from one another when their 84% confidence intervals did not overlap, while we considered them not statistically different from one another when their confidence intervals did overlap.
Measurement of on-screen urban agglomerations using high quality satellite imagery.‡ Estimated within the urban agglomerations (United Nations 2015). †

Table 2 .
Built cover (average ± SD) of the sites in which we surveyed House Sparrows (Passer domesticus) by city and "land class."

Table 3 .
General linearized mixed models for House Sparrow (Passer domesticus) abundances in the three studied urban-agricultural landscapes.The most parsimonious model was selected based on the lowest AIC score.

Table 4 .
Pairwise Tukey comparisons for House Sparrow (Passer domesticus) abundances in the five land classes of the three studied urban-agricultural landscapes.

Table 5 .
Generalized additive mixed model for House Sparrow (Passer domesticus) abundances in the three studied urbanagricultural landscapes, considering "city" as a random factor.
Compared to Barcelona, we recorded greater average House Sparrow densities in Mexico City (~4 times more), while we found lower average House Sparrow densities in Los Angeles.Although there is a limited knowledge of House Sparrow density patterns along the Mexican territory, our density results are similar to those from another study conducted in Mexico City (~7 sparrows/ point-count; Ortega-Álvarez and MacGregor-Fors 2011).However, other studies have reported different House Sparrow densities in other Mexican urban areas, with higher values in Morelia (~20 sparrows/point-count; MacGregor-Fors et al. 2010) and lower values in Xalapa (~0.5 sparrows/point-count; J. F. Escobar-Ibáñez, personal communication).Given that the negative relationship between House Sparrow numbers and native bird species richness has only been assessed in a city where the sparrow is highly abundant, i.e., Morelia, its potential adverse effects could be density dependent, and this observed pattern may not necessarily hold across its nonnative range, although future work must be done to feasibly demonstrate this.