Área de estudio

El estudio se llevó a cabo en las áreas excluidas (EPE) y con pastoreo (EPP) del rancho “El Palmar”, del municipio de Charcas, SLP, así como en el rancho “Los Garza” (LG), Charcas, SLP, y el ejido “El Tepetate” (ET) del municipio de Santo Domingo, SLP. Al iniciar la instalación del parque eólico los propietarios del rancho El Palmar decidieron excluir temporalmente al pastoreo de ganado en sus 1300 ha, y por ello constituye una excelente oportunidad para evaluar el impacto de los diferentes sistemas de producción ganadera en la condición de un mismo sitio de agostadero con base en la composición florística, densidad y producción de biomasa y la sección transversal al pie, por forma vital y valor forrajero de las especies perennes presentes, además de la cobertura actual del suelo. El Palmar colinda al oeste con el ejido “El Tepetate” del municipio de Santo Domingo, SLP, en el que desde su creación en 1932 (^{PHINA, 2024}) como en todos los ejidos y comunidades agrarias se practica el pastoreo irrestricto y continuo por bovinos, ovinos, caprinos y equinos, tanto del propio ejido como de otros ejidos en sus agostaderos comunales que normalmente carecen de cercos divisorios; asimismo, colinda al norte con el rancho Los Garza de propiedad privada y con antecedente de aprovechamiento mediante el pastoreo exclusivo de ovinos por más de 60 años. Por su parte, dentro del propio rancho el Palmar en 2019 sus propietarios decidieron reiniciar la explotación mediante el pastoreo de bovinos para carne en un predio de 1000 ha, y dejaron 300 ha completamente excluidas. El disponer de áreas excluidas al pastoreo en México es poco frecuente, por lo que contar con dicha área, es una excelente oportunidad para contrastar el impacto de cuatro sistemas de manejo completamente distintos en la condición actual de agostadero del zacatal de Bouteloua gracilis, así como en su grado de riesgo a erosión mediante la evaluación de la vegetación y la cobertura actual del suelo.

Las cuatro unidades de producción colindan entre sí y se ubican en una meseta a casi 2450 m de altitud. (Figura 1). Su clima es seco templado (BS0kw), con temperatura media anual entre 12 y 18 °C, con abril como el mes más caliente y seco; junio es el mes más húmedo con una precipitación media de 88 mm; además, registra lluvias de verano (300 mm) y porcentaje de precipitación invernal entre 5 y 10.2%; la estación seca abarca de noviembre a mayo (^{INEGI, 2017}). El tipo de vegetación observado en las cuatro unidades de producción correspondió a un zacatal con predominio de Bouteloua gracilis (willd. Ex kunth) lag. Ex griffiths, acompañado principalmente de Aristida adsencionis L., B. scorpioides. Lag., Erioneuron avenaceum (Humb., Bonpl. & Kunth) Tateoka, Muhlenbergia villiflora Hitchc., Buddleja scordioides Kunth, Dalea bicolor Humb. & Bonpl. ex Willd., Cylindropuntia imbricata (Haw.) F.M. Knuth, Flourencia sernua DC., Larrea tridentata (Sessé & Moc. ex DC.) Cov., Opuntia rastrera F.A.C. Weber, O. cantabrigensis Lynch, Parthenium incanum Kunth, Talinopsis frutescens A. Gray y Yucca filifera hort EX Engelm.

Figura 1: Localización geográfica del área de estudio. 

Figure 1: Geographic location of the study area.  

Caracterización de la vegetación

Para caracterizar la vegetación se utilizó primero, una imagen de satélite LANDSAT de febrero de 2021, a la cual mediante ARC-GIS v. 2010 (^{ESRI, 2011}), se le realizó un proceso de clasificación no supervisada el cual generó un índice de vegetación con cuatro clases de vegetación presentes en el área de estudio; posteriormente, con el apoyo de un GPS Garmin Montana 680 de alta precisión, se acudió a campo para realizar la supervisión correspondiente y cotejar las cuatro variantes de vegetación identificadas previamente. En cada una de las unidades de producción, se recolectaron todas las especies perennes; las especies anuales se descartaron por carecer de efectos significativos en la condición del agostadero, pues su presencia temporal es de apenas unas pocas semanas (^{Aldrete y Aguirre, 1982}). Los especímenes recolectados se identificaron y depositaron en el Herbario Isidro Palacios (SLPM) del Instituto de Investigación de Zonas Desérticas de la Universidad Autónoma de San Luis Potosí.

