Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Does Fire Influence the Landscape-Scale Distribution of an Invasive Mesopredator?

  • Catherine J. Payne,

    Affiliation Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Melbourne, Victoria, Australia

  • Euan G. Ritchie,

    Affiliation Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Melbourne, Victoria, Australia

  • Luke T. Kelly,

    Affiliation Australian Research Council Centre of Excellence for Environmental Decisions, School of Botany, University of Melbourne, Melbourne, Victoria, Australia

  • Dale G. Nimmo

    dale@deakin.edu.au

    Affiliation Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Melbourne, Victoria, Australia

Abstract

Predation and fire shape the structure and function of ecosystems globally. However, studies exploring interactions between these two processes are rare, especially at large spatial scales. This knowledge gap is significant not only for ecological theory, but also in an applied context, because it limits the ability of landscape managers to predict the outcomes of manipulating fire and predators. We examined the influence of fire on the occurrence of an introduced and widespread mesopredator, the red fox (Vulpes vulpes), in semi-arid Australia. We used two extensive and complimentary datasets collected at two spatial scales. At the landscape-scale, we surveyed red foxes using sand-plots within 28 study landscapes – which incorporated variation in the diversity and proportional extent of fire-age classes – located across a 104 000 km2 study area. At the site-scale, we surveyed red foxes using camera traps at 108 sites stratified along a century-long post-fire chronosequence (0–105 years) within a 6630 km2 study area. Red foxes were widespread both at the landscape and site-scale. Fire did not influence fox distribution at either spatial scale, nor did other environmental variables that we measured. Our results show that red foxes exploit a broad range of environmental conditions within semi-arid Australia. The presence of red foxes throughout much of the landscape is likely to have significant implications for native fauna, particularly in recently burnt habitats where reduced cover may increase prey species’ predation risk.

Introduction

Predators shape ecosystems worldwide [1]. They can exert top-down regulation of lower trophic levels [2] and induce trophic cascades which flow through entire ecosystems [3]. Predators introduced to areas outside of their native range can have a particularly strong effect on native species [4], and have caused population declines and extinctions in a range of ecosystems [5]. Many invasive predators are ‘mesopredators’: smaller predator species that increase in abundance or activity following the removal of apex predators [6]. For example, in Australia, persecution of the native apex predator, the dingo (Canis dingo), has led to increases in the density or activity of invasive mesopredators (e.g. the red fox [Vulpes vulpes]) throughout large portions of the continent [3].

Fire is another globally significant process that affects environments worldwide [7]. Fire influences ecosystems via bottom-up control by altering the availability of key resources for biota. Fire incinerates plant matter, altering vegetation structure [8], [9], which in turn affects the distribution and abundance of animals [10].

Invasive mesopredators and fire share an important characteristic from a conservation perspective: both can be manipulated through management interventions. Invasive mesopredators are managed using lethal control and exclusion fencing, and fire using suppression or prescribed burning. However, management of mesopredators and fire usually occurs in isolation, without consideration of the potential effects of fire on mesopredators [11]. It is important to rapidly address this significant knowledge gap because some fire regimes may exacerbate the effects of invasive mesopredators by simplifying vegetation and amplifying predation risk [12], [13]. For example, interactions between fire regimes and invasive mesopredators have been hypothesised as a cause of lower survival of reptile species in recently-burned areas [14], and a contributor to the collapse of small mammal communities in northern Australia [15].

The red fox is one of the world’s most widely distributed mesopredators. It is common in both the northern and southern hemispheres. Foxes, and a second introduced mesopredator, the feral cat (Felis catus), are widely regarded as the primary cause of extinctions and declines of Australia’s marsupial fauna [5]. Evidence for the negative impact of foxes has been demonstrated through predator-control experiments that have shown that prey species increase in both range and activity when foxes are removed [16], [17]. Further evidence comes from dietary studies showing foxes eat a wide range of native mammal, reptile, bird, and invertebrate prey [18][20].

Despite indications that foxes may inhibit the recovery of native species following fire [12], [21], whether foxes are themselves influenced by fire remains poorly known. This knowledge gap limits the ability of land managers to consider the effects of fire management on red foxes, which could have negative ramifications for native biodiversity. While foxes are widely considered as habitat generalists, they do display local variability in occurrence related to habitat or landscape structure [22]. For example, in some regions, foxes prefer heterogeneous landscapes [22], as they are able to use multiple landscape elements on a daily or seasonal basis [23], [24]. Fire management in many regions seeks to maximise landscape heterogeneity by creating mosaics of fire ages (i.e. ‘patch mosaic burning’; [25]). Does such management inadvertently favour invasive mesopredators?

The few studies that have explored the topic have focused on relatively short temporal scales (<30 years and often <10 years post fire) or small spatial scales (but see [26]). However, in some ecosystems, post-fire vegetation recovery continues for a century or more after fire [27]. Consequently, animal species respond to fire over similarly long time-frames [28]. The effects of fire can also occur across multiple spatial scales [29]; while time since fire may affect a species’ occurrence at any point in the landscape, the area and composition of fire-ages within a ‘whole’ landscape can play a critical role in affecting species’ landscape-level distributions [30]. This is likely to be especially true for large, mobile species, such as the red fox.

