Introducción

El margay (Leopardus wiedii, Schinz 1821) es el más pequeño de los felinos silvestres en México y está en peligro de extinción debido a la destrucción y fragmentación de su hábitat por la conversión a tierras agrícolas y ganaderas y a desarrollos urbanos; también está amenazado por el tráfico y comercio ilegal de pieles y ejemplares vivos y el sacrificio incidental de individuos que cazan aves de corral (^{Oliveira et al., 2015}). A nivel mundial, el margay aparece como especie casi amenazada en la Lista Roja de la Unión Internacional para la Conservación de la Naturaleza (^{Oliveira et al., 2015}) y como especie sobre la que no se pueden comercializar individuos o sus partes de manera legal según el Apéndice I de la Convención sobre el Comercio Internacional de Especies Amenazadas de Fauna y Flora Silvestres (^{CITES, 2015}). En nuestro país, el tamaño de las poblaciones de margay se desconoce y es posible que se esté reduciendo continuamente por lo que se incluyó dentro de la Norma Oficial Mexicana NOM-059-SEMARNAT-2010 en la categoría de especie en peligro de extinción (^{Semarnat, 2010}).

La conservación del margay es importante ya que este pequeño felino posee características morfológicas particulares, en comparación con otros mesocarnívoros, que le permiten jugar una variedad de funciones que serían de difícil sustitución en los ecosistemas donde habita (^{Farias y Svensson, 2014}). Es un cazador solitario que captura a sus presas tanto en el suelo como en los árboles porque puede desplazarse en el dosel y quedar suspendido de las ramas por las extremidades traseras al flexionar la articulación anterior al tarso con una supinación de hasta 180° (^{Oliveira, 1998}; ^{Oliveira et al., 2015}; ^{Tewes y Schmidly, 1987}). Sus hábitos arbóreos le permiten prevenir encuentros agonísticos desfavorables con especies competidoras (^{Oliveira et al., 2015}). Entre sus presas registradas a lo largo de su distribución, se incluyen anfibios, lagartijas, roedores, ardillas y pequeños primates (^{Bianchi et al., 2011}; ^{Cinta-Magallón et al., 2012}; ^{Oliveira, 1998}; ^{Oliveira et al., 2015}; ^{Oliveira-Calleia et al., 2009}).

Aun cuando el margay muestra un área de distribución histórica desde Sonora y Tamaulipas en México hasta el norte de Argentina y noroeste de Uruguay (^{Hall, 1981}; ^{Oliveira, 1998}), el área ocupada actualmente por esta especie debe ser mucho menor debido a su preferencia por la vegetación nativa con cobertura densa y en buen estado de conservación (^{Aranda et al., 2012}; ^{Aranda-Sánchez y Valenzuela-Galván, 2015}; ^{Briones-Salas et al., 2016}; ^{Domínguez-Castellanos y Ceballos, 2005}; ^{Oliveira, 1998}; ^{Oliveira et al., 2015}). Se ha registrado la presencia del felino desde el nivel del mar hasta alrededor de los 1,000 m snm, aunque hay algunos registros a 3,000 m (^{Oliveira, 1998}; ^{Oliveira et al., 2015}). Las estimaciones de densidad poblacional del margay son escasas y se han utilizado diferentes métodos para obtenerlas; en general se reportan de 1 a 5 individuos/100 km², pero en algunos sitios de la región amazónica y en ausencia del ocelote (Leopardus pardalis) se han hecho estimaciones de hasta 15 a 25 individuos/100 km² (^{Di Bitteti et al., 2010}; ^{Oliveira et al., 2010}, ²⁰¹⁵). Además de la baja densidad poblacional y de su reducida área de ocupación, se trata de una especie difícil de registrar por sus hábitos nocturno-crepusculares y arborícolas, a pesar de su amplia distribución geográfica (^{Oliveira et al., 2015}). En México, la distribución histórica del margay se registró a lo largo de ambas costas, por el océano Pacífico, desde el sur de Sonora y por el golfo de México, desde Tamaulipas hasta unirse en Oaxaca y continuar al istmo de Tehuantepec, Chiapas y la península de Yucatán (^{Hall, 1981}; ^{Leopold, 1959}; ^{Oliveira et al., 2015}; ^{Tewes y Schmidly, 1987}). Los reportes recientes de fotocapturas de margay documentaron su presencia en sitios anteriormente no incluidos dentro del área de distribución histórica en los estados de Guanajuato (^{Iglesias et al., 2008}); Morelos (^{Aranda y Valenzuela-Galván, 2015}; ^{Valenzuela-Galván et al., 2013}); Puebla (^{Farías et al., 2015}) y San Luis Potosí (^{Martínez-Calderas et al., 2012}), y con ello aumentaron el área de distribución conocida para el felino.

