Evaluación de Criterios del Anexo Normativo II, Método de evaluación del riesgo de extinción de plantas en México

Criterio A. Características de la distribución geográfica. Subcriterio A-1), extensión de la distribución. Se consultaron los registros de presencia de la especie en los herbarios. Los registros de presencia de la especie en la península de Yucatán fueron superpuestos al mapa de uso de suelo y vegetación Serie VI del ^{INEGI (2016)}, se descartaron aquellos que se encuentran en zonas urbanas o de uso agropecuario. La superficie ocupada por B. pliabilis fue estimada con base en la extensión de los remanentes de vegetación primaria que son hábitat de esta especie. También se consideró para la estimación el trabajo en campo. Subcriterio A-2), número de poblaciones o localidades conocidas existentes. Las localidades de los sitios de presencia se superpusieron a un mapa 1:4,000,000 para discernir las poblaciones. También, los registros se superpusieron a las cartas de uso de suelo y vegetación del INEGI Serie II (^{INEGI, 2001}) y Serie VI (^{INEGI, 2016}) para identificar las poblaciones remanentes. Los períodos fueron seleccionados considerando la fecha de inclusión de la especie en la ^{NOM-059-ECOL-1994} como especie amenazada y una fecha reciente para su comparación. Para estimar el aislamiento geográfico entre las poblaciones, se calculó espacialmente la distancia euclidiana entre las 6 poblaciones silvestres muestreadas. Subcriterio A-3), número de provincias biogeográficas en las que se encuentra el taxón (distribución histórica). Los registros conocidos se superpusieron al mapa de provincias biogeográficas (^{Conabio, 1997}). La representatividad de la distribución de B. pliabilis en el territorio mexicano se estimó considerando el área de distribución biogeográfica conocida para esta especie.

Criterio B. Características del hábitat. Subcriterio B-1), ¿en cuántos tipos de vegetación se presenta? Los registros conocidos de B. pliabilis se superpusieron al mapa de vegetación potencial de ^{Rzedowski (1990)} para identificar los tipos de vegetación donde se presenta. Los subcriterios B-2), ¿el taxón tiene un hábitat especializado? B-3), ¿la permanencia de la población es dependiente de un hábitat primario? y B-4), ¿la permanencia de la población requiere de regímenes de perturbación particulares o está asociada a etapas transitorias en la sucesión? fueron respondidos con base en la revisión bibliográfica y mediante la caracterización biofísica de las poblaciones conocidas, empleando mapas temáticos de suelo, clima y características geomorfológicas de la península de Yucatán (^{Bautista et al., 2005}; ^{INEGI, 2008}; ^{Inifap-Conabio, 1995}).

Criterio C. Vulnerabilidad biológica intrínseca. Subcriterio C-1), demografía. El número total de individuos, reclutamiento y densodependencia fueron estimados en las 6 poblaciones silvestres caracterizadas de B. pliabilis. Los datos recabados fueron analizados estadísticamente para estimar la densidad poblacional y la relación entre la altura y el diámetro base de las poblaciones. Subcriterio C-2, genética. Para este subcriterio no se cuenta con estudios genéticos por lo que sólo se estimó el nivel de diferenciación entre poblaciones, con base en caracteres morfológicos de las 6 poblaciones caracterizadas. Subcriterio C-3), interacciones bióticas especializadas. La información biológica de B. pliabilis relacionada con su reproducción y propagación fue obtenida de ^{Hernández- Sandoval et al. (2012)}.

