Introduction
Biological corridors are important to connect populations that would otherwise be isolated, and thus mitigate the effects of fragmentation, one of the main threats to biodiversity (Bennett et al. 2004). Corridors are generally identified from spatial analyses in combination with inferences about movement patterns (i. e. "functional corridors") and habitat requirements (i. e. "structural corridors"; Cushman et al. 2009).
Large carnivores often function as keystone species in terms of their value for wildlife conservation. When extirpated, ecosystem functions may be destabilized and ecosystems could be driven to a disturbance state (Cramer and Portier 2001). The jaguar (Panthera onca) is an apex predator with a historical distribution that ranged from the southern United States to the Río Negro in Argentina and included much of México (Seymour 1989). Presently, the distribution of the jaguar in México is much reduced as compared to its historic range (Sanderson et al. 2002). Moreover, due to its large home range, this species is susceptible to the influences of habitat fragmentation (Dixon et al. 2006). Therefore, it is considered endangered in México (SEMARNAT 2010).
The Sierra Madre Oriental (SMO), the eastern mountain range of México, has been classified as having only medium priority for the long-term conservation and viability of the jaguar in México (Sanderson et al. 2002). This classification was given because the SMO was considered isolated from the rest of the distributional range of jaguars in México (Rabinowitz and Zeller 2010). However, recent studies indicate that the SMO may have significant potential for the conservation of jaguars (Villordo-Galván et al. 2010; Rodríguez-Soto et al. 2011). Thus, characterizing the SMO as a region of limited conservation potential for jaguars may be premature, and could result in well-intentioned conservation priorities that may actually be detrimental for the species, as has happened to other organisms (e. g.Morrison et al. 2012). For example, we are aware of no evidence that indicates that jaguar populations in the SMO are genetically isolated from other populations, which could justify marginalization of the SMO as a viable conservation unit for jaguars in México. Because of the presence of jaguar populations in the SMO, we suggest that the SMO deserves attention to evaluate both the real status of jaguar habitat and the connectivity among populations in this region.
Distributional studies have reported jaguars in areas where they were considered extirpated, suggesting that these areas may contain habitat that could sustain the species (Sanderson et al. 2002). These studies have also identified the need to promote connectivity among these areas to increase both dispersal potential and maintenance of these populations (Grigione et al. 2009; Villordo-Galván et al. 2010; Rodríguez-Soto et al. 2011; Briones-Salas et al. 2012). This is the present case in the SMO of México. Therefore, our objectives were to (1) evaluate the current habitat conditions for jaguars within the SMO and (2) to model habitat linkages that could facilitate exchange of jaguars with established populations, particularly to the south of the SMO.
Material and Methods
Study area. The flora and fauna of the SMO are mainly of tropical origin, and the SMO in northeastern México is a key region for biodiversity conservation (Ortega-Huerta and Peterson 2004). It is located at the northern limits of low and medium tropical moist forests, as well as tropical oak (Quercus spp.) forest. Cloud forests are also present in this region. Dominant climates include warm sub-humid (AW2), sub-temperate humid (A) Cf, and temperate with continuous precipitations C(f). The topography of this area range from very rugged to flat to undulating, with a of humid tropical climate. The altitude ranges from 0 to 2,300 m and the annual rainfall varies from 600 to 2,500mm (INEGI 2011).
Our evaluation area was centered at 21° 4´ 46´´ N and -98° 47´ 8´´ W and included 30,837 km2 in the southern SMO. This region includes seven priority conservation areas that are managed under different categories (Ortega-Huerta and Peterson 2004). These areas include the Abra-Tanchipa Biosphere Reserve which harbors one of the most important jaguar populations in the SMO (Figure 1; Villordo-Galván et al. 2010).
Increasing human populations and resultant habitat loss and increased fragmentation of remaining natural habitats is the primary threat for long term conservation of the SMO. In addition to increasing urbanization, rural areas are influenced by frequent land conversions from natural habitats to agriculture and cattle pasturing. In all cases, road development threatens the connectivity of wildlife populations and natural habitats throughout the mountain range. Numerous rural agrarian communities are found throughout the SMO, and these communities require services, environmental products, and food as they grow; this often results in overuse of the natural resources (including water, timber, forage for livestock, game, etc.) in the region. Current policies affecting habitat conservation, restoration, and management are mostly inadequate in the region, as are policies aimed at conservation of biodiversity (Villordo-Galván et al. 2010; Sahagún-Sánchez et al. 2011).
Identification of potential habitat patches. We developed a jaguar occurrence database for the SMO that included only photographic records that (1) were collected within the previous 12 years and (2) were verified as jaguars by experts. We used maximum entropy modeling (Maxent 3.3.3; Phillips et al. 2006) to predict potential jaguar habitat patches from presence records. Maximum entropy modeling uses only known locations to determine key correlates of a species presence in a region. Because it uses only presence data, Maxent eliminates the need for pseudo-absence data which can bias wildlife-habitat models and thus provides a less biased alternative to other approaches that require the generation of known non-use areas (i. e. discriminant analysis, logistic regression, etc; Phillips et al. 2006; Baldwin and Bender 2008). We identified variables that likely influence jaguar presence in Mexican habitats, including elevation, slope, land use and vegetation type, distance to roads, distance to water sources, and human population density (Ortega-Huerta and Medley 1999; Monroy-Vilchis et al. 2008; Rosas-Rosas et al. 2010; Colchero et al. 2011; Rodríguez-Soto et al. 2011). We obtained baseline geospatial digital maps for these variables from INEGI (r) (Mexican National Institute of Statistics and Geography). Digital map resolution was 30 m at 1:250,000.
We modeled all possible 1-6 variable candidate models and compared resultant models using receiver operating characteristic (ROC) plots and the critical ratio test (Pearce and Ferrier 2000 as modified by Baldwin and Bender 2008). We compared models using ROC plots to assess relative performance and to establish thresholds for identifying the likelihood of a site being occupied by jaguar (Phillips et al. 2006). We then used the area-under-the-curve (AUC) to assist in selecting the most-appropriate model (Fielding and Bell 1997; Phillips et al. 2006). AUC provides an index of accuracy of models; values range from 0.5 to 1.0 with values of 0.5 indicating no fit greater than that expected by chance alone, whereas models with AUC > 0.7 indicate good fit (Swets 1988). We used 30 % of locations as test data to approximate SEs around AUC scores.
We compared all possible models, and reported models with the highest AUC value for each subset of habitat variables (i. e. 1-6 variable models). Because higher dimensioned models often have greater AUC even if some variables contribute little to the model (Baldwin and Bender 2008), we (1) used AUC to select the model with the best fit (highest AUC) from each set of 1-6 variable models (to ensure that only models that fit the data well were considered) and then (2) compared whether the highest AUC model differed statistically from lower (more parsimonious) models. If they did not differ, we selected the most parsimonious (fewer variable) model. For the latter, we used the critical ratio test (Pearce and Ferrier 2000; Baldwin and Bender 2008) to compare the highest AUC model to simpler models to determine if the increase in explanatory value was significant at α = 0.05. We also derived thresholds for probability of jaguar presence for test data by maximizing sensitivity and minimizing specificity (Fielding and Bell 1997; Phillips et al. 2006). We used these thresholds to convert probabilities to binary response (presence-absence) and used the equal test sensitivity and specificity threshold values to calculate successful classification percentages to corroborate results from ROC curves. We corroborated model selection using concordance (percent successful classification of sites) because a model that poorly classifies the data it was built from is unlikely to have any true predictive ability (Hosmer and Lemeshow 1989).
Once a preferred model was determined, we used the binary predictions to generate a model of 2 habitat quality classifications, high (probability of jaguar presence > 0.390) and low (probability of presence < 0.390). We discarded low quality habitat patches and used only high quality patches in further analyses. We calculated the area of each patch in order to differentiate between fundamental patches and stepping-stone patches. Fundamental patches encompassed at least 100 km2 of optimal jaguar habitat, which was the minimal area reported capable of holding a population of 5 individual jaguars (Núñez-Pérez 2011). In contrast, the minimal area for a stepping-stone was 1 km2 (Bunn et al. 2000; Bennett et al. 2004).
Identification of potential corridors. We used a functional connectivity model (Rabinowitz and Zeller 2010) utilizing a least cost-path method (LCP; LaRue and Nielsen 2008) to analyze jaguar movement opportunities. This method involves analyzing landscape structure and the response of the target species to the landscape (Dickson et al. 2005; Cushman et al. 2009). The result is a parameters value matrix that reflects the difficulty or ease for jaguar movements through the landscape, termed a permeability or resistance matrix. The permeability matrix in turn allows identification of the least-cost paths of movement (Rabinowitz and Zeller 2010).
We used the same landscape characteristics described for maximum entropy modeling to evaluate landscape structure and its permeability to jaguar movements (Figure 2). However, because habitat attributes that correlate with jaguar presence may not be the same as those that facilitate or inhibit long range movements of jaguars, we used Analytical Hierarchy Analysis (AHP; Saaty 1980) to develop the landscape permeability matrix rather than the maximum entropy modeling results. AHP is a process that utilizes pair-wise comparison matrices to clarify the relative importance of criteria that measure each elements (habitat variables) contribution to movements. Experts then compare every possible pairing and enter ratings, which are based on a continuous scale of nine points (Saaty 1980). We interviewed jaguar experts to quantify how landscape characteristics may influence jaguar movements, and analyzed their responses using Expert Choice 11v (Expert Choice Inc., Pittsburgh, PA; http://expertchoice.com). We selected the interviews with the best consistency index (IC) for each criterion to develop inputs for each variable in the permeability matrix. This evaluation requires a minimum value for data acceptance of IC < 1 (Malczewski 2004).
We then used the Multi-Criteria-Evaluation (MCE) module in IDRISI Andes (Eastman 2004) to assign the values derived from AHP using the Weighted Linear Combination (WCL) method (Eastman 2004). This allowed production of a raster reflecting the condition of the evaluation area with respect to the evaluated landscape elements (Malczewski 2004). The WCL is a weighted average, in which the landscape variables are input in raster format with weights assigned to each from the AHP. The decision maker directly assigns the weights of "relative importance" to each alternative by multiplying the importance weight assigned for each attribute by the scaled value given to the alternative on that attribute, and summing the products over all attributes. When the overall scores are calculated for all of the alternatives, the alternative with the highest overall score is chosen (Malczewski 2004). The output is a map of suitable habitat scaled from 0 to 1, where 1 = no resistance to movement. We reclassified the resultant map into 9 categories using the inverse of the score to create the permeability matrix; higher values thus indicated a major barrier to jaguar movements. We also calculated the percentage of the area occupied by each category.
Last, we used the permeability matrix to identify LCPs for jaguar movements using the "Cost Path" tool of ArcGIS 9.3 (ESRI, Redlands, California, USA). LCPs are represented as lines that result from cell selection for the shortest distance and the lowest costs of movements in the landscape permeability matrix. The LCPs can be interpreted as potential corridors (LaRue and Nielsen 2008; Cushman et al. 2009).
For LCPs, we selected as movement origins each of the jaguar records from the SMO, and used the potential habitat patches identified from maximum entropy modeling as destinations. We calculated the length of each route, and termed routes ≥ 20 km (maximum distance traveled regularly by a jaguar; Sollmann et al. 2011) dispersal routes. In contrast, we classified shorter routes (0.5 km to 19.9 km) as daily travel routes (Schaller and Crawshaw 1980). We interpreted the resultant biological routes regionally from north to south using each Mexican state as a separate region. This regionalization was necessary because all Mexican states are independent administrative units within the country, and therefore each state has its own administrative policy. Finally we described the relationship between route type and habitat patches relative to their location with respect to natural protected areas and Priority Terrestrial Regions within the study area.
Results
We obtained and classified 68 jaguar locations (Rosas-Rosas et al. 2002; Villordo-Galván et al. 2010; Hernández Saint-Martin et al 2013; Hernández Saint-Martin et al. 2014; Martinez-Calderas, per. comm.) that met the Class I credibility criterion of Tewes and Everett (1986; Figure 3). Maximum entropy modeling indicated that distribution of jaguars in the SMO was most strongly associated with slope and land use-vegetation type (Table 1). The best supported model contained slope and land use-vegetation type and fit data well (AUC = 0.881; SE = 0.056; Concordance = 79 %); slope accounted for 53% of the proportional variance in the model and likelihood of jaguar presence was > 0.70 in areas with slope < approximately 9o, then decreased rapidly until likelihood of presence was < 0.50 with slopes > 11o. Land use-vegetation type accounted for 47 % of the proportional variance in the model; presence of jaguar was most strongly associated with tropical forest (Likelihood > 0.70; Figure 4). Collectively, jaguar presence was most likely in tropical and oak forests with moderate slope.
a S = slope; L = landuse-vegetation type; R = distance to road; W = distance to water; E = elevation; P = human population density.
We identified 581 patches of suitable habitat ranging from 1 km2 to 740.5 km2 which totaled 23% of our evaluation area. Of these, 3 were > 100 km2 and thus met the criteria for fundamental patches. The remaining (578) patches were classified as stepping-stone (Figure 5). Of our evaluation landscape, 1 % was classed as fundamental patches and 99 % as stepping-stones (Figure 5).
We surveyed 9 jaguar researchers from universities and research centers in Latin America for our AHP analysis. Results of surveys suggested that areas with population densities < 7 inhabitants/ km2; with slopes from 0 - 3°; elevations of > 1000 m; presence of permanent water within 2 km; and > 16 km from roads were the most likely to be used by jaguar. Further, tropical forests, oak forests, and tropical deciduous forests were considered the most suitable vegetation types for jaguar in the study area (Table 2). Respondents considered the type of land use- vegetation and the human population density to be the most important landscape characteristics associated with jaguar presence (Table 3).
The resultant permeability matrix indicated that approximately 49 % and 3 % of the SMO constituted areas of medium and low resistance to jaguar movements, respectively (Figure 6). The remaining 48 % was characterized by attributes that probably represented a barrier for jaguar movements. Most of the areas of higher resistance and therefore the higher cost areas for jaguar movements were located in the western portion of the SMO (Figure 6).
Results of the LCP analysis showed 61 dispersal routes totaling > 2,000 km, and > 200 travel routes (Figure 7). The largest number of both route types was in Hidalgo (Table 4). Potentially significant movement corridors for jaguar in the SMO included:
In the northern region of our study area, routes originate in habitat patches in Queretaro Jalpan de Serra, Landa de Matamoros and Arroyo Seco in Querétaro and reach Tepehuacán de Guerrero and Jacala de Ledezma in San Luis Potosí and Hidalgo, respectively. A fourth route originates in San Ciro de Acosta and Lagunillas in San Luis Potosi and reaches Xichú municipality in Guanajuato.
In the middle region, numerous dispersal routes originate in Hidalgo (e. g. in localities including Pisaflores, Chapulhuacan, Tlanchinol, Yahualica, San Felipe Orizatlan, Hujutla de Reyes, and Atlapexco, Huatla) providing for movement throughout Hidalgo and into other states, including a route towards Veracruz that passes through through Chiconamel.
In the south, two potential dispersal routes interconnected by numerous travel routes could function as major jaguar corridors to northern areas of the state of Puebla. These corridors may be particularly important for jaguar dispersal south to other jaguar populations in southern México (Figure 8). The Sierra Madre Oriental corridors originate in the state of Veracruz through the Cuenca Hidrográfica del Río Necaxa including municipalities of Z. Mena, Venustiano Carranza, Xicotepec and the priority región of Cuetzalan through minicipalties such as Tuzamapa de Galeana, Jonotla and Zoquiapan.
Discussion
Our regional analysis of the southern SMO identified habitat patches sufficient to maintain multiple populations of jaguars (fundamental patches), as well as numerous smaller habitat patches that provided the resources necessary to facilitate movements of jaguars (stepping-stones). These results affirm the importance of the SMO for the conservation of jaguars in México in both the short and long term (Grigione et al. 2009; Rodríguez-Soto et al. 2011). Although the fundamental patches we identified were separated by relatively long distances (Figure 7) with respect to the movement patterns of the jaguar, this was mitigated by the number and distribution of stepping-stones as intermediates that provided connectivity among the fundamental patches. Because of the number, function, and importance of the stepping-stone patches, the fundamental patches were likely not isolated despite the relatively high degree of fragmentation. Stepping-stones help maintain the populations that reside in the fundamental patches by facilitating the flow of individuals and serving as forage and cover stations during dispersal (Bennett et al. 2004).
The biological corridor concept has evolved into a more comprehensive paradigm from its origin of simple linear features (Bennett and Mulongoy 2006). Currently, corridors are viewed as a mosaic of land use types and other landscape features that allow connectivity among patches of suitable habitat across the landscape (Bennett and Mulongoy 2006). Understanding the role of the landscape matrix as habitat for wildlife, particularly in fragmented landscapes such as the SMO, is important because the distribution of habitat components influences the movement patterns of the target species (Riley et al. 2003). Our results corroborated the presence of suitable habitat for the jaguar where jaguars have been located and studied in the states of Tamaulipas, San Luis Potosi, and Querétaro (Ortega-Huerta and Medley 1999; Bravo and González 2006; Villordo-Galván et al. 2010). Moreover, our model predictions indicated several additional areas with high likelihood of suitable habitat in little studied areas, including Veracruz, Hidalgo, and Puebla. A priority should be to survey these areas for jaguar presence, and to analyze with greater detail the potential connectivity between these patches and the known populations of jaguar that occur to the south.
Our permeability matrix indicated a high proportion (48 %) of areas with high resistance to jaguar movement, likely due to the degradation of the landscape, particularly changes in land use to agriculture and cattle ranching (Williams-Linera et al. 2003; Rivera et al. 2010). Ongoing development results in only remnants of original habitat being available for jaguars, and reaffirms the need to protect and restore habitat that connects core jaguar populations (Villordo-Galván et al. 2010; Ávila Nájera et al. 2011). However, despite high disturbance and fragmentation in this region of México, our results and previous studies indicate that biodiversity is still high in the Huasteca Region of San Luis Potosí (Villordo-Galván et al. 2010). In addition, jaguars are resident in this region of San Luis Potosí and likely juveniles or adults are dispersing towards southern areas through the proposed corridors (Villordo-Galván et al. 2010; Ávila Nájera et al. 2011).
Moreover, it is possible that areas classified as having high resistance to jaguar movements could result simply from the characteristics of the habitat variables proposed by experts from other countries exceeding the standards of suitability for jaguars in México. This would result in an apparent lower level of habitat suitability in the SMO than is actually present. This could also have occurred when the Jaguar Conservation Units were delineated by Sanderson et al. (2002), which declared the SMO as likely isolated from other jaguar populations in México and therefore as an area of secondary importance for long-term jaguar conservation. This conclusion was made despite a lack (at that time) of recent field surveys, and was likely premature. A similar case was presented by Morrison et al. (2012). Therefore, it is important to continue monitoring studies to better understand the current condition of the jaguar habitat and movement patterns throughout their range in México.
In our evaluation we analyzed the suitability of the SMO to provide corridors for movements of jaguars within the SMO and among other regions. The routes described were based on the best available data and knowledge, although they require further validation in field. With respect to existing landscape policy in the SMO, Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO; http://www.conabio.gob.mx) identified 4 Priority Terrestrial Regions (RTP) in our evaluation area. Three of these - the Sierra del Abra Tanchipa, Sierra Gorda-Rio Moctezuma, and Bosque Mesofilos de la Sierra Madre Oriental - are postulated to act as biological corridors, whereas the fourth region (RTP-Cuetzalan) had no information to support this biological function (CONABIO; Fig 1). Our results confirm the biological functions of three of the four areas defined by CONABIO, while providing further information on the function of the RTP-Cuetzalan. We identified several travel or dispersal routes which converged or originated in this RTP, providing potential routes of dispersal for jaguars into the Sierra Norte de Puebla. For this reason, the RTP-Cuetzalan is important for connectivity not only for the jaguar but also for other species that inhabit this region. Consequently, we recommend that the RTP-Cuetzalan be considered for designation as a Natural Protected Area, because it still retains significant high quality habitat for wildlife which facilitates its ability to support individuals dispersing to other areas, as well as providing for many resident species.
In addition to the existing data on each RTP, vertebrate studies in the SMO further indicate the presence of the main prey of the jaguar (Núñez et al. 2000; Rosas-Rosas et al. 2008) in San Luis Potosí (Ávila-Nájera et al. 2011; Rueda et al. 2013), Hidalgo (Sánchez-Rojas et al. 2009), and Puebla (Villarreal-Espino-Barros et al. 2008; Bravo et al. 2010). Because the distribution and movements of jaguars depend largely on the presence of food, these records further indicate the potential of the SMO to support the long-term needs and movements of the jaguar (González and Miller 2002).
Last, the accuracy of model predictions such as our maximum entropy results and permeability matrix are determined by the quality of data. We used an appropriate number of presence records and landscape variables that functionally relate to jaguar movements, which resulted in a good model based on model fit criteria (Guisan et al. 2006; Anderson and Gonzalez 2011). However, it is also necessary to consider the reliability of the records (Mateo et al. 2011). Despite limiting our modeling to only highly reliable jaguar records within the region, and thus limiting our overall sample size, the predictive ability of our preferred maximum entropy model was good, possibly due to the precise demarcation of the study area (Soberón 2010) and the reliability of the records that we used. In addition, we consider at least one record in five of the seven States studied (i. e. 70 % of the area contained information for input into models).
Further, the landscape attributes we found best associated with suitable habitat for the jaguar were also frequently identified in the literature (Quigley and Crawshaw 1992; Monroy Vilchis et al. 2008; Conde et al. 2010; Villordo-Galván et al. 2010; Colchero et al. 2011). Habitat quality and distribution of jaguars is most frequently associated with specific vegetation types and land use (Monroy-Vilchis et al. 2008; Rosas-Rosas and Valdez 2010; Colchero et al. 2011). The variables we modeled thus allowed us to locate specific areas with habitat features highly conducive to the presence of jaguars.
Implications. Biological corridors of the SMO initiate at fundamental patches in the States of Tamaulipas and San Luis Potosí, with continuity towards the State of Queretaro. In Queretaro, historic and present reports of jaguar presence and the presence of primary prey suggests that suitable habitat conditions exist to allow jaguars to meet their basic requirements (Téllez and López 1995). Although there are no historical or recent records of jaguars in Hidalgo, we identified routes that connect directly with the state of Puebla, with further possibilities of connection to northeastern Oaxaca where the presence of the jaguar was recently documented (Figures 7 and 8; Briones-Salas et al. 2012).
In Puebla, routes apparently diverge and may indicate the presence of a barrier (Figure 6). However, the more westerly route is promising because of proximity among multiple small fragments of habitat, presence of multiple travel routes, and proximity to the Necaxa River hydrographic basin. The Necaxa River watershed is important because of its similarity to the Sierra del Abra Tanchipa in San Luis Potosi. Like the Sierra del Abra Tanchipa, which is also surrounded by fragmented landscapes of secondary vegetation and agricultural and livestock areas, this region has a core of habitat conducive to prey and vegetation coverage that can provide suitable habitat for jaguars and feasibly allow dispersal (as well as support permanent populations) in these highly fragmented landscapes.
We conclude that the SMO provides multiple biological corridors for connectivity with more southern populations of jaguar in México (Figure 8). Based on this, we recommend development of management strategies to facilitate exchange of individuals from the SMO with other known jaguar populations in southern states (e. g. Oaxaca, Tabasco, Guerrero). Principal strategies may focus on habitat management of fundamental patches and habitat restoration between or around stepping-stones. The quality of the habitat patch is important because it determines availability of prey. Finally we suggest continued field research in order to seasonally monitor prey populations, jaguar presence, and movements. Working with local communities to improve habitat, enhance conservation of landscapes, and conduct educational programs on minimizing human-jaguar conflicts is also needed.
This information can also help to identify areas where management policies should be revised and also where future development should be reviewed and possibly minimized to mitigate any potential negative effects to critical patches or corridors. Because corridors overlap multiple jurisdictions, it is also important to coordinate local and regional conservation efforts so that individual actions can be linked across larger regions to form cohesive connectivity networks.