Context

Regime shifts are well known for driving penetrating ecological change, yet we do not recognise the consequences of these shifts much beyond species diversity and productivity. Sound represents a multidimensional space that carries decision-making information needed for some dispersing species to locate resources and evaluate their quantity and quality.

Objectives

Here we assessed the effect of regime shifts on marine soundscapes, which we propose has the potential function of strengthening the positive or negative feedbacks that mediate ecosystem shifts.

Methods

We tested whether biologically relevant cues are altered by regime shifts in kelp forests and seagrass systems and how specific such shifted soundscapes are to the type of driver; i.e. local pollution (eutrophication) vs. global change (ocean acidification).

Results

Here, we not only provide the first evidence for regime-shifted soundscapes, but also reveal that the modified cues of shifted ecosystems are similar regardless of spatial scale and type of environmental driver. Importantly, biological sounds can act as functional cues for orientation by dispersing larvae, and observed shifts in soundscape loudness may alter this function.

Conclusions

These results open the question as to whether shifted soundscapes provide a functional role in mediating the positive or negative feedbacks that govern the arrival of species associated with driving change or stasis in ecosystem state.

Context

While remote sensing imagery is effective for quantifying land cover changes across large areas, its utility for directly assessing the response of animals to disturbance is limited. Soundscapes approaches—the recording and analysis of sounds in a landscape—could address this shortcoming.

Objectives

In 2011, a massive wildfire named “the Horseshoe 2 Burn” occurred in the Chiricahua National Monument, Arizona, USA. We evaluated the impact of this wildfire on acoustic activity of animal communities.

Methods

In 2013, soundscape recordings were collected over 9 months in 12 burned and 12 non-burned sites in four ecological systems. The seasonal and diel biological acoustic activity were described using the “Bioacoustic Index”, a detailed aural analysis of sound sources, and a new tool called “Sonic Timelapse Builder” (STLB).

Results

Seasonal biophony phenology showed a diurnal peak in June and a nocturnal peak in October in all ecological systems. On June mornings, acoustic activity was lower at burned than at non-burned sites in three of four ecological systems, due to a decreased abundance of cicadas directly impacted by the death of trees. Aural analyses revealed that 55% of recordings from non-burned sites contained insect sounds compared to 18% from burned sites. On October nights, orthopteran activity was more prevalent at some burned sites, possibly due to post-fire emergence of herbaceous.

Conclusions

Soundscape approaches can help address long-term conservation issues involving the responses of animal communities to wildfire. Acoustic methods can serve as a valuable complement to remote sensing for disturbance-based landscape management.

Context

Scale dependence of bat habitat selection is poorly known with few studies evaluating relationships among landscape metrics such as class versus landscape, or metrics that measure composition or configuration. This knowledge can inform conservation approaches to mitigate habitat loss and fragmentation.

Objectives

We evaluated scale dependence of habitat associations and scaling patterns of landscape metrics in relation to bat occurrence or capture rate in forests of southwestern Nicaragua.

Methods

We captured 1537 bats at 35 locations and measured landscape and class metrics across 10 spatial scales (100–1000 m) surrounding capture locations. We conducted univariate scaling across the 10 scales and identified scales and variables most related to bat occurrence or capture rate.

Results

Edge and patch density, at both landscape and class levels, were the most important variables across species. Feeding guilds varied in their response to metrics. Certain landscape and configuration metrics were most influential at fine (100 m) and/or broad (1000 m) spatial scales while most class and composition metrics were influential at intermediate scales.

Conclusions

These results provide insight into the scale dependence of habitat associations of bat species and the influence of fine and broad scales on habitat associations. The effects of scale, examined in our study and others from fine (100 m) to broad (5 km) indicate habitat relationships for bats may be more informative at larger scales. Our results suggest there could be general differences in scale relationships for different groups of landscape metrics, which deserves further evaluation in other taxonomic groups.

Context

Scale is the lens that focuses ecological relationships. Organisms select habitat at multiple hierarchical levels and at different spatial and/or temporal scales within each level. Failure to properly address scale dependence can result in incorrect inferences in multi-scale habitat selection modeling studies.

Objectives

Our goals in this review are to describe the conceptual origins of multi-scale habitat selection modeling, evaluate the current state-of-the-science, and suggest ways forward to improve analysis of scale-dependent habitat selection.

Methods

We reviewed more than 800 papers on habitat selection from 23 major ecological journals published between 2009 and 2014 and recorded a number of characteristics, such as whether they addressed habitat selection at multiple scales, what attributes of scale were evaluated, and what analytical methods were utilized.