Un método eficiente y robusto para elaborar hipótesis acerca de la distribución geográfica potencial de una especie es el modelado del nicho ecológico, a partir de registros georreferenciados y de la información de variables ambientales como predictores (^{Sánchez-Cordero et al., 2001}). Los modelos de distribución potencial son un recurso importante en la biología de la conservación debido a que son utilizados como representaciones del espacio geográfico idóneo para la presencia de la especie en estudio (^{Mateo et al., 2011}; ^{Sánchez-Cordero et al., 2001}). Sin embargo, también hay que tomar en consideración que los sitios con condiciones climáticas favorables no necesariamente estarán ocupados por la especie bajo estudio debido a que la distribución geográfica de los organismos también está determinada por las interacciones entre especies y por su historia biogeográfica, y estos factores no pueden ser incorporados en los algoritmos (^{Guisan y Zimmermann, 2000}; ^{Soberón y Peterson, 2005}).

Aún con esta desventaja, el modelado del nicho ecológico se ha convertido en información fundamental para enfocar los esfuerzos de conservación (^{Sánchez-Cordero et al., 2001}; ^{Soberón y Peterson, 2005}). El algoritmo de MaxEnt (Maximum Entropy, ^{Phillips et al. 2006}, ²⁰⁰⁹) posee robustez y alta capacidad predictiva basado en la entropía máxima; por lo tanto, es de los más utilizados para estimar la idoneidad de las condiciones ambientales que a su vez permiten la presencia de una especie (^{Elith et al., 2011}; ^{Merow et al., 2013}).

Las áreas naturales protegidas (ANP) pueden ofrecer espacios geográficos con mayor probabilidad de cumplir con el objetivo de conservación de los felinos silvestres (^{Bezaury-Creel y Gutiérrez-Carbonell, 2009}; ^{Estes et al., 2011}; ^{Woodroffe y Ginsberg, 1998}). Los felinos silvestres se benefician por la protección, y la abundancia de sus poblaciones está relacionada positivamente con el tamaño de las reservas (^{Di Bitteti et al., 2010}; ^{Massara et al., 2015}, ²⁰¹⁶). Debido a que las ANP son áreas estáticas y que el tamaño es uno de los principales factores en cuanto a su eficiencia, la Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (Conabio) hizo una propuesta para la expansión de varias ANP, así como para la creación de una nueva figura de conservación, a lo que llamó las Regiones Terrestres Prioritarias y que es un instrumento de planeación territorial (^{Arriaga et al., 2000}).

Las áreas que mantienen altos niveles de integridad ecológica, como muchas de las ANP, son importantes para la conservación del margay debido a que este felino no se adapta fácilmente a hábitat degradado debido a su preferencia por la vegetación nativa densa y a que busca refugio entre las copas de los árboles y en cañones rocosos con pendiente pronunciada (^{Oliveira, 1998}; ^{Oliveira et al., 2015}). Nuestro objetivo fue generar un modelo de distribución potencial de L. wiedii en México para identificar ANP desprovistas de registros georreferenciados de margay y con condiciones ambientales favorables para su presencia. El modelo de distribución potencial fue nuestra hipótesis sobre la distribución geográfica actual del margay, y nuestra predicción fue que podríamos identificar ANP donde sería factible verificar la presencia del margay.

Materiales y métodos

Elaboramos una base de datos biológica con la información de los registros de presencia de L. wiedii para México con coordenadas geográficas de 1944 a 2015. Se obtuvo la información disponible de los registros de presencia de L. wiedii en México en los portales de redes de bases de datos interconectadas “Mammal Networked Information System” (MaNIS), “Global Biodiversity Information Facility” (GBIF), Irekani del Instituto de Biología de la Universidad Nacional Autónoma de México (IBUNAM) y Unidad de Informática para la Biodiversidad (UNIBIO) del IBUNAM. También se realizó una búsqueda de registros georreferenciados en artículos científicos (^{Almazán-Catalán et al., 2013}; ^{Aranda-Sánchez y Valenzuela-Galván, 2015}; ^{Botello et al., 2006}, ²⁰⁰⁸; ^{Faller-Menéndez et al., 2005}; ^{Farías et al., 2015}; ^{Horváth et al., 2008}; ^{Iglesias et al., 2008}; ^{Lira-Torres et al., 2005}, ²⁰¹²; ^{Martínez-Calderas et al., 2012}, ²⁰¹⁶; ^{Meave-del Castillo y Luis-Martínez, 2000}; ^{Morales-Pérez, 2005}; ^{Navarrete-Gutiérrez et al., 1996}; ^{Pérez-Irineo y Santos-Moreno, 2012}; ^{Valenzuela-Galván et al., 2013}) y se revisaron los ejemplares depositados en la Colección Nacional de Mamíferos (CNMA) del Departamento de Zoología del IBUNAM. No utilizamos registros de observaciones directas ni de registros fotográficos que no estuvieran incluidos en alguna colección científica o citados en algún artículo científico. Se incluyeron 12 campos, 3 con información taxonómica (género, especie y subespecie), 5 con información geográfica (latitud, longitud, estado, municipio, localidad) y 4 de la colección biológica (nombre de la colección, número de catálogo, tipo de ejemplar, fecha de recolecta).

El modelo del nicho ecológico de L. wiedii se generó con el software MaxEnt versión 3.4.0 (Maximum Entropy, ^{Phillips et al., 2017}) a partir de los registros de presencia y 19 coberturas digitales de variables climáticas elaboradas por ^{Cuervo-Robayo et al. (2014)} con datos climáticos de 1910 a 2009 y resolución espacial de 1 km². Para eliminar la correlación geográfica de los datos utilizamos la herramienta SDMtoolbox v2.4 para ArcGis 10.0 (^{Brown, 2020}). El software eliminó los registros múltiples de presencia dentro de un área determinada y seleccionó registros independientes espacial y ambientalmente a partir de una matriz de registros coordenados y sus perfiles bioclimáticos. La herramienta “spatially rarefy occurrence data” filtró los datos de las localidades y redujo los registros de presencia a un solo punto dentro de una distancia Euclideana especificada. Esta herramienta se usó para aplicar rarefacción a los registros de margay a varias distancias de acuerdo con la heterogeneidad topográfica y climática y fueron filtrados espacialmente a 10 km², 20 km² y 25 km² en áreas con heterogeneidad alta, media y baja, respectivamente. Este método de filtrado gradual es muy útil para estudios con pocos registros ya que puede maximizar el número de localidades independientes.

