Introducción

El bosque mesófilo de montaña (BMM) es una comunidad vegetal que se distribuye en las montañas del territorio mexicano, con una estructura, afinidad y composición florística diversa, por lo que su clasificación es ambigua y complicada (Gual-Díaz & González-Medrano, 2014). Este bosque contiene una mezcla de especies de procedencia biogeográfica holártica y neotropical, un componente importante de epífitas (Jardel-Peláez, Excurra, Cuevas-Guzmán, Santiago-Pérez, & Cruz-Cerda, 2004). Su estructura se asemeja a la selva húmeda (Challenger, 1998), aunque en las zonas más secas se han registrado comunidades con fisonomía, estructura y composición florística diferente (Sánchez-Rodríguez, López-Mata, García-Moya, & Cueva-Guzmán, 2003).

Estas comunidades juegan un papel crucial en el mantenimiento de la diversidad biótica (Connell, 1978; Christensen et al., 1989), sin embargo están expuestas a perturbaciones directas e indirectas (Pickett, Kolasa, Armesto, & Collins, 1989) que pueden inducir procesos de cambio continuo a diferentes escalas espaciales y temporales (White & Pickett, 1985; Condit, Hubell, & Foster, 1992; Pickett, Wu, & Cadenasso, 1999). Ejemplos de dichas perturbaciones incluyen la fragmentación del hábitat, la caída de árboles, la incidencia de la ganadería, incendios, vientos fuertes, tormentas relacionadas con huracanes y deslizamientos de tierra (Garwood, Janos, & Brokaw, 1979; Denslow, 1987; Pickett et al., 1989; Phillips & Sheil, 1997; Pickett et al., 1999). Estas alteraciones afectan el ambiente local e influyen en el desarrollo y reproducción de los individuos de las diferentes especies (Williams-Linera, 2007; Santana et al., 2014), además generan cambios espacio-temporales que influyen en la dinámica de poblaciones y especies (Hubbell & Foster, 1990; Condit et al., 1992; Condit, 1995; Chazdon, Redondo-Brenes, & Vilchez-Alvarado, 2005).

No siempre se pueden determinar con exactitud las perturbaciones, porque el desarrollo de las comunidades depende de muchos factores interactuantes que promueven cambios abruptos en la estructura, composición y dinámica de la vegetación (Chazdon, Redondo-Brenes, & Vilchez-Alvarado, 2005), llegando a ser negativos, incluso a escalas espaciales relativamente pequeñas (Homeier, Breckle, Günter, Rollenbeck, & Leuschner, 2010; López & Duque, 2010). El estudio de la estructura del BMM se ha limitado principalmente a trabajos realizados en un tiempo específico (Meave, Soto, Calvo-Irabién, Paz-Hernández, & Valencia Ávalos, 1992; Santiago-Pérez & Jardel-Peláez, 1993; Ruiz-Jiménez, Meave, & Contreras, 2000; Mejía-Domínguez, Meave, & Ruiz-Jiménez, 2004). En los trópicos, la dinámica estructural de bosques, que incluyen diferentes tipos de ensambles, han sido estudiados a través de parcelas permanentes, incluyendo la búsqueda de los factores ambientales y bióticos que influencian los procesos de cambio (Eggeling, 1947; Lugo & Brown, 1981; Lieberman, Lieberman, Hartshorn, & Peralta, 1985; Veillon, 1985; Condit et al., 2005). En la presente investigación, se utiliza información proveniente de ocho parcelas permanentes de 2,500 m² cada una, en un periodo de 12 años, para investigar la dinámica de la estructura del BMM en la Sierra de Manantlán, a través de la pregunta: ¿Existen cambios en la composición florística y atributos estructurales en el tiempo y el espacio?

Método

Área de estudio

La Sierra de Manantlán es un área montañosa localizada al SO de México, en los estados de Jalisco y Colima, entre los 103°51’12’’ y 104°27’05’’ de longitud oeste y 19°26’47’’ y 19°42’05’’ de latitud norte, con una superficie de 139,575 ha, en un gradiente elevacional de 400 a 2 860 m (Jardel-Peláez, 1992). El área específica de estudio corresponde a la comunidad indígena de Cuzalapa localizada en la porción centro-sur de la Sierra de Manantlán (fig. 1). El clima corresponde al cálido Aw2(w)(i´), con temperatura media anual mayor de 22°C y temperatura del mes más frío mayor de 18°C y al semicálido (A)Ca(w2)(w)(i´)g, con temperatura media anual mayor de 18°C y con poca oscilación térmica (García, 2004). El mes más caliente corresponde a mayo y el más húmedo a julio, con porcentaje de lluvia invernal menor de 5%. La precipitación anual oscila entre 1 500 y 1 700 mm. Los suelos corresponden a Entisoles del grupo Lithic Udorthent, con laderas pronunciadas, cuyas pendientes oscilan entre 40 a > 100% (Martínez-Rivera & Sandoval, 1993). El ambiente fisiográfico y antrópico donde se localizaron las ocho parcelas permanentes presentaron una amplitud elevacional de 1 150 a 1 410 m, con pendientes de 18 a 48%, la geoforma dominante fue la ladera media con topografía regular y exposiciones Oeste, Suroeste y Noroeste, información adicional sobre cada parcela se consigna en el cuadro 1.

Establecimiento de las parcelas y registro de datos

Se establecieron ocho parcelas permanentes en áreas con BMM, en una amplitud elevacional de 1 000-1 500 m en 1998. Las parcelas corresponden a cuadros de 2 500 m² (50×50 m), subdivididas en 25 subparcelas de 10×10 m, con la finalidad de facilitar la toma de datos. La localización de cada parcela fue georreferenciada utilizando coordenadas UTM, con un sistema de posicionamiento global marca Garmin, con el datum WGS84 13-N.

Se registró en cada parcela la topografía, geoforma, pendiente, rocosidad, así como algunos factores de perturbación como intensidad de pastoreo, corta de árboles e incidencia de incendios (Olvera-Vargas, Moreno-Gómez, & Figueroa-Rangel, 1996). Todos los individuos arbóreos con diámetro a la altura del pecho (DAP₁₃₀) ≥ 2.5 cm se etiquetaron y se midió su diámetro y altura, esta última fue estimada con un clinómetro Haga. Cada individuo fue ubicado con coordenadas cartesianas. Las parcelas se volvieron a censar en el 2004 y en el 2010. Se recolectó material botánico, el cual fue identificado con el apoyo de especialistas y los ejemplares fueron incorporados al herbario ZEA del Centro Universitario de la Costa Sur de la Universidad de Guadalajara, como material de referencia de la investigación.