Criterio D. Impacto de la actividad humana. Subcriterio D-1), ¿cómo afecta al taxón la alteración antrópica del hábitat?, D-2), ¿cuál es el nivel de impacto de las actividades humanas sobre el hábitat del taxón?, y D-3), ¿existe evidencia que indique un deterioro en la calidad o extensión del hábitat como efecto de cambios globales o se prevé un cambio drástico en el uso del suelo? Para actualizar el nivel de impacto de las actividades humanas sobre el hábitat de B. pliabilis se generaron espacialmente áreas de influencia a distancias de 0.1, 1, 2 y 3 km de cada una de las 6 poblaciones caracterizadas en campo. En las áreas de influencia se estimó la superficie de uso de suelo y la vegetación presente, reportada en cartografía del INEGI Serie II (^{INEGI, 2001}) y Serie VI (^{INEGI, 2016}); también se estimó la distancia a carreteras, densidad de vías de comunicación y distancia a localidades urbanas y rurales. Los subcriterios D-4), ¿cuál es el impacto del uso sobre el taxón? y D-5), ¿el taxón es cultivado o propagado ex situ? Fueron respondidos con base en una revisión bibliográfica.

Índice de Rareza

Criterio A. Características de la distribución geográfica. Subcriterio A-1), extensión de la distribución. Los datos de herbario y de campo confirmaron que Beaucarnea pliabilis es endémica de la provincia biogeográfica península de Yucatán, cuya principal área de distribución se encuentra en los estados mexicanos de Campeche, Quintana Roo y Yucatán, en el Petén guatemalteco y al norte de Belice. En lo que respecta a México, los 3 estados representan 7.2% de la superficie total del país.

Los tipos de selvas donde habita B. pliabilis, en las últimas décadas, han reducido su cobertura en el país y en la península de Yucatán, donde su extensión es menor a 1% con respecto al territorio nacional (tabla 1). Solo en caso de la selva baja caducifolia, se presentó un ligero aumento del 0.001%, correspondiente a la regeneración de la selva en ese periodo. La mayoría de las localidades donde fueron colectados u observados ejemplares de B. pliabilis han sido transformados en uso de suelo de tipo agropecuario y urbano. Actualmente, los remanentes de las poblaciones se restringen a parches de vegetación secundaria arbórea y de selva mediana subperennifolia. En la península de Yucatán, la vegetación secundaria arbórea ocupaba, en 2001, una extensión de 2.63% con respecto a la superficie continental del país y, recientemente (datos del 2016), se estima que ocupa 1.92%. Con esta información, se le asigna un valor de 3 al criterio A1 (tabla 2).

Subcriterio A-2), número de poblaciones o localidades conocidas existentes. Con base en las colecciones de herbarios, la especie se ha distribuido en 17 municipios de la península de Yucatán, que podrían corresponder a alrededor de 15 poblaciones. Pero el cambio de uso de suelo ha impactado las poblaciones y actualmente se pueden encontrar en aproximadamente 10 localidades en la península (fig. 1). La pérdida y transformación de las selvas primarias ha ocasionado la fragmentación y aislamiento de poblaciones de B. pliabilis. Es importante resaltar que es una especie dioica, que sus poblaciones han sido severamente saqueadas y que se desconoce en las poblaciones remanentes la proporción de sexos. Para estimar la distancia entre las poblaciones remanentes, seleccionamos 6 sitios que se encuentran en parches de vegetación secundaria arbórea y en un parche de selva mediana subperennifolia, la distancia más corta entre ellas es de 51.9 km, con una distancia promedio entre las poblaciones de alrededor de 200 km (tabla 3). Dada esta densidad de las poblaciones, que corresponde a 4-8 poblaciones visibles en un mapa a una escala 1:4,000,000, se le asignó un valor de 2 al subcriterio A-2 (tabla 2).

Subcriterio A-3), número de provincias biogeográficas en las que se encuentra el taxón. Históricamente, las poblaciones de esta especie se distribuían ampliamente en 2 de las 19 provincias biogeográficas de México (^{Conabio, 1997}): Yucatán y Petén que representan 6.8% de la superficie del país (fig. 2). El mayor número de recolectas se registró en la selva baja caducifolia y en la selva mediana subcaducifolia de la provincia Yucatán, sobre suelos rocosos y formando asociaciones con especies características de estas selvas como las cactáceas candelabriformes, Actualmente, ya no se encuentran las poblaciones registradas en la selva mediana subcaducifolia, debido a la transformación de su hábitat para uso agropecuario. Por lo anterior, se le asignó un valor de 2 al subcriterio A-3 (tabla 2).

