Housing or habitat: what drives patterns of avian species richness in urbanized landscapes?

Conservation of biodiversity in cities depends on ecologically sensitive urban planning, informed by an understanding of patterns of species distributions and richness. Because urbanized landscapes are heterogeneous mosaics, and many species move between different land-cover types, it is valuable to compare ‘whole landscapes’ (broad-scale spatially heterogeneous areas) that systematically differ in landscape structure. We tested the relative influence of housing cover and canopy tree cover on avian species richness, to identify the components of landscape structure that most strongly influence landscape-scale richness (i.e., the pooled richness of multiple sites within a whole landscape). We selected 30 residential landscapes (each 1 km2) in Melbourne, Australia, stratified to represent concurrent gradients of housing and canopy tree cover. Five point-count surveys were conducted at each of 10 sites per landscape (for a total of 50 surveys per landscape) and the data pooled to represent the whole landscape mosaic. Up to 82% of variation in avian richness was explained by properties of the whole landscape. Housing cover was most dominant and a strong predictor for multiple response groups including native, terrestrial, forest, and aquatic birds. As housing cover increased, the richness of all groups decreased. Tree cover, primarily comprised of scattered trees in residential areas, had less influence on richness. Nonetheless, for forest birds, the extent of native vegetation surrounding a landscape had an important positive influence, indicating the value of potential source habitat for urban bird populations. Cities can be home to a diverse avifauna. The strong influence of landscape structure on species richness indicates a scope to plan and manage urbanized areas to support a diversity of birds that require natural habitat elements. We conclude that urbanizing environments can best be designed to benefit native birds by protecting patches of native vegetation (particularly large source areas) combined with localized higher housing cover, rather than uniform (lower) housing cover across the entire landscape.

richness (i.e., the pooled richness of multiple sites within a whole landscape). Methods We selected 30 residential landscapes (each 1 km 2 ) in Melbourne, Australia, stratified to represent concurrent gradients of housing and canopy tree cover. Five point-count surveys were conducted at each of 10 sites per landscape (for a total of 50 surveys per landscape) and the data pooled to represent the whole landscape mosaic. Results Up to 82% of variation in avian richness was explained by properties of the whole landscape. Housing cover was most dominant and a strong predictor for multiple response groups including native, terrestrial, forest, and aquatic birds. As housing cover increased, the richness of all groups decreased. Tree cover, primarily comprised of scattered trees in residential areas, had less influence on richness. Nonetheless, for forest birds, the extent of native vegetation surrounding a landscape had an important positive influence, indicating the value of potential source habitat for urban bird populations. Conclusions Cities can be home to a diverse avifauna. The strong influence of landscape structure on species richness indicates a scope to plan and manage urbanized areas to support a diversity of birds that require natural habitat elements. We conclude that urbanizing environments can best be designed to benefit native birds by protecting patches of native vegetation (particularly large source areas) combined with localized higher housing cover, rather than uniform (lower) housing cover across the entire landscape.