Análisis de los datos

Para el análisis de los cambios estructurales de los individuos arbóreos se estimó la densidad, como un indicador de abundancia; el área basal, para estimar la dominancia; y la frecuencia. Se determinó el valor de importancia relativo de cada especie (VIR) mediante la sumatoria de los valores relativos de área basal + densidad + frecuencia (índice de importancia de Cottam) (Mueller-Dombois & Ellenberg, 1974; Matteucci & Colma, 1982). Se estimó, además, en cada uno de los censos, la riqueza de especies, géneros y familias para cada parcela, para conocer los cambios en la composición florística en el espacio y el tiempo.

Debido a que los datos no cumplieron los supuestos de homocedasticidad (Zar, 2010), aplicamos la prueba no-paramétrica de Kruskall-Wallis para conocer si existían diferencias significativas entre parcelas, para las variables estructurales y de composición florística.

Para conocer las diferencias entre el tiempo de estas mismas variables a través de los diferentes censos (C1,C2 y C3), se aplicó una prueba no-paramétrica de Friedman; ésta corresponde a una prueba para muestras relacionadas (Zar, 2010). Ambos procedimientos se efectuaron en el paquete estadístico SPSS v.22.0 (IBM Corp. Released, 2011).

Resultados

Composición florística

Se registraron 115, 128 y 132 especies leñosas, 88, 98 y 100 géneros y 52, 54 y 53 familias para el C1, C2 y C3, respectivamente. Las localidades con mayores cambios en la composición de especies fueron ALP, PMII y ET (cuadro 2). Arroyo La Paloma presenta 10 nuevas incorporaciones de especies en el C2, pero para el C3 se pierden nueve. Patrones similares se registran para PMII y ET (cuadro 2). La riqueza de géneros varió de 88 a 100 del C1 al C3, los géneros con mayor número de especies fueron Quercus (6) y Solanum (5), mientras que Ficus, Cinnamomum, Piper y Prunus presentaron tres cada uno. La riqueza de familias osciló entre 52 y 54 del C1 al C3.

Densidad

La dendidad promedio para los C1, C2, y C3 fue de 1 190, 1 239 y 1 164, respectivamente. El Tigre, BNII y ALP presentaron las densidades mayores en el C1 con 1,524, 1,488 y 1,464 individuos, respectivamente. Las menores las registraron PMII con 896 y LM con 956. En el C2, las densidades mayores fueron para BNII con 1,728 y ALP con 1,520 y la menor para LM con 996. En el C3, BNII tuvo la mayor densidad con 1,680 y PMI la menor con 816. Los cambios en densidad no presentaron un patrón común, pero durante los 12 años estudiados la densidad disminuyó en BNI, PMI y ET; con un incremento en LM y PMII; en BNII, ALP y ASS se presentó un incremento hacia el C2 y una disminución en los dos primero para el C3 (cuadro 3).

Styrax radians fue la especie con mayor densidad en las dos hectáreas, la cual duplica a Conostegia xalapensis, la segunda en importancia. Otras especies con altos valores de densidad fueron Piper villiramulum, Cordia prunifolia, Croton gossypiifolius, Quercus magnoliifolia, Zinowiewia concinna, Schoepfia flexuosa y Oreopanax peltatus (cuadro 4). La dinámica de la densidad entre censos no varió mucho para las especies que presentaron altos valores en este atributo, aunque fue más notorio para las que tuvieron bajas densidades, ya que algunas de ellas murieron de un censo a otro (anexo 1). Las especies que registraron mayores densidades en una o más parcelas fueron Styrax radians, Guarea glabra, Cordia prunifolia, Myriocarpa longipes y Ardisia compressa (anexo 1).

Frecuencia

Styrax radians se presentó en las ocho parcelas durante lo tres censos, mientras que Siparuna tecaphora se encontró en siete de ellas en el C1 y en todas en el C2 y C3. Otras especies con frecuencia entre el 75 y 87% durante los tres censos fueron Cecropia obtusifolia, Conostegia xalapensis, Cordia prunifolia, Inga vera subsp. eriocarpa, Magnolia iltisiana, Piper villiramulum, Symplococarpon purpusii y Xylosma flexuosa. Algunas especies (Toxicodendron sp., Liabum sp., Guazuma ulmifolia, Machaerium salvadorense, Lippia umbellata, Acacia cochliacantha, Aphananthe monoica, Casimiroa watsonii y Alchornia latifolia) sólo se registraron para una parcela y sólo para uno, dos o tres censos.

Área basal

El área basal total para los tres censos fue de 72.3, 72.4 y 75.1 m². El valor más bajo por parcela se registró en el C1 en BNI con 22.33 m² ha^-1 y el mayor en ALP con 62.91 m² ha^-1. En el C2, PMI tuvo el valor menor con 21.06 m² ha^-1 y el mayor para ALP con 65.6 m² ha^-1. En el C3 el valor más bajo se registró en PMI con 16.88 m² ha^-1 y el mayor en ALP con 66.05 m² ha^-1 (cuadro 3). El comportamiento observado en el área basal no permite detectar un patrón común, pues hay parcelas que tuvieron aumento a través del tiempo, otras disminuyeron y luego aumentaron (cuadro 3). Styrax radians, Zinowiewia concinna, Quercus candicans, Myrsine coriacea, Magnolia iltisiana, Lonchocarpus salvadorensis, Clusia salvii aumentaron sus valores de área basal con el tiempo, mientras que Conostegia xalapensis, Podachaenium eminens, Schoepfia flexuosa y Xylosma flexuosa registraron una disminución en este atributo (Anexo 1) (Cuadro 4). Agonandra racemosa, Guarea glabra y Piper villiramulum aumentaron del C1 al C2 y luego disminuyeron en el C3, mientras que Quercus uxoris se redujo del C1 al C2 con una recuperación en el C3.

Valor de importancia relativa (VIR)

Styrax radians fue la especie de mayor VIR en los tres censos, con valores mayores al 10%, seguida por Ficus aurea, Quercus iltisii, Inga vera subsp. eriocarpa, Magnolia iltisiana y Lonchocarpus salvadorensis, lo cual parece estar influenciado por sus valores de área basal. Algunas especies como Conostegia xalapensis, Guarea glabra, Piper villiramulum y Cordia prunifolia presentaron altos VIR influenciados por sus densidades (anexo 1).

