Next Article in Journal
The Permeability of Natural versus Anthropogenic Forest Edges Modulates the Abundance of Ground Beetles of Different Dispersal Power and Habitat Affinity
Next Article in Special Issue
Facilitation or Competition? Effects of Lions on Brown Hyaenas and Leopards
Previous Article in Journal
A National Survey of Managed Honey Bee Colony Winter Losses (Apis mellifera) in China (2013–2017)
Previous Article in Special Issue
Exotic Prey Facilitate Coexistence between Pumas and Culpeo Foxes in the Andes of Central Chile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inferring Species Interactions from Long-Term Monitoring Programs: Carnivores in a Protected Area from Southern Patagonia

by
Francisco Díaz-Ruiz
1,*,
Alejandro Rodríguez
2,
Diego Procopio
3,
Sonia Zapata
3,
Juan Ignacio Zanón-Martínez
4,5 and
Alejandro Travaini
3
1
Biogeography, Diversity and Conservation Research Team, Department of Animal Biology, Faculty of Sciences, University of Málaga, 29071 Málaga, Spain
2
Department of Conservation Biology, Estación Biológica de Doñana, CSIC, Américo Vespucio 26, 41092 Sevilla, Spain
3
Centro de Investigación de Puerto Deseado, Universidad Nacional de la Patagonia Austral, CONICET, Avenida Prefectura Naval s/n, 9050 Puerto Deseado, Santa Cruz, Argentina
4
Centro para el Estudio y Conservación de las Aves Rapaces en Argentina (CECARA), Facultad de Ciencias Exactas y Naturales, Universidad Nacional de La Pampa, Avda. Uruguay 151, 6300 Santa Rosa, La Pampa, Argentina
5
Instituto de las Ciencias de la Tierra y Ambientales de La Pampa (INCITAP), Consejo Nacional de Investigaciones Científicas y Técnicas, Avda. Uruguay 151, 6300 Santa Rosa, La Pampa, Argentina
*
Author to whom correspondence should be addressed.
Diversity 2020, 12(9), 319; https://doi.org/10.3390/d12090319
Submission received: 11 July 2020 / Revised: 14 August 2020 / Accepted: 18 August 2020 / Published: 21 August 2020

Abstract

:
Protected areas recently created in Argentina often include previously degraded lands, such as sheep ranches in the Patagonian deserts. We show the results of a 14-year monitoring program of three formerly persecuted carnivores, the culpeo fox (Lycalopex culpaeus), the South American grey fox (Lycalopex griseus) and the puma (Puma concolor), in two abandoned sheep ranches that were incorporated into a Patagonian national park approximately 25 years ago. The culpeo fox population underwent an average annual decline of 10–23%, whereas the grey fox and puma populations increased at an average annual rate of 7% and 19%, respectively. The grey fox’s increasing trends were strongly correlated with the decline of the culpeo fox, whereas the correlations between the fox and puma trends were weaker. Culpeo fox decline was stronger in the ranch where sheep and predator controls had been removed earlier. These relationships between species trends support the competitive release hypothesis, assuming that puma competition with the culpeo fox for trophic resources is stronger than competition with the grey fox, and that the puma can exclude culpeo foxes through interference. Species trends suggest a competitive hierarchy between fox species, with grey fox being the inferior competitor. However, mechanisms other than competition could not be discounted. Our study illustrates how long-term monitoring of interacting species allows a better understanding of ecological processes and wildlife ecology.

Graphical Abstract

1. Introduction

Few areas across the planet remain where natural systems have not been damaged, directly or indirectly, by human activities, which in many cases pose a threat for biodiversity conservation [1,2,3]. The establishment of protected areas is a common conservation strategy to reduce the intensity of detrimental human activity or to regulate land use to better conserve biodiversity and restore altered ecological processes [4,5,6]. Over the last century, the area covered by nature reserves has increased globally, especially in developing countries harbouring high biodiversity [6], and currently makes up approximately 15% of Earth’s land area [7]. Once a protected area is created, it is assumed that critical ecological processes, resources, or species of conservation concern will be preserved or restored, because major conservation threats (e.g., deforestation, habitat fragmentation or poaching) are expected to be largely mitigated within their limits. However, the benefits of protected areas for biodiversity conservation are often unknown, mainly due to the lack of long-term monitoring programs allowing an objective evaluation of their effectiveness [8].
Long-term monitoring provides reliable data on the status and trends of key natural resources and/or species populations targeted for conservation or harvesting [9]. When a monitoring program is well designed, for example when it is based on power analysis, monitoring may detect in a timely fashion significant deviation from the desired conservation conditions [10]. This allows the application of appropriate management measures to reverse negative trends and achieve conservation goals. Further, an adaptive management framework should be adopted when there is a lack of knowledge about the best management option [10,11]. Long-term monitoring programs are often viewed as a tool exclusively applied to management. However, monitoring data could also be used to address research questions [12,13]. Thus, many long-term data sets collected for management and conservation purposes may have high scientific value, contributing a spatial and temporal depth that is rarely seen in research projects, and allowing the proposal and, sometimes, testing of ecological hypotheses [13].
Nowadays, the creation of new protected areas is not always done on pristine or untransformed natural areas. Rather, protection is often applied to large areas degraded by human activity [14]. In these cases, long-term monitoring programs may provide valuable information about the ecological processes involved in the ecosystem restoration that would be expected as a result of land protection [14,15]. During the last two decades, several protected areas have been created in different landscapes of Patagonia, both in Chile and Argentina. Some of these new protected areas are good examples of land use conversion of formerly degraded ecosystems as a consequence of unsustainable livestock farming, e.g., Monte León National Park in Argentina [16] and Patagonia National Park in Chile [17].
In the Argentine Patagonian steppes, sheep (Ovis aries) husbandry has been the main land use and economic activity since European settlement at the end of the 19th century (1890) [18]. The pressure exerted on the Patagonian ecosystems by this activity has been high throughout the 20th century, affecting wildlife communities both directly and indirectly. Direct persecution of species perceived as a threat to livestock and the introduction of exotic species have modified the composition and structure of mammal assemblages and the strength of interspecific interactions. For example, large native herbivores such as the guanaco, Lama guanicoe, or the lesser rhea, Pterocnemia pennata, underwent significant population declines due to a combination of intense hunting, habitat degradation and direct competition with mammalian herbivores introduced into Patagonia by the first European settlers [19,20], in particular sheep and the European hare, Lepus europaeus [21,22,23]. In addition to the competition with native species for grasses, the European hare modified food web interactions, becoming a very important resource for many Patagonian native predators before the decrease of native prey [22,24].
Among native predators, mammalian carnivores have been heavily persecuted in sheep ranches throughout Patagonia, either as a measure to control livestock predation [25] or as a way for rural people to get extra economic income [26]. Pumas, Puma concolor, are the largest carnivore species in Patagonian ecosystems [27]. Perceived by many ranchers as one of the main threats to livestock [28], lethal control has been widely implemented across Patagonia, where pumas were extirpated from many localities and regions during the second half of the 20th century [28,29]. Smaller carnivores, including the culpeo fox, Lycalopex culpaeus, and the South American grey fox, Lycalopex griseus (henceforth “grey fox”) have also been heavily persecuted throughout Patagonia, where they are simultaneously considered a threat to livestock and a valuable species for their pelts [30]. As culpeo foxes regularly prey on lambs and young sheep [22,31], intensive culpeo fox control has been a common activity in sheep ranches [25], reaching annual reductions of over 75% in some local populations [32]. Although the grey fox does not pose a risk to livestock, it has also been culled, either for the value of its fur or indirectly due to the extensive use of nonselective predator control methods [25]. Patagonian sheep stock peaked in 1952 at 25 million head and decreased continuously thereafter, as a consequence of overgrazing effects, desertification and wool price devaluation [18,33]. During the second half of the 20th century, livestock activity was completely abandoned in large regions of Patagonia [34], which has allowed the recovery of some native species such as the puma and the guanaco [23,35]. However, predator persecution of foxes and pumas has persisted in some areas [29,36]. Accordingly, the persecution exerted on carnivores over time might have altered the competitive interactions between these species, as well as their relative abundances.
Changes in the relative abundance of sympatric carnivores can have far-reaching ecological consequences, including the precipitation of trophic cascades and species declines [37,38]. Carnivores play an important role in preserving ecological function and dynamics of ecosystems, driving top-down effects such as controlling populations of native or introduced species [39,40]. Larger predators may regulate populations of smaller predators through exploitative competition upon shared resources, or through interference competition by harassment and intraguild predation [41]. One consequence of losing larger predators is an increase of mesopredator abundance [38,42,43]. As a rule, larger species tend to dominate competitive interactions over smaller species [44]. The competitive release hypothesis states that a reduction in the abundance of a dominant competitor triggers an increase in the numbers of subordinate competitors [45,46], whereas the recovery of the superior competitor would produce the opposite effect.
The puma plays a dominant role over both fox species, including high potential to kill them [47]. Going down in the competitive hierarchy, a dominant role of the larger culpeo fox (weighing ~8 kg) over the smaller grey fox (weighing ~2.6 kg) [48] would be expected. Although there is evidence of the recovery of the puma throughout many regions of Patagonia after the reduction of predator control associated with the abandonment of sheep ranches [35,36], it is unknown how puma resurgence has affected other carnivores such as the culpeo fox and the grey fox and their competitive interactions. The diets of these three species overlap. Pumas mainly prey on larger prey such as guanacos, but European hares are an important diet item [24,49]. Grey foxes mainly feed on small mammals and invertebrates, although European hares are also included in their diet [48,50]. Culpeo foxes feed mainly on European hares and small mammals [48,50].
Twenty years ago, our research team started a long-term monitoring programme of the populations of the two fox species present in Patagonia [30], whose main objective was to document population trends in order to guide conservation management [51]. This monitoring programme was carried out in a protected area which had recently incorporated abandoned sheep ranches. Adopting an adaptive monitoring framework [11], in 2002, we included the puma as a monitored species in order to better understand changes in fox populations. In this work, we present population trends of the three species and explore whether the observed trends support the predictions of dominance hierarchy and the competitor release effects. We also highlight the relevance of long-term monitoring programs in recently created protected areas for a better understanding of ecological processes triggered after the cessation of human activities.

