Área de estudio

La RBSAT se ubica en la porción media Noroeste de la Gran Sierra Plegada o provincia fisiográfica de la Sierra Madre Oriental, en los municipios Ciudad Valles y Tamuín en el estado de San Luis Potosí. Cubre una superficie de 21 464 ha que se extienden entre las coordenadas 22°05’ y 22°24’ N y 98°52’ y 99°01’ O (^{Semarnap, 1994}) (Figura 1). Su área de influencia la conforman nueve localidades: Laguna del Mante, Los Sabinos Número Dos, Gustavo Garmendia (La Unión), León García, El Sabino del Obispo, El Aguaje, Las Palmas, Estación Tamuín y Fracción la Tima (^{Durán, 2018}).

Figura 1 Reserva de la Biosfera Sierra del Abra Tanchipa y su zona de influencia. 

El clima predominante en la región corresponde al cálido subhúmedo con lluvias en verano Aw1(x’) con una precipitación media anual de 965 mm y una temperatura promedio anual de 24.5 °C (^{Inegi, 2002}; ^{García, 2004}). La Reserva se localiza en la región hidrológica 26 (RH26) Río Panuco, correspondiente a la cuenca del Río Tamuín, aunque las corrientes superficiales son escasas, posee ríos subterráneos conectados a través de diversos sótanos propiciados por el fracturamiento e infiltraciones de agua en las rocas calcáreas (^{Inegi, 2002}). Los suelos predominantes corresponden a los Feozems y Leptosoles (^{Inegi, 1983}).

La vegetación dominante corresponde a la selva baja caducifolia y mediana subcaducifolia. La primera se caracteriza por tener individuos en el estrato arbóreo con alturas de 4 a 15 m que pierden sus hojas en la temporada seca y la segunda por alcanzar alturas de 25 a 30 m con un follaje denso (^{Rzedowski, 2006}; ^{Miranda y Hernández, 2016}).

Distribución actual

Se llevó a cabo la consulta con taxónomos y expertos en vegetación, así como una revisión y verificación de los registros de Quercus oleoides para el estado de San Luis Potosí, en el herbario Isidro Palacios (SLPM) del Instituto de Investigación de Zonas Desérticas (IIZD) de la Universidad Autónoma de San Luis Potosí. Se corroboraron los registros de la especie en la RBSAT y se georreferenciaron los sitios de colecta.

Con estos datos y con el apoyo de guías locales y expertos conocedores de la flora regional, se realizaron siete recorridos de campo en la Reserva y su área de influencia para verificar la presencia de la especie y obtener información complementaria de los sitios donde se desarrolla. Adicionalmente, se aplicaron 15 entrevistas semiestructuradas a algunos habitantes de los ejidos Laguna del Mante, Gustavo Garmendia y Los Sabinos, para conocer los usos principales de la especie, antiguas zonas de distribución, su importancia para los lugareños y acciones que podrían ser implementadas para su conservación.

Los puntos identificados durante los recorridos de trabajo de campo, así como los datos de los registros de colectas fueron representados en el Sistema de Información Geográfica Arcgis 10.3., para realizar un análisis espacial y determinar las principales características ambientales de los sitios donde se desarrolla la especie. Para ello se consideró: la altitud, la pendiente, la orientación de ladera, el tipo de suelo, la geología y los usos de la tierra. Estas capas de información con una resolución de 30 m de tamaño de píxel fueron consultadas del servicio digital de Inegi.

Distribución potencial

Para generar el modelo se usaron siete registros geográficos de las especies procedentes del herbario, así como los datos de 56 sitios georreferenciados en campo. Las principales variables ambientales consideradas en el análisis fueron: i) temperatura media anual; ii) precipitación anual (1970-2010), obtenidas de la base de datos Worldclim (^{Fick y Hijmans, 2017}); iii) geología; iv) edafología (^{Inegi, 2007}); v) altitud; y vi) pendiente, generadas a partir del modelo digital de elevación del terreno (CEM 3.0) (^{Inegi, 2013}). Para estimar las contribuciones de cada variable, así como sus respectivos parámetros, se utilizó el algoritmo de máxima entropía del programa Maxent 3.4.1.

Durante la ejecución de Maxent se determinaron 35 puntos aleatorios de entrenamiento de sensibilidad máxima. La evaluación y contribución de cada variable ambiental a la predicción se basó en la estimación de los porcentajes de contribución relativa de Jackknife que definen las variables más importantes para el modelo (^{Phillips et al., 2006}). Asimismo, se utilizaron la prueba de evaluación de curva- característica de funcionamiento del receptor (ROC) y el área bajo la curva (AUC), que establecen la capacidad de discriminación del modelo y los valores de aptitud y rendimiento (^{Hanley y McNeil, 1982}). Los valores del modelo se interpretan como insuficientes en un intervalo de 0.50-0.70; pobre de 0.71-0.80; bueno de 0.81-0.90 y excelente de 0.91-1.0 (^{Araujo y Guisan, 2006}).

Áreas con aptitud para la conservación

Se realizó una evaluación espacial multicriterio (EMC) por lo que se seleccionaron variables continuas y categóricas como elementos de medición y evaluación. Lo anterior implicó la multiplicación del valor de cada criterio estandarizado a una escala común y un rango de ponderación, para contar con alternativas múltiples a partir de los principios previos (^{Malczewski, 1999}; ^{Eastman, 2012}).

Las variables elegidas para el análisis fueron agrupadas en tres componentes: ambiental, de presión antrópica y subsidios y rezago social (Figura 2) y se basaron en los criterios por ^{Moffett y Sarkar (2006)}, ^{Fabbro et al. (2014)} y ^{Hernández (2014)}. Éstas fueron: conectividad, pendiente, elevación, distancia a la Reserva, uso de suelo, susceptibilidad a incendios, grado de degradación, distancia a zonas agropecuarias (^{Inegi, 2017}), distancia a carreteras, distancia a localidades (^{Inegi, 2011}), subsidios (pago por servicios ambientales, PROCAMPO y PROGAN) e índice de rezago social.

Figura 2 Esquema conceptual utilizado para la evaluación de áreas prioritarias para la conservación. 

En el módulo de evaluación multi-criterio, los componentes se integraron en una sumatoria lineal ponderada de acuerdo a la siguiente fórmula:

Donde:

S = Adecuación de la celda

W _i = Peso asignado al factor i

X _i = Valor del factor i

La estimación de los pesos y la ponderación de cada criterio implicó la consulta a expertos, a personal técnico de la Dirección de la RBSAT y a ejidatarios. Se aplicaron quince matrices de opinión y valor agregado para cada criterio. Los datos obtenidos se analizaron y se comprobaron mediante el Analytic Hierarchy Process (AHP), que consiste en ponderar valores directos o comparaciones por pares e inconsistencias (^{Saaty, 2008}). Los valores representados en una escala numérica (1 a 9) de acuerdo a su importancia, se incorporaron a una matriz de comparación (^{Eastman, 2012}). Los pesos se integraron y calcularon en el módulo Decision support Weight, para promediar los resultados de todos los criterios en un nivel de consistencia menor a 0.1 (^{Clark Labs, 2015}).

La estandarización consistió en trasformar cada criterio en unidades comparables en el módulo Fuzzy a escala byte (0-255), en la que 0 representa el valor mínimo de la alternativa y 255 el máximo para la conservación. Los valores de la matriz de comparación (AHP) derivados de cada cuestionario se consensuaron para calcular los valores ponderados de importancia relativa. A partir de ello, se derivó una composición cartográfica de valores <44 no óptimos y > 195 muy óptimos para la conservación. Posteriormente la información fue estandarizada a cuatro niveles de importancia (no apto, bajo, medio y alto).

Distribución y uso de Quercus oleoides

En la RBSAT y su área de influencia se identificaron 56 sitios con Quercus oleoides, principalmente en los ejidos Laguna del Mante, Gustavo Garmendia y Los Sabinos Dos (Figura 3). En dichos espacios predominan los suelos someros (Leptosoles réndzicos), altitudes inferiores a 350 m, terrenos planos o con suaves pendientes ubicados en las partes bajas de la sierra. Los encinares de la RBSAT se desarrollan en su mayoría, sobre estratos rocosos de calizas y caliza-lutitas.