En general las especies mantuvieron su VIR a través del tiempo o lo incrementaron ligeramente, con la excepción de las tres especies de Quercus, las cuales redujeron su VIR del C1 al C3. De forma opuesta, Styrax radians se presentó en todas las parcelas con un incremento de su VIR a lo largo del tiempo (anexo 1) (cuadro 4).

Cambios espacio-temporales de la estructura y la composición florística

El BMM de la Sierra de Manantlan presentó diferencias significativas entre parcelas con relación a la densidad (p < 0.0001), mientras que para el área basal no se encontraron. Con relación a la riqueza de especies, géneros y familias, las ocho parcelas presentaron diferencias significativas con p < 0.0001 para las tres variables.

De acuerdo con la prueba no-paramétrica de Friedman se encontraron diferencias significativas (p<0.0001) entre los tres censos para todas las variables estructurales y de la composición florística.

Discusión

Cambios espaciales en la composición florística

Los resultados sugieren al BMM de la sierra de Manantlán como una comunidad espacialmente heterogénea. Alrededor del 70% de los géneros (Clethra, Dendropanax, Meliosma, Oreopanax, Persea, Quercus, Styrax, Symplocos y Zinowewia) corresponden a la flora vascular exclusiva o preferente del BMM (Rzedowski, 1996).

La diferencia entre las parcelas en términos de riqueza de especies, podría atribuirse en primer lugar a la presencia de especies raras, especies turistas (aquellas que se encuentran por periodos breves, no se reproducen y no mantienen poblaciones estables), fenómenos demográficos, al área restringida que ocupa el BMM y su heterogeneidad ambiental (Halffter & Moreno, 2005), así como a factores que cambian o alteran las condiciones biológicas y/o ambientales del hábitat (Halffter & Moreno, 2005; Williams-Linera, 2007; Ruiz-Jiménez, Téllez-Valdés, & Luna-Vega, 2012). La menor riqueza encontrada en BNI, se atribuye a su cercanía con asentamientos humanos, áreas de cultivo de maíz y de agostadero, como ha sido registrado por otros autores (Klug, B., Scharfetter-Lehrl G. y Scharfetter E., 2002; Cascante-Marín, A., V. Meza-Picado y A. Estrada-Chavarría, 2010). La mayor riqueza en las parcelas El Tigre, arroyo La Paloma y Puerto Marínez II, puede ser consecuencia de una topografía más accidentada, laderas escarpadas localizadas en cañadas con mayor humedad, en ecotono con otros tipos de vegetación y con menos incidencia de ganado y actividades agrícolas (Zamudio, Rzedowski, Carranza, & Calderón de Rzedowski, 1992; Cartujano, Zamudio, & Luna-Vega, 2002), o simplemente a la distribución azarosa de las especies (Halffter & Moreno, 2005). La razón de que la riqueza de géneros y familias se mantenga mas o menos constante a través del tiempo se atribuye a que conforme se incrementa el nivel taxonómico, la sensibilidad en los cambios es menor.

Cambios espacio-temporales en la estructura

Cambios en densidad

El comportamiento registrado para la densidad, en las unidades muestreadas, en las cuales la mortalidad (Sánchez-Rodríguez, en proceso), se incrementa después de perturbaciones que alteran el medio físico-biológico de la comunidad, o esta aminora por la ausencia de perturbaciones o la disminución gradual de ellas, lo que permite una recuperación de la densidad por la incorporación de nuevos individuos en las categorías diamétricas inferiores, es un comportamiento ya descrito por otros autores para otras comunidades vegetales (Condit, Hubell, & Foster, 1992; Vílchez-Alvarado, Braulio Chazdon & Milla-Quesada, 2008).

Las diferencias en densidad entre las parcelas y los cambios registrados en el tiempo se relacionan con el comportamiento descrito en el párrafo previo y se atribuye a las condiciones particulares de cada una de las parcelas. Especies heliófilas como Croton gossypiifolius, responden a la apertura de claros, logrando una densidad relativa del 18% con una disminución en el C2 y C3 al 5 y 3%, respectivamente, lo cual podría atribuirse a la mortalidad por denso dependencia. La parcela LM con la menor densidad en 1998, podría ser resultado de un incendio ocurrido antes del muestreo en la primavera de 1998 y la caída de un árbol en 2003, lo que generó claros de dimensiones similares a las reportadas en otras regiones (Durán, Reif, & Hernández, 2011), permitiendo una recuperación en el número de individuos (1224 en el C3). Una perturbación similar se presentó en PMI, lo cual, junto con la influencia de la ganadería observada en el año 2004, podría ser la causa de una disminución de la densidad en un 28% para el C3 y reflejo de las diferencias estadísticas según la prueba de Friedman entre los censos C1 y C2. La parcela ALP, fue la segunda con mayor densidad en el C1, pero para el C3 perdió el 30% de su densidad, pérdidas atribuidas a la ganadería y la corta selectiva de individuos con diámetros menores a 10 cm, factores responsables de las diferencias estadísticas significativas.

Especies con cambios en densidad como Clusia salvinii, Montanoa andersonii, Critonia hebebotrya, Podachaenium eminens, Verbesina fastigiata, Arachnothryx leucophylla, Schoepfia flexuosa, Oreopanax peltatus, Solanum schlechtendalianum, Solanum umbellatum, Myriocarpa longipes, Phenax sp., Piper villiramulum y Conostegia xalapensis respondieron a la apertura de claros (Castillo-Campos & Laborde, 2004; Albuquerque, Velázquez, & Mayorca-Saucedo, 2006; Saldaña-Vázquez, Sosa, Hernández-Montero, & López-Barrera, 2010; Mortero-Montiel, 2011; Morales-Arias, 2016), los cuales se produjeron por la incidencia de incendios forestales, tala clandestina y deslaves provocados por huracanes. Los factores descritos se consideran como un riesgo de pérdida del BMM de la zona de estudio y coinciden con los registrados en otras regiones del país (Toledo-Aceves, Meave, González-Espinosa, & Ramírez-Marcial, 2011; González-Espinosa et al., 2012) y casi idénticas con los BMM de Michoacán (Santana et al., 2014) con las que comparte, además, tala clandestina, actividades productivas agrícolas y ganaderas (Ramírez-Marcial, 2001), así como el cambio de uso del suelo para el establecimiento de huertas de aguacate (Mendoza, Granados, Geneletti, Pérez-Salicrup, & Salinas, 2011).