Introduction
Landscape change through urbanization typically involves the systematic replacement of native vegetation or farmland with built infrastructure, impervious surfaces, and high human population densities (Hahs and McDonnell 2006;Banville et al. 2017). The resulting urban areas represent complex, heterogeneous mosaics of different land-uses, including industrial land, residential neighborhoods, private gardens, open green spaces, riparian corridors, and small patches of native vegetation (Breuste et al. 2008). Despite such changes, cities can form valuable habitats for a variety of fauna, including some threatened species (Ives et al. 2016;Soanes and Lentini 2019).
Birds are a particularly ubiquitous, abundant, and diverse taxonomic group in urban areas, and play a key role in connecting urban residents with nature (Cox and Gaston 2016). Strong positive associations have been demonstrated between avian species richness and human mental health (Cox et al. 2017;Methorst et al. 2021a), life satisfaction (Methorst, Rehdanz, Methorst et al. 2021a, b) and happiness (Cameron et al. 2020). Clearly, there are ecological and social benefits to maintaining avian diversity in cities. Successful conservation, however, will depend on ecologically sensitive urban planning informed by an understanding of what drives the distribution of bird species and patterns of avian richness. The challenge lies in ensuring that cities can accommodate both a growing human population and diverse bird communities.
Much research to date has compared the value of different types of urban land-cover for birds: for example, gardens (Daniels and Kirkpatrick 2006;van Heezik et al. 2013), urban green spaces (Sandström et al. 2006;Rega et al. 2015), wooded streetscapes (White et al. 2005;Ikin et al. 2013), riparian areas (Banville et al. 2017), and patches of native vegetation (Sewell and Catterall 1998). Attention has also been given to specific habitat features, such as the abundance of significant large trees (Barth et al. 2015;Threlfall et al. 2017), diversity of tree species (Ferenc et al. 2014), and the extent of understorey vegetation (Daniels and Kirkpatrick 2006;Threlfall et al. 2017). Collectively, such work has enhanced understanding of how urbanization influences ecological patterns and processes (McDonnell and Hahs 2008;Alberti and Wang 2022), especially with regard to what makes urban habitats valuable for birds, and how habitat can be improved on a local scale.
These approaches to studying birds in urbanizing regions frequently share two characteristics. First, birds are surveyed along gradients of human influence, often an urban-rural gradient (McDonnell and Pickett 1990), which uses a space for time substitution to represent the process of urbanization (Pickett et al. 2011). Studies generally focus on a single broad measure of human influence, such as housing density (Pidgeon et al. 2014), human population density (Geschke et al. 2018), or distance from the city center (Connelly et al. 2020). However, the urban-rural gradient is complex (McDonnell and Pickett 1990), often inconsistently defined (Padilla and Sutherland 2019), and difficult to capture with a single measure (Padilla and Sutherland 2021). Recent studies (e.g., Padilla and Sutherland 2021;Alberti and Wang 2022) have demonstrated the benefits of considering multidimensional gradients of urbanization intensity.
Second, studies are typically site or patch-based, and data analysis and inference relate to the occurrence of bird species at the site or patch scale (e.g., White et al. 2005;Barth et al. 2015). Often such studies also include the influence of surrounding land-uses by establishing a buffer zone around individual sites, transects or patches (e.g., Mayorga et al. 2020;Malekian et al. 2021;Curzel et al. 2021). This approach, termed a 'patch-landscape' (McGarigal and Cushman 2002) approach, provides valuable evidence for the importance of landscape context in shaping avian communities. However, birds are mobile, and individuals often move between multiple sites or patches of the same, or different, habitat type within a broader heterogeneous landscape. Further, some land-uses may pose a barrier to movement within a landscape, or are of lower quality because they harbor predators, competitors, or human-induced disturbances (Shoffner et al. 2018).
Alternatively, in a 'landscape' approach (McGarigal and Cushman 2002;Bennett et al. 2006), a whole landscape is the unit of study (and of inference), with data pooled from multiple sites, transects or patches to represent the landscape. This approach recognises that landscape structure-the composition and configuration of different types of land-uses-has a profound effect on the ecological processes that shape the occurrence, abundance and richness of species (Wu and Hobbs 2002;Fahrig 2005). Using a 'whole of landscape' approach to examine how landscape structure influences faunal assemblages has provided important insights on the distribution of birds in agricultural (e.g., Haslem and Bennett 2008;Zingg et al. 2018;Sirami et al. 2019) and forested (e.g., McGarigal and McComb 1995;Sitters et al. 2014) environments. A landscape approach has also proven valuable in urban areas because it captures the inherent heterogeneity of cities and has clear applications in urban design (e.g., Belaire et al. 2014;Geschke et al. 2018;Souza et al. 2019). Further, the landscape as a spatial unit of inference is likely to be more relevant for urban planners, policy makers and local government management authorities than individual gardens, parks or streets; and therefore, offers novel insights for conserving birds in cities.
In this study, we used a landscape approach, whereby bird data were pooled from multiple sites within 100 ha residential landscapes, to investigate the influence of urban development on avian species richness in Melbourne, Australia. By studying birds in spatially heterogeneous landscapes (n = 30), explicitly stratified to sample concurrent gradients of housing cover and canopy tree cover, we aimed to: (i) investigate change in the biophysical properties of whole landscapes across an urbanization gradient; (ii) test the utility of a landscape approach with multiple gradients to study patterns of avian species richness in urbanized areas; (iii) determine the relative influence of housing cover, canopy tree cover, and other landscape properties on the richness of bird species; and (iv) compare the responses of different groups of birds to identify assemblages that are more (or less) vulnerable to urbanization. We predicted that: (i) the biophysical properties of landscapes would differ substantially across the urbanization gradient sampled; (ii) bird assemblages would respond strongly to the structure of 'whole' landscapes; (iii) avian species richness would decline with increasing housing cover and increase with canopy tree cover; and (iv) bird assemblages would respond to urbanization in different ways, with forest birds being most vulnerable due to their strong association with native vegetation.

Study area
The study area encompassed ~ 1300 km 2 in the northern and eastern suburbs of Greater Melbourne, a city of 4.92 million people in south-eastern Australia (Fig. 1). The selected area formerly supported native eucalypt forests and woodlands (Department of Environment Land Water and Planning 2021), but now comprises mostly residential land-use, interspersed with industry and commercial zones, urban parklands, riparian corridors, and intact patches of native vegetation. The climate is temperate: mean annual rainfall varies from 670 to 1260 mm across the region, with mean daily maximum temperatures of 13.4 °C in winter (July) and 26.7 °C in summer (January) (Bureau of Meteorology 2021).