Study Area

The study was conducted in the Bosques Petrificados de Jaramillo National Park (hereafter denoted as BPJNP) located in the central north of the Santa Cruz province, Argentine Patagonia (47.66° S 65.89° W; Figure 1). Currently, the BPJNP extends over 78,543 ha of Patagonian shrub-steppe, an ecoregion of special interest for global conservation [52]. In 1997, two neighbouring sheep ranches were incorporated into the national park, increasing its size by about 30,000 ha. Prior to its incorporation, Cerro Horqueta ranch (CH) had ~1000 sheep, while El Cuadro ranch (EC) had ~500, but the livestock load for both fields together exceeded 5000 sheep in the 1950s. Livestock raising and predator persecution went on at least until 1999 in EC, while these activities ceased in CH around 1992. In both ranches, the vegetation is dominated by tussock grasses and low, dome-shaped spiny shrubs [53], whose cover ranges from less than 10% in the most arid sites to 60%, suggesting a marked degradation that could be attributed to many decades of overgrazing, as observed in other Patagonian steppes [54,55]. Summer temperatures average 17 °C, and winters are relatively mild with an average of ~12 freezing days/year. Annual rainfall ranges between 100 and 300 mm, and snow is rare.

2. Materials and Methods

2.1. Field Surveys

The monitoring program was designed and launched in 2000 [30,51]. To assess the distribution and relative abundance of foxes, we used bait stations, a noninvasive method to survey carnivores at large scales [56]. A bait station consisted of a 1-m-diameter circular surface of sifted, smoothed earth, with a 30 g bait buried at its centre [57]. Bait composition consisted of a mixture of minced meat, hydrogenous oil, cornstarch and the commercial lure Cat Passion (O´Gorman Enterprise Inc., Broadus, Montana, USA) [58].
We set 8 lines of 6 stations in each abandoned sheep ranch. Stations were set on alternate sides along secondary, unpaved roads. Separation of adjacent stations within each line was 0.5 km, and lines were at least 1 km apart. We assumed that the probability of any individual fox visiting more than one line per survey approached zero. Therefore, we considered each line as an independent sampling unit. Lines were surveyed twice per year, at six-month intervals, for 14 years (2000–2013) in autumn and 13 years (2000–2012) in spring. For each survey, we operated bait stations during three consecutive nights. Each morning, we recorded whether fox tracks were printed on the smoothed soil. Fox tracks were easily distinguishable between the species, since grey fox tracks are narrower and much shorter than culpeo fox tracks [58]. Tracks for which accurate species identification was not possible due to strong winds and/or poor printing were excluded from the analyses. In order to maximise accurate track identification, all surveys were conducted in dry weather.
Bait stations were designed to monitor fox populations, and they did not perform well in detecting puma trends in our study area (unpublished data). Therefore, between 2002 and 2009, we surveyed the puma population using fresh signs (scats, tracks and scrape marks) detected along foot transects as indices of relative abundance [59,60]. The study area was divided into 134 4-km2 square cells (72 in CH and 62 in EC; Figure 1). With the aid of a GPS (Garmin 3Plus; Lenexa, Kansas, USA), a mean distance (± SD) of 3.01 ± 0.71 km was erratically travelled in each cell, starting at one of the cell edges and passing through the cell’s centre. All fresh signs of pumas within ~2 m of each side of the transect were recorded. Transects were mainly performed by the same observer (95%; D. Procopio), who was occasionally supported by A. Travaini and J.I. Zanón. Logistical difficulties from financial constraints and a lack of human resources prevented us from sampling all cells each year. Despite this, on average, more than 50% of cells from each ranch were visited annually (mean% ± SE for CH: 58% ± 8%; for EC: 60% ± 7%). Surveys were carried out in autumn, except for the first two years, when they were done in spring. By pooling spring and autumn data, we assumed that season had a small effect on the index of puma abundance.

2.2. Analyses