Aphanante monoica, Casearia arguta, Ficus insipida, Coussapoa purpusii, Calatola laevigata, Beilschmiedia manantlanensis, Alchornia latifolia, Saurauia serrata, Turpinia occidentalis son especies con densidades bajas, lo cual podría atribuirse a que algunas de ellas se encuentran en su límite superior de distribución elevacional como Aphananthe monoica, Ficus insipida y Coussapoa purpusii, las cuales encuentran su óptimo de distribución ecológica en elevaciones menores a 1200 m (Ibarra-Manríquez y Wendt, 1992; Berg, 1989; Cuevas-Guzmán, 2002), mientras que Calatola laevigata, Beilschmiedia manantlanensis, Alchornia latifolia y Turpinia occidentalis son especies de áreas bien conservadas, umbrófilas y con dispersión zoocora (Cuevas-Guzmán, 2002), condiciones que en conjunto podrían estar limitando su densidad. Casearia arguta es una especie que se ve favorecida en espacios abiertos, por lo cual se esperaría que tuviera frecuencias mayores en las localidades estudiadas, entonces lo que podría estar limitando su distribución es la temperatura, ya que se le observa como un elemento común de áreas abiertas en los bosques tropicales subcaducifolios (Cuevas-Guzmán, 2002), en elevaciones de preferencia por debajo de los 400 m (Nee, 1999).

Cambios en frecuencia

Styrax radians, Siparuna thecaphora y Conostegia xalapensis fueron las únicas especies con frecuencia del 100%. Styrax radians debe su importancia en la comunidad a que su abundancia y alta frecuencia se asocia con hábitats perturbados, así como a su alta capacidad de reproducción y establecimiento (Hernández & Carreón, 1987; Fritsch, 1997; Santiago-Pérez, Jardel-Peláez, Cuevas-Guzmán, & Huerta-Martínez, 2009) y por encontrarse en su óptimo elevacional (Carranza, 2008). Comportamiento similar al de S. radians se registra en la Sierra de Cacoma (Morales-Arias, 2016).

La importancia en frecuencia y densidad de Conostegia xalapensis podría atribuirse a su preferencia por espacios abiertos dentro del bosque (Puig, 1976) y su relación favorable con factores antrópicos como la incidencia de fuego y la ganadería (Almeda, 1993; Cartujano et al., 2002). En otros sitios de México, con indicios de perturbación similares al área de estudio, también se registró la presencia de Conostegia (Acosta-Castellanos, 1997; Cartujano et al., 2002).

Otra especie que responde a la dinámica de claros producidos por la caída de árboles y a temperaturas más altas que las prevalecientes en dosel cerrado es Siparuna thecaphora; dichas condiciones son indispensables para la presencia de semillas (Martínez-Garza, Flores-Palacios, De La Peña-Domene, & Howe, 2009), su germinación y establecimiento (Williams-Linera, 1997). Asimismo prefiere lugares húmedos para mantenerse y pasar a etapas de desarrollo más avanzadas (Morales-Arias, 2016), condiciones que se registran en BNII y ET, donde la especie manifiesta sus mayores densidades (Santana et al., 2014).

La importancia de Inga vera respecto a su frecuencia parece responder a su localización en laderas medias de ambientes rocosos y pendientes pronunciadas (Morales-Arias, 2016). Por otra parte, Magnolia iltisiana y Symplococarpon purpusii son parecidas en su importancia tanto en frecuencia como en densidad; ambas se registraron en la parcela LM con un promedio de 20 individuos durante los tres censos. Sin embargo, en el resto de las parcelas en las que coinciden, se encuentran en bajas densidades, lo que sugiere que este comportamiento podría atribuirse a una posible asociación y una dependencia de ambas a algún factor microambiental particular. Este patrón se ha registrado en otros estudios en Minatitlán, Colima (Padilla-Velarde, Cuevas-Guzmán, & Koch, 2008), Manantlán (Muñoz-Mendoza, 1992; Santiago-Pérez & Jardel-Peláez, 1993) y Cacoma (Morales-Arias, 2016).

Piper villiramulum, Cordia prunifolia, Cecropia obtusifolia, Solanum schlechtendalianum y Ardisia compressa deben sus altas frecuencias a que son especies de áreas abiertas, lo cual es coincidente con lo registrado en otras investigaciones (Alvarez-Buylla & Martinez-Ramos, 1992; Martínez-Ramos & Alvarez-Buylla, 1995; Sánchez-Rodriguez et al., 2003; Padilla-Velarde et al., 2008), además parecen requerir de microhábitats húmedos para mantenerse y pasar a etapas más avanzadas (Morales-Arias, 2016). Otras de sus características, como tasas de crecimiento rápido, estrategias eficientes de dispersión de propágulos (Nee, 1993; Hernández-Montero, Rojas-Soto, & Saldaña-Vázquez, 2011), abundante producción de semillas con dominancia en los reservorios, son rasgos típicos de especies pioneras que en parte puede explicar su dominancia estructural en las primeras etapas de la sucesión (Muñiz-Castro, Williams-Linera, & Martínez-Ramos, 2012), lo que las perfila como especies importantes en la dinámica de claros.

Cambios en área basal

El área basal de las parcelas estudiadas, 16.9 a 66 m² ha^-1, coincide con lo encontrado en otros BMM de México (Catalán-Heverástico, López-Mata, & Terrazas, 2003; Santana et al., 2014) y con bosques tropicales de Sudamérica (Poorter et al., 2001; Lewis et al., 2004; Samper & Vallejo, 2007). El área basal, en la mayoría de los sitios estudiados, (como se registró para BNI y BNII, ALP, ASS, PMII y ET) incrementó con la edad del bosque (Lewis et al., 2004; Vílchez-Alvarado et al., 2008). Sin embargo, los cambios entre períodos fueron diferentes; por ejemplo, en ET el cambio en el C2 fue de 3.64 m², mientras que para el C3 la diferencia fue de 1.34 m². Los cambios mayores se atribuyen a la eliminación por competencia de individuos de diámetros pequeños, que al liberar espacio permiten el desarrollo de otros (Condit et al., 1992). LM presentó disminución hacia 2004 y un incremento posterior. Este evento pudo deberse a la muerte y caída de un individuo de Quercus iltisii ocurrida en años previos a 2004. Este individuo era de grandes dimensiones (158.1 cm de diámetro y 28 m de altura) y su derribo pudo ser el responsable de la pérdida de área basal, junto con la eliminación de otras especies durante su caída. A esta situación se agrega la ocurrencia de un incendio superficial previo al primer muestreo de 1998 (Jardel-Peláez et al., 2006).

