Área de estudio

El estudio comprendió las ANP estatales La Trinidad, La Hediondilla y Llano La Soledad, con una extensión de 33.02 km², 43.73 km² y 76.07 km², respectivamente, ubicadas en el municipio de Galeana, Nuevo León (Figura 1). El clima es seco con una precipitación anual de 300 mm (^{Guadarrama et al., 2003a}, ^2003b, ^2003c). Las ZSCE representan un mosaico biológico de más de 200 especies de vertebrados y 284 especies de plantas (^{Conabio, 2008}; ^{Guadarrama et al., 2003a}, ^2003b, ^2003c). El pastizal halófito, dominado por Muhlenbergia villiflora es la principal comunidad vegetal y está rodeado de matorral xerófilo-micrófilo, con Larrea tridentata y Flourencia cernua como componentes principales (^{Estrada et al., 2010}).

Fuente: Elaboración propia con imagen de relieve obtenida de la plataforma de Google Earth (22 de junio, 2021).

Figura 1 Ubicación y zonas de manejo de las zonas sujetas a conservación ecológica en el municipio de Galeana, Nuevo León. 1) Zona núcleo de uso restringido, 2) Zona de amortiguamiento de aprovechamiento sustentable de agroecosistemas, 3) Zona de amortiguamiento de aprovechamiento sustentable de los recursos naturales. 

Cada área protegida cuenta con una zona núcleo de uso restringido (ZN) que busca mantener las condiciones naturales del ecosistema. Además, La Hediondilla y La Trinidad contemplan dos zonas de manejo (Figura 1): zona de amortiguamiento de aprovechamiento sustentable de agroecosistemas (ZASA) y zona de amortiguamiento de aprovechamiento sustentable de los recursos naturales (ZARN), orientadas a regular las actividades agropecuarias (^{Guadarrama et al., 2003a}, ^2003b, ^2003c; ^{Semarnat, 2004}).

Evaluación de la salud del ecosistema

Se evaluó mensualmente (de febrero de 2010 a enero de 2011) la condición de la salud del ecosistema de las tres ZSCE. De forma sistemática, se establecieron transectos y parcelas permanentes en las que se seleccionaron y evaluaron ocho indicadores que determinan los atributos Estabilidad del sitio/suelo, Funcionalidad hidrológica e Integridad biótica (Tabla 1), de acuerdo con el método cualitativo sugerido por ^{Pellant et al. (2005)}. Los datos de cada indicador se obtuvieron a través de la observación del grado de diferencia de la descripción del sitio ecológico, mediante las categorías descriptivas extrema (E), moderada a extrema (Mo), moderada (M), ligera a moderada (L) y nula a ligera (N) descritas para cada indicador, donde nula a ligera representa una categoría sin deterioro, mientras que extrema significa mayor deterioro (^{Pellant et al., 2000}).

Considerando el tamaño de la superficie de cada área protegida, el número de transectos y parcelas varió de la siguiente manera: en La Trinidad se establecieron 30 parcelas sobre cinco transectos de 2.6 km cada uno, en La Hediondilla 32 parcelas en cuatro transectos de 4 km cada uno y en el Llano La Soledad 60 parcelas en cinco transectos de 6 km cada uno. En cada transecto se establecieron parcelas permanentes de 10 m × 10 m (con una equidistancia de 500 m), y dentro de estas se ubicaron mensualmente al azar dos subparcelas de 1.0 m²; en las primeras se midieron los indicadores que requirieron de un horizonte amplio para su apreciación (indicadores 1 y 5) y dentro de las subparcelas el resto de los indicadores que demandaron mayor exactitud en su percepción (indicadores del 2 al 8). La mayoría de las parcelas se situaron dentro del área de pastizal halófito, sin embargo, algunas quedaron en la zona de matorral xerófilo-micrófilo (donde predominan las plantas arbustivas). Por consiguiente, se interpretó de forma separada la condición de la salud de ambos ecosistemas.

Producción de biomasa

Cada mes se estimó la producción de biomasa aérea en las tres ZSCE. De acuerdo con las características de la vegetación de las parcelas, se calculó la producción de biomasa del pastizal dentro de las subparcelas de 1.0 m² mediante el corte de la vegetación a 3 cm sobre la superficie del suelo. Asimismo, se estimó la biomasa de arbustos en las parcelas de 100 m² y la biomasa de vegetación invasora dentro de las subparcelas de 1.0 m², ambos con el método no destructivo de Adelaide, que consta de una muestra (rama) que representa la especie de interés y se emplea para medir el número de veces que la muestra de referencia cabe en dicha especie (^{Andrew et al., 1979}; ^{Andrew et al., 1981}). Todas las muestras vegetales se depositaron en bolsas de papel, posteriormente, fueron secadas en estufa (Felisa FE-293D) a una temperatura de 75 °C durante 72 horas y, finalmente, se pesaron en una balanza digital (Sartorius LC 620 S) para determinar su biomasa.

Análisis estadísticos

Con los datos categóricos de cada indicador se generaron tablas de frecuencia con las que se estableció el estado funcional (categoría) que mejor describió cada atributo, esto con base en la elección de las categorías descriptivas que se observaron con mayor frecuencia durante la evaluación de los indicadores (^{Pellant et al., 2005}). Esta evaluación incluyó el análisis de todas las observaciones recabadas durante la medición de cada atributo (^{Pellant et al., 2005}) e información recopilada de los programas de manejo de cada ANP (^{Guadarrama et al., 2003a}, ^2003b, ^2003c), lo que permitió llegar al criterio que definió la condición de salud del ecosistema de cada zona de manejo en las tres ZSCE.

Para comprobar la normalidad de los datos de la producción de biomasa aérea se utilizó la prueba de Kolmogorov-Smirnov, y para cada tipo de biomasa se calculó la estadística descriptiva (^{Zar, 2010}). Por medio de gráficos se visualizó la producción promedio mensual de biomasa aérea entre las zonas de manejo de cada área evaluada y entre ANP, esto con base en las condiciones hídricas del año evaluado a partir de los datos de precipitación mensual acumulada, obtenidos de las estaciones climatológicas 19182-San Roberto y 5136-Las Hormigas de la Comisión Nacional del Agua (^{Conagua, 2020}). Mediante el análisis de Kruskall-Wallis se evaluó si existieron diferencias en la producción mensual de biomasa dentro de cada área protegida y en la producción anual de biomasa entre ellas (^{Zar, 2010}). Para determinar cuál de las ANP fue diferente con respecto a su producción anual, se realizó un análisis post hoc de comparación múltiple, con la función “kruskalmc” del paquete “pgirmess” (^{Giraudoux, 2022}) para el lenguaje de programación R. Los análisis estadísticos se realizaron con el software Statistica versión 7 (^{StatSoft, Inc., 2004}) y el software R versión 4.0.2 (^{R Core Team, 2020}).

Salud del ecosistema

En La Trinidad, la condición del pastizal correspondió a la categoría de nula a ligera en ambas zonas de manejo (Tabla 2). En contraste, la ZN presentó el estado funcional moderado en los dos primeros atributos (Tabla 2); asimismo, los tres atributos obtuvieron la categoría de nula a ligera, lo que indicó que dentro de esta misma zona de manejo existieron dos tipos de condición de salud definidos por la Estabilidad del sito/suelo y la Funcionalidad hidrológica del pastizal. Situación similar ocurrió para el área de matorral (que únicamente estuvo presente en la ZASA), ya que presentó dos estados funcionales en los atributos S y H; no obstante, de manera general, la condición de su salud fue extrema (Tabla 2).

La salud del pastizal de la ZN de La Hediondilla fue extrema para los atributos S y H, en tanto que la salud del matorral (solo presente en la ZN) tuvo el estado de extremo en todos sus atributos (Tabla 2). El resto de las zonas de La Hediondilla presentaron una condición de salud sin deterioro (Nula a ligera). En el Llano La Soledad, los tres atributos determinaron una salud del ecosistema de pastizal sin deterioro (Nula a ligera), mientras que la evaluación del matorral indicó una condición de salud extrema en cuanto a todos los atributos (Tabla 2).

Producción de biomasa

La Trinidad presentó una producción media anual de biomasa total de 0.40 Mg/ha^-1 (Tabla 3). La ZN tuvo una mayor producción promedio en julio con 0.44 Mg/ha^-1 (Figura 2), en tanto que la ZASA presentó dos picos de producción: uno en julio (3.09 Mg/ha^-1) y otro en noviembre (0.97 Mg/ha^-1), además fue la zona de manejo con mayor producción anual de biomasa, promediando 0.65 Mg/ha^-1 (Tabla 3). La ZARN mostró mayor producción de biomasa en el mes de julio con una media de 0.74 Mg/ha^-1 (Figura 2). De acuerdo con las diferentes comunidades vegetales presentes en La Trinidad, la biomasa de plantas invasoras fue la de mayor producción anual dentro del área (2.01 Mg/ha^-1) (Tabla 3). Con relación a zonas de manejo, la ZN de La Trinidad registró un promedio anual de 0.20 Mg/ha^-1 de biomasa de pastizal, en tanto que la ZASA tuvo una producción anual de plantas invasoras de 2.56 Mg/ha^-1 (Tabla 3). La ZARN produjo una media anual de biomasa de pastizal de 0.98 Mg/ha^-1, seguida de 0.89 Mg/ha^-1 de biomasa de vegetación invasora (Tabla 3).

Nota. Datos de precipitación de la estación climatológica 19182-San Roberto, Galeana, Nuevo León. ZN: zona núcleo de uso restringido; ZASA: zona de amortiguamiento de aprovechamiento sustentable de agroecosistemas; ZARN: zona de amortiguamiento de aprovechamiento sustentable de los recursos naturales.

Fuente: Elaboración propia.

Figura 2 Produccción mensual de biomasa de las zonas de manejo en la zona sujeta a conservación ecológica La Trinidad. 

El ANP La Hediondilla registró una producción promedio anual de biomasa total de 0.53 Mg/ha^-1 (Tabla 3). Con respecto a la producción media mensual, la ZN alcanzó en enero 0.44 Mg/ha^-1, en tanto que la ZARN mostró en julio su máxima producción mensual con 0.22 Mg/ha^-1 (Figura 3). La ZASA destacó con tres picos de producción de biomasa aérea, promediando tanto en julio como en diciembre 2.5 Mg/ha^-1 y en octubre 1.05 Mg/ha^-1 (Figura. 3); por lo tanto, fue la zona de manejo que produjo mayor cantidad de biomasa aérea con una media anual de 0.88 Mg/ha^-1 (Tabla 3). Considerando las diferentes comunidades vegetales en La Hediondilla, la producción anual de plantas invasoras fue la mayor durante el año evaluado (2.23 Mg/ha^-1) (Tabla 3). En cuanto a sus zonas de manejo, la ZN tuvo una mayor producción promedio anual de 0.15 Mg/ha^-1 de biomasa aérea de pastizal, en tanto que la ZASA y ZARN produjeron la mayor cantidad anual de biomasa de vegetación invasora con una media de 2.42 Mg/ha^-1 y 0.51 Mg/ha^-1, respectivamente (Tabla 3).

Nota. Datos de precipitación de la estación climatológica 5136-Las Hormigas, Saltillo, Coahuila. ZN: zona núcleo de uso restringido; ZASA: zona de amortiguamiento de aprovechamiento sustentable de agroecosistemas; ZARN: zona de amortiguamiento de aprovechamiento sustentable de los recursos naturales.

Fuente: Elaboración propia.

Figura 3 Produccción mensual de biomasa de las zonas de manejo en la zona sujeta a conservación ecológica La Hediondilla. 

La producción promedio anual de biomasa total en el Llano La Soledad fue de 0.11 Mg/ha^-1 (Tabla 3). Durante los meses evaluados se presentaron diversos incrementos de biomasa aérea: en marzo con una media de 0.17 Mg/ha^-1, en julio 0.21 Mg/ha^-1 y en diciembre 0.16 Mg/ha^-1 (Figura 4). En esta área protegida se registraron únicamente biomasa de arbustos y de pastizal, de las cuales destacó la segunda con 0.14 Mg/ha^-1 de producción media anual (Tabla 3).

Nota. Datos de precipitación de la estación climatológica 5136-Las Hormigas, Saltillo, Coahuila. ZN: zona núcleo de uso restringido.

Fuente: Elaboración propia.

Figura 4 Produccción mensual de biomasa en la zona sujeta a conservación ecológica Llano La Soledad. 

La producción anual de biomasa difirió entre áreas protegidas (H _{(2, N = 2322)} = 23.62, p < 0.001). De acuerdo al análisis post hoc de comparación múltiple, el ANP que registró la mayor producción anual de biomasa total fue La Hediondilla (Tabla 3). Para cada área de estudio la producción mensual de biomasa fue diferente entre los meses de muestreo (La Trinidad: H _{(11, N = 630)} = 46.08, p < 0.001; La Hediondilla: H _{(11, N = 612)} = 38.32, p = 0.001; Llano La Soledad: H _{(11, N = 1080)} = 104.88, p < 0.001). En las tres áreas protegidas, julio fue el mes con el mayor promedio de productividad de biomasa aérea: La Trinidad 1.82 Mg/ha^-1, La Hediondilla 1.57 Mg/ha^-1 y Llano La Soledad 0.21 Mg/ha^-1 (Figura 5). Particularmente, La Trinidad presentó en noviembre un segundo pico de producción con una media de 0.60 Mg/ha^-1, mientras que para La Hediondilla y el Llano La Soledad un segundo incremento considerable ocurrió en diciembre con 1.50 Mg/ha^-1 y 0.16 Mg/ha^-1, respectivamente (Figura 5).

Fuente: Elaboración propia.

Figura 5 Produccción mensual de biomasa en la zona sujeta a conservación ecológica La Trinidad (T), La Hediondilla (H) y Llano La Soledad (S). 

Salud del ecosistema

Las dos zonas de manejo de La Trinidad en el área de pastizal no presentaron ninguna afectación en su salud. Por el contrario, en la ZN se identificó un moderado deterioro en los atributos funcionales Estabilidad del sitio/suelo y Funcionalidad hidrológica, a pesar de ser un sitio de uso restringido. Esto puede ser consecuencia de la heterogeneidad de la zona a causa de los cambios que ha sufrido su cobertura y composición florística debido a la agricultura (^{Hernández, 2010}; ^{Yen, 2006}), que comenzó a intensificarse en la década de 1990 (^{Zaragoza-Quintana et al., 2012}), además del sobrepastoreo (^{Guadarrama et al., 2003b}; ^{Hernández, 2010}) y la creación de nuevos caminos (^{Guadarrama et al., 2003b}). ^{Ochoa (2006)} observó resultados similares a los de la ZN de La Trinidad, ya que una evaluación de la salud del pastizal realizada en colonias de Cynomys ludovicianus, en general, mostró un estado de degradación moderado a causa de la presencia o ausencia que C. ludovicianus tienen sobre el pastizal, así como por la actividad ganadera que, en Janos, Chihuahua, es significativamente mayor a la del resto del país.

De forma conjunta, la ZN de La Trinidad presentó un estado funcional sin alteración en su Estabilidad del sitio/suelo y Funcionalidad hidrológica. Obtener esta condición en ambos atributos indica que dentro de la ZN se encontraron parcelas con escasa afectación en los recursos del suelo (incluyendo nutrientes y material orgánico) y en su dinámica hídrica (almacenamiento y liberación de agua) (^{Pellant et al., 2005}; ^{Printz et al., 2014}). Con respecto a la salud del matorral de esta ANP, esta reflejó un estado extremadamente modificado debido a que también ha sufrido el impacto agrícola al ser transformado a cultivos de papa y alfalfa (^{Estrada et al., 2010}; ^{Hernández, 2010}). Simultáneamente, sobre este tipo de vegetación ocurrió el apacentamiento de ganado caprino, bovino y equino, que junto con la agricultura ocasionaron entre 1971 y 2003 una reducción del 35.6% de su superficie original en La Trinidad (^{Hernández, 2010}), lo que en conjunto pudo afectar severamente la salud de la ZSCE. No obstante, el matorral también mostró una condición de salud prácticamente sin deterioro en los atributos S y H, lo que indica la presencia de sitios que mantuvieron tanto las propiedades del suelo como su capacidad hidrológica (^{Pellant et al., 2005}; ^{Printz et al., 2014}).

En el área protegida La Hediondilla, la ZN se caracterizó por presentar un estado crítico o deterioro extremo en la salud de sus ecosistemas (pastizal y matorral); el resto de las zonas de manejo no mostraron deterioro. Históricamente, La Hediondilla es uno de los pastizales halófilos del Altiplano Mexicano que ha perdido gran parte de su cobertura vegetal original (^{Estrada et al., 2010}; ^{Yen, 2006}). Sin embargo, a pesar de que dentro de su ZN no se observó ningún cambio de uso de suelo, al menos hasta el año 2011 (^{Zaragoza-Quintana et al., 2012}), los resultados del presente estudio indican que existen algunos factores que podrían estar repercutiendo considerablemente en su salud. Por ejemplo, la superficie que comprende esta zona de uso restringido es mucho menor (menos del 10%) con respecto a la extensión de las otras zonas de manejo (^{Guadarrama et al., 2003c}). Además, la ZN ha sido rodeada por caminos y parcelas circulares para uso agrícola, algunas de las cuales posteriormente se recolonizaron con vegetación invasora (^{Zaragoza-Quintana, 2011}). Estos disturbios podrían afectar la resistencia y resiliencia de la zona a través de efectos indirectos, ya que, de acuerdo con ^{Duniway et al. (2010)}, la salud de los pastizales puede afectarse negativamente ante la presencia de carreteras y caminos en áreas adyacentes, alterando así la función hidrológica y la estabilidad del suelo y del sitio, lo que ocasionaría un impacto en la función ecológica del ecosistema.

El atributo Integridad biótica no presentó afectación alguna en el estado funcional de la zona núcleo de La Hediondilla y de La Trinidad, lo cual indica que en ambas zonas se mantuvieron tanto las características funcionales del pastizal como la estructura de sus comunidades en un contexto de variabilidad normal y una capacidad de resiliencia presente en los sitios (^{Blair et al., 2014}; ^{Pellant et al., 2005}). Caso contrario ocurre en un pastizal mediano abierto de Zacatecas, en el que este atributo es el más deteriorado a causa del sobrepastoreo que ocasiona la disminución de cobertura vegetal, además, el ganado pastorea libremente, repercutiendo en la producción de materia seca y en la función productiva de ese ecosistema (^{Echavarría-Cháirez et al., 2015}).

El Llano La Soledad es el área protegida que en los últimos 40 años ha presentado menor cambio de uso de suelo y pérdida de cobertura vegetal nativa (^{Estrada et al., 2010}; ^{Zaragoza-Quintana et al., 2012}), debido a que toda la superficie de la ZSCE es de uso restringido, lo que limita las actividades agrícolas en su interior (^{Guadarrama et al., 2003a}). Esta escasa actividad pudo influir sobre la salud del pastizal, ya que no presentó ningún tipo de deterioro. Sin embargo, el área de matorral reflejó una salud extremadamente alterada. Al respecto, en el Llano La Soledad se tiene registro de que en este tipo de vegetación la riqueza de especies es alta, la cual está relacionada con los disturbios que ocasiona principalmente el sobrepastoreo de ganado caprino, vacuno y equino, que ha sido una actividad común dentro de la ANP (^{Estrada et al., 2010}; ^{González et al., 2017}; ^{Yen, 2006}). La presencia de este tipo de ganado podría estar relacionada con la deteriorada salud que mostró esta comunidad vegetal en el presente estudio. Los indicadores que resultaron con las categorías extremas en el matorral (tales como área de suelo desnudo, compactación del suelo y fragmentación de la costra biológica) están asociados a los impactos negativos que generan las actividades ganaderas (^{Duniway et al., 2010}; ^{Guzmán-Aranda et al., 2012}; ^{Menegaz & Arellano, 2016}).

A partir del presente estudio se pudo detectar que específicamente las ZN de La Hediondilla y La Trinidad manifestaron importantes alteraciones en la salud del hábitat del perrito llanero mexicano, por lo que se considera apropiada una revisión y actualización de sus planes de manejo, ya que aparentemente estas zonas no están cumpliendo con el objetivo de mantener las condiciones naturales del ecosistema y su conservación a corto, mediano y largo plazo (^{Guadarrama et al., 2003b}, ^2003c). Particularmente, se sugiere como estrategia reconsiderar un aumento de superficie de la ZN de La Hediondilla o la eliminación de caminos que cruzan a través de esta, así como la implementación de programas sobre el uso del suelo de las áreas circundantes para mejorar la salud del ANP. Respecto a La Trinidad, debe atenderse con especial énfasis aquellos atributos funcionales que definieron tanto el deterioro de la ZN como de la ZASA, por ejemplo, mediante protocolos y acciones a corto plazo (como el manejo del ganado) que contribuyan en la mejora y mantenimiento de la calidad del sitio, con el fin de disminuir el deterioro de estos pastizales y matorrales, los cuales están sufriendo una rápida transformación con preocupantes consecuencias sobre sus funciones ecológicas (^{Echavarría-Cháirez & Rubio-Aguirre, 2020}; ^{Souther et al., 2019}).

Debido a que las tres ANP albergan las colonias más extensas que existen del perrito llanero mexicano en Nuevo León (^{Guadarrama et al., 2003a}, ^2003b, ^2003c), es sustancial atender las sugerencias señaladas en beneficio de C. mexicanus y asegurar su papel ecológico que como especie clave representa para la comunidad del pastizal. Un hábitat de buena calidad es necesario para la presencia del perrito llanero mexicano y en un momento determinado para favorecer el incremento de sus poblaciones, así como de las especies de fauna relacionadas con la distribución de este mamífero, tales como la zorra norteña Vulpes macrotis, el tejón Taxidea taxus, el águila real Aquila chrysaetos y el gorrión de Worthen Spizella wortheni (^{Conabio, 2008}; ^{Contreras-Balderas et al., 2008}).

Producción de biomasa

La producción de biomasa de La Trinidad y La Hediondilla estuvo integrada de plantas invasoras, pastizal y matorral natural. Particularmente, la ZN de La Trinidad destacó por presentar únicamente biomasa de pastizal; en contraste, la ZASA, aparte de ser la zona que registró la mayor producción anual de biomasa, tuvo los mayores promedios anuales de biomasa de plantas arbustivas e invasoras. Esto se debió a que en su interior se ha desarrollado con más intensidad la agricultura (^{Hernández, 2010}; ^{Zaragoza-Quintana, 2011}), lo que ha generado la formación de considerables parches de vegetación invasora, algunos presentes desde el año 2000 (^{Zaragoza-Quintana et al., 2012}). La ZARN tuvo el mayor promedio anual de biomasa de pastizal debido a que es una zona de aprovechamiento con menos actividades antrópicas, lo que ha dejado una mayor superficie de pastizal nativo (^{Zaragoza-Quintana, 2011}).

La Hediondilla también produjo grandes cantidades de biomasa de plantas invasoras, principalmente la ZASA, la cual además reportó el mayor promedio anual de biomasa de pastizal. Por tanto, fue la zona con más productividad primaria debido, en parte, al tamaño de su superficie y a que en esta se ha concentrado gran actividad agrícola, lo que ha ocasionado la presencia de numerosos sitios cubiertos de vegetación invasora y otros en proceso de recuperación natural (^{Yen, 2006}; ^{Zaragoza-Quintana et al., 2012}). No obstante que la ZARN también cuenta con diversos fragmentos en recuperación y con vegetación invasora, se presentó una menor producción de biomasa en el pastizal que puede estar relacionada a la intensidad del pastoreo en la zona (^{Zaragoza-Quintana, 2011}). Con relación a la ZN (la cual es de tamaño considerablemente menor en comparación con las otras zonas de manejo), esta tuvo el mayor promedio anual de biomasa arbustiva y estuvo ausente de plantas invasoras, lo cual puede ser consecuencia de la restricción de actividades antrópicas en su interior (^{Guadarrama et al., 2003c}).

La productividad primaria fue diferente entre las tres ANP. La mayor producción anual de biomasa fue en La Hediondilla, seguida de La Trinidad y el Llano La Soledad, lo cual es consecuencia de las condiciones particulares de cada ZSCE. Por ejemplo, La Trinidad y La Hediondilla han sido pastizales sometidos a un fuerte impacto antropogénico por las actividades agrícolas y la apertura de caminos (^{Estrada et al., 2010}; ^{Zaragoza-Quintana et al., 2012}), lo que ocasionó que su producción de biomasa estuviera fuertemente ligada con los cambios en sus coberturas. Especialmente La Hediondilla, además de ser el sitio que resultó con la salud más deteriorada, su productividad primaria estuvo influenciada por la presencia de sitios con vegetación invasora, que son resultado del cambio de uso de suelo que caracteriza a esta área protegida.

^{Medrano (2010)} y ^{Hernández-Gómez et al. (2013)} mencionan que la máxima productividad de biomasa que registraron en La Hediondilla y en el Llano La Soledad fue entre los meses de julio y septiembre, de acuerdo con la temporada de lluvias de la región. Por otro lado, en el presente trabajo la vegetación reflejó el máximo incremento de producción en el mes de julio, es decir, poco después del paso del huracán “Alex” en el estado de Nuevo León, por lo que en general la dinámica de la biomasa en las tres ANP parece relacionarse con la precipitación del año evaluado (^{Conagua, 2020}).

Las actividades antrópicas en el Llano La Soledad son escasas, por lo que la productividad primaria que se encontró en esta zona es similar a los reportados para esta zona protegida. El promedio anual de biomasa obtenido para la comunidad del pastizal del Llano La Soledad (0.14 Mg/ha^-1) fue similar a los registrados en otros trabajos desarrollados en la misma área (0.11 Mg/ha^-1 en ^{Medrano, 2010}; 0.16 Mg/ha^-1 y 0.13 Mg/ha^-1 en ^{Hernández-Gómez et al., 2013}), a pesar de que hubo diferencias en los tamaños de muestra y tratamientos de las parcelas evaluadas en cada estudio. En el caso de La Trinidad y La Hediodilla, la productividad media anual obtenida en este trabajo (0.49 Mg/ha^-1 y 0.19 Mg/ha^-1, respectivamente) fue superior al estimado por ^{Medrano (2010)} para esta última área protegida (0.05 Mg/ha^-1). Pese a que el esfuerzo de muestreo entre ambos estudios fue similar, esta diferencia puede deberse a las características del fragmento evaluado por ^{Medrano (2010)}, el cual era frecuentemente utilizado como agostadero para ganado caprino. En cambio, para el presente trabajo, la distribución de los puntos de muestreo al azar dentro de todas las zonas de La Trinidad y La Hediodilla permitió una mayor representación de la heterogeneidad de la comunidad vegetal de cada ANP.

Dado que la ganadería caprina es importante en las áreas de distribución del perrito llanero mexicano y la competencia por alimento con este es baja (^{Mellado et al., 2005}), se considera a partir de los presentes resultados que el uso de suelo debe ser enfocado (al menos en las ANP destinadas a esta especie) a impulsar la ganadería bajo un marco de pastoreo sostenible. Adicionalmente, la implementación de un monitoreo de los pastizales que incluya la cercanía a zonas urbanas (^{Hopson et al., 2015}), parámetros biofísicos (productividad secundaria, fertilidad del suelo, ciclo de nutrientes, temperatura acumulada) y el uso de herramientas tecnológicas como sensores e imágenes de detección remota (^{Xu & Guo, 2015}), podrían ser de utilidad para comprender mejor la dinámica de las poblaciones de C. mexicanus y el manejo del ecosistema que habita.