Visitation indices denote the possible ways of converting field data into a quantitative indicator of species activity at a bait station. Different visitation indices are built using different combinations of daily records throughout the three-day operation period. Travaini et al. [30] used power analysis to examine the effect of choosing each of seven possible visitation indices on the ability of the monitoring program to detect changes in fox abundance. They found that the proportion of zeroes (no visits) had a large effect and that recording fox visits within 72 h after activation showed the lowest variation and improved the expected power to detect population trends of both fox species monitored simultaneously. We used this visitation index, defined as the number of bait stations visited by each fox species per line during any of the three consecutive nights [30]. Thus, our index of relative abundance could take integer values between zero and six visited stations per line.
Population trends were estimated using the R package rtrim [61], a reimplementation of the original TRIM software [62]. TRIM was specifically developed for the analysis of wildlife monitoring data and has been widely used in the study of population trends of several vertebrate taxa, e.g., [63,64]. TRIM allows dealing with common issues in long-term time series: missing counts, overdispersion and serial correlation. Missing counts are replaced with predicted counts estimated by log-linear Poisson regressions models, allowing for analysing data sets with up to 50% of missing counts [62,65]. In our data set, missing values were between 6% and 12% for foxes, and 41% for puma. TRIM includes a correction parameter when overdispersion occurs and also controls for serial correlation among time series counts by using a Generalised Estimating Equations (GEE) approach [62]. TRIM computes two types of indices: (i) model indices, based entirely on the model; and (ii) imputed indices based on the observations plus, for missing counts, estimated values based on the model. These indices represent change between years in terms of relative variation in total population size. Imputed indices are employed to estimate a mean annual change rate, since they show a more realistic course in time and a trend category is assigned according to its statistical significance and the estimated magnitude of change [62].
For each fox species and season, we fitted a linear trend model for the visitation index with all years as change points (time-effect model), using the complete data set. We set 2000 as the base year, so the index for 2000 took the value 1. We used a stepwise procedure to identify significant changes in the slope using Wald tests to remove nonsignificant change points from the model (default significance threshold: 0.2) according to the parsimony principle [64]. This procedure provides a high confidence on true variations in the studied populations [63]. In order to check this, we built a model where all change points were included as mandatory, and we compared it with the model built under a stepwise procedure using the Akaike Information Criterion (AIC) and considering differences between models in ΔAIC > 2 to reject the less parsimonious model [66]. The same procedure was used to evaluate the puma population trend for the 2002–2009 period, using the number of fresh signs as count data and setting 2002 as the base year. All models were run with overdispersion and serial correlation taken into account (default TRIM thresholds: >3.0 and >0.4, respectively). We estimated overall trends using the multiplicative trend slope (β), which reflects the annual average percent change [62]. When β < 1, the population decreases; when β = 1, the population remains stable; and when β > 1, the population increases. In addition, we used β to estimate the overall population size change (PSC) for the whole time series using the formula: (PSC) = (β (n − 1) − 1) × 100, where n is the total number of years for the whole monitoring period, and change is expressed as a percentage [67].
Seasonal effects on population trends of foxes were tested using the goodness of fit of models and 95% confidence intervals for the PSC. Goodness of fit of models was tested by the likelihood ratio (LR) and chi-squared tests. Nonsignificant values (i.e., p > 0.05) for one or both tests indicate that the data fitted a Poisson distribution [62].
We also evaluated how differences between sheep ranches in the cessation time of livestock activities (including predator control) affected the slopes of both fox species and puma overall trends, assuming a similar management model in both ranches. To do this, we included ranch identity as a covariate in a simple linear trend model for each species and season. The data structure prevented the fit of a time-effect model with a stepwise procedure [62]. Differences on trend slopes were evaluated by the Wald test, using a significant level value of p < 0.05 [62].
We used simple regression models to evaluate how species population trends were related to each other according to the hypothesis of size-based dominance hierarchy [37,44]. We used the trend index of the putative dominant species at year “t” as the explanatory variable for the observed trend index for the putative subordinate species at year “t + 1”, following a procedure to assess top-down regulation relationships [37,68].
Trend graphs were built using imputed trend indices, which were more realistic and more easily visualised than model indices [62]. Analyses were developed using “rtrim”, “stats” and “graphics” R packages [69].

3. Results

Stepwise procedure models (SWM) were the most parsimonious for every species (Supplementary Material Table S1). All models fitted a log-linear distribution and were not affected by overdispersion and serial correlation, except the spring model for the culpeo fox, for which both the chi-squared and the likelihood ratio tests indicated a lower model fit (Supplementary Material Table S1), which, however, should not affect the reliability of this trend [62].
According to the best model, the overall population trend for the culpeo fox was significantly decreasing in both seasons (Figure 2; Table 1), with an estimated average annual population decline of 23% in autumn and 10% in spring, which were classified by TRIM as a “steep decline” and “moderate decline”, respectively (Table 1). In terms of population size, this means a population reduction of 97% (95% CI: 92–99%) in autumn and of 75% (95% CI: 42–89%) in spring at the end of the study period. These models included five significant change points, both for autumn and spring, which were similar between seasons (Supplementary Material Table S2; Figure 2).
On the contrary, the overall population trend for grey foxes was significantly positive and similar between seasons, with an estimated average annual population increase of 7%, which was classified by TRIM as a “moderate increase” (Table 1; Figure 2). The grey fox population size doubled during the study period (a 125% increase; 95% CI: 80–180% in autumn and 44–247% in spring). These models included four significant change points, both for autumn and for spring, which were similar between seasons (Supplementary Material Table S2; Figure 2).
The puma population trend significantly increased for the studied time period (2002–2009), with an estimated annual population increase of 19%, classified by TRIM as a “steep increase” (Table 1; Figure 3). This means a tripling of its population size in the eight years of monitoring, with an estimate increase of 238% (95% CI: 43–625%). The model included three significant change points with two periods of significant increase and one with a significant decrease (Supplementary Material Table S2; Figure 3).
Observed trends were similar between abandoned ranches, following the same pattern of overall trends previously described, i.e., decreasing trends for the culpeo fox and increasing trends for the grey fox and puma. Despite this, trend slopes varied significantly between ranches for the culpeo fox in autumn (Wald test = 13.2, df = 1, p < 0.001), showing a decrease rate for CH (model slope ± SE: −0.49 ± 0.09; average annual population decline: 39%) higher than that of EC (model slope ± SE: −0.16 ± 0.08; average annual population decline: 15%; Figure 4). On the contrary, no significant differences were found during spring (Wald test = 0.059, df = 1, p > 0.1) (Supplementary Material Figure S1). Similarly, no significant differences were found for the grey fox slopes between ranches both in autumn (Wald test = 0.02, df = 1, p > 0.1; CH: 0.06 ± 0.01; EC: 0.06 ± 0.02) and spring (Wald test = 0.24, df = 1, p > 0.1; CH: 0.05 ± 0.02; EC: 0.07 ± 0.03), showing for both ranches similar average annual population increases of ~7% (Supplementary Material Figure S1). Finally, no differences were found for puma populations at the ranch level (Wald test = 9.39, df = 1, p > 0.1), showing a similar slope (CH: 0.12 ± 0.08; EC: 0.13 ± 0.08) and similar average annual population increases of ~13% (Supplementary Material Figure S1).
The grey fox showed an increasing trend with the decreasing trend of the culpeo fox population (Figure 5). Simple regression models showed significant negative relationships between the trends for the two fox species, both in autumn (R2 = 0.70, F1,10 = 23.94, p < 0.001 ) and spring (R2 = 0.79, F1,9 = 33.87, p < 0.001). A similar relationship was found between trends of the culpeo fox and puma for both seasons, with a decreasing trend of the culpeo fox mirroring the increasing trend of the puma population. These relationships were not significant either in autumn (R2 = 0.40, F1,5 = 3.37, p = 0.12) or spring (R2 = 0.36, F1,5 = 2.83, p = 0.15; Figure 5). Finally, a weak positive relationship was found between the grey fox and puma trends, both in autumn (R2 = 0.27, F1,5 = 1.90, p = 0.22) and spring (R2 = 0.40, F1,5 = 2.07, p = 0.21; Figure 5).

4. Discussion

