Introducción

Los helechos y licófitos son un grupo de plantas vasculares que tienen su mayor representación taxonómica en los bosques templados húmedos de las montañas tropicales (^{Kessler et al., 2011}, ^{Kreft et al., 2010}, ^{Kreft et al., 2004}). Uno de los tipos de bosque de esta zona ecológica en México es el denominado bosque mesófilo de montaña (BMM); se trata de la formación vegetal con mayor riqueza de flora con relación al área que ocupa (aprox. 0.4% del territorio nacional) (^{Gual-Díaz y Gonzáles-Medrano, 2014}, ^{Rzedowski, 1996}, ^{Villaseñor, 2010}). Los helechos y licófitos constituyen el 10% (630 spp.; ^{Tejero-Díez, Torres-Díaz y Gual-Díaz, 2014}) del total de las plantas vasculares consideradas para este tipo de vegetación (6,163 spp.; ^{Villaseñor y Gual-Díaz, 2014}). En una escala más reducida, el BMM en el estado de Veracruz es uno de los 3 más diversos en México, junto con los de Oaxaca y Chiapas. En él habitan al menos 375 especies (50 de manera exclusiva) de helechos y licófitos. Esta cifra representa el 66% de las 564 especies de estas plantas registradas para esta entidad, y el 59% de las consideradas a nivel nacional para este tipo de vegetación (^{Krömer et al., 2015}, ^{Tejero-Díez et al., 2011}).

Debido a la elevada contribución de la pteridoflora en las regiones montañosas del trópico, así como a su sensibilidad a las condiciones ambientales e impacto antrópico, este grupo ha sido referido como indicador para medir aspectos relacionados con la biodiversidad y sus cambios por efecto de la transformación de los ecosistemas (^{Paciencia y Prado, 2005}, ^{Williams-Linera et al., 2005}), incluso en el BMM de Veracruz (^{Carvajal-Hernández y Krömer, 2015}, ^{Carvajal-Hernández et al., 2014}). En el BMM de la región central de Veracruz, ^{Williams-Linera et al. (2005)} observaron que la riqueza de helechos y licófitos no está correlacionada positivamente con la de los árboles, sino que más bien su composición florística fue altamente complementaria. Posteriormente, un estudio en la zona capital de Veracruz constató que las diferencias en la composición de las especies arbóreas influyen en la composición de las especies de helechos y licófitos, aunque en términos de riqueza presentaron valores similares (^{Hernández-Rojas, 2010}).

En el BMM de la región de Los Tuxtlas, en el sureste de Veracruz, se presentan diferentes asociaciones de especies arbóreas dominantes (p. ej., bosques de Podocarpus-Alfaroa), las cuales varían con relación a la orientación de las vertientes de los volcanes San Martín Tuxtla y Santa Marta (^{Conabio, 2010}, ^{Ramírez, 1999}). A la fecha, no existen estudios comparativos de la pteridoflora del BMM de estos 2 volcanes, a pesar de la gran importancia de este hábitat para estas plantas (^{Acebey, Krömer, Vázquez-Torres y Tejero-Díez, 2015}). Se ha mencionado que el BMM es el tipo de vegetación mejor conservado de esta región (^{Castillo-Campos y Laborde, 2004}), sin embargo, está amenazado y en peligro de desaparecer como en otras regiones de México (^{Conabio, 2010}, ^{Williams-Linera, 2002}). El presente estudio pretende contribuir a una mejor comprensión de cómo la heterogeneidad vegetal influye en los patrones de diversidad y composición de los helechos y licófitos en las distintas asociaciones del BMM. Para ello, se plantearon las siguientes preguntas: i) ¿Cómo cambia la riqueza y composición de especies en 3 diferentes áreas de BMM en los 2 volcanes principales de Los Tuxtlas? ii) ¿Cuáles especies contribuyen a estas diferencias y por qué?

Materiales y métodos