Figura 2 Distribución conocida de Beaucarnea pliabilis en las provincias biogeográficas (^{Conabio, 1997}). 

En cuanto a la representatividad de la distribución de B. pliabilis en México, la especie es endémica del país, por lo que las poblaciones forman parte de una distribución periférica o extralimital. Se le asigna un valor de 0 al subcriterio A-4 (tabla 2). Al sumar los valores de impacto de los subcriterios, obtenemos que para el criterio A, suman 7 de un máximo de 11.

Criterio B. Características del hábitat. Subcriterio B-1), ¿en cuántos tipos de vegetación se presenta? La especie se puede encontrar en 3 tipos de vegetación potencial reconocidas por ^{Rzedowski (1990)}, en los bosques tropical caducifolio, tropical subperennifolio y en el bosque espinoso de la península de Yucatán (fig. 3); estos tipos de vegetación representan una extensión potencial de 4.6%, 0.15% y 4.6%, respectivamente, con respecto al territorio nacional. Históricamente, la distribución de B. pliabilis abarcaba también el bosque tropical subcaducifolio, pero éste ha sido transformado a coberturas de uso agropecuario y urbano. Se le asigna un valor de 1 al subcriterio B-1 (tabla 2).

Figura 3 Distribución de Beaucarnea pliabilis superpuesta al mapa de vegetación potencial de ^{Rzedowski (1990)}. Históricamente, la especie se encontraba en el bosque tropical subcaducifolio, pero actualmente solo quedan individuos aislados. 

Subcriterio B-2), ¿el taxón tiene un hábitat especializado? ^{Miranda (1959)} señaló que era frecuente encontrar a esta especie en 2 variantes de la selva baja decidua (renombrada como selva baja caducifolia en 1963 por este autor), en la selva baja decidua con cactáceas candelabriformes que se extendía en una franja de 10 a 15 km y paralela a la costa, desde el sur de Telchac Puerto hasta Sisal, y en la selva baja decidua con Pseudophoenix sargentii H. Wendl. ex Sarg. que se extendía paralela a la costa, cerca de Puerto Juárez, Quintana Roo y cerca de El Cuyo, Yucatán. En la década pasada, se señaló su presencia de modo preferente en la asociación “guayacán” y de modo dominante en la asociación “despeinada” de selva mediana subcaducifolia, donde comparten la dominancia con especies como el Guaiacum sanctum L., Lysiloma latisiliqua A. Gray ex Sauvalle y Metopium brownei (Jacq.) Urb. (^{Martínez y Galindo-Leal, 2002}).

Las selvas bajas inundables también han sido reportadas como hábitat de B. pliabilis. En una matriz de selva baja inundable de Celestún, Yucatán, se han encontrado individuos de B. pliabilis en altillos con selva baja caducifolia, donde comparten el hábitat con Plumeria obtusa L., Nopalea gaumeri Britton et Rose y Sabal yapa C. Wright ex Becc., y en entre las localidades Sierra Papacal y Chuburná Puerto. En Yucatán se han reportado en parches de selva baja caducifolia localizadas en medio de pastizales y ciperáceas inundables, donde comparte el hábitat con otras especies endémicas como Furcraea cahum Trel. y Nopalea gaumeri Britton et Rose (^{Hernández-Sandoval et al., 2012}). También, se ha registrado de modo notable en la asociación mixta de la vegetación de “Bajos” de Calakmul (^{Martínez y Galindo, 2002}).