Our results for foxes highlight the importance of designing monitoring programs on the basis of pilot studies and power analysis for the early detection of significant changes in population size [9,10]. Well-designed monitoring programs may be able to detect population change over a short time, but others may require a decade or two to show, for example, a 25% population change with 80% power at α = 0.10 [9]. Our monitoring program was designed to detect significant changes in the trends of both fox species in the first 5 years of monitoring (50% decrease and even a 30% increase in three years) with an estimated power of ~76% [30]. Thus, our monitoring program was effective because it was able to detect quick significant changes in population trends for both fox species during the first 5 years, in autumn as well as in spring. Further, the sensitivity of the monitoring program was greater than expected in detecting both significant increases and decreases [30,51].
The estimated trends were concurrent between seasons for both fox species, showing a significant decreasing trend for the culpeo fox population and a significant increasing trend for the grey fox population. The similarity between seasonal trends was maximal for the grey fox, obtaining the same mean annual change rates (7% of annual increase). Conversely, the mean annual decrease rate estimated for the culpeo fox in autumn was twice the corresponding spring value (a 23% vs 10% annual decrease, respectively). Overall, autumn trends produced more accurate estimates with narrower confidence intervals for population change, thus reducing uncertainty. Red fox (Vulpes vulpes) populations reach maximum abundance in autumn during juvenile dispersal [70,71], increasing the probability of detection and improving the accuracy of population parameter estimates during this season as compared with other seasons, including spring [72,73]. Patagonian foxes have a similar phenology and annual peaks of abundance usually occur in autumn, during juvenile dispersal [74], resulting in higher visitation rates with lower variability and more precise estimates. On the contrary, minimum population abundances are found in spring, after winter mortality, coupled with a reduction in fox activity during the breeding season [74]. Using annual troughs in fox abundance might have increased variability in visitation rate as well as uncertainty in estimates of population change. This effect was especially pronounced for the culpeo fox population, which was apparently involved in a drastic population decline. Decreasing numbers and reduced activity could explain seasonal differences in the estimated magnitude and rate of population collapse and the lack of fit in spring models [75]. Conversely, the higher abundance of the grey fox during the study period as a result of steady population increase might have diluted seasonal differences in population trend. Likewise, the sensitivity for detecting population increases is usually greater than the sensitivity for detecting declines when the sampling effort is kept constant [76].
Once fox trends were detected, the next step was inferring possible processes causing them. Population trends of culpeo and grey foxes showed a specular antagonistic pattern, where the continued increase of the grey fox population size was significantly related to the decline of the culpeo fox. Culpeo and grey foxes are opportunistic predators with similar ecological requirements and compete for the same resources [77]. With a larger body size (more accentuated in our study area; body size ratio > 3) [48], the culpeo fox could play a dominant role in competitive interactions similar to those described in other carnivore guilds [44]. Culpeo foxes may exclude grey foxes from high-quality habitats with abundant prey [78,79]. Under unfavourable conditions (e.g., reduced prey diversity and availability), competitive interactions could be more intense [50,77] and culpeo foxes could even kill grey foxes [47]. According to the competitor release effect, a reduction in the abundance of a dominant competitor allows subordinate competitors to take advantage of released resources, increasing their abundance [45,46]. An example of this ecological process has been shown with mesocarnivores in the UK, where the intense culling of the Eurasian badger Meles meles, a dominant competitor, was associated with an increase in densities of a subordinate competitor, the red fox [80]. Thus, the reduction in the population size of the dominant culpeo fox would have favoured grey fox numbers once competitive stress was relaxed.
In general, both fox species show considerable resilience to extraction from predator control, maintaining fairly stable populations over time [32,81,82]. Furthermore, when the intensity of predator persecution is reduced, culpeo foxes can respond by numerically increasing their density [81]. Hence, the deep decline observed in the culpeo fox population was an unexpected result in a protected area such as BPJNP where predator control was forbidden years before monitoring began and poaching rarely occurs. The first hypothesis that we considered to explain culpeo fox trends was the recovery of the puma population in BPJNP triggering top-down regulation processes.
On the one hand, the significant increase in puma population observed during our study suggests a population recovery similar to that described in other Patagonian regions where livestock activity and predator persecution has ceased [35]. In the absence of human persecution, large carnivores self-regulate their populations through social mechanisms and density dependence [83]. Accordingly, a significant increase in abundance would not be expected in a medium–high-density puma population, in which case most of the territories would be occupied [84,85]. Therefore, it is likely that pumas in BPJNP occurred at very low densities. Unlike foxes, predator persecution has decimated puma populations throughout Argentina, which were extirpated in many regions of Patagonia [28]. Thus, the status of the puma population at the creation of the BPJNP could have been at low density or even absent from some areas of the BPJNP, which could have favoured the release of subordinate competitors [38,42], such as the culpeo fox.
On the other hand, population trends of culpeo foxes (both autumn and spring) were negatively related to puma population trends, supporting the recovery of the regulation exerted by pumas on the culpeo fox population. Conversely, grey fox population trends were positively related to the increasing population trend of pumas, showing no regulation of pumas over grey foxes. Although these relationships were not statistically significant, the pattern observed in our study was similar to that described for the interaction between three canid species in North America [37]. This study shows how the recovery of the dominant predator, the grey wolf (Canis lupus) suppresses coyote (Canis latrans) populations, which in turn releases red foxes from regulation by coyotes. Levi and Willmers [37] point out how these interactions are determined by the trophic overlap between species, in turn related to their body size. Thus, wolves find a greater potential competitor in the larger coyote than in the smaller red fox. Similarly, pumas are more likely to focus their efforts on displacing or even killing culpeo foxes (we have direct evidences of that in the BPJNP; authors’ unpublished data) than grey foxes because they might perceive culpeo foxes as stronger competitors due to their larger overlap on food resources. Trophic niche overlap between both species can be high. In a protected area located in the high mountains of central Argentina, where native prey increased after livestock cessation, pumas and culpeo foxes showed a high diet overlap which was maintained over the 6 years of study (Pianka’s overlap index values range: 0.92–0.99) [86]. In southwestern Patagonia, pumas and culpeo foxes select the exotic European hare even in protected areas where native large-bodied guanacos occur [22,24]. In our study area, pumas mainly consume guanacos (60% of consumed biomass) and European hare (30% of consumed biomass) [49]. In contrast, European hare (36% of consumed biomass) and small rodents and armadillos (30% of consumed biomass) make up the bulk of the culpeo fox diet in BPJNP [50]. Despite the fact that the guanaco is the main prey of pumas, the European hare is an important prey for both species in BPJNP. So, it is expected that the increase in puma abundance intensified competitive interactions between both species.
Although our results would support these hypotheses, other factors not considered here could have been involved in the observed ecological processes. For example, our results partially support the hypothesis of a regulation process of the culpeo fox population by pumas. However, accepting this as the sole mechanism provoking the culpeo fox decline would imply that puma interference would be more efficient in depressing fox numbers than human persecution, which is questionable. Furthermore, marked changes in the availability of food resources advantageously exploited by each species could have benefited the grey fox or harmed the culpeo fox [48,50], explaining the observed trends in both fox species populations. On the other hand, our results in autumn support how the differences in time since the abandonment of livestock activity in ranches have significantly affected the population reduction rate of culpeo foxes. The regression was much more accentuated in the ranch where sheep were previously removed, i.e., CH, than in the ranch where about 500 sheep along with predator persecution remained in the area seven years longer, i.e., EC. Therefore, the removal of sheep from the BPJNP would be a limiting factor for the culpeo fox. Sheep have been described as an important food resource for culpeo foxes in active Patagonian ranches, both as live prey (mainly lambs and young sheep) and carrion [77]. Further, an important prey for the culpeo fox in Patagonia, the European hare [22,77], could have been harmed by the disappearance of sheep grazing that favours the proliferation of introduced lagomorphs by reducing the height of the grass and shrub cover [87,88]. However, we do not know whether management in both ranches was different, and how this could have influenced the results. Unfortunately, we do not have enough data to explain the dramatic population reduction observed for the culpeo fox in BPJNP, so we can only present our discussion as hypotheses that could have acted: puma competition, reduction in key food resources for culpeo foxes or a synergy between them. All these alternative explanations deserve further investigation to achieve a better understanding of the competitive interactions between these species.
Our long-term study highlights the importance of wildlife monitoring programs not only as a tool for an effective management and conservation of biodiversity [8], but also for proposing hypotheses about ecological processes and interactions [12,13,15]. The implementation of monitoring programs for carnivores in newly protected areas created on human-altered landscapes can be a key tool to assess the relevance of predators to restore degraded ecosystems and confer resilience against globally threatening processes [40]. Our long-term monitoring may provide insight into the poorly understood competitive interactions among the three most common carnivore species of Patagonian ecosystems, and further specific studies are required to assess the hypotheses proposed here.

Supplementary Materials

The following materials are available online at https://www.mdpi.com/1424-2818/12/9/319/s1; Figure S1: Non-significant ranch effect (EC: El Cuadro; CH: Cerro Horqueta) on population trends of culpeo fox, grey fox and puma for each season, Table S1: Trend models fitted for each species and season, and ordered according to their AIC values, Table S2: Model slope and standard error (SE) for change points selected in the stepwise modelling process for the three species and seasons.

Author Contributions