PMI presentó una pérdida en área basal de 7.79 a 4.2 m² del C1 al C3, lo cual se atribuye a la incidencia de incendios y ganadería, que generan alteraciones, que deben de estar influyendo para las diferencias estadísticas significativas (T² = 10.97 p = 0.0001) encontradas entre censos para la parcela

Las áreas basales mayores se registran para Styrax radians, Inga vera subsp. eriocarpa, Lonchocarpus salvadorensis, Magnolia iltisiana y Quercus iltisii, lo cual se atribuye a sus diámetros de grandes dimensiones y sus densidades altas. La importancia en área basal de Ficus insipida y F. aurea podría deberse a sus características biológicas, como el tener diámetros y alturas grandes (Banack, Horn, & Gawlicka, 2002); además para F. aurea, su sistema multitallos y su comportamiento como hemiepífito (Cuevas-Guzmán, 2002), podrían contribuir al incremento en su área basal.

Cambios en VIR

Styrax radians fue la especie que registró el mayor valor de importancia en todo el BMM del área de estudio, lo cual se atribuye a que es una especie favorecida por la apertura de claros, su mecanismo de dispersión por zoocoría, y su capacidad para rebrotar (Fritsch, 1997). Los valores de importancia mayores registrados en los tres censos para S. radians (10.5 a 11.4%), coinciden con lo encontrado en el BMM de Michoacán donde S. argenteus, un congénere con atributos biológicos y requerimientos ecológicos semejantes, se encontró como la especie más importante (Santana et al., 2014). Las perturbaciones generadas por las plantaciones de aguacate, cambio de uso para agricultura y ganadería, descritas por Santana y colaboradores (2014) y Pérez-Salicrup y colaboradores (por Pérez-Salicrup, Peñaloza-Guerrero, & Aguiar-Eleutério, 2006), y que afectan los bosques mesófilos de montaña de Michoacán, se consideran factores que influyen en la estructura actual, en la cual S. argenteus se presenta como la especie dominante. Un patrón similar, pero con fuentes de perturbación diferentes, entre ellas ganadería, cortas clandestinas e incendios que afectan a los bosques de Jalisco, particularmente en el área de estudio, parecen afectar la estructura de S. radians como ocurre en Michoacán con S. argenteus. La similitud ambiental del área de estudio con algunos bosques mesófilos del estado de Michoacán, como condiciones climáticas cálidas y secas con una estacionalidad de lluvias bien marcada, podrían estar generando las similitudes estructurales en los BMM de ambas localidades, además de su distribución restringida a barrancas o cañadas, donde se conserva una mayor humedad, debido a que las pendientes pronunciadas limitan la radiación solar (Challenger, 1998; Rzedowski, 2006) y la evapotranspiración.

Conclusiones

Cambios espacio-temporales fueron observados, con un incremento en la riqueza de especies, géneros y familias, en las parcelas que tuvieron una recuperación hacia una fase más conservada después de que ocurrió una perturbación.

La tendencia del área basal a través del tiempo, en la mayoría de los sitios estudiados, fue a incrementarse conforme madura el bosque, con una disminución cuando se presenta alguna perturbación.

La densidad de las especies es dinámica, un grupo de ellas como Conostegia xalapensis, Croton gossypifolius, Montanoa andersonii, Piper villiramulum, Podachaenium eminens, Solanum schlechtendalianum, S. umbellatum y Verbesina fastigiata, la incrementan después de que se presenta una perturbación, mientras que otras como Aphananthe monoica, Calatola laevigata, Beilschmiedia manantlanensis, Alchornia latifolia y Turpinia occidentalis responden en forma positiva después de cierto tiempo de que ocurrió la perturbación.

Styrax radians se presentó como la especie más importante a nivel comunidad y una de las más importantes por parcela, cuando se expresó en términos de su valor de importancia relativa, posiblemente por su preferencia por hábitats con cierto grado de perturbación y su alta capacidad de establecimiento.

Los bosques mesófilos del área de estudio muestran ser dinámicos, aún en periodos de tiempos cortos, como el estudiado, respondiendo en cambios en su estructura, diversidad y en importancia de las especies que están presentes, siendo elementos de mucha influencia en la configuración del bosque, los factores de perturbación como la actividad ganadera, la extracción selectiva de especies y la incidencia de incendios forestales. Se considera que mientras se mantengan las perturbaciones a la escala que están ocurriendo, el bosque persistirá y permitirá su dinámica en el largo plazo.

Literatura citada

Acosta-Castellanos, S. (1997). Afinidades fitogeográficas del bosque mesófilo de montaña de la zona de Pluma hidalgo, oaxaca, México. Polibotánica, 6, 25-39.

Albuquerque, L. B., Velázquez, A., y Mayorca-Saucedo, R. (2006). Solanaceae composition, pollination and seed dispersal syndromes in neotropical Mountain Cloud Forests. Interciencia, 31(11), 807-816.

Almeda, F. (1993). Melastomataceae. Fascículo 10. In Flora del Bajío y de regiones adyacentes (pp. 1-36). Instituto de Ecología, A.C. http://www1.ecologia.edu.mx/publicaciones/ resumeness/FLOBA/Flora 78.pdf

Alvarez-Buylla, E., y Martinez-Ramos, M. (1992). Demography and allometry of Cecropia obtusifolia, a Neotropical pioneer tree - an evaluation of the climax-pioneer paradigm for tropical rain forests. Journal of Ecology, 80(2), 275-290.

Banack, S. A., Horn, M. H., y Gawlicka, A. (2002). Disperser- vs . Establishment-Limited Distribution of a Riparian Fig Tree ( Ficus insipida ) in a Costa Rican Tropical Rain Forest 1. Biotropica, 34(2), 232-243. https://doi.org/10.1111/j.1744-7429.2002.tb00534.x

Carranza, E. (2008). Styracaceae. Fascículo 22. Flora de Jalisco y áreas colindantes. Universidad de Guadalajara.

Cartujano, S., Zamudio, S., y Luna-Vega, I. (2002). El bosque mesófilo de montaña en el municipio de Landa de Matamoros, Querétaro, México. Boletín de la Sociedad Botánica de México, 70, 13-43.