Study design
We surveyed birds in 30 circular landscapes, each 1 km² in size (100 ha; 564 m radius), a landscape size widely used for studying birds (e.g., Zingg et al. 2018;Sambell et al. 2019;Sirami et al. 2019). Landscapes were selected to target heterogeneous residential areas (including small open green spaces, native vegetation patches and riparian areas), and avoided major roads, industrial zones, shopping centers, golf courses and public transport hubs. The 30 study landscapes each exhibited fine-scale variation, or heterogeneity, in land-cover types and biophysical properties (see Table 1 for examples). They were positioned ≥ 1 km apart and stratified to represent two main gradients of urbanization: housing cover and canopy tree cover.
For the housing cover gradient, we selected landscapes to represent three categories of housing density (see Theobald 2005): (i) Peri-urban (or Exurban, low density: 1-2 residential dwellings per ha (dph)); (ii) Suburban (medium density: 2-10 dph); and (iii) Urban (high density: >10 dph) (Appendix A). Our focus was on residential environments, so we did not sample inner City or outer Rural categories because they are more heavily dominated by commercial uses, or mixed lifestyle and agricultural land-uses, respectively. Sampling the full urbanization gradient (i.e., as depicted in Appendix A) may have reduced our capacity to detect trends across residential areas specifically, as subtle changes could have been diluted or masked by the significant differences regularly observed between the city center and undeveloped natural areas. For the canopy tree cover gradient, we selected landscapes representing a continuous gradient in the cover of canopy trees, spanning the maximum range available within each housing category (Peri-urban: 24-63%; Suburban: 20-51%; Urban: 13-37%). Tree cover was comprised mostly of canopy trees among residential land and urban parks. The number of residential dwellings per landscape was derived from the Vicmap Address-Address Point spatial layer (Department of Environment Land Water and Planning 2018) in ArcMap 10.7.1 (ESRI 2016), and canopy tree cover was initially visually estimated by observing aerial imagery (captured in 2017; sourced from Google Earth Pro 7.3.3.7786).
We selected an equal number of Peri-urban, Suburban and Urban landscapes (n = 10 each), ensuring half in each housing category (n = 5) contained a waterway, such as a natural flowing creek or drainage line. Waterways in Urban landscapes tended to be more modified than those in Peri-urban. Four landscapes were later reassigned once housing and canopy tree cover were systematically quantified (Table 1 and below). The final set of landscapes included 10 Periurban, 8 Suburban and 12 Urban (n = 30), spanning a housing gradient of 1.03-11.78 dph, and an estimated tree cover gradient of 13-63% cover (Appendix B).
Ten survey sites for birds were positioned in each landscape, distributed across four quadrants (≥ 2 sites per quadrant). The 10 sites were allocated among four land-cover types: residential, open green space, native vegetation patches, and riparian areas (see Fig. 2). In landscapes with waterways, two sites were allocated to riparian vegetation: one at the widest point, the second at a point where riparian width was considered representative of the landscape. For all other landcover types, sites were allocated in proportion to their cover in the landscape, with one site allocated per 10 ha (or ~ 10% cover). All sites were spaced > 200 m apart, and > 50 m from the landscape boundary.
Potential sites were marked on a map of each landscape, then checked in the field for suitability.

Bird surveys
Bird surveys, conducted by the same observer (JEH), consisted of a 10-min point count covering 0.5 ha (40 m radius) at each site. All birds seen or heard within the 40 m radius were considered to be 'on-site'; birds recorded > 40 m from the observer but within the landscape area were considered 'off-site'.
Surveys were conducted in fair weather, avoiding high winds (≥ 20 km/hour), rain and high temperatures (> 30 °C). All surveys within a landscape for a single round were completed on the same day (except in unsuitable weather). Five rounds of bird surveys were undertaken (1500 surveys in total) from 2018 to 2019, including three rounds in the breeding season (September to December) and two in the non-breeding season (March to August). During each season, sites were surveyed at least once in the early morning (0-3 h after sunrise) and once in the late morning (3-6 h after sunrise). The order in which sites and landscapes were visited was randomized for each round.
Bird data were pooled from 50 surveys per landscape (including both on-site and off-site records) to calculate landscape-level species richness for five response groups: all native species, terrestrial species, forest species, aquatic species, and exotic bird species. Terrestrial and aquatic species were subsets of all native species. Aquatic species included taxa typically associated with streams, wetlands and lakes. Forest species, a subset of terrestrial species, included taxa that depend on native vegetation. Species were classified as forest species if included in ≥ 55% of lists analysed by Fraser et al. (2015) and were also considered woodland, forest, or heathland-dependent by . Taxonomy followed Menkhorst et al. (2017).