Conceptualization, A.T., A.R. and F.D.-R.; methodology, A.T., A.R. and F.D.-R.; formal analysis, F.D.-R. and D.P.; Investigation (field work), A.T., D.P., S.Z. and J.I.Z.-M.; writing—original draft preparation, F.D.-R. and A.T.; writing—review and editing, F.D.-R. and A.R.; project administration and funding acquisition, A.T. and A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FONCYT through the projects PICT (grant PICT-2012-0262), the BBVA Foundation through a grant under the Conservation Biology Program (to A. Rodríguez), Universidad Nacional de la Patagonia Austral through project B-013, CONICET (PEI-6065), the Rufford Foundation (grant RSG-30.06.09) and the Spanish Ministry of Science and Innovation (grant CGL2011-27469), and co-funded by the European Regional Development Fund (FEDER).

Acknowledgments

We thank the support received from the staff of Bosques Petrificados de Jaramillo National Park. We also thank the collaboration in the fieldwork from R. Martínez Peck, M. Santillán, C. Zoratti, G. Soria, F. Escobar, G. Aguilera, E. Aguilar, P. Collavino, E. Daher, M. Yaya, D. Brescia and M. Bronfman. D. Procopio was supported by a CONICET predoctoral fellowship. Francisco Díaz-Ruiz enjoyed a postdoctoral research contract, “Juan de la Cierva” (ref: FJCI-2015-24949), from the Spanish Ministry of Economy, Industry and Competitiveness, and his research stay in UNPA (Argentina) was financed by a research grant from the program “Becas Iberoamérica-Santander Universidades” (2015 Call) and a postdoctoral contract financed by the European Social Fund (ESF) and the Junta de Comunidades de Castilla-La Mancha (Operational Programme FSE 2007/2013). David McNitt (Virginia Polytechnic Institute and State University, USA) carefully revised the English text.

Conflicts of Interest

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