Mediciones estructurales y funcionales de la vegetación

Con base en recorridos de campo, y el plano generado por la clasificación no supervisada, entre agosto y septiembre de 2021 (segunda mitad del periodo de lluvias o de crecimiento), se ubicaron cuatro áreas de muestreo, una en cada unidad de producción en estudio. Luego, en cada unidad de producción se realizó el levantamiento de cuatro transectos de 100 m con un total de 80 puntos de muestreo, 320 cuadrantes y 1280 plantas medidas (Figura 2). Se aplicó el método sin parcela de cuadrantes centrados en un punto sobre transectos desarrollado por ^{Cottam y Curtis (1956)}, con las modificaciones propuestas por ^{Aldrete y Aguirre (1982)} para evaluar vegetación de diversos estratos como ocurre en los matorrales áridos. Las especies registradas con el muestreo se clasificaron según su valor forrajero y su reacción sucesional al pastoreo, como decrecientes, crecientes e invasoras, acorde con ^{Vallentine (2000)} y ^{Bolaños y Aguirre (2000)}.

Figura 2: Ubicación de los puntos de muestreo en las cuatro áreas de estudio.  

Figure 2: Location of sampling points in the four study areas. 

Para el cálculo de la densidad absoluta de cada especie estudiada se empleó la fórmula utilizada por ^{Aldrete y Aguirre (1982)} para determinar el recíproco del cuadrado de la distancia media de dicha especie, con la fórmula siguiente:

En donde la unidad de área es la superficie de referencia para expresar la densidad, en las mismas unidades que el área media de las especies.

Por otra parte, se utilizó la fórmula del cono truncado invertido para estimar la biomasa volumétrica instantánea de cada individuo:

Donde: π = 3.1416; h = altura o distancia entre los dos radios; R = radio de la copa; r = radio basal.

Diámetro basal

El diámetro basal de cada una de las plantas medidas se determinó con una cinta diamétrica (Forestry supliers, Mod. 283D) al momento de realizar las estimaciones de densidad y biomasa; para el caso de las herbáceas y zacates se midió todo el diámetro de la planta o del macollo, mientras en las arbustivas y arborescentes fue la base del fuste.

Estado de la superficie del suelo

Para evaluar el estado actual de la superficie del suelo en cada unidad de producción y con ello el grado de riesgo inherente a la erosión y posterior desertificación, se utilizaron cinco Líneas ^{Canfield de intercepción (Canfield, 1941)} de 5 m de longitud, colocadas de forma paralela y alterna al transecto principal, georreferenciando sus puntos de inicio y final para facilitar estudios posteriores; sobre estas líneas se midieron las porciones de suelo desnudo y las cubiertas con vegetación, mantillo, piedras o heces. El valor medio de cada tipo de cobertura fue expresado en centímetros.

Análisis estadístico

Se realizó la comparación de la cobertura media actual del suelo (suelo desnudo, cubierto con vegetación, con mantillo, con piedras y con heces), así como la correspondiente a las secciones transversales al pie obtenidas para cada una de las unidades de producción y clasificadas de acuerdo con su valor forrajero, mediante un diseño experimental completamente aleatorio con el procedimiento GLM de SAS (^{SAS Institute, 1999}). Asimismo, a las medias aritméticas de las secciones transversales al pie de las plantas evaluadas por el muestreo se les realizó un ANOVA unidireccional, con un diseño experimental completamente aleatorio y, un nivel de significancia de (P ≤ 0.05), después en ambos casos se hizo una comparación de medias por el método de Tukey (P ≤ 0.05), (^{Tukey, 1949}).