Results

Our results show that despite widespread recognition of the importance of multi-scale analyses of habitat relationships, a large majority of published habitat ecology papers do not address multiple spatial or temporal scales. We also found that scale optimization, which is critical to assess scale dependence, is done in less than 5 % of all habitat selection modeling papers and less than 25 % of papers that address “multi-scale” habitat analysis broadly defined.

Conclusions

Our review confirms the existence of a powerful conceptual foundation for multi-scale habitat selection modeling, but that the majority of studies on wildlife habitat are still not adopting multi-scale frameworks. Most importantly, our review points to the need for wider adoption of a formal scale optimization of organism response to environmental variables.

Context

In the Rio Grande Plains of southern Texas, subtropical savanna vegetation is characterized by a two-phase pattern consisting of discrete woody patches embedded within a C₄ grassland matrix. Prior trench transect studies have suggested that, on upland portions of the landscape, large woody patches (groves) occur on non-argillic inclusions, while small woody patches (clusters) are dispersed among herbaceous vegetation where the argillic horizon is present.

Objective

To test whether spatial heterogeneity of subsurface soil texture drives the landscape-scale pattern of woody patches in this subtropical savanna.

Methods

Landscape-scale spatial patterns of soil texture were quantified by taking spatially-specific soil samples to a depth of 1.2 m in a 160 m × 100 m plot. Kriged maps of soil texture were developed, and the locations of non-argillic inclusions were mapped.

Results

Visual comparison of kriged maps of soil texture to a high resolution aerial photograph of the study area revealed that groves were present exclusively where the non-argillic inclusions were present. This clear visual relationship was further supported by positive correlations between soil sand concentration in the lower soil layers and total fine root biomass which mapped the locations of groves.

Conclusions

Subsurface non-argillic inclusions may favor the establishment and persistence of groves by enabling root penetration deeper into the profile, providing greater access to water and nutrients that are less accessible on those portions of the landscape where the argillic horizon is present, thereby regulating the distribution of grove vegetation and structuring the evolution of this landscape.

Context

Management of wintering waterfowl in North America requires adaptability because constant landscape and environmental change challenges existing management strategies regarding waterfowl habitat use at large spatial scales. Migratory waterfowl including mallards (Anas platyrhynchos) use the lower Mississippi Alluvial Valley (MAV) for wintering habitat, making this an important area of emphasis for improving wetland conservation strategies, while enhancing the understanding of landscape-use patterns.

Objectives

We used aerial survey data collected in the Arkansas portion of the MAV (ARMAV) to explain the abundance and distribution of mallards in relation to variable landscape conditions.

Methods

We used two-stage, hierarchical spatio-temporal models with a random spatial effect to identify covariates related to changes in mallard abundance and distribution within and among years.

Results

We found distinct spatio-temporal patterns existed for mallard distributions across the ARMAV and these distributions are dependent on the surrounding landscape structure and changing environmental conditions. Models performing best indicated seasonal surface water extent, rice field, wetland and fallow (uncultivated) fields positively influenced mallard presence. Rice fields, surface water and weather were found to influence mallard abundance. Additionally, the results suggest weather and changing surface water affects mallard presence and abundance throughout the winter.

Conclusions

Using novel datasets to identify which environmental factors drive changes in regional wildlife distribution and abundance can improve management by providing managers additional information to manage land over landscapes spanning private and public lands. We suggest our analytical approach may be informative in other areas and for other wildlife species.

Context

Climate change alters the vegetation composition and functioning of ecosystems. Measuring the magnitude, direction, and rate of changes in vegetation composition induced by climate remains a serious and unmet challenge. Such information is required for a predictive capability of how individual ecosystem will respond to future climates.

Objectives

Our objectives were to identify the relationships between 20 climate variables and 39 ecosystems across the southwestern USA. We sought to understand the magnitude of relationships between variation in vegetation composition and bioclimatic variables as well as the amount of ecosystem area expected to be affected by future climate changes.

Methods

Bioclimatic variables best explaining the plant species composition of each ecosystem were identified. The strength of relationships between beta turnover and bioclimate gradients was calculated, the spatial concordance of ecosystem and bioclimate configurations was shown, and the area of suitable climate remaining within the boundaries of contemporary ecosystems under future climate projections was measured.

Results