Cascante-Marín, A., V. Meza-Picado y A. Estrada-Chavarría. (2011). Tree turnover in a premontane tropical forest (1998-2009) in Costa Rica. Plant Ecol. 212: 1101-1108.

Castillo-Campos, G., y Laborde, J. (2004). La vegetación. In S. Guevara, J. Laborde, & G. Sánchez (Eds.), Los Tuxtlas: el paisaje de la sierra. (pp. 235-261). Xalapa, Veracruz: Instituto de Ecología, A.C. y Unión Europea.

Catalán-Heverástico, C., López-Mata, L., y Terrazas, T. (2003). Estructura, composición florística y diversidad de especies leñosas de un bosque mesófilo de montaña de Guerrero, México. Anales del Instituto de Biología. Serie Botánica, 74(2), 209-230. http://www.redalyc.org/articulo.oa?id=40074204

Challenger, A. (1998). Utilización y conservación de los ecosistemas terrestres de México. Pasado, presente y futuro. México, D. F.: Comisión Nacional para el Conocimiento y Uso de la Biodiversidad, Instituto de Biología, UNAM, Agrupación Sierra Madre.

Chazdon, R. L., Redondo-Brenes, A., y Vilchez-Alvarado, B. (2005). Effects of climate and stand age on annual tree dynamics in tropical second-growth rain forests. Ecology, 86(7), 1808-1815. https://doi.org/10.1890/04-0572

Christensen, N. L., Agee, J. K., Brussard, P. F., Hughes, J., Dennis, H., Minshall, G. W., Williams, S. (1989). the Interpreting Yellowstone fires of ecosystem responses and management implications. Sciences-New York, 39(10), 678-685.

Condit, R. (1995). Research in Large, Long Term Tropical Forest Plots. Tree, 10(1), 18-22.

Condit, R., Ashton, P., Balslev, N., Brokaw, S., Bunyavejchewin, G., y Chuyong Co, L. (2005). Tropical tree a-diversity: Results from a worldwide network of large plots. Biolokiske Skrifter, 55, 565-582.

Condit, R., Hubell, S. P., y Foster, R. B. (1992). Stability and change of a neotropical forest over a decade. BioScience, 42, 822-828.

Connell, J. H. (1978). Diversity in tropical rain forests and coral reefs High diversity of trees and corals is maintained. Science, 199(4335), 1302-1310. https://doi.org/10.2307/1745369

Cuevas-Guzmán, R. (2002). Análisis de gradientes de la vegetación de la cañada El Tecolote, en la Sierra de Manantlán, Jalisco, México. Colegio de Postgraduados.

Denslow, J. S. (1987). Tropical rainforest gaps and tree species diversity. Annual Review of Ecology and Systematics, 18(1), 431-451. https://doi.org/10.1146/annurev.es.18.110187.002243

Durán, C., Reif, A., y Hernández, L. (2011). Understanding small-scale disturbances in the Guyana´s Montane Forests: Gap characterization in Sierra de Lema, Venezuela. Interciencia, 36(4), 272-280.

Eggeling, W. . J. (1947). Observations on the ecology of the Budongo rain forest, Uganda. Journal of Ecology , 34(1), 20-87. https://doi.org/10.2307/2256760

Fritsch, P. W. (1997). A Revision of Styrax (Styracaceae) for Western Texas, Mexico, and Mesoamerica. Annals of the Missouri Botanical Garden, 84(4), 705-761. https://doi.org/10.2307/2992026

García, E. (2004). Modificaciones al sistema de clasificación climática de Köppen. Quinta edición. Serie Libros Núm. 6. México, D.F. Instituto de Geografía, Universidad Nacional Autónoma de México.

Garwood, N. C., Janos, D. P., y Brokaw, N. V. L. (1979). Earthquake-caused landslides: a major disturbance to tropical forests. Science, 205(4410), 997-999. https://doi.org/205/4410/997[pii]\r10.1126/science.205.4410.997

González-Espinosa, M., Meave, J. A., Ramírez-Marcial, N., Toledo-Aceves, T., Lorea-Hernández, F. G., e Ibarra-Manríquez, G. (2012). Los bosques de niebla de México: conservación y restauración de su componente arbóreo. Ecosistemas, 21(1-2), 36-52. https://doi.org/10.7818/re.2014.21-1-2.00

Gual-Díaz, M., y González-Medrano, F. (2014). Los bosques mesófilos de montaña en México. In F. Gual-Díaz, y M. González-Medrano (Ed.), Bosques mesófilos de montaña de México, diversidad, ecología y manejo (pp. 21-68). México, D.F.: CONABIO.

Halffter, G., y Moreno, C. E. (2005). Significado biológico de las Diversidades Alfa, Beta y Gamma. (pp. 5-18). En: Sobre Diversidad Biológica: El significado de las Diversidades alfa, beta y gamma. m3m-Monografias 3ercer Milenio, vol 4. SEA, CONABIO, Grupo DIVERSITAS & CONACYT, Zaragoza. IV+ 242 pp.

Hernández-Montero, J. R., Rojas-Soto, O. R., y Saldaña-Vázquez, R. A. (2011). Consumo y dispersión de semillas de Solanum schlechtendalianum (Solanaceae) por el murciélago frugívoro Sturnira ludovici (Phyllostomidae). Chiroptera Neotropical, 17, 1017-1021.

Hernández, H. M., y Carreón, Y. A. (1987). Notas sobre la ecología reproductiva de árboles en un bosque mesófilo de montaña en Michoacán, México. Boletín de la Sociedad Botánica de México , 47, 5-35.

Homeier, J., Breckle, S., Günter, S., Rollenbeck, R. T., y Leuschner, C. (2010). Tree diversity , forest structure and productivity along altitudinal and topographical gradients in a species-rich Ecuadorian montane rain forest. Biotropica, 42(2), 140-148.

Hubbell, S. P., y Foster, R. B. (1990). Structure, dynamics, and equilibrium status of old-growth forest on Barro Colorado Island. Four Neotropical Rainforests, 522-541.

IBM Corp. Released. (2011). IBM SPSS Statistics for Windows, Version 22.0. 2011.

Jardel-Peláez, E. J. (1992). Estrategia para la conservación de la Reserva de la Biosfera Sierra de Manantlán. Guadalajara, Jalisco: Universidad de Guadalajara.