Measures of human disturbance
Prior studies have indicated that the number of human pedestrians, domestic dogs, and domestic cats in urban areas can have a negative impact on avian communities (Villegas and Garitano-Zavala 2010; MacGregor-Fors and Schondube 2011; Fig. 1 a The study area in Melbourne, Australia, indicating the extent of land (light gray), water (dark gray), major roads (white lines), and the location of the 30 study landscapes (circles). b Aerial images of example landscapes from each housing category (Peri-urban in green, Suburban in blue, Urban in black) showing the upper (top row) and lower (bottom row) limits of canopy tree cover in each (indicated by the percentage value in the centre) ◂ Belaire et al. 2014). We incorporated these variables by counting the number of pedestrians, dogs and cats observed within a landscape during each visit. These direct counts were then averaged across the five rounds to give a mean value per landscape, per survey round.

Landscape properties
We collated variables to describe four broad properties of the study landscapes: human infrastructure, tree cover, landscape composition, and landscape context ( Table 1).
Many of these variables were based on the extent of different land-cover types, which we quantified by using a grid of points overlayed onto an aerial image of each landscape (captured in December  Fig. 3). Where multiple land-cover types were recorded for a single point (e.g., vegetation from a garden overhanging a road; a riparian area located within an open green space), the most dominant land-cover type was noted.
The origin of canopy trees may influence avian communities in urban areas, as native trees provide important feeding resources for many indigenous species (Young et al. 2007;Wood and Esaian 2020). For each point that fell on a canopy tree, we recorded whether it was native or exotic by comparing the colour of the foliage, the shape and density of the crown, and looking for evidence of leaf loss in autumn or winter. Most native canopy trees were Eucalypts which tend to have duller, grey-green foliage and more irregular and patchier crowns. The most The extent of each land-cover type was then visually estimated-here the land-cover types are digitised to show residential land (red), open green space (yellow), native vegetation (green), and riparian areas (blue). c Each site (yellow points) covered an area of 0.5 ha (40 m radius; solid yellow lines). Aerial images are provided for some example sites: c1 and c2 depict residential sites where the observer stood on the footpath; c3 shows an open green space site in a grassy area with a single eucalypt tree; c4 depicts a riparian site on the edge of a natural flowing creek common exotic canopy trees included pine (e.g., Pinus radiata), plane (Platanus spp.), and elm trees (Ulmus spp.) which have more vibrant, green foliage, denser crowns and were more likely to be deciduous. Where canopy trees could not be accurately assessed from aerial imagery alone, we used Google Street View or ground-truthing to identify trees to family level.
In addition, riparian zones in landscapes, and all patches of native vegetation within a landscape (> 1 ha) and within 2000 m of the landscape boundary (> 5 ha) were digitised and measured (Appendix C).
We represented human infrastructure at the landscape-scale by two spatially derived variables: total impervious cover (%) and housing cover (%) ( Table 1), both of which were highly correlated with the number of dwellings per landscape (r = 0.92 and r = 0.91, respectively). Tree cover was represented by two measures: total canopy tree cover (%) and the proportion of canopy tree cover comprised of exotic tree species (Table 1). Landscape composition was assessed by five metrics: open green space cover (%), riparian cover (ha), compositional diversity (heterogeneity; H), maximum patch size (native vegetation; ha), and mean patch size (native vegetation; ha) ( Table 1). We assessed landscape compositional diversity by using the Shannon-Wiener diversity index, based on the relative proportion of dominant land-cover types in each 100 ha landscape (residential, gardens, open green space, native vegetation patches, and riparian areas; Table 1). Five variables relating to landscape  (Table 1). We measured surrounding native vegetation because it may serve as a source area for forest bird species, the taxa of greatest conservation concern.
Data analysis

Assessment of survey effort
To assess bird survey effort, we generated a series of landscape-level species accumulation curves using the function 'specaccum' from the package 'vegan' (Oksanen et al. 2020). We used both the random permutation and collector methods to visualize the accumulation of species over 50 surveys per landscape. For many Suburban and Urban landscapes, species richness plateaued after 30 to 40 surveys; Peri-urban landscapes approached a plateau by 50 surveys (see Appendix D for example curves).

Landscape properties across an urbanization gradient
To examine variation in landscape properties across the gradient sampled, we carried out an ordination using non-metric multidimensional scaling (NMDS), based on a Bray-Curtis dissimilarity matrix of the landscape variables in Table 1. Landscape variables were fitted as environmental vectors using the function 'envfit' from the package 'vegan' (Oksanen et al. 2020). We also compared landscape properties between Peri-urban, Suburban and Urban landscapes using one-way ANOVAs and Tukey HSD post-hoc tests.