El área de estudio se localiza en la sierra de Los Tuxtlas, la cual emerge sobre la llanura costera del golfo de México, ubicada al sureste del estado de Veracruz (^{Conabio, 2010}). La región tiene una superficie de 315,525 ha, de las cuales 155,122 forman parte de la Reserva de la Biosfera Los Tuxtlas (^{Laborde, 2004}). Esta sierra está constituida por 2 grandes macizos volcánicos separados en la parte central por una caldera volcánica que alberga al lago de Catemaco. El macizo noreste está constituido por el volcán San Martín Tuxtla (1,680 m) y el macizo sureste es conocido como la sierra de Santa Marta (1,680 m) (^{Guevara, Laborde y Sánchez-Ríos, 2004}).

El clima de esta región está fuertemente influido por el pronunciado gradiente altitudinal, su topografía compleja y la proximidad al mar del golfo de México. Estos volcanes interceptan los vientos dominantes que provienen del norte, los cuales propician precipitación sobre las cuestas de barlovento y la producción de una sombra de lluvia sobre el lado opuesto (^{Soto, 2004}). Por lo tanto, la precipitación es de 3,000-4,500 (7,000 en las partes altas) mm/año en las vertientes de exposición barlovento, mientras que en las de sotovento es de 1,500-3,500 mm/año (^{Gutiérrez-García y Ricker, 2011}, ^{Soto, 2004}). A pesar de que llueve durante casi todo el año, el régimen lluvioso es marcadamente estacional, con una época húmeda de junio a febrero, y una seca entre marzo y mayo. De acuerdo con la clasificación climática de Köppen modificada por ^{García (1981)}, el clima del BMM, por encima de los 1,000 m de altitud, es semicálido húmedo, con una temperatura media superior a los 18 °C ((A) C (fm)); mientras que, por encima de los 1,500 m, el clima es templado con temperaturas inferiores a los 18 °C (C (fm)). En ambos casos el clima es muy húmedo y en el mes más seco caen al menos 40 mm de lluvia (^{Soto, 2004}). Las frecuentes neblinas características de este tipo de vegetación, y la consiguiente alta humedad atmosférica, suplen la escasez de lluvia en el periodo seco del año (^{Gual-Díaz y González-Medrano, 2014}).

Según ^{Campos (2004)}, a nivel regional, los volcanes San Martín Tuxtla y Santa Marta presentan tipos de suelo distintos; los de la sierra de Santa Marta son geológicamente más antiguos: luvisol (739 km²), acrisol (360 km²) y cambisol (65 km²). En la vertiente norte y noreste, que es más húmeda, la textura de los suelos es franco arcillo limosa y franco limosa (^{Conanp, 2006}). Los suelos del área de influencia del volcán San Martín Tuxtla, de geología más reciente, son predominantemente de tipo andosol (689 km²), la textura de los suelos es de franco arcillo arenosa a franco arcillosa y franca, poseen una capacidad alta para retener agua y tienen altos contenidos de materia orgánica (^{Conanp, 2006}).

Castillo-^{Campos y Laborde (2004)} caracterizaron 9 tipos de vegetación para la región de Los Tuxtlas, donde tanto el BMM como la selva mediana perennifolia son los más conservados, por encontrarse en sitios poco accesibles en las laderas de los volcanes. Sin embargo, el área actual que ocupa el BMM es muy reducida (11,518 ha en la sierra de Santa Marta y 1,072 ha en San Martín Tuxtla; ^{Conabio, 2010}), y se encuentra encima de los 1,000 m de altitud (^{Castillo-Campos y Laborde, 2004}). Se ha observado que entre los (750−) 800 m y 1,100 m de altitud existe una zona de transición entre la selva alta y mediana perennifolia y el BMM (^{Acebey et al., 2015}, ^{Castillo-Campos y Laborde, 2004}), cuya amplitud varía dependiendo de la orientación y del suelo. Para evitar las transiciones, en este estudio se eligieron 3 áreas de BMM por encima de los 1,000 m en los 2 volcanes principales, y en diferente orientación.