Se generó un modelo utilizando los parámetros predeterminados por MaxEnt y con las opciones de ignorar datos duplicados, considerar siembra de datos al azar y no considerar acotación ni extrapolación. El 75% de los registros se utilizó para el entrenamiento del modelo y 25% para la evaluación del modelo. El número máximo de puntos de fondo fue de 1,000. La región M (^{Soberón y Peterson, 2005}) fue recortada de la cobertura digital de ecorregiones de la Wildlife World Fund (WWF, ^{Olson et al., 2001}) y el recorte se basó en las ecorregiones con registros de margay por lo que la península de Baja California no se incluyó. El análisis jackknife de MaxEnt estimó la contribución relativa (%) de cada una de las 19 variables al modelo generado y a partir de esta información, se seleccionó el subconjunto de variables que presentaron los mayores porcentajes de contribución y que sumados explicaron de 90 a 95% de la variabilidad. Tras la selección de las variables: (i) precipitación del cuatrimestre más lluvioso, (ii) oscilación anual de la temperatura, (iii) temperatura mínima promedio del periodo más frío, (iv) estacionalidad de la temperatura, (v) estacionalidad de la precipitación, (vi), temperatura promedio del cuatrimestre más frío (vii) precipitación del cuatrimestre más seco y (viii) precipitación anual, repetimos el ejercicio de modelado con estas 8 variables para generar 15 réplicas y su promedio utilizando la opción de submuestra que replicó conjuntos de muestras elegidos y eliminó el porcentaje de prueba aleatoria sin reemplazo. Se obtuvo la curva ROC/AUC (acrónimo en inglés de Receiver Operating Characteristic / Area Under the Curve). La curva ROC evaluó la habilidad predictiva del modelo y AUC midió la capacidad del modelo para clasificar correctamente los registros de la especie como presencia o ausencia (^{Fawcett, 2006}; ^{Phillips et al., 2006}). Utilizamos la ROC, y no la ROC parcial, debido a que no contamos con muchos registros considerando la extensión del país y la distribución de la especie (^{Peterson, 2011}). Utilizamos la clasificación de ^{Araújo y Guisan (2006)} del valor de AUC que indica la precisión del modelo como insuficiente (0.50 - 0.60), pobre (0.60 - 0.70), promedio (0.70 - 0.80), bueno (0.80 - 0.90) o excelente (0.90 - 1.00).

El mapa de distribución potencial logístico se desplegó en el sistema de información geográfica ArcView 3.3 (Environmental Systems Research Institute [ESRI], Redlands, California, EUA). El criterio de punto de corte que utilizamos fue el del umbral logístico de presencia de entrenamiento de percentil 10 (^{Phillips et al., 2006}) y en nuestro modelo tuvo el valor de ≥ 0.4125. El mapa de distribución potencial del margay fue proyectado mediante ArcView sobre coberturas digitales de mapas de a) ANP federales de México (^{Conanp, 2016}), b) ANP estatales, municipales, ejidales y privadas de México (^{Conabio, 2015}) y c) ecorregiones terrestres de México (^{INEGI-Conabio-INE, 2008}). Determinamos las ANP que mostraron las condiciones ambientales favorables para la presencia de L. wiedii como aquellas dentro de los límites del modelo logístico de distribución potencial. Identificamos a las ANP con valores de idoneidad más altos (0.8 - 1.0) de acuerdo con el modelo logístico. Utilizamos el mapa de ecorregiones terrestres en México para identificar los tipos de vegetación presentes en las ANP y en los sitios donde los registros fueron obtenidos (^{Inegi-Conabio-INE, 2008}).

Resultados