Predictors of avian species richness at the landscape-scale
We tested correlations between all pairs of landscape variables; where r s >0.61, one variable was excluded from subsequent analyses (Appendix E) (Dormann et al. 2013). Despite efforts to select landscapes along independent, concurrent gradients of housing cover and canopy tree cover, there was a moderate correlation between these variables (r = − 0.58); however, this was not considered strong enough to distort models (Dormann et al. 2013). Variables retained for modelling were those critical to the study design (i.e., housing cover, canopy tree cover) and those considered likely to have ecological influence (e.g., extent of surrounding native vegetation). The final set of variables included: housing cover, canopy tree cover, exotic tree proportion, open green space cover, riparian cover, rainfall long term, and surrounding native vegetation within 1000 m. All variables were standardised prior to further analysis.
Generalized linear models, assuming a Poisson distribution, were used to test hypotheses concerning the relative influence of landscape properties on avian species richness at the landscape-scale. We treated each of the four groups of variables as independent hypotheses which may explain variation in bird species richness: (i) Human infrastructure (Housing); (ii) Tree cover (Tree Cover + Tree Exotic); (iii) Landscape composition (Open Green + Riparian); and (iv) Landscape context (Rain Long + Vege 1000) (see Table 1) We predicted that avian species richness would be highest in landscapes with: (i) less housing cover; (ii) more canopy tree cover, but fewer exotic trees; (iii) more open green space and riparian areas; and (iv) higher annual rainfall and a greater extent of surrounding native vegetation.
We used a global model and the function 'dredge' from the package 'MuMIn' (Barton 2020), to compare all possible combinations of these hypotheses, for a total of 16 candidate models for each of the five bird response groups (native, terrestrial, forest, aquatic, exotic). The function 'subset' was used to ensure predictor variables representing the same hypothesis were included together (e.g., Tree Cover + Tree Exotic) (Barton 2020). We then ranked candidate models using Akaike's information criterion corrected for small sample sizes (AIC c ), to identify the most parsimonious model for each bird group (Hurvich and Tsai 1989). Models were deemed similarly plausible when the difference between the AIC c value (∆ i ) and the 'best' model (smallest AIC c value) was less than two (∆AIC c < 2.0) (Burnham and Anderson 2002). Where multiple models had ∆AIC c < 2.0, all associated variables were included in a single final model. For each final model, we plotted the residuals, calculated the dispersion ratio, and used Cook's distance to check the leverage of individual observations (Cook and Weisberg 1984); there was no evidence of overdispersion and no landscapes with high leverage. For each bird response group, we plotted the modelled relationships (using 'ggplot2' (Wickham 2016)) for all variables with moderate to strong evidence for an effect on species richness (i.e., P < 0.05). All data analyses were carried out in R version 4.0.1 (R Core Team 2016) and R Studio version 1.1.419 (RStudio Team 2016). We interpret the analyses using the language of evidence, as described by Muff et al. (2022).

Landscape properties across an urbanization gradient
An NMDS ordination of the 30 landscapes (Fig. 4a) revealed that the main ordination axis (NMDS1) broadly represents a gradient in urban intensification. Low values of NMDS1 represent landscapes with greater cover of impervious surfaces and housing, a higher proportion of exotic tree cover, and more people observed per visit (Fig. 4a). In contrast, high values of NMDS1 represent landscapes with greater canopy tree cover, larger size (mean and maximum) of native vegetation patches, and a greater extent of native vegetation surrounding the landscape (Fig. 4a). All landscape variables, except the direct count of domestic dogs, explained a proportion of the variation observed among the 30 landscapes (Appendix F).
Most landscape variables showed clear differences across housing categories, with the strongest evidence for differences found between the low (Periurban) and high (Urban) ends of the gradient (Fig. 4b, Appendix G). Impervious surface cover, housing cover, people, and exotic tree proportion increased across the housing cover gradient; canopy tree cover, maximum native patch size, and the extent of surrounding native vegetation (within 1000 m, 1500 m, and 2000 m) decreased (Fig. 4b; Appendix G).

The bird community
A total of 76 bird species was recorded, including 71 native and five exotic species (Appendix H). Ten native species were recorded in all 30 landscapes, including Australian magpie (Gymnorhina tibicen) and noisy miner (Manorina melanocephala). Two exotic species, common myna (Acridotheres tristis) and spotted dove (Spilopelia chinensis), also occurred in all landscapes. In contrast, 10 species were recorded in just one landscape: these included aquatic species such as Australasian darter (Anhinga novaehollandiae) and forest birds such as scarlet robin (Petroica boodang).