La primera área de estudio (SME) se localizó sobre la ladera noreste del volcán San Martín Tuxtla (18°33′34″ N, 95°10′12″ O), abarcando un intervalo altitudinal entre 1,060-1,670 m (Fig. 1). Aquí el BMM se caracteriza por tener un estrato arbóreo de 20-30 m con dosel cerrado y otro con árboles emergentes de cobertura abierta que sobrepasan los 40 m. Los géneros arbóreos más comunes son Carpinus, Ulmus y Quercus. Un tercer estrato que crece por debajo de los árboles de dosel o sotobosque está caracterizado por palmas del género Chamaedorea y helechos arborescentes. Se destaca también la gran biomasa de epífitas vasculares (orquídeas, helechos, bromelias, peperomias y aráceas), y musgos sobre los troncos y ramas más gruesas. En este estudio, de acuerdo con la clasificación de Castillo-^{Campos y Laborde (2004)}, se consideró como parte del BMM la vegetación de la cima del volcán, señalada por otros autores (^{Ramírez, 1999} sensu ^{Miranda y Hernández, 1963}, ^{Sousa, 1968}) como selva baja perennifolia (SBP) o bosque enano (elfin forest), debido a que está dominada por Oreopanax xalapensis y Clusia salvinii, las cuales son comúnmente aceptadas como indicadoras de BMM (^{Rzedowski, 1978}). Esta última asociación del BMM se caracteriza por presentar una altura de dosel de 8-10 m, una gran abundancia de epífitas vasculares y no vasculares, así como palmeras y helechos en el sotobosque a causa de la presencia de neblina persistente por su posición altitudinal encima de 1,500 m.

Figura 1 Mapa de la región de Los Tuxtlas, Veracruz, México, mostrando el área de estudio (rectángulos cuadriculados). SME: ladera noreste San Martín Tuxtla; SMO: ladera suroeste San Martín Tuxtla; SMS: ladera suroeste Santa Marta. 

La segunda área (SMO) fue establecida sobre la ladera suroeste del volcán San Martín Tuxtla (18°33′ N, 95°12′ O) entre 1,060-1,660 m. Por debajo de 1,000 m no se observaron fragmentos de BMM naturales, ya que la vegetación natural ha sido transformada en potreros. Debido a la exposición de sotavento, principalmente en la parte inferior, el BMM se aprecia más seco comparado con el de la primera área, y con una menor abundancia de epífitas no vasculares. Sin embargo, la estructura y composición del bosque es similar.

La tercera área (SMS) fue ubicada sobre la ladera suroeste del volcán de Santa Marta (18°18′19″ N, 94°49′50″ O) entre 1,000-1,570 m. La estructura del BMM de esta área es similar a la de las anteriores, el bosque es poliestratificado denso con árboles de dosel entre 20-25 m y árboles emergentes de hasta 35 m, con 2-3 estratos inferiores. La pteridoflora se destaca por su riqueza y abundancia de especies, principalmente en los estratos herbáceo y arbustivo (^{Ramírez, 1999}). El epifitismo alcanza un alto grado de exuberancia, donde los troncos y las ramas de los árboles están cubiertos por musgos, helechos, bromelias, aráceas y orquídeas. Los géneros arbóreos forman asociaciones que están relacionadas con el intervalo altitudinal, por ejemplo: alrededor de los 1,000 m dominan Liquidambar y Quercus, por encima de 1,200 m los géneros dominantes son Podocarpus, Alfaroa, Quercus, y Liquidambar en las cañadas (^{Ramírez, 1999}). A partir de 1,500 m se presenta, al igual que en el volcán San Martín Tuxtla, una asociación de baja altura (bosque enano), aunque florísticamente la dominancia está dada por especies como Hedyosmum mexicanum, Podocarpus matudae y Weinmannia pinnata.

El muestreo fue realizado entre abril y diciembre de 2005 en el volcán San Martín Tuxtla, y entre agosto y septiembre de 2006 en el volcán de Santa Marta. En el área 1 (SME) se muestrearon en total 13 parcelas, en el área 2 (SMO) 10 parcelas, y en el área 3 (SMS) 12 parcelas. Las parcelas fueron georreferenciadas (coordenadas y elevación) y tuvieron mayormente una forma cuadrada (20 × 20 m); ocasionalmente la forma se adecuó a la topografía de los sitios a muestrear. Este tamaño corresponde al área mínima recomendada por ^{Kessler y Bach (1999)}, que mantiene un criterio de homogeneidad tanto fisionómico como ecológico de las parcelas. En cada parcela se registraron (presencia/ausencia) todas las especies de helechos y licófitos terrestres y epífitos del sotobosque (troncos, arbustos y árboles jóvenes hasta una altura de 8 m en el forófito) (^{Gradstein, Nadkarni, Krömer, Holz y Nöske, 2003}). Para cada especie se indicó la forma de crecimiento de acuerdo con el sustrato según ^{Grayum y Churchill (1987)}: terrestre (T), epífita (E), hemiepífita (HE), rupícola (R). Se recolectaron especímenes por triplicado o tetraplicado y se herborizaron siguiendo las recomendaciones en ^{Lorea y Riba (1990)}. Los ejemplares se clasificaron de acuerdo con ^{Christenhusz, Zhang y Schneider (2011)} y los nombres de las especies se modernizaron de acuerdo con estudios taxonómicos recientes. Se confirmó la adecuada determinación por comparación con especímenes del Herbario Nacional de México (MEXU) y del herbario del Instituto de Ecología, A.C. (XAL), así como con la consulta de especialistas. Un juego completo de ejemplares fue depositado en MEXU y duplicados fueron donados a los herbarios XAL y UC de la Universidad de California.

Para analizar el esfuerzo de muestreo se realizaron curvas de acumulación de especies (^{Colwell, Mao y Chang, 2004}) mediante los estimadores de riqueza no paramétricos ICE y Chao2, los cuales están basados en datos de presencia-ausencia y toman en cuenta las especies raras (únicas y duplicadas). Estos estimadores son considerados más confiables por mantener un crecimiento inicial rápido en la curva y una asíntota bien definida, además de presentar un menor sesgo (^{López-Gómez y Williams-Linera, 2006}, ^{Villaseñor et al., 2005}). Los valores promedio de riqueza observados se ajustaron mediante el modelo de Clench (^{Soberón y Llorente, 1993}) para evaluar la representatividad del muestreo; este método es el más utilizado y ha demostrado un buen ajuste en la mayoría de las situaciones reales y con la mayoría de los taxones (^{Soberón y Llorente, 1993}). Dado que el número de parcelas de muestreo en cada sitio fue distinto, se utilizó el método de rarefacción para compararlas (^{Gotelli y Colwell, 2001}).

La representación gráfica de la similitud florística entre las 3 áreas de estudio se realizó mediante un análisis de escalamiento multidimensional no-paramétrico (NMDS), basado en los registros de incidencia (presencia-ausencia) de especies en cada parcela por sitio de muestreo. La similitud de la composición de especies de helechos y licófitos en las áreas se comparó mediante el análisis no-paramétrico de similitudes (Anosim) (^{Clarke, 1993}), cuya prueba para determinar las diferencias significativas se basó en el índice de Jaccard como una medida de distancia. Este índice es el más utilizado en análisis ecológicos cuando los datos son de incidencia (^{Legendre y Legendre, 1998}). Además, se utilizó la corrección secuencial de Bonferroni para la prueba post-hoc (^{Zar, 1996}). Para identificar cómo cada especie contribuye a la disimilitud se utilizó el análisis no-paramétrico SIMPER, que usa generalmente el valor de Bray-Curtis como medida de similitud (^{Clarke, 1993}). Si bien este tipo de análisis emplea de forma general una variable cuantitativa para determinar la importancia, en este trabajo se utilizó la distancia euclidiana, y lo que se comparó fueron las frecuencias de las especies. Finalmente, se calculó el índice de biodiversidad taxonómica (IB), definido por el número de especies dividido entre el logaritmo natural del área en km² (IB = S/lnA, donde S es el número de especies registradas y A el tamaño del área) (^{Ponce et al., 2002}, ^{Squeo et al., 1998}). Este índice fue empleado tanto para comparar la riqueza de la pteridoflora, entre los fragmentos de BMM estudiados, como con otros en la zona central de Veracruz. Los análisis se realizaron mediante los programas EstimateS Win 9.00 (^{Colwell, 2009}), Statistica (StatSoft, versión 7) y la aplicación Past.exe (^{Hammer, Harper y Ryan, 2001}).

Resultados

