Exotic Prey Facilitate Coexistence between Pumas and Culpeo Foxes in the Andes of Central Chile

: Coexistence between species with similar ecological niches implies species must segregate along one or more niche axes to survive. Space, time, and trophic resources are regarded as the principal axes upon which species segregate. We examined segregation along these niche axes to determine mechanisms underlying coexistence between the two main predators, puma ( Puma concolor ) and culpeo foxes ( Lycalopex culpaeus ) in the Andes of Central Chile. We used occupancy modeling to examine space use and overlap, Kernel Density Estimation to determine temporal activity patterns and overlap, and analysis of prey remains in feces to assess diet breadth and similarity. We found high spatial overlap and positive associations between detection of the carnivores lending little support for spatial segregation. Similarly, we found high nocturnal, temporal overlap between pumas and foxes that matched peaks in activity of prey. In contrast, we found relatively low dietary overlap indicating niche segregation likely occurs along the dietary axis. The Puma diet was dominated by introduced, exotic hares and foxes appeared to shift away from hares to rabbits, small mammals, and seeds. Given that lagomorphs are the main dietary resource for pumas in particular, management decisions regarding the control or eradication of such exotic species could negatively a ﬀ ected puma survival. Author Contributions: Conceptualization, C.O., M.K., and C.B.; methodology, C.O. and M.K.; formal analysis, C.O. and M.K.; investigation, C.O., A.M. and N.G.; resources, C.O., C.B.; data curation, C.O., N.G.; writing—original draft preparation, C.O.; writing—review and editing, M.K.; visualization, C.O.; supervision, M.K. and C.B.; project administration, A.M. and C.O.; funding acquisition, A.M., N.G., C.B.

Previous studies on spatial associations between carnivores have reported spatial segregation. For instance, red foxes (Vulpes vulpes) and Iberian lynxes (Lynx pardinus) show habitat segregation in Europe and are sensitive to prey abundance, suggesting that, when active, foxes avoided Iberian lynxes using habitats not used by them [12]. Others found that carnivores co-occur to a high degree, The study area is located in the mountain range of the Andes, in the Mediterranean-steppe transition that belongs to the Chilean Winter Rainfall Mediterranean region-Valdivian Forests Biodiversity Hotspot [34]. Weather corresponds to a Mediterranean climate, with cold rainy winters (June-September) and a summer dry season (December-March) [35]. DiCastri and Hajek [36] reported that the annual averages for maximum, mean, and minimum temperatures recorded in a climate station near the study area (i.e., the town of Sewell) are 14.5 °C, 9.5 °C, and 5.0 °C, respectively. The relative humidity (annual mean) is 50% with a precipitation annual mean of 1052 mm.
Topography is dominated by glacier basins embedded in steep mountain chains [32]. Vegetation comprises three broad types forming an altitudinal gradient between a dense sclerophyllous forest dominated by litre (  The study area is located in the mountain range of the Andes, in the Mediterranean-steppe transition that belongs to the Chilean Winter Rainfall Mediterranean region-Valdivian Forests Biodiversity Hotspot [34]. Weather corresponds to a Mediterranean climate, with cold rainy winters (June-September) and a summer dry season (December-March) [35]. DiCastri and Hajek [36] reported that the annual averages for maximum, mean, and minimum temperatures recorded in a climate station near the study area (i.e., the town of Sewell) are 14.5 • C, 9.5 • C, and 5.0 • C, respectively. The relative humidity (annual mean) is 50% with a precipitation annual mean of 1052 mm.
RCNR belongs to the national protected area system of Chile, was created in 1985, and protects recovering Mediterranean vegetation and diverse wildlife including the puma, the culpeo fox, and a few remaining guanacos (about 100 individuals, and decreasing, N. Guarda, comm. pers. 2020), which are the largest native herbivore of Chile ( Figure A2).
Native prey at our study site include guanacos (Lama guanicoe) and plain viscachas (Lagidium viscacia) [38][39][40], (which occur at very low abundances at the study area [32]) and small mammals (mostly rodents). European hares (Lepus europaeus) and rabbits (Oryctolagus cuniculus) are both exotic lagomorphs well adapted to the Mediterranean ecosystems of central Chile, introduced at the end of the 19th century [41,42]. European hares reach 3-5 kg while rabbits are smaller (1.5-2.5 kg) [43][44][45]. Hares are an important prey item for both pumas and culpeo foxes at other sites in Chile and Argentina [38][39][40]46,47], whereas rabbits are an important prey item for culpeo foxes [48][49][50]. Previous studies have assessed the contribution of exotic species as prey for pumas in other biomes of Chile [38,39,51,52]. However, the ecological and trophic relationships between pumas and culpeo foxes, and the contribution of exotic lagomorphs to their diets, have not been assessed in the Chilean Andean ecoregion.
The viscacha is a medium size rodent that was previously extant but now is extremely rare in the study area. Introduced rabbits and hares are abundant at the study site ( Figure A2). A rich biodiversity of small mammals from central Chile is well represented. A minimum puma density between 0.3 and 0.75 adult individuals/100 km 2 has been estimated based on individual recognition from camera trap records in RCNR and VV&LC [29] and culpeo fox density remains unknown. VV&LC are used as summer ranges for livestock with the animals arriving in spring and leaving in the early autumn. Although cattle and horses are grazed in the area with little supervision, sheep and goats are managed more intensively (i.e., goats are milked daily and spend the night in enclosures near houses) [33].

Camera-Trap Surveys and Trap Rates
We deployed a total of 50 single-camera stations in the study area: 36 in RCNR in 2015, and 14 in VV&LC in 2012-2013 ( Figure 1). We set cameras along valley bottoms and occasionally on ridges, using a line-sampling scheme [29] that has been reported as a successful approach for snow leopards (Panthera uncia) in mountainous environments [53]. Cameras were deployed along pathways separated 1 to 1.5 km apart [7,9,54], prioritizing local placement where pumas were more likely to be detected [29]. Cameras were mounted on the base of trees or rocks approximately 0.4-0.6 m above ground and were programmed to capture three photos when triggered, with 5-s delay between triggers. Cameras were active 24 h/day and set to record local (GMT minus 4:00) date and time when triggered. We used two camera-trap models: UNIT (Stealth Cam LLC ® Digital Scouting Cameras, Grand Prairie, TX, USA) and Trophy Cam HD (Bushnell ® Outdoor Products, Kansas, MO, USA), with trigger times of 0.8 and 0.7 s, respectively.
To increase detection probability, especially in open areas without an obvious trail, camera-stations were baited with a combination of wildcat urine and wildcat lure (e.g., Hawbakers Wild Cat Lure No. 2 and No. 1; Calvin Klein Obsession perfume [55]) placed 2-3 m from the camera in perforated plastic containers, at approximately 30 cm above ground. Understory growth in front of cameras was cleared to avoid false triggering and provide clear view of detected animals [56]. Camera-traps were in the field for 113 ± 38 (Mean ± SD) days and were inspected monthly to replace batteries and memory cards and re-bait sites.
To gain understanding of whether activity rates of each species varied seasonally or across the two study areas, we calculated 30-min capture rates (CR 30min ) and 24-h capture rates (CR 24h ) for each species, divided by the number of traps nights at each station, multiplied by 100 camera-nights. We considered the total period elapsed between the camera setting date and the last record (usually camera retrieval date) as the effective sampling period (i.e., number of trap-nights). Photographs of Diversity 2020, 12, 317 5 of 26 the same species within a 30-min interval were considered as a single event, unless animals were clearly different individuals based on sex, scars, tail color patterns, or other permanent or transitory marks [7,9,54,56]. A 24-h capture rate is considered a low biased estimator of trap rates [57], and thus we also tallied photographic events on a 24-h time-frame. Due to the difficulty in determining small mammal species (e.g., rodents, marsupials, and insectivores) from camera-trap images, we pooled independent detections of these species into the category of small mammals. We calculated CR 24h for each species or grouping and compared them between study sites and between seasons (fall-winter vs. spring-summer) using Student's t-tests to determine whether study area or weather might influence capture rates.

Single-Species, Single-Season Models, and Co-Occurrence Modeling
We used single-season occupancy modeling in Program PRESENCE [58], incorporating covariates we thought would influence distribution and detection of pumas and culpeo foxes [59,60]. It should be noted that, due to the large home range size of our carnivores relative to camera spacing, we consider occupancy to be equivalent to site use [61]. We created detection histories for both carnivores using data from the 36 camera-stations within RCNR only, because the VV&LC site was surveyed two years before RCNR, with most of the sampling effort focused on the opposite climate season. To obtain fine-scale information on detection probabilities for the carnivores, each encounter occasion was represented by three days of sampling (~10 occasions per month), yielding a total of 42 sampling occasions over the 4-month survey (May-August 2015).
We used a small number of habitat and biotic variables-due to our small sample size once we removed VV&LC stations-we thought it to be most influential in this system and likely explain carnivores site use and detection probabilities [62]: habitat (forest, shrubland, or Andean steppe), elevation, distance to the main river (Cipreses river), distance to the border of the reserve, capture rates (both CR 30min and CR 24h ) of prey (hares, rabbits, small mammals, and guanacos) and of the competing carnivore at each camera-site. We expected occupancy probability for both species to be negatively influenced by elevation because topography is extremely rugged at higher elevations and elevation has been regarded as an influential variable to carnivore spatial distribution elsewhere [63,64]. We also expected that detection would increase for both species closer to rivers because water can be scarce, and we expected carnivores to seek out water in areas they already occupied, as found elsewhere [65]. Since human disturbance (e.g., livestock grazing and tourism) increases near the edges of the reserve, we expected occupancy for pumas to increase farther away from the reserve boundary, which would limit contact with livestock ranches and tourism activities. In contrast, foxes might be more tolerant of disturbance due to their smaller size and may not be affected by distance to border, especially considering the reserve borders are somewhat "soft" containing a gradient marginal habitat. Additionally, foxes might take advantage of puma absence near the border and have higher hunting rates. In addition, we expected occupancy and detection of both species to be positively influenced by hare capture rate since hares occur in the diets of both carnivores, but that only rabbits and small mammals would positively influence fox occupancy and detection because these prey items occur in the fox diet but only have a marginal contribution to puma diet [38][39][40]66]. Finally, we did not expect puma occupancy or detection to be influenced by fox capture rate, but that fox occupancy and particularly detection could be negatively influenced by puma capture rate.
We tested variables for correlation (|r s | ≥ 0.7) and did not include correlated variables in the same model. We also standardized continuous variables to center them on zero so that the beta values were comparable [20,61]. Vegetation GIS-layers were obtained from Luebert and Pliscoff [35], and elevation data were obtained from Digital Elevation Models (DEMs) available online from the US Geological Survey (USGS) and the National Aeronautics and Space Agency (NASA) [67].
We first modeled detection while holding occupancy constant to find the best detection model, we then modeled occupancy while holding detection constant. We then combined our best detection and occupancy models to find the best model as determined by Akaike Information Criterion corrected for small sample size (AIC c ) with ∆AIC c < 2.0 denoting competing models [68,69].
To further examine potential species interactions between pumas and foxes, we fit two species co-occurrence models using the psiBa parameterization [58,61,70] within program PRESENCE [58], which requires designating a dominant and subordinate species. Due to their larger size, we considered pumas as the dominant species (A) and the smaller culpeo foxes as subordinate (B). We estimated eight parameters: Ψ P -probability that a site is occupied by puma, the dominant species, Ψ FP -the probability that area is occupied by foxes (subordinate species) given pumas are present; Ψ Fp -the probability that area is occupied by foxes given pumas are not present; p P -probability of detecting pumas given foxes are not present; p F -probability of detecting foxes given pumas are not present; r P -probability of detecting pumas given both are present; r FP -probability of detecting foxes given both are present and pumas were detected; and r Fp -probability of detecting foxes given both are present and pumas were not detected [61].
We tested eight types of different candidate models, representing different hypotheses regarding the dependence/independence of occupancy and detection parameters and between species (Table A1 in Appendix A).
We estimated the Species Interaction Factor (Φ) as follows: which measures spatial interaction between two species [70] to determine if they use habitat independently (Φ equal to 1.0 or 95% CIs overlaps 1.0), if they co-occur more than expected by chance (Φ > 1.0 and 95% CIs do not overlap 1.0), or co-occur less than expected by chance (Φ < 1.0 and 95% CIs do not overlap 1.0) [71]. We used the 'deltamethod' [61,72] in program R [73] to calculate the SIF and its standard error from the occupancy estimates using the variance-covariance matrix from the top model [61]. We also used the derived Detection Interaction Factor (δ) to evaluate if detection probability for both species was independent (δ = 1.0 or 95% CIs overlap 1.0) or not (δ 1.0 and 95% CIs do not overlap 1.0) [74]: To determine the relative importance of habitat features versus species interactions, we included the important habitat variables from the single season models for pumas and foxes in our co-occurrence models. We compared models with and without these predictor variables. If models with habitat or biotic covariates rose to the top while species interactions were deemed independent, then habitat or biotic factors are more important in influencing occupancy and detection. If species interactions without habitat rose to the top, then species interactions drive state variables. If both are included in the top model, then species interactions occur, but are mediated by habitat or biotic features.
We ranked models according to AIC c denoting competing models as those with ∆AIC c < 2.0 [68] and report top models and the constant model.

Temporal Activity Overlap Estimation
Because there is little quantitative information on timing of activity of predators and prey in central Chile, we used data from all 50 camera stations to assess strength of selection for a certain diel period. We classified all 30 min independent events into four periods of the diel cycle: dawn-from 1 h before sunrise to 1 h after sunrise; day-from 1 h after sunrise to 1 h before sunset; dusk-from 1 h before sunset to 1 h after sunset; and night-from 1 h after sunset and 1 h before sunrise [9]. Species that chose dawn and dusk were considered crepuscular, those that chose daylight hours were diurnal, nighttime hours were nocturnal, and those that occurred evenly throughout the day and night were Diversity 2020, 12, 317 7 of 26 consider cathemeral. Each event was classified considering the sunrise and sunset times on the study area for each record's date. Only species with 10 or more independent events were considered for the activity pattern diel period selectivity analysis. We calculated a 'density of detections' [9] by dividing the total number of independent detections in each diel period by the duration of that period per 100 camera nights. Species preference/strength of selection for each diel period was assessed using a modified Ivlev's selectivity index [75] proposed by Jacobs [76]; hereafter, JSI, broadly used in ecological studies to evaluate selection over many types of resources (e.g., [6,77]). JSI is formulated as follows [76]: where: D is the Jacob's Selectivity Index (JSI), u indicates the used proportion of the resource, and a indicates the available proportion of the resource. A 5000-replicates bootstrap resampling was performed and JSI was calculated for each bootstrap sample to determine the mean JSI and 95% confidence intervals [78]. We considered that each diel period was positively or negatively selected when the 95% CI was positive or negative, respectively, and did not overlap zero. JSI formula was coded in R software [73], and bootstrap resamples were run using boot package [79] in R [73]. This analysis was conducted by pooling all records for each species and each species was considered as diurnal (JSI day ≥ 0.8), strictly nocturnal (JSI night ≥ 0.8), or facultative nocturnal (0.4 ≥ JSI night > 0.8) [9], crepuscular (JSI dusk ≥ 0.8; JSI dawn ≥ 0.8). We expect both pumas, foxes, and their main prey (lagomorphs and small mammals) to show nocturnal activity, but lagomorphs could exhibit a diurnal to near-cathemeral activity pattern, as it has been reported by previous studies [8,18].
Each species' independent events are considered a random sample of the underlying continuous temporal distribution describing the probability of detecting that species within a particular interval of the day [4,80]. We used a chi-square to evaluate whether the observed differences at RCNR and VV&LC were explained by differences in sampling effort under random conditions, or not (α = 0.05). We also used a two-way, circular ANOVA [81,82] to assess significant differences (α = 0.05) in each species' activity patterns between seasons (spring-summer vs. fall-winter) or areas (RCNR vs. VV&LC). This analysis was performed using CircStats Toolbox [83] for MATLAB Release 2016a (The MathWorks, Inc., Natick, MA, USA).
We also estimated overlap in temporal activity patterns between pairs of species non-parametrically using Kernel Density Estimates (KDE) [80], only when ≥10 independent records per species were obtained [84]. We estimated the coefficient of overlap (∆) between species pairs using the ∆ 1 estimator when the smallest sample was under 50 detections, and ∆ 4 estimator otherwise [4,80,84]. ∆ ranges between 0 (no temporal overlap) to 1 (full overlap). We conducted a 10,000 sample smoothed bootstrap to generate confidence intervals for ∆ [85]. Activity pattern analysis and estimation of coefficients of overlap were performed using Overlap package [85] for R Software [73].

Scat Collection and Diet Analysis
Carnivore scats were collected between 2012 and 2015 along 3-km long transects (n = 8 to 10) per year in areas where carnivores are likely to transit. Scats found during the fieldwork outside those transects were also collected. Scat host species was identified based on size and shape, which is reliable at this study site given that wild carnivores' body sizes do not overlap and their scats have different morphology ( Figure A1 in Appendix A). Additionally, only felid scats large enough to be positively identified as pumas were collected.
We dried and de-wormed scats on a stove at 70 • C for 3 h and re-hydrated them in soapy water for at least 24 h [66]. We then washed scats and strained them on a fine mesh (1 mm) strainer and dried them again on stove for 5-10 h [66]. Hairs present in scats were analyzed under a microscope to describe Diversity 2020, 12, 317 8 of 26 the medulla-cortex and scale patterns [86] to identify prey species consumed by carnivores. Teeth in scats were identified using keys available in literature [39,87]. We excluded scats with unidentifiable remains from analysis and we pooled rodents and insectivores into a single category as small mammals. Given that most lagomorph remains in puma scats could not be identified at the species level, we pooled all the lagomorphs (i.e., rabbits and hares) into a single category.
We calculated frequency of occurrence (FO) as the number of scats positive for a prey item divided by all scats collected (with identifiable prey) and the relative percent occurrence (RPO) as the percentage of a specific prey item relative to all prey items found in scats (because there could be more than one prey item per scat) [66]. Given pumas consume large prey, we also calculated relative prey biomass (RBM) as the percentage of the consumed biomass of each prey item corrected for the total biomass consumed by the puma [66]: where BM i is the consumed biomass of the i prey species; FO i is the frequency of occurrence of the i prey species; X i is the average body mass of the i prey species, extracted from bibliographical data; n indicates the total number of prey items (i = 1, 2, 3, ..., n). We calculated RBM only for pumas because the correction factor to derive RBM from RPO is not available for foxes and, given that foxes do not consume large prey, this biomass correction is not needed [66]. To assess which species had a more generalist or specific diet, we calculated the trophic niche breadth for each carnivore, using a version of the Gladfelter-Johnson's index [88], modified to measure trophic niche breadth based on frequency occurrence, which ranges from 0 to 1, is more sensitive to resource use patterns, and is suitable to measure trophic niche breadth in a guild context [89]: where: B is the modified Gladfelter-Johnson's index; FO i is the frequency of occurrence of the i prey species; σ is the standard deviation of occurrence frequencies; and R is the number of prey-species consumed by the carnivore guild. We also calculated the Pianka's dietary niche overlap index [90] as a measure of trophic niche overlap between pumas and culpeo foxes Equation (7), given its frequent use in diet studies (e.g., [30,31,91,92]): where α is the Pianka's dietary niche overlap index between species j and species k; p i,j is the proportion of the prey item i in species j diet; p i,k is the proportion of the prey item i in species k diet; n indicates the total number of prey items (i = 1, 2, 3, ..., n).

Camera-Trap Records and Trap Rates
The total sampling effort for both study areas was 9100 camera nights from 50 camera-stations, 7106 from 36 camera-stations at RCNR, and 1994 from 14 camera-stations at VV&LC. We obtained a total of 5953 capture events of all species combined. Carnivores were captured in 3070 photos. The most frequent carnivore captured was the culpeo fox ( Figure 2), almost eight times more frequent than the puma, whereas the least frequent carnivore was the pampas cat, detected only four times and only inside the RCNR. Given the small sample size for pampas cats, this species was excluded from statistical analysis (a minimum of 10 independent records is required to conduct temporal activity overlap estimation). No species exhibited statistically significant seasonal differences in VV&LC, but in RCNR, culpeo foxes, and small mammals had significantly lower capture rates in spring summer than in winter-fall (p > 0.01). Lagomorphs were most frequently captured with hares exhibiting similar 24-h capture rates (CR 24h ) in RNRC and VV&LC (CR 24h,RCNR = 4.62, SE = 1.24, CR 24h,VV&LC = 4.42, SE = 1.78), while rabbits were more frequently captured in VV&LC (CR 24h,RCNR = 2.66, SE = 1.33; CR 24h,VV&LC = 6.85, SE = 3.30). Small mammals had similar photographic capture rates in both study sites, and were nearly absent in the spring-summer. Guanacos were the least frequently captured species and were only found in RCRN (CR 24h,RCNR = 2.18; SE = 0.81). The southern mountain viscacha, hog-nosed skunk (Conepatus chinga), and lesser grison (Galictis cuja) were not recorded during the camera-trap surveys.

Single Season Occupancy Models
Only two pairs of covariates were correlated: distance to the reserve border with elevation (rs = 0.73), and 24-h puma capture rates with 24-h hare capture rates (rs = 0.74). Therefore, we did not combine these correlated variables in the same model.
Pumas were recorded at 16 out of 36 sites (i.e., 44.0% naïve occupancy), and estimated occupancy (i.e., site use) from the top-ranked, single season model was 0.51 (SE = 0.12). The top model (with no competing models) suggests that puma site use strongly increases with hare activity (based on CR24h); however, the 95% CIs on the beta estimate overlapped zero suggesting high variability across the landscape (Tables 1 and 2, Figure 3). Puma detection also increases with hare activity (based on CR24h), and this effect was supported (95% CIs did not overlap zero) but was fairly weak. Puma detection increased with proximity to the river (m), and this effect was supported and strong, particularly at high hare capture rates (Table 1, Figure 3).

Single Season Occupancy Models
Only two pairs of covariates were correlated: distance to the reserve border with elevation (r s = 0.73), and 24-h puma capture rates with 24-h hare capture rates (r s = 0.74). Therefore, we did not combine these correlated variables in the same model.
Pumas were recorded at 16 out of 36 sites (i.e., 44.0% naïve occupancy), and estimated occupancy (i.e., site use) from the top-ranked, single season model was 0.51 (SE = 0.12). The top model (with no competing models) suggests that puma site use strongly increases with hare activity (based on CR 24h ); however, the 95% CIs on the beta estimate overlapped zero suggesting high variability across the landscape (Tables 1 and 2, Figure 3). Puma detection also increases with hare activity (based on CR 24h ), and this effect was supported (95% CIs did not overlap zero) but was fairly weak. Puma detection increased with proximity to the river (m), and this effect was supported and strong, particularly at high hare capture rates (Table 1, Figure 3).  β (95% CI) Puma Culpeo foxes were recorded at 33 out of 36 sites (91.6% naïve occupancy) and the two top competing models estimated occupancy, Ψ (95% CI), as constant at 0.920 (0.770-0.975) indicating foxes are ubiquitous across the landscape. The top model indicated strong support for a decrease in fox detection probability with increasing elevation, and the competing model, which also included elevation, indicated that fox detection decreased with increasing puma detection, but this effect was weak and not supported (Table 2, Figure 3).   Culpeo foxes were recorded at 33 out of 36 sites (91.6% naïve occupancy) and the two top competing models estimated occupancy, Ψ (95% CI), as constant at 0.920 (0.770-0.975) indicating foxes are ubiquitous across the landscape. The top model indicated strong support for a decrease in fox detection probability with increasing elevation, and the competing model, which also included elevation, indicated that fox detection decreased with increasing puma detection, but this effect was weak and not supported (Table 2, Figure 3). (a)

Two-Species Co-Occurrence Models
All but three out of the 44 candidate, 2-species co-occurrence models reached successful numerical convergence. There was only one top-ranked model with no competing models ( Table 3). As with the single season models, puma site use probability increased with hare activity, whereas fox site use probability was constant.
Fox occupancy was the same whether with pumas present or not (ΨFP(.) = ΨFp(.), Table 3), since foxes used almost every site. The Species Interaction Factor (Φ (SE)) was 1.0 (0.0) indicating independent site use by the two species. Interestingly, puma detection probabilities were greater when foxes were present (rP > pP), and increased as the hare CR24h increased (Figure 4a). Contrary to our predictions, fox detection probability was higher when pumas were present and detected at the sites (rFP > rFp > pF), and was higher at lower elevations (Figure 4b). The Detection Interaction Factor δ (95% CI) was 1.55 (1.06-2.05) indicating strong support for co-detection between the two species.

Two-Species Co-Occurrence Models
All but three out of the 44 candidate, 2-species co-occurrence models reached successful numerical convergence. There was only one top-ranked model with no competing models ( Table 3). As with the single season models, puma site use probability increased with hare activity, whereas fox site use probability was constant.
Fox occupancy was the same whether with pumas present or not (Ψ FP (.) = Ψ Fp (.), Table 3), since foxes used almost every site. The Species Interaction Factor (Φ (SE)) was 1.0 (0.0) indicating independent site use by the two species. Interestingly, puma detection probabilities were greater when foxes were present (rP > pP), and increased as the hare CR 24h increased (Figure 4a). Contrary to our predictions, fox detection probability was higher when pumas were present and detected at the sites (r FP > r Fp > p F ), and was higher at lower elevations (Figure 4b). The Detection Interaction Factor δ (95% CI) was 1.55 (1.06-2.05) indicating strong support for co-detection between the two species. r P -probability of detecting pumas, given both are present; r FP -probability of detecting foxes, given both are present and pumas were detected; and r Fp -probability of detecting foxes, given both are present, and pumas were not detected; Ψ P -probability that a site is occupied by pumas, Ψ FP -the probability that area is occupied by foxes, given pumas are present, Ψ Fp -probability that area is occupied by foxes, given pumas are not present.    AICc-Akaike Information Criterion corrected for small samples; ΔAICc-difference between each model's AICc and the top-model AICc; wi-model weight; L-model likelihood; K-number of parameters; pP-probability of detecting pumas, given foxes are not present; pF-probability of detecting foxes, given pumas are not present; rP-probability of detecting pumas, given both are present; rFP-probability of detecting foxes, given both are present and pumas were detected; and rFp-probability of detecting foxes, given both are present, and pumas were not detected; ΨPprobability that a site is occupied by pumas, ΨFP-the probability that area is occupied by foxes, given pumas are present, ΨFp-probability that area is occupied by foxes, given pumas are not present. 1 Le30-capture rate of hares from 30-min independent events per 100 trap night (CR30min); 2 Le24capture rate of hares from 24-h independent events per 100 trap night (CR24h); 3 Elev-Elevation

Diel Selectivity and Activity Pattern Overlap
We found that only pumas (F = 10.20, p = 0.002) and guanacos (F = 19.83, p = 1.6 × 10 −5 ) showed significantly higher activity during spring-summer than in fall-winter. Therefore, we plotted overlap separately for each season. No species showed significant differences in temporal activity patterns between RCNR and VV&LC, thus we combined sites.
The Jacob's Selectivity Index indicated that pumas, foxes, hares, and rabbits are facultatively nocturnal (0.4 ≥ JSI night > 0.8), whereas small mammals are strictly nocturnal (JSI night = 0.99), guanacos are near diurnal (JSI day = 0.77), and hares also show some selection for crepuscular time periods (Table 4).  The highest activity pattern overlap for pumas occurred with culpeo foxes (∆ = 0.84) in fall-winter (Figure 5a). Pumas showed a peak in activity around midnight in the fall-winter, but shifted to more activity after midnight in the spring-summer, leading to slightly lower overlap with foxes (∆ = 0.81) in that season, while foxes exhibited nocturnal activity peaks at~10:00 p.m. and~6:00 a.m. in both seasons (Figure 5a). We did not find seasonal differences for foxes-or any of their prey species-activity patterns; thus, we combined seasons in the overlap analysis for foxes and prey. The culpeo fox had an activity pattern highly synchronized with both hares (Δ = 0.89) and rabbits (Δ = 0.88) year-round and had higher synchrony with small mammals (Δ = 0.81) than pumas ( Figure 6). Foxes also showed lowest temporal overlap with guanacos (Δ = 0.71), but we note that foxes are not known to prey on guanacos. Hare activity in fall-winter occurred mostly from~8:00 p.m. to 6:00 a.m., but, interestingly, in spring-summer hares had crepuscular peaks just after sunset and before sunrise (Figure 5b). This led to higher overlap with pumas, particularly in the spring-summer (∆ = 0.82) when pumas displayed the largest peak in activity just before sunrise, matching the pre-dawn crepuscular peak in hare activity (Figure 5b). Rabbit activity in fall-winter had a single peak at~4:00 a.m. and in spring-summer exhibited multiple peaks at~8:00 a.m.,~2:00 a.m. and just before 6:00 a.m., and had similar overlap with pumas in both seasons (∆ = 0.81, Figure 5c). Small mammals displayed the highest peaks in nocturnal activity just past midnight in fall-winter and at~3:00 a.m. in spring-summer, with low relative overlap with puma activity with ∆ = 0.71 and ∆ = 0.64 in those seasons, respectively (Figure 5d). Guanacos exhibited more density of activity during night-time in the fall-winter, which led to higher overlap with pumas (∆ = 0.82) compared to the spring-summer when they had peaks of activity in morning and late afternoon resulting in lower overlap with pumas (∆ = 0.54).
We did not find seasonal differences for foxes-or any of their prey species-activity patterns; thus, we combined seasons in the overlap analysis for foxes and prey. The culpeo fox had an activity pattern highly synchronized with both hares (∆ = 0.89) and rabbits (∆ = 0.88) year-round and had higher synchrony with small mammals (∆ = 0.81) than pumas ( Figure 6). Foxes also showed lowest temporal overlap with guanacos (∆ = 0.71), but we note that foxes are not known to prey on guanacos. We did not find seasonal differences for foxes-or any of their prey species-activity patterns; thus, we combined seasons in the overlap analysis for foxes and prey. The culpeo fox had an activity pattern highly synchronized with both hares (Δ = 0.89) and rabbits (Δ = 0.88) year-round and had higher synchrony with small mammals (Δ = 0.81) than pumas ( Figure 6). Foxes also showed lowest temporal overlap with guanacos (Δ = 0.71), but we note that foxes are not known to prey on guanacos.