Predictors of avian species richness at the landscape-scale
The best performing models for native species (∆AICc < 2.0) included measures of human infrastructure (Housing) and landscape context (Rain Long, Vege 1000) (Appendix I). The final model explained 82% of the variation in native species richness (Table 2). There was strong evidence (P < 0.01) for a negative effect of housing cover on native species richness at the landscape-scale (Fig. 5a).
For terrestrial species richness, all seven landscape properties were represented in the best performing models (Appendix I). The final (global) model explained 80% of the variation in terrestrial species richness (Table 2), with strong evidence (P < 0.01) that landscapes with greater housing cover have lower richness of terrestrial bird species (Fig. 5b).
The best performing models for forest species included measures of human infrastructure (Housing), landscape composition (Open Green, Riparian) and landscape context (Rain Long, Vege 1000) Fig. 4 a A non-metric multidimensional scaling ordination of the 30 study landscapes based on the full complement of biophysical landscape properties, including measures of human infrastructure, tree cover, landscape composition, and landscape context. Environmental vectors are shown as grey lines. b Dot plots of the mean values for Peri-urban, Suburban and Urban landscapes for all variables where ANOVAs and Tukey HSD post-hoc tests showed moderate to strong evidence for a difference. The top row displays variables that increased across the housing cover gradient; the bottom row displays variables that decreased. The letter codes indicate differences between housing categories for which post-hoc tests showed moderate to strong evidence for a difference. See Table 1 for a description of variables and abbreviations of variable names ◂ (Appendix I). This combination of landscape properties explained 82% of the variation in forest bird richness (Table 2). There was strong evidence (P < 0.01) for a negative influence of housing cover on the richness of forest species in urbanizing landscapes, and moderate evidence (P < 0.05) for a positive relationship with the extent of surrounding native vegetation (Fig. 5c).
For aquatic species, the best performing models included measures of human infrastructure (Housing), tree cover (Tree Cover, Tree Exotic) and landscape composition (Open Green, Riparian) (Appendix I). The final model explained 63% of the variation (Table 2), with strong (P < 0.01) and moderate evidence (P < 0.05), respectively, that housing cover and canopy tree cover negatively affect the richness of aquatic birds in urbanized landscapes (Fig. 5d).
The best performing models for exotic species richness at the landscape scale included only human infrastructure (Housing) (Appendix I); however, the null model was equally parsimonious. There was little evidence for an effect of housing cover on the richness of exotic birds, and this model explained only 5% of the variation (Appendix I).

Discussion
This study investigated change in avian species richness along gradients of housing cover and canopy tree cover in a large, urbanized region. Importantly, we examined how the properties of 'whole' landscapes (1 km 2 ) affect bird communities. There were marked differences in the spatial patterns of land cover between Peri-urban, Suburban and Urban landscapes. Species richness of several bird habitat groups also varied greatly between landscapes. Housing cover, the proportion of the residential landscape dominated by buildings, negatively influenced the richness of all bird groups, except for exotic species which were ubiquitous across the region. Forest species were also positively influenced by the extent of native vegetation surrounding the landscape, suggesting that landscape context plays a key role for this assemblage. Landscape properties explained > 80% of the variation in richness for native, terrestrial and forest birds, demonstrating the value of a landscape approach for understanding patterns of avian species richness in urbanized areas. Our findings indicate there is scope to enhance habitat for avifauna in urban landscapes via improved management and planning.