Las curvas de acumulación de especies observadas se aproximaron a la asíntota esperada, con excepción de la de SMO (Fig. 2). De acuerdo con el análisis de rarefacción, el traslape de los límites de confianza (95%) en SME y SMS muestra que las diferencias en la riqueza no fueron significativas, mientras que la disparidad del valor para SMO sí lo fue. De acuerdo con los estimadores de riqueza ICE y Chao2, el déficit de muestreo para SME y SMS fue del 8-14%, mientras que para SMO fue del 33-53% (Tabla 1). Tanto para SME y SMS, los valores observados son muy similares a los de los límites de confianza inferiores, a diferencia de SMO. Sin embargo, se obtuvo un buen ajuste de las curvas al modelo de Clench para las 3 áreas (con un coeficiente de determinación R² = 0.99). De acuerdo con este modelo, para que la curva alcance la asíntota, se estima la presencia de 56 especies en SME, 54 en SMS y 48 en SMO, de las cuales se registraron 47, 46 y 38, respectivamente (Tabla 1). Por lo tanto, el déficit de muestreo correspondería al 16% para SME y SMS, y al 21.5% para SMO.

Figura 2 Curvas de rarefacción para las 3 áreas de estudio incluyendo los intervalos de confianza del 95% y valores de riqueza estimados por ICE y Chao2. SME: ladera noreste San Martín Tuxtla; SMO: ladera suroeste San Martín Tuxtla; SMS: ladera suroeste Santa Marta. 

En total se registraron 71 especies de helechos y 3 de licófitos, distribuidas en 37 géneros y 15 familias de acuerdo con la clasificación de ^{Christenhusz et al. (2011)} (Apéndice). En la ladera noreste del volcán San Martín Tuxtla (SME) se encontraron 47 especies, distribuidas en 23 géneros y 11 familias; las familias con 5 o más especies fueron: Polypodiaceae (13 spp., 27%), Dryopteridaceae (9, 19%), Aspleniaceae e Hymenophyllaceae (6, 13% cada una). Los géneros mejor representados fueron Elaphoglossum (8, 17%) y Asplenium (6, 13%). La mayor parte de las especies fueron registradas como epífitas (32, 68%), y terrestres (6, 13%) (Apéndice). En la ladera suroeste del volcán San Martín Tuxtla (SMO) se encontraron 38 especies, distribuidas en 20 géneros y 9 familias. Las familias con 5 o más especies fueron: Polypodiaceae (11, 29%), Dryopteridaceae (8, 21%) y Aspleniaceae (6, 16%). Los géneros mejor representados fueron Elaphoglossum y Asplenium (6, 16% cada uno). Las especies en su mayoría fueron epífitas (22, 58%), y terrestres (8, 21%). En la ladera suroeste de la sierra de Santa Marta (SMS) se encontraron 46 especies, distribuidas en 27 géneros y 12 familias. Las familias con 5 o más especies fueron: Polypodiaceae (21, 45%), Dryopteridaceae (8, 17%) e Hymenophyllaceae (5, 11%). Los géneros mejor representados fueron Elaphoglossum (7, 15%) y Pleopeltis (6, 13%). Las especies en su mayoría fueron epífitas (30, 65%), y terrestres (10, 22%). El IB, es decir, el número de especies con relación al área, osciló entre 4.58 a 5.49 para las 3 áreas de estudio (Tabla 2).

El análisis NMDS de la presencia/ausencia de las especies por parcela en las 3 áreas de estudio (Fig. 3) muestra una marcada separación entre el ensamblaje de especies en el BMM en Santa Marta (SMS) con respecto a las áreas del BMM estudiadas en el volcán San Martín Tuxtla (SME y SMO) (stress: 0.1411). Esta separación es confirmada por el análisis Anosim, donde se indica una diferencia significativa entre las 3 áreas de BMM estudiadas (R = 0.4014 global, p < 0.001). Las comparaciones entre pares de áreas fueron significativamente distintas (p < 0.001), tanto para SME-SMS (R = 0.593) como para SMO-SMS (R = 0.5325); mientras que entre las laderas opuestas del volcán San Martín Tuxtla (SME y SMO) la diferencia no fue significativa (R = 0.073, p = 0.3423). De forma particular para cada área de BMM estudiada, en SME y SMO las parcelas con mayor similitud son esencialmente aquellas dispuestas en las laderas altas entre 1,450 y 1,675 m, mientras que las de los límites altitudinales inferiores son relativamente diferentes. Dicha relación entre altitud y ensamblaje de especies también se refleja en SMS, donde las parcelas entre 1,000 y 1,100 m tienen mayor similitud, en segundo lugar las ubicadas entre 1,200 a 1,400 m y finalmente las que se encuentran entre 1,480 y 1,570 m (Fig. 3).

Figura 3 Ordenación de las parcelas en el NMDS proyectado con el primer y el segundo eje. Los grupos marcados comprenden las parcelas de las 3 áreas: cruces rojas (SME) ladera noreste San Martín Tuxtla; cuadrados azules (SMO) ladera suroeste San Martín Tuxtla y círculos verdes (SMS), ladera suroeste Santa Marta, stress = 0.1411. 

El análisis de similitud (SIMPER) mostró una disimilitud total entre las 3 áreas de BMM del 19.5%. Estos valores están dados por la contribución de especies con frecuencias muy distintas y algunas exclusivas de cada área (Tabla 3). En la comparación de áreas por pares encontramos que SME-SMO obtuvo un promedio de disimilitud del 15.1%, SME-SMS del 22.5%, y SMO-SMS del 20.8% (Tabla 3). Estas diferencias se explican principalmente por la contribución porcentual de especies exclusivas (Pleopeltis polypodioides, Polypodium fraternum, Elapholossum pringlei y Blechnum occidentale registradas solo para SMS), o con frecuencias contrastantes entre áreas (Tabla 3). El área que presentó el mayor número de especies exclusivas fue SMS (22), seguida de SME (10) y SMO (4). Por otro lado, las áreas de BMM estudiadas compartieron entre sí 19 especies (25.7%), de las cuales 8 fueron epífitas, 7 epífitas/terrestres, 3 epífitas/rupícolas y 1 terrestre.

Discusión

Para las áreas SME y SMS se obtuvieron altos porcentajes de completitud (> 84%), por lo cual el muestreo fue satisfactorio. Incluso para SMO, que presentó valores más bajos con relación a los estimadores de riqueza, se registraron el 79% de especies estimadas por el modelo de Clench (Tabla 1). Los menores valores obtenidos por los estimadores no paramétricos se relacionan con el cálculo de la riqueza total pronosticada mediante la estimación de las especies no observadas, que se basa en el número de especies raras únicas (singletons) o duplicadas (doubletons) registradas dentro de las muestras (^{Chao, Chazdon, Colwell y Shen, 2005}). Por lo tanto, en SMO existe una alta incidencia de especies únicas, lo que influye en un bajo porcentaje de especies pronosticadas. Esto podría estar relacionado con un menor tiempo de colonización y con un ambiente más seco en SMO, que promueve la escasez de hábitats favorables para los helechos, sobre todo los higrófilos como grammitidáceas e himenofiláceas (^{Krömer et al., 2014}, ^{Mehltreter, 2010}). En este sentido, solo pocos individuos de ciertas especies con alta especificidad pudieron establecerse, y por lo tanto son raras.

El número de especies de helechos y licófitos registrados entre SME y SMS es muy similar, mientras que en SMO fue significativamente menor. Esto podría deberse al microclima más seco de SMO por el tipo de exposición de la ladera en sombra de lluvia. Este clima más seco se evidencia también por un menor número de especies epífitas, cuya diversidad y abundancia están fuertemente correlacionadas con la precipitación y el aumento de la humedad relativa (^{Kreft et al., 2004}, ^{Kreft et al., 2010}, ^{Krömer et al., 2013a}, ^{Krömer et al., 2005}). El hecho de que los BMM del volcán Santa Marta (SMS), que presentan también una exposición a barlovento, tengan una riqueza mayor a la ladera oeste del San Martín Tuxtla (SMO) podría deberse a su mayor estabilidad histórica relacionada con una mayor edad de sus suelos y bosques. Teóricamente, esta mayor edad permite la colonización de los nichos disponibles por parte de especies provenientes de fuentes cercanas, lo que ocurre de una forma paulatina a través de una escala de tiempo mayor (^{Emerson y Gillespie, 2008}). Se sabe que las erupciones más recientes del volcán Santa Marta datan de hace 2.4 y 1.0 millones de años atrás, y están cubiertas por suelos de hasta 2 m de espesor (^{Martín-del Pozzo, 1997}). En contraste, las últimas erupciones del volcán San Martín Tuxtla ocurrieron en 1664 y 1793 (^{Lozoya, 1984}). Especialmente esta última destruyó la vegetación de los flancos del volcán. En este sentido, a pesar de las condiciones relativamente secas de SMS, las especies han tenido mayor tiempo para establecerse.