De un total de 97 registros de L. wiedii para México con coordenadas únicas, 66 fueron registros independientes; es decir, sitios sin correlación espacial ni ambiental y se utilizaron para generar el modelo de nicho ecológico. El total de los 97 registros de presencia de margay se distribuyeron en 18 estados de la República Mexicana. El 74% se localizaron en San Luis Potosí (n = 23), Oaxaca (15), Colima (10), Chiapas (9), Jalisco (9) y Guerrero (6). Los estados con menos de 5 registros fueron: Campeche (4), Morelos (4), Veracruz (4), Puebla (3), Durango (2), Guanajuato (2), Chihuahua (1), Nayarit (1), Quintana Roo (1), Sinaloa (1), Sonora (1) y Yucatán (1). Asimismo, del total de los 66 registros independientes, la mayoría se concentraron en: San Luis Potosí (15), Oaxaca (11), Chiapas (8), Jalisco (6) y Guerrero (5), y los estados con menos de 5 registros fueron: Campeche (3), Colima (3), Morelos (3), Veracruz (3), Durango (2), Puebla (2), Chihuahua (1), Guanajuato (1), Quintana Roo (1), Sinaloa (1) y Yucatán (1). De acuerdo con el mapa de ecorregiones terrestres en México (^{Inegi-Conabio-INE, 2008}), los 97 registros georreferenciados de margay se ubicaron casi en su totalidad en 3 ecorregiones que fueron las sierras templadas (n = 50), las selvas cálido-húmedas (25) y las selvas cálido-secas (21), y 1 sitio se ubicó en la ecorregión de elevaciones semiáridas meridionales. Los 66 registros independientes también se ubicaron en las sierras templadas (27), las selvas cálido-húmedas (20) y las selvas cálido-secas (18), y 1 registro se mantuvo en las elevaciones semiáridas meridionales. Los tipos de vegetación donde se localizaron los 97 registros fueron bosque de pino-encino (n = 20), selva baja caducifolia (18), bosque de encino (16), selva alta y mediana perennifolia (9), selva mediana (9), selva alta y mediana subperennifolia (6), bosque de pino (5), selva baja (5), bosque mesófilo de montaña (4), pastizal cultivado (3) y matorral xerófilo (2).

Elegimos como modelo de distribución potencial del margay al modelo promedio de las 15 réplicas (Fig. 1). La eficiencia de predicción del modelo lo ubicó en el intervalo del valor de AUC de 0.8574 con una precisión clasificada como buena de acuerdo con los valores considerados por ^{Araújo y Guisan (2006)}. El modelo predijo mayor probabilidad de presencia de condiciones favorables para el margay en las sierras Madre Oriental, Norte de Oaxaca y Madre Occidental. La predicción tuvo menor probabilidad en la depresión del Balsas, la planicie costera de Michoacán, Guerrero, Oaxaca y Chiapas y en las regiones de manglares y humedales de Tabasco y las regiones norte y este de la península de Yucatán. El subconjunto de 8 variables ambientales que presentaron los mayores porcentajes de contribución y que sumados explicaron de 90 a 95% de la variabilidad estuvieron relacionados con la estacionalidad de la temperatura y de la precipitación, los periodos con mayor precipitación y los periodos más fríos.

Figura 1 Modelo de distribución potencial del margay (Leopardus wiedii) en México con los polígonos de las áreas naturales protegidas federales. Las estrellas identifican a los 66 registros de margay utilizados para generar el modelo. Los colores cálidos (rojo y anaranjado) indican los valores de idoneidad ambiental más altos (0.8 - 1.0) para la presencia del margay. Los tonos de colores representan valores de un continuo entre 0 y 1 de idoneidad desde el color más frío (azul) hacia el color más cálido (rojo). 

De los 97 sitios con registros georreferenciados de margay, 39 se localizaron dentro de un ANP federal y 1 en ANP privada (Tabla 1). No hubo registros georreferenciados en ANP estatales, municipales y ejidales. Mediante el modelado del nicho ecológico del margay identificamos 40 ANP federales con condiciones climáticas favorables para la presencia del felino de las cuales no encontramos registros con coordenadas geográficas en las colecciones científicas ni en los artículos científicos (Tabla 2). El modelado predijo que 17 de las 40 ANP contaron con condiciones ambientales de alta idoneidad (0.8 - 1.0) de acuerdo con el modelo logístico (Fig. 2) para la presencia de este felino (Tabla 2). Los 40 registros de margay que se localizaron dentro de algún ANP se ubicaron en bosque de pino-encino (n = 12), bosque de encino (8), selva baja caducifolia (4), selva alta y mediana subperennifolia (4), pastizal cultivado (3), bosque mesófilo de montaña (2), matorral xerófilo (2) y bosque de pino (1).

Figura 2 Modelo de distribución potencial del margay (Leopardus wiedii) en México con los polígonos de las áreas naturales protegidas federales. El color rojo indica los valores de idoneidad ambiental más altos (0.8 - 1.0) para la presencia del margay. 