Finalmente, mediante una matriz de 16 × 16 se realizó la ordenación y clasificación mediante los módulos DECORANA y TWINSPAN del programa PC-ORD v 7.0. (^{McCune y Mefford, 2018}). La matriz incluyó las 16 variables de vegetación y cobertura de suelo resultantes de la evaluación de la vegetación de las cuatro unidades de producción, así como los 16 transectos utilizados en dicha evaluación.

Sitio de agostadero

Los agostaderos son los sitios de alimentación del ganado mediante vegetación espontánea, además, sostienen a la fauna silvestre y son fuente de bienes y servicios ecosistémicos. ^{Holechek et al. (2011)}, describieron que un sitio de agostadero es una unidad de paisaje cuyo potencial de producción está determinado por su clima, suelo y topografía, que lo normal es que haya heterogeneidad y por ello existan muchos sitios de agostadero distintos. Derivado de ello, en la presente investigación y con la información generada por el proceso de clasificación no supervisada y los recorridos de campo, se pudo comprobar que en las cuatro unidades de producción se tenía en realidad un único sitio de agostadero compartido, esto es, con el mismo clima, mismo suelo e igual topografía, y se le identificó como zacatal de Bouteloua gracilis. Con base en ello se pudo determinar que las diferencias mostradas por el índice de vegetación generado de inicio con el proceso de clasificación no supervisada aplicado a la imagen de satélite, resultaron ser más bien distintas condiciones del agostadero, debidas al aprovechamiento pecuario histórico al que el zacatal ha estado sometido (Figuras 3 y 4).

Figura 3: Aspecto de la colindancia entre el ejido El Tepetate (A), rancho Los Garza (B) y rancho El Palmar excluido (C), en el zacatal de Bouteloua gracilis. Las diferencias en la vegetación son resultado del aprovechamiento individual al que cada unidad de producción ha estado sometida. Imagen capturada el 15 de septiembre de 2021.  

Figure 3: Boundary aspect between the common land Ejido El Tepetate (A), Los Garza Ranch (B) and El Palmar Ranch (C), excluded in the Bouteloua gracilis grassland. The differences in vegetation are the results of the individual use to which each production unit has been subjected. Image captured in September 15, 2021. 

Figura 4: Aspecto aéreo de la colindancia entre el ejido El Tepetate ( A), rancho Los Garza (B ) y rancho El Palmar excluido (C). Imagen capturada el 02 de octubre de 2021. 

Figure 4: Aerial aspect of the boundary between Ejido El Tepetate (A), Los Garza Ranch (B) and El Palmar ranch (C) excluded. Image captured on October 2, 2021. 

Composición florística

La mayor riqueza florística de especies perennes medidas por el muestreo se encontró en la porción pastoreada del rancho El Palmar (EPP) con 22 especies, seguida por su porción excluida al pastoreo (EPE) con 21; por su parte en el ejido El Tepetate (ET) hubo 18 especies y finalmente en el rancho Los Garza (LG) se registraron 17 (Cuadro 1). Además, para cada unidad de producción se registraron otras 35 especies perennes respectivamente, pero fuera de las unidades muestreadas con el transecto. La clasificación del valor forrajero de las plantas presentes en cada unidad de producción arrojó que la unidad de producción con mayor proporción de especies deseables fue EPP con 52.4%, seguida de EPE y ET con 50% respectivamente; por su parte, LG tuvo el menor porcentaje de especies deseables con 46.7%, por lo que no se registraron diferencias significativas entre las cuatro unidades de producción (P = 0.058); en contraste, respecto a las especies forrajeras indeseables sí se registraron diferencias significativas (P < 0.003), pues ET registró la proporción más alta con 27.8%, luego LG con 20%, mientras EPP tuvo 19.0% y por último EPE con 10 por ciento.