Las similitudes en términos de riqueza de especies en SME y SMS, a pesar de las diferencias climáticas y con asociaciones distintas del BMM, coinciden con los resultados obtenidos por ^{Hernández-Rojas (2010)}, quien también obtuvo valores de riqueza iguales o muy similares al comparar bosques diversos y mixtos con bosques dominados por Oreomunnea mexicana y Fagus, respectivamente. Así lo demuestra también el IB para las 3 áreas de este estudio, el cual oscila entre 4.58 y 5.49, muy similar a la variación de este mismo índice (4.92 a 5.79) calculado para algunas localidades de BMM natural en el municipio de Tlalnelhuayocan, región central de Veracruz (Tabla 2), mientras que los valores más bajos (3,04 a 4,19) de dicho índice corresponden a 3 localidades (Parque Ecológico, Xolostla y Banderilla) cercanas a la ciudad capital de Xalapa (^{Williams-Linera, 2002}), donde el BMM está sometido a un mayor impacto antrópico. Estudios recientes han demostrado que la riqueza de helechos y licófitos se pierde entre el 40 y el 70% en la vegetación perturbada o secundaria derivada de BMM (^{Carvajal-Hernández y Krömer, 2015}, ^{Carvajal-Hernández et al., 2014}), por lo que su evaluación es un indicador confiable del estado de conservación del bosque.

El ordenamiento NMDS por parcelas para cada área distingue claramente entre la pteridoflora de SME y SMO del San Martín Tuxtla y SMS del Santa Marta, lo cual es corroborado por el análisis de similitud Anosim (R = 0.4014 global, p < 0.001). La diferencia de similitud florística más obvia es que SMS tiene casi el doble de especies de Polypodiaceae que las otras 2 áreas. Especialmente el género Pleopeltis fue el más numeroso después de Elaphoglossum en SMS, mientras que en SME y SMO fueron Elaphoglossum y Asplenium (Apéndice). En su mayor parte Polypodiaceae se caracteriza por contener especies epífitas de dosel adaptadas a tolerar condiciones de sequía estacional (^{Hietz, 2010}, ^{Kessler, 2010}, ^{Watkins et al., 2006}). Parece ser que estas diferencias en la composición entre los volcanes San Martín Tuxtla y Santa Marta son reflejo de las diferencias climáticas a las que están sometidas las 3 áreas estudiadas. De acuerdo con el estudio de ^{Krömer, Acebey, Kluge et al. (2013)}, la precipitación en SMS y SMO es menor comparada con SME, principalmente a menores altitudes, donde notaron además una menor cobertura de briófitas epífitas, la cual sirve como indicador de humedad atmosférica (^{Karger et al., 2012}). Además, durante el invierno, la vertiente continental suroeste del volcán Santa Marta es más caliente y seca que la vertiente del golfo de México debido a que los «nortes» fríos y húmedos son interceptados por la sierra (^{Soto, 2004}).

