Next Article in Journal
Comparison of 5 Versus 7-Day Ovsynch + Progesterone Releasing Intravaginal Device Protocols (PRID) and a Modified G7G with an Option of Heat Detection Protocol for 1st Service in Lactating Dairy Cows
Next Article in Special Issue
In Its Southern Edge of Distribution, the Tawny Owl (Strix aluco) Is More Sensitive to Extreme Temperatures Than to Rural Development
Previous Article in Journal
Management of Low Birth Weight in Canine and Feline Species: Breeder Profiling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 10
10.3390/ani11102954
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Colonization of Urban Habitats: Tawny Owl Abundance Is Conditioned by Urbanization Structure

by
Nerea Pagaldai
1,*,
Juan Arizaga
1,
María V. Jiménez-Franco
2 and
Iñigo Zuberogoitia
1,3
1
Department of Ornithology, Aranzadi Sciences Society, Zorroagagaina 11, 20014 Donostia-San Sebastian, Spain
2
Centro de Investigación e Innovación Agroalimentaria y Agroambiental (CIAGRO-UMH), Department of Applied Biology, Miguel Hernández University of Elche, 03202 Elche, Spain
3
Estudios Medioambientales Icarus S.L. C/San Vicente, 8. 6 Planta, Dpto 8. Edificio Albia I., 48001 Bilbao, Spain
*
Author to whom correspondence should be addressed.
Animals 2021, 11(10), 2954; https://doi.org/10.3390/ani11102954
Submission received: 6 August 2021 / Revised: 10 October 2021 / Accepted: 11 October 2021 / Published: 13 October 2021
(This article belongs to the Special Issue Owls' Responses to Environmental Challenges)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Alteration of natural habitats due to urbanization is an increasing concern worldwide. Some species, including owls, can exploit this novel environment, although the consequences at the population level have not been described. In this study, we analyze the effect of different urban variables on tawny owl (Strix aluco) population abundance. At the local scale, forest and urban cover, as well as the clumpiness index, affected tawny owl abundance. At the landscape scale, its abundance decreased in complex-shaped urban patches and when distance between them was greater. Urban habitats do not substitute natural habitats in terms of abundance, but the species can easily colonize patchy urban habitats.

Abstract

Natural habitats are being altered and destroyed worldwide due to urbanization, leading to a decrease in species abundance and richness. Nevertheless, some species, including tawny owls, have successfully colonized this novel habitat. Consequences at the population level have not been described; thus, our main objective was to describe the effects that urban structure have on the tawny owl population at local and landscape levels. Data were obtained from 527 survey points over 7 months in a large-scale owl survey in the Basque Country (northern Spain) in 2018. At the local scale, the interaction between forest and urban cover affected tawny owl abundance, the optimum being in medium forested areas. The interaction between urban cover and clumpiness index (urban patch distribution) showed a generally negative effect. At the landscape scale, its abundance decreased in complex-shaped urban patches and when distance between them was greater. In conclusion, at the local scale, when a minimal forest structure is present in urbanized areas, the species can exploit it. At the landscape scale, it prefers smaller urban towns to cities. Thinking ahead, the current tendency toward “green capitals” should benefit tawny owl populations.

1. Introduction

The process of urbanization has led to the alteration and destruction of the natural habitat and a decrease in species abundance and richness [1]. As a result [2,3], only a few species are able to successfully adapt to this new urban environment. Although habitat loss is a major consequence of the conversion from native to non-native habitats, habitat fragmentation does have distinct and sometimes influential effects (e.g., [4]). The fragmentation of habitats due to the expansion of cities into their urban surroundings affects the habitats available for birds in all cities across the world [5]. However, the ability to cope with the challenges of living in urban areas depends not only on the characteristics of each species but also on local conditions [6]. Despite a large number of negative influences on urban birds [7], cities still allow ecological niches to occur, meaning that those species able to exploit them perform well [8]. While this leads to an increase in avian biomass, the downside is a reduction in richness [9].
Differences have been observed in bird species along rural–urban gradients, with sensitive species generally tending to locate at the outskirts of the city which still have a relatively high proportion of “natural” habitat. However, with increasing expanses of urbanization within the “core” area of the city, the avian community tends to be dominated by four or five “urban” bird species [10]. Raptors are considered “sentinels” of different local and large-scale environmental change, being sensitive to changes in land use and highly susceptible to local extinctions [11]. For instance, raptor richness tends to be negatively affected by urban development [12]. Indeed, urbanization is leading to systematic raptor population declines due to the reduction in the area of suitable habitat available [13,14]. Nevertheless, some studies have revealed that this group of birds has a certain capacity to colonize urban areas [15,16]. This could be due to the fact that raptors usually have large home ranges, which can extend well beyond the urban boundary and, therefore, do not need to meet all their ecological requirements within urban spaces [9]. In addition, moderate levels of landscape alteration may also permit the co-occurrence of many species, mostly due to edge effects [17]. In such conditions, community structure changes along the transition zone, or ecotone, between different land uses [18]. Some generalist species may in fact even show higher densities and reproduction and survival levels in urban areas [19].
Tawny owl (Strix aluco) is a ubiquitous and resilient raptor that successfully exploits heterogeneous patchy habitats [20]. Due to its abundance and broad geographic range, it is an appropriate bird species model with which to test how raptors colonize urban environments and which drivers from these urban habitats play a role in explaining the occurrence of these species. It is necessary to understand how urbanization affects the tawny owl at the landscape scale [21,22]. The effect of habitat fragmentation (i.e., forest fragmentation) on tawny owl populations has been studied previously [20,21]; however, to our knowledge, the effect of urban habitat fragmentation remains unknown. In addition, it is not common to study habitat selection through a multiscale approach, even though different variables can be important to a species at different spatial scales. As such, studies conducted at arbitrarily defined scales may suffer from serious limitations [23].
In this work, we aimed to study how the composition of urban areas and their distribution and structure affect tawny owl abundance at two different spatial scales (landscape and local). Various studies [14,24] have shown the negative effect of urbanization on raptors. However, the characteristics that urban habitats need to possess for them to be exploited by tawny owls and whether such characteristics vary at different scales have not been examined. It is, however, known that tawny owls prefer heterogeneous habitats; hence, we might predict that variables such as urban patch distribution (clumpiness index) and the number of urban patches in an area may influence species abundance, although the effect might vary at the local and at the landscape scale. Our specific objectives were to investigate the influence of the quantity, fragmentation, and complexity of urban areas on tawny owl abundance at the landscape and the local scale. In this way, if management measures need to be taken in the future, the data from this work will enable the best scale for the intervention to be selected and provide an understanding of how it will affect tawny owls.

2. Materials and Methods

2.1. Study Area

This study was carried out in the Basque Country, in northern Spain (7234 km2, Figure 1). It comprises two principal biogeographic areas, the Cantabric and the Mediterranean regions, along with a transition zone (subCantabric) between them. The Cantabric region (subprovince Cantabro-Atlantic) lies along the coast of the Bay of Biscay (northern half of the Basque Country). It has an Atlantic climate and mild temperatures (annual mean temperature is 14 °C), and the mean precipitation is 1200–2000 mm. The landscape is mountainous and densely populated, with extensive, dense urban areas located in valleys. Forestry plantations (Pinus radiata and Eucalyptus spp.) became widespread in the 1980s, gradually replacing grasslands and traditional agricultural activities, but a few remnants of native forest still exist. The Mediterranean region occupies the southern half of the Basque Country. It is drier and has a more continental climate (mean temperature is 20 °C during summer and 4 °C during winter). Mean precipitation is low around the year, approximately 600 mm [24]. The landscape is less mountainous than in the north, being dominated by arable lands, vineyards, Mediterranean scrub, and holm-oak woods, and it is much less densely populated.

2.2. Data Collection

The data used in this work were obtained from a large-scale owl survey conducted in 2018 (see details in [25,26]). First, we randomly selected a total of 65 Universal Transverse Mercator (UTM) squares measuring 5 × 5 km that represented all the vegetation types of the Basque Country (Figure 1 and Table 1), each of which was considered to be a survey unit (SU) of 25 km2. Second, we selected survey points (SPs) in each SU. These were chosen to represent the main habitat types present in each SU. Each SP was surrounded by a 1 km2 buffer zone with a radius 564.2 m, and there was a minimum distance of 1 km between SPs. Overall, we surveyed 527 SP (Table A1). From January to July 2018, we conducted 2507 surveys in total, in the 527 SPs distributed in 65 SUs.
We surveyed owls for seven consecutive months (one survey per month from January to July). As the survey area was very large and the work was carried out over a long period, 24 observers were involved, each of them responsible for surveying a certain number of SUs. In this way, it was possible to have all the SUs surveyed each month, with 4–8 SPs per SU surveyed each time (Table A1; 6.5 ± SD: 0.7; range: 6–8). In each SP, we applied a three-period survey protocol: 5 min listening for spontaneous calls, 5 min of playing recorded calls, and 5 min listening for calls, making a total of 15 min per survey (for more details, see [25,26]). In the recordings, which were created to our own design based on our previous experience (e.g., [27]), we included a mix of territorial and mating calls of both sexes of particular species in tracks lasting 5 min. Calls were downloaded from xeno-canto (https://xeno-canto.org/, accessed on 10 December 2017). The volume of the recording was loud enough to be heard by an observer at a distance of 300 m but not so loud that it caused distortion in the area close to the speakers. All the surveys were conducted during the first few hours after dusk (3.04 ± SD: 1.73). Since surveying eight SPs required on average 180 min, only 1–2 SUs were surveyed in a single night by each observer. Surveys were conducted on dry nights and were suspended if it rained or if it became windy, because these factors significantly reduce owl detectability [25]. The recordings of tawny owl calls were used in February because it is when the species demonstrates high vocal and territorial activity [28,29]. When a tawny owl responded to the recording in any other month, this was also taken into account in the data and analysis, even though the probability of a response was lower. It made it necessary to work with detectability. As a result of this variability in the likelihood of tawny owl responsiveness, the concept of detectability is important to this study. Detectability, or the probability of detection (p), is defined as the probability of detecting at least one individual at a specific site during a survey, assuming that individuals are present at the site during the sampling period [30]. Detectability is an important source of variation in monitoring programs and can often account for false negatives. For example, nocturnal species usually have low detectability rates, which means that, even though tawny owls are present, they often do not call and, thus, are not detected [25]. It is, therefore, important to employ a hierarchical modeling framework for tawny owl populations in order to estimate abundance that is corrected for imperfect detection.

2.3. Environmental Variables

We selected 13 site-specific covariates (Table 2): (1) REG (categorical), relating to the biogeographic zone (Cantabric, SubCantabric, Mediterranean); (2) ALT (linear), altitude (m a.s.l.) of the SP/SU; (3) ALT2 (quadratic), quadratic altitude (m a.s.l.) of the SU. The 527 circular SP buffer areas and 65 square SU areas were overlaid onto digital vegetation and urban maps (www.geoeuskadi.eus, accessed on 15 March 2019) to obtain the percentage of vegetation types; (4) FOR (linear), forested area which incorporated various plantation types; (5) FOR2 (quadratic), the quadratic function of the percentage of forested area; (6) URB (linear), percentage of urban area. We also obtained various indices of the urban habitat structure in each SP and SU using the ‘sample_lsm’ function of the ‘landscapemetrics’ package in R. These were as follows: (7) CAI (linear), the core area index, which refers to the central area of the urban patch; (8) CAI2 (quadratic), the quadratic effect of CAI for each SU; (9) CLU (linear), the clumpiness index, which refers to the distribution of urban patches, having a value between −1 and 1, and indicates whether the urban area is totally disaggregated (−1), randomly distributed (0), or totally aggregated (1); (10) ENN (linear), the mean Euclidean nearest neighbor, which refers to the distance between the closest two patches having the value of 0 when the nearest neighbor is physically adjoining (aggregated patches) and an unlimited upper limit (isolated patches); (11) SHAPE (linear), the mean shape index of the SP, which is a gradient of perimeter dependent shape complexity and is 0 if all the urban patches are squares, while it increases without limit as the complexity of the patch shape increases; (12) NP (linear), the number of urban patches inside the SP; (13) PAF, the urban perimeter–area fractal dimension for each SU, a value that reflects the complexity of the shape of the SU but is not perimeter dependent and gives a value between 1 and 2, where 1 indicates a shape with a very simple perimeter and 2 indicates a shape with highly convoluted perimeter.
Regarding survey covariates, we used the same ones used in our previously published work surveying the entire owl assemblage [25]. These were (1) BROAD (binomial), whether the individual owl response was before (spontaneous, 0) or during/after the recording (1), (2) EXPER (binomial), level of observer’s experience surveying owls (0 for observers with no previous experience, 1 for those with experience) (see, e.g., [31,32]), (3) DATE (linear), Julian date, used to control for seasonal effects, (4) HOUR (linear), time in the day survey, measured as hours after sunset, and (5) WIND (linear), wind speed (km/h). This variable was obtained from the meteorological station closest to each SP (source: http://www.euskalmet.euskadi.eus, accessed on 15 March 2019).
For the analyses at the local scale (1 km2), we took into account the five survey covariates and 10 of the 13 site-specific covariates. CAI2 and ALT2 were not included because the AIC value increased when adding them into the model. PAF was also dropped because a minimum of 10 patches per SP is needed to measure this index, and most of our SPs did not meet this requirement, although it was taken into account in the analyses at the landscape scale (25 km2). The fieldwork at the landscape scale was not designed to consider detectability; therefore, the covariates for the various SPs were not joined as there could be considerable variability from one to another. Thus, in the landscape-scale models, only the 13 site-specific covariates were considered. The survey covariates were not considered (none were included).

2.4. Statistical Analysis

The analyses were run at the two spatial scales. The predictive variables recorded for the analysis can be classified into two groups: those that might affect species detectability (survey covariates) and those that might affect tawny owl abundance (site covariates).
We performed binomial N-mixture occupancy models where abundance and detectability were estimated as a function of site-specific and survey-specific covariates using the log link function [33]. Observations were generated through a combination of (a) a state process to determine abundance (i.e., counts) at each site, and (b) a detection process that yields observations conditional on the state process [30,34,35]. To control for these two sources of variation, models were built with the “pcount’=” function from the “unmarked” package [33] in R [36]. At the local scale, the abundance of tawny owls per SP was considered as the response variable, with a zero-inflated Poisson distribution. Second, at the landscape scale, the maximum number of tawny owls detected in all the SPs in an SU was considered as the response variable, with a Poisson distribution of errors.
Given the differences between local and landscape scale characteristics, model fitting was conducted in a slightly different way in each case. At the local scale, as the detectability model was previously obtained [24], we kept these variables (DATE, HOUR, WIND, BROAD, and EXPER) constant when performing the abundance model selection, which was performed in a hierarchical process that considered all the possible combinations of the site-specific covariates [36]. We first standardized the variables and fitted both linear and quadratic functions for each variable, and then ran an ANOVA for each pair of models per variable and chose to include the quadratic function when it provided fit better than the linear function. Consequently, we considered the quadratic function of forest (FOR2) for the local scale and the quadratic functions of forest (FOR2), altitude (ALT2), and urban core area index (CAI2) for the landscape scale. In addition, in case tawny owl abundance is affected by the interaction of different variables with the urban variable, we tried different saturated models (at both local and landscape scale) with and without interactions. In the case of local scale, the interactions of FOR, FOR2, CLU, and NP lowered AIC values; hence, these interactions were included in the model. Model selection was carried out using the AIC criterion [37]. Models differing by <2 AIC units were considered to fit the data equally well, and, in this case, we used model averaging procedures to calculate estimates of regression parameters. The relative importance of each predictor (RVI) was assessed by summing the AIC weights across the highest-ranked models that include that predictor [38].

3. Results

Overall, 680 tawny owl calls were noted. Our models showed an evident negative effect of URB on species abundance at both spatial scales (Table 3), although it was strongly influenced by other factors that also impacted species abundance.

3.1. Local Scale

Tawny owls were detected in 371 SPs (70.4%), at an incidence of 1–5 individuals per SP. Models revealed an effect on abundance for five variables (REG, ALT, FOR2, URB, and CLU), as well as the interaction of FOR × URB and CLU × URB (Table 3). These variables were always included as main effects in the top-ranked models (Table A2 and Table A3). There was a positive influence on tawny owl abundance related to Cantabric region, ALT, and FOR (Table 3 and Figure 2). The FOR variable showed a quadratic effect, with the species showing a preference for intermediate percentages of forest cover, rather than highly or poorly forested areas (Figure 2). In contrast, the proportion of urbanized area (URB) showed a negative effect (Figure 2), and CLU revealed tawny owls’ predilection for sites where the urban space is disaggregated (Figure 2). In addition, the FOR × URB interaction showed a negative effect with respect to lowly and highly forested areas (Figure 3) and was optimal in intermediate forested areas, where species abundance was almost constant throughout the urbanization gradient (Figure 3). Models also showed how the effect of urbanization changed depending on its aggregation level (CLU × URB; Figure 4). Urban patch aggregation had almost no effect in lowly urbanized areas but showed a negative effect as urbanization increased (Figure 4).

3.2. Landscape Scale

The maximum number of tawny owls detected in a single SU was 22, and the minimum was 0. The mean of the maximum number of individuals detected was 10.49, and the median was 10. Abundance models revealed clear support for five factors (REG, FOR2, ENN, SHAPE, and NP), two of them already detected in the local-scale analysis (REG and FOR2). No interactions were included at this scale. These factors were always included as the main effect in the top ranked models (Table A2 and Table A3), as well as in the most parsimonious model (Table 3). REG had a very slight effect in favor of the Cantabric region. FOR2 showed a negative effect (Figure 5), reaching its highest values at intermediate levels of forested area. ENN and SHAPE also showed a negative effect on species abundance (Figure 5 and Figure 6). Tawny owl abundance decreased when the shape of urban patches was more complex and when the distance between urban patches increased. In contrast, NP had a slight positive effect (Figure 6), with species abundance increasing as the number of urban patches increased.

4. Discussion

In a previous study, we analyzed the effects of survey related factors (i.e., the observer’s experience and wind speed) and the principal site-related habitat factors on tawny owl detection and abundance in the study area [25]. We demonstrated that the presence of forested habitats positively affected its abundance, whereas urban habitats had a negative effect. In this study, which complements the previous one, we focused on analyzing the effects the distribution and structure of urban habitats have on the species abundance. Our results show that the abundance of tawny owls in urban habitats strongly depends on the forested areas that surround them, as well as on the structure of the urban area itself. As such, the best sites for tawny owls in cities are periurban areas or areas where an urban–forested matrix is present. Therefore, our results also show new biodiversity potential associated with urban areas as an urban ecosystem service [39] through increasing numbers of tawny owls.
When a local-scale approach was used, the species was found to be more abundant in places with a higher proportion of forests, such as the Cantabric region. This confirms results from our previous study which analyzed the main distribution pattern of this owl in the region [26]. One novel finding, however, is that we detected a quadratic effect, indicating that both extensive forest areas and small patches of forest negatively affected tawny owl abundance (for similar results, see also [40]), this being different to results obtained in another work [41]. It seems, therefore, that the species may benefit from using woodland areas with some degree of forest fragmentation, with the non-forested areas comprising small patches of open habitats such as grasslands, meadows, or pastures [42]. Urbanization showed a forest-dependent effect in addition to its overall negative effect. When the forested area was low (<30%), the percentage of urban area had a slightly negative effect on tawny owl abundance. At intermediate forest percentages (30–60%), urbanization had almost no effect, with abundance being more or less constant throughout the urban range. However, when forest covered extensive areas (>60%), the percentage of urbanized area showed the worst effect. This is probably due to the permanent food supplies that urban environments provide [43], as well as the availability of suitable new nesting and roosting sites [44]. In addition, it has been demonstrated that home range size decreases at the intermediate level of forested area [42], with higher densities being found in urban areas mixed with patchily distributed forests than in purely forested areas [20,45], but always depending on the proportion of suitable habitats available [19]. Regarding the aggregation level of urban patches (CLU), something similar happens in areas of low urbanization (<10%), where aggregation level has a slight, nearly constant, effect on tawny owl abundance. Nevertheless, the negative effect becomes stronger as urbanization level increases, probably due to the lack of suitable habitats since tawny owl is a forest-dwelling raptor and, therefore, needs trees to be able to exploit urban areas. Moreover, in highly urbanized areas (>20%), where urban patches tend to be more aggregated, there are not so many trophic resources to exploit due to the reduced availability of green areas. Tawny owls have been shown to prefer heterogeneous patchy habitats [46], although urban habitat was not taken into account in this study. Our results suggest that this habitat could in fact be included as yet another habitat type that the species can exploit. At this point, we should note that the degree of urbanization within our region ranged from moderate to low since our SUs and SPs were never surrounded by more than 20% of urban land use. This is a good reflection of our study area, which mainly comprises spread out rural villages with few cities.
At the landscape scale, our results also showed that the highest abundance was at intermediate levels of forest cover, again contrasting with previous results [41]. Nevertheless, the amount of urban habitat did not appear to affect abundance, which could be due to the habitat variability encompassed at this scale. However, two urban indices showed significant effects, which means that urbanization affects tawny owl abundance differently when considered at different scales, showing that multiple scales need to be taken into account in future studies [23]. The distance to the nearest neighboring patch of the same class (ENN) resulted in a quadratic effect. That is, when urban patches were very close to each other, the space for forest or crop patches was not big enough or did not contain sufficient diversity to ensure high tawny owl abundances. In contrast, when urban patches were very far apart, the forest or crop habitats between them were larger. As a consequence of tawny owl preference for an urban–rural mix, with high levels of forest or cropland, the abundance of the species will decrease. It is, thus, clear that there is an optimal point where the gap between urban areas is sufficient to maintain a mixture of different crops and forested areas capable of ensuring a high trophic resource, which likely results in an increase in tawny owl abundance. Regarding the urban patch shape index (SHAPE), our results do not match with other researchers [47,48], who observed that the shape index was not relevant to species abundance. In the present study, however, an increase in shape complexity, which is related to a lower urban core area, showed a negative effect on species abundance. This could be attributable to the low urbanization rate in the study area, as well as the narrow range of the variable. The shape index starts at 0 and can increase indefinitely, but our shape index range was only from 1.2 to 1.7. As our shape complexity ratio was so small, we probably did not observe the real effect of this variable. Thus, a wider range of shape complexity should be considered in future studies. The number of urban patches showed a positive effect on species abundance here, which may have been a consequence of the previously mentioned suitable habitats [19] and resource availability, as towns and villages offer high levels of trophic resources and nesting niches.

5. Conclusions

Tawny owls have the capacity to colonize new environments [49,50] depending on the surrounding habitat and urban habitat structure [18,20]. This ability is known to be associated with success in the colonization of urban habitats, as well as with ecological and evolutionary consequences in terms of consistent changes in distribution, abundance, behavior, and life history [51]. At the local scale, when minimum forest structure is available in urbanized areas or its surroundings, and when urbanization is patchier, the species can easily exploit it, which has a positive effect on its abundance [18,20]. At the landscape scale, we observed that this species prefers smaller urban villages split in more patches than a large, urbanized city. In this sense, we think that the current tendency toward “green capitals” (the Green Capital Award is a European award for cities based on their environmental records) and the environmental changes involved [52,53] may benefit tawny owl populations, at least in terms of abundance in urban areas. From a biodiversity [54] and pollution [55] point of view, it is valuable to have an urban bioindicator such as an urban tawny owl population, as well as the possibility of monitoring it over time. Nevertheless, we observed that tawny owls do not exploit highly urbanized environments. Our study contributes to the building of friendly cities through adding to our knowledge of how new urban areas should be constructed. It is recommended to maintain small forest patches, transforming the urban matrix into heterogeneous areas, where tawny owls can exploit them as foraging areas in this new urban ecosystem [20]. In the context of global change where landscapes are increasingly anthropized, we want to highlight the importance of assessing multiscale relationships between urban biodiversity and the generation of new urban ecosystem services [39].

Author Contributions

Conceptualization, I.Z. and J.A.; methodology, I.Z.; validation, I.Z., M.V.J.-F. and J.A.; formal analysis, N.P.; data curation, N.P.; writing—original draft preparation, N.P., I.Z., M.V.J.-F. and J.A.; writing—review and editing, N.P., I.Z., M.V.J.-F. and J.A.; visualization, M.V.J.-F.; supervision, I.Z. and J.A.; funding acquisition, J.A. All authors have read and agreed to the published version of the manuscript.

Funding

The fieldwork for this research was funded by the Department of the Environment, Territorial Planning, and Housing of the Basque Government. N.P. was supported by a Basque Government predoctoral fellowship grant (PRE_2018_1_0004). M.V.J.-F. was supported by a postdoctoral grant “Juan de la Cierva-Incorporación” (reference IJC2019-039145-I) financed by the MCIN/AEI/ 10.13039/501100011033).

Institutional Review Board Statement

Not applicable. Not humans or animals were involved. The study was conducted with permission of the three Basque Councils and the Basque Government for using playbacks during the fieldwork.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://doi.org/10.5061/dryad.dncjsxkwg.

Acknowledgments

The fieldwork was funded by Department of the Environment, Territorial Planning, and Housing, Basque Government. Fieldwork was carried out by I.Z., N.P., Izaskun Aguirre, Ainara Azkona, Gorka Belamendia, Sergio de Juan, Rosa de Siria Apaolaza, Alexandra Egunez, Nober Fuente, Jordi Gómez, Carlos González de Buitrago, Ana Gracianteparaluceta, Maite Laso, Emilio López, Alfonso López de Armentia, Javier López de Luzuriaga, Diego Malo, Jabi Manzano, I. Martínez, Aitor Montes, Jose Ignacio Ocariz, Silvia Prats, Rakel Santos, and Brian Webster. Ronnie Lendrum helped by improving the written English.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table
Table A1. A summary table of the number of survey points (SPs) in each survey unit (SU) each month. In total, 2507 SPs were surveyed.
Table A1. A summary table of the number of survey points (SPs) in each survey unit (SU) each month. In total, 2507 SPs were surveyed.
(Month) JanuaryFeburaryMarchAprilMayJune–JulyMean
(Playback Sp.) Strix aluco Values
SUTotal SPs per SUSPs SurveyedSPs SurveyedSPs SurveyedSPs SurveyedSPs SurveyedSPs SurveyedSPs Surveyed
30TVN08NW11010000
30TVN68SE85666866
30TVN78NW84656776
30TVN84NW88866887
30TVN88NW85876666
30TVN93NW88866877
30TVN96NE88866667
30TVN96NW88786667
30TVN98NW84666776
30TVN98SW84865666
30TVN99NW87468676
30TWN02NE88876667
30TWN03NE86886867
30TWN04NW88866667
30TWN07NW88766656
30TWN07SW88876647
30TWN09SW95557666
30TWN13SE88868888
30TWN14NE88888888
30TWN16NE88865556
30TWN16NW85745866
30TWN18NW85666766
30TWN18SE85667666
30TWN19NW88877777
30TWN21SE11100000
30TWN24NE87536666
30TWN24SW86626676
30TWN25NW1210986467
30TWN26SE77777777
30TWN30NE88866687
30TWN30SW 88876767
30TWN31NE33300001
30TWN31SW 88586667
30TWN33SW87656566
30TWN34NE86847846
30TWN36SE86666666
30TWN37SW86677877
30TWN38SW87776867
30TWN39NE87788777
30TWN42SW88866667
30TWN43NW88456666
30TWN46SE86676746
30TWN48NW87778888
30TWN52SW88866667
30TWN54NW66666666
30TWN55NE86676656
30TWN56NW86655686
30TWN57NE84666686
30TWN58SE84676676
30TWN59SW86456686
30TWN67SW86566766
30TWN69SW86556676
30TWN75NW86664776
30TWN76NE85764776
30TWN76NW87564776
30TWN76SW85765776
30TWN77NE84657776
30TWN78NE84677756
30TWN78SW84686666
30TWN89NE86566666
30TWN89SW85766676
30TWN99NW86554786
30TWN99SW86466776
30TWP00NE87788878
30TWP10SW87788878
30TWP20SE88778888
30TWP90SE86456886
30TWP90SW86464786
Total527422426402399434424
Table
Table A2. Local-scale model selection. In the right-hand columns, the AIC, deltaic, and AIC weights (wAIC; used for calculating the RVI index) for each model are shown. The variables used for detectability were as follows: P (DATE + HOUR + WIND + BROAD + EXPER).
Table A2. Local-scale model selection. In the right-hand columns, the AIC, deltaic, and AIC weights (wAIC; used for calculating the RVI index) for each model are shown. The variables used for detectability were as follows: P (DATE + HOUR + WIND + BROAD + EXPER).
ModelsFormulaAICdeltaAICwAIC
m.25REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + FOR:URB + FOR2:URB + CLM:URB3235.880.000.28
m.23REG + ALT + FOR + FOR2 + URB + CAI + CLU + SHAPE + NP + FOR:URB + FOR2:URB + CLM:URB + NP:URB3236.911.030.17
m.29REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + NP + FOR:URB + FOR2:URB + CLM:URB3239.941.060.16
m.21REG + ALT + FOR + FOR2 + URB + CLU + ENN + SHAPE + NP + FOR:URB + FOR2:URB + CLM:URB + NP:URB3237.601.730.12
m.24REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + NP + FOR:URB + FOR2:URB + CLM:URB + NP:URB3238.092.210.092
m.15 (sat)REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + NP + FOR:URB + FOR2:URB + CLM:URB + NP:URB3238.812.930.064
m.22REG + ALT + FOR + FOR2 + URB + CAI + ENN + SHAPE + NP + FOR:URB + FOR2:URB + NP:URB3239.663.780.042
m.27REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + NP + FOR:URB + CLM:URB + NP:URB3240.104.220.034
m.26REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + NP + CLM:URB + NP:URB3241.185.300.020
m.28REG + ALT + FOR + FOR2 + URB + CAI + CLU + ENN + SHAPE + NP + FOR:URB + FOR2:URB + NP:URB3241.415.530.018
m.19REG + ALT + FOR + URB + CAI + CLU + ENN + SHAPE + NP + FOR:URB + CLM:URB + NP:URB3245.029.150.003
Table
Table A3. Landscape-scale model selection. In the columns on the right, the AIC, deltaAIC, and AIC weights (wAIC; used for calculating the RVI index) for each model are shown. No detectability variables were used at this scale: P(.).
Table A3. Landscape-scale model selection. In the columns on the right, the AIC, deltaAIC, and AIC weights (wAIC; used for calculating the RVI index) for each model are shown. No detectability variables were used at this scale: P(.).
ModelsFormulaAICdeltaAICwAIC
mo.11REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + ENN + PAF + SHAPE + NP1234.390.000.22
mo.8REG + ALT + ALT2 + FOR + FOR2 + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1234.740.350.18
mo.13REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + SHAPE + NP1235.551.150.12
mo.4REG + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1235.811.410.11
mo.2 (sat)REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1236.351.960.081
mo.9REG + ALT + ALT2 + FOR + FOR2 + URB + CLU + ENN + PAF + SHAPE + NP1236.361.970.081
mo.3ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1236.782.390.065
mo.5REG + ALT + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1236.812.410.065
mo.10REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CLU + ENN + PAF + SHAPE + NP1237.733.330.041
mo.15REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE1237.873.480.038
mo.12REG + ALT + ALT2 + FOR + FOR2 + URB + CAI + CAI2 + CLU + PAF + SHAPE + NP1243.048.650.0029
mo.7REG + ALT + ALT2 + FOR + URB + CAI + CAI2 + CLU + ENN + PAF + SHAPE + NP1243.999.600.0018

References

  1. McKinney, M.L. Urbanization, biodiversity, and conservation. BioScience 2002, 52, 883–890. [Google Scholar] [CrossRef]
  2. Sol, D.; González-Lagos, C.; Moreira, D.; Maspons, J.; Lapiedra, O. Urbanization tolerance and the loss of avian diversity. Ecol. Lett. 2014, 17, 942–950. [Google Scholar] [CrossRef]
  3. Lopucki, R.; Klich, D.; Kitowski, I.; Kiersztyn, A. Urban size effect on biodiversity: The need for a conceptual framework for the implementation of urban policy for small cities. Cities 2020, 98, 102590. [Google Scholar] [CrossRef]
  4. De Lima Filho, J.A.; Vieira, R.J.A.; de Souza, C.A.M.; Ferreira, F.F.; de Oliveira, V.M. Effects of habitat fragmentation on biodiversity patterns of ecosystems with resource competition. Physica A 2021, 564, 125497. [Google Scholar] [CrossRef]
  5. Hedblom, M.; Murgui, E. Urban Bird Research in a Global Perspective. In Ecology and Conservation of Birds in Urban Environments, 1st ed.; Murgui, E., Hedblom, M., Eds.; Springer: Cham, Switzerland, 2017; pp. 3–10. [Google Scholar] [CrossRef]
  6. Shanahan, D.F.; Stronhbach, M.W.; Warren, P.S.; Fuller, R.A. The challenges of urban living. In Avian Urban Ecology, 1st ed.; Gil, D., Brumm, H., Eds.; OUP Oxford: Oxford, UK, 2014; pp. 3–20. [Google Scholar] [CrossRef]
  7. Santos, S.M.; Lourenco, R.; Mira, A.; Beja, P. Relative Effects of Road Risk, Habitat Suitability, and Connectivity on Wildlife Roadkills: The Case of Tawny Owls (Strix aluco). PLoS ONE 2013, 8, e79967. [Google Scholar] [CrossRef] [Green Version]
  8. Donázar, J.A.; Cortés-Avizanda, A.; Fargallo, J.A.; Margalida, A.; Moleón, M.; Morales-Reyes, Z.; Moreno-Opo, R.; Pérez-García, J.M.; Sánchez-Zapata, J.A.; Zuberogoitia, I.; et al. Roles of Raptors in a Changing World: From Flagships to Providers of Key Ecosystem Services. BioOne 2016, 63, 181–234. [Google Scholar] [CrossRef] [Green Version]
  9. Chace, J.F.; Walsh, J.J. Urban effects on native avifauna: A review. Landsc. Urban Plan. 2006, 74, 46–69. [Google Scholar] [CrossRef]
  10. Palomino, D.; Carrascal, L.M. Urban influence on birds at a regional scale: A case study with the avifauna of northern Madrid province. Landsc. Urban Plan. 2006, 77, 276–290. [Google Scholar] [CrossRef]
  11. Rullman, S.; Marzluff, J.M. Raptor presence along an urban-wildland gradient: Influences of prey abundance and land cover. J. Raptor Res. 2014, 48, 257–272. [Google Scholar] [CrossRef]
  12. Palomino, D.; Carrascal, L.M. Habitat associations of a raptor community in a mosaic landscape of Central Spain under urban development. Landsc. Urban Plan. 2007, 83, 268–274. [Google Scholar] [CrossRef]
  13. Martínez, J.A.; Zuberogoitia, I. Habitat references and causes of population decline for barn owls Tyto alba: A multiscale approach. Ardeola 2004, 51, 303–317. [Google Scholar]
  14. Ivajnšič, D.; Denac, D.; Denac, K.; Pipenbaher, N.; Kaligarič, M. The Scops owl (Otus scops) under human-induces environmental change pressure. Land Use Policy 2020, 99, 104853. [Google Scholar] [CrossRef]
  15. Solonen, T. Timing of breeding in rural and urban Tawny Owls Strix aluco in southern Finland: Effects of vole abundance and winter weather. J. Ornithol. 2014, 155, 27–36. [Google Scholar] [CrossRef]
  16. Solonen, T.; Lokki, H.; Sulkava, S. Diet and brood size in rural and urban Northern Goshawks Accipiter gentilis in southern Finland. Avian Biol. Res. 2019, 12, 3–9. [Google Scholar] [CrossRef]
  17. Anderson, D.L. Landscape Heterogeneity and Diurnal Raptor Diversity in Honduras: The Role of Indigenous Shifting Cultivation. Biotropica 2001, 33, 511–519. [Google Scholar] [CrossRef]
  18. Fröhlich, A.; Ciach, M. Nocturnal noise and habitat homogeneity limit species richness of owls in an urban environment. Environ. Sci. Pollut. Res. 2019, 26, 17284–17291. [Google Scholar] [CrossRef] [Green Version]
  19. Ranazzi, L.; Manganaro, A.; Salvati, L. Density fluctuations in urban population of Tawny Owl Strix aluco: A long-term study in Rome, Italy. Ornis Svecica 2002, 12, 63–67. [Google Scholar]
  20. Fröhlich, A.; Ciach, M. Noise pollution and decreased size of wooded areas reduces the probability of occurrence of Tawny Owl Strix aluco. Ibis 2018, 160, 634–646. [Google Scholar] [CrossRef]
  21. Redpath, S.M. Habitat Fragmentation and the Individual: Tawny Owl Strix aluco in Woodland Patches. J. Anim. Ecol. 1995, 64, 652–661. [Google Scholar] [CrossRef]
  22. Sunde, P. What do we know about territorial behaviour and its consequences in Tawny Owls? In Ecology and Conservation of European Forest-Dwelling Raptors, 1st ed.; Zuberogoitia, I., Martínez, J.E., Eds.; Departamento de Agricultura de la Diputación Foral de Bizkaia: Bilbao, Spain, 2011; pp. 253–260. [Google Scholar]
  23. Martínez, J.A.; Serrano, D.; Zuberogoitia, I. Predictive models of hábitat preferences for the Eurasian Eagle Owl Bubo bubo: A multiscale approach. Ecography 2003, 26, 21–28. [Google Scholar] [CrossRef] [Green Version]
  24. Euskalmet: Agencia Vasca de Meteorología. El Clima en Euskadi: Zonas Climáticas; Gobierno Vasco: Vitoria-Gasteiz, Spain. Available online: https://www.euskalmet.euskadi.eus/clima/euskadi/ (accessed on 15 March 2019).
  25. Zuberogoitia, I.; Martínez, J.E.; González-Oreja, J.A.; González de Buitrago, C.; Belamendia, G.; Zabala, J.; Laso, M.; Pagaldai, N.; Jiménez-Franco, M.V. Maximizing detection probability for effective large-scale nocturnal bird monitoring. Divers. Distrib. 2020, 26, 1034–1050. [Google Scholar] [CrossRef]
  26. Zuberogoitia, I.; Martínez, J.E.; González-Oreja, J.A.; González de Buitrago, C.; Belamendia, G.; Zabala, J.; Laso, M.; Pagaldai, N.; Jiménez-Franco, M.V. Data from: Maximizing detection probability for effective large-scale nocturnal bird monitoring. Divers. Distrib. 2020. Dryad, Dataset. [Google Scholar] [CrossRef]
  27. Zuberogoitia, I.; Campos, L.F. Censusing owls in large areas: A comparison between methods. Ardeola 1998, 45, 47–53. [Google Scholar]
  28. Zuberogoitia, I.; Martínez, J.A. Methods for surveying Tawny Owl Strix aluco populations in large areas. Biota 2000, 1, 137–146. [Google Scholar]
  29. Zuberogoitia, I.; Burgos, G.; González-Oreja, J.A.; Morant, J.; Martínez, J.E.; Zabala, J. Factors affecting spontaneous vocal activity of Tawny Owls Strix aluco and implications for surveying large areas. Ibis 2019, 161, 495–503. [Google Scholar] [CrossRef]
  30. MacKenzie, D.I.; Nichols, J.D.; Lachman, G.B.; Droege, S.; Royle, J.A.; Langtimm, C.A. Estimating site occupancy rates when detection probabilities are less than one. Ecology 2002, 83, 2248–2255. [Google Scholar] [CrossRef]
  31. Zabala, J.; Zuberogoitia, I.; Martínez-Climents, J.A.; Martínez, J.E.; Azkona, A.; Hidalgo, S.; Iraeta, A. Occupancy and abundance of Little Owl Athene noctua in an intensively managed forest area in Biscay. Ornis Fennica 2006, 83, 97–107. [Google Scholar]
  32. Zuberogoitia, I.; Zabala, J.; Martínez, J.E. Bias in little owl population estimates using playback techniques during surveys. Anim. Biodivers. Conserv. 2011, 34, 395–400. [Google Scholar]
  33. Fiske, I.; Chandler, R. Unmarked: An R package for fitting hierarchical models of wildlife occurrence and abundance. J. Stat. Softw. 2011, 43, 1–23. [Google Scholar] [CrossRef] [Green Version]
  34. Kéry, M.; Royle, J.A.; Schmid, H. Modeling avian abundance from replicated count using binomial mixture models. Ecol. Appl. 2005, 15, 1450–1461. [Google Scholar] [CrossRef]
  35. Royle, J.A. N-Mixture Models for Estimating Population Size from Spatially Replicated Counts. Biometrics 2004, 60, 108–115. [Google Scholar] [CrossRef]
  36. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 23 July 2021).
  37. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  38. Giam, X.; Olden, J.D. Quantifying variable importance in a multimodel inference framework. Methods Ecol. Evol. 2016, 7, 388–397. [Google Scholar] [CrossRef]
  39. Krener, P.; Hamstead, Z.; Haase, D.; McPhearson, T.; Frantzeskaki, N.; Andersson, E.; Kabisch, N.; Larondelle, N.; Rall, E.L.; Voigt, A.; et al. Key insights for the future of urban ecosystem services research. Ecol. Soc. 2016, 21, 29. [Google Scholar] [CrossRef] [Green Version]
  40. Michel, V.T.; Jiménez-Franco, M.V.; Naef-Daenzer, B.; Grüebler, M.U. Intraguild predator drives forest edge avoidance of a mesopredator. Ecosphere 2016, 7, e01229. [Google Scholar] [CrossRef] [Green Version]
  41. Ranazzi, L.; Manganaro, A.; Ranazzi, R.; Salvati, L. Woodland cover and Tawny Owl Strix aluco density in a Mediterranean urban area. Biota 2000, 1, 27–34. [Google Scholar]
  42. Burgos, G.; Zuberogoitia, I. A telemetry study to discriminate between home range and territory size in Tawny Owls. Bioacustics 2018, 29, 109–121. [Google Scholar] [CrossRef]
  43. Rodewald, A.D.; Shustack, D.P. Consumer resource matching in urbanizing landscapes: Are synanthropic species over-matching? Ecology 2008, 89, 515–521. [Google Scholar] [CrossRef] [Green Version]
  44. Sacchi, R.; Gentilli, A.; Razzetti, E.; Barbieri, F. Effects of building features on density and flock distribution of Feral Pigeons Columba livia var. domestica in an urban environment. Can. J. Zool. 2002, 80, 48–54. [Google Scholar] [CrossRef] [Green Version]
  45. Millsap, B.A. Survival of Florida Burrowing Owls along an urban-development gradient. J. Raptor Res. 2002, 36, 3–10. [Google Scholar]
  46. López-Peinado, A.; Lis, A.; Perona, A.M.; López-López, P. Habitat preferences of the tawny owl (Strix aluco) in a special conservancy area of eastern Spain. J. Raptor Res. 2020, 54, 402–413. [Google Scholar] [CrossRef]
  47. Whytock, R.C.; Fuentes-Montemayor, E.; Watts, K.; Macgregor, N.A.; Call, E.; Mann, J.A.; Park, K.J. Regional land-use and local management create scale-dependent ‘landscapes of fear’ for a common woodland bird. Landsc. Ecol. 2020, 35, 607–620. [Google Scholar] [CrossRef] [Green Version]
  48. Alaniz, A.J.; Carvajal, M.A.; Fierro, A.; Vergara-Rodríguez, V.; Toledo, G.; Ansaldo, D.; Moreira-Arce, D.; Rojas-Osorio, A.; Vergara, P.M. Remote-sensing estimates of forest structure and dynamics as indicators of habitat quality for Magellanic woodpeckers. Ecol. Indic. 2021, 126, 107634. [Google Scholar] [CrossRef]
  49. Mak, B.; Francis, R.A.; Chadwick, M.A. Living in the concrete jungle: A review and socio-ecological perspective of urban raptor habitat quality in Europe. Urban Ecosyst. 2021, 1–21. [Google Scholar] [CrossRef]
  50. Gryz, J.; Krauze-Gryz, D. Influence of habitat urbanization on time of breeding and productivity of tawny owl (Strix aluco). Pol. J. Ecol. 2018, 66, 153–161. [Google Scholar] [CrossRef]
  51. Møller, A.P.; Díaz, M. Avian preference for close proximity to human habitation and its ecological consequences. Curr. Zool. 2018, 64, 623–630. [Google Scholar] [CrossRef]
  52. European Commision. European Commision, Europe 2020: A Trategy for Smart, Sustainable and Inclusive Growth. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:52010DC2020&from=EN (accessed on 23 July 2021).
  53. Gulsrud, N.M.; Ostoić, S.K.; Faehnle, M.; Maric, B.; Paloniemi, R.; Pearlmutter, D.; Simson, A.J. Challenges to Governing Urban Green Infrastructure in Europe—The Case of the European Green Capital Award. In The Urban Forest, 1st ed.; Pearlmutter, D., Calfapietra, C., Samson, R., O’Brien, L., Ostoić, S.K., Sanesi, G., Alonso del Amo, R., Eds.; Springer: Cham, Switzerland, 2017; Volume 7, pp. 235–258. [Google Scholar] [CrossRef]
  54. Sergio, F.; Newton, I.; Marchesi, L. Top predators and biodiversity. Nature 2005, 436, 192. [Google Scholar] [CrossRef]
  55. Lodenius, M.; Solonen, T. The use of birds of prey as indicators of metal pollution. Ecotoxicology 2013, 22, 1319–1334. [Google Scholar] [CrossRef] [PubMed]
Animals 11 02954 g001 550
Figure 1. Distribution of the 527 survey points (SPs) and 65 survey units (SUs) in the study area. Forest and urban habitats are shown, which were the land cover types present in the top-ranked model.
Figure 1. Distribution of the 527 survey points (SPs) and 65 survey units (SUs) in the study area. Forest and urban habitats are shown, which were the land cover types present in the top-ranked model.
Animals 11 02954 g001
Animals 11 02954 g002 550
Figure 2. Expected abundance of tawny owl (Strix aluco) at the local scale and based on the top-ranked model in relation to (a) forested gradient, (b) urban gradient, and (c) clumpiness index. The clumpiness index represents the aggregation level of habitat patches (in our case, urban patches). It is equal to −1 when the urban patch is maximally disaggregated, 0 when it is randomly distributed, and 1 when maximally aggregated. The blue lines show a cubic smoothing spline fitted to the data in order to show the general trend, while shaded areas show standard error boundaries.
Figure 2. Expected abundance of tawny owl (Strix aluco) at the local scale and based on the top-ranked model in relation to (a) forested gradient, (b) urban gradient, and (c) clumpiness index. The clumpiness index represents the aggregation level of habitat patches (in our case, urban patches). It is equal to −1 when the urban patch is maximally disaggregated, 0 when it is randomly distributed, and 1 when maximally aggregated. The blue lines show a cubic smoothing spline fitted to the data in order to show the general trend, while shaded areas show standard error boundaries.
Animals 11 02954 g002
Animals 11 02954 g003 550
Figure 3. Expected tawny owl (Strix aluco) abundance at the local scale and based on the top-ranked model in three forest levels (a–c) along an urban gradient. The blue lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Figure 3. Expected tawny owl (Strix aluco) abundance at the local scale and based on the top-ranked model in three forest levels (a–c) along an urban gradient. The blue lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Animals 11 02954 g003
Animals 11 02954 g004 550
Figure 4. Expected tawny owl (Strix aluco) abundance at the local scale and based on the top-ranked model in three levels of urbanization (a–c) and at different aggregation levels of the urban patches (clumpiness index). The clumpiness index represents the aggregation level of habitat patches (in our case, urban patches). It is equals to −1 when the urban patch is maximally disaggregated, 0 when it is randomly distributed, and 1 when maximally aggregated. The blue lines show a cubic smoothing spline fitted to the data in order to show the general trend, while shaded areas show standard error boundaries.
Figure 4. Expected tawny owl (Strix aluco) abundance at the local scale and based on the top-ranked model in three levels of urbanization (a–c) and at different aggregation levels of the urban patches (clumpiness index). The clumpiness index represents the aggregation level of habitat patches (in our case, urban patches). It is equals to −1 when the urban patch is maximally disaggregated, 0 when it is randomly distributed, and 1 when maximally aggregated. The blue lines show a cubic smoothing spline fitted to the data in order to show the general trend, while shaded areas show standard error boundaries.
Animals 11 02954 g004
Animals 11 02954 g005 550
Figure 5. Expected tawny owl (Strix aluco) abundance at the landscape scale, based on the top-ranked model along (a) a forested gradient, and (b) different Euclidean nearest neighbor distances (ENNs) of urban patches. The orange lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Figure 5. Expected tawny owl (Strix aluco) abundance at the landscape scale, based on the top-ranked model along (a) a forested gradient, and (b) different Euclidean nearest neighbor distances (ENNs) of urban patches. The orange lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Animals 11 02954 g005
Animals 11 02954 g006 550
Figure 6. Expected tawny owl (Strix aluco) abundance at the landscape scale, based on the top-ranked model with respect to (a) different shape index of urban patches, and (b) the number of urban patches. The shape index represents the degree of irregularity of habitat patches (in our case, urban habitat). It is equal to 1 when the urban patch is maximally compact (almost square) and increases without limit as the patch shape becomes more irregular. The orange lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Figure 6. Expected tawny owl (Strix aluco) abundance at the landscape scale, based on the top-ranked model with respect to (a) different shape index of urban patches, and (b) the number of urban patches. The shape index represents the degree of irregularity of habitat patches (in our case, urban habitat). It is equal to 1 when the urban patch is maximally compact (almost square) and increases without limit as the patch shape becomes more irregular. The orange lines show a cubic smoothing spline fitted to the data to show the general trend, while shaded areas show standard error boundaries.
Animals 11 02954 g006
Table
Table 1. Comparison between the percentage distribution of land types by survey unit (SU, 25 km2), by survey point (SP, 1 km2), and across the whole study area. Values were obtained by calculating the mean habitat (forest and urban) percentage and standard deviation for each SU and SP area, as well as for the whole study area.
Table 1. Comparison between the percentage distribution of land types by survey unit (SU, 25 km2), by survey point (SP, 1 km2), and across the whole study area. Values were obtained by calculating the mean habitat (forest and urban) percentage and standard deviation for each SU and SP area, as well as for the whole study area.
HabitatSP (Mean % ± SD)SU (Mean % ± SD)Study Area (%)
Forest42.8 ± 26.345.1 ± 23.154.6
Urban6 ± 11.78.2 ± 11.73.1
Table
Table 2. Description of the variables used for the top-ranked models at both spatial scales, local and landscape. The “SCALE” column shows the spatial scale used for each variable in the models (LO: local and LA: landscape).
Table 2. Description of the variables used for the top-ranked models at both spatial scales, local and landscape. The “SCALE” column shows the spatial scale used for each variable in the models (LO: local and LA: landscape).
VariablesContractionDescriptionTypeRangeUnitScale
RegionREGClimatic region (Cantabric, Subcantabric, Mediterranean)Binary0–1–2NoneLO + LA
AltitudeALTAltitude above sea levelContinuous0–1500MetersLO + LA
ForestFORForested area percentageContinuous0–100%LO + LA
UrbanURBUrbanized area percentageContinuous0–100%LO + LA
Mean core area indexCAICAI is the percentage that the core area (interior area) takes in a patch. The mean value of the CAI of all the urban patches in each SP or SU is calculated.Continuous0–100%LO + LA
Clumpiness indexCLUDescribes how the entire group of urban patches is distributed. Equal to −1 for maximally disaggregated, 0 for randomly distributed, and 1 for maximally aggregated urban patches.Continuous−1–1NoneLO + LA
Euclidean nearest mean neighborENNENN measures the distance between one urban patch and its nearest urban patch (does not take into account the whole patch group). It is calculated as the mean of ENN of all the urban patches in each SP or SU.Continuous>0MetersLO + LA
Shape indexSHAPESHAPE describes the ratio between the actual perimeter of the urban patch and its hypothetical minimum perimeter. The mean SHAPE value of all urban patches is calculated.Continuous≥1NoneLO + LA
Number of patchesNPThe number of urban patchesDiscrete≥1 LO + LA
Perimeter Area Fractal DimensionPAFDescribes the complexity of an urban patch. Approaches 1 for those with simple shapes (i.e., like a square) and approaches 2 for those that are very irregular.Continuous1 ≤ PAF ≤ 2 LA
Table
Table 3. Conditional model averaging of the top-ranked models (delta AICc < 2), examining the effect of predicted factors on the abundance of tawny owls. Model average parameter estimates, adjusted standard errors, 95% confidence intervals, and relative variable importance (RVI, the sum of Akaike weights across the set of models in which the variable appears) are shown. Effects with 95% CI not overlapping zero are shown in bold.
Table 3. Conditional model averaging of the top-ranked models (delta AICc < 2), examining the effect of predicted factors on the abundance of tawny owls. Model average parameter estimates, adjusted standard errors, 95% confidence intervals, and relative variable importance (RVI, the sum of Akaike weights across the set of models in which the variable appears) are shown. Effects with 95% CI not overlapping zero are shown in bold.
Parameters NBEstimateSE (Estimate)Lower 95% CIUpper 95% CIRVI
Local scale: 1 km2
Intercept4.770.443.915.63
REG−0.760.11−0.98−0.540.99
ALT0.450.080.280.610.99
FOR0.070.08−0.080.230.99
FOR2−0.250.08−0.40−0.090.99
URB1.010.390.241.770.99
CAI0.070.07−0.060.200.89
CLU−0.510.26−1.010.000.99
ENN0.020.05−0.080.120.83
SHAPE−0.080.06−0.210.050.92
NP−0.080.09−0.250.100.68
FOR × URB−0.060.16−0.370.250.99
FOR2 × URB−0.220.11−0.44−0.010.96
URB × CLU−0.920.48−1.850.010.97
URB × NP0.030.06−0.080.140.51
Landscape scale: 25 km2
Intercept3.81.061.725.87
REG−0.140.14−0.420.130.89
ALT0.060.11−0.160.280.90
ALT2−0.130.07−0.260.010.86
FOR0.220.080.060.371
FOR2−0.210.08−0.36−0.061
URB0.070.1−0.120.260.84
CAI0.310.22−0.120.750.86
CAI2−0.190.13−0.430.060.78
CLU00.09−0.180.180.80
ENN−0.60.22−1.03−0.171
PAF0.080.09−0.110.261
SHAPE−0.530.14−0.8−0.251
NP−0.130.08−0.290.030.96
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pagaldai, N.; Arizaga, J.; Jiménez-Franco, M.V.; Zuberogoitia, I. Colonization of Urban Habitats: Tawny Owl Abundance Is Conditioned by Urbanization Structure. Animals 2021, 11, 2954. https://doi.org/10.3390/ani11102954

AMA Style

Pagaldai N, Arizaga J, Jiménez-Franco MV, Zuberogoitia I. Colonization of Urban Habitats: Tawny Owl Abundance Is Conditioned by Urbanization Structure. Animals. 2021; 11(10):2954. https://doi.org/10.3390/ani11102954

Chicago/Turabian Style

Pagaldai, Nerea, Juan Arizaga, María V. Jiménez-Franco, and Iñigo Zuberogoitia. 2021. "Colonization of Urban Habitats: Tawny Owl Abundance Is Conditioned by Urbanization Structure" Animals 11, no. 10: 2954. https://doi.org/10.3390/ani11102954

APA Style

Pagaldai, N., Arizaga, J., Jiménez-Franco, M. V., & Zuberogoitia, I. (2021). Colonization of Urban Habitats: Tawny Owl Abundance Is Conditioned by Urbanization Structure. Animals, 11(10), 2954. https://doi.org/10.3390/ani11102954

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop