SciELO - Scientific Electronic Library Online

vol.9 número1Spatial variation of bat diversity between three floodplain-savanna ecosystems of the Colombian LlanosPopulation density of Saguinus leucopus (Mammalia: Primates) and landscape characteristics in eastern Antioquia, Colombia índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados




Links relacionados

  • Não possue artigos similaresSimilares em SciELO



versão On-line ISSN 2007-3364

Therya vol.9 no.1 La Paz Jan. 2018 


Ecological niche and occupation by gray fox (Urocyon cinereoargenteus) at Las Barajitas Canyon, Sonora

Lucila Armenta-Méndez1  * 

Juan Pablo Gallo-Reynoso1 

Alberto Macías-Duarte2 

Marcelino Montiel-Herrera3 

Leopoldo Villarruel-Sahagún2 

1 Centro de Investigación en Alimentación y Desarrollo A. C., Laboratorio de Ecofisiología, Carretera al Varadero Nacional km. 6.6, Colonia Las Playitas, CP. 85480, Guaymas. Sonora, México. E-mail: and (JPG).

2 Universidad Estatal de Sonora, Ley Federal del Trabajo, Col. Apolo, CP. 83100, Hermosillo. Sonora, México. Email: alberto.macias@ (AMD) y (LVS).

3 Universidad de Sonora, Departamento de Medicina y Ciencias de la Salud, Edificio 7K Blvd. Luis Donaldo Colosio y Reforma Colonia Centro, CP. 83000, Hermosillo. Sonora, México. Email: (MMH).


Hutchinson defines ecological niche considering the relevance of the environment in the n-dimensional space, according to the patterns of species presence in it. For Peterson and collaborators, these patterns set the basis for niche modeling, since they directly affect the proportion of sites occupied by species. To determine the ecological variables that would favor the conservation of biodiversity at the landscape level, the ecological niche and site occupation by the gray fox (Urocyon cinereoargenteus) were simulated using photo-trapping monitoring at Las Barajitas Canyon, Sonora, Mexico. The potential habitat of the fox was modeled with the program MaxEnt, and the occupation (ψ) and probability of detection (p) were estimated with the program MARK. Vegetation was sampled using the quadrants and Braun-Blanquet methods to estimate importance values for the main plants for inclusion in the presence-absence models, as well as to know the state of the vegetation, using the species accumulation function. The sampling effort over one-year was 124 days/night-trap, obtaining 575 photographs of the gray fox. The analysis of the activity pattern revealed that the gray fox is a nocturnal-twilight species. The evaluation of the vegetation sampled evidenced the dominance of the Ruellia-Jatropha-Bursera association. The plant species accumulation function showed a positive linear relationship (r² = 0.97). In accordance with the information criteria (AICc), positive correlations were found between the sympatric species cacomixtle and coyote, and seasonality (i.e. spring and autumn), resulting in ψ = 38.5 % and p = 7.7 %. The spatial dynamics of plant communities and the relationships between animal species observed in this study indicate that, Las Barajitas Canyon is a suitable environment for the gray fox, since the Area Under the Curve (AUC) obtained was close to 1 (0.873).

Key words: Coyote; gray fox; Las Barajitas; occupation-abundance; Sonora; Urocyon cinereoargenteus


Hutchinson define el nicho ecológico incluyendo la importancia del medio ambiente en un espacio n-dimensional, en función de los patrones de presencia de las especies. Para Peterson y colaboradores, estos patrones constituyen la base para su modelación, ya que afectan directamente en la proporción de sitios ocupados por las especies. Para determinar las variables ecológicas que pudieran favorecer la conservación de la biodiversidad a nivel de paisaje, se simuló el nicho ecológico y la ocupación de la zorra gris (Urocyon cinereoargenteus) mediante un monitoreo de foto-trampeo en el Cañón de Las Barajitas, Sonora, México. Se modeló su hábitat potencial con el programa MaxEnt y se estimaron la ocupación (ψ) y la probabilidad de su detección (p) con el programa MARK. Se muestreó la vegetación por los métodos de cuadrantes y de Braun-Blanquet, para generar valores de importancia de las principales plantas e incluirlos en los modelos de presencia-ausencia, así como conocer el estado de la vegetación, usando la función de acumulación de especies. El esfuerzo de un año de muestreo fue de 124 días/noche trampa, obteniendo 575 fotografías de zorra gris. El análisis de los patrones de actividad colocó a la zorra gris como especie nocturna-crepuscular. La evaluación de la vegetación muestreada presentó una dominancia de la asociación Ruellia-Jatropha-Bursera. La función de la acumulación de especies vegetales mostró una relación lineal positiva (r² = 0.97). De acuerdo con los criterios de información (AICc), se encontraron correlaciones positivas entre especies simpátricas, el cacomixtle y el coyote y la estacionalidad, primavera y otoño, resultando en ψ = 38.5 % y p = 7.7 %. La dinámica espacial de comunidades vegetales y relaciones de fauna observadas en este estudio indican que, el Cañón de Las Barajitas es un ambiente adecuado para la zorra gris, ya que el valor del área bajo la curva (AUC) resultó cercano a 1 (0.873).

Key words: Coyote; gray fox; Las Barajitas; occupation-abundance; Sonora; Urocyon cinereoargenteus


The understanding and prediction of the responses of species to climate change at different spatial and temporal scales is a key challenge in contemporary ecology. A useful tool to assess the sensitivity of species to climate change is niche modeling, which involves determining the climate-habitat interactions that limit their distribution range (Soberon and Nakamura 2009). However, these niche models are unsuitable for characterizing demographic or ecological processes that are intrinsic to each species (Fordham et al. 2013). At the local level, the area occupied by a species (defined as the proportion of sampling units where a species is detected; Noon et al. 2012) is a parameter that has been used as an alternative to abundance. The area occupied by a species can be used to infer habitat selection and allows proposing hypothesis on species population dynamics and distribution. The occupation model (Gastón and Warren 1997; Mackenzie and Nichols 2004; Buckley and Freckleton 2010) may be useful to analyze the variables that define an ecological niche and favor the conservation of biodiversity in a given area, as some of these variables directly influence the proportion of sites occupied by a given species (Royle and Nichols 2003; Benayas 2009). In this way, occupation also assumes an n-dimensional niche based on presence-absence models.

Most studies on species distribution are based on models that assume a “global” ecological relationship (i. e., water availability, precipitation, vegetation). These ecological relationships may be affected by processes and conditions in the local habitat, and give rise to a heterogeneous pattern that is more accurately described through the values of parameters in a local model (Miller 2012). To test the above, we chose Las Barajitas Canyon, located in the arid coastal region of Sonora where anthropogenic disturbance is low. Nevertheless, due to its location and proximity to the central Gulf of California coast, this site is vulnerable to natural phenomena (tropical storms and hurricanes). Therefore, it is an ideal location to study to identify the conditions associated with ecological wildlife parameters, such as occupation, species’ activity patterns and the identification of a particular ecological niche. The importance of studying niches and occupation in this canyon lies in that these will allow relating the role and likelihood of persistence of a species in a given area (Fordham et al. 2013). In addition, it may represent a model for monitoring changes in the presence-absence of the species in the canyon related to potential natural and anthropogenic disturbances at a local scale.

Mammals are an interesting group to study given their wide distribution range and the diversity of times of the day when they carry out their activities (Bennie et al. 2014). Within mammals, the gray fox (Urocyon cinereoargenteus) is a small-to-medium-sized canid (3 to 5 kg) that is a generalist carnivore with a distribution range spanning from temperate forests to xeric shrublands. It is seemingly more abundant in sites with abundant rocky outcrops and dense vegetation cover (Fritzell and Haroldson 1982; Cypher 2003). In Mexico, it is distributed throughout the continental (Hall 1981); however, there is little information about its ecology (Servin et al. 2014). The gray fox is a species that, for being omnivorous, could function as a stabilizer of ecological networks (Fagan 1997); Due to its broad distribution and biological plasticity, it may also serve as an umbrella species (Brown 2003; Botello 2015).

By using photo trapping — a proven effective non-invasive method to construct occupation models (O’Connell et al. 2010) — the objective was to obtain the baseline conditions and requirements associated with a little disturbed habitat, which influence the ecological niche. In turn, these determine biological plasticity by modeling the distribution of this fox species in Sonora and the Sonoran Desert, as well as the occupation-abundance relationship of U. cinereoargenteus in Las Barajitas Canyon.

Materials and Methods

Study Area: Las Barajitas Canyon (28° 04’ N and -111° 09’ W) is part of a series of canyons in Sierra El Aguaje, belonging to the Sonoran Desert subdivision, in the central Gulf Coast. It is located in the southern edge of the Sonoran Desert, near Guaymas. It includes a broad transition from thorny shrubland to tropical deciduous forest (Felger 1999), which has led to a remarkable biodiversity. The local vegetation is dominated by a Bursera-Jatropha association (Shreve and Wiggins 1951). The local climate is characterized by a transition from tropical arid to tropical semi-arid conditions (Burquez et al. 1999; Peinado et al. 2011; Bogan et al. 2014). Sierra El Aguaje is located within the Cajón del Diablo hunting reserve, corresponding to the Priority Terrestrial Region RTP-18 (CONABIO 2000). Gallo-Reynoso and González-Martínez (2003) reported 799 invertebrate species (9 of which are endemic), as well as 67 mammals, 20 reptiles, 10 amphibian and 250 bird species.

An area of 5.5 km² in Las Barajitas Canyon was sampled, comprising several types of habitats that included riparian, closed xeric shrubland, caves, permanent and intermittent watering holes, open and closed streams, and areas dominated by palm trees (Figure 1).

Figure 1 Location of the study area, northwestern Mexico. Sampled area shown in green in the satellite image. Las Barajitas Canyon, Sonora. 

Photo trapping: Fourteen photo trapping stations (TS) were established separated by at least 500 m, from May 2013 to April 2014. At each TS, one camera trap (10 Wild-viewXtreme 4, two Bushnell Trophy Cam, and two Bushnell Trophy Cam HD cameras) was placed on tree trunks or wooden stakes at 50 cm above ground and facing north (Chavez et al. 2013). Given the heterogeneity of the habitat, TSs were relocated in each subsequent sampling event (season of the year) to cover most of the study area and identify the most visited sites, thereby setting the criteria for medium- or long-term sampling. In all TSs, except those located adjacent to watering holes, baits (mixture of canned sardine in tomato sauce, raisins and apples; Harrison 1997) were placed during each sampling event. Motion-activated cameras were set to take three photographs per second at 1-min intervals every time motion was detected. Photographs separated by 30 minutes of inactivity between images were considered as separate records (Srbek-Araujo and Chiarello 2005). The value of the Relative Abundance Index (RAI) that was used was the number of separate images per 100 daystrap (Monroy-Vilchis et al. 2009), and population density was estimated as RAI/area studied (Chavez et al. 2013).

Activity Patterns: We analyzed the records of U cinereoargenteus and other sympatric species such as the coyote (Canis latrans), bobcat (Lynx rufus) and cacomixtle (ringtailed cat; Bassariscus astutus), at 2-h intervals throughout the day, thus determining the frequency of occurrence (FO) for each interval. According to Monroy-Vilchis et al. (2009), for this study the twilight period was considered between 06:00 h and 08:00 h at dawn, and between 18:00 h and 20:00 h in the evening.

Vegetation Description and Plant Species-Accumulation Function: To estimate plant species richness, 14 sites close to TSs were sampled, 11 using the quadrant method (an area of 25 m × 50 m) and, on sites with very high vegetation density and cover, three were conducted through vegetation relevés (Braun-Blanquet et al. 1932). In each quadrant, the presence and abundance of each perennial plant species were recorded to calculate species diversity and density. The height of five random individuals of each species was also recorded with measuring tape to determine the vertical structure of the plant community. Photographs of species were taken for subsequent identification in the Herbarium at the University of Sonora. For sites sampled by the quadrant method, an importance value (IV) was estimated for vegetation in each quadrant, as the sum of relative density, relative frequency and dominance of the relative vertical structure, as per Wilder (pers. comm. 2016). For vegetation relives, the percentage of species cover was used (Table 1). These parameters were used as covariates for the occupation models.

Table 1 Importance Values (IV) of vegetation and percent plant coverage of relevés (Braun-Blanquet 1932) for the sites sampled. 

Quadrant IV/ % Method
1 190 Quadrant
2 147 Quadrant
3 221 Quadrant
4 189 Quadrant
5 269 Quadrant
6 205 Quadrant
7 175 Quadrant
8 209 Quadrant
9 309 Quadrant
10 210 Quadrant
11 85 Br-Bl
12 203 Quadrant
13 95 Br-Bl
14 90 Br-Bl

The vegetation was assessed through the species-accumulation function (Soberón and Llorente 1993). A matrix of the increase in the number of species was constructed with rows corresponding to cumulative species a), columns to each site sampled b), and each cell to the number of individuals of each species in each sampling site. For linear regression models, the a/b quotient was calculated, where the response variable is obtained by adding a new species (Soberón and Llorente 1993). These species-accumulation functions based on an appropriate model allow the prediction of species richness (Brown 2003). In this study, this function was used for integrating two vegetation sampling methods (quadrants and Br-Bl), thereby determining the state of the vegetation at a given time, and hence discuss the influence of this variable on the presence-absence of the gray fox. Data were analyzed with the program Sigma-Plot version 11.0 (Systat Software Inc., Chicago).

Modeling of Occupation: The occupation (ψ) and probability of detection (p) models were designed by selecting categorical explanatory variables and quantitative variables likely associated with the presence of the gray fox. The categorical variables selected were TS and season of the year (SY, with spring, summer, autumn and winter as categories). The quantitative variables selected were VI (vegetation), presence of sympatric species and height above mean sea level of each TS. Capture histories by TS were produced based on photographic records of U. cinereoargenteus. According to Macias-Duarte and Panjabi 2013, seven occupation models were produced (log-linear relationship model), the parameters of which were estimated with the program MARK (White et al. 2000). The best model to infer the key factors for the detection of and occupation by U. cinereoargenteus was selected using the Akaike (1987) criterion for small samples (AICc; Burnham and Anderson 2002). The statistical significance of the occupation model was determined by a Chi-square test with 95% confidence (White et al. 2000). The occupation-abundance relationship was explained applying a linear regression model with occupation as the TSs where the gray fox was captured, and abundance as the independent records of the gray fox in each TS by season.

Ecological Niche Modeling: The ecological niche and the potential distribution of the gray fox were modelled using records of the presence of this species from the “Cuatro Gatos” project of the Sky Island Alliance non-governmental organization, and from the Global Biodiversity Information Facility (, and considering the environmental variables associated with these points, cut in ArcGIS 10.1, according to the distribution of this species (Phillips et al. 2006). The bioclimatic variables were obtained from World-clim ( The data were analyzed with the program MaxEnt, which is based on a maximum-entropy algorithm (a higher probability of presence represents less environmental restrictions; Petterson et al. 2011). The bases of the species distribution models (or their ecological niche) are cartographic representations of mathematical or statistical relationships between the actual known distribution and a set of independent variables within a space for the presence of a species, depending on the variables used to derive such distribution. This value represents a probability of presence of the species (Mateo et al. 2011).


Camera trapping: The sampling effort was 164 days/camera. It captured a total of 7,459 photographs of fauna, with 13% corresponding to the gray fox and the rest to sympatric species including coyote (16%), ring-tailed cat (6%) and bobcat (2%), as well as other animal species (Procyon lotor, Tamias dorsalis, Tayassu tajacu, Odocoileus virginianus, Macrotus californicus, Neotoma albigula, Otospemorphilus variegatus, Nasua narica, Lepus californicus, Sylvilagus audobonii and several bird species). Ninety five photographs were considered as individual records of gray fox, with 2 including 2 individuals (female and offspring; female and male), resulting in a RAI = 0.95 for subtropical xeric schrubland. By applying the formula of Chavez et al. (2013), the density of gray foxes is 1.73 individuals per km² in the study area. The highest number of detections occurred in the spring. The sites with the highest Frequency of Occurrence (FO) were those where vegetation cover was >75% according to the method of Braun-Blanquet (1932). Activity Patterns: The activity patterns of the gray fox revealed a higher FO during the night and at twilight hours. The coyote is most active throughout the day, but both this species and the gray fox were found to prefer the darkest hours. In addition, the level of activity of the coyote and the gray fox remained constant throughout the year, although it decreases in the summer (Figure 2).

Figure 2 Activity patterns of the gray fox (Urocyon cinereoargenteus) and other sympatric species. L. rufus (red), B. astutus (orange), U. cinereoargenteus (ligth green), and C. latrans (dark green). 

Vegetation Description: A total of 3,233 individuals were recorded in 14 sampling sites, belonging to 59 different species distributed in 32 families. The family with the highest species richness was Fabaceae (10), but families with the highest frequency were Euphorbiaceae, Burseraceae and Acanthaceae, being present in all sampling sites. This assessment of the flora of Las Barajitas Canyon found that the dominant plant association is Ruellia-Jatropha-Bursera. With regard to species-accumulation function, a positive linear relationship was observed between species richness and number of sampling sites. The number of sites sampled was insufficient for the function to reach an asymptote of its distribution (Figure 3a). The sites sampled by the Braun-Blanquet method revealed a more pronounced convergence between desert and subtropical species (Table 2).

Figure 3 A) Species-accumulation function with the vegetation sampled in Las Barajitas Canyon, Sonora, Mexico. B) Association between occupation and abundance of the gray fox (Urocyon cinereoargenteus) (r² = 0.88) in Las Barajitas Canyon, Sonora, Mexico, in photo trapping stations in May 2013-August 2014. 

Table 2 Evaluation of the plant community with higher Frequency of Occurrence of gray fox, by the method of Braun-Blanquet (1932)

Species Family Index Trait
Watering hole
Washingtonia robusta Arecaceae 3 characteristic sp
Bonellia macrocarpa Primulaceae 1 preferential sp.
Carnegia gigantea Cactaceae 2 selective sp.
Pachycereus pringlei Cactaceae 1 preferential sp.
Stenocereus thurberi Cactaceae 1 preferential sp.
Bursera microphylla Burseraceae 2 selective sp.
Lycium sp. Solanaceae P associated sp.
Maytenus phyllanthoides Celastraceae R incidental sp.
Justicia californica Acanthaceae P associated sp.
Macroptilium atropurpureum Fabaceae P associated sp.
Gossypium davidsonii Malvaceae P associated sp.
Hibiscus denudatus Malvaceae R incidental sp.
Hofmeisteria crassifolia Asteraceae P associated sp.
Bebbia juncea Asteraceae P associated sp.
Tephrosia vicioides Fabaceae 1 preferential sp.
Stegnosperma halimifolium Stegnospermataceae P associated sp.
False mangrove
Maytenus phyllanthoides Celastraceae 3 characteristic sp.
Simmondsia chinensis Simmondsiaceae 1 preferential sp.
Jatropha cuneata Euphorbiaceae 2 selective sp.
Lycium sp. Solanaceae 1 preferential sp.
Brahea brandegeei Arecaceae P associated sp.
Atriplex barclayana Amaranthaceae P associated sp.
Forchhammeria watsonii Resedaceae R incidental sp.
Stenocereus thurberi Cactaceae R incidental sp.
Prosopis glandulosa vartorreyana Fabaceae 2 selective sp.
Echinocereus scopulorum Cactaceae R incidental sp.
Dense shrubland
Randia thurberi Rubiaceae 1 preferential sp.
Simmondsia chinensis Simmondsiaceae 2 selective sp.
Maytenus phyllanthoides Celastraceae 1 preferential sp.
Forchammeria watsoni Resedaceae 1 preferential sp.
Bonellia macrocarpa Primulaceae 1 preferential sp.
Lycium sp. Solanaceae 1 preferential sp.
Olneya tesota Fabaceae 1 preferential sp.
Krameria sonorae Krameriaceae 1 preferential sp.
Ruellia californica Acanthaceae 2 selective sp.
Sideroxylon occidentale Sapotaceae 3 characteristic sp.
Sapindus saponaria Sapindaceae 3 characteristic sp.
Marsdenia edulis Apocynaceae 1 preferential sp.
Prosopis glandulosa vartorreyana Fabaceae 3 characteristic sp.
Ziziphus amole Rhamnaceae 3 characteristic sp.
Phaulothamnus spinescens Achatocarpaceae 3 characteristic sp.
Ziziphus obtusifolia Rhamnaceae 2 selective sp.

Index Values for species:

3 = > 75 %, 2 = 50 to 75 %, 1 = 25 to 50 %, P = 5 % and R = < 5 %.

Occupation: The best model (lowest ΔAICc) for occupation (ψ) and probability of detection (p) of the gray fox included the variables climate season and sympatric species (Tables 3 and 4), with ψ = 38.5% (DT = 1.02) p = 7.7% (DT = 0.36). There is a positive linear correlation between occupation and abundance (R² = 0.88, Figure 3b).

Table 3 Models classified by AICc for occupation by gray fox. 

Model AICc ΔAICc Wi
p (station+sympatric species)ψ(station+sympatric species) 447.77 0 0.52
p(station)ψ(station) 447.99 0.22 0.47
p (sympatric species)ψ(sympatric species) 477.21 29.44 0
p(station+IV+alt) ψ(station+IV+alt) 486.55 38.77 0
p(station+IV+alt+caco+lynx+coyote) ψ(station+IV+alt+caco+lynx+coyote) 492.51 44.73 0

p = probability of detection, ψ= occupation. IV= importance value of vegetation, alt= altitude above sea level of the photo trapping station. Sympatric species= coyote, bobcat, and cacomixtle.

Table 4 Best model for the occupation of gray fox (lower ΔAICc). 

Parameter β SE Lower Limit
of 95% CI
Upper Limit
of 95% CI
Autumn 1.9725 1.4729 -0.9143 4.8595
Winter -1.0615 1.3196 -3.6479 1.5249
Spring 1.4667 1.2979 -1.0771 4.0105
Cacomixtle 1.6695 0.2970 1.0874 2.2515
Bobcat -0.0488 0.3449 -0.7249 0.6272
Coyote 1.6142 0.3351 0.9573 2.2711

Estimators of regression coefficient (β), standard error (SE) and confidence interval (95% CI).

Ecological Niche Modeling: The models of distribution and probability of occurrence produced by MaxEnt for the State of Sonora and the Sonoran Desert indicate that the probability of presence of U. cinereoargenteus in the study area ranges from 0.46 to 0.62 (Figures 4a, 4b). The bioclimatic variables of highest importance to generate the models for Sonora were precipitation in the coldest trimester, medium diurnal range, and mean temperature in the driest trimester; for the Sonoran Desert, seasonal precipitation, minimum temperature of the coldest season, and precipitation in the warmest trimester, with an AUC = 0.873.

Figure 4 Potential distribution model in Sonora, Mexico, with MaxEnt for the gray fox (Urocyon cinereoargenteus). A) More intense colors (red) indicate a higher probability of presence. B) More intense colors (red) indicate a higher probability of presence. Presence spots of gray fox were used for modeling. 


Camera trapping: Under the assumption that all TSs have the same probability of capture (Chavez et al. 2013), U cinereoargenteus could be anywhere within the habitat at Las Barajitas Canyon. It was found that TSs with the highest FO were the sites with the most dense vegetation cover, which is considered to be protective of daytime movements across the range of the gray fox (Haroldson and Fritzell 1984; Farias et al. 2005; Sollmann et al. 2013). Fritzell and Haroldson (1982), considering different locations, seasons and estimation methods, found that the density of the gray fox varied from 1.2 to 2.1 individuals per km², which is consistent with the density found in this study (1.73). The value of this parameter indicates that the habitat allows maintaining stable gray fox populations.

Activity Patterns: The coyote (8 to 20 kg), bobcat (5 to 15 kg), and gray fox (3 to 5 kg; Fedriani et al. 2000) are sympatric species across most of North America (Nicholson and Hill 1981; Farias et al. 2005). This convergence is observed in Las Barajitas Canyon, where cacomixtle (0.87 to 1.1 kg) adds to the above species (Poglayen-Neuwall and Toweill 1988). Sargeant et al. (1987) suggest that the sympatric relationship between gray fox and coyote occurs due to the semi-arboreal habit of the former, which provides an effective mechanism to escape from predators. This feature may be considered as a response to pressure from predation by larger species such as bobcat and coyote. Our results differ from those reported by Farias et al. (2005), who mention that the gray fox avoids areas of activity of the potential predators, in both time and space. Similar to the coyote, however, the bobcat is a potential predator of this species. By observing that these species display the same activity patterns in twilight and evening hours, their presence at similar times of the day may be associated with a lower risk of predation. These findings suggest a particular detection in sites frequently visited, where individuals can hide, climb trees or move fast enough to avoid confrontations with these predators.

Vegetation Description: The regularity of flowering and fruiting of plant species is key for sustaining the fauna in a given locality (León de la Luz et al. 1996). Las Barajitas Canyon is a valuable site for the convergence of desert and subtropical species, which provide spaces with greater coverage throughout the year. This is a safer habitat for the gray fox, in accordance with the sites sampled and the frequency of occurrence of this species in these sites. Thus, in terms of vegetation structure in the niche of the gray fox, this canyon is a stable ecosystem, given the relationship between vegetation and fauna to maintain their processes.

Occupation: This depends on the effective or realized niche, which fosters a successful adaptation (Soberon and Nakamura 2009). For an area of 5.5 km² (our field study), occupation values can be considered as high, given the climatic conditions in Las Barajitas Canyon. The model selected to determine and model the presence of U. cinereoargenteus indicates a negative influence in summer, and that the abundance of this species is not directly related to plant diversity. However, floral phenology triggers biological activity and interactions (frugivory), with deciduous shrubs providing a protective dense cover that offers a suitable space for survival.

The biotic conditions associated with climate seasonality are likely to provide the gray fox with a suitable habitat to settle as a closed population, according to the occupation model (Gastón and Warren 1997). This has led to determine the occupation of the habitat and the ecological niche, as measured with a log-linear model, making possible to identify significant effects between resources and environmental conditions.

Ecological Niche: The maps produced with MaxEnt for the study area are within an acceptable value range for classification (Mateo et al. 2011), since U cinereoargenteus is a generalist species with a preference for forested habitats (Fritzell and Haroldson 1982; Chamberlain and Leopold 2000).

The ecological niche of U. cinereoargenteus in Las Barajitas Canyon, as described based on occupation models and correlated with categorical and quantitative habitat variables, indicate that this species does not show niche overlap in a virtually undisturbed habitat (Gastón and Warren 1997). It would therefore be appropriate to consider that the presence of omnivorous generalist species is key for the self-regulation of natural systems (Fagan 1997; Brown 2003; Botello 2015). The ecological plasticity of the gray fox (Lord 1961), coupled with the heterogeneity and floristic diversity of its habitat, allow the diversification of niches between species, contributing to a lower competition for resources; as a result, biodiversity is favored at the local level.

As the existence of an inverse relationship between the density of gray fox and coyote has been previously suggested (Farias et al. 2005). Sollman et al. (2012) state that species with similar feeding habits show a division of habitat and niche, as well as differences in their respective activity patterns, hence facilitating a stable coexistence. In Las Barajitas Canyon, under conditions of little disturbance, the gray fox and the coyote seemingly do not show niche or occupation overlap, and their populations coexist. The analytical techniques used in this study were suitable to identify local factors that may either favor or restrain the particular occupation and ecological niche of the gray fox. Hence, landscape complexity provides environments that favor species and their interspecific relationships, leading to ecological niches based on the local landscape factors in this canyon. For omnivorous species such as the gray fox, the season of the year is seemingly not associated with a drastic change in niche and occupation, suggesting that that this species serves as a major stabilization factor in ecological networks (Fagan 1997).


To M. Mackenzie, for allowing access to Las Barajitas and for the enthusiasm in its conservation. To R. Felix, for the unconditional support in the edition of maps in this work. To Sky Island Alliance, for providing records of the gray fox in northern Sonora. To the Herbarium at University of Sonora, for its important archive and dedication to the native flora. To R. Lara, for the critical review of the manuscript. Last, to the responsible citizens that pay their taxes, allowing CONACYT to grant scholarships to conduct research in Mexico. María Elena Sánchez-Salazar translated the manuscript into English.

Literature cited

Akaike, H.. 1987. Factor analysis and AIC. Psychometrika 52:317-332. [ Links ]

Benayas, J. M. R.. 2009. La rareza de las especies. Investigación y ciencia 392:63. [ Links ]

Bennie, J. J., J. P., Duffy, R., Inger, and K. J., Gaston. 2014. Biogeography of time partitioning in mammals. Proceedings of the National Academy of Sciences 111:13727-13732. [ Links ]

Bogan, M. T., N., Noriega-Félix, S. L., Vidal-Aguilar, L. T., Findley, D. A., Lytle, O. G., Gutiérrez-Ruacho, and A., Varela-Romero. 2014. Biogeography and conservation of aquatic fauna in spring-fed tropical canyons of the southern Sonoran Desert, México. Biodiversity and conservation 23:2705-2748. [ Links ]