References

  1. Caro, T.; Darwin, J.; Forrester, T.; Ledoux-Bloom, C.; Wells, C. Conservation in the Anthropocene. Conserv. Biol. 2012, 26, 185–188. [Google Scholar] [CrossRef] [PubMed]
  2. Newbold, T.; Hudson, L.N.; Phillips, H.R.P.; Hill, S.L.L.; Contu, S.; Lysenko, I.; Blandon, A.; Butchart, S.H.M.; Booth, H.L.; Day, J.; et al. A global model of the response of tropical and sub-tropical forest biodiversity to anthropogenic pressures. Proc. R. Soc. B Biol. Sci. 2014, 281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Venter, O.; Sanderson, E.W.; Magrach, A.; Allan, J.R.; Beher, J.; Jones, K.R.; Possingham, H.P.; Laurance, W.F.; Wood, P.; Fekete, B.M.; et al. Sixteen years of change in the global terrestrial human footprint and implications for biodiversity conservation. Nat. Commun. 2016, 7, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Balmford, B.; Green, R.E.; Onial, M.; Phalan, B.; Balmford, A. How imperfect can land sparing be before land sharing is more favourable for wild species? J. Appl. Ecol. 2019, 56, 73–84. [Google Scholar] [CrossRef] [Green Version]
  5. Le Saout, S.; Hoffmann, M.; Shi, Y.; Hughes, A.; Bernard, C.; Brooks, T.M.; Bertzky, B.; Butchart, S.H.M.; Stuart, S.N.; Badman, T.; et al. Protected areas and effective biodiversity conservation. Science 2014, 342, 803–805. [Google Scholar] [CrossRef]
  6. Naughton-Treves, L.; Holland, M.B.; Brandon, K. The role of protected areas in conserving biodiversity and sustaining local livelihoods. Annu. Rev. Environ. Resour. 2005, 30, 219–252. [Google Scholar] [CrossRef] [Green Version]
  7. UNEP-WCMC; IUCN Protected Planet. The World Database on Protected Areas (WDPA), April 2020 Version; UNEP-WCMC: Cambridge, UK, 2020. [Google Scholar]
  8. Kiffner, C.; Binzen, G.; Cunningham, L.; Jones, M.; Spruiell, F.; Kioko, J. Wildlife population trends as indicators of protected area effectiveness in northern Tanzania. Ecol. Indic. 2020, 110. [Google Scholar] [CrossRef]
  9. Gitzen, R.A.; Millspaugh, J.J.; Cooper, A.B.; Licht, D.S. Design and Analysis of Long-term Ecological Monitoring Studies; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  10. Lindenmayer, D.B.; Piggott, M.P.; Wintle, B.A. Counting the books while the library burns: Why conservation monitoring programs need a plan for action. Front. Ecol. Environ. 2013, 11, 549–555. [Google Scholar] [CrossRef] [Green Version]
  11. Westgate, M.J.; Likens, G.E.; Lindenmayer, D.B. Adaptive management of biological systems: A review. Biol. Conserv. 2013, 158, 128–139. [Google Scholar] [CrossRef]
  12. Geupel, G.R.; Humple, D.; Roberts, L.J. Monitoring decisions: Not as simple as they seem? Trends Ecol. Evol. 2011, 26, 107, author reply 108-9. [Google Scholar] [CrossRef]
  13. Gitzen, R.A.; Millspaugh, J.J. Ecological monitoring: The heart of the matter. In Design and Analysis of Long-Term Ecological Monitoring Studies; Gitzen, R.A., Millspaugh, J.J., Cooper, A., Licht, D.S., Eds.; Cambridge University Press: Cambridge, UK, 2012; pp. 3–22. [Google Scholar]
  14. Pringle, R.M. Upgrading protected areas to conserve wild biodiversity. Nature 2017, 546, 91–99. [Google Scholar] [CrossRef] [PubMed]
  15. Cross, S.L.; Bateman, P.W.; Cross, A.T. Restoration goals: Why are fauna still overlooked in the process of recovering functioning ecosystems and what can be done about it? Ecol. Manag. Restor. 2020, 21, 4–8. [Google Scholar] [CrossRef]
  16. National Parks Agency, Argentine Government. Parque Nacional de Monte León. Available online: https://www.argentina.gob.ar/parquesnacionales/monteleon (accessed on 20 May 2020).
  17. National Parks Agency, Chile Government. Parque Nacional Patagonia. Available online: https://www.conaf.cl/parques/parque-nacional-patagonia/ (accessed on 20 May 2020).
  18. Cibils, A.; Borrelli, P.R. Grasslands of Patagonia. In Grasslands of the World; Suttie, J.M., Reynolds, S.G., Batello, C., Eds.; Food and Agriculture Organization of the United Nations: Rome, Italy, 2005; pp. 121–170. [Google Scholar]
  19. Grigera, D.E.; Rapoport, E.H. Status and distribution of the European hare in South America. J. Mammal. 1983, 64, 163–166. [Google Scholar] [CrossRef]
  20. Bonino, N.; Coss, D.; Menegheti, J. Dispersal of the European hare, in South America. Folia Zool. 2010, 59, 9–15. [Google Scholar] [CrossRef]
  21. Barri, F.R.; Martella, M.B.; Navarro, J.L. Nest-site habitat selection by Lesser Rheas (Rhea pennata pennata) in northwestern Patagonia, Argentina. J. Ornithol. 2009, 150, 511–514. [Google Scholar] [CrossRef]
  22. Novaro, A.J.; Funes, M.C.; Walker, R.S. Ecological extinction of native prey of a carnivore assemblage in Argentine Patagonia. Biol. Conserv. 2000, 92, 25–33. [Google Scholar] [CrossRef]
  23. Travaini, A.; Zapata, S.C.; Bustamante, J.; Pedrana, J.; Zanón, J.I.; Rodríguez, A. Guanaco abundance and monitoring in southern patagonia: Distance sampling reveals substantially greater numbers than previously reported. Zool. Stud. 2015, 54. [Google Scholar] [CrossRef] [Green Version]
  24. Buenavista, S.; Palomares, F. The role of exotic mammals in the diet of native carnivores from South America. Mamm. Rev. 2018, 48, 37–47. [Google Scholar] [CrossRef]
  25. Travaini, A.; Zapata, S.; Martinez Peck, R.; Delibes, M. Percepción y actitud humanas hacia la prelación de ganado ovino por el zorro colorado (Pseudalopex culpaeus) en Santa Cruz, Patagonia Argentina. Mastozoología Neotrop. 2000, 7, 117–129. [Google Scholar]
  26. Funes, M.; Novaro, A. Rol de la fauna silvestre en la economia del poblador rural, provincia del Neuquen, Argentina. Rev. Argent. Prod. Anim. 1999, 19, 265–271. [Google Scholar]
  27. Elbroch, L.M.; Kusler, A. Are pumas subordinate carnivores, and does it matter? PeerJ 2018, 6, e4293. [Google Scholar] [CrossRef] [PubMed]
  28. Guerisoli, M.; de las, M.; Luengos Vidal, E.; Franchini, M.; Caruso, N.; Casanave, E.B.; Lucherini, M. Characterization of puma–livestock conflicts in rangelands of central Argentina. R. Soc. Open Sci. 2017, 4, 170852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Llanos, R.; Travaini, A.; Montanelli, S.; Crespo, E. Estructura de edades de pumas (Puma concolor) cazados bajo el sistema de remoción por recompensas en Patagonia. ¿Selectividad u oportunismo en la captura? Ecol. Austral. 2014, 24, 311–319. [Google Scholar]
  30. Travaini, A.; Rodríguez, A.; Procopio, D.; Zapata, S.C.; Zanón, J.I.; Martínez-Peck, R. A monitoring program for Patagonian foxes based on power analysis. Eur. J. Wildl. Res. 2010, 56, 421–433. [Google Scholar] [CrossRef]
  31. Bellati, J.; von Thungen, J. Lamb predation in Patagonian ranches. In Proceedings of the 14th Vertebrate Pest Conference, Sacramento, CA, USA, 6–8 March 1990; Volume 6, pp. 263–268. [Google Scholar]
  32. Novaro, A.J. Sustainability of harvest of culpeo foxes in Patagonia. Oryx 1995, 29, 18–22. [Google Scholar] [CrossRef]
  33. Mueller, J.P. Introducción a la producción ovina Argentina. Available online: http://www.produccion-animal.com.ar/ (accessed on 3 May 2020).
  34. Andrade, L. Sociología de la Desertificación. Los Productores Ovinos Extensivos de la Patagonia Austral; Minio & Dávila Editores: Buenos Aires, Argentina, 2005. [Google Scholar]
  35. Walker, S.; Novaro, A. The world’s southernmost pumas in Patagonia and the southern Andes. In Cougar: Ecology and Conservation; Hornocker, M., Negri, S., Eds.; University of Chicago Press: Chicago, IL, USA, 2010; pp. 91–99. [Google Scholar]
  36. Caro, J.; Zapata, S.C.; Zanón, J.I.; Rodríguez, A.; Travaini, A. Ganadería ovina y usos alternativos del suelo en la Patagonia austral, Argentina. Multequina 2017, 26, 33–50. [Google Scholar]
  37. Levi, T.; Wilmers, C.C. Wolves-coyotes-foxes: A cascade among carnivores. Ecology 2012, 93, 921–929. [Google Scholar] [CrossRef] [Green Version]
  38. Prugh, L.R.; Stoner, C.J.; Epps, C.W.; Bean, W.T.; Ripple, W.J.; Laliberte, A.S.; Brashares, J.S. The rise of the mesopredator. Bioscience 2009, 59, 779–791. [Google Scholar] [CrossRef]
  39. Estes, J.A.; Terborgh, J.; Brashares, J.S.; Power, M.E.; Berger, J.; Bond, W.J.; Carpenter, S.R.; Essington, T.E.; Holt, R.D.; Jackson, J.B.C.; et al. Trophic downgrading of planet earth. Science 2011, 333, 301–306. [Google Scholar] [CrossRef] [Green Version]
  40. Ritchie, E.G.; Elmhagen, B.; Glen, A.S.; Letnic, M.; Ludwig, G.; McDonald, R.A. Ecosystem restoration with teeth: What role for predators? Trends Ecol. Evol. 2012, 27, 265–271. [Google Scholar] [CrossRef]
  41. Polis, G.A.; Holt, R.D. Intraguild predation: The dynamics of complex trophic interactions. Trends Ecol. Evol. 1992, 7, 151–154. [Google Scholar] [CrossRef]
  42. Ritchie, E.G.; Johnson, C.N. Predator interactions, mesopredator release and biodiversity conservation. Ecol. Lett. 2009, 12, 982–998. [Google Scholar] [CrossRef]
  43. Ripple, W.J.; Wirsing, A.J.; Wilmers, C.C.; Letnic, M. Widespread mesopredator effects after wolf extirpation. Biol. Conserv. 2013, 160, 70–79. [Google Scholar] [CrossRef] [Green Version]
  44. Monterroso, P.; Díaz-Ruiz, F.; Lukacs, P.M.; Alves, P.C.; Ferreras, P. Ecological traits and the spatial structure of competitive coexistence among carnivores. Ecology 2020, 101, e03059. [Google Scholar] [CrossRef] [PubMed]
  45. Caut, S.; Casanovas, J.G.; Virgos, E.; Lozano, J.; Witmer, G.W.; Courchamp, F. Rats dying for mice: Modelling the competitor release effect. Austral. Ecol. 2007, 32, 858–868. [Google Scholar] [CrossRef] [Green Version]
  46. Larsen, E. Competitive release in microhabitat use among coexisting desert rodents: A natural experiment. Oecologia 1986, 69, 231–237. [Google Scholar] [CrossRef]
  47. De Oliveira, T.G.; Pereira, J.A. Intraguild predation and interspecific killing as structuring forces of carnivoran communities in South America. J. Mamm. Evol. 2014, 21, 427–436. [Google Scholar] [CrossRef]
  48. Zapata, S.C.; Procopio, D.; Martínez-Peck, R.; Zanón, J.I.; Travaini, A. Morfometría externa y reparto de recursos en zorros simpátricos (Pseudalopex culpaeus y P. griseus) en el sureste de la Patagonia Argentina. Mastozoología Neotrop. 2008, 15, 103–111. [Google Scholar]
  49. Zanón Martínez, J.I.; Travaini, A.; Zapata, S.; Procopio, D.; Santillán, M.Á. The ecological role of native and introduced species in the diet of the puma Puma concolor in southern Patagonia. Oryx 2012, 46, 106–111. [Google Scholar] [CrossRef] [Green Version]
  50. Zapata, S.C.; Travaini, A.; Delibes, M.; Martínez-Peck, R. Food habits and resource partitioning between grey and culpeo foxes in southeastern Argentine Patagonia. Stud. Neotrop. Fauna Environ. 2005, 40, 97–103. [Google Scholar] [CrossRef]
  51. Travaini, A.; Zapata, S.C.; Zoratti, C.; Soria, G.; Escobar, F.; Aguilera, G.; Collavino, P. Diseño de un programa de seguimiento de poblaciones de cánidos silvestres en ambientes esteparios de la Patagonia, Argentina. Acta Zool. Mex. 2003, 90, 1–14. [Google Scholar]
  52. Olson, D.M.; Dinerstein, E. The Global 200: Priority ecoregions for global conservation. Ann. Mo. Bot. Gard. 2002, 89, 199–224. [Google Scholar] [CrossRef]
  53. Soriano, A. Deserts and semi-deserts of Patagónia. In Temperate Deserts and Semideserts; West, N.E., Ed.; Elsevier: Amsterdam, The Netherlands, 1983; pp. 423–459. [Google Scholar]
  54. Ares, J. Structural and dynamic characteristics of overgrazed lands of northern Patagonia, Argentina. In Managed Grasslands; Bremeyer, A., Ed.; Elsevier: Amsterdam, The Netherlands, 1990; pp. 149–175. [Google Scholar]
  55. Bertiller, M.B.; Bisigato, A. Vegetation dynamics under grazing disturbance. The state-and-transition model for the Patagonian steppes. Ecol. Austral. 1998, 8, 191–199. [Google Scholar]
  56. Ray, J.C.; Zielinski, W.J. Track stations. In Noninvasive Survey Methods for Carnivores; Long, R.A., MacKay, P., Zielinski, W.J., Ray, J.C., Eds.; Island Press: Washington, DC, USA, 2008; pp. 75–190. [Google Scholar]
  57. Linhart, S.B.; Knowlton, F.F. Determining the Relative Abundance of Coyotes by Scent Station Lines. Wildl. Soc. Bull. 1975, 3, 119–124. [Google Scholar]
  58. Travaini, A.; Peck, R.M.; Zapata, S.C.; Delibes, M. Selection of odor attractants and meat delivery methods to control culpeo foxes (Pseudalopex culpaeus) in Patagonia. Wildl. Soc. Bull. 2001, 29, 1089–1096. [Google Scholar] [CrossRef]
  59. Shaw, H.G.; Beier, P.; Culver, M.; Grigione, M. Puma Field Guide. The Cougar Network; Cougarnet.org: Boston, MA, USA, 2007. [Google Scholar]
  60. Harmsen, B.J.; Foster, R.J.; Gutierrez, S.M.; Marin, S.Y.; Doncaster, C.P. Scrape-marking behavior of jaguars (Panthera onca) and pumas (Puma concolor). J. Mammal. 2010, 91, 1225–1234. [Google Scholar] [CrossRef] [Green Version]
  61. Bogaart, P.; Van Der Loo, M.; Pannekoek, J. Rtrim: Trends and Indices for Monitoring Data. R Package Version 2.0.6. 2018. Available online: https://CRAN.R-project.org/package=rtrim (accessed on 1 August 2020).
  62. Pannekoek, J.; Van Strien, A. TRIM 3.53 Manual. Trends and Indices for Monitoring Data; Voorburg: Statistics Netherlands: Den Haag, The Netherlands, 2006. [Google Scholar]
  63. Loman, J.; Andersson, G. Monitoring brown frogs Rana arvalis and Rana temporaria in 120 south Swedish ponds 1989-2005. Mixed trends in different habitats. Biol. Conserv. 2007, 135, 46–56. [Google Scholar] [CrossRef]
  64. Underwood, J.G.; Silbernagle, M.; Nishimoto, M.; Uyehara, K. Managing conservation reliant species: Hawai’i’s endangered endemic waterbirds. PLoS ONE 2013, 8, e67872. [Google Scholar] [CrossRef] [Green Version]
  65. Kasahara, S.; Koyama, K. Population trends of common wintering waterfowl in Japan: Participatory monitoring data from 1996 to 2009. Ornithol. Sci. 2010, 36, 23–36. [Google Scholar] [CrossRef]
  66. Burnham, K.P.; Anderson, D.R. Model. Selection and Multimodel Inference; Burnham, K.P., Anderson, D.R., Eds.; Springer: New York City, NY, USA, 2002. [Google Scholar]
  67. Gómez-Catasús, J.; Pérez-Granados, C.; Barrero, A.; Bota, G.; Giralt, D.; López-Iborra, G.M.; Serrano, D.; Traba, J. European population trends and current conservation status of an endangered steppe-bird species: The Dupont’s lark Chersophilus duponti. PeerJ 2018, 2018, 1–19. [Google Scholar] [CrossRef] [Green Version]
  68. Fernandez-de-Simon, J.; Díaz-Ruiz, F.; Rodríguez-de la Cruz, M.; Delibes-Mateos, M.; Villafuerte, R.; Ferreras, P. Can widespread generalist predators affect keystone prey? A case study with red foxes and European rabbits in their native range. Popul. Ecol. 2015, 57, 591–599. [Google Scholar] [CrossRef] [Green Version]
  69. R CORE TEAM. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  70. Saunders, G.; Coman, B.J.; Kinnear, J.L.; Braysher, M.L. Managing Vertebrate Pests: Foxes; Saunders, G., Coman, B.J., Kinnear, J.L., Braysher, M.L., Eds.; Australian Government Publishing Service: Canberra, Australia, 1995. [Google Scholar]
  71. Macdonald, D.W.; Reynolds, J.C. Red fox Vulpes vulpes Linnaeus, 1758. In Canids: Foxes, Wolves, Jackals and Dogs. Status Survey and Conservation Action Plan; Sillero-Zubiri, C., Hoffmann, M., Macdonald, D.W., Eds.; IUCN/SSC Canid Specialist Group: Cambridge, UK, 2004; pp. 129–136. [Google Scholar]
  72. Field, S.A.; Tyre, A.J.; Thorn, K.H.; O’Connor, P.J.; Possingham, H.P. Improving the efficiency of wildlife monitoring by estimating detectability: A case study of foxes (Vulpes vulpes) on the Eyre Peninsula, South Australia. Wildl. Res. 2005, 32, 253–258. [Google Scholar] [CrossRef]
  73. Vine, S.J.; Crowther, M.S.; Lapidge, S.J.; Dickman, C.R.; Mooney, N.; Piggott, M.P.; English, A.W. Comparison of methods to detect rare and cryptic species: A case study using the red fox (Vulpes vulpes). Wildl. Res. 2009, 36, 436–446. [Google Scholar] [CrossRef]
  74. Sillero-Zubiri, C.; Hoffmann, M.; Macdonald, D.W. Canids: Foxes, Wolves, Jackals and Dogs. Status Survey and Conservation Action Plan; IUCN/SSC Canid Specialist Group: Cambridge, UK, 2004. [Google Scholar]
  75. Bart, J.; Beyer, H.L. Analysis options for estimating status and trends in long-term monitoring. In Design and Analysis of Long-term Ecological Monitoring Studies; Gitzen, R.A., Millspaugh, J.J., Cooper, A.B., Licht, D.S., Eds.; Cambridge University Press: Cambridge, UK, 2012; pp. 253–278. [Google Scholar]
  76. Gibbs, J.P. Monitor: User’s Manual. Available online: https://www.esf.edu/efb/gibbs/monitor/ (accessed on 5 January 2020).
  77. Novaro, A.J.; Funes, M.C.J.; Jiménez, J. Selection for introduced prey and conservation of culpeo and chilla zorros in Patagonia. In The Biology and Conservation of Wild Canids; Macdonald, D., Sillero-Zubiri, C., Eds.; Oxford University Press: Oxford, UK, 2004; pp. 243–254. [Google Scholar]
  78. Johnson, W.E.; Franklin, W.L. Spatial resource partitioning by sympatric grey fox (Dusicyon griseus) and culpeo fox (Dusicyon culpaeus) in southern Chile. Can. J. Zool. 1994, 72, 1788–1793. [Google Scholar] [CrossRef]
  79. Jiménez, J.E.; Yáñez, J.L.; Tabilo, E.L.; Jaksic, F.M. Niche-complementarity of South American foxes: Reanalysis and test of a hypothesis. Rev. Chil. Hist. Nat. 1996, 69, 113–123. [Google Scholar]
  80. Trewby, I.D.; Wilson, G.J.; Delahay, R.J.; Walker, N.; Young, R.; Davison, J.; Cheeseman, C.; Robertson, P.A.; Gorman, M.L.; Mcdonald, R.A. Experimental evidence of competitive release in sympatric carnivores. Biol. Lett. 2008, 4, 170–172. [Google Scholar] [CrossRef] [Green Version]
  81. Jiménez, J.E.; Novaro, A.J. Culpeo Pseudalopex culpaeus (Molina, 1782). In Canids: Foxes, Wolves, Jackals and Dogs. Status Survey and Conservation Action Plan; Sillero-Zubiri, C., Hoffmann, M., Macdonald, D.W., Eds.; IUCN/SSC Canid Specialist Group: Gland, Switzerland; Cambridge, UK, 2004; pp. 44–49. [Google Scholar]
  82. González del Solar, R.J.R. Chilla Pseudalopex griseus (Gray, 1837). In Canids: Foxes, Wolves, Jackals and Dogs. Status Survey and Conservation Action Plan; Sillero-Zubiri, C., Hoffmann, M., Macdonald, D.W., Eds.; IUCN/SSC Canid Specialist Group: Gland, Switzerland; Cambridge, UK, 2004; pp. 53–63. [Google Scholar]
  83. Wallach, A.D.; Izhaki, I.; Toms, J.D.; Ripple, W.J.; Shanas, U. What is an apex predator? Oikos 2015, 124, 1453–1461. [Google Scholar] [CrossRef]
  84. Stoner, D.C.; Wolfe, M.L.; Choate, D.M. Cougar exploitation levels in Utah: Implications for demographic structure, population recovery, and metapopulation dynamics. J. Wildl. Manag. 2006, 70, 1588–1600. [Google Scholar] [CrossRef]
  85. Cooley, H.S.; Wielgus, R.B.; Koehler, G.; Maletzke, B. Source populations in carnivore management: Cougar demography and emigration in a lightly hunted population. Anim. Conserv. 2009, 12, 321–328. [Google Scholar] [CrossRef]
  86. Pia, M.V. Trophic interactions between puma and endemic culpeo fox after livestock removal in the high mountains of central Argentina. Mammalia 2013, 77, 273–283. [Google Scholar] [CrossRef]
  87. Simonetti, J.A. Effect of goats upon native rodents and european rabbits in chilean matorral. Rev. Chil. Hist. Nat. 1983, 56, 27–30. [Google Scholar]
  88. Meyers, K.; Parer, I.; Wood, D.; Cooke, B.D. The rabbit in Australia. In The European Rabbit: The History and Biology of a Successful Colonizer; Thompson, H.V., King, C.M., Eds.; Oxford University Press: Oxford, UK, 1994; pp. 108–157. [Google Scholar]