Across the southwestern USA, four climate variables account for most of the climate related variation in vegetation composition. Twelve ecosystems are highly sensitive to climate change. By 2070, two ecosystems lose about 4000 (15 %) and 7000 (31 %) km² of suitable climate area within their current boundaries (the Western Great Plains Sandhill Steppe and Sonora-Mojave Creosotebush-White Bursage Desert Scrub ecosystems, respectively). The climatic areas of riparian ecosystems are expected to be reduced by half.

Conclusions

Results provide specific climate and vegetation parameters for anticipating how, where and when ecosystem vegetation transforms with climate change. Projecting the loss of suitable climate for the vegetation composition of ecosystems is important for assessing ecosystem threats from climate change and for setting priorities for ecosystem conservation and restoration across the southwestern USA.

Context

Wildfire activity in boreal forests is projected to increase dramatically in response to anthropogenic climate change. By altering the spatial arrangement of fuels, land-cover configuration may interact with climate change to influence fire-regime dynamics at landscape and regional scales.

Objectives

We evaluate how land cover interacts with weather conditions to influence boreal-forest burning from 2012 to 2014 in Alaska.

Methods

Using geospatial fire and land-cover data, we quantify relationships between area burned and land cover, and test whether observed patterns of burning differ from random under varying weather conditions and fire sizes.

Results

Mean summer moisture index was correlated with annual area burned (ρ = ?0.78, p < 0.01), the total number of fires (ρ = ?0.68, p = 0.01), and the number of large fires (>500 km²; ρ = ?0.58, p = 0.04). Area burned was related positively to percent cover of coniferous forest and woody wetlands, and negatively to percent cover of shrub scrub, dwarf scrub, and open water and barren areas. Fires preferentially burned coniferous forest, which represented 50.1 % of the area burned in warmer/drier summers and 40.3 % of area burned in cooler/wetter summers, compared to the 34.5 % (±4.2 %) expected by random selection of land-cover classes. Overall vegetation tended to burn more similarly to random in warmer/drier than cooler/wetter years.

Conclusions

Land cover exerted greater influences on boreal fire regimes when weather conditions were less favorable for forest burning. Reliable projections of boreal fire-regime change thus require consideration of the interactions between climate and land cover, as well as feedbacks from land-cover change.

Purpose

Wildlife conservation requires understanding how landscape context influences habitat selection at spatial scales broader than the territory or habitat patch.

Objectives

We assessed how landscape composition, fragmentation, and disturbance affected occurrence and within-season site-fidelity of a declining grassland songbird species (Henslow’s Sparrow, Ammodramus henslowii).

Methods

Our study area encompassed eastern Kansas (USA) and North America’s largest remaining tracts of tallgrass prairie. We conducted 10,292 breeding-season point-count surveys over 2 years, and related occurrence and within-season site-occupancy dynamics of sparrows to landscape attributes within 400-, 800-, and 1600-m radii.

Results

Sparrows inhabited < 1% of sites, appearing and disappearing locally within and between breeding seasons. Early in spring, sparrows responded to landscape attributes most strongly within 400-m radii, settling in areas containing > 50% unburned prairie. Later in summer, sparrows responded to landscape attributes most strongly within 800-m radii, settling in areas containing > 50% unfragmented prairie, including sites burned earlier the same year. Sparrows avoided landscapes containing woody vegetation, disappeared from hayfields after mowing, and were most likely to inhabit landscapes containing Conservation Reserve Program (CRP) fields embedded within rangeland.

Conclusions

Landscape context influenced habitat selection at spatial scales broader than both the territory and habitat patch. Protecting contiguous prairies from agricultural conversion and woody encroachment, promoting CRP enrollment, and maintaining portions of undisturbed prairie in working rangelands each year are critical to reversing the conservation crisis in North America’s remaining grasslands. As landscape change alters natural areas worldwide, effective conservation requires suitable conditions for threatened species at multiple spatial scales.

Context

The problem of how ecological mechanisms create and interact with patterns across different scales is fundamental not only for understanding ecological processes, but also for interpretations of ecological dynamics and the strategies that organisms adopt to cope with variability and cross-scale influences.

Objectives

Our objective was to determine the consistency of the role of individual habitat patches in pattern-process relationships (focusing on the potential for dispersal within a network of patches in a fragmented landscape) across a range of scales.

Methods

Network analysis was used to assess and compare the potential connectivity and spatial distribution of highland fynbos habitat in and between protected areas of the Western Cape of South Africa. Connectivity of fynbos patches was measured using ten maximum threshold distances, ranging from five to 50 km, based on the known average dispersal distances of fynbos endemic bird species.