Landscape change across an urbanization gradient
Landscapes sampled in this study showed complex patterns of variation in multiple landscape properties. The primary gradient contrasted landscapes with greater housing cover against those with more native vegetation cover. Landscapes with greater housing cover had greater impervious surface cover and more pedestrians observed during survey rounds. They also tended to have a greater proportion of exotic tree cover, reflecting the composition of residential gardens in older, more densely settled suburbs. Landscapes with greater canopy tree cover were characterised by larger patches of native vegetation, both within the landscape and at buffer distances of up to 2 km. This primary urbanization gradient aligns with trends observed in other regions, such as cities in North America (Melles et al. 2003;Padilla and Sutherland 2021), South America (Souza et al. 2019), and Europe (Suarez-Rubio and Krenn 2018). Housing cover values for Greater Melbourne (for 2018) were comparable to many cities in the conterminous United States in 2000 (Theobald 2005).
The second gradient identified by the ordination of landscape attributes, was represented by measures of landscape composition, notably the cover of open green space, cover of riparian vegetation and landscape heterogeneity. In residential environments, green spaces such as sports ovals, local parks and undeveloped cleared areas, along with retained streamside vegetation add visible heterogeneity to the landscape, a factor previously identified as an important feature of urban landscapes for birds (Souza et al. 2019). These findings, that urban landscapes encompass variation in multiple landscape properties, underscore the value of a landscape approach for understanding drivers of avian biodiversity.
Housing cover, as a key determinant of the structure and composition of urbanized landscapes, Fig. 5 Predicted relationships between the richness of bird response groups and landscape properties for which there was strong evidence: a native species, b terrestrial species, c forest species, and d aquatic species. For exotic birds, there was no evidence for a relationship between species richness and any of the landscape variables included in this study ◂ provides a useful measure for quantifying urban intensification (e.g., Pidgeon et al. 2014), as it is strongly correlated with several other landscape attributes (e.g., impervious surfaces, density of people), which have the potential to influence species richness at the landscape-scale. Despite deliberately designing our study to sample independent gradients of both housing cover and tree cover across urbanized landscapes, there was a moderate association between these variables (r = − 0.58). This suggests it is difficult to maintain high levels of canopy tree cover (e.g., > 40% landscape cover) in areas with greater housing cover. Conversely, at least in our study region, it was rare that landscapes with low housing cover also had low canopy tree cover (< 25%). This is likely because the study region was formerly covered by forest vegetation (Department of Environment Land Water and Planning 2021), and much land clearing has occurred due to urban expansion. This contrasts with other areas where urban development has expanded into cleared agricultural land, where prior tree cover is relatively low. Land-use history can have a strong influence on spatial patterns and conservation outcomes in urbanizing landscapes (Padilla and Sutherland 2021;Alberti and Wang 2022) and other environments (e.g., Lunt and Spooner 2005). Relationships between avian species richness and landscape properties Changes in the biophysical properties of landscapes influence the extent and quality of habitat available for species in urban areas (McKinney 2002). This occurs at multiple spatial scales: both local and landscape-scale habitat variables are influential (Melles et al. 2003;Mayorga et al. 2020;Curzel et al. 2021).
Here, the extent of landscape-scale housing cover was a primary driver of the availability and suitability of habitat for native birds: as housing cover increased, fewer bird species were recorded in urbanized landscapes. These findings are consistent with other reports of a negative relationship between bird species richness and housing density (Pidgeon et al. 2014;Wood et al. 2015), or the extent of impervious surfaces (Chamberlain et al. 2019;Souza et al. 2019;Lerman et al. 2021;Suárez-Castro et al. 2022). The response to housing cover is not necessarily a direct relationship; rather, that increasing housing cover serves as a proxy for a range of other processes. Built structures such as houses, shops, roads and carparks all displace natural habitat and reduce the land available for land-cover types such as gardens, parks and patches of native vegetation which collectively provide more suitable and diverse habitat for most birds in urban landscapes (Sewell and Catterall 1998;Daniels and Kirkpatrick 2006;Sandström et al. 2006). Greater housing cover is also associated with increases in other human-induced disturbances, such as vehicle traffic and road strikes, domestic cats, air pollution, anthropogenic noise, and artificial light at night (Ciach and Fröhlich 2017;Shoffner et al. 2018). Relative to housing cover, landscape-scale canopy tree cover did not have a strong influence on the richness of native, terrestrial or forest birds across the urbanization gradient sampled. This contrasts with other studies in urban environments which have reported a positive relationship between avian species richness and canopy tree cover, either for landscape units (Chamberlain et al. 2019;Callaghan et al. 2020) or sites with surrounding tree cover (Callaghan et al. 2018;Lerman et al. 2021). The lack of a relationship between landscape-scale tree cover and bird richness here may be due to several factors. First, it is likely that the form of tree cover in these predominantly residential landscapes affects the degree to which it influences bird species. Tree cover mostly consisted of street trees and trees scattered among residential gardens and open parkland, rather than substantial intact patches of native vegetation. Whilst scattered trees undoubtedly provide habitat and feeding resources in urban areas (Wood and Esaian 2020), they are potentially of lesser value than intact patches of native forest (e.g., Haslem and Bennett 2008).
The origin of tree species in the urban environment is also important. In many landscapes, exotic trees contributed substantially to the overall tree cover. Whilst some bird species use resources provided by a combination of native and exotic plants (French et al. 2005;Daniels and Kirkpatrick 2006), a prevalence of introduced tree species reduces the amount of habitat for certain assemblages (Wood and Esaian 2020), such as insectivores and nectarivores that feed in eucalypt canopies (White et al. 2005;Young et al. 2007;Ikin et al. 2013). Third, this study sampled a gradient focussed primarily on residential landscapes, avoiding high-density inner city and nearnatural landscapes at the urban fringe. It is likely that if a longer gradient was sampled, including the more extreme situations at either end, landscape-level tree cover would have had a greater influence on avian species richness. Finally, the positive influence of canopy tree cover may simply have been masked by the moderate negative correlation with housing cover.
Bird species that depend on natural habitats are considered 'urban avoiders' and are among the most sensitive to changes associated with urbanization (McKinney 2002;Conole and Kirkpatrick 2011;Ferenc et al. 2014;Shoffner et al. 2018). In southern Australia, woodland-dependent birds are a group of conservation concern due to the extensive loss and fragmentation of wooded habitats throughout their range Ford 2011). Here, our results show that in addition to a negative response to increasing housing cover, forest-dependent birds were positively associated with surrounding tree cover (within 1000 m of the landscape boundary) comprising native vegetation in intact patches (≥ 5 ha). These results are consistent with the view that forest species are frequently lost from peri-urban suburbs and that forest remnants, either within or close to residential areas, are essential for the persistence of many species in this assemblage (Sewell and Catterall 1998;Palmer et al. 2008;Geschke et al. 2018;Shoffner et al. 2018). It is likely that intact forest areas adjacent to the study landscapes act as source populations for forest species from which dispersal and spill over of individuals occurs into neighboring residential areas (see also van Heezik et al. 2010;van Heezik and Adams 2016).
The value of a landscape approach to understanding avian diversity in urban environments We sampled birds at multiple sites and pooled the data to represent the richness of avifauna for 'whole' landscapes. There was much variation in richness between these residential landscapes (e.g., 5-25 species per landscape for forest birds); and importantly, such variation was strongly related to properties of the landscape, indicating that this is a valuable scale for seeking understanding. Between 80 and 82% of the variation in richness of all native species, terrestrial species, and forest species, respectively, was explained by landscape-scale properties. Less variation was explained for aquatic species (63%), which were not a primary focus, or for exotic species (5%) which were ubiquitous throughout the region.
There are multiple benefits of a landscape approach. First, it provides a means of measuring and understanding the consequences of the inherent heterogeneity of cities, which are complex mosaics of different land-uses (Breuste et al. 2008). Here, even primarily residential landscapes showed variation in multiple landscape properties (see also Padilla and Sutherland 2021), and our results provide an understanding of how bird groups respond to different aspects of habitat modification and urbanization. Second, a landscape approach explicitly recognises that bird species may move between different parts of the landscape and that ecological processes operate at spatial scales that incorporate multiple habitats or land-use types. For example, flowering of eucalypt trees in these landscapes provides a seasonal nectar source that some bird species exploit, and to do so they move between trees in residential gardens, along streets or in remnant patches of vegetation. Third, the findings from studies at this scale are highly relevant to urban planners, policy makers and local government management authorities, who must plan at broad scales rather than individual sites. Finally, the strong response of bird species to the overall properties of a landscape implies that individual actionsby local residents, community groups or government -matter, because they contribute to change in the wider landscape. For example, setting aside a local area from housing to restore native bushland will alter the landscape structure, and influence the avifauna of the broader landscape.