Botello, F., V., Sánchez-Cordero, and M. A., Ortega-Huerta. 2015. Disponibilidad de hábitats adecuados para especies de mamíferos a escalas regionales (estado de Guerrero) y nacional (México). Revista Mexicana de Biodiversidad 86:226-237. [ Links ]

Buckley, H. L., and R. P., Freckleton. 2010. Understanding the role of species dynamics in abundance–occupancy relationships. Journal of Ecology 98:645-658. [ Links ]

Búrquez, A., A., Martínez-yrízar, R.S., Felger, and D., Yetman. 1999. Vegetation and Habitat Diversity at the Southern Edge of the Sonoran Desert. Ecology of Sonoran Desert Plants and Plant Communities (Robichaux R. H. ed). University of Arizona Press. Tucson, U. S. A. [ Links ]

Braun-Blanquet, J., Fuller, G., and H., Conrad. 1932. Plant sociology: the study of plant communities. Mc Graw Hill Book Co. NewYork, U. S. A. [ Links ]

Brown, J. H.. 2003. Macroecología, Primera edición. Fondo de cultura económica. Ciudad de México, México. [ Links ]

Burnham, K. P., and D. R., Anderson. 2002. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2 ed. Springer-Verlang, New York, U. S. A. [ Links ]

Conabio. 2000. Listado de regiones terrestres prioritarias de la Comisión Nacional para uso de la Biodiversidad (2000) acceso en: ]

Cypher, B. L.. 2003. Foxes. Wild Mammals of North America: Biology, Management, and Conservation. Segunda edición. The Johns Hopkins University Press, Baltimore, U. S. A. [ Links ]

Chamberlain, M. J., and B. D, Leopold. 2000. Spatial use patterns, seasonal habitat selection, and interactions among adult gray foxes in Mississippi. The Journal of wildlife management 742-751. [ Links ]

Chávez, C., A, de la Torre, H., Bárcenas, R. A., Medellín, H., Zarza, and G., Ceballos. 2013. Manual de fototrampeo para estudio de fauna silvestre. El jaguar en México como estudio de caso. Alianza WWF-Telcel, Universidad Nacional Autónoma de México. Ciudad de México, México. [ Links ]

Farías, V., T. K., Fuller, R. K., Wayne, and R. M, Sauvajot. 2005. Survival and cause specific mortality of gray foxes (Urocyon cinereoargenteus) in southern California. Journal of Zoology 266:249-254. [ Links ]

Felger, S. R.. 1999. The Flora of Cañón de Nacapule: A Desert Bounded Tropical Canyon near Guaymas, Sonora. México. Proceedings of the San Diego Society of Natural History. San Diego, U. S. A. [ Links ]

Fordham, D. A., H. R., Akcakaya, M. B., Araújo, D. A., Keith, and B. W., Brook. 2013. Tools for integrating range change, extinction risk and climate change information into conservation management. Ecography 36:956-964. [ Links ]

Fagan, W. F.. 1997. Omnivory as a stabilizing feature of natural communities. The American Naturalist 150:554-567. [ Links ]

Fedriani, J. M., T. K., Fuller, R. M., Sauvajot, and E. C., York. 2000. Competition and intraguild predation among three sympatric carnivores. Oecologia 125:258-270. [ Links ]

Fritzell, E. K., and K. J., Haroldson. 1982. Urocyon cinereoargenteus. Mammalian species 189:1-6 [ Links ]

Gallo-Reynoso, J. P., and A., González-Martínez. 2003. Estudio previo justificativo para proponer el establecimiento de la región Sierra del Aguaje, Bahía de San Francisco e Isla San Pedro Nolasco y sus aguas aledañas como una nueva Área Natural Protegida. CIAD, A.C., Grupo de Estudios Ambientales, A. C. Guaymas, México. [ Links ]

Gaston, K. J., and P. H, Warren. 1997. Interspecific abundance-occupancy relationships and the effects of disturbance: a test using microcosms. Oecologia 112:112-117. [ Links ]

Hall, R. E.. 1981. The mammals of North America. John Wiley and Sons. New York, U. S. A. [ Links ]

Haroldson, K. J., and E. K., Fritzell. 1984. Home ranges, activity, and habitat use by gray foxes in an oak-hickory forest. The Journal of wildlife management 48:222-227. [ Links ]

Harrison, R. L.. 1997. A comparison of gray fox ecology between residential and undeveloped rural landscapes. The Journal of wildlife management 61:112-122. [ Links ]

Hutchinson, G. E.. 1957. Cold spring harbor symposium on quantitative biology. Concluding remarks 22:415-427. [ Links ]

León de la Luz, J. L., L. R. C., Benet, and M. C., Estrada. 1996. Fenología floral de una comunidad árido-tropical de Baja California Sur, México. Acta Botánica Mexicana 35: 45-64. [ Links ]

Lord, Jr. R. D.. 1961. A population study of the gray fox. American Midland Naturalist 66:87-109. [ Links ]

Macías-Duarte, A., and A., Panjabi. 2013. Association of habitat characteristics with winter survival of a declining grassland bird in Chihuahuan Desert grasslands of Mexico. The Auk 130:141-149. [ Links ]