Results

Network connectivity increased predictably with scale. More interestingly, however, the relative contributions of individual protected areas to network connectivity showed strong scale dependence.

Conclusions

Conservation approaches that rely on single-scale analyses of connectivity and context (e.g., based on data for a single species with a given dispersal distance) are inadequate to identify key land parcels. Landscape planning, and specifically the assessment of the value of individual areas for dispersal, must therefore be undertaken with a multi-scale approach. Developing a better understanding of scaling dependencies in fragmenting landscapes is of high importance for both ecological theory and conservation planning.

Context

Past human land use has received increasing attention as an important driver of ecosystem change also in seemingly natural landscapes. Quantification of historical land use is therefore critical for assessing the degree of human impact and requires integration of ecology, history and archaeology.

Objective

This study aims to assess and compare levels of resource use by different actors during 355 years across a large landscape of northern Sweden.

Method

Data on resource use derived from case studies were extrapolated using demographic data to estimate harvested resources at the landscape scale. Here, we examined the use of the key-specie Scots pine by native Sami peoples and farmers and through commercial logging, and reconstructed historical forest conditions in order to interpret harvest levels and sustainability.

Results

We show that (1) the pre-industrial use of Scots pine resources in Pite Lappmark was sustainable from a landscape perspective, and (2) that the early commercial logging, in contrast, was not sustainable. Large and old Scots pine trees were logged at a very high rate, reaching up to 300 % of the annual ingrowth.

Conclusion

We suggest that historical landscape studies should incorporate analysis at different spatial scales, as such an approach can mirror the overall use of resources. Only then can land use data be applied across larger spatial scales, function as reference values and be compared to those of other regions, time-periods and types of human impact.

Context

Wildfires play a crucial role in maintaining ecological and societal functions of North American boreal forests. Because of their contagious way of spreading, using statistical methods dealing with spatial autocorrelation has become a major challenge in fire studies analyzing how environmental factors affect their spatial variability.

Objectives

We aimed to demonstrate the performance of a spatially explicit method accounting for spatial autocorrelation in burn rates modelling, and to use this method to determine the relative contribution of climate, physical environment and vegetation to the spatial variability of burn rates between 1972 and 2015.

Methods

Using a 482,000 km² territory located in the coniferous boreal forest of eastern Canada, we built and compared burn rates models with and without accounting for spatial autocorrelation. The relative contribution of climate, physical environment and vegetation to the burn rates variability was identified with variance partitioning.

Results

Accounting for spatial autocorrelation improved the models’ performance by a factor of 1.5. Our method allowed the unadulterated extraction of the contribution of climate, physical environment and vegetation to the spatial variability of burn rates. This contribution was similar for the three groups of factors. The spatial autocorrelation extent was linked to the fire size distribution.

Conclusions

Accounting for spatial autocorrelation can highly improve models and avoids biased results and misinterpretation. Considering climate, physical environment and vegetation altogether is essential, especially when attempting to predict future area burned. In addition to the direct effect of climate, changes in vegetation could have important impacts on future burn rates.

Context

A recent hypothesis, the habitat amount hypothesis, predicts that the total amount of habitat in the landscape can replace habitat patch size and isolation in studies of species richness in fragmented landscapes.

Objectives

To test the habitat amount hypothesis by first evaluating at which spatial scale the relationship between species richness in equal-sized sample quadrats and habitat amount was the strongest, and then test the importance of spatial configuration of habitat—measured as local patch size and isolation—when habitat amount was taken into account.

Methods

A quasi-experimental setup with 20 habitat patches of dry calcareous grasslands varying in patch size, patch isolation and habitat amount at the landscape scale was established in the inner Oslo fjord, Southern Norway. We recorded species richness of habitat specialists of vascular plants in equal-sized sample quadrats and analysed the relationship between species richness, habitat amount in the landscape and patch size and isolation.

Results

Although the total amount of habitat in a 3 km-radius around the local patch was positively related to species richness in the sample quadrats, local patch size had an additional positive effect, and the effect of patch size was higher when the amount of habitat within the 3 km-radius was high than when it was low.

Conclusions

In our study system of specialist vascular plants in dry calcareous grasslands, we do not find support for the habitat amount hypothesis.

Context

Land use change and forest degradation have myriad effects on tropical ecosystems. Yet their consequences for low-order streams remain very poorly understood, including in the world´s largest freshwater basin, the Amazon.

Objectives