Jardel-Peláez, E. J., Excurra, E., Cuevas-Guzmán, R., Santiago-Pérez, A. L., y Cruz-Cerda, P. (2004). Vegetación y patrones del paisaje. In R. Cuevas-Guzmán & E. J. Jardel-Peláez (Eds.), Flora y vegetación de la Estación Científica Las Joyas (65-115). Autlán de Navarrro, Jalisco, México: Universidad de Guadalajara.

Jardel-Peláez, E. J., Morfín-Ríos, J. E., Vargas-Jaramillo, S., Michel-Fuentes, J. M., Cuevas-Guzmán, R., Castillo-Navarro, F., Quiñónez, H. (2006). Fire regime and fire effects on Wester Mexico subtropical montane forest ecosystems. In 3rd International Congress on Fire Ecology and Management (pp. 13-17). San Diego, California.

Klug, B., Scharfetter-Lehrl G. y Scharfetter E. (2002). Effects of trampling on vegetation above the timberline in the Eastern Alps, Austria. Arctic and Alpine Research 34(4) 377-388

Lewis, S. L., Phillips, O. L., Baker, T. R., Lloyd, J., Malhi, Y., Almeida, S., Vinceti, B. (2004). Concerted changes in tropical forest structure and dynamics: evidence from 50 South American long-term plots. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 359(1443), 421-36. https://doi.org/10.1098/rstb.2003.1431

Lieberman, D., Lieberman, M., Hartshorn, G. S., y Peralta, R. (1985). Growth rates and age-sizes relationships of a tropical wet forest trees in Costa Rica. Journal TropicalEcology, 1, 97-109.

López, W., y Duque, Á. (2010). Patrones de diversidad alfa en tres fragmentos de bosques montanos en la región norte de los Andes, Colombia. Revista de Biologia Tropical, 58(1), 483-498. https://doi.org/10.15517/rbt.v58i1.5223

Lugo, A., y Brown, S. (1981). Ecological monitoring in the Luquillo Forest Service. Ambio, 10, 541-550.

Martínez-Garza, C., Flores-Palacios, A., De La Peña-Domene, M., y Howe, H. F. (2009). Seed rain in a tropical agricultural landscape. Journal of TropicalEcology , 25(5), 541-550. https://doi.org/10.1017/S0266467409990113

Martínez-Ramos, M., y Alvarez-Buylla, E. R. (1995). Ecología de poblaciones de plantas en selvas húmedas de México. Boletín de la Sociedad Botánica de México , 56, 121-153.

Martínez-Rivera, L. M., y Sandoval L., J. de J. (1993). Levantamiento taxonomico de suelo de la subcuenca de Cuzalapa Sierra de Manantlán, Jal. Terra, 11(1), 3-11.

Matteucci, S. y Colma, A. (1982). Metodología para el estudio de la vegetación. Washington, D.C.: Programa Regional de Desarrollo Cinetífico y Tecnológico de la Secretaría General de los Estados Americanos. http://bases.bireme.br/cgi-bin/wxislind.exe/iah/online/?IsisScript=iah/iah.xis&src=google&base=REPIDISCA&lang=p&nextAction=lnk&exprSearch=137297&indexSearch=ID

Meave, J. A., Soto, M. A., Calvo-Irabién, L. M., Paz-Hernández, H., y Valencia Ávalos, S. (1992). Análisis sinecológico del bosque mesófilo de montaña de Omiltemi, Guerrero. Boletín de la Sociedad Botánica de México , (52), 31-77.

Mejía-Domínguez, N. R., Meave, J. A., y Ruiz-Jiménez, C. A. (2004). Análisis estructural de un bosque mesófilo de montaña en el extremo sur oriental de la Sierra Madre del Sur (Oaxaca), México. Boletín de la Sociedad Botánica de México , 74, 13-29. http://www.socbot.org.mx/publicaciones/boletin/2000/contenido/2004/1/pdf/3.pdf

Mendoza, M. E., Granados, E. L., Geneletti, D., Pérez-Salicrup, D. R., y Salinas, V. (2011). Analysing land cover and land use change processes at watershed level: A multitemporal study in the Lake Cuitzeo Watershed, Mexico (1975-2003). Applied Geography, 31(1), 237-250. https://doi.org/10.1016/j.apgeog.2010.05.010

Morales-Arias, J. G. (2016). Patrones estructurales y heterogeneidad ambiental de un bosque mesófilo de monatña en el Centro Occidente de México. Universidad de Guadalajara.

Mortero-Montiel, O. (2011). Abandono de tierras y el desarrollo de la vegetación secundaria en dos ejidos de La Sierra de Santa Marta, Veracruz. Universidad Veracruzana.

Mueller-Dombois, D., y Ellenberg, H. (1974). Aims and methods of vegetation ecology. Wiley, (August), 93-135. https://doi.org/10.1016/0304-3746(75)90017-7

Muñiz-Castro, M. A., Williams-Linera, G., y Martínez-Ramos, M. (2012). Dispersal mode, shade tolerance, and phytogeographical affinity of tree species during secondary succession in tropical montane cloud forest. Plant Ecology, 213(2), 339-353. https://doi.org/10.1007/s11258-011-9980-5

Muñoz-Mendoza, M. E. (1992). Distribución de especies arbóreas del bosque mesófilo de montaña en la Reserva de la Biosfera Sierra de Manantlán. Universidad de Guadalajara.

Nee, M. (1993). Solanaceae II. Fascículo 72. Flora de Veracruz. Instituto de Ecología, A.C.

Nee, M. (1999). Flacourtiaceae. Fascículo 111. Flora de Veracruz. Instituto de Ecología, A.C.

Olvera-Vargas, M., Moreno-Gómez, S., y Figueroa-Rangel, B. L. (1996). Sitios permanentes para la investigación silvícola. Guadalajara, Jal., México: Universidad de Guadalajara.

Padilla-Velarde, E., Cuevas-Guzmán, R., y Koch, S. (2008). Plantas vasculares y vegetación de la parte alta del Arroyo Agua Fría, municipio de Minatitlán, Colima, México. Acta Botanica Mexicana, 84, 17-23.

Pérez-Salicrup, D. R., Peñaloza-Guerrero, C., y Aguiar-Eleutério, A. A. (2006). Regeneration of Styrax argenteus in natural forest and in plantations of Cupressus lindleyi in Michoacan, Mexico. New Forests, 32(2), 231-241. https://doi.org/10.1007/s11056-005-6104-x