Diet Analysis
We found that exotic European lagomorphs made up the majority of the diet for pumas at 76% frequency occurrence (FO), 63% relative percent occurrence (RPO), and 70% by biomass (RBM), with the majority of these items being hares (Table 5). Small mammals made up the next most frequent diet item for pumas, but their biomass contribution was small (2% RBM). Although we never photographed viscachas in the study area, we did find evidence of this diet item in two puma scats. We only found ungulates occasionally in the puma diet, but they made up the second largest diet item by biomass (26%), and this was mostly driven by horse remains found in scats (5% RPO; 15% RBM) and cattle (3% RPO; 8% RBM). Other puma diet items included one culpeo fox, four birds, one reptile, but no seeds (Table 5).

Diet Analysis
We found that exotic European lagomorphs made up the majority of the diet for pumas at 76% frequency occurrence (FO), 63% relative percent occurrence (RPO), and 70% by biomass (RBM), with the majority of these items being hares (Table 5). Small mammals made up the next most frequent diet item for pumas, but their biomass contribution was small (2% RBM). Although we never photographed viscachas in the study area, we did find evidence of this diet item in two puma scats. We only found ungulates occasionally in the puma diet, but they made up the second largest diet item by biomass (26%), and this was mostly driven by horse remains found in scats (5% RPO; 15% RBM) and cattle (3% RPO; 8% RBM). Other puma diet items included one culpeo fox, four birds, one reptile, but no seeds (Table 5).
For foxes, the main prey item was small mammals at 75% FO, but the relative percent occurrence was lower at 36% (Table 5). Lagomorphs were the next largest animal prey species accounting for 25% FO and 12% RPO, of which nearly all were rabbits as opposed to hares. We did not find viscachas or guanacos, but did find a 2% RPO each for horses and cattle within fox scats. We found one scat with a hog nosed skunk, four with birds, but no reptiles in fox scats. Unlike pumas, fox scats contained seeds which made up the highest frequency of occurrence (88% FO) and the highest relative frequency of occurrence (43% RPO) of all diet items.

Carnivores Trophic Niche Breadth and Overlap
The Gladfelter-Johnson's index of trophic niche breadth was smaller for pumas (B Puma = 0.06) than for foxes (B Fox = 0.07), indicating that foxes consume a slightly broader variety of prey items.
The Pianka's dietary niche overlap index between the puma and culpeo fox was 0.30, a low to moderate value-considering both carnivores are opportunistic predators, which indicates low trophic overlap between the species.
Finally, we plotted the central Chilean food web across the three trophic niches (species interaction factor and detection interaction factor, temporal overlap ∆, and diet RPO) to visualize the strength of interactions and to assess factors leading to co-occurrence between pumas and foxes (Figure 7). moderate value-considering both carnivores are opportunistic predators, which indicates low trophic overlap between the species.
Finally, we plotted the central Chilean food web across the three trophic niches (species interaction factor and detection interaction factor, temporal overlap Δ, and diet RPO) to visualize the strength of interactions and to assess factors leading to co-occurrence between pumas and foxes ( Figure 7). Relative percent occurrence of prey in predator diet (RPO%), the temporal overlap coefficient (∆) for each species pair, and, for the two carnivores, the species interaction factor (Φ) and detection interaction factor (δ) are presented beside each arrow. The width of the arrows (black component) is proportional to the RPO% as a proxy of the intensity of the trophic interaction, whereas the width of the gray component is proportional to the coefficient of temporal overlap (Δ). Potential prey are arranged from smallest to largest body size. Although not included in the diagram, hog-nosed skunk made up 1% RPO for foxes, reptiles made up 1% RPO for pumas, and birds made up 4% and 2% RPO for foxes and pumas, respectively.

Camera-Trap Records and Trap Rates
To our knowledge, this is the first study to describe and characterize the predator-prey community in the Central Chilean Biodiversity hotspot. We provide new information on the distribution of native and exotic prey in a protected area and an area with livestock ranching, across two seasons. We evaluate factors driving carnivores coexistence by investigating three niche dimensions (space, time, and diet) simultaneously for the main carnivores, pumas and culpeo foxes. Surprisingly, we found few differences in the photographic trapping rates between the protected area (RCNR) and area with livestock grazing (VV&LC) with two notable exceptions. We captured no guanacos in VV&LC, which points to the importance of protected areas for this species, and we captured significantly fewer foxes in RNRC than VV&LC in the spring-summer season only. However, this lower fox capture rate seems to be driven by a seasonal difference within the protected area, since the winter trapping rate was similarly high across the two sites. This may be an artifact of our sampling regime because we could only effectively put cameras in valley bottoms due to steep and rugged terrain. Foxes were potentially limited to these valley bottoms due to high snow in winter thus increasing trapping rates, whereas, once snow melts, they could potentially expand and move up to higher elevations where other vegetation types occur and prey gather in spring-summer [33]. Relative percent occurrence of prey in predator diet (RPO%), the temporal overlap coefficient (∆) for each species pair, and, for the two carnivores, the species interaction factor (Φ) and detection interaction factor (δ) are presented beside each arrow. The width of the arrows (black component) is proportional to the RPO% as a proxy of the intensity of the trophic interaction, whereas the width of the gray component is proportional to the coefficient of temporal overlap (∆). Potential prey are arranged from smallest to largest body size. Although not included in the diagram, hog-nosed skunk made up 1% RPO for foxes, reptiles made up 1% RPO for pumas, and birds made up 4% and 2% RPO for foxes and pumas, respectively.

Camera-Trap Records and Trap Rates
To our knowledge, this is the first study to describe and characterize the predator-prey community in the Central Chilean Biodiversity hotspot. We provide new information on the distribution of native and exotic prey in a protected area and an area with livestock ranching, across two seasons. We evaluate factors driving carnivores coexistence by investigating three niche dimensions (space, time, and diet) simultaneously for the main carnivores, pumas and culpeo foxes. Surprisingly, we found few differences in the photographic trapping rates between the protected area (RCNR) and area with livestock grazing (VV&LC) with two notable exceptions. We captured no guanacos in VV&LC, which points to the importance of protected areas for this species, and we captured significantly fewer foxes in RNRC than VV&LC in the spring-summer season only. However, this lower fox capture rate seems to be driven by a seasonal difference within the protected area, since the winter trapping rate was similarly high across the two sites. This may be an artifact of our sampling regime because we could only effectively put cameras in valley bottoms due to steep and rugged terrain. Foxes were potentially limited to these valley bottoms due to high snow in winter thus increasing trapping rates, whereas, once snow melts, they could potentially expand and move up to higher elevations where other vegetation types occur and prey gather in spring-summer [33]. The VV&LC site has a similar elevation level throughout and we found similar fox trapping rates year-round. Alternatively, foxes might be deterred in the summer by high human use of the area for summer recreation. Small mammals also had a significantly lower trapping rate in RCNR in summer, which also could be related to the higher human use. Although pets are not allowed in the reserve, some visitors bypass regulations and domestic dogs, either as pets, rancher's working dogs, or free-ranging dogs (Muñoz and Osorio, personal observation), are present in the valley bottom during the spring-summer.

Space Use, Detection, and Overlap
Overall, we found little evidence of segregation between foxes and pumas along the spatial niche axis. In fact, evidence pointed more towards facilitation and co-detection. Foxes were recorded at nearly every camera station in RCNR in 2015, and thus it was no surprise that fox occupancy (i.e., site use) was estimated as high (0.92) and constant across the study site. Puma occupancy was lower (0.51) and increased with hare capture rate, although this effect was highly variable across the landscape with CIs overlapping 0. Thus, we found little evidence of habitat segregation between the species. However, detection rate for pumas increased with proximity to the main river, especially in areas with high hare activity, while fox detection strongly declined with elevation, indicating higher activity rates at lower elevation, as has been reported for culpeo foxes in other studies elsewhere in Chile [93] and Colombia [94].
Our two-species, co-occurrence models indicated independent occurrence across the study site with foxes having the same occupancy at sites with or without pumas present. Detection, however, revealed interesting insights that were contrary to our predictions. The detection interaction factor was strongly positive, indicating a predator was more likely to have higher detection rates at stations with the opposite carnivore present. In fact, for pumas, this was mediated by hare capture rate; puma detection increased with hare capture rate and was higher, and increased more strongly, when foxes were present. For foxes, detection interacted with puma capture rate and elevation, such that foxes had high detection rates at low elevation especially at sites with higher (vs. medium and low) puma capture rates, and fox detection declined with increasing elevation.
We predicted carnivores to be less active at sites with the other carnivore present to avoid potential aggressive encounters, yet we found the opposite, indicating that other factors (i.e., hares and distance to rivers for pumas) and habitat (elevation for foxes) were important, but so was the presence or detection rate of the other carnivore. This could indicate that the presence or activity rate of the opposite predator itself was an attractant, or that there was some other habitat/biotic variable attractive to both species that we failed to account for (e.g., see [62]). Other studies have shown co-occurrence to be common for carnivores across a wide-range of studies [14,20], but here we show that co-detection across sites, perhaps a subtler interaction, also may be common in carnivores.
We note that sample size (i.e., camera traps) was small in VV&LC, and it was surveyed in a different year; this likely limited our ability to make comparisons across sites. Given the large number of private ranches in our region, we suggest that future efforts focus on such areas to provide more insight into how predators utilize such areas across the mosaic landscape. Future research could expand trapping arrays to enlarge sample size and collect micro-habitat variables on-site as it has been done in other regions [26,62] to further explore factors responsible for positive associations in detection rates of pumas and foxes.

Temporal Activity and Overlap
We found that most species were either nocturnal (small mammals) or facultatively nocturnal (pumas, foxes, hares, rabbits), meaning they can shift diel activity according to environmental/biological cues. Guanacos were the only diurnal species, as reported previously [95] and hares showed some selection for crepuscular periods, which has been described elsewhere [96,97].
Similar to findings for the spatial niche axis, we found little evidence of temporal avoidance as a means to facilitate coexistence between pumas and foxes. In fact, temporal overlap between pumas and foxes almost always was higher than between pumas and prey. Other studies have reported temporal partitioning among carnivores. For example, felids in Sumatra with similar body sizes had low temporal overlap suggesting temporal avoidance [20], jaguars and pumas in Belize avoided using the same site simultaneously [22], and cheetahs (Acynonix jubatus) in Tanzania avoided lions (Panthera leo) and hyenas both spatially and temporally [98,99]. Although it has been reported that interference competition is likely to result in temporal partitioning [21], our lack of evidence for temporal segregation indicates that interference competition is not likely an important driver of co-existence between foxes and pumas.
Pumas exhibited more nocturnal behavior in the spring-summer, likely due to the high heat in that season [100]. In the spring-summer, pumas shifted activity to near dawn, matching the peak in hare activity and, to a lesser extent, rabbit activity. Previous studies in Bolivia and Mexico have shown that pumas did not match the daily activity patterns of any particular prey species [2,8], perhaps due to the complexity of those systems, which harbor many more prey species than our study site.
Pumas had low temporal overlap with small mammals, especially in spring-summer. We predicted pumas to have high overlap with guanacos as a potential main prey item, but instead we found similar overlap to hares in fall-winter and the lowest overlap of all species combinations in the spring-summer. Interestingly, guanacos, the only other species besides pumas that exhibited significant season differences, had more nocturnal activity in fall-winter leading to the higher overlap with pumas, but switched to being strongly diurnal in the spring-summer. Given that guanacos are a primary prey species for pumas in other areas, the lack of association between pumas and guanacos in our study is likely due to low guanaco density.
Foxes exhibited that the highest of all temporal pairwise species overlap with hares and rabbits and their overlap with small mammals was substantially higher than that of pumas. Foxes had lowest overlap with guanacos, which is to be expected considering they are not considered a main prey and that foxes were more strongly nocturnal than pumas.

Trophic, Dietary Overlap
Our diet analysis provided the strongest evidence that co-existence between pumas and foxes of the Chilean central Andes is facilitated by trophic niche axis segregation. Surprisingly, we found that puma diet was dominated by hares and that, when rabbits and hares were combined, the puma diet consisted almost entirely of these exotic lagomorphs (63% relative percent occurrence and 70% relative biomass). This is in stark contrast to other Andean areas where pumas select guanacos as their main prey item up to 88.5% of time [31,101] even though hares are available [102][103][104][105]. As guanacos are present in RNRC, we expected them to be a main prey item for pumas, yet guanaco consumption was nearly absent, only documented twice in the study area. Low guanaco density, as a result of past over-hunting and habitat change [33], may prevent pumas from hunting them in the study area.
In an area of Patagonia where guanacos are scarce (under 5 ind/km 2 ), pumas have been shown to switch prey and consume guanacos only as secondary prey items [106]. Although we have no data on guanaco abundance in our study area, the fact that they were rarely consumed indicates potentially even lower density than in Patagonia.
Interestingly, after lagomorphs, small mammals made of the next highest percent relative proportion in the puma diet, but their biomass contribution was small. Also noteworthy was that we found evidence of the native vischacha in puma diet even though we never photographed this species, either indicating there could be a few remaining on the landscape or that we did not sample the proper habitat where this species could be easier to detect, if present on the landscape. We also found that culpeo fox remains in one puma scat, suggesting that, if intraguild predation does occur, it is rare. Finally, we did find a small contribution of cattle and horse in puma scats (8% relative percent occurrence) similar to other areas [31,38], but, due to the larger body sizes of livestock, this made up a fairly high relative biomass (25%). Given the absence of livestock within RCNR, the predation or scavenging events associated to these findings should have occurred outside RCNR, either in VV&LC or elsewhere.
The majority of culpeo fox diet consisted of seeds and small-mammals (combined relative frequency of occurrence of 79%). Lagomorphs were the next largest diet component, making up 12% relative frequency of concurrency, but most of those were rabbits as opposed to hares. In contrast to our observations, studies of culpeo fox diet in other areas reveal that lagomorphs are often an important prey item [30,47,50], indicating that culpeo foxes may be exhibiting a diet shift away from hares, when pumas are present and guanacos are scarce. The high proportion of fruits and seeds present in the diet of culpeo foxes could also play a role in trophic niche segregation, as it has been previously reported for other opportunistic generalist predators that are more carbohydrate tolerant, like the pine marten (Martes martes) [107], and the red fox [108]. Trophic niche shifts appear to be a key mechanism underlying coexistence between large felids (i.e., jaguars and pumas [17,22]), but the coexistence mechanisms between the largest felid and largest canid of the Andean Mediterranean Region of Chile have not been previously explored.
Dietary niche breadth was higher for foxes than for pumas, as has been reported previously in Andean areas of Argentina [30,31]. Despite the opportunistic feeding behavior described for both pumas [109][110][111][112] and culpeo foxes [30], overlap in diet between the two carnivores was low (α = 0.30) compared to other studies reporting Pianka's dietary niche overlap indices between pumas and culpeo foxes over 0.8 (e.g., α = 0.97 at Quebrada del Condorito National Park, mountains of central Argentina [30]; α = 0.82 at northern Patagonian ranches, Neuquén Province, Argentina [31]). In addition, based on the scats in which lagomorphs were determined to the species level (hares and rabbits), we observed that pumas consume mainly hares, while foxes consume mainly rabbits. Given that we computed the Pianka's dietary niche overlap index considering rabbits and hares pooled within a single category (lagomorphs), this difference in lagomorphs consumption would decrease the observed trophic overlap to an even lower value.

Conclusions
We found high spatial overlap including co-detection, high temporal overlap, and relatively low dietary similarity (despite low prey species richness) between pumas and foxes. Rather than spatial or temporal exclusion, we found dietary segregation with exotic, introduced lagomorphs and native small mammals and seeds most likely to be the strongest factor promoting coexistence between pumas and foxes. We found evidence that (i) puma space use was positively associated with hare trapping rates, (ii) pumas shifted temporal activity seasonally to match peaks in hare activity, a phenomenon not previously observed, and (iii) lagomorphs (especially hares) made up the majority of puma diet. Foxes were ubiquitous on the landscape, yet their detection increased with puma detection, but also was influenced negatively by elevation. Foxes had more overlap with small mammals, their main mammalian prey, than pumas, but their temporal overlap with lagomorphs also was high and their diet also contained lagomorphs (mostly rabbits). Thus, we report that exotic species (hares and rabbits) are now key trophic resources in Mediterranean central Chile. This is highly relevant for wildlife management and conservation planning, given that, currently, eradication of introduced, exotic lagomorphs is the most common policy response in Chile in an effort to restore the native ecosystem. In this case, a removal plan for hares is likely to severely, negatively affect survival of the top predator. The prey base has changed dramatically in the study area since the 15th century [29] and, given how widespread guanacos were throughout South America historically, it is likely that guanacos previously were an important prey item for pumas in the study area. Thus, without guanaco population augmentation before hare removal, pumas would be unlikely to survive lagomorph removal, while foxes likely could survive on seeds and small mammals. Conservation planning and management actions should also focus on identifying and mitigating threats to guanacos in the study area, so they can naturally increase their population size.
Our work in this relatively low complexity community with introduced hares and rabbits as dominant prey species and only pumas and culpeo foxes competing for them gave us insight into carnivore coexistence and the value of exotic prey in a highly altered ecosystem. In this conservation conundrum, exotic species have become important prey items for native carnivores, and management decisions regarding the control and/or eradication of such species in the Mediterranean area of central Chile could benefit from more research on restoring native prey before large-scale, lagomorph removal campaigns.