Determine the degree to which physical and chemical characteristics of the instream habitat of low-order Amazonian streams change in response to past local- and catchment-level anthropogenic disturbances.

Methods

To do so, we collected field instream habitat (i.e., physical habitat and water quality) and landscape data from 99 stream sites in two eastern Brazilian Amazon regions. We used random forest regression trees to assess the relative importance of different predictor variables in determining changes in instream habitat response variables.

Results

Multiple drivers, operating at multiple spatial scales, were important in determining changes in the physical habitat and water quality of the sites. Although we found few similarities in modelled relationships between the two regions, we observed non-linear responses of specific instream characteristics to landscape change; for example 20 % of catchment deforestation resulted in consistently warmer streams.

Conclusions

Our results highlight the importance of local riparian and catchment-scale forest cover in shaping instream physical environments, but also underscore the importance of other land use changes and activities, such as road crossings and upstream agriculture intensification. In contrast to the property-scale focus of the Brazilian Forest code, which governs environmental regulations on private land, our results reinforce the importance of catchment-wide management strategies to protect stream ecosystem integrity.

Context

How do young birds achieve spatial knowledge about the environment during the initial stages of their life? They may follow adults, so gaining social information and learning; alternatively, young birds may acquire knowledge of the environment themselves by experiencing habitat and landscape features. If learning is at least partially independent of adults then young birds should respond to landscape composition at finer spatial scale than adults, who possess knowledge over a larger area.

Objectives

We studied the responses of juvenile, immature and adult Caspian Gull Larus cachinnans to the same habitat and landscape variables, but at several spatial scales (ranging from 2.5 to 15 km), during post-breeding period.

Methods

We surveyed 61 fish ponds (foraging patches) in southern Poland and counted Caspian gulls.

Results

Juvenile birds responded at finer spatial scales to the factors than did adults. Immature birds showed complicated, intermediate responses to spatial scale. The abundance of juvenile birds was mostly correlated with the landscape composition (positively with the cover of corridors and negatively with barriers). Adult abundance was positively related to foraging patch quality (fish stock), which clearly required previous spatial experience of the environment. The abundance of all age classes were moderately correlated with each other indicating that social behaviour may also contribute to the learning of the environment.

Conclusions

This study shows that as birds mature, they respond differently to components of their environment at different spatial scales. This has considerable ecological consequences for their distribution across environments.

Context

Forest insect outbreaks are influenced by ecological processes operating at multiple spatial scales, including host-insect interactions within stands and across landscapes that are modified by regional-scale variations in climate. These drivers of outbreak dynamics are not well understood for the western spruce budworm, a defoliator that is native to forests of western North America.

Objectives

Our aim was to assess how processes across multiple spatial scales drive western spruce budworm outbreak dynamics. Our objective was to assess the relative importance and influence of a set of factors covering the stand, landscape, and regional scales for explaining spatiotemporal outbreak patterns in British Columbia, Canada.

Methods

We used generalized linear mixed effect models within a multi-model interference framework to relate annual budworm infestation mapped from Landsat time series (1996–2012) to sets of stand-, landscape-, and regional-scale factors derived from forest inventory data, GIS analyses, and climate models.

Results

Outbreak patterns were explained well by our model (R ² = 93%). The most important predictors of infestation probability were the proximity to infestations in the previous year, landscape-scale host abundance, and dry autumn conditions. While stand characteristics were overall less important predictors, we did find infestations were more likely amongst pure Douglas-fir stands with low site indices and high crown closure.

Conclusions

Our findings add to growing empirical evidence that insect outbreak dynamics are driven by multi-scaled processes. Forest management planning to mitigate the impacts of budworm outbreaks should thus consider landscape- and regional-scale factors in addition to stand-scale factors.

Context

Landscape spatio-temporal heterogeneity is regarded as an important driver of biodiversity. In agricultural landscapes, the composition and configuration of cultivated fields and their multi-year dynamics should be considered. But the habitat-matrix paradigm in landscape ecology has resulted in little consideration of cropped areas.

Objectives

The main objective of our study was to determine the influences of spatial and multi-year temporal heterogeneity of the crop mosaic on carabid beetle assemblages of agricultural landscapes.

Methods

Carabids were sampled in 40 cereal fields in western France, and their species richness, total abundance and abundance of species groups with different dispersal abilities were measured. For each sampling site, we computed different metrics that characterized crop mosaic spatial and temporal heterogeneity. We quantified relationships between carabid assemblages and heterogeneity metrics and tested their significance.

Results