Phillips, O. P., y Sheil, D. (1997). Forest turnover, diversity and CO2. Trends in Ecology and Evolution, 12, 249-250.

Pickett, S. T., Kolasa, J., Armesto, J. J., y Collins, S. L. (1989). The ecological concept of disturbance and its expression at various hierarchical levels. Oikos, 54(2), 129. https://doi.org/10.2307/3565258

Pickett, S. T., Wu, J., y Cadenasso., M. L. (1999). Patch dynamics and the ecology of disturbed ground: a framework for synthesis. Ecosystems of Disturbed Ground. Amsterdam: Elsevier.

Poorter, L. R., Boot, G. A., Hayashida-Oliver, Y., Leigue-Gomez, L., Peña-Claros, M., y Zuidema, P. A. (2001). Tropical humid forest stand dynamics in the Northern Bolivian Amazon. In B. Mostacedo, y T. S. Fredericksen (Ed.), Regeneración y silvicultura de bosques tropicales en Bolivia (pp. 55-76). Santa Cruz, Bolivia: Mostacedo, F. y Fredericksen, T.S.

Puig, H. (1976). Vegetation de la Huasteca, Mexique. Mission Archeologique et Ethnologique Française au Mexique. Etudes Mesoamericaines, 5, 214-222.

Ramírez-Marcial, N. (2001). Diversidad florística del bosque mesófilo en el norte de Chiapas y su relación con México y Centroamérica. Boletín de la Sociedad Botánica de México , 69, 63-76.

Ruiz-Jiménez, C. A., Meave, J. A., y Contreras, J. . L. (2000). El bosque mesófilo de la región de Huautla de Jiménez (Oaxaca), México: análisis estructural. Boletín de la Sociedad Botánica de México , 65, 23-37.

Ruiz-Jiménez, C. A., Téllez-Valdés, O., y Luna-Vega, I. (2012). Clasificación de los bosques mesófilos de montaña de México: afinidades de la flora. Revista Mexicana de Biodiversidad, 83(4), 1110-1144. https://doi.org/10.7550/rmb.29383

Rzedowski, J. (1996). Analsis preliminar de la flora vascular de los bosques mesofilos de montaña de Mexico. Acta Botanica Mexicana , 35, 25-44.

Saldaña-Vázquez, R. A., Sosa, V. J., Hernández-Montero, J. R., y López-Barrera, F. (2010). Abundance responses of frugivorous bats (Stenodermatinae) to coffee cultivation and selective logging practices in mountainous central Veracruz, Mexico. Biodiversity and Conservation, 19(7), 2111-2124. https://doi.org/10.1007/s10531-010-9829-6

Samper, C. K., y Vallejo, M. I. (2007). Estructura y dinámica de poblaciones de plantas de un bosque Andino. Revista Academia Colombiana Científica, 31, 57-68.

Sánchez-Rodríguez, E. (En proceso). Dinámica temporal de un bosque mesófilo demontaña en el Centro Occidente de México. Tesis Doctoral. Universidad de Guadalajara.

Sánchez-Rodríguez, E. V., López-Mata, L., García-Moya, E., y Cueva-Guzmán, R. (2003). Estructura, composición florística y diversidad de especies leñosas de un bosque mesófilo de montaña en la sierra de Manantlán, Jalisco. Boletín de la Sociedad Botánica de México , 73, 17-34.

Santana, G., Mendoza, M., Salinas, V., Pérez-Salicrup, D., Martínez, Y., y Aburto, I. (2014). Análisis preliminar de la diversidad y estructura arbórea-arbustiva del bosque mesófilo en el Sistema Volcánico Transversal de Michoacán, México. Revista Mexicana de Biodiversidad, 85(4), 1104-1116. https://doi.org/10.7550/rmb.41519

Santiago-Pérez, A. L., y Jardel-Peláez, E. J. (1993). Composición y estructura del bosque mesófilo de montaña en la sierra de Manantlán, Jalisco-Colima. Biotam, 5, 13-26.

Santiago-Pérez, A. L., Jardel-Peláez, E. J., Cuevas-Guzmán, R., y Huerta-Martínez, F. M. (2009). Vegetación de bordes en un bosque mesófilo de montaña del Occidente de México. Boletín de la Sociedad Botánica de México , 85, 31-49.

Toledo-Aceves, T., Meave, J. A., González-Espinosa, M., y Ramírez-Marcial, N. (2011). Tropical montane cloud forests: Current threats and opportunities for their conservation and sustainable management in Mexico. Journal of Environmental Management, 92(3), 974-981. https://doi.org/10.1016/j.jenvman.2010.11.007

Veillon, J. P. (1985). El crecimiento de algunos bosques naturales de Venezuela en relación con los patrones del medio ambiente. Revista Forestal Venezolana, 29, 5-120.

Vílchez-Alvarado, B., Chazdon, R. L., y Milla-Quesada, V. M. (2008). Dinámica de la regeneración en cuatro bosques secundarios tropicales de la región Huetar Norte, Costa Rica; su valor para la conservación o uso comercial. Recursos Naturales y Ambiente, 55, 118-128.

White, P. S., y Pickett, S. T. A. (1985). Natural disturbance and patch dynamics: an introduction. (S. T. A. Pickett & P. S. White, Eds.), The ecology of natural disturbance and patch dynamics. Orlando, Florida: Academic Press, Inc. https://doi.org/citeulike-article-id:13497046\rdoi: 10.1016/b978-0-12-554520-4.50006-x

Williams-Linera, G. (1997). Phenology of deciduous and broadleaved-evergreen tree species in a Mexican tropical lower montane forest. Global Ecology and Biogeography Letters, 6(2), 115-127. https://doi.org/10.2307/2997568

Williams-Linera, G. (2007). El bosque de niebla del centro de Veracruz: ecología, historia y destino en tiempos de fragmentación y cambio climático. Xalapa, Veracruz: Comisión Nacional para el Conocimiento y Uso de la Biodiversidad e Instituto de Ecología, A.C.

Zamudio, S., Rzedowski, J., Carranza, E., y Calderón de Rzedowski, G. (1992). La vegetación del estado de Querétaro. Querétaro, Qro., México: Consejo de Ciencia y Tecnología del estado de Querétaro e Instituto de Ecología, A.C.

Zar, J. H. (2010). Biostatistical analysis (5ta ed.). New Jersy: Pearson Prentice Hall.

Notas de autor

^*Autor para correspondencia: valentes@cucsur.udg.mx