Las poblaciones se encuentran en sitios con alto porcentaje de pedregosidad, exclusivamente de roca caliza, en suelos escasos de tipo litosol y de rendzina principalmente, y en paisajes con cierta homogeneidad morfogenética y ambiental (^{Bautista et al., 2005}; ^{Inifap-Conabio, 1995}) (tabla 4). Los sitios se encuentran en climas secos muy cálidos (BS₀) con precipitación anual promedio de 550 mm, en los semi secos muy cálidos (BS₁) con precipitación anual promedio de 650 mm, y en los cálidos subhúmedos (Aw₀, Aw₁) con precipitación anual promedio entre 1,000 y 1,300 mm (tabla 4). La temperatura promedio en todos los sitios oscila entre 24.5 y los 26 °C. Con base en la diversidad de climas, asociaciones vegetales y hábitats donde ha sido reportada B. pliabilis, se considera que no tiene un hábitat especializado y se le asigna un valor de 0 al subcriterio B-2.

Subcriterio B-3), ¿la permanencia de la población es dependiente de un hábitat primario? La especie fue registrada en las selvas de la península de Yucatán, actualmente ya no se encuentra en las localidades reportadas en las colectas de Schott 892 (MO), Lundell 8128 (MO), Gaumer 23520 (MO), Goldman 1901 (MO), Miranda 8234 (MEXU) y Matuda 37487 (MEXU); los remanentes de poblaciones aún existentes sólo se encuentran en los parches de selva aún conservados. De acuerdo con datos históricos, las poblaciones eran más abundantes en sitios conservados (^{Miranda, 1959}), por lo que el número poblacional encontrado, escasez de plántulas y juveniles y ausencia de estructuras reproductivas en los muestreos, indican la dependencia de un hábitat primario, valor 1 punto.

Subcriterio B-4), ¿la permanencia de la población requiere de regímenes de perturbación particulares o está asociada a etapas transitorias en la sucesión? Las poblaciones evaluadas en este trabajo son remanentes de poblaciones alguna vez más abundantes, donde permanecen individuos en parches de vegetación secundaria arbórea y arbustiva de selva baja caducifolia, selva mediana subcaducifolia y de selva mediana subperennifolia, donde comparten hábitat con otras especies endémicas (tabla 5). No se cuenta con información que indique una asociación con regímenes de perturbación transitorios, valor 0 puntos (tabla 2).

Subcriterio 5), amplitud del intervalo altitudinal que ocupa el taxón, la especie se distribuye en la península de Yucatán, la cual es predominantemente plana, en la porción central y al sur puede alcanzar una elevación de hasta 380 m snm; la población de XMA es la que se encuentra a mayor altitud, en una zona de lomeríos menores a 200 m snm, valor 3 puntos (tabla 2; ^{Lugo, 1999}). Los valores de impacto para este criterio suman 5 puntos (5/9).

Criterio C. Vulnerabilidad biológica intrínseca. Subcriterio C-1.1), demografía. En los transectos de las 6 poblaciones muestreadas se encontró un total de 223 individuos de diferentes etapas de desarrollo (fig. 4), de los cuales 28% fueron adultos, 19% juveniles y 53% plántulas. En las poblaciones, excepto en CHI, se encontraron individuos de por lo menos 2 etapas de desarrollo, pero la proporción fue variable. La densidad varió de 0 a 118 plántulas/0.1ha, 0 a 43 juveniles/0.1ha y 0 a 18 adultos/0.1ha. El cálculo del número total de individuos que realizamos con base en los censos y registros de poblaciones es de 501-5,000 individuos, por lo que le asignamos a este subcriterio (C-1.1) el valor de 2 puntos (tabla 2).

Figura 4 Distribución de las etapas de desarrollo en las poblaciones de Beaucarnea pliabilis. 