Total carabid abundance increased with increase in temporal heterogeneity of the crop mosaic. However, all species were not influenced in the same way by spatial and temporal heterogeneity metrics. Some species with high dispersal power such as Trechus quadristriatus were more abundant in landscapes with high spatial heterogeneity, whereas the abundance of less mobile species such as Poecilus cupreus were only positively influenced by temporal crop dynamics.

Conclusions

Our results suggest that both the spatial and temporal heterogeneity of the crop mosaic affects farmland biodiversity, at least for species that use crops during their life cycle or disperse through fields. We highlight the importance of taking this heterogeneity into account in further ecological studies on biodiversity in agricultural landscapes.

Context

Multi-objective management can mitigate conflicts among land-use objectives. However, the effectiveness of a multi-objective solution depends on the spatial scale at which land-use is optimized. This is because the ecological variation within the planning region influences the potential for site-specific prioritization according to the different objectives.

Objectives

We optimized the allocation of forest management strategies to maximize the joint production of two conflicting objectives, timber production and carbon storage, at increasing spatial scales. We examined the impacts of the extent of the planning region on the severity of the conflict, the potential for its mitigation, and the strategies that were identified as optimal.

Methods

Using forecasted data from a forest simulator, we constructed Pareto frontiers optimizing the joint provision of the objectives in production forests in Finland. Optimization was conducted within increasing hierarchical spatial scales and outcomes were compared in terms of the severity of the conflict and the solution to mitigate it.

Results

The trade-offs between timber production and carbon storage appeared less severe and could be mitigated more effectively the larger the planning regions were, but the improvements became minor beyond the scale of ‘large forest holding’. The results thus indicate that this scale, approximately 100 stands or 200 ha, is large enough to effectively mitigate the conflict between timber production and carbon storage.

Conclusions

Management planning over relatively small forest areas (200 ha) can mitigate ecosystem service trade-offs effectively. Thus the effective use of multi-objective optimization tools may be feasible even in small-scale forestry.

Context

Livestock predation by tiger and leopard in Bhutan is a major threat to the conservation of these felids. Conflict mitigation planning would benefit from an improved understanding of the spatial pattern of livestock kills by the two predators.

Objectives

We aimed to identify the landscape features that predict livestock kills by tiger and leopard throughout Bhutan. Our goals were to: (1) identify the predictors that have the largest influence in determining livestock kills, (2) assess the influence of scale across the different predictors evaluated and identify the scale at which each was most important.

Methods

We used livestock kills obtained from compensation records of tiger (n = 326) and leopard (n = 377) across Bhutan between 2003 and 2012 to run predation risk models with MaxEnt algorithm, using a multi-scale modeling approach (1, 2, 4, 8 and 16 km).

Results

Human-presence (density of settlements and roads) and land-cover (percentage of tree cover and meadow patches) were the main variables contributing to livestock kills by both species. Livestock kills were likely driven by a trade-off between livestock density and predator ecology, and the balance of this trade-off varied with scale. Risk maps revealed different hotspots for tiger and leopard kills, and analysis showed both species preferentially killed equids over other livestock types.

Conclusions

Our results highlight the importance of evaluating scale when investigating the spatial attributes of livestock kills by tiger and leopard. Our findings provide guidance for reducing conflict between humans and large felids throughout the country.

Context

Although multi-scale approaches are commonly used to assess wildlife-habitat relationships, few studies have examined selection at multiple spatial scales within different hierarchical levels/orders of selection [sensu Johnson’s (1980) orders of selection]. Failure to account for multi-scale relationships within a single level of selection may lead to misleading inferences and predictions.

Objectives

We examined habitat selection of the federally threatened eastern indigo snake (Drymarchon couperi) in peninsular Florida at the level of the home range (Level II selection) and individual telemetry location (Level III selection) to identify influential habitat covariates and predict relative probability of selection.

Methods

Within each level, we identified the characteristic scale for each habitat covariate to create multi-scale resource selection functions. We used home range selection functions to model Level II selection and paired logistic regression to model Level III selection.

Results

At both levels, EIS selected undeveloped upland land covers and habitat edges while avoiding urban land covers. Selection was generally strongest at the finest scales with the exception of Level II urban edge which was avoided at a broad scale indicating avoidance of urbanized land covers rather than urban edge per se.

Conclusions

Our study illustrates how characteristic scales may vary within a single level of selection and demonstrates the utility of multi-level, scale-optimized habitat selection analyses. We emphasize the importance of maintaining large mosaics of natural habitats for eastern indigo snake conservation.