In addition to the effects that fire may have on species’ occurrence, other environmental factors may be locally important. With regard to foxes, this includes climate [26], the distribution of vegetation types [31], and the distance to roads [24] and agricultural land [22]. Foxes rely on free standing water for drinking, particularly when temperatures are high (>30°C), as is common in many semi-arid environments. Hence, as annual rainfall decreases (aridity intensifies) permanent water may be reduced in its availability and limit fox occurrence. Foxes are often thought of as edge specialists [22]. They often prefer to hunt in open areas such as resource-rich agricultural fields or structurally simple vegetation types adjacent to more complex vegetation which provides cover during the day [22], [32]. Their ability to hunt may be further enhanced where roads create easy access and increased visibility in otherwise structurally complex habitats [24], [33].

Here, we examine what drives the occurrence (reporting rate) of red foxes in semi-arid Australia at multiple spatial scales, with a particular emphasis on the role of fire. We conducted two large-scale natural experiments. First, we explored landscape-scale patterns of fox occurrence in relation to the properties of fire mosaics; namely, the amount and diversity of fire age-classes within each of 28 study landscapes (each 12.6 km2). Second, we explored site-scale patterns of fox occurrence in relation to fire history at 108 sites stratified along a century-long post-fire chronosequence. In both cases, we also quantified the influence of other environmental variables such as vegetation type and distance to agricultural land. Our aims were: 1) to determine the drivers of fox distribution in semi-arid Australia; and 2) to understand the specific role of fire in influencing fox occurrence at large scales relevant to fire and mesopredator management.

Materials and Methods

Study region