El reclutamiento, tomando como evidencia la presencia de plántulas, se registró en 2 transectos, aunque la mayor densidad se registró en la población de SP (117 plántulas de un total de 118). También, la densidad de juveniles fue relativamente baja en todas las poblaciones, excepto para la población de FCP donde se registraron 32 individuos. Aunque, la densidad de adultos fue variable, 3 poblaciones (FCP, RL y XMB) presentaron una densidad similar (tabla 6). En la población de CHI estuvieron ausentes las plántulas y los adultos, y sólo se encontró un individuo de la etapa juvenil. Por lo tanto, consideramos para el subcriterio C-1.2 que hay reclutamiento en algunas poblaciones, lo que corresponde a un valor de 2 puntos (tabla 2).

En cuanto a los atributos demográficos, se consideraron las siguientes preguntas (tabla 2). a) ¿Hay evidencia de denso-dependencia en la reproducción? no se contó con evidencia de densodependencia en la reproducción, valor 0 puntos. b) ¿Hay clonalidad? la especie no tiene la capacidad de generar nuevos individuos independientes por medio de reproducción asexual, valor 1 punto. c) ¿Hay evidencia de decrecimiento de las poblaciones en el país? las observaciones históricas señalan poblaciones alguna vez más densas en la península de Yucatán (^{Miranda, 1959}), patrón que se puede observar en los sitios más conservados como el de FCP, o en especies hermanas del género Beaucarnea que localmente, y en sitios conservados son abundantes (^{Castillo-Gómez, 2011}; ^{Hernández-Sandoval et al., 2012}; ^{Pérez-Farrera et al., 2012}; ^{Rubio Méndez et al., 2018}), valor 1 punto. d) ¿Hay evidencia de una varianza muy grande en la fecundidad? de acuerdo con nuestras observaciones y experiencia en campo, la floración año con año es muy variable; en las ciudades de la península de Yucatán se han presentado años con floraciones masivas y años con floraciones escasas; en condiciones silvestres ha sido raro encontrar individuos en etapa reproductiva, tal es el caso de la presencia de un gran número de plántulas en la localidad de SP, donde las plántulas se encontraban cerca de un solo individuo adulto progenitor, valor 1. e) ¿El taxón es dioico, los individuos son dicógamos o autoincompatibles? la especie presenta sexos separados, inflorescencias con flores pistiladas o estaminadas (^{Hernández-Sandoval, 1993}), valor 1 punto. f) ¿La floración es sincrónica o gregaria? las poblaciones silvestres florecen de abril a julio (^{Cabrera-Cano et al., 1982}; ^{López-Segovia, 1985}), en la actualidad es raro encontrar individuos silvestres en floración, en las ciudades florecen de manera sincrónica, valor 1 punto; g) ¿El taxón produce pocos propágulos? los individuos producen de 1 a 6 inflorescencias, y cada individuo puede producir miles de frutos provistos por una sola semilla (Orellana, datos sin publicar), las cuales tienen una alta tasa de germinación (^{López-Segovia, 1985}), valor 0 puntos.

Subcriterio C-2. Genética. No se cuenta con datos genéticos para B. pliabilis, por lo que los subcriterios 1) cantidad de variación genética molecular, 2) estructura genética molecular y 3) variación molecular, recibieron un valor de 0 (tabla 2). En cuanto a la variabilidad morfológica entre poblaciones, se observaron diferencias entre poblaciones de las selvas bajas y las selvas medianas. En la población con mayor precipitación, en la selva mediana de FCP, se encontraron adultos con alturas máximas de hasta 10.5 m, mientras que las demás poblaciones no superan, en general, los 8 m. En el caso del área basal, esta fue mayor en las poblaciones de las selvas bajas caducifolias del norte de Yucatán, SP, DZE y RL, con valores entre 0.67 y 4.33 m² e incluso un valor extremo en 15.5 m², en comparación con la de las poblaciones de las selvas medianas, FCP y XMA, cuyos valores se ubicaron entre 0.02 y 1.0 m2 (fig. 5). Las poblaciones de las selvas medianas son más numerosas, sin embargo, presentan un área basal por hectárea menor al de las poblaciones de las selvas bajas caducifolias. Dada esta variabilidad entre las poblaciones, se le asignó el valor de 1 punto al criterio C-2.4 (tabla 2).