Discusión

En nuestra búsqueda de registros de presecia de margay, encontramos el mayor número en San Luis Potosí y Oaxaca porque fue en estos estados de la República Mexicana donde se desarrollaron la mayoría de los trabajos de investigación que publicaron registros georreferenciados de margay en años recientes (^{Botello et al., 2006}, ²⁰⁰⁸; ^{Cinta-Magallón et al., 2012}; ^{Lira-Torres et al., 2005}, ²⁰¹²; ^{Martínez-Calderas et al., 2012}, ²⁰¹⁶; ^{Pérez-Irineo y Santos-Moreno, 2012}). En contraste, las investigaciones realizadas en el suroeste del Estado de México, el sur de Michoacán y el sureste de Tamaulipas verificaron la presencia del margay por medio de fotocapturas y capturas de ejemplares, pero no publicaron las coordenadas geográficas de los sitios de colecta (^{Carvajal-Villareal et al., 2012}; ^{Charre-Medellín et al., 2015}; ^{Monroy-Vilchis et al., 2011}). En la península de Yucatán los registros fueron muy escasos y no obtuvimos ningún registro en Tabasco (^{Faller-Menéndez et al., 2005}). De esta manera, se comprende que el mayor número de registros en San Luis Potosí y Oaxaca, así como el menor número de registros en el sureste mexicano no estuvieron relacionados con una mayor o menor abundancia de las poblaciones del margay.

Es notable mencionar que la mayor parte de los registros colectados en años recientes fueron fotocolectas. En comparación con métodos de colecta anteriormente utilizados como la trampa de pelo, las estaciones olfativas y la captura de ejemplares vivos, el fototrampeo posee ventajas para la captura de registros del margay ya que es una especie con actividad nocturna principalmente, sigilosa y con baja densidad poblacional (^{Downey, 2005}; ^{Goldman, 1943}; ^{Iglesias et al., 2008}; ^{Vanderhoff et al., 2011}). El fototrampeo es un método no invasivo que permite la identificación del margay con relativa facilidad y la probabilidad de detección es confiable con un esfuerzo de muestreo apropiado aun cuando la abundancia relativa del felino sea baja (^{Tobler et al., 2008}).

El modelo generado en nuestro estudio permitió obtener una hipótesis bien sustentada sobre distribución geográfica del felino en la que se puede apreciar la información reciente sobre el área de distribución en el centro de México (Fig. 1). A pesar de que en la península de Yucatán el modelo resultó en un subestimado en comparación con el resto del país debido al número tan reducido de registros y a que no se capturó la variación ambiental, no se presentaron omisiones.

La probabilidad de encontrar las condiciones ambientales que permiten la presencia de una especie se interpreta como un índice de similitud ambiental entre una región determinada y los sitios donde la especie ha sido registrada (^{Guisan y Zimmermann, 2000}; ^{Peterson y Soberón, 2012}). El clima en los puntos de registro de la especie es representante de la gama ambiental del organismo; sin embargo, en caso de que no se cuente con una muestra representativa de la gama ambiental completa pueden presentarse omisiones en el modelo (^{Guisan y Thuiller, 2005}; ^{Guisan y Zimmermann, 2000}; ^{Soberón y Peterson, 2005}; ^{Phillips et al., 2006}). Nuestro modelo no presentó omisiones y obtuvo un valor aceptable del AUC, por lo que asumimos que es un modelo robusto para dirigir el diseño de muestreos hacia sitios específicos en las ANP para verificar la presencia del margay. Asimismo, aun cuando se contó con un número reducido de registros, el modelo de distribución potencial se mostró robusto, ya que coincidió con la distribución actual basada en los registros y la potencial, de una manera que parece consistente considerando la distribución histórica (^{Hall, 1981}).