Figure 1. Left: Location of the Bosques Petrificados de Jaramillo National Park (BPJNP) in the Santa Cruz province (grey region with oblique bars), southern Argentine Patagonia (grey shaded regions). Right: spatial distribution of bait station lines designed for fox population monitoring (black dots), and sampled grids for monitoring the puma population in each abandoned sheep ranch.
Figure 1. Left: Location of the Bosques Petrificados de Jaramillo National Park (BPJNP) in the Santa Cruz province (grey region with oblique bars), southern Argentine Patagonia (grey shaded regions). Right: spatial distribution of bait station lines designed for fox population monitoring (black dots), and sampled grids for monitoring the puma population in each abandoned sheep ranch.
Diversity 12 00319 g001
Figure 2. Culpeo (black dots) and grey fox (grey triangles) population trends based on the imputed trend indices estimated by the most parsimonious time-effect model (stepwise time effect model; Supplementary Material Table S1) for each season: autumn (2000–2013) and spring (2000–2012). Time-points identified by the model as significant change points on the population trends are marked with asterisks (*) for the culpeo fox and hashes (#) for the grey fox. Vertical bars represent the SE of trend indices. Horizontal dashed lines represent the index value for the base year (2000).
Figure 2. Culpeo (black dots) and grey fox (grey triangles) population trends based on the imputed trend indices estimated by the most parsimonious time-effect model (stepwise time effect model; Supplementary Material Table S1) for each season: autumn (2000–2013) and spring (2000–2012). Time-points identified by the model as significant change points on the population trends are marked with asterisks (*) for the culpeo fox and hashes (#) for the grey fox. Vertical bars represent the SE of trend indices. Horizontal dashed lines represent the index value for the base year (2000).
Diversity 12 00319 g002
Figure 3. Puma population trends based on the imputed trend indices estimated by the most parsimonious time-effect model (stepwise time effect model; Supplementary Material Table S1) for autumn (2002–2009). Time points incorporated in the model as significant change points on population trend are marked with asterisks (*). Vertical bars represent the SE of trend indices. The horizontal dashed line represents the index value for the base year (2002).
Figure 3. Puma population trends based on the imputed trend indices estimated by the most parsimonious time-effect model (stepwise time effect model; Supplementary Material Table S1) for autumn (2002–2009). Time points incorporated in the model as significant change points on population trend are marked with asterisks (*). Vertical bars represent the SE of trend indices. The horizontal dashed line represents the index value for the base year (2002).
Diversity 12 00319 g003
Figure 4. A significant ranch effect (EC: El Cuadro ranch; CH: Cerro Horqueta ranch) on the population trends of the culpeo fox in autumn. Population trends are based on the imputed trend indices estimated by a simple linear trend model (without change points; Supplementary Material Table S1). Vertical bars represent the SE of trend indices. The horizontal dashed line represents the index value for the base year (2000).
Figure 4. A significant ranch effect (EC: El Cuadro ranch; CH: Cerro Horqueta ranch) on the population trends of the culpeo fox in autumn. Population trends are based on the imputed trend indices estimated by a simple linear trend model (without change points; Supplementary Material Table S1). Vertical bars represent the SE of trend indices. The horizontal dashed line represents the index value for the base year (2000).
Diversity 12 00319 g004
Figure 5. Relationships between culpeo fox and grey fox population trends (a) and puma and both fox species’ population trends (b,c) for autumn (black dots) and spring (grey diamonds). Goodness of fit for each season model is represented by the coefficient of determination (R2).
Figure 5. Relationships between culpeo fox and grey fox population trends (a) and puma and both fox species’ population trends (b,c) for autumn (black dots) and spring (grey diamonds). Goodness of fit for each season model is represented by the coefficient of determination (R2).
Diversity 12 00319 g005
Table 1. Population trends of the culpeo fox, grey fox and puma in the BPJNP from count data spanning 2000–2013. The model slope and the mean annual percent change (+/−) are shown.
Table 1. Population trends of the culpeo fox, grey fox and puma in the BPJNP from count data spanning 2000–2013. The model slope and the mean annual percent change (+/−) are shown.
SpeciesSeasonModel Slope (SE)Annual Change (%)Long-Term Trends
Culpeo foxAutumn−0.26 (0.04)−23Steep decline ***
Spring−0.1 (0.03)−10Moderate decline *
Grey foxAutumn0.06 (0.01)7Moderate increase ***
Spring0.06 (0.01)7Moderate increase *
PumaSpring/Autumn0.17 (0.05)19Steep increase *
* (p < 0.05); *** (p < 0.001).

Share and Cite

MDPI and ACS Style

Díaz-Ruiz, F.; Rodríguez, A.; Procopio, D.; Zapata, S.; Zanón-Martínez, J.I.; Travaini, A. Inferring Species Interactions from Long-Term Monitoring Programs: Carnivores in a Protected Area from Southern Patagonia. Diversity 2020, 12, 319. https://doi.org/10.3390/d12090319

AMA Style

Díaz-Ruiz F, Rodríguez A, Procopio D, Zapata S, Zanón-Martínez JI, Travaini A. Inferring Species Interactions from Long-Term Monitoring Programs: Carnivores in a Protected Area from Southern Patagonia. Diversity. 2020; 12(9):319. https://doi.org/10.3390/d12090319

Chicago/Turabian Style

Díaz-Ruiz, Francisco, Alejandro Rodríguez, Diego Procopio, Sonia Zapata, Juan Ignacio Zanón-Martínez, and Alejandro Travaini. 2020. "Inferring Species Interactions from Long-Term Monitoring Programs: Carnivores in a Protected Area from Southern Patagonia" Diversity 12, no. 9: 319. https://doi.org/10.3390/d12090319

APA Style

Díaz-Ruiz, F., Rodríguez, A., Procopio, D., Zapata, S., Zanón-Martínez, J. I., & Travaini, A. (2020). Inferring Species Interactions from Long-Term Monitoring Programs: Carnivores in a Protected Area from Southern Patagonia. Diversity, 12(9), 319. https://doi.org/10.3390/d12090319

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

Article Metrics

Back to TopTop