Las áreas del volcán San Martín Tuxtla (SME y SMO) comparten el 63% de las especies registradas. Las diferencias más notables en composición entre ambas áreas se dan hacia los límites inferiores del intervalo altitudinal (Fig. 3). Esto podría deberse a que en SMO el clima es más contrastante por un mayor grado de intervención humana a menor altitud, mientras que en SME todavía existe un continuo de vegetación natural hacia las tierras bajas de la Estación de Biología Tropical, (^{Krömer, Acebey, Kluge et al., 2013}). La mayor similitud de composición entre ambas áreas se da a partir de los 1,400 m en las laderas altas del volcán, lo que está relacionado probablemente con su cercanía física, como lo mencionan ^{Ruiz-Jiménez, Téllez-Valdés y Luna-Vega (2012)} en su clasificación de los BMM de México a nivel de floras. No obstante, los helechos y licófitos por su alta capacidad de dispersión difieren en sus patrones de distribución comparados con los de las plantas con semilla, por lo cual no se pueden hacer generalizaciones a nivel regional (^{Sanginés-Franco, Luna-Vega, Alcántara-Ayala y Contreras-Medina, 2011}). Incluso para algunas localidades de BMM no se ha encontrado una mayor similitud florística con relación a la cercanía (^{Ruiz-Jiménez et al., 2012}, ^{Carvajal-Hernández et al., 2014}). Es decir, a pesar de su dispersión efectiva por esporas, factores climáticos e históricos, así como el disturbio antropogénico, tienen un papel preponderante para permitir su establecimiento (^{Karst et al., 2005}, ^{Ruiz-Jiménez et al., 2012}).

En concordancia con las diferencias en composición de especies, el análisis de SIMPER mostró disimilitud en los promedios de abundancia relativa de las especies. Las especies que marcan la disparidad fueron aquellas que reflejan las condiciones climáticas de cada área; por ejemplo, Asplenium auriculatum y Scoliosorus ensiformis en SME podrían tener altos requerimientos de humedad debido a que se encontraron asociadas a lugares umbrosos y húmedos (A.R. Acebey, obs. pers.). En SMO destacan especies terrestres probablemente relacionadas con preferencias específicas del suelo, por ejemplo, Asplenium miradorense y Asplenium achilleifolium. En SMS destacan especies generalistas (Hymenophyllum polyanthos y Melpomene xiphopteroides) y xerotolerantes (P. polypodioides), así como indicadoras del grado de perturbación antropogénica como B. occidentale, que es una especie de lugares abiertos y secundarios (^{Carvajal-Hernández et al., 2014}, ^{Mickel y Smith, 2004}).

Las 3 áreas comparten una sola especie terrestre, Pteris orizabae, mientras que Polystichum mickelii se comparte entre SMO y SMS. Esta disimilitud se explica por la heterogeneidad en la composición de los suelos a nivel local y regional (^{Poulsen et al., 2006}, ^{Tuomisto y Poulsen, 1996}, ^{Tuomisto y Ruokolainen, 1994}, ^{Tuomisto et al., 2002}), dado que las partes altas de los volcanes San Martín Tuxtla y Santa Marta presentan tipos de suelos distintos (^{Campos, 2004}, ^{Flores-Delgadillo et al., 1999}).

Finalmente, el grado de intervención humana afecta además de a la riqueza de especies a la composición de los helechos y licófitos en el BMM (^{Carvajal-Hernández y Krömer, 2015}, ^{Carvajal-Hernández et al., 2014}), ya que las modificaciones en la estructura del bosque causan un microclima más seco y caliente por la mayor entrada de luz (^{Köster et al., 2009}, ^{Paciencia y Prado, 2005}). Por la larga historia de uso de la tierra y mayor densidad poblacional (^{Guevara et al., 2004}), la sierra de Santa Marta, a diferencia del volcán San Martín Tuxtla, presenta una mayor cantidad de especies consideradas secundarias, tales como Pteridium caudatum, Sticherus palmatus, Macrothelypteris torresiana y Gleichenella pectinata. Asimismo, aunque algunas especies higrófilas estaban presentes en Santa Marta, su frecuencia fue más baja, confirmando así la sensibilidad de los helechos y licófitos a las condiciones microclimáticas (^{Krömer et al., 2014}, ^{Larrea y Werner, 2010}).

Se concluye que la riqueza de especies de las asociaciones del BMM en los 2 volcanes fue similar, aunque esta puede variar considerablemente entre áreas de acuerdo con la orientación de la vertiente y sus condiciones ambientales. Estos factores afectan principalmente a la disimilitud en la composición entre ambos volcanes, lo que indica una diversidad beta alta, que enfatiza la importancia de la región de Los Tuxtlas para la conservación de estas plantas y la biodiversidad en general.