Figura 3 Presencia de Quercus oleoides Schltdl. & Cham. en la RBSAT y su área de influencia. 

En la mayoría de los sitios estudiados es común encontrar pocos individuos adultos (3 a 5). Los más longevos se localizaron en el ejido Los Sabinos con alturas promedio de 35 m y DAP de 331 cm, mientras que en Gustavo Garmendia algunos sitios reúnen ejemplares de hasta 25 m de altura y DAP de 220 cm. Los más jóvenes y de menor altura están presentes en los ejidos Laguna del Mante y Gustavo Garmendia donde los rebrotes son numerosos. Se identificó un gradiente ambiental relacionado con las características fisonómicas de los individuos, es decir, los más desarrollados están agrupados en la porción sur de la RBSAT, la cual tiene mayor humedad y precipitación (Figura 3).

De acuerdo con la información recabada en las entrevistas, además de las grandes superficies cubiertas por los encinares que fueron removidos para dar paso a los campos agrícolas y praderas cultivadas entre 1960 y 1970, esta especie fue ampliamente aprovechada por los pobladores de la zona, para proveerse de carbón vegetal y de madera para la construcción de viviendas, lienzos y la elaboración de muebles. Los incendios forestales es otra de las razones que reconocen los pobladores como causas de su desaparición, lo que contra- argumentan algunos autores como ^{Zavala (2000)}, ^{Peña y Bonfil (2003)} y ^{Encina et al. (2007)} quienes sostienen que el fuego es un factor que favorece la presencia de ciertas especies de encinos.

Una vez que se desmontan estos bosques, la presencia de incendios o quemas periódicas y el ramoneo del ganado en el sotobosque hacen que su recuperación sea más difícil (^{Flores, 2007}), ya sea mediante el crecimiento o la regeneración de los tocones, ya que el ritmo de la especie es más bien lento.

Lo anterior pudiera explicar por qué los encinos de interés solo persisten como relictos dispersos en los terrenos poco accesibles y más restringidos del área de influencia de la RBSAT. La mayor parte de las comunidades humanas desconoce los efectos asociados a la deforestación, pero reconoce que es una especie nativa, que se le debe proteger.

Actualmente, la madera de encino, se emplea en labores de construcción de casas, elaboración de muebles en general, como leña, además de la elaboración de carbón y fabricación de utensilios de uso doméstico, lo que evidencia su importancia en las regiones tropicales del país (^{Pérez et al., 2000}; ^{Pennington y Sarukhán, 2005}; ^{Rzedowski, 2006}; ^{Pérez y Dávalos, 2008}; ^{Montoya, 2009}).

Distribución potencial

Entre los diferentes algoritmos de predicción de nichos ecológicos, Maxent es uno de los más utilizados y sus principales aplicaciones permiten la identificación de zonas de protección y conservación de especies amenazadas, la distribución potencial de flora y fauna, los efectos del cambio climático y la cobertura potencial de plagas y enfermedades infecciosas (^{Guisan y Zimmermann, 2000}; ^{Dudík et al., 2004}; ^{Elith et al., 2006}; ^{Phillips et al., 2006}; ^{Phillips y Dudík, 2008}; ^{Austin, 2007}; ^{Mateo et al., 2011}; ^{Cruz et al., 2014}). Por ello, fue elegido para establecer la distribución potencial de Quercus oleoides.

La prueba Jackknife, que determina las variables más importantes para el modelo (^{Phillips et al., 2006}), mostró que la probabilidad de presencia del taxón de interés aumenta en áreas con suelos poco profundos, lo mismo que en terrenos de pendiente suave o planos y cuando la altitud es inferior a los 600 m. Del total de los sitios georreferidos en campo, 76 % se ubican en altitudes entre 300 y 325 m. Asimismo, 89 % de los registros de campo indica que la especie se desarrolla en terrenos planos o casi planos, en pendientes no mayores a 14° donde los suelos suelen ser más profundos por la acumulación de materiales. En ambos casos, coinciden con lo señalado por ^{Rzedowski (2006)} y ^{Pennington y Sarukhán (2005)}.

Para el modelo, la edafología hizo la mayor contribución con 47.3 %, lo que coincide con lo señalado por ^{Márquez et al. (2005)}, ^{Pennington y Sarukhán (2005}) y ^{Rzedowski (2006)}, seguido de la pendiente (17.8 %) y la altitud (12.7 %). El resto de las variables aportaron menos de 10 %, y aunque la temperatura y la precipitación lo hicieron todavía menos, se considera que los climas cálido húmedos, con una precipitación anual de 965 mm y una temperatura media de 24.7 °C, influyen en la presencia de los encinares tropicales y explican su similitud fisionómica con las selvas (^{Pennington y Sarukhán, 2005}).

La capacidad de predicción de los datos en el modelo de distribución mostró un AUC (curva ROC) de 0.988, con una desviación estándar de 0.002, interpretándose como una predicción de prueba óptima (^{Phillips et al., 2006}). Como resultado, Quercus oleoides podría haberse extendido sobre una superficie de 15 528 ha de forma paralela a la sierra, desde los límites de Ciudad Valles hasta la presa La Lajilla (Figura 4). Las áreas más favorables para su eventual recuperación, se ubican al noroeste de la RBSAT, que coincide con lo referido por ^{Sabás et al. (2015)}, quienes reconocen su presencia en la zona.

Figura 4 Distribución potencial de Quercus oleoides Schltdl. & Cham. en la RBSAT y su área de influencia. 

En la distribución potencial por ejido, la mayor superficie corresponde a Laguna del Mante, seguido de Gustavo Garmendia y Los Sabinos. De acuerdo con ^{Zavala (2000)}, los factores antrópicos son determinantes para la presencia de la especie. Las observaciones y registros en campo permiten inferir que la presión humana y las actividades agropecuarias condicionan, en gran medida, la ausencia de esta especie en algunos lugares, a pesar de existir el ambiente propicio para su desarrollo.

Áreas con aptitud para la conservación

La identificación de las áreas con aptitud para la conservación se fundamenta en una serie de etapas que incluye la selección de indicadores cuantificables y variables para discriminar opciones, cuyo propósito es identificar zonas con condiciones y características que deban ser conservadas (^{Hernández, 2014}).

El procedimiento implica seleccionar criterios que puede ser medidos y evaluados para tomar una decisión (^{Eastman, 2012}). Se derivó una composición cartográfica de valores <44 no óptimos y > 195 muy óptimos para la conservación fuera del polígono de protección de la RBSAT. Posteriormente la información se clasificó en cuatro niveles de importancia (no apto, bajo, medio y alto). Se descartaron las áreas poco favorables para la conservación, por corresponder a usos de la tierra agrícola, ganadero o urbano (Figura 5).

Figura 5 Áreas con mayor aptitud para la conservación de Quercus oleoides Schltdl. & Cham. 

Las áreas con valores de 155-216 son zonas con alto potencial para la conservación y se asocian con ejidos apoyados por pago por servicios ambientales, mayor conectividad del paisaje y distancia al área protegida.

El siguiente nivel (115-155) indica una probabilidad media y se asocia principalmente con los aspectos de conectividad y distancia al área protegida, fundamentales para conservar áreas prioritarias (^{Moffett y Sarkar, 2006}; ^{Fabbro et al., 2014}; ^{Hernández, 2014}). Los valores inferiores a 114 son áreas no aptas o poco favorables para la conservación, por su cercanía a zonas con actividades agropecuarias y de crecimiento demográfico (^{Challenger y Caballero, 1998}; ^{Flores, 2007}; ^{FAO, 2016}).

Las áreas con mayor aptitud para la conservación se localizan al suroeste de la Reserva, en los ejidos Laguna del Mante, Gustavo Garmendia, Los Sabinos y Las Palmas con una superficie estimada de 198 ha. Los componentes que contribuyen con mayor peso en la matriz de comparación (AHP) son: el componente ambiental en 37 %, los subsidios y el rezago social, con 34 %.