nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
Introducción
Los humedales son ecosistemas altamente productivos que aportan diversos servicios ecosistémicos (p. ej. provisión y filtración de agua, protección de inundaciones y huracanes, fertilización y aporte de sedimentos, captura de carbono) (M.E.A., 2005; Costanza et al., 2008). Sin embargo, a nivel mundial son uno de los más amenazados debido a actividades humanas como la agricultura, la ganadería (Mitsch et al., 2009) y más recientemente la urbanización. México ha perdido 62% de sus humedales y específicamente los de la costa centro oeste del Golfo de México (estado de Veracruz) son de los más amenazados debido a la fuerte presión ganadera, ya que esta entidad ocupa el primer lugar a nivel nacional en producción de carne de bovino (Contreras-Espinosa y Warner, 2004; Landgrave y Moreno-Casasola, 2012; Moreno-Casasola et al., 2014).
Diversos estudios han asociado el pastoreo del ganado bovino con la pérdida de biodiversidad (Verdú et al., 2000; Watkinson y Ormerod, 2001; Durant et al., 2008) y con la disminución de la riqueza de especies nativas. Por ejemplo Dorrough et al. (2006) encontraron que en bosques con pastizales en Australia, la riqueza de plantas nativas está correlacionada negativamente con la intensidad agrícola y de pastoreo, mientras que varias especies exóticas se veían favorecidas por estas actividades. Sin embargo, existen estudios que mencionan que la riqueza de especies vegetales en ciertos ecosistemas puede verse beneficiada con el pastoreo (Dupré y Diekmann, 2001; Pykälä, 2005).
En algunos humedales del estado de California (EUA) y de Brasil se ha observado que zonas sin pastoreo presentan mayor cobertura de pastos exóticos y menor cobertura de especies nativas que áreas con pastoreo, y se ha hecho énfasis en que la ganadería no es necesariamente dañina para los humedales, siempre y cuando el pastoreo sea de baja intensidad (Marty, 2005; Junk y da Cunha, 2012). En México, el único trabajo que existe al respecto es el de Travieso-Bello et al. (2005), quienes reportaron que la ganadería puede disminuir la riqueza de especies nativas de humedales, y también la biomasa; sin embargo, en este estudio no se realizaron tratamientos excluyendo al ganado y solo se muestreó una época del año, por lo que no presentaron comparaciones temporales.
Es necesario empezar a generar conocimiento acerca del nivel de impacto que tiene la ganadería en diferentes tipos de humedales (manglares, selvas inundables, popales, tulares), principalmente en regiones con tendencia al aumento de esta actividad o en donde la ganadería ocupa grandes extensiones de estos ecosistemas. Esta información podría ayudar a mejorar las políticas de manejo ganadero, logrando una ganadería sustentable en la zona costera, y garantizando la conservación de las funciones y servicios ecosistémicos. En este trabajo se evaluó el efecto de la ganadería y la estacionalidad sobre la composición florística, cobertura vegetal, biomasa aérea (BA) y biomasa subterránea (BS) de los humedales herbáceos (popales) de la costa centro oeste del Golfo de México, específicamente los de la laguna de Alvarado, Veracruz, para responder las siguientes preguntas: 1) ¿La ganadería tiene un impacto negativo en estos humedales? y 2) ¿La variación estacional determina el nivel de impacto ganadero y la composición de la vegetación? La hipótesis de trabajo plantea que una densidad baja de ganado mantendrá la biodiversidad florística del humedal, presentándose cambios estacionales en la composición de las especies.
Materiales y métodos
Área de estudio
El área de estudio se localiza cerca de Alvarado, Veracruz, en la costa centro oeste del Golfo de México, entre los 18°46'N, 95°34'E y los 18°42'N, 95°58'E (Fig. 1). La altitud es de 10 m y el clima es cálido con abundantes lluvias en verano. La temperatura es de 22-26 ºC y la media de la precipitación pluvial es de 1748.3 mm (INEGI, 2005). Los datos fueron obtenidos para el proyecto Nº 48247 CONACYT-CONAGUA “Inventario, delimitación, caracterización y uso sustentable de los humedales de la cuenca del Río Papaloapan, México”. Otros datos derivaron de entrevistas con los propietarios de los terrenos, observaciones de campo y mediciones.
Figura 1: Localidades de Río Blanco (RB1, RB2) y Río Limón (RL1, RL2) en los humedales de Alvarado, Veracruz, México.
En los sitios de muestreo la vegetación natural que domina es el popal, término asignado por Faustino Miranda (Miranda y Hernández-X, 1963). Este ecosistema corresponde a una comunidad vegetal sobre superficies pantanosas permanentemente estancadas en la planicie costera de Tabasco, sur de Veracruz, Chiapas y Campeche. Se establece sobre llanuras aluviales, prácticamente sin declives, atravesadas por ríos que por medio de filtraciones e inundaciones cubren la superficie en zonas con climas tropicales (Rzedowski, 2006). La vegetación se caracteriza por presentar herbáceas emergentes enraizadas de 1 a 3 m de alto con hojas anchas. Algunas especies que se han reportado para los popalestulares del estado de Veracruz son Typha domingensis Pers., Thalia geniculata L., Pontederia sagittata C. Presl, Sagittaria lancifolia L., Sagittaria latifolia Willd., Cyperus articulatus L., Cyperus giganteus Vahl, Eleocharis cellulosa Torr., Hymenachne amplexicaulis (Rudge) Nees, Heliconia latispatha Benth., Leersia ligularis Trin., Hydrocotyle bonariensis Lam. y Bacopa monnieri (L.) Wettst. (Moreno-Casasola et al., 2010).
La economía del municipio de Alvarado está basada fundamentalmente en la ganadería bovina, en donde las razas más empleadas para esta actividad son el ganado Cebú para aprovechamiento de carne y el ganado Suizo y Holandés para aprovechamiento de leche (Portilla Ochoa, 2005). Es común que el pastoreo del ganado sea extensivo y también rotatorio, ya que se saca a los animales del terreno en el período de mayor inundación llevándolos a zonas más altas. La carga animal permitida en el municipio es de una cabeza de ganado por hectárea; sin embargo, en muchos ejidos esto no se respeta, excediendo la carga permitida (INEGI, 2007).
Para este estudio se seleccionaron cuatro sitios, todos impactados por la ganadería, ya que después de una búsqueda extensiva y conversaciones informales con los habitantes de la localidad y la Asociación Ganadera, no se encontraron humedales herbáceos totalmente libres de ganado para considerarlos como sitios control. Se seleccionaron dos localidades de trabajo, una fue Río Blanco, con dos sitios de muestreo (a partir de aquí llamados RB1 y RB2). La segunda localidad fue Río Limón, también con dos sitios de muestreo (RL1 y RL2). En RB1 el ganado es introducido aproximadamente en febrero o marzo y lo sacan aproximadamente en septiembre. En RB2 el ganado entra en enero y lo sacan en septiembre. En RL1 y RL2 el ganado es introducido entre enero y febrero y lo sacan entre octubre y noviembre. Ver Cuadro 1 para una descripción más detallada. Esta variación en el momento de extracción del ganado y el tiempo que permanece fuera del humedal depende del nivel de inundación de los predios particulares.
Cuadro 1: Características de los diferentes tipos de manejo en los cuatro sitios de muestreo localizados en el municipio de Alvarado, Veracruz, México.
| RB1 | RB2 | RL1 | RL2 | |
|---|---|---|---|---|
| Coordenadas | 18°44'49.66"N | 18°45'23.51"N | 18°37'25.49"N | 18°37'39.39"N |
| 95°52'42.24"O | 93°53'36.38"O | 95°53'24.95"O | 95°52'30.87"O | |
| Manejo | 1 cabeza ha-1 | 2 o 3 cabezas ha-1 | 2 o 3 cabezas ha-1 | 2 cabezas ha-1 |
| Influencia humana | Muy poca, en los cinco muestreos nunca se observó a personas trabajando en el terreno. | Muy poca, en los cinco muestreos solo se observó dos veces al capataz cuidando el terreno. Solo iba de entrada por salida. | Se encuentra la casa del dueño a 120 metros en línea recta al oeste del sitio. El dueño va entre una y dos veces por día al sitio en caballo para arrear a las vacas. | En el terreno se encuentra la casa de los dueños, el señor al menos una vez al día arrea al ganado en caballo. |
| Hidrología | Sin modificar | Sin modificar | Sin modificar | Sin modificar |
| Humedad del suelo e inundación | Permanece inundado 9 meses/año, y el resto del año el suelo se encuentra saturado en agua. | Permanece inundado 8 meses/año; el suelo está saturado en agua 2 meses y seco el resto del año. | Permanece inundado 6 meses y seco el resto de año. | Permanece inundado 9 meses/año, y el resto del año el suelo se encuentra saturado en agua. |
| Permanencia de ganado en el sitio | 6-10 meses | 8-10 meses | 10 meses | 10 meses |
Diseño experimental
Debido a que no se encontró ningún popal no impactado por la ganadería, no se tuvo un sitio control. De tal manera no se pudo medir la verdadera magnitud del impacto, y éste solo fue estimado de manera parcial, en relación con el sitio que presentó más características de un humedal herbáceo conservado.
RB1 y RB2 se encuentran separados entre sí por 2.5 km, mientras que RL1 y RL2 por 2 km. Los sitios de Río Blanco se encuentran a 14 km de los de Río Limón (Fig. 1). En cada sitio se colocaron ocho cuadros al azar de 1 m2 cada uno. Para evitar la entrada del ganado, cuatro de estos cuadros fueron cercados con alambre de púas de la siguiente manera: dos cuadros se colocaron dentro de un cerco de alambre de púas, cuidando que haya un espacio entre cuadros y el alambre exterior de al menos 1 m, y en otra área del lugar se colocaron los otros dos cuadros de la misma manera. Se realizaron tres muestreos en temporada de secas (enero, marzo y mayo de 2009) y dos en la de inundación (septiembre y noviembre de 2009).
Estructura y composición florística
En cada cuadro se tomaron datos de composición de especies vegetales, porcentaje de cobertura vegetal para cada especie por medio del método de cobertura-abundancia de Westhoff y Van der Maarel (1978). Se estimó la frecuencia y cobertura relativa, y con esto se obtuvo el Valor de Importancia Relativa (VIR) (Moreno-Casasola y López Rosas, 2009) para cada una de las especies. Para obtener los valores de diversidad y de equidad por sitio, temporalidad y tratamiento, se utilizó el índice de Shannon-Wiener y el índice de equidad de Pielou (Liang et al., 2007). El reemplazo de especies únicamente se realizó por sitio y se calculó por medio del índice de Whittaker (1972) para RB2, RL1 y RL2, en relación con el sitio RB1, ya que éste es el que recibe menos influencia humana y tiene menor carga animal.
La biomasa áerea (BA) se obtuvo cortando la vegetación al ras del suelo dentro de cada cuadro de muestreo. Se colectó toda la materia incluyendo la que estaba enraizada a ese nivel y la fitomasa seca y muerta. Esto fue al inicio (enero 2009) y al final del experimento (noviembre 2009), es decir en noviembre se colectó la biomasa que creció durante once meses en presencia o ausencia del ganado. Después de cada corte, el material se secó durante 120 horas a 65 °C y posteriormente se pesó. La biomasa subterránea (BS) se obtuvo por medio de un núcleo de 15 cm de largo por 10 cm de diámetro. En cada esquina (por fuera) de cada cuadro de 1 m2 se tomaron cuatro muestras, durante los cinco muestreos. Estas muestras se lavaron en laboratorio hasta eliminar la mayor cantidad de tierra posible. Posteriormente se pusieron en flotación para obtener la mayor cantidad de material vivo (Cejudo y Capistrán, 2009). Para separar este material se utilizaron dos tamices, uno de 0.1 mm y otro de 0.039 mm de abertura para raíces finas. Una vez separadas las raíces, se secaron en un horno a 70 °C por 24 horas, para obtener el peso seco.
Análisis estadísticos
Para conocer si el índice de diversidad, la riqueza media de especies, la BA, la BS y el VIR por especie varían entre los sitios de estudio, se aplicó un ANOVA de una vía y la prueba de rangos múltiples de Tukey. Cuando las variables no pudieron ser normalizadas se utilizó una ANOVA por rangos de una vía de Kruskal-Wallis y la prueba de comparaciones múltiples de Dunn. La riqueza media por temporadas (secas-lluvias) se calculó tomando en cuenta los valores de cada uno de los ocho cuadros en cada sitio. La riqueza media, BA y BS entre temporadas (lluvias-secas) y tratamientos (presencia de ganado y exclusión de ganado) se compararon mediante una prueba de t pareada o si los datos no pudieron ser normalizados se utilizó una prueba de Rango de Wilcoxon (Zar, 1999). No se realizaron comparaciones estadísticas en los casos en donde la especie solo se presentaba en un solo tratamiento o una sola temporada, solo se menciona el porcentaje de VIR en donde la especie estuvo presente.
Para visualizar la similitud entre los sitios de muestreo y cuadros con y sin ganado en función de las variables medidas se realizaron análisis de componentes principales (ACP). Es importante mencionar que en los ACPs no se incluyeron todas las especies. Las especies que quedaron fuera del análisis son todas ruderales, a excepción de Lippia nodiflora (L.) Michx., y solo se presentaron en un sitio, en un muestreo (generalmente fue el de septiembre) y en un solo cuadro (de alguno de los ocho del sitio), por lo que fueron especies que no aportaron información relevante al análisis. El nivel de significancia utilizado en todos los análisis fue de 0.05 y se utilizaron los paquetes estadísticos R versión 2.5.1 (Crawley, 2014), Statgraphics plus 5.1 y MVSP 3.1. (Manugistics, 1998).
Resultados
Variación por sitio
Se registró un total de 29 especies (Anexo) y con base en los resultados obtenidos se encontró que había especies clave con VIRs arriba de 40 (Cuadros 2). Con respecto a la diversidad, riqueza, BA y BS, se observó un gradiente de perturbación entre los cuatro sitios, siendo RB1 el menos impactado por la ganadería, seguido de RB2, RL1 y RL2 (Cuadro 3).
Cuadro 2 Se muestran únicamente las especies que presentaron los valores más altos de importancia relativa (VIR) para cada sitio de muestreo en la costa centro oeste del Golfo de México. El valor indica la media ± 1 EE del VIR de los cinco muestreos del año. Se aplicó una ANOVA y la prueba de rangos múltiples de Tukey (*) solo en los casos en donde las especies fueron encontradas en más de dos sitios. Solo se encontraron diferencias significativas en Sagittaria lancifolia L. (F=17.04, p=0.00) que se encontró en los cuatro sitios, y en Thalia geniculata L. (F=11.88, p=0.008) presente en dos sitios. Valores con la misma letra son estadísticamente iguales con p>0.05.
| RB1 | RB2 | ||
| Sagittaria lancifolia L. | 77.3±2.79*a | Echinochloa colona (L.) Link | 86.3±4.32 |
| Eleocharis cellulosa Torr. | 69.3±7.47 | Eleocharis celulosa Torr. | 42±12.04 |
| Hydrocotyle verticillata Thunb. | 68.5±2.02 | Bacopa monnieri (L.) Wettst. | 40.9±11.15 |
| Echinochloa colona (L.) Link | 68.5±17.57 | Sagittaria lancifolia L. | 21.7±6.3*b |
| Bacopa monnieri (L.) Wettst. | 56.1±5.84 | Lemna minor L. | 21.1±13.51 |
| Sporobolus virginicus (L.) Kunth | 25.7±16.77 | Thalia geniculata L. | 18.3±7.72*a |
| Nymphaea ampla (Salisb.) DC. | 16.5±10.13 | ||
| RL1 | RL2 | ||
| Eleocharis cellulosa Torr. | 67.8±8.16 | Hymenachne amplexicaulis (Rudge) Ness | 73.1±3.43 |
| Hymenachne amplexicaulis (Rudge) Ness | 66.1±17.38 | Thalia geniculata L. | 48.4±4.09*b |
| Leersia hexandra Sw. | 65.9±7.92 | Ludwigia octovalvis (Jaq.) P.H. Raven | 54.2±2.34 |
| Sagittaria lancifolia L. | 33.3±8.01*b | Eleocharis cellulosa Torr. | 44.7±5.54 |
| Cyperus articulatus L. | 27.6±11.95 | Pontederia sagittata C. Presl | 36.3±1.36 |
| Nymphaea ampla (Salisb.) DC. | 3.47±3.47 | Ambrosia psilostachya DC. | 24.4±15.44 |
| Sagittaria lancifolia L. | 24.0±6.84*b | ||
| Paspalum sp. | 26.3±16.17 |
Cuadro 3: Diversidad de especies (índice de Shannon - H’ e índice de equidad de Pielou - J’), riqueza de especies, biomasa aérea (BA) y biomasa subterránea (BS) para los cuatros sitios de muestreo en el municipio de Alvarado, Veracruz, México. En el caso del índice de Shannon se indica el valor de éste (H’) y el tamaño de muestra (S). Se indica la media ± 1 EE en los casos donde se aplicó un ANOVA y la prueba de rangos múltiples de Tukey (*) y se muestra la mediana y los percentiles (25-75%) en los casos donde se aplicó un ANOVA por rangos de Kruskal-Wallis y una prueba de Dunn (**). Valores con la misma letra son estadísticamente iguales con p>0.05.
| Característica | RB1 | RB2 | RL1 | RL2 | Valor estadístico | Valor de p |
|---|---|---|---|---|---|---|
| Índice de diversidad | S=9, H’=1.81 b | S=10, H’=1.63 b | S=9, H’=1.6 b | S=19, H’=2.32 a | H=6.28 | 0.04 |
| Índice de equidad | J’=0.95a | J’=0.61bc | J’=0.85a | J’=0.77ac | F=6.72 | 0.00 |
| Riqueza de especies* | 5±0.31 b | 5.6±1.36 b | 5.2±0.37 b | 9.2±1.28 a | F=4.71 | 0.01 |
| BA (g m-²)* | 660±106.2 a | 352.9±86.7 b | 600±84.3 a | 410.37±65.8 ab | F=2.87 | 0.04 |
| BS (g m-²)** | 23.5 (19.83-24.95)a | 16.99 (12.65-20)b | 13.91 (9.72-17.25)b | 9.38 (6.3-11.88)c | H=40.59 | 0.00 |
| BA:BS** | 24.62(19.11-32.62)ab | 15.81(9.56-20.56)b | 39.4(27.31-54.84)a | 43.11(18.55-59.8)a | H=13.58 | 0.00 |
Los cuatro sitios compartieron solo dos especies: Sagittaria lancifolia y Eleocharis cellulosa, dos especies características de humedales. El VIR de la primera especie fue significativamente mayor en RB1 en comparación con los demás sitios, mientras que los valores para E. cellulosa no fueron significativamente diferentes entre sitios. Se encontraron otras especies compartidas en dos o tres sitios, aunque la única especie que presentó diferencias significativas en el VIR fue Thalia geniculata, en donde fue mayor el valor de éste en RL2 que en RB2 (Cuadro 2). Es importante mencionar que hay especies que solo se presentaron en un sitio y en un muestreo y por tanto sus VIR fueron muy bajos, por lo que no se presentan en el Cuadro 2.
Los valores del índice de reemplazo de especies de Whittaker de cada sitio con respecto a RB1 fueron de 1.66 para RB2, 1.74 para RL1 y de 1.91 para RL2. Por tanto, RB2 es más parecido a RB1 y RL2 es el que más difiere de éste. En cuanto a la diversidad de especies, RL2 fue significativamente mayor con respecto a los de más sitios, seguido de RB1, RB2 y RL1. De acuerdo al valor del índice de equidad de Pielou, el sitio RB1 fue en donde la proporción de la diversidad fue más equitativa, mientras que en RB2 y RL2, esta proporción no es homogénea y hay especies dominantes. En cuanto a la riqueza de especies, el sitio RL2 fue significativamente mayor al resto de los sitios (Cuadro 3). Los resultados de BA y BS por sitio serán descritos más adelante.
Variación por temporalidad de secas-lluvias y asociación con especies y sitios de muestreo
En el ACP se observa que los dos primeros componentes explicaron 46% de la varianza acumulada. El primero muestra un gradiente de intensidad de forrajeo y las variables de especies de mayor peso fueron Bacopa monnieri, Echinochloa colona (L.) Link e Hydrocotyle verticillata Thunb. (Asociadas positivamente), así como Hymenachne amplexicaulis, Pontederia sagittata, Ludwigia octovalvis (Jacq.) P.H. Raven y Paspalum sp. (asociadas negativamente). Esta última especie apareció solo en RL2 y en un muestreo. En el segundo el mayor peso correspondió a Cyperus articulatus asociada de manera positiva y Lemna minor L. y Nymphaea ampla (Salisb.) DC. Asociadas negativamente (Fig. 2). Ambas especies solo aparecieron en la temporada de lluvias, y la primera únicamente en RB2.
Figura 2: Asociación de las cuatro localidades de estudio (casos) con 18 especies (variables) con respecto a los dos componentes principales. Los dos componentes principales explican 46% de la variación. Sobre el primer componente se establece un gradiente de intensidad de forrajeo con los sitios de mayor intensidad en el extremo izquierdo (RL2). Dentro de cada sitio las épocas del año tienden a separarse sobre el eje dos en función de la temporalidad de secas-lluvias, excepto para el sitio RL1.
Según el índice de Shannon, la diversidad de especies fue mayor en lluvias en todos los sitios, aunque las diferencias solo fueron significativas en RB2 (Secas H’=0.98, Lluvias H’=1.83) y en RL2 (Secas H’=1.8, Lluvias H’=2.37). El índice de equidad de Pielou reflejó que todos los sitios presentaron un mayor índice en lluvias; sin embargo, solo en RB1 y RB2 las diferencias fueron significativas. La riqueza media de especies fue mayor en lluvias en todos los sitios, aunque solo en RB1 las diferencias no fueron significativas (Cuadro 4).
Cuadro 4: Diversidad de especies (índice de Shannon - H’), riqueza de especies, biomasa aérea (BA) y biomasa subterránea (BS) por temporadas (Secas-Lluvias) en cada uno de los sitios de muestreo en el municipio de Alvarado, Veracruz, México. En el caso del índice de Shannon, se indica el valor de éste (H’) y la muestra (S). Se indica la media ± 1EE al aplicarse la prueba de t pareada y se muestra la mediana y los percentiles (25-75%) en los casos donde se aplicó la prueba de Rango de Wilcoxon. Valores con la misma letra son estadísticamente iguales con p>0.05.
| RB1 | RB2 | |||||
| Característica | Secas | Lluvias | Secas | Lluvias | ||
| Índice de diversidad | S=7, H’=1.57a | S=7, H’=1.81a | z=0.88 p=0.37 | S=5, H’=0.98a | S=10, H’=1.83b | z=2.67 p= 0.00 |
| Índice de equidad | J’=0.92a | J’=0.98b | z=2.53 p=0.01 | J’=0.68a | J’=0.98b | z=2.66 p=0.00 |
| Riqueza de especies | 4.5(3.25-5)a | 5(4-5) a | z=1.7 p=0.08 | 2(1-2) a | 4(3-4)b | z=3.32 p=0.00 |
| BA (g m-²) | 820±196a | 500±48.3a | t=1.52 p=0.18 | 185.37±62.7a | 520.5±142b | t=-2.39 p=0.04 |
| BS (g m-²) | 23.7±1.7a | 21.22±1.3a | t=2.17 p=0.06 | 13.65±2.62a | 19.41±0.76b | t=-3.21 p=0.01 |
| BA:BS | 34.06±7.55a | 24.13±3.34a | t=1.19 p=0.27 | 15.21(0-18)a | 19.25(11.8-57.1)b | z=2,24 p=0.02 |
| RL1 | RL2 | |||||
| Característica | Secas | Lluvias | Secas | Lluvias | ||
| Índice de diversidad | S=8, H’=1.52a | S=6, H’=1.54a | z=1.01 p=0,311 | S=10, H’=1.8a | S=15, H’=2.37b | z=2.09 p=0.03 |
| Índice de equidad | J’=0.89a | J’=0.99a | z=1.19 p=0.23 | J’=0.78a | J’=0.97a | z=1.75 p=0.07 |
| Riqueza de especies | 2(2-2)a | 3(3-4)b | z=3.03 p=0.00 | 2(2-3.75)a | 3.5(3-4)b | z=2.81 p=0.00 |
| BA (g m-²) | 730.7±108.3a | 470.6±118a | t=4.74 p=0.00 | 242.25±84.6a | 578±57.8b | t=-3.54 p=0.00 |
| BS (g m-²) | 14.78±1.43a | 13.46±1a | t=0.55 p=0.59 | 7.82±1.16a | 9.77±1a | t=-1.54 p=0.16 |
| BA:BS | 54.09±8.64a | 35.21±7.38a | t=2.23 p=0.06 | 29.49±9.96a | 62.92±11.51b | t=-2.41p=0.04 |
Las especies que incrementaron su VIR y se vieron beneficiadas en los cuatro sitios durante la temporada de lluvias fueron Sagittaria lancifolia y Eleocharis cellulosa (Cuadro 5). En los sitios RB1 y RB2, la especie que incrementó su VIR en esta temporada fue Bacopa monnieri (RB1 Secas 47.6, Lluvias 68.92; RB2 Secas 30.55, Lluvias 56.36), al igual que lo hicieron Hymenachne amplexicaulis (Secas 52.31, Lluvias 86.9), Cyperus articulatus (Secas 11.34, Lluvias 52.08) en RL1, y Thalia geniculata (Secas 41.89, Lluvias 58.33) en RL2. Se encontraron especies que solo estuvieron presentes durante la temporada de lluvias como el pasto Sporobolus virginicus (L.) Kunth en RB1, y Lemna minor, Nymphaea ampla y Pontederia sagittata en RB2. Durante la temporada de secas la especie más importante en RB1 fue Echinochloa colona (Secas 90.27, Lluvias 35.87), mientras que en RB2 ésta misma se presentó constante durante todo el año con un VIR alto (Secas 85, Lluvias 88). El pasto Leersia hexandra Sw. en RL1 presentó un VIR alto durante todo el año (Secas 67.36, Lluvias 63.64), y en RL2 la gramínea Hymenachne amplexicaulis se encontró asociada a la temporada de secas con VIR más alto que en lluvias (Secas 76.54, Lluvias 68) (Cuadro 5).
Cuadro 5: Valor de importancia relativa (VIR) de cada especie por temporalidad de secas-lluvias para los cuatro sitios de muestreo en el municipio de Alvarado, Veracruz, México. Se indica la media ± 1 EE del Valor de Importancia Relativa. Las especies con asterisco (*) solo se presentaron en una temporada.
| Secas | Lluvias | Secas | Lluvias | ||
|---|---|---|---|---|---|
| Especie | RB1 | Especie | RL1 | ||
| Sagittaria lancifolia L. | 78.04±3.46 | 76.28±6.34 | Eleocharis cellulosa Torr. | 61.88±12.63 | 76.6±7.62 |
| Eleocharis cellulosa Torr. | 59.84±8.38 | 83.53±3.27 | Hymenachne amplexicaulis (Rudge) Ness | 52.31±27.7 | 86.9±1.02 |
| Hydrocotyle verticillata Thunb. | 65.83±1.71 | 72.71±1.98 | Sagittaria lancifolia L. | 21.06±1.85 | 51.59±8.54 |
| Bacopa monnieri (L.) Wettst. | 47.6±2.93 | 68.92±5.56 | Cyperus articulatus L. | 11.34±11 | 52.08±6.9 |
| Sporobolus virginicus (L.) Kunth* | 0 | 51±16.67 | Leersia hexandra Sw. | 67.36±11.46 | 63.64±15.02 |
| Echinochloa colona (L.) Link | 90.27±3.02 | 35.87±30 | |||
| Especie | RB2 | Especie | RL2 | ||
| Echinochloa colona (L.) Link | 85.09±7.39 | 88.19±4.16 | Hymenachne amplexicaulis (Rudge) Ness | 76.54±1.31 | 68.05±8.33 |
| Thalia geniculata L. | 19.21±10.06 | 17.01±16 | Ludwigia octovalvis (Jacq.) P.H. Raven | 50.15±0.9 | 41.66±1.38 |
| Eleocharis cellulosa Torr. | 27.54±14.4 | 64±5.9 | Thalia geniculata L. | 41.89±1.15 | 58.33±1.38 |
| Bacopa monnieri (L.) Wettst. | 30.55±15.36 | 56.36±11.68 | Eleocharis cellulosa Torr. | 38.88±4.86 | 53.47±10.41 |
| Sagittaria lancifolia L. | 13.42±6.9 | 34.02±0.1 | Sagittaria lancifolia L. | 20.83±11.47 | 28.81±5.9 |
| Lemna minor L. * | 0 | 52±12.5 | Pontederia sagittata C. Presl | 35.87±1.85 | 36.8±2.77 |
| Nymphaea ampla (Salisb.) DC. * | 0 | 41.31±1.73 | |||
| Pontederia sagittata C. Presl * | 0 | 31.25±8.33 | |||
Variación entre tratamientos y asociación con las especies y sitios de muestreo
El segundo ACP que se realizó para los tratamientos explicó 32% de la varianza entre los dos primeros componentes siendo éstos los más importantes. Las variables de mayor peso en el primer componente fueron las especies Bacopa monnieri, Echinochloa colona, Hydrocotyle verticillata y Sagittaria lancifolia (asociadas positivamente) y Ludwigia octovalvis, Hymenachne amplexicaulis y Paspalum sp. (Asociadas negativamente). En el segundo componente los mayores pesos correspondieron a Thalia geniculata y Gomphrena hispida L., asociadas positivamente, e Hymenachne amplexicaulis y Eleocharis cellulosa asociadas negativamente (Fig. 3).
Figura 3: Asociación de las cuatro localidades de estudio (casos) con respecto a los dos componentes principales. Dentro de cada localidad hay tratamientos “sin ganado” SG (color gris) y “con ganado” CG (color negro). Como variables se utilizaron 22 especies. Los dos componentes principales explican 32% de la variación.
La riqueza media de especies fue significativamente mayor en el tratamiento “sin ganado” en todos los sitios excepto RB2. No se observaron diferencias en cuanto al Índice de diversidad de Shannon e Índice de equidad de Pielou entre los tratamientos (Cuadro 6). En el ACP se observa que las especies que se asocian al tratamiento “sin ganado” en el sitio RB1 fueron H. verticillata, S. lancifolia y B. monnieri, presentando diferencias significativas en cuanto al VIR solo las dos últimas especies mencionadas; en el sitio RB2 fueron B. monnieri y L. minor; en el sitio RL1 fueron L. hexandra, C. articulatus y E. cellulosa, presentando diferencias significativas en el VIR solo la primera especie; y en el sitio RL2 se observa la presencia exclusiva de H. amplexicaulis en este tratamiento. En cuanto al tratamiento “con ganado”, en RB1 se observa como la gramínea S. virginicus es exclusiva de este tratamiento, mientras que en RB2 E. colona se asocia a la presencia del ganado presentando un VIR mayor, sin embargo, las diferencias no fueron significativas. En RL1, H. amplexicaulis se encuentra asociada a este tratamiento, y en RL2 las especies asociadas al ganado fueron especies como L. octovalvis, G, hispida, C. anguria y Paspalum sp. (Fig. 3, Cuadro 7).
Cuadro 6: Diversidad y riqueza de especies, biomasa aérea (BA) y biomasa subterránea (BS) tomadas en las parcelas de cada uno de los sitios de estudio, de acuerdo al tratamiento ganadería (“sin ganado” y ”con ganado”) en el municipio de Alvarado, Veracruz, México. En el caso del índice de Shannon, se indica el valor de éste (H’) y la muestra (S). Se indica la media ± 1EE al aplicarse la prueba de t pareada (*) y se muestra la mediana y los percentiles (25-75%) en los casos donde se aplicó la prueba de Rango de Wilcoxon (**). Valores con la misma letra son estadísticamente iguales con p>0.05.
| Característica | RB1 | RB2 | ||||
| Sin ganado | Con ganado | Sin ganado | Con ganado | |||
| Índice de diversidad | S=8, H’=1.8a | S=8, H’=1.8a | z=1.19 p=0.23 | S=8, H’=1.5a | S=8, H’=1.5a | z=0.28 p=0.77 |
| Índice de equidad | J’=0.97a | J’=0.96a | z=1.64 p=0.1 | J’=0.95a | J’=0.95a | z= 0.8 p=0.42 |
| Riqueza de especies** | 5(4-5)a | 4(3-5)b | z=2.37 p=0.01 | 2(1-4)a | 2(1-3)a | z=0.08 p=0.93 |
| BA (g m-²) * | 724.2±212a | 595.8±44.3a | t=0.59 p=0.57 | 240.12±52.7a | 365.12±104.78a | t=-1.99 p=0.08 |
| BS (g m-²)* | 26.71±1.99a | 20.29±1.13b | t=4.04 p=0.00 | 18.02±2.48a | 16.15±1.12a | t=0.73 p=0.47 |
| BA:BS** | 19.25(14.2-37.7)a | 28.33(21.7-41)a | z=1.26 p=0.2 | 9.56(0-15.8)a | 20.56(16.4-57)b | z=2.52 p=0.01 |
| Característica | RL1 | RL2 | ||||
| Sin ganado | Con ganado | Sin ganado | Con ganado | |||
| Índice de diversidad | S=7, H’=1.6a | S=7, H’=1.4a | z=1.24 p=0.21 | S=15, H’=2.1a | S=11, H’=2a | z=0.72 p=0.46 |
| Índice de equidad | J’=0.98a | J’=0.983a | t=1 p=0.31 | J’=0.95a | J’=0.95a | z=1.32 p=0.18 |
| Riqueza de especies** | 3(2-3)a | 2(1.2-3)b | z=2.54 p=0.01 | 3.5(2-4)a | 2.5(1-3)b | z=2.88 p=0.00 |
| BA (g m-²) * | 628.8±142.9a | 572.5±99.2a | t=0.27 p=0.78 | 462.12±75.42a | 358.6±110.18a | t=0.89 p=0.40 |
| BS (g m-²)* | 14.77±1.02a | 12.42±1.37a | t=1.45 p=0.16 | 8.3±0.8a | 10.4±1.3a | t=-1.47 p=0.15 |
| BA:BS* | 45.05±11.18a | 53.5±13.93a | t=-0.4 p=0.69 | 56.16±10.16a | 36.25±13.42a | t=1.12 p=0.29 |
Cuadro 7: Valor de importancia relativa (VIR) de las especies más representativas para los cuatro sitios de muestreo en el municipio de Alvarado, Veracruz, México. Estos valores se presentan para el tratamiento ganadería (“sin ganado” y ”con ganado”). Se indica la media ± 1EE al aplicarse la prueba de t pareada y se muestra la mediana y los percentiles (25-75%) en los casos donde se aplicó la prueba de rango de Wilcoxon. Valores con la misma letra son estadísticamente iguales con p>0.05. Las especies que presentan asterisco (*), solo aparecieron en uno de los tratamientos y no se pudo realizar prueba estadística.
| Especie | RB1 | Estadístico | |
|---|---|---|---|
| Sin ganado | Con ganado | ||
| Sagittaria lancifolia L. | 83.33(81.9-83.3)a | 63.42(63.42-81.9)b | z=2.03, p=0.04 |
| Eleocharis cellulosa Torr. | 76.85±7.3a | 60.83±9.31a | t=2.5, p=0.06 |
| Hydrocotyle verticillata Thunb. | 71.66±3.33a | 65.55±3.17a | t=1.23, p=0.28 |
| Echinochloa colona (L.) Link | 79.16(32-93)a | 90.27(39-93)a | z=1.46, p=0.14 |
| Bacopa monnieri (L.) Wettst. | 63.51±5.8a | 48.05±6.86b | t=3.47, p=0.02 |
| Sporobolus virginicus (L.) Kunth* | 0 | 28±39.31 | |
| RB2 | |||
| Echinochloa colona (L.) Link | 80.55±7.77a | 91.94±1.83a | t=-1.61, p=0.18 |
| Eleocharis cellulosa Torr. | 39.72±17.7a | 31.94±13.5a | t=0.41, p=0.69 |
| Bacopa monnieri (L.) Wettst. | 61.11(0-72)a | 0(0-54)a | z=1.09, p=0.27 |
| Thalia geniculata L. | 12.77±8.25a | 12.77±8.25a | t=0, p=1 |
| Pontederia sagittata C. Presl* | 15±9.55 | 0 | |
| Lemna minor L.* | 20.5±12.59 | 0 | |
| RL1 | |||
| Eleocharis cellulosa Torr. | 67.91±6.42a | 45.27±18.69a | t=1.2, p=0.29 |
| Hymenachne amplexicaulis (Rudge) Ness | 47.5±20.02a | 67.22±17.81a | t=-1.16, p=0.31 |
| Leersia hexandra Sw. | 73.05±9.51a | 40±17.1b | t=3.45, p=0.02 |
| Cyperus articulatus L. | 22.45±14.18a | 18.05±11.16a | t=0.31, p=0.76 |
| Sagitaria lancifolia L.* | 36.8±6.06 | 0 | |
| RL2 | |||
| Hymenachne ampllexicaulis (Rudge) Ness | 93(81.25-97)a | 0(0-59.7)b | z=2.02, p=0.04 |
| Ludwigia octovalvis (Jacq.) P.H. Raven | 52.77(0-55.5)a | 56.94(51.3-67.36)a | z=1.46, p=0.14 |
| Ambrosia psilostachya DC. | 23.33±14.52a | 26.29±16.73a | t=-1.27, p=0.32 |
| Paspalum sp. | 23.33±14.31a | 29.44±18.16a | t=-1.12, p=0.32 |
| Pontederia sagittata C. Presl * | 42.5±1.36 | 0 | |
| Sagitaria lancifolia L.* | 30.27±8.5 | 0 | |
| Eleocharis cellulosa Torr.* | 51.11±4.5 | 0 | |
Biomasa por sitio, por temporalidad de secas-lluvias y tratamiento
En el sitio RB2 la BA (352.9 g m-²) fue significativamente menor que en RB1 (660 g m-²) y en RL1 (600g m-²). Con respecto a la BS, en RB1 (23.5g m-²) hubo mayor cantidad de ésta con respecto a los demás sitios y RL2 (9.385g m-²) resultó menor estadísticamente a los sitios restantes (Cuadro 3). Con respecto a la relación BA:BS, RB2 (15.81) fue significativamente menor en comparación con RL1 (39.4) y RL2 (43.11). En cuanto a la temporalidad, en el sitio RB2 en lluvias la cantidad de BS, BA y la relación BA:BS fueron mayores. En RL2 solo BA y la relación BA:BS fueron mayores estadísticamente en esta misma temporada (Cuadro 4). En los tratamientos “con ganado” y “sin ganado”, no hubieron diferencias significativas en cuanto a la BA; sin embargo, en cuanto a la BS, ésta fue mayor en las parcelas “sin ganado” en RB1 (Cuadro 6).
Discusión
La composición florística de los cuatro humedales herbáceos de Alvarado en la costa centro oeste del Golfo de México, estudiados en el presente trabajo, se está viendo afectada negativamente por la carga ganadera (dos cabezas de ganado por ha o más). Los resultados indican que el gradiente de transformación encontrado entre los sitios de muestreo se debe a las diferencias de manejo, tanto en carga animal como en el tiempo de permanencia en el predio, y esto a su vez causa alteración en la estructura vegetal. Se observa que RB1, sitio con menor carga animal, fue el más conservado debido a que las especies presentes son principalmente nativas del humedal estuvieron constantes a lo largo del año con altos VIRs (arriba de 50) y presentaron equidad en cuanto a la proporción de su abundancia según el índice de Pielou (0.95) (Cuadro 3). Las especies a las que nos referimos y que se han documentado como propias de humedales son Bacopa monnieri, Hydrocotyle verticillata, Sagittaria lancifolia y Eleocharis cellulosa (Moreno-Casasola et al., 2010). Estas últimas dos fueron las únicas especies compartidas en los cuatro sitios de muestreo, y principalmente Sagittaria lancifolia se encontró en mayor abundancia en las localidades menos perturbadas. Se ha reportado que esta especie es sensible a la presencia de algunas gramíneas exóticas como Echinochloa pyramidalis (Lam.) Hitchc. & Chase (López Rosas et al., 2010), la cual no se encuentra en nuestros sitios de muestreo, pero sí se reportó a Echinochloa colona que podría estar jugando el mismo papel, ya que está considerada como una especie exótica e invasora (Hanan y Mondragón, 2009). Durante la temporada de secas el índice de equidad de Pielou fue significativamente menor que en lluvias, lo que indica que hay una especie dominante la cual es Echinochloa colona según los resultados del VIR por temporalidad (Cuadro 4 y 5). Sin embargo, el crecimiento de esta especie puede estar siendo controlado por la entrada de agua durante la temporada de lluvias y el régimen de inundación que permanece constante durante el año, y por tanto no se ve afectado el crecimiento de Sagittaria lancifolia. Esto coincide con los resultados del trabajo de López Rosas y Moreno-Casasola (2012), quienes encontraron que la competencia entre Echinochloa pyramidalis y Sagittaria lancifolia era asimétrica en diferentes condiciones de inundación. Por ejemplo, en cuadrantes en donde el suelo se encontraba seco se veía favorecido el crecimiento de Echinochloa P. Beauv. y en cuadrantes en donde el suelo se mantenía inundado, el crecimiento de la gramínea era limitado y dominaba Sagittaria L.
El sitio RB2 fue el más conservado después de RB1. La mayoría de las especies de este sitio son propias de un humedal herbáceo; sin embargo, el pasto exótico Echinochloa colona fue el único que presentó un VIR alto durante todo el año (86%) en comparación a todas las demás que no obtuvieron más de 42%. Es probable que las condiciones de disturbio de este sitio (pisoteo del ganado y sequía) estén favoreciendo el crecimiento de esta especie (Mejía-Saulés y Aranda, 1992) lo que limita el crecimiento de las especies propias del humedal, que durante la época de lluvias claramente se ven favorecidas (López Rosas et al., 2010). Los sitios RL1 y RL2 son los que más difieren de un humedal herbáceo conservado. En ambos sitios la carga animal es mayor y además el ganado permanece por más tiempo. Esto podría tener como consecuencia una mayor compactación del suelo y modificaciones en la hidrología y por tanto se favorece el crecimiento de especies exóticas y se inhibe el de las especies nativas del humedal. Además, los sitios RB1 y RB2 contienen más materia orgánica debido al tipo de suelo, mientras que RL1 y RL2 carecen de ésta (Rodríguez-Medina y Moreno-Casasola, 2013).
Además del manejo ganadero, otro factor que influye en el nivel de impacto y en el cambio de la vegetación es la variación estacional. En dos de los cuatro sitios se observa mayor diversidad y en tres sitios mayor riqueza vegetal en lluvias (Cuadro 4), lo que sugiere que el régimen de inundación está determinando el crecimiento de especies propias del humedal (Hess et al., 2003). Sin embargo, algunos autores mencionan que, tanto en humedales como en otros ecosistemas, la riqueza de especies aumenta al inicio de la temporada de lluvias, debido al crecimiento de especies facultativas que aprovechan situaciones medias de disturbio, cambios en el hidroperíodo y nuevas condiciones de humedad. Esto podría estar sucediendo en RB2 y RL1, pero sobre todo en RL2, que fue el que presentó mayor riqueza de especies. En estos sitios, al inicio de la temporada de lluvias se registra un incremento en el número de especies no nativas del humedal; sin embargo, conforme avanza la temporada de lluvias, los niveles de inundación se elevan y la riqueza de especies disminuye debido a que muchas plantas ruderales no soportan las condiciones de inundación, y generalmente se mantienen especies que son tolerantes a estos ambientes, como son las hidrófitas (Pollock et al., 1998; Weiher, 2003; Travieso-Bello et al., 2005). En humedales perturbados en donde la hidrología ha sido modificada y el sitio permanece menor tiempo inundado (como el sitio RL2), la riqueza de especies nativas decrece e incrementa la de especies facultativas que suelen estar presentes durante todo el año, únicamente con variaciones en su abundancia (Grace y Pugesek, 1997; Travieso-Bello et al., 2005).
Algunas de las especies nativas de los humedales y otras introducidas desarrollan estrategias para tolerar las condiciones de inundación o sequía. En el presente trabajo es importante conocer algunas de estas adaptaciones para entender la permanencia o ausencia de las especies en los sitios. Por ejemplo, los estudios realizados por Busch et al. (2004) y Baksh y Richards (2006) mencionan que, bajo condiciones críticas de inundación, los brotes jóvenes de Eleocharis cellulosa al entrar en contacto con el agua crecen más rápido y producen tallos más gruesos y altos, generando una mayor biomasa aérea en relación con la biomasa subterránea. En secas, estos brotes mueren, pero la especie puede soportar largos períodos de estiaje y producir nuevos brotes, aunque con crecimiento más lento. El hecho de que Eleocharis cellulosa haya presentado un VIR bajo en todos los sitios en la temporada de secas comparado con lluvias puede deberse a lo mencionado anteriormente y también a que ésta es una especie de mayor palatabilidad para el ganado, el cual tiene mayor presencia en secas (Rosen et al., 2008). También se encontraron en los sitios tres especies de gramíneas dominantes: Hymenachne amplexicaulis, Leersia hexandra y Echinochloa colona. La primera es una especie nativa pero invasora (especie translocada) (Hill, 1996; Houston y Duivenvoorden, 2003; Gordon et al., 2007) y se utiliza como forrajera de alta calidad. Su éxito reproductivo se basa en que es estolonífera con culmos que se dispersan sobre el suelo húmedo o flotan en la superficie del agua, a partir de los cuales producen raíces y forman densas colonias (Mejía-Saulés y Aranda, 1992), requiriendo de la temporada de secas para establecerse (Hill, 1996). Los sitios RL1 y RL2 presentan ambos escenarios, y en secas la planta logra establecerse ya que el sitio no se encuentra saturado en agua debido a la compactación del suelo originada por la ganadería, mientras que en lluvias el sitio se inunda logrando que la planta se disperse y se reproduzca con éxito. Esto explica por qué Hymenachne amplexicaulis fue tan exitosa en RL2 y RL1. En este último sitio, además de esta especie también se presentó Leersia hexandra, especie nativa pero invasora, que además es un pasto forrajero de primera calidad utilizado en el pastoreo, muy común en zonas inundables y puede tolerar hasta un metro de inundación (Mejía-Saulés y Aranda, 1992). Esta especie presentó un VIR alto durante todo el año, lo que sugiere que es muy plástica, con capacidad de adaptación a la sequía e inundación.
La tercera gramínea presente en los sitios fue Echinochloa colona, un pasto exótico que habita en lugares húmedos de regiones tropicales y subtropicales, y que a diferencia de Hymenachne amplexicaulis solo tolera parcialmente las inundaciones. Prefiere la temporada de estiaje y tiene la capacidad de adaptarse a suelos pobres en nutrientes (Mejía-Saulés y Aranda, 1992; Calderón y Rzedowski, 2004). Esta especie estuvo presente en RB1 y RB2 siendo más importante en el segundo sitio, en donde no se presentaron diferencias en el VIR entre secas y lluvias. Esto se puede deber a que es un sitio que no está inundado o saturado todo el año y además la carga animal es elevada, lo que se traduce en mayor pisoteo del ganado y mayor compactación. Esto podría indicar que es un suelo pobre en nutrientes, lo cual favorece el establecimiento de Echinochloa colona durante todo el año. A diferencia de RB2, en RB1 el VIR de esta especie fue mucho más alto en secas que en lluvias, lo que coincide con lo mencionado anteriormente por Mejía-Saulés y Aranda (1992) de que tolera parcialmente las inundaciones y prefiere la época de estiaje.
Se ha reportado que en los humedales de agua dulce del sureste de México la introducción de gramíneas y las prácticas ganaderas han causado pérdida de la biodiversidad, ya que muchas de las especies nativas son reemplazadas por las introducidas y las prácticas ganaderas no les permiten mantenerse en el sitio (Travieso-Bello et al., 2005; López Rosas y Moreno-Casasola, 2012). De acuerdo a nuestros resultados esta situación puede estar presentándose en RB2, RL1 y RL2, ya que la ganadería ha favorecido que las tres gramíneas antes mencionadas (una de ellas exótica) dominen sobre las especies propias de los humedales.
Variación por tratamiento (“con ganado” y “sin ganado”)
Se encontraron diferencias significativas entre sitios en cuanto a la carga animal, resultando ser RB1 (sitio más conservado y con poca carga animal), el sitio que presentó mayores características propias a las de un humedal conservado. Sin embargo, en cuanto al tratamiento las diferencias solo fueron significativas en la riqueza de especies y en BS en RB1. Por ejemplo, en las localidades de Río Limón en donde la riqueza fue significativamente mayor al excluir al ganado, se puede observar que las especies que se vieron favorecidas fueron básicamente las gramíneas y algunas ruderales como Leersia hexandra, Hymenachne amplexicaulis, Ambrosia psilostachya DC. y Paspalum sp., esto se debe a que, en áreas muy impactados por esta actividad, los humedales han perdido características intrínsecas del ecosistema, sobre todo en el componente suelo que es primordial para que la composición florística y la hidrología se mantengan. Las características del suelo que son afectadas por la ganadería son la capacidad de retención de agua, el espacio aéreo (porosidad), cantidad de materia orgánica y de micro y macro nutrientes (Rodríguez-Medina y Moreno-Casasola, 2013). Estas alteraciones traen como consecuencia que, aunque se excluya al ganado de la zona, las condiciones inicialmente son más favorables para especies dominantes como son las gramíneas y ruderales, limitando el crecimiento de la vegetación propia del humedal (Bantilan-Smith et al., 2009).
Biomasa
El hecho de que la cantidad de BA en RB1 fuera estadísticamente mayor en comparación con los otros sitios se debe a que muchas especies propias de humedales herbáceos, al verse sometidas a situaciones de anaerobiosis por inundación constante (sitio RB1 casi inundado todo el año), aumentan su tamaño (por ejemplo, tallos más anchos y altos, y/u hojas más anchas), con lo que generan más aerénquima y pueden captar más oxígeno y transportarlo hacia las raíces (Moreno-Casasola y Wagner, 2009; Mitsch et al., 2009). Por ejemplo, López-Rosas y Moreno-Casasola (2012) encontraron que especies propias del humedal, como Sagittaria lancifolia y Typha domingensis, generaban una cantidad estadísticamente mayor de biomasa en condiciones de inundación que durante la temporada de secas, y observaron lo contrario con el pasto invasor Echinochloa pyramidalis. En el presente trabajo no se midió la biomasa por especie; sin embargo, en RB1 la especie de mayor VIR fue Sagittaria lancifolia (probablemente debido a las condiciones de inundación), lo que podría explicar las altas cantidades de BA (Cuadro 2).
En el caso de RL1 y RL2, donde también se obtuvieron grandes cantidades de BA, las especies nativas de los humedales no fueron dominantes, lo que se explica por la variación del hidroperíodo y a la elevada carga animal. Sin embargo, en estas localidades el pasto Hymenachne amplexicaulis fue abundante y se ha reportado que esta especie genera gran cantidad de biomasa aérea (Gordon y Feo, 2007), al igual que Leersia hexandra, gramínea que también estuvo presente en RL1 (Calderón y Rzedowski, 2004).
En cuanto a los tratamientos, se observa que en RB1, RL1 y RL2, la cantidad de BA fue mayor en el tratamiento “sin ganado”; sin embargo, las diferencias no son significativas estadísticamente. El incremento de BA en ausencia de ganado se debe a que la ausencia de herbivoría acumula mayor biomasa vegetal seca y en descomposición, y a que la fitomasa viva no es consumida por los herbívoros (Olff y Ritchie, 1998; Travieso-Bello et al., 2005; Menghi et al., 2011). Es importante mencionar que en el primer muestreo del tratamiento “sin ganado” del sitio RB2 hubo vandalismo y se robó el alambre de púas que rodeaba los cuadros, permitiendo el acceso de ganado al tratamiento y la presencia de herbivoría. Esto hecho pudo influir en la cantidad de BA encontrada en este sitio y en los resultados entre tratamientos. Otro aspecto importante a considerar es que la BA que se obtuvo al inicio del experimento fue cortada dentro de los cuadros a muestrear, y esto pudo tener un efecto inicial de herbívora, e influir en los resultados entre tratamientos. Sin embargo, en cuanto a temporadas (lluvias-secas) no afectó ya que sí se observaron diferencias estadísticamente significativas. La biomasa subterránea fue significativamente mayor en RB1 (sitio permanentemente saturado o inundado). Esto es común que suceda en humedales herbáceos en buen estado de conservación, ya que cuando algunas plantas de estos ecosistemas se encuentran en situaciones de anaerobiosis, la disponibilidad de nutrientes es limitada, estimulando que las plantas generen mayor cantidad de biomasa en las raíces, aumentando así su área y la capacidad de absorción (Ponnamperuma, 1972; Chapin, 1980; Trought y Drew, 1980). Además, en condiciones de anaerobiosis la descomposición de la materia orgánica es lenta, generando un incremento en la acumulación de raíces muertas en proceso de descomposición en el suelo (Richardson y Vepraskas, 2001). El hecho de que RL1 y RL2 fueran los sitios con menor cantidad de BS posiblemente se deba a que son sitios con mucha carga animal y por consecuencia con suelos compactados, por lo que las raíces tienen limitado su crecimiento (Unger y Kaspar, 1994) (Cuadro 3).
Como ya se mencionó anteriormente, la ganadería es una actividad que va en aumento debido a la demanda productiva que existe, al aumento del precio de la carne de res y a las necesidades económicas de los pobladores, lo que se ve reflejado en un incremento de terrenos destinados a esta actividad. Ecosistemas como los humedales que antes no se veían afectados, ahora lo están. Esta actividad de manera intensiva, además de perturbar negativamente la composición florística, afecta otros componentes sumamente importantes para el buen funcionamiento de los humedales, como son las propiedades fisicoquímicas del suelo y la hidrología (Rodríguez-Medina y Moreno-Casasola, 2013) y por tanto también se afectan servicios ecosistémicos dependientes de la porosidad del suelo como son el almacenamiento de carbono y agua en este componente (Campos et al., 2011). Sin embargo, existen estudios que mencionan que en diversos ecosistemas, incluyendo los humedales, la ganadería moderada, específicamente la herbivoría, funciona como control de especies dominantes (principalmente gramíneas) y favorece el crecimiento de las especies nativas. Con los resultados observados en este trabajo, coincidimos con esta idea de que los humedales herbáceos de agua dulce pueden conservar sus funciones ecológicas en presencia de la ganadería, siempre y cuando ésta sea de bajo impacto. Además, creemos que esta actividad de manera controlada puede favorecer a la diversidad de especies nativas de los humedales y limitar el crecimiento de los pastos forrajeros introducidos (Marty, 2005; Junk y da Cunha, 2012; Junk et al., 2014). Un ejemplo de esto es el humedal del sitio RB1. La idea sería moderar dos aspectos: 1) el tiempo en el que el ganado se encuentra en los humedales y 2) la carga animal. Podría considerarse también rotar el ganado más de una vez al año con el objetivo de dejar “descansar” el terreno y dar oportunidad a que la vegetación nativa se recupere. Otra práctica de manejo que se podría combinar con la ganadería de bajo impacto son pequeñas quemas controladas de bajas temperaturas (asegurando que no hay fuerte acumulación de material seco), ya que se ha reportado que se utilizan para frenar la reproducción vegetativa de especies invasoras (López Rosas et al., 2010). Sin embargo, habría que profundizar más en esta área y evaluar qué tanto estas pequeñas quemas podrían contribuir más a la recuperación del humedal o en su defecto perjudicarlo. En municipios como Tecolutla, Veracruz, se ha reportado de manera informal por los pobladores, que anualmente hay fuegos espontáneos en tulares, que son atendidos por los habitantes o autoridades correspondientes para evitar se propaguen a cultivos o caseríos, y que estas quemas no han afectado hasta el momento la presencia de las especies dominantes del humedal.
Concluimos de acuerdo a nuestros resultados obtenidos que: 1) la ganadería está teniendo un efecto negativo en tres de los humedales herbáceos muestreados y esto está relacionado con la carga animal y el tiempo de inundación. Con respecto a los tratamientos “con y sin ganado” por sitio en general no se observaron diferencias estadísticamente significativas; 2) La variación estacional está determinando la composición florística de cada sitio e influyendo en el nivel de impacto ganadero. Por tanto, retomamos la hipótesis y consideramos que una ganadería con una carga de una cabeza por hectárea permite conservar la composición de especies nativas del humedal y un crecimiento en biomasa aérea y subterránea de esta vegetación comparable a la de un humedal sin impacto ganadero. Los autores están conscientes de las limitaciones que el diseño experimental tiene para los resultados, por lo que se sugiere que para futuros estudios el tiempo de muestreo se extienda a dos o más años, para que las diferencias entre tratamientos con respecto a la composición florística y la biomasa sean más evidentes. También se recomienda que los cuadros cercados para excluir al ganado sean de mayor tamaño que 1 m², ya que con parcelas tan pequeñas hubo una fuerte presión de la ganadería externa. Cabe decir que este diseño podrá implementarse gracias a la información obtenida en el presente estudio. El hecho de no poder encontrar parcelas sin ganado, y que haya sido complicado encontrar parcelas sin menos de dos cabezas de ganado por ha, aun cuando en muchas zonas rurales del país con vegetación terrestre hay abandono de terrenos por ser improductivos, muestra la presión que está ejerciendo la ganadería en estos humedales. En la zona, a pesar de ser un sitio Ramsar, no hay Áreas Naturales Protegidas ni estaciones de campo donde se puedan montar experimentos a más largo plazo. Los resultados obtenidos pueden ser la base para convencer a algunos ganaderos locales de permitir usar una mayor superficie de su terreno por varios años, y contar con su apoyo para la protección de materiales asegurando que se impida la entrada del ganado.
Estos resultados, junto con los obtenidos en Brasil, son los primeros que se enfocan a obtener información sobre las implicaciones de la actividad ganadera en los humedales tropicales, una actividad ampliamente difundida en todos los continentes. Nuestros datos indican la necesidad de llevar a cabo experimentos más extensos, en distintos tipos de humedales y bajo diversos tipos de manejo, para poder dar recomendaciones que aseguren la sustentabilidad del funcionamiento del humedal, así como de una actividad económica fundamental para las zonas rurales del trópico húmedo.