This study was undertaken in the Murray Mallee region of south-eastern Australia (Fig. 1). The climate in the region is semi-arid, with mean annual rainfall of 200–350 mm and average daily maximum temperatures are 30–33°C in summer and 15–18°C in winter (Australian Bureau of Meteorology; http://www.bom.gov.au). The vegetation is predominantly ‘tree mallee’ characterised by an overstorey of Eucalyptus species (<5–8 m) with a multi-stemmed growth form [34]. Two vegetation types are common throughout region [35]. ‘Triodia Mallee’ has a canopy of Eucalyptus dumosa and E. socialis with an understorey of Triodia scariosa and mixed shrubs, and occurs mainly on sandier soils typical of dunes. ‘Chenopod Mallee’ has a canopy of E. oleosa and E. gracilis with an open understorey of chenopod species, and occurs on heavier soils typical of swales.

thumbnail
Figure 1. Map of study area showing all landscapes (circles) considered in this study (grey shading indicates mallee vegetation; majority of white areas indicates agricultural land used for grazing and cereal crops).

The dashed box shows the spatial extent of the site-scale study. An inset shows an example of a study landscape including the position of 10 sites within where site-scale data were collected. Within the inset, different hatching represents different fire ages.

https://doi.org/10.1371/journal.pone.0107862.g001

Mallee vegetation is fire-prone with large fires (i.e. >100,000 ha) occurring somewhere in the region on a bidecadal basis [36], although individual sites can go long periods without fire (i.e. >100 years; [37]). Fire is actively managed in the region through prescribed burning and suppression for both asset protection and conservation objectives [36]. Most wildfires are ignited by lightning strikes and are stand-replacing, essentially resetting vegetation succession to ‘year-zero’ (Fig. 2; [8]).

thumbnail
Figure 2. Examples of mallee vegetation with differing fire histories.

(a) A recently burned site; (b) A long unburned site.

https://doi.org/10.1371/journal.pone.0107862.g002

Site selection

We refer to two datasets in this study derived from two different natural experiments that differed in both their spatial grain and extent. We refer to these as ‘landscape-scale’ and ‘site-scale’ datasets throughout, in reference to the spatial scale of the response and predictor variables (i.e. the spatial grain) of the respective datasets.

Landscape-scale data.

The landscape-scale dataset consists of 28 study landscapes, each with a 4 km diameter circle (12.6 km2; Fig. 1), distributed throughout a 104, 000 km2 study area. These landscapes were selected as part of a broad-scale natural experiment: the Mallee Fire and Biodiversity Project. Study landscapes were selected to allow a comparison of the effects of different approaches to patch mosaic burning on biodiversity, with a particular emphasis on the role of the area and diversity of fire-ages (‘pyrodiversity’, see [25]). Thus, landscapes were stratified according to number and spatial extent of fire-age classes within the landscape [29]. The fire history of the region was mapped using the ENVI package [38] and then converted to shape files for use in ArcMap version 9.2 [39]. Only fires that occurred post-1971 were mapped due to limited availability of Landsat imagery prior to this time (see [36]).

Site-scale data.

We collected site-scale data within a subset of the 28 study landscapes located within the region’s largest national park; Murray Sunset National Park (6630 km2; Fig. 1). Ten sites were established within each of 10 of the original study landscapes. Sites were distributed to incorporate a range of fire-age classes (range = 7–105 years), as well as capturing geographic and topographic variation. We established an additional landscape, containing 12 sites, following an experimental burn during the study (fire age = 0 years), resulting in 11 landscapes containing 112 sites. We omitted four sites to comply with ethics permits due to their close proximity to active nesting sites of the endangered malleefowl (Leipoa ocellata). This resulted in a total of 108 sites being surveyed. All sites were a minimum of 200 m apart and typically >100 m from the edge of a fire-age class.

Predator surveys

Landscape-scale data.

We surveyed large mammalian predators using track surveys from three sand-plots within each study landscape (n = 84 in total). Each sand-plot was a 100 m×2 m area smoothed out by dragging a weight along an unsealed vehicle track. The locations of sand-plots within landscapes were chosen to incorporate variation in the topography (dunes and swales) within each landscape. Sand-plots were typically >500 m apart. We checked each sand-plot for tracks once per day by walking along the transect and identifying tracks to species level for three consecutive days in spring (October–November 2007), and again in summer (January–March 2008), resulting in six survey nights for each sand-plot and thus 18 survey nights per landscape. Following checking, sand-plots were smoothed over in preparation for the following day. If the sand-plot was heavily disrupted on one day (due to weather or vehicle disturbance), it was surveyed for an additional day.

Site-scale data.

We used camera traps (Passive ScoutGuard 550; ScoutGuard IR Cameras, Australia) to survey for mammalian predators at the site-scale during April–July 2012. We installed one camera per site and deployed each for a minimum of 15 nights. We attached cameras to a post at a height of 0.5 m and positioned them facing southward. Vegetation was removed within the immediate area of the camera to reduce false triggering. A 15 second video was taken each time the camera motion sensor was triggered. To attract predators to the front of the camera from the local vicinity, we placed a scent lure of tuna oil soaked into chemical wadding inside a bait holder made from PVC piping with steel mesh at one end. We positioned the lure 3 m from the base of the camera post, and secured it to the ground with a peg.

Predictor variables

Landscape-scale data.

Six predictor variables were chosen to represent the properties of the study landscapes (Table 1). Three of these variables represent the fire history of the landscape: (1) the extent of recently burnt vegetation in the landscape (<10 years since fire; ‘recently burned’); (2) the extent of long unburnt vegetation (unburned since 1972; ‘long unburned’); and (3) the diversity of fire-ages within a landscape (‘fire diversity’). Fire diversity was calculated as the Shannon-Wiener diversity index of the proportional cover of fire age-classes within each landscape.

thumbnail
Table 1. Predictor variables included in models using the landscape-scale and the site-scale datasets.

https://doi.org/10.1371/journal.pone.0107862.t001

Three predictor variables were chosen to describe properties of the study landscapes other than fire history. We used a measure of mean solar radiation (‘solar radiation’) as a surrogate for aridity across the region. The solar radiation variable represents the total amount of solar energy falling on a horizontal space per day (MJ/m2). We derived these values from a gridded data set (5 km resolution) extending over 18 years (1990–2008; Australian Bureau of Meteorology http://www.bom.gov.au, 2009). Solar radiation was the mean of the 18 yearly averages of the grids that overlaid each landscape. Solar radiation is negatively correlated with annual rainfall and positively correlated with temperature. We used the proportional extent of Triodia mallee vegetation (‘Triodia Mallee’) within the landscape to capture differences in vegetation types. The extent of mallee vegetation in the study area was mapped in previous work (see [35]). Finally, we used the distance from the centre point of each landscape to the closest area of contiguous non-mallee vegetation (‘distance to agricultural land’), to capture the context of landscapes with respect to landscape modification. The area surrounding each reserve is comprised almost entirely of grazing land and grain crops. We calculated distance to agricultural land using ArcGIS [39].

Site-scale data.

Eight predictor variables were chosen at the site level (Table 1). The fire history of sites was represented by the time since the last fire (‘time since fire’; range: 0–105 years). This was determined using two methods. Recent fire history (since 1972) was calculated using the fire history maps (see [36]). Fire-ages for sites burnt prior to the availability of satellite imagery (i.e. before 1972) were estimated using regression models of the relationship between stem diameter and tree age, and then using stem diameter to estimate the age of trees in areas where fire history was unknown (see [37] for detailed methods). This extended the time since fire axis from 0–32 years to 0–105 years.

Vegetation type was considered as a categorical variable with two levels: Triodia Mallee or Chenopod Mallee (‘vegetation type’). We again considered the effects of landscape modification by including the distance of sites to both the border of the National Park (1.72–21.28 km; ‘distance to edge’) and dirt roads (range: 28–1044 m; ‘distance to road’). We used park boundary as a proxy for an edge habitat because the park forms abrupt boundaries with cleared agricultural land and other non-mallee vegetation. We calculated distance variables using ArcGIS [39]. Aridity (solar radiation) was not considered at this scale as the data were collected from a single reserve.

Four additional predictor variables were included to describe vegetation structure at the sites. We established vegetation transects in representative areas 15 m from each camera location. We recorded substrate type and vegetation structure at 1 m intervals along a 50 m transect using a 2 m structure pole (2 cm diameter) held vertically above the ground. The four variables considered in the analysis represent the cover of open, bare ground (‘bare ground cover’), spinifex (‘Triodia cover’), eucalypt shrubs (defined as Eucalypt trees <3 m in height ‘eucalypt cover’) and non-eucalypt shrubs (‘shrub cover’). Bare ground was included because it gives an approximation of the ‘openness’ of the vegetation at the ground level. The cover of spinifex, eucalypt shrubs, and non-eucalypt shrubs were included as they form the majority of the ground and understorey structural complexity, and are known to drive fauna in the region [34], [40].

Response variables

For both datasets, the response variable was the ‘reporting rate’ of foxes. At the landscape-scale, we defined reporting rate as the number of nights that a fox was recorded as ‘present’ and ‘absent’, respectively, at a sand pad over the 18 nights of sampling per landscape (i.e. three sand-plots surveyed for six nights in each landscape). Likewise, at the site-scale, reporting rate is the number of nights that foxes were and were not detected at the site, respectively, over the course of sampling (i.e. 15 nights).

Statistical analysis

We used generalised linear mixed models (GLMMs) with the Laplace approximation [41] to examine the relationship between response and predictor variables at both landscape and site-scales. In landscape-scale models, we included ‘reserve’ as a random effect to account for spatial clustering of landscapes in conservation reserves (Fig. 1). Similarly, in the site-scale models, we included ‘landscape’ as a random effect to account for potential spatial correlation due to the clustering of sites into landscapes. Because we were studying the reporting rate of red foxes, a proportion, we modelled the response variable (at both scales) using a binomial distribution of errors and a logit link function.

For the landscape-scale dataset, we developed a set of candidate models that included all combinations of the six landscape-scale predictor variables. At the site-level, we developed two separate sets of models. As fire affects the variables used to describe vegetation structure (e.g. Triodia cover, bare ground cover; [8]), including both fire and vegetation structure variables in the same model could result in unreliable parameter estimates due to colinearity between predictor variables [42]. Thus, one model set (model set 1) included time since fire, vegetation type, distance to edge and distance to road, and a second model set (model set 2) included the vegetation structure variables (bare ground cover, Triodia cover, eucalypt cover and shrub cover). All combinations of predictors within the two sets of models were considered, meaning all eight site-level variables were in the same number of models overall. All variables included within a model set had low levels of colinearity (i.e. r <0.5). We tested both datasets for overdisperson using Pearson’s residuals [43], and found no evidence of overdispersion.

We compared each set of candidate models using Akaike’s Information Criterion corrected for small sample sizes (AICc; [44]). To compare the level of support for each model relative to the most parsimonious model, we calculated the difference (Δi) between the AICc value of the best model (lowest AICc value) and the AICc value of each candidate model [44]. We considered models with Δi<2 to have substantial support [44]. We also calculated the Akaike weight (wi) for each model. By summing these weights to calculate predictor weights (∑wi) for each variable, we were able to explore the influence of individual predictor variables at both the landscape and site level.

When there was no clear ‘best model’ (i.e. the most parsimonious model was not strongly weighted [wi<0.9]), we used model averaging to determine the direction and magnitude of the effect of each predictor variable [44]. We considered a variable as important when the associated 95% confidence interval of the averaged estimate did not overlap with zero. We performed all statistical analyses in R version 2.15.1 [45] using the lme4 package [41] and the MuMIn package [46].

Ethics statement

The landscape-scale data were collected with approval from animal ethics committees at La Trobe University (approval number AEC06/07[L]V2) and Deakin University (approval number A41/2006), and permits from the Department of Sustainability and Environment, Victoria (permit 10003791), the Department of Environment and Heritage, South Australia (permit 13/2006), and the National Parks and Wildlife Service, NSW (license number S12030). The site-scale data were collected in accordance with the regulations of the Deakin University Animal Ethics Committee (approval number B10-2012) and in accordance with Department of Sustainability and Environment, Victoria (approval number 10006279).

Results

At the landscape-scale, we recorded fox tracks in 24 of 28 (86%) study landscapes. We detected foxes on 3.32±0.49 (mean ± standard error) of 18 nights per landscape over the total sampling period. Other large-bodied, mammalian predators were uncommon: we detected cats at only 7 of 28 landscapes (25%). At the site-scale, we observed foxes at 62 of 102 (61%) sites (six cameras failed to reach the full 15 day survey period due to fault and were excluded from further analysis i.e. n = 102) and found the species to be widely distributed across the study area. We did not detect any cats at the site-scale over the 15 night sampling period.

At the landscape-scale, all models were a poor fit for the data and explained <6.5% of the variation in the data (% deviance explained). At the site-scale, all models explained <3.5% of the variation in the data. For both datasets, model selection indicated there was a similar level of support for several models (Δi<2; Table 2), including the intercept-only model (i.e. only an intercept terms, no predictor variables), which received substantial support at both scales. As no single model was supported as being clearly best (i.e. wi>0.9; Table 2), we employed multi-model inference using model averaging to estimate the size, direction and uncertainty of parameter effects for fox explaining reporting rate in both datasets.

thumbnail
Table 2. Model selection results for red fox reporting rate for landscape-scale and sits-scale datasets.

https://doi.org/10.1371/journal.pone.0107862.t002

The model-averaged coefficients for each predictor variable, in both datasets, were small and uncertain. The 95% confidence intervals of all predictor variables overlapped with zero (Fig. 3). The ∑wi for all predictor variables was low: <0.5 and <0.6 for the landscape- and site-scale datasets respectively.

thumbnail
Figure 3. Model-averaged regression coefficients and 95% confidence intervals of models describing the reporting rate of foxes at both the landscape-scale (a) and site-scale (b and c).

https://doi.org/10.1371/journal.pone.0107862.g003

Graphical exploration of the data further highlights that fox activity was not strongly linked to key predictor variables (Fig. 4). In summary, the data shows that neither fire, nor any other predictor variable measured, affected the reporting rate of foxes at either the landscape- or site-scale.

thumbnail
Figure 4. Relationships between the reporting rate of the red fox and the properties of fire mosaics.

Circles are raw data points.

https://doi.org/10.1371/journal.pone.0107862.g004

Discussion

Introduced mesopredators and fire are two processes that shape ecosystems around the world [4], [7]. Here, we have shown that a widespread and ecologically devastating mesopredator, the red fox [5], is largely unaffected by fire and is an extreme habitat generalist in semi-arid Australia. This result was confirmed using two large, complementary datasets, collected at different times and characterised by differing spatial scales and sampling strategies.

Fire and the red fox

Our findings show that fire does not exert a strong influence on the distribution of the red fox in semi-arid mallee ecosystems. Despite conducting two intensive natural experiments across a broad geographic region, we did not detect a relationship between the reporting rate of foxes and fire history at either the landscape- or site-scale. At the landscape-scale, the red fox was recorded equally often in landscapes dominated by recently burned or long unburned vegetation, and in landscapes with a single fire age-class as those with a diversity of fire ages. At the site-scale, the red fox has a similar reporting rate in recently burned sites as in sites unburned for over a century. The post-fire preferences of the red fox are thus extremely broad, both spatially and temporally (also see [12], [47]).

Fire causes significant changes to vegetation structure over century-long time frames in mallee ecosystems [8]. In doing so, fire affects the distribution of a large range of fauna species [34]. Indeed, work conducted within the same study landscapes has shown the large and long-term effects fire has on birds, reptiles, and small mammals [28][30]. The lack of a response to fire by foxes is therefore not typical of native fauna in the region. It also suggests that foxes are not restricted to areas with particular soil or vegetation attributes for denning. This is consistent with foxes not being affected by any of the vegetation attributes measured (e.g. Triodia cover, shrub cover etc.).

A related way that fire could influence foxes is by altering the distribution of prey resources. As mentioned above, the distribution of many prey species are significantly affected by fire in the study region (e.g. birds, mammals, reptiles). Thus, foxes occupy a range of post-fire ages despite the strong influence of fire on the type and abundance of prey available. Red foxes have a broad and generalist diet [48], being able to consume a wide range of prey including both vertebrates and invertebrates, and even vegetation [19], [20]. Furthermore, foxes are capable of prey switching to capitalize on the most abundant prey source available [18], [49], thereby reducing their reliance on any particular prey item. This flexibility in their diet is likely to be a key component of their life history that allows them to occur within such a broad range of post-fire conditions.

One objection to our findings at the site-scale may be that the local site is not a relevant spatial scale to characterize the effects of fire, as foxes are a relatively large and mobile species. Given the large estimated home ranges of foxes in other parts of arid Australia (e.g. 8–33 km2; [31]), foxes may select broader areas (i.e. kms2) that capture their resource requirements across entire landscapes, and this might include a large area of a particular fire-age, or multiple fire ages. Such use of multiple habitat types by foxes has been demonstrated in other systems [23], [24]. Our landscape-scale study characterized land mosaics at a large scale relevant to the home range of foxes (12.6 km2), and still failed to detect any relationship between fox activity and fire history. Therefore, our results suggest that the lack of relationships between fox reporting rate and fire history does not stem from spatial scaling issues. Instead, foxes are resilient towards the effects of fire at multiple temporal and spatial scales.

Climate and distance to modified land

In addition to fire, we examined other variables that could influence the distribution of the red fox. Here, we again found red foxes to be flexible to a broad range of ecological conditions. Foxes displayed no response to an aridity gradient across the study region. This lack of response to aridity is unsurprising, as the geographic range of the red fox spans the northern hemisphere and much of Australia, suggesting the species is capable of coping with a range of climatic conditions.

Despite foxes occupying a broad climatic niche in space, fluctuations in populations do occur in response to extreme weather events. For example, fox populations in arid areas rise rapidly following high rainfall events, in response to increased prey availability [50]. Our site-scale study was carried out during a year of record high rainfall (Australian Bureau of Meteorology, Ouyen Station). Considered in isolation, this may suggest that the wide distribution of the fox was partly due to a productivity-related increase in food resources (predominantly populations of native and introduced rodents; [40]). However, the landscape-scale data were collected near the end of a severe, decade-long drought. Foxes were widely distributed across the region despite the drought. This indicates that, in semi-arid Australia, foxes can be widespread during a broad range of climatic conditions and despite fluctuations in their prey populations which accompany climactic extremes [40].

Some studies have found that foxes are positively associated with edges between fragments and modified land (e.g. agricultural land) [22], [32]. Our results indicate that foxes do not show a preference for edge habitats in mallee ecosystems, despite our sites and landscapes capturing a broad gradient of distances to agricultural land, from <2 km to >30 km. Edge habitats may be more important for foxes in highly fragmented landscapes, where they occur with small remnant patches of wooded vegetation which provide the only available cover [22]. While the mallee region has been subject to large amounts of land clearing, there are still relatively large intact areas of native vegetation. Edges may be less important in this region because the interior mallee vegetation provides sufficient shelter and prey. Nevertheless, it is also possible that edge effects occur closer to the agricultural boundary than we sampled (i.e. <2 km).

The use of roads and tracks by foxes is also well documented [51], [52]. Foxes have been found to be more abundant along roadsides [33]. In the mallee system, however, we found similar reporting rates at varying distances (28–1044 m) from roads, indicating foxes use areas well away from roads equally as often as sites close to roads. One hypothesis for the use of roads by foxes is that they provide ‘runways’ which facilitate movement and allow access to foraging areas that would be otherwise difficult to reach [51], [52]. In contrast to environments with a dense understory, mallee vegetation is relatively open, and is unlikely to limit the movement of foxes to roadsides. This may explain the lack of preference for sites near roads in the current study.

Implications

Fire is used as a conservation tool in Australia and around the world [25]. This study suggests it is unlikely that any particular approach to fire management will alter the reporting rate of the red fox in the semi-arid mallee systems of Australia. However, the presence of foxes in recently burned sites and landscapes is a concern. Predation by invasive mesopredators has been hypothesized as a cause of low post-fire survival in reptiles [14] and mammals [12], [13], due to the reduced cover available in burnt habitats. Although we found no effect of fire history on red fox occurrence, it is possible that predation pressure differs across fire ages due to increased predation risk in recently burned areas. Thus, assessing predation pressure directly across a range of post-fire ages is an important area for further research.

The loss of apex predators can cause smaller predators to increase in abundance, expand their range, and change their temporal activity; this is known as ‘mesopredator release’ [2], [6]. Red foxes have been shown to select particular habitats which may allow them to avoid dominant predators (e.g. coyotes; [53]). As such, one further explanation for the lack of obvious habitat selection by foxes in this system may be the lack of regulating predators. In other Australian systems, the presence of the dingo, Australia’s largest terrestrial apex predator, has been shown to affect fox distributions [3]. Dingoes are largely extinct from the study area but were once common, and as such there is no direct regulation of the abundance or distribution of foxes via biotic interactions. Thus, one potential way to control red foxes in mallee communities is by reinstating dingoes as the apex predator. As this is likely to be a controversial idea owing to the proximity of mallee vegetation to agricultural land and livestock, trialing reintroductions in a controlled and experimental way would be an important first step towards a proof of concept, and a potential solution to this complex conservation issue.

Supporting Information

Table S1.

Landscape-scale data on the distribution of red foxes in semi-arid land mosaics.

https://doi.org/10.1371/journal.pone.0107862.s001

(DOCX)

Acknowledgments

Thanks to all members of the Mallee Fire and Biodiversity Team, particularly Andrew Bennett, Mike Clarke, Lisa Farnsworth and Lauren Brown, and to the many volunteers, agency staff, and land owners who assisted with logistics and field work.

Author Contributions

Conceived and designed the experiments: DGN CJP EGR LTK. Performed the experiments: DGN CJP LTK. Analyzed the data: DGN CJP. Contributed to the writing of the manuscript: CJP DGN LTK EGR.

References

  1. 1. Estes JA, Terborgh J, Brashares JS, Power ME, Berger J, et al. (2011) Trophic downgrading of planet Earth. Science 333: 301–306.
  2. 2. Ritchie EG, Johnson CN (2009) Predator interactions, mesopredator release and biodiversity conservation. Ecology Letters 12: 982–998.
  3. 3. Letnic M, Ritchie EG, Dickman CR (2012) Top predators as biodiversity regulators: the dingo Canis lupus dingo as a case study. Biological Reviews 87: 390–413.
  4. 4. Salo P, Korpimäki E, Banks PB, Nordström M, Dickman CR (2007) Alien predators are more dangerous than native predators to prey populations. Proceedings of the Royal Society B: Biological Sciences 274: 1237–1243.
  5. 5. Johnson C (2006) Australia’s Mammal Extinctions: a 50000 year history. New York: Cambridge University Press.
  6. 6. Crooks KR, Soule ME (1999) Mesopredator release and avifaunal extinctions in a fragmented system. Nature 400: 563–566.
  7. 7. Bowman DMJS, Balch JK, Artaxo P, Bond WJ, Carlson JM, et al. (2009) Fire in the Earth system. Science 324: 481–484.
  8. 8. Haslem A, Kelly LT, Nimmo DG, Watson SJ, Kenny SA, et al. (2011) Habitat or fuel? Implications of long-term, post-fire dynamics for the development of key resources for fauna and fire. Journal of Applied Ecology 48: 247–256.
  9. 9. Smit IPJ, Asner GP, Govender N, Kennedy-Bowdoin T, Knapp DE, et al. (2010) Effects of fire on woody vegetation structure in African savanna. Ecological Applications 20: 1865–1875.
  10. 10. Letnic M, Tamayo B, Dickman CR (2005) The responses of mammals to La Nina (El Nino Southern Oscillation) - associated rainfall, predation and wildfire in central Australia. Journal of Mammalogy 86: 689–703.
  11. 11. Driscoll DA, Lindenmayer DB, Bennett AF, Bode M, Bradstock RA, et al. (2010) Fire management for biodiversity conservation: Key research questions and our capacity to answer them. Biological Conservation 143: 1928–1939.
  12. 12. Arthur AD, Catling PC, Reid A (2012) Relative influence of habitat structure, species interactions and rainfall on the post-fire population dynamics of ground-dwelling vertebrates. Austral Ecology 37: 958–970.
  13. 13. Sutherland EF, Dickman CR (1999) Mechanisms of recovery after fire by rodents in the Australian environment: a review. Wildlife Research 26: 405–419.
  14. 14. Smith AL, Bull CM, Driscoll DA (2012) Post-fire succession affects abundance and survival but not detectability in a knob-tailed gecko. Biological Conservation 145: 139–147.
  15. 15. Woinarski JCZ, Legge S, Fitzsimons JA, Traill BJ, Burbidge AA, et al. (2011) The disappearing mammal fauna of northern Australia: context, cause, and response. Conservation Letters 4: 192–201.
  16. 16. Kinnear JE, Onus ML, Bromilow RN (1988) Fox control and rock-wallaby population dynamics. Wildlife Research 15: 435–450.
  17. 17. Risbey DA, Calver MC, Short J, Bradley JS, Wright IW (2000) The impact of cats and foxes on the small vertebrate fauna of Heirisson Prong, Western Australia. II. A field experiment. Wildlife Research 27: 223–235.
  18. 18. Catling PC (1988) Similarities and contrasts in the diets of foxes, Vulpes vulpes, and cats, Felis catus, relative to fluctuating prey populations and drought. Wildlife Research 15: 307.
  19. 19. Glen AS, Fay AR, Dickman CR (2006) Diets of sympatric red foxes Vulpes vulpes and wild dogs Canis lupus in the Northern Rivers Region, New South Wales. Australian Mammalogy 28: 101–104.
  20. 20. Risbey DA, Calver MC, Short J (1999) The impact of cats and foxes on the small vertebrate fauna of Heirisson Prong, Western Australia. I. Exploring potential impact using diet analysis. Wildlife Research 26: 621–630.
  21. 21. Letnic M, Dickman CR (2005) The responses of small mammals to patches regenerating after fire and rainfall in the Simpson Desert, central Australia. Austral Ecology 30: 24–39.
  22. 22. Graham CA, Maron M, McAlpine CA (2012) Influence of landscape structure on invasive predators: feral cats and red foxes in the brigalow landscapes, Queensland, Australia. Wildlife Research 39: 661–676.
  23. 23. Lucherini M, Lovari S, Crema G (1995) Habitat use and ranging behaviour of the red fox (Vulpes vulpes) in a Mediterranean rural area: is shelter availability a key factor? Journal of Zoology 237: 577–591.
  24. 24. Meek PD, Saunders G (2000) Home range and movement of foxes (Vulpes vulpes) in coastal New South Wales, Australia. Wildlife Research 27: 663–668.
  25. 25. Parr CL, Andersen AN (2006) Patch mosaic burning for biodiversity conservation: a critique of the pyrodiversity paradigm. Conservation Biology 20: 1610–1619.
  26. 26. Southgate R, Paltridge R, Masters P, Ostendorf B (2007) Modelling introduced predator and herbivore distribution in the Tanami Desert, Australia. Journal of Arid Environments 68: 438–464.
  27. 27. Gosper CR, Yates CJ, Prober SM (2013) Floristic diversity in fire-sensitive eucalypt woodlands shows a ‘U’-shaped relationship with time since fire. Journal of Applied Ecology 50: 1187–1196.
  28. 28. Watson SJ, Taylor RS, Nimmo DG, Kelly LT, Haslem A, et al. (2012) Effects of time since fire on birds: How informative are generalized fire response curves for conservation management? Ecological Applications 22: 685–696.
  29. 29. Nimmo DG, Kelly LT, Spence-Bailey LM, Watson SJ, Taylor RS, et al. (2013) Fire mosaics and reptile conservation in a fire-prone region. Conservation Biology 27: 345–353.
  30. 30. Kelly LT, Nimmo DG, Spence-Bailey LM, Taylor RS, Watson SJ, et al. (2012) Managing fire mosaics for small mammal conservation: a landscape perspective. Journal of Applied Ecology 49: 412–421.
  31. 31. Moseby KE, Stott J, Crisp H (2009) Movement patterns of feral predators in an arid environment – implications for control through poison baiting. Wildlife Research 36: 422–435.
  32. 32. Catling P, Burt R (1995) Why are red foxes absent from some eucalypt forests in eastern New South Wales? Wildlife Research 22: 535–545.
  33. 33. Towerton AL, Penman TD, Kavanagh RP, Dickman CR (2011) Detecting pest and prey responses to fox control across the landscape using remote cameras. Wildlife Research 38: 208–220.
  34. 34. Bradstock RA, Cohn JS (2002) Fire regimes and biodiversity in semi-arid mallee ecosystems. In: R. A Bradstock, J. E Williams and M. A Gill, editors. Flammable Australia: The fire regimes and biodiversity of a continent. Cambridge: Cambridge University Press. 238–258.
  35. 35. Haslem A, Callister KE, Avitabile SC, Griffioen PA, Kelly LT, et al. (2010) A framework for mapping vegetation over broad spatial extents: A technique to aid land management across jurisdictional boundaries. Landscape and Urban Planning 97: 296–305.
  36. 36. Avitabile SC, Callister KE, Kelly LT, Haslem A, Fraser L, et al. (2013) Systematic fire mapping is critical for fire ecology, planning and management: A case study in the semi-arid Murray Mallee, south-eastern Australia. Landscape and Urban Planning 117: 81–91.
  37. 37. Clarke MF, Avitabile SC, Brown L, Callister KE, Haslem A, et al. (2010) Ageing mallee eucalypt vegetation after fire: insights for successional trajectories in semi-arid mallee ecosystems. Australian Journal of Botany 58: 363–372.
  38. 38. ITT (2005) ENVI. Version 4.2. Boulder, Colorado: ITT Industries.
  39. 39. Environmental Systems Research Institute (2007) Arc View. Version 9.2. Redlands, California: ESRI.
  40. 40. Kelly LT, Dayman R, Nimmo DG, Clarke MF, Bennett AF (2013) Spatial and temporal drivers of small mammal distributions in a semi-arid environment: The role of rainfall, vegetation and life-history. Austral Ecology 38: 786–797.
  41. 41. Bates D, Maechler M, Bolker B (2012) lme4: Linear mixed-effects models using S4 classes. R package (Version 0.999999-0). Available: http://CRAN.R-project.org/package=lme4.
  42. 42. Quinn GP, Keough MJ (2002) Experimental design and data analysis for biologists. New York: Cambridge University Press.
  43. 43. Zuur AF, Ieno EN, Walker N, Saveliev AA, Smith GM (2009) Mixed effects models and extensions in ecology with R. New York: Springer.
  44. 44. Burnham KP, Anderson DR (2002) Model selection and multi-model inference. New York: Springer.
  45. 45. R Development Core Team (2012) R: A language and environment for statistical computing. Vienna, Austria. Available: http://www.R-project.org/: R Foundation for Statistical Computing.
  46. 46. Bartoń K (2012) MuMIn: Multi-model inference. R package (Version 1.7.11). Available: http://CRAN.Rproject.org/package=MuMIn.
  47. 47. Catling PC, Coops N, Burt RJ (2001) The distribution and abundance of ground-dwelling mammals in relation to time since wildfire and vegetation structure in south-eastern Australia. Wildlife Research 28: 555–565.
  48. 48. White JG, Gubiani R, Smallman N, Snell K, Morton A (2006) Home range, habitat selection and diet of foxes (Vulpes vulpes) in a semi-urban riparian environment. Wildlife Research 33: 175–180.
  49. 49. Leckie FM, Thirgood SJ, May R, Redpath SM (1998) Variation in the diet of red foxes on Scottish moorland in relation to prey abundance. Ecography 21: 599–604.
  50. 50. Pavey CR, Eldridge SR, Heywood M (2008) Population dynamics and prey selection of native and introduced predators during a rodent outbreak in arid Australia. Journal of Mammalogy 89: 674–683.
  51. 51. Carter A, Luck GW, McDonald SP (2012) Ecology of the red fox (Vulpes vulpes) in an agricultural landscape. 2. Home range and movements. Australian Mammalogy 34: 175–187.
  52. 52. Frey SN, Conover MR (2006) Habitat use by meso-predators in a corridor environment. The Journal of Wildlife Management 70: 1111–1118.
  53. 53. Gosselink TE, Deelen TRV, Warner RE, Joselyn MG (2003) Temporal habitat partitioning and spatial use of coyotes and red foxes in east-central Illinois. The Journal of Wildlife Management 67: 90–103.