Los hábitos sigilosos del margay no son compatibles con los métodos de captura de ejemplares, lo que pudo ser un impedimento para que la especie se pudiera registrar de manera regular en un pasado, ya que debido a su amplia distribución podría haber habido más registros en las colecciones científicas. Sin embargo, el número y distribución de los registros conocidos reflejó los ambientes en donde habita la especie. Finalmente, este estudio permitió identificar áreas y ANP en donde se podría específicamente aplicar el muestreo con cámaras trampa con el fin de obtener registros georreferenciados y disponer de más y mejores datos para hacer evaluaciones y recomendaciones para la conservación del margay.

Nuestra predicción se cumplió; identificamos 17 ANP federales con condiciones climáticas con alta idoneidad para la presencia del margay, pero en las cuales no encontramos registros con coordenadas geográficas de ejemplares resguardados en colecciones científicas o reportados en artículos científicos (Tabla 2). Es importante mencionar que en los planes de manejo de las 17 ANP sin registros y otras más, se incluyó la presencia de la especie. Es posible que se haya incluido al felino en los listados faunísticos de las ANP que están localizadas dentro de la distribución histórica del margay (^{Hall, 1981}) aún sin tener ejemplares en colecciones o fotocolectas científicas debido a que el margay es una especie críptica muy difícil de observar en vida silvestre y de capturar (^{Domínguez-Castellanos y Ceballos, 2005}). Asimismo, existen registros y menciones en la literatura científica de la presencia del margay en áreas adyacentes a ciertas ANP. Por ejemplo, San Luis Potosí fue el estado del que encontramos el mayor número de registros georreferenciados, pero ninguno estuvo localizado dentro de algún ANP.

En nuestros resultados, la Reserva de la Biosfera (RB) Sierra del Abra Tanchipa, en San Luis Potosí, presentó valores de idoneidad altos de acuerdo al modelo logístico y concordó con el trabajo de ^{Martínez-Calderas et al. (2016)}, quienes generaron un modelo de distribución potencial del margay para el noreste de México y reportaron que el área con mayor potencial para la presencia del margay se encontró en las regiones montañosas de la sierra Madre Oriental entre la Gran Sierra Plegada y el Karso Huasteco; región que incluye a la RB Sierra del Abra Tanchipa. Otro caso es El Parque Nacional (PN) Cañón del Sumidero, Chiapas; fue otro sitio que nuestro modelo logístico predijo con alta idoneidad climática y coincidió con ^{Arroyo-Chacón et al. (2013)}, quienes mencionaron avistamientos de margay por parte de los pobladores, a pesar de que no documentaron fotocapturas o capturas de ejemplares. Por otro lado, también encontramos comentarios sobre la presencia del margay en áreas adyacentes a ANP predichas por nuestro modelo logístico que no contaron con alta idoneidad de condiciones climáticas. ^{Quijano-Hernández y Calmé (2002)} reportaron que cazadores del ejido de Felipe Carrillo Puerto que colinda con la RB Sian Ka’an, Quintana Roo, confieren valor comercial a la piel del margay. La Reserva Privada El Zapotal, en Yucatán, cuenta con un registro de margay y colinda con la RB Ría Lagartos (^{Faller-Menéndez et al. 2005}). ^{Buenrostro-Silva et al. (2012)} presentaron un listado de los mamíferos del PN Lagunas de Chacahua, Oaxaca y aunque no obtuvieron registros del margay, señalaron que previo a su estudio, ^{Lira-Torres et al. (2005)} registraron al felino en el cerro de La Tuza de Monroy, adyacente al PN Lagunas de Chacahua.

A partir de la información disponible en literatura y colecciones científicas, y de nuestro modelo de distribución potencial, resultó evidente la necesidad de obtener un mayor número de registros confiables de presencia del felino en México. Por lo tanto, sugerimos realizar investigaciones en campo para comprobar la presencia del margay en las ANP, identificar la ubicación precisa de sus poblaciones, así como los tipos de vegetación en que habitan y proponemos resaltar los sitios de alta idoneidad reportados en este trabajo (Tabla 2). Esta información será útil para tomar medidas para la protección del margay y de su hábitat a corto, mediano y largo plazo.