Conclusion
Urbanizing landscapes in major cities are complex mosaics which have differing landscape properties, including varying extents of human infrastructure (e.g., housing cover) and tree cover (all tree cover, exotic trees). In residential landscapes across Melbourne, the extent of housing cover had the strongest influence on avian species richness at the landscape-scale. Residential housing offers little value as direct habitat for birds, and its presence restricts the amount of space available for other land-uses, such as gardens, parks or native vegetation, which offer more suitable habitat. Additionally, the limited canopy tree cover that remains in highly developed residential landscapes is often dominated by introduced tree species, which offer fewer resources for native avifauna. Heavily urbanized areas often lack nearby large patches of native vegetation which are essential for the persistence of forest birds.
Continued expansion of urbanized areas will come at a cost to biodiversity. The amount of land available is finite, and so there must be trade-offs between land allocated to housing and to other landcover types such as native vegetation. Our findings suggest that benefits for native birds could best be achieved by protecting native vegetation, particularly source areas of intact habitat, and compensating with local areas of higher housing cover, rather than having uniform lower housing cover across the entire landscape. feedback which substantially improved this manuscript. The study was conducted with approval of the La Trobe University Animal Ethics Committee (Approval No. AEC18-11) and the Department of Environment, Land, Water and Planning (Permit No. 10008805 Funding Open Access funding enabled and organized by CAUL and its Member Institutions. This study was funded by the Holsworth Wildlife Research Endowment and the Ecological Society of Australia, BirdLife Australia, the Australian Wildlife Society, and the Field Naturalists Club of Victoria. JEH was supported by an Australian Government Research Training Program (RTP) Stipend Scholarship.

Data availability
The datasets analysed in the present study are available from the corresponding author upon reasonable request.

Declarations
Competing interests Financial interests: The authors declare they have no known financial interests. Non-financial interests: Andrew F. Bennett is a current editor with the journal Landscape Ecology.

Ethical approval
The study was conducted with approval of the La Trobe University Animal Ethics Committee (Approval No. AEC18-11) and the Department of Environment, Land, Water and Planning (Permit No. 10008805).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.