Figura 5 Relación observada y estimada entre la altura (m) y el área basal (m²) de los individuos de las poblaciones de Beaucarnea pliabilis registrados en los transectos de las localidades DZE (Dzemul), FCP (Felipe Carrillo Puerto), RL (Río Lagartos), SP (Sierra Papacal) y XMA (Xmabén), estimada mediante un modelo de regresión múltiple (r² = 0.27; p < 0.0000). 

Subcriterio C-3. Interacciones bióticas especializadas. Subcriterio 1), ¿el taxón requiere una “nodriza” para su establecimiento? las observaciones en campo indican que las plántulas requieren en sus primeras etapas la sombra de otros árboles, como ocurre con otras especies del género, valor 1 punto (^{Cardel et al., 1997}; ^{Hernández-Sandoval et al., 2012}). 2), ¿El taxón requiere un hospedero o forofito específico? la especie no requiere forofito, valor 0 puntos. 3), ¿El taxón requiere un polinizador específico? la polinización en la especie es cruzada y mediante vectores bióticos o abióticos como el viento, valor 0 puntos. 4), ¿El taxón tiene un dispersor específico? las semillas son dispersadas principalmente por el viento, valor 0 puntos. 5), ¿El taxón presenta mirmecofilia obligada? no, valor 0 puntos. 6), ¿El taxón presenta dependencia estricta de la micorriza? en Beaucarnea gracilis se ha encontrado asociación con micorrizas (^{Camargo-Ricalde et al., 2003}), para B. pliabilis no contamos con esta información, pero en ninguno de los casos hay evidencia de que la dependencia sea estricta, valor 0 puntos. 7), ¿El taxón sufre una afectación importante por depredadores, patógenos? en las poblaciones silvestres se observaron individuos con daños en tronco y ramas (presentan oquedades y rupturas del tronco principal original), también se encontraron excretas de ganado en las poblaciones que colindan con ranchos ganaderos, y se ha observado depredación de plántulas por ganado, valor 1 punto. Los valores de impacto para este criterio suman 10 puntos (10/23, tabla 2).

Criterio D. Impacto de la actividad humana. Subcriterios D-1), ¿cómo afecta al taxón la alteración antrópica del hábitat? D-2) ¿Cuál es el nivel de impacto de las actividades humanas sobre el hábitat del taxón? D-3) ¿Existe evidencia que indique un deterioro en la calidad o extensión del hábitat como efecto de cambios globales o se prevé un cambio drástico en el uso del suelo? en la década de los noventa, las poblaciones analizadas, en un área de 100 m², se encontraban en una matriz de parches de vegetación secundaria, de tipo agropecuario y de selva baja caducifolia; una década después, la selva fue reemplazada por la cobertura de tipo agropecuario (fig. 6). Conforme nos alejamos de las poblaciones muestreadas, a distancia de 1, 2 y 3 km, se registra una reducción en la superficie ocupada por la selva baja caducifolia y por la vegetación secundaria, en tanto que aumenta la superficie de la cobertura de uso agropecuario, se registran la presencia de las coberturas de uso urbano y zonas desprovistas de vegetación. También, se registran las coberturas manglar (por la cercanía de la población de la localidad de río Lagartos al mar) y selva mediana subperennifolia a distancias de 2 y 3 km². El incremento de la cobertura de uso agropecuario y la pérdida de selva baja caducifolia en las inmediaciones de las poblaciones muestran una tendencia hacia la pérdida de hábitat.

Figura 6 Cobertura de uso de suelo y vegetación presentes alrededor de las poblaciones de Beaucarnea pliabilis en radios de 0.1, 1, 2 y 3 km utilizando las series II (a) y V (b) del INEGI. SBC: Selva baja caducifolia, SMSC: selva mediana subperennifolia, VS: vegetación secundaria, AGP: agropecuario, ZU: zona urbana, MG/PH: manglar/pastizal halófilo, DV: desprovisto de vegetación, CA: cuerpos de agua. 

