text new page (beta)
Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Similars in
SciELO 
Permalink
Introducción
La abundancia y distribución espacial de los seres vivos en la tierra generalmente no es aleatoria ni estable sino que se rige por la amplitud ecológica o capacidad de tolerancia que presenta cada especie a factores ambientales (Chapman, 1976; Tsiftsis, Tsiripidisa, Karagiannakidoua y Alifragisb, 2008; Thuiller, Gassó, Pino y Vilá, 2012). En el caso de las plantas se ha encontrado que su abundancia y distribución se asocia principalmente a factores climáticos, edáficos y geográficos (Kimmins, 2003), así como bióticos, entre los que destaca la interacción competitiva con otras especies y la transformación y pérdida del hábitat que realiza el ser humano debido al cambio de uso de suelo y a las actividades productivas (Meier, Lischke, Schmatz y Zimmermann, 2011) que también generan cambios significativos en las condiciones bióticas y abióticas del ecosistema (Murcia, 1995; Lavergne, Molina y Debussche, 2005).
En ecosistemas templados de México, la abundancia y distribución del género Pinus está influenciada por el clima, el tipo de suelo y la historia biogeográfica de la región (Perry, 1991; Challenger, 1998; Sánchez-González, 2008). Este género es el grupo más importante en el área forestal debido a su historia de aprovechamiento y manejo, sus poblaciones y comunidades ecológicas han variado en una escala espacial y temporal (Richardson et al., 2007). En este contexto, un ecosistema forestal puede ver modificada la cantidad y distribución de especies con potencial productivo debido a la práctica silvícola de seleccionar especies con buenas características morfológicas y de alto valor comercial, por ejemplo, de P. patula a nivel nacional (Castellanos-Bolaños et al., 2008; Ramírez-Santiago, 2015) y de P. halepensis, a nivel internacional (Richardson et al., 2007).
Recientemente, se ha señalado que la abundancia y distribución de los organismos puede cambiar drásticamente en el futuro (Gómez-Mendoza y Arriaga, 2007; van Zonneveld, Koskela, Vinveti y Jarvis, 2009). En México, Gómez-Mendoza y Arriaga (2007) analizaron la sensibilidad de diversas especies de pino al cambio climático bajo dos escenarios, uno severo (22.7 °C y 660 mm) y otro conservador (21.8 °C y 721 mm) para el año 2050, concluyen con drásticas reducciones en el área de distribución actual de este taxa. A una escala regional, Cruz-Cárdenas et al. (2016) mencionan que el cambio climático afecta a especies sensibles como P. leiophylla, P. teocote y P. hartwegii y su área de distribución puede disminuir en más de 70%.
En este contexto, comprender la amplitud ecológica y tolerancia de un organismo a diversos factores ambientales a lo largo de un gradiente, es un aspecto esencial en ecología aplicada, ya que no solo permite describir su nicho ecológico (Dolédec, Chessel y Gimaret-Carpentier, 2000), sino conocer su grado de especificidad, es decir si son generalistas o especialistas, de acuerdo con los intervalos de tolerancia y distribución de los organismos (Hulbert y White, 2007; Thuiller et al., 2012). El uso del análisis multivariado de marginalidad o outlying mean index (OMI) (Dolédec et al., 2000), permite determinar la amplitud, marginalidad y tolerancia del nicho ecológico de las especies (Tsiftsis et al., 2008; Thuiller et al., 2012), lo cual constituye una herramienta importante para predecir potenciales efectos del cambio climático sobre una especie y planificar labores de reforestación en áreas de distribución potencial de estas.
Este estudio se realizó a una escala de paisaje utilizando el análisis de marginalidad y diversas variables ambientales que afectan la abundancia y distribución de pinos en el municipio de Capulálpam de Méndez, Oaxaca, México, comunidad certificada en el manejo forestal de su bosque. Se conocen las especies presentes en su territorio pero se desconocen los factores ambientales que inciden en su abundancia y distribución.
Objetivos
El objetivo del trabajo fue determinar la abundancia y distribución de las especies de pinos y evaluar la sensibilidad de su nicho ecológico con la finalidad de realizar recomendaciones para la conservación y manejo forestal en el municipio de Capulálpam de Méndez.
Materiales y métodos
Descripción del área de estudio
El área de estudio se limitó al municipio de Capulálpam de Méndez, ubicado en la Sierra Juárez al Norte del estado de Oaxaca. Se localiza entre las coordenadas 96° 27' de longitud oeste y 17° 18' latitud Norte. Cuenta con una superficie de 3849 ha (Fig. 1), de las cuales 2557.51 ha se encuentran bajo aprovechamiento forestal (Unión de Comunidades Productoras Forestales Zapotecos-Chinantecos de la Sierra Juárez [Uzachi], 2003). Limita al norte con los terrenos comunales de Ixtlán de Juárez y San Juan Tepanzacoalcos, al sur con terrenos comunales de San Juan Juquila Vijanos, Santiago Xiacui y pueblos mancomunados (Lachatao, Amatlán y Yavesía); al este con la comunidad de San Miguel Yotao y al oeste con terrenos comunales de Santa María Yahuiche (Uzachi, 2003).
Los tipos de suelo que predominan son Cambisol, Litosol y Redzina (Plan de Desarrollo Municipal [PDM], 2009). El clima, según la clasificación de Köppen modificada por E. García, es C(w2) (w) big, definido como un clima templado subhúmedo, el más húmedo de los subhúmedos con lluvias abundantes en verano. La precipitación total es de 1115 mm anuales y los meses con mayor intensidad van de junio hasta octubre, mientras que la temperatura promedio anual es de 15.2 °C (PDM, 2009). La altitud varía de los 2000 m a los 3000 m snm, con una altitud promedio de 2120 m snm (PDM, 2009). Los tipos de vegetación que se encuentran son bosque de pino (P. leiophylla-P. lawsonii), pino-encino (Pinus lawsonii-Quercus urbanii-Quercus crassifolia), latifoliadas, asociaciones de bosques de encino con matorrales secundarios (Quercus-Acacia), y en algunas partes se encuentran pinos con elementos de bosque mesófilo de montaña donde ocurren los géneros Quercus, Oreopanax, Pinus y Dendropanax (Uzachi, 2003; PDM, 2009).
Manejo forestal de la comunidad
En la década de 1950 el bosque de Capulálpam de Méndez fue concesionado 25 años, por decreto presidencial, a la empresa Fábricas de Papel Tuxtepec (Fapatux) para el aprovechamiento de madera del género Pinus, quien utilizó el Método Mexicano de Ordenación de Montes (MMOM) (Uzachi, 2003) también conocido por Método Mexicano de Ordenación de Bosques Irregulares (MMOBI) (Hernández-Díaz, Corral-Rivas, Quiñones-Chávez, Bacon-Sobbe y Vargas-Larreta, 2008). A principios de la década de 1980 la concesión llegó a su fin y las comunidades indígenas recuperaron el control de sus bosques, y a partir de entonces desarrollaron un manejo forestal comunitario (Bray y Merino, 2004). Este tipo de manejo forestal consistió en el Método Desarrollo Silvícola (MDS) con la aplicación de cortas de regeneración de árboles padre, en combinación con el MMOBI, utilizando como método de regeneración el de selección grupal. A partir del año 2003 el método de regeneración de árboles padre fue sustituido por el sistema de corta de regeneración de matarrasa en franjas (Uzachi, 2003).
Obtención de datos de abundancia
Se utilizó información de abundancia proporcionada por la Unión de Comunidades Productoras Forestales Zapotecos-Chinantecos de la Sierra Juárez (Uzachi). Esta consistió en una base de datos con información de abundancia de especies de 636 sitios de muestreo ubicados en áreas de aprovechamiento forestal, obtenidos de censos en parcelas circulares de 1000 m2, que forman parte del inventario forestal del Plan de Manejo Forestal 2015-2022. Adicionalmente, se incorporaron 354 sitios de muestreo de áreas con uso diferente al forestal (agrícola, urbano, turístico), con la finalidad de tener representado todo el municipio. En total se trabajó con una base de datos con 990 sitios con datos de abundancia total y promedio por especie; así como cuatro variables del sitio: coordenadas geográficas (latitud, X y longitud, Y), altitud (ALT) y pendiente (%P). Para analizar las diferencias en la abundancia entre parcelas se realizó un análisis Anova Kruskall-Wallis en rangos y prueba múltiple de contraste de Tukey usando el programa XLSTATTM (AddinSoft, 2014).
Distribución espacial
El análisis espacial fue realizado usando los datos de coordenadas geográficas de cada parcela en el programa PASTTM 2.17 (Hammer, Harper y Ryan, 2001), para analizar la distribución espacial de Ripley's K de cada especie. Cuando la función empírica alcanza valores superiores al modelo nulo a una determinada distancia, implica la existencia de un patrón espacial en agregados significativos a esa distancia; por el contrario, cuando aparece por abajo, implica la existencia de un patrón espacial regular. Adicionalmente, se realizó una prueba de uniformidad de Rayleigh con la que se estimó la orientación cardinal preferente de cada especie, usando el programa PASTTM 2.17 (Hammer et al., 2001).
Obtención del modelo digital de elevación
La ubicación específica de cada especie en Capulálpam de Méndez se situó en un modelo digital de elevación (MDE) del municipio, generado a partir de las curvas de nivel (equidistancia 40 m) del conjunto de datos vectoriales de las cartas topográficas E14D38, escala 1:50 000 (Instituto Nacional de Estadística, Geografía e Informática [Inegi], 2001). Este procedimiento se desarrolló mediante el sistema de información geográfica ArcGis 9.3TM (Esri, 2000) en la extensión 3D Analyst Tools. El MDE se realizó con el sistema de coordenadas UTM, zona 14N, con datum de referencia WGS84 y resolución de 15 m.
Análisis de datos del nicho ecológico
Para realizar el análisis de marginalidad (Outlying Mean Index, OMI) se empleó el software ADE-4TM (Thioulouse, Chessel, Dolédec y Olivier, 1997), donde se utilizaron datos de las variables del sitio y nueve variables climáticas obtenidas de la base de datos WorldClim de la página electrónica http://www.worldclim.org (Hijmans et al., 2005), debido a que el municipio carece de datos climáticos: precipitación máxima (Pmax), precipitación máxima del cuartil más caluroso (Pmts), precipitación mínima (Pmin), precipitación máxima del cuartil más lluvioso (Pmlv), temperatura promedio anual (Tpm), temperatura máxima del cuartil más frío (Tmf), temperatura máxima promedio (Tmaxp), temperatura mínima promedio (Tminp) y desviación estándar de la temperatura (Tdv) y c) los datos de la abundancia de cada especie.
El análisis utiliza dos matrices: una botánica que contiene la abundancia de la especie (en renglones) y los sitios (en columnas); y otra matriz ambiental (climática más geográfica) con valores de las variables ambientales (en renglones) y mismos n sitios (en columnas). Como primer paso, la matriz ambiental fue analizada en sus componentes principales (ACP) con la finalidad de ordenar los sitios en función de las variables ambientales. El OMI organiza la inercia asociada con la distribución de una especie en tres términos: 1) el OMI, que estima la marginalidad de una determinada especie, es decir, la distancia entre las condiciones encontradas en la zona donde se encuentra y las condiciones medias del hábitat en toda el área de estudio. Valores altos de OMI indican que las especies tienen nichos marginales, lo que corresponde a especies que ocurren en hábitats atípicos; 2) la amplitud del nicho o tolerancia de la especie corresponde a la dispersión de los sitios de muestreo que contienen una especie a lo largo de un gradiente ambiental. Valores bajos de tolerancia significa que la especie cuenta con un intervalo limitado de condiciones (e.g. especies especialistas), mientras que los valores altos implican que una especie se distribuye a través de condiciones ambientales muy variables (e.g. especies generalistas); finalmente, 3 la tolerancia residual, que representa la variación en el hábitat no explicada por las variables explicativas del estudio. Valores bajos de este parámetro, indican que las variables analizadas explican la distribución de una especie, mientras que valores altos indican una débil relación (Dolédec et al., 2000; Tsiftsis et al., 2008).
Se utilizó una prueba de Monte Carlo para comparar la marginalidad observada de la especie versus su distribución a través de 10 000 permutaciones al azar y una significancia de P = 0.05, bajo la hipótesis nula de que el ambiente no afecta la abundancia de la especie.
Resultados
Abundancia de las especies
Se encontró un total de 10 especies de pinos, de las cuales una tiene dos variedades, P. pseudostrobus (Tabla 1). Pinus patula presentó la mayor abundancia (4814 individuos) y presencia en los puntos de muestreo (58%), seguida en abundancia de P. pseudostrobus var. pseudostrobus (1253 individuos) y P. douglasiana (1385 individuos), mientras que P. chiapensis y P. leiophylla fueron las especies con menor abundancia y ocurrencia (Tabla 1; Fig. 2).
Tabla 1 Abundancia total (At) y promedio (Ap) de las especies de Pinus en Capulálpam de Méndez.
At-suma de todos los individuos de la especie en los sitios muestreados y Ap-número estimado con base en el promedio por 1000 m2.
*Existe controversia de si P. pseudostrobus Lindl. var. oaxacana (≈P. oaxacana Mirov) es o no una especie. El nombre aceptado en trópicos (www.tropicos.org) una base de datos taxonómica internacional, es Pinus pseudostrobus var. apulcensis (Lindl.) Shaw, no obstante, Conafor (2007a) la reconoce con el nombre brindado por Uzachi.
** A P. rudis autores como Farjon y Styles (1997) y la base de datos trópicos consideran a P. rudis y P. hartwegii como sinónimos. Sin embargo, de la misma forma que en el caso anterior, Conafor (2007b) y la Uzachi la reconocen como P. rudis, por lo que en este estudio se usará el nombre proporcionado por Uzachi, dejando la discusión a futuros trabajos taxonómicos o moleculares, lo cual no es objetivo de este estudio.
En general, las parcelas presentaron una baja diversidad de coníferas. El número promedio de especies presentes por sitio fue de 1.38 ± 0.85. El mayor porcentaje de las parcelas registró, de manera predominante, la presencia de una, dos o tres especies (48%, 30% y 20.94%, respectivamente); lo que representó 98.9% del total analizado. De manera poco frecuente, se ubicaron parcelas con presencia de cuatro a cinco especies (0.91% y 0.15% respectivamente). Este resultado se relaciona con la baja abundancia total y promedio de pinos presentes en las parcelas (Tabla 1) en donde predominaron sitios de muestreo con ausencia de pinos (49% de los sitios analizados) y sitios con menos de 10 individuos en 1000 m2 (27%), siendo escasas las parcelas con más de 50 individuos (< 3% de los sitios analizados).
Distribución de las especies de pino en Capulálpam de Méndez
De manera general, la mayoría de las especies mostró una distribución espacial significativa (R = 1.45; P < 0.016) dispersa a cortas distancias (< de 300 m) y agregada a grandes distancias (de 300 m en adelante), con excepción de P. leiophylla y P. rudis que mostraron una distribución espacial al azar (R = 1.01; P < 0.007).
Utilizando como base el modelo de elevación digital, se definieron tres estratos altitudinales para Capulálpam de Méndez: a) la zona baja que comprende de los 1680 m a 2205 m snm; b) la zona media de 2205 m a 2730 m snm y c) la zona alta con altitudes de 2730 m a 3080 m snm, dentro de los cuales se ubicó la distribución y presencia de cada especie. Para brindar una mejor interpretación de los resultados se describe el patrón de distribución a nivel de cada una de ellas (Fig. 3a-j).
*Se indica a P. pseudostrobus var. pseudostrobus con el símbolo de color rojo, a P. pseudostrobus var. oaxacana de color azul y la presencia de ambas en el mismo sitio se indica de color amarillo.
Figura 3 Distribución de las especies de Pinus en Capulálpam de Méndez.
P. ayacahuite tuvo una mayor presencia en la zona media del municipio (Fig. 3a), con mayor preferencia por la orientación NE (81.4°, R = 0.78, P < 0.005). La especie se registró en zonas húmedas localizadas a una altitud que va de los 2205 m a los 2730 m snm. No se encontró en la zona baja, mientras que su presencia en la zona alta es mínima. Las especies que mostraron una asociación con P. ayacahuite fueron P. patula y P. pseudostrobus var. pseudostrobus presentes en 65% y 23% de los sitios respectivamente y en menor asociación con P. douglasiana (9%), P. leiophylla (0.7%) y P. rudis (0.5%), mientras que en 11% de los sitios se presentó de manera exclusiva.
P. douglasiana mostró un patrón similar a P. ayacahuite, al presentarse con mayor frecuencia en la zona media (2205 m - 2730 m snm) (Fig. 3b). La especie mostró preferencia por laderas y la orientación NE (82.3°, R = 0.77, p < 0.002). Este pino se encontró solo en 69% de las parcelas y cuando se encuentra en asociación es con P. patula.
P. leiophylla se registró solo en 12 sitios que corresponden a laderas de la zona media del municipio (Fig. 3c). Su distribución espacial fue significativa aleatoria y mostró una orientación preferente por el E (93°, R = 0.61, p < 0.009). Su asociación con otros pinos fue baja, menor a 0.5% de las parcelas con P. patula, P. ayacahuite, P. teocote, P. pseudostrobus var. pseudostrobus y P. douglasiana.
P. pseudostrobus var. oaxacana (Pinus oaxacana Mirov) mostró alta presencia y abundancia en la zona baja y media (Fig. 3d) y una orientación NE (78.9°, R = 0.72, p < 0.005). Este pino se asoció principalmente con P. teocote con quién estuvo presente en 45% de los sitios, con P. pseudostrobus var. pseudostrobus en 37% de los sitios y con P. patula en 30% de las parcelas. En menor porcentaje se asoció con P. douglasiana y P. ayacahuite en 3% y 1.7% de los sitios, respectivamente. En 33% de las parcelas muestreadas la especie se presentó de manera exclusiva.
P. patula fue la especie más abundante de todas y su distribución se concentró principalmente en la parte media, en la parte baja se presentó en pequeñas porciones y en la parte alta se le encontró alrededor de los 2900 m (Fig. 3e), ubicándose en las laderas y sitios con una orientación preferente NE (38°, R = 0.53, p < 0.003). P. patula se encontró asociada preferentemente con P. pseudostrobus var. pseudostrobus (42%) y P. ayacahuite (13%).
P. pseudostrobus var. pseudostrobus se registró en las tres zonas (Fig. 3d), con excepción de las altitudes máxima (3080 m snm) y mínima (1680 m snm). Esta especie mostró preferencia por sitios con buena exposición, orientados preferentemente al NE (73°, R = 0.77, p < 0.003). Se encontró asociada en mayor medida con P. patula (64%), P. teocote (25%), P. pseudostrobus var. oaxacana (22%) y P. ayacahuite (10%).
P. rudis, al igual que P. leiophylla, presentó baja abundancia (Fig. 3f). Esta especie se encontró con mayor frecuencia en altitudes elevadas ya que a pesar de encontrarse desde los 2200 m snm con unos pocos individuos, su abundancia se incrementó a partir de los 2500 m - 2800 m snm donde mostró una orientación preferente NE (82°, R = 0.71, p < 0.001). P. rudis se encontró asociada con P. teocote (37.5%), P. patula (29%), P. pseudostrobus var. pseudostrobus (17%), así como con P. ayacahuite y P. douglasiana en 12.5% de los sitios.
P. teocote se registró en solo 18% de los sitios analizados (Fig. 3g); no obstante, es de las especies más abundantes en el municipio, particularmente en la parte baja y media, con presencia menos abundante en la parte alta donde mostró una orientación preferente NE (66°, R = 0.43, p < 0.05). Se le encontró asociada con P. pseudostrobus var. pseudostrobus (45%), P. patula (33%), y en menor medida con P. rudis (8%). Esta especie estuvo sola en 41% de los sitios analizados.
P. lawsonii presentó baja abundancia (Fig. 2) y ocurrencia (Fig. 3h) en el municipio. Esta especie se encontró con mayor frecuencia a altitudes medias, ya que a pesar de encontrarse desde los 2100 m snm con unos pocos individuos, su abundancia máxima se encontró entre los 2250 m y 2400 m snm, donde mostró una orientación preferente NE (85.5°, R = 0.83, p < 0.001). P. lawsonii se encontró asociada con P. patula (60%), P. pseudostrobus var. pseudostrobus (35%), P. teocote (30%) y en menor ocurrencia con P. douglasiana (10%).
P. devoniana solo se registró en 16 parcelas al SE del municipio (Fig. 3i), donde predomina en laderas de la parte baja a una altitud promedio de 1989 m snm, con una distribución espacial significativa aleatoria a cortas y grandes distancias y uniforme a distancias intermedias y una orientación promedio NE (73°, R = 0.88, p < 0.001). Se asoció principalmente con P. pseudostrobus var. pseudostrobus (75%) y en menor medida con P. patula (19%) y P. teocote, (35%).
P. chiapensis mostró presencia solamente en dos sitios del municipio con abundancia muy baja de individuos (Fig. 3j), a una altitud promedio de 2210 m snm, en un patrón agregado y con orientación SO. P. patula fue la única especie asociada en ambos sitios.
Análisis de nicho ecológico
El análisis del índice de marginalidad indicó un efecto directo de las variables ambientales sobre la abundancia y la distribución de las especies de pinos en Capulálpam de Méndez. Ochenta por ciento de las especies son influenciadas significativamente por las variables ambientales analizadas. Las coníferas con mayores valores del OMI fueron P. ayacahuite, P. chiapensis, P. devoniana, P. lawsonii y P. rudis mientras que las cinco especies restantes mostraron un valor menor a la unidad (Tabla 2). Los valores del índice de tolerancia mostraron que los pinos con mayor tolerancia o amplitud de nicho en el municipio fueron P. teocote, P. patula y P. rudis, mientras que las menos tolerantes fueron P. chiapensis, P. lawsonii y P. devoniana (Tabla 2). Sesenta y cuatro por ciento de las coníferas presentaron una tolerancia baja (valores menores a la media), lo que indicó un nicho muy sensible a cambios, es decir, un nicho con baja amplitud o muy especializado. Se encontró que las coníferas mostraron un bajo porcentaje de tolerancia residual, lo cual indica que los parámetros de hábitat analizados explicaron satisfactoriamente su distribución (Tabla 2).
Tabla 2 Parámetros de nicho para las especies de Pinus en Capulálpam de Méndez.
Nombres de especies y abreviaturas se encuentran en Tabla 1. OMI = Índice promedio de marginalidad; TT1 = Índice de tolerancia; T2 = Porcentaje de tolerancia residual; P = representa el número de simulaciones (1000) que representan un mayor valor del observado en el OMI (NS = no significativo a un valor de P > 0.05).
Los primeros dos ejes del análisis OMI explicaron 90% de la varianza obtenida (Fig. 4), en el inciso a de la Figura 4 se muestra que los pinos responden de manera diferente a las variables geográficas y climáticas de la Figura 4 inciso b, mostrando una correlación positiva o negativa significativa. La posición de las especies a lo largo de los dos primeros ejes del OMI permitió separar claramente dos grupos de coníferas; el primero (parte izquierda de la Figura 4a) que consistió en pinos asociados a variables geográficas (latitud, longitud y altitud) y las variables climáticas de precipitación, mientras que el segundo grupo (parte derecha de la Figura 4a) mostró a especies asociadas con las variables climáticas de temperatura.
a) Distribución de las especies en los primeros dos ejes del análisis de marginalidad y b) peso canónico de las variables ambientales. La longitud de las flechas describe la importancia relativa de cada variable analizada y la dirección de la flecha indica la correlación entre variables. Ver Tabla 1 para abreviatura de las especies y sección datos de nicho ecológico en métodos para la abreviatura de las variables.
Figura 4 Valores canónicos de las variables y su efecto sobre la distribución de las especies de Pinus en Capulálpam de Méndez.
Discusión
Las diez especies estudiadas en Capulálpam de Méndez, representan 83% de los pinos documentados para la Sierra Norte de Oaxaca y 71% a nivel estatal (14 especies de acuerdo con Del Castillo, Pérez de la Rosa, Vargas-Amado y Rivera-García, 2004). Este es un porcentaje de especies alto si se compara con regiones de bosques templados de mayor extensión tanto en el estado de Oaxaca (Del Castillo et al., 2004), como a nivel nacional (e.g.Márquez-Linares, González-Elizondo y Álvarez-Zagoya, 1999; Enríquez-Enríquez, Koch y González-Elizondo, 2003; Díaz, Sosa-Ramírez y Pérez-Salicrup, 2012). De las especies presentes, P. patula mostró la mayor abundancia y distribución en el municipio, tres veces más que P. pseudostrobus var. pseudostrobus, que es la segunda especie en abundancia. Este resultado es congruente con su uso al ser P. patula y P. pseudostrobus las especies más aprovechadas forestalmente en México (Perry, 1991; Bermejo y Pontones, 1999; Castellanos-Bolaños et al., 2008).
En Oaxaca, estas dos especies son sumamente apreciadas en el sector forestal, particularmente desde el punto de vista económico, producto de las buenas características morfológicas de sus fustes (Del Castillo et al., 2004; Comisión Nacional Forestal [Conafor], 2007a). En Ixtlán de Juárez, predio comunal colindante con Capulálpam de Méndez, se ha encontrado que P. patula y P. pseudostrobus presentan una mayor abundancia y distribución con respecto a otras especies (Hipólito Vázquez datos no publicados). De igual manera, estas dos especies son registradas por Castellanos-Bolaños et al. (2008) como las de mayor distribución en la Sierra Juárez. Esto puede ser debido a que se impulsa un proceso de selección o reemplazo de especies no comerciales por aquellas con potencial económico en las comunidades forestales de la sierra.
La modificación de la estructura forestal es una consecuencia directa del manejo forestal (Castellanos-Bolaños et al., 2008; Ramírez-Santiago, 2015), lo cual no solo puede favorecer la abundancia y distribución de las coníferas de interés comercial, sino afectar a otras especies (Solís-Moreno et al., 2006; Richardson et al., 2007; Jardel, 2012). La alta presencia de P. patula en el municipio con fines de manejo forestal es ejemplo de una práctica común en el ámbito forestal que conduce a la modificación de la estructura del bosque (Perry, Graham y Richardson, 1998; Richardson et al., 2007). Por consiguiente, una posible respuesta a la baja abundancia y distribución de las especies de menor valor comercial (e.g. P. teocote, P. leiophylla y P. rudis) en la región, sea producto de prácticas inadecuadas de manejo forestal en el pasado cuando Fábricas de Papel Tuxtepec (Fapatux) utilizaba el Método Mexicano de Ordenación de Bosques Irregurales (MMOBI), cuya corta selectiva propició un cambio en la composición de los bosques al reducir la dominancia de pinos y su reemplazo por encinos de acuerdo con Uzachi (2003). No obstante, un análisis histórico de los últimos 50 años realizado por Gorgonio (2015), señala un incremento en las poblaciones de encino en la Sierra Norte a partir de la década de 1980, período en el que termina la concesión a Fapatux e inicia el manejo forestal comunitario. Esta evidencia, en conjunto con la encontrada por Castellanos-Bolaños et al. (2008) y Ramírez-Santiago (2015) exponen el efecto del manejo forestal en la estructura del dosel, sin embargo, se requiere una mayor investigación para concluir qué tipo de manejo forestal (Fapatux o Comunitario) causó la reducción de pinos no comerciales. No obstante, en el futuro se puede causar una reducción de especies de pinos no comerciales, debido a que actualmente los programas de manejo del municipio y la región utilizan como tratamiento silvícola la matarrasa en franja en donde se extraen todas las especies presentes y posteriormente se realizan plantaciones con especies comerciales como P. patula o P. pseudostrobus. Con dichas coníferas se realizan también los programas de reforestación de áreas plagadas o afectadas por incendios, lo que incrementa la abundancia y distribución de estas especies.
Adicionalmente, los procesos históricos de cambio de uso de suelo en el municipio (zonas agrícolas y asentamientos humanos) (Marcial-Juárez, 2011) son la causa de ausencia de pinos en 42% de las parcelas estudiadas, y puede ser otra razón de la baja abundancia de pinos en zonas agrícolas abandonadas, donde se reforestó con especies comerciales. Esta práctica, en conjunto con un manejo forestal inadecuado constituyen factores que restan superficie a los bosques (Challenger, 1998; Sáenz-Romero, Snively y Lindig-Cisneros, 2003; Richardson et al., 2007), que en sinergia con factores ambientales modifican la composición y estructura forestal (Perry, 1991; Sánchez-González, 2008).
Por otra parte, la compleja topografía del territorio nacional, unida a las grandes diferencias a nivel regional en latitud y altitud han dado forma a un mosaico climático y ecológico favorable para la gran diversidad de pinos y distribución de las especies que existe en México (Perry, 1991; Bermejo y Pontones, 1999). Capulálpam de Méndez, cuenta con una accidentada topografía y altitudes que varían desde los 1600 m hasta los 3000 m snm, intervalo altitudinal en el que se sitúa la distribución de los pinos en el país (Bermejo y Pontones, 1999; del Castillo et al., 2004; Sánchez-González, 2008).
Diversos estudios demuestran que tanto las variables geográficas (altitud y latitud) como climáticas (precipitación y temperatura) influyen en la abundancia y distribución de las especies forestales (e.g. Álvarez-Moctezuma, Ochoa-Gaona, de Jong y Soto-Pinto, 1999; Poulos y Camp, 2005; Naranjo-Luna, 2014). Los valores de especificidad del hábitat y amplitud ecológica de las especies de pinos indican de manera general un efecto de las variables ambientales en la abundancia y distribución de las especies de pinos en Capulálpam de Méndez. En algunas especies como el caso de P. patula, los valores pueden estar sobreestimados debido a que su distribución ha sido incrementada debido a las plantaciones en áreas de bajo aprovechamiento, y su marginalidad y amplitud ecológica puede ser más reducida, similar a las demás especies.
La afinidad ecológica e intervalo de tolerancia de los pinos a variables ambientales puede ser relevante al momento de seleccionar a las especies para plantaciones o restauración con especies forestales para establecerlas en su nicho adecuado, particularmente si a futuro se elige un sistema de manejo forestal cercano a lo natural, en donde se predomina la diversidad de especies, y no monocultivos forestales (Bannister y Pyttel, 2010). Adicionalmente, es importante saber cuáles especies son sensibles a factores climáticos para tomar decisiones sobre estrategias de mitigación contra el cambio climático. Un estudio realizado por Gómez-Mendoza y Arriaga (2007) concluye que la distribución a nivel nacional P. rudis, especie presente en Capulálpam de Méndez, es uno de los pinos más vulnerables en México, del que se pronostica una disminución de su distribución de entre 43% y hasta 64%, que coincide con la alta marginalidad y baja amplitud obtenida en este estudio.
La reducción de la distribución de especies forestales descrita a nivel nacional por Gómez-Mendoza y Arriaga (2007); revela un escenario menos radical en la disminución en superficie de otras especies (< 36%), por ejemplo, P. douglasiana, cuya reducción sería de 1% a 2% de su área actual. Sin embargo, este escenario nacional debe tomarse de manera precautoria en el municipio de Capulálpam de Méndez, dada la sensibilidad de las especies a factores climáticos y donde variaciones en temperatura y precipitación a una escala pequeña, puede tener efectos más severos que a nivel nacional. En este contexto Cruz-Cárdenas et al. (2016) encontraron, a nivel regional, decrementos drásticos mayores a 60% en el área de distribución de coníferas sensibles como P. douglasiana, P. rudis, P. leiophylla, P. teocote y P. devoniana.
Las especies anteriores se encuentran también presentes en el área de estudio, donde muestran una vulnerabilidad mediana a alta a procesos de cambio ambiental. En Sierra Juárez, región donde se localiza Capulálpam de Méndez se registran cambios en el régimen pluviométrico e incremento en la temperatura promedio (≈ 0.3 °C) en las últimas tres décadas (Clark-Tapia, Suárez-Mota, Matías-Nolasco y Velasco-Hipólito, 2016). Este cambio ambiental asociado con la alta marginalidad de coníferas incrementa la vulnerabilidad de las especies forestales sensibles, un aspecto que debe considerarse en los planes de manejo forestal del municipio y de la región, de donde se extraen todas las especies del área de aprovechamiento, pero solo se promueve la reforestación con especies comerciales. Lo anterior puede incidir de manera directa en la reducción de las poblaciones y su área de distribución, por lo que se sugiere retener especies de bajo valor comercial, para conservar el nicho de especies asociadas y evitar la pérdida a nivel local de una especie forestal.
Conclusiones
El municipio de Capulálpam de Méndez presenta una alta riqueza de especies de pinos, dentro las cuales Pinus patula presenta la mayor abundancia y distribución, en contraste con P. chiapensis y P. leiophylla que es menor. El manejo forestal, ha generado cambios en la estructura de la vegetación, reflejándose en la dinámica poblacional de las especies forestales en el municipio, lo cual favorece a las especies de alto valor comercial como P. patula y P. pseudostrobus. De manera general, las especies presentaron requerimientos ambientales específicos y alta sensibilidad a procesos de cambio ambiental, por lo que se recomienda la retención de especies de bajo valor comercial en áreas de aprovechamiento forestal, y la utilización de un sistema silvícola cercano a lo natural, donde se promueva la reforestación mixta de especies.