Mateo, R. G., Á. M., Felicísimo., and J., Muñoz. 2011. Modelos de distribución de especies: Una revisión sintética. Revista chilena de Historia Natural 84:217-240. [ Links ]

MacKenzie, D. I., and J. D., Nichols. 2004. Occupancy as a surrogate for abundance estimation. Animal biodiversity and conservation 27:461-467. [ Links ]

Miller, J. A.. 2012. Species distribution models: Spatial autocorrelation and non-stationarity. Progress in Physical Geography 36:681-692. [ Links ]

Monroy-Vilchis, O., C., Rodríguez-Soto, M., Zarco-González, and V., Urios. 2009. Cougar and jaguar habitat use and activity patterns in central Mexico. Animal Biology 59:145-157 [ Links ]

Moreno, C. E.. 2001. Métodos para medir la biodiversidad. M&T–Manuales y Tesis SEA 1:84 [ Links ]

Nicholson, W. S., and E., Phill. 1981. An ecological study of the gray fox in east-central Alabama. Alabama Cooperative Wildlife Project Final Report. Project W-44-5. CIUDAD, U. S. A. [ Links ]

Noon, B. R., L. L., Bailey, T. D., Sisk, and K. S., McKelvey. 2012. Efficient Species- Level Monitoring at the Landscape Scale. Conservation Biology 26:432-441. [ Links ]

O’Connell, A. F., J. D., Nichols, and K. U., Karanth (eds.). 2010. Camera traps in animal ecology: methods and analyses. Springer Science and Business Media. Tokyo, Japón. [ Links ]

Peinado, M., M. Á., Macías, F. M., Ocaña-Peinado, J. L., Aguirre, and J., Delgadillo. 2011. Bioclimates and vegetation along the Pacific basin of Northwestern Mexico. Plant ecology 212:263-281. [ Links ]

Peterson, A. T., J., Soberón, R. G., Pearson, R. P., Anderson, E., Martínez-Meyer, M., Nakamura, and M.B., Araújo. 2011. Ecological niches and geographic distributions. Princeton, New Jersey: Princeton University Press. Monographs in Population Biology. New Jersey, U. S. A. [ Links ]

Phillips, S. J., Anderson, R. P., and R. E., Schapir. 2006. Maximum entropy modeling of species geographic distributions. Ecological modelling 190:231-259. [ Links ]

Poglayen-Neuwall, I., and D. E., Toweill. 1988. Bassariscus astutus. Mammalian Species 327:1-8. [ Links ]

Royle, J. A., and J. D., Nichols. 2003. Estimating abundance from repeated presence–absence data or point counts. Ecology 84: 777–790. [ Links ]

Servín, J., A., Bejarano, N., Alonso-Pérez, and E., Chacón. 2014. El tamaño del ámbito hogareño y el uso de hábitat de la zorra gris (Urocyon cinereoargenteus) en un bosque templado de Durango, México. Therya 5:257-269. [ Links ]

Sargeant, A. B., S. H., Allen, and J. O., Hastings. 1987. Spatial relations between sympatric coyotes and red foxes in North Dakota. The Journal of Wildlife Management 51:285-293. [ Links ]

Soberon, J., and M., Nakamura. 2009. Niches and distributional areas: concepts, methods, and assumptions. Proceedings of the National Academy of Sciences106:19644-19650. [ Links ]

Soberón, J., and Llorente, B. J. 1993. The Use of Species Accumulation Functions for the Prediction of Species Richness. Conservation Biology 7:480–488. [ Links ]

Sollmann, R., M., Furtado, M. H., Hofer, A. T., Jácomo, N. M., Tórres, and L., Silveira. 2012. Using occupancy models to investigate space partitioning between two sympatric large predators, the jaguar and puma in central Brazil. Mammalian Biology-ZeitschriftfürSäugetierkunde 77:41-46. [ Links ]

Sollmann, R., A., Mohamed, H., Samejima, and A., Wilting. 2013. Risky business or simple solution–Relative abundance indices from camera-trapping. Biological Conservation 159:405-412. [ Links ]

Srbek-Araujo, A. C., and A. G, Chiarello. 2005. Is camera-trapping an efficient method for surveying mammals in Neotropical forests? A case study in south-eastern Brazil. Journal of Tropical Ecology21:121-125. [ Links ]

Shreve, F., and L., Wiggins. 1951. Vegetation and flora of the Sonoran Desert. Vegetation of the Sonoran Desert. Carnegie Institution of Washington. Washington, U. S. A [ Links ]

White, G. C., K., Burnham, and D., Anderson. 2000. Advanced features of program MARK. In Wildlife, land, and people: priorities for the 21st century. Proceedings of the second international wildlife management congress. The Wildlife Society. Bethesda, U. S. A. [ Links ]

Associated editor: Rafael Avila

Received: September 14, 2017; Accepted: December 22, 2017

* Corresponding author: Lucila Armenta-Méndez, e-mail:

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License