El impacto negativo de las actividades humanas en el área de distribución de B. pliabilis se ve reforzado por la presencia de vías de comunicación, asentamientos rurales y urbanos en las inmediaciones de sus poblaciones silvestres. Cinco de 6 poblaciones de B. pliabilis se encuentran a menos de 200 m de la carretera más cercana en línea recta; las poblaciones de DZE, SP y XMA son las más cercanas a alguna carretera, están a menos de 50 m; en tanto que las de RL y FCP están a menos de 116 y 177 m de la carretera más cercana, respectivamente (fig. 7a). La densidad de vías de comunicación (calles, terracería y carreteras) en un área de 100 m² alrededor de las poblaciones silvestres muestreadas es relativamente alta, 31.7 m por m², en tanto que la densidad en las superficies de 1, 2 y 3 km², es entre 160 y 170.8 m por m² (fig. 7c).

Figura 7 Distancia de las localidades de muestreo de Beaucarnea pliabilis a las vías de comunicación, asentamientos rurales y urbanos. a) Distancia a carreteras, b) distancia a asentamientos rurales y urbanos. Localidades CHI (Chiquilá), DZE (Dzemul), FCP (Felipe Carrillo Puerto), RL (Río Lagartos), SP (Sierra Papacal) y XMA (Xmabén). 

La cercanía a asentamientos rurales y urbanos también es un factor que contribuye a la pérdida de hábitat, sitios como DZE, FCP, RL, y XMA se encuentran a menos de 1 km de algún asentamiento (fig. 7b).

Con base en esta información podemos evaluar los subcriterios, D-1, la pérdida de cobertura vegetal de la península de Yucatán ha ocasionado la pérdida y/o reducido el tamaño de las poblaciones, las cuales son relativamente más abundantes en sitios mejor conservados como en las localidades de SP y FCP, valor 1 punto. Subcriterio D-2), el impacto de las actividades humanas es fuerte y afecta a todas las poblaciones (figs. 1, 6, 7), valor 3 puntos. Subcriterio D-3), en la tabla 1 podemos ver la reducción en las selvas que habita B. pliabilis (y un muy pequeño aumento en la selva baja caducifolia) del 2001 al 2016, y en la figura 6 los cambios en el uso de suelo, en fechas recientes. Lo cual da sustento a considerar una fuerte pérdida de hábitat por cambio de uso de suelo, debido a actividades como la ganadería, la agricultura comercial mecanizada, el crecimiento urbano y de infraestructura impulsados por el crecimiento demográfico y económico de la región, valor 1 punto.

En el subcriterio D-4) ¿cuál es el impacto del uso sobre el taxón? la especie tiene un alto valor ornamental, por lo que sus poblaciones silvestres han sufrido una extracción indiscriminada, hecho que condujo a su inclusión en la lista de especies amenazadas de la NOM-059 desde el año 1994, el impacto de uso es fuerte en algunas o moderado en todas las poblaciones, valor 2 puntos. 5) ¿La especie ha sido cultivada y propagada ex situ? su inclusión en la NOM-059 en 1994 promovió la creación de unidades de manejo de vida silvestre (UMA), aunque en el 2011 se reportaron hasta 34 UMA que propagaban la especie, muchas de estas no continúan a la fecha; no obstante, es frecuente encontrarla en viveros de las principales ciudades, valor -1 punto. Los valores de impacto para este criterio suman 6 puntos (6/10, tabla 2).

Finalmente, los criterios A, B, C del índice de rareza suman 1.718, y el criterio D del índice de impacto antropogénico tuvo un valor total de 0.6, dando una suma total del MER de 2.312 (tabla 2) lo cual ubica a la especie en la categoría de en peligro de extinción (P), por tener un valor igual o mayor a 2.