La composición florística total de las especies perennes registrada en el zacatal de B. gracilis para cada unidad de producción de este estudio fue superior a lo descrito por ^{Lara-Juárez, Castillo, Tristán, Rendón y Aguirre (2016)} en estudios florísticos de especies perennes, en Charcas SLP para dos unidades de producción colindantes, el ejido Francisco I. Madero con 22 especies y el potrero “siete vueltas” del rancho Laguna Seca con 35 especies muy probablemente porqué aquellos sitios de agostadero correspondieron a matorrales semiáridos con menor cantidad de precipitación media anual y además en el caso del ejido con pastoreo irrestricto por diversas especies animales y ganado incluso de otros ejidos, mientras que el rancho privado tampoco cuenta con cercos divisorios por sitio de agostadero que obliguen al ganado a pastorear cada sitio en forma ordenada, lo que demuestra que las malas prácticas no son exclusivas del usufructo comunal, sino que ocurren también por falta de conocimiento en los ranchos privados. Por otra parte, los resultados de ambas porciones del rancho El Palmar fueron muy similares a las 55 especies perennes registradas en las evaluaciones de 18 parcelas ganaderas de tres sitios de agostadero del ejido El Castañón, Catorce, SLP ^{Aguirre-Rivera, Negrete y Castro (2020)}., en donde a partir de su parcelación de los agostaderos comunales en 1993, ahora existe la posibilidad de controlar la carga animal, mayor número de aguajes, saladeros y de cercos divisorios para distribuir y uniformar el pastoreo.

Densidad y biomasa

La densidad florística en general fue más alta en los ranchos con mayor deterioro aparente (Cuadro 2), en contraste, la producción de biomasa estimada para cada unidad de producción fue superior hasta cinco a seis veces para los ranchos mejor conservados (Cuadro 3). Un efecto más ocasionado por el pastoreo diferenciado en las cuatro unidades de producción fue que diversas plantas preferidas y muy consumidas por el ganado ya no existen en ET y LG, las dos unidades de producción en peor condición visual aparente lo que evidencia su poca resistencia a la defoliación intensa y frecuente que ahí se permite.

Respecto a la densidad de las diversas plantas forrajeras deseables evaluadas en las cuatro UP, las correspondientes a ET y LG, fueron mayores, sin embargo, esto ocurre por sus diámetros y vigor disminuidos que ocupan poco espacio en el terreno; en contraste, sus evaluaciones de biomasa volumétrica instantánea y de la sección transversal al pie fueron significativamente menores con respecto a las de las mismas plantas en EPP y EPE. Lo anterior coincide con ^{Abate, Ebro y Nigatu (2012)} en agostaderos áridos del sur de Etiopía quienes sugieren que se genera una tendencia regresiva en la condición de los agostaderos conforme se da el deterioro de las unidades de producción por sobrepastoreo.

Las densidades relativas de especies forrajeras deseables registradas en ET, EPP y LG del presente estudio (99.9, 98.8 y 94.9%) fueron en general superiores a las registradas por ^{Aguirre-Rivera et al. (2020)} en seis parcelas ganaderas con buena condición de agostadero de tres sitios de matorrales áridos, aprovechados como agostaderos comunales por más de 50 años, pero ya bajo control privado en los recientes 25 años en el ejido. El Castañón, Catorce, SLP, (97.8, 93.6 y 90.9%) posiblemente porqué en ese ejido se registra menor precipitación media anual (300 mm), y aún quedan remanentes del maltrato sufrido en los agostaderos por el pastoreo comunal al que fue sometido. Por otra parte, la menor densidad de especies forrajeras deseables del presente estudio se registró en el área excluida del rancho El Palmar (85.2%), y fue causada principalmente por la propagación de O. rastrera, especie considerada como forrajera menos deseable tras 12 años de exclusión del ganado, lapso en el que los efectos inherentes del pastoreo como el pisoteo y el consumo ocasional de dicha planta no se han registrado; esto coincide con lo señalado por ^{Flores y Yeaton (2000)} en cuanto a la dinámica estructural de los zacatales áridos, que por sucesión son sustituidos por especies suculentas leñosas.}

Con respecto a la cantidad de biomasa volumétrica instantánea de las especies forrajeras deseables de EPP (7.6 mil m³ ha^-1) y EPE (5.9 mil m³ ha^-1) éstas fueron muy superiores a las registradas en LG (1.3 mil m³ ha^-1) y ET (1.2 mil m³ ha^-1); sin embargo, la biomasa de especies forrajeras deseables de ambas porciones de rancho El Palmar, fueron incluso poco inferiores a lo reportado por ^{Aguirre-Rivera et al. (2020)} en agostaderos correspondientes al ecotono entre el matorral desértico micrófilo-matorral desértico rosetófilo en buena condición (8.2 mil m³ ha^-1) del ejido El Castañón y Anexos, Catorce, San Luis Potosí, aunque en ese caso mucha de esa biomasa fue aportada principalmente por arbustivas forrajeras y no por zacates como en el presente estudio. Así, hay que considerar que en las áreas que dejan de ser defoliadas por mucho tiempo como en el caso de EPE con 12 años sin presencia de ganado, se registra una reducción en la producción de biomasa, pues las plantas dejan de producir al no tener el estímulo de la defoliación. En cuanto a la cantidad producida de biomasa volumétrica instantánea de las especies forrajeras deseables de LG y ET, las dos unidades de producción con la peor condición de agostadero de este estudio, coincide con lo señalado por ^{Schlesinger et al. (1996)}, respecto a la disminución de la producción de forraje de buena calidad en sitios de agostadero con intensidad de pastoreo pesada y con condición de agostadero pobre, en el sur de Nuevo México.

Sección transversal al pie

Ésta es considerada una medida más confiable para determinar la producción de materia viva que la correspondiente a la biomasa volumétrica instantánea, ya que ésta permanece en el terreno a pesar de que la planta haya sido parcial o totalmente defoliada. Los mayores valores registrados para el promedio de este atributo en las plantas deseables y menos deseables correspondieron a las dos unidades de producción visualmente más conservadas (EPE y EPP), en tanto que para el de las plantas indeseables los mayores promedios los registraron LG y ET, las dos unidades visualmente más deterioradas por el sobrepastoreo (Cuadro 4). Las diferencias encontradas en las medias aritméticas de las secciones transversales al pie para las especies deseables (P = 0.001) y de las especies indeseable (P = 0.02) de las cuatro unidades de producción resultaron significativas entre las cuatro unidades de producción, no así en lo correspondiente a las especies menos deseables (P > 0.05). Este comportamiento coincide con lo mencionado por ^{Abate et al. (2012)} en agostaderos del sureste de Etiopía, que señalan que conforme incrementa la carga animal, la sección transversal al pie de las especies forrajeras deseables disminuye y la condición de agostadero se vuelve regresiva, tal y como ocurre en la presente investigación, pues en ET y LG, la sección transversal al pie de sus plantas forrajeras deseables está notoriamente disminuida respecto a las dos áreas del rancho El Palmar, debido principalmente al desequilibrio registrado entre la capacidad de carga y la carga animal, además de los otros factores negativos del manejo del pastoreo como ocurre principalmente en ET donde se registra pastoreo sin restricciones por diferentes especies animales, e incluso animales de otros ejidos, ya que no cuenta con cercos que delimiten y dividan los agostaderos y que puedan obligar a los animales a consumir el forraje de manera uniforme y ordenada.

Evaluación de la cobertura actual de la superficie del suelo

En este rubro, EPE fue la unidad de producción con la más alta estabilidad potencial del sitio al registrar la mayor proporción media de suelo cubierto con vegetación, y la menor proporción de suelo desnudo; por otro lado, ET registró la mayor proporción media de suelo desnudo y la menor proporción media de suelo cubierto con vegetación, lo que conlleva menor estabilidad potencial del sitio, mayor grado de riesgo a erosión y de desertificación de las cuatro UP estudiadas (Cuadro 5). Estadísticamente las diferencias en las proporciones medias de suelo desnudo (P < 0.0001) y, cubierto con vegetación (P < 0.0002), fueron significativas entre las cuatro unidades de producción en estudio; no así las correspondientes a suelo cubierto con mantillo (P = 0.2344), con heces (P = 0.0948) y con rocas (P = 0.1445). ^{Parker (1951)} y luego ^{Holechek et al. (2011)} reconocieron que la porción de suelo desnudo de un sitio de agostadero nos permite determinar su grado de riesgo a erosión y el inherente peligro de desertificación; este punto de vista se considera también en los nuevos enfoques para la determinación de la condición de sitios, pues permite establecer con respecto a la cantidad de suelo desnudo, sí la condición de sitio es satisfactoria o insatisfactoria con base en su estabilidad potencial (^{Smith, 2003}). entre más alta sea la porción con suelo desnudo, mayor es la posibilidad de erosión durante las tormentas que dispersan las partículas ya que no existe barrera física ni biológica que las proteja, seguida de la formación casi inmediata de la escorrentía y el posterior arrastre de sedimentos fértiles fuera de la cuenca, aunque es posible que en las zonas bajas se encuentren elementos como piedras, heces y vegetación que pueden incrementar la rugosidad del suelo y con ello disminuir la velocidad del escurrimiento superficial.

En cuanto a la cobertura actual del suelo en las cuatro UP estudiadas, las proporciones medias de suelo cubierto por vegetación de este estudio fueron superiores a 55 puntos, las cuales pueden considerarse buenas y con menor grado de riesgo a erosión y desertificación en el corto plazo, pues resultaron superiores a las registradas con igual tamaño de muestra (cinco líneas Canfield de 500 cm por ^*transecto de 100 m) por ^{Lara et al. (2016)}, en un estudio en tres sitios de agostadero con buena, regular y pobre condición visual del rancho Laguna Seca, Charcas, SLP (con 47, 44 y 42 puntos porcentuales), así como las del ejido Francisco I. Madero Charcas, SLP (con 12, 14 y 12 puntos porcentuales respectivamente). Por otra parte, también superan al 47% de cobertura de suelo de un sitio de vegetación identificado como matorral desértico micrófilo con muchos izotes y al 28% de otro sitio de vegetación identificado como matorral desértico micrófilo con pocos izotes, registrados por ^{Aguirre-Rivera et al. (2020)} en las parcelas ganaderas del ejido El Castañón, Catorce, SLP, el cual es un ejido administrado desde hace 25 años nuevamente como propiedad privada; por su parte, el porcentaje de suelo cubierto con vegetación de ET fue similar al 54% registrado en seis parcelas de un sitio de vegetación denominado ecotono entre el matorral desértico micrófilo y matorral desértico rosetófilo, del mismo estudio en Castañón, sólo que a diferencia de ET que al parecer está en deterioro, el sitio de agostadero del Castañón parece encontrarse en recuperación tras 25 años de manejarse como pequeña propiedad.

Análisis multivariable de las variables evaluadas

El módulo DECORANA del programa PC-ORD v 7.0 (^{McCune y Mefford, 2018}), utilizado para la ordenación de las 16 variables estudiadas indicó que sólo los primeros dos ejes proporcionaron información significativa para la interpretación con un valor de varianza acumulada de 0.5981 (Figura 5). Así, en el eje 1 se observó un primer gradiente correspondiente a la condición de agostadero de buena a regular, con la biomasa y densidad de especies forrajeras deseables en el extremo izquierdo asociadas a la buena condición del agostadero; en tanto que en el extremo derecho y como indicadoras de peor condición, se presenta la biomasa de especies forrajeras indeseables y menos deseables, la proporción de suelo cubierto con piedras y la proporción de suelo desnudo. En el eje 2 por su parte, se puede inferir un gradiente correspondiente a la intensidad de pastoreo que va de pesada a ligera; así, en la parte inferior se tuvo la mayor carga animal con la presencia de las variables correspondientes a la biomasa de especies forrajeras menos deseables, la proporción de suelo desnudo, la densidad de especies indeseables, la densidad de especies forrajeras menos deseables y la proporción de suelo cubierto con piedras; por su parte, en la parte superior se registró la menor cobertura de heces en el suelo y el diámetro basal de especies forrajeras deseables, variables que indican menor carga animal y por ende menor intensidad y frecuencia de defoliación. Dichos gradientes encontrados en la ordenación de las variables y transectos de nuestro estudio, coinciden con el comportamiento de los gradientes derivados de la ordenación de variables realizadas por ^{Aguirre-Rivera et al. (2020)} y Aguirre-Rivera, Negrete y Castro (2023), de sus trabajos de investigación realizados primero en parcelas seleccionadas del ejido El Castañón y Anexos, Catorce, San Luis Potosí, para determinar el efecto de la parcelación de los agostaderos comunales, en sus parcelas ganaderas administradas de forma privada entre 1993 y 2013, y luego en la implementación de una propuesta actualizada para evaluar agostaderos de tres ranchos privados de zonas áridas del altiplano potosino-zacatecano.

Figura 5: Ordenación bidimensional de las variables evaluadas en el estudio.  

Figure 5: Bi-dimensional ordering of the variables evaluated in the study.. 

Asimismo, mediante el módulo TWINSPAM se realizó la clasificación de los 16 transectos evaluados los cuales formaron seis grupos de entre uno y hasta cuatro transectos (Figura 6). La variable que generó la dicotomía en el primer nivel de clasificación fue la proporción de suelo cubierto con heces, la cual separó a los grupos 1 (tres transectos) de los grupos 2, 3, 4, 5 y 6 (13 transectos). El grupo 1 se conformó por los transectos 2, 3 y 4 de EPE en condición de agostadero excelente, sin carga animal, ni intensidad, ni frecuencia de pastoreo al encontrarse excluidos de ganado. En el segundo nivel de clasificación las variables indicadoras de una nueva partición fueron el diámetro basal de las especies forrajeras deseables y el número total de especies de plantas presentes que ocasionaron que los grupos 2, 3 y 4 se separaran de los grupos 5 y 6; en el primero se incluyeron los transectos 1, 2, 3 y 4 de EPP, el 1 de EPE y el 2 de ET con condición de agostadero buena, con base en los diámetros de sus especies forrajeras deseables y la mayor diversidad de plantas presentes; por su parte en el segundo grupo quedaron los transectos 1, 2, 3 y 4 de LG, y 1, 3 y 4 de ET, debido posiblemente a su peor condición de agostadero. Posteriormente en el tercer nivel de clasificación los grupos 2 y 3 que incluyó a los transectos 1, 2, 3 y 4 de EPP y el 1 de EPE se separaron del grupo 4 conformado sólo por el transecto 2 de ET debido a la variable indicadora correspondiente al diámetro basal de especies forrajeras indeseables, por lo que la diferencia se originó en la condición de agostadero buena de los grupos 2 y 3, respecto a la condición de agostadero regular del grupo 4; en este mismo nivel de clasificación también se dividieron los grupos 5 y 6 todos ellos correspondientes tanto al rancho pastoreado por ovinos durante más de 60 años y como al ejido con pastoreo comunal irrestricto, divididos por la variable indicadora correspondiente a la biomasa de especies deseables, por lo que la diferencia entre ellos, correspondió a la condición de agostadero regular del grupo 5 respecto a la condición de agostadero pobre de los tres transectos que conformaron el grupo 6. Por último, en el cuarto nivel de clasificación la variable indicadora correspondiente a la proporción de suelo cubierto con piedras separó a los grupos 2 y 3, el primero incluyó a los transectos 2, 3 y 4 de EPP, mientras el otro a los transectos 1 de EPE y el 1 de EPP, la diferencia fue posiblemente la condición de agostadero buena a excelente de los transectos del grupo 2, respecto a la condición de agostadero buena de los transectos del grupo 3. La formación de los grupos derivados de la clasificación coincide con el comportamiento en la conformación de grupos registrada en la investigación realizada en los agostaderos de tres ranchos privados productores de bovinos para carne en el altiplano potosino-zacatecano por ^{Aguirre-Rivera et al. (2023)} en su propuesta para la implementación de una nueva metodología para evaluar agostaderos de matorrales semiáridos.

Figura 6: Clasificación de los transectos evaluados en las cuatro unidades de producción.  

Figure 6: Classification of the transects evaluated in the four-production unit.