INTRODUCCIÓN
El río Zahuapan es receptor de aguas residuales de origen urbano e industrial en el estado de Tlaxcala. Como consecuencia, la contaminación del río es notable (García-Nieto et al. 2011). Desde una perspectiva ecológica y de manejo de cuencas, la contaminación es solamente un aspecto del problema del deterioro del río. Se debe considerar también a la pérdida del hábitat ribereño y su biodiversidad, así como la del acuático, y la alteración del sistema hidrológico (caudales). Los conceptos de condición e integridad del río permiten trascender la visión de centrar la atención solamente en la contaminación. Este estudio tuvo como objetivo evaluar la condición e integridad del Río Zahuapan. Se tomó como base la condición de hábitat y parámetros físico - químicos de calidad del agua, mediante el protocolo de biomonitoreo rápido para ríos y corrientes basado en la presencia de macroinvertebrados bentónicos (Barbour et al. 1999).
En una evaluación integral del deterioro de un río son centrales los conceptos de cuenca y condición e integridad ecológica del río. La cuenca permite considerar la interacción de los componentes físicos, biológicos y socioeconómicos, así como tomar en cuenta relaciones espaciales de los procesos, desde los sitios de las partes altas a las partes bajas (Frissell et al. 1986). La integridad ecológica de los ecosistemas acuáticos se define en relación con una condición de referencia, donde los impactos de las actividades humanas son mínimos o imperceptibles. A partir de ahí se evalúa el grado en el cual los componentes físicos, químicos y biológicos de un ecosistema y sus relaciones están presentes, funcionan y se mantienen (Schallenberg et al. 2011). Integridad ecológica, salud del río y condición ecológica son términos que se consideran a veces como equivalentes en estudios de monitoreo de los ríos con el propósito de medir una condición total de los ecosistemas acuáticos, aunque para algunos esto es controversial (Oliveira y Cortes 2006).
El biomonitoreo o monitoreo biológico utiliza organismos bioindicadores o sus respuestas para determinar la condición o cambio del ambiente o de un ecosistema. A través de él se evalúan los efectos de la influencia humana sobre un ecosistema acuático (Karr 2006). Un bioindicador es un organismo (o parte de un organismo) que contiene información de la calidad del ambiente (Li et al. 2010). Para la evaluación de ríos, los macroinvertebrados bentónicos, las algas y los peces son los más utilizados.
Los métodos de biomonitoreo incluyen el uso de índices de diversidad, índices bióticos, métodos multimétricos, métodos multivariados y enfoques funcionales (Richards y Horn 2005, Li et al. 2010). El índice de integridad biótica (IBI, por sus siglas en inglés) ha sido ampliamente utilizado para evaluar la integridad o condición de los ríos (Harrison y Whitfield 2006). Barbour et al. (1999), desarrollaron un protocolo de biomonitoreo rápido para usarse en ríos y corrientes, basado en la presencia de perifiton, macroinvertebrados bentónicos y peces. Es un método multimétrico que evalúa también la calidad de la estructura del hábitat e incluye aspectos físicos y químicos del agua del río.
En México, los primeros estudios sobre bioindicadores de ríos se realizaron en los años 90. Lyons et al. (1995) desarrollaron un índice de integridad biótica para la región centro-noroeste de México, basado en comunidades de peces, a partir de 10 parámetros considerados en índices y que han sido utilizadas en regiones de Norteamérica. Ellos evaluaron también la calidad ambiental de los sitios a partir de variables del hábitat y de la calidad de agua, estandarizando los datos a una escala de 10 puntos y considerando el valor más bajo obtenido como indicador de la calidad ambiental menor. Para la misma región, Weigel et al. (2002) utilizaron ocho medidas para desarrollar un índice de integridad biótica basado en macroinvertebrados. Para la condición ambiental, evaluaron la perturbación como mínima, moderada y severa. Otros estudios orientaron el monitoreo a la evaluación de índices de diversidad (Godínez-Domínguez y González-Sansón 1999). Los índices bióticos utilizados inicialmente en México identificaron taxonómicamente a los macroinvertebrados hasta el nivel de género, aunque recientemente se ha considerado la identificación a nivel de familia como una opción, por la facilidad de identificación (Gamboa et al. 2008).
MATERIALES Y MÉTODOS
Área de estudio
El estudio se realizó en la subcuenca del río Zahuapan (Fig. 1), la cual tiene una superficie de 1632.7 km2. El clima predominante es templado subhúmedo y solamente en las áreas con mayor altitud se presenta un clima semifrío subhúmedo. En la zona de estudio los tipos de vegetación que predominan son el pastizal, el bosque de táscate, el bosque de pino y el bosque de encino. En la subcuenca está asentada el 61.18 % de la población del estado de Tlaxcala (INEGI 2010) y se tiene la presencia de una zona industrial importante.
Identificación de factores de perturbación
Las zonas de agricultura de riego, consideradas como fuentes no puntuales de contaminación, se delimitaron con la cartografía digital de uso de suelo vegetación, serie IV del Instituto Nacional de Estadística y Geografía. Se precisó la delimitación de uso de suelo con fotografías aéreas y se realizaron recorridos de campo para validarla. Los sitios de descarga de aguas residuales (fuentes puntuales de contaminantes) se ubicaron con el Localizador del Registro Público de Derechos de Agua (LOCREPDA),un sistema de información de la Comisión Nacional del Agua (CONAGUA 2012). Adicionalmente, a través de un sistema de información geográfica se relacionaron los sitios de descarga con el origen de las mismas: industrial, público-urbano, servicios y pecuario.
Sitios de muestreo
Para evaluar la condición del río se seleccionaron once sitios de muestreo (Fig. 1, S1 a S11), dos de ellos fueron sitios de referencia, localizados en ecosistemas forestales (S1 y S2), donde la perturbación por actividades humanas es mínima. Los demás sitios se ubicaron en áreas con influencia de alguno de los factores de perturbación de la condición del río: agrícola-pastizal (S3), agrícola (S4 y S11), urbano (S5, S9 y S10) e industrial y urbano-industrial (S6, S7 y S8). En todos los sitios se evaluaron en la primer quincena del mes de mayo de 2013 tres componentes: a) macroinvertebrados presentes, b) condición de hábitat y c) parámetros de calidad del agua. Se seleccionó un periodo de muestreo previo al establecimiento del periodo lluvioso para evitar alteraciones del caudal y de la presencia de macroinvertebrados.
Muestreo de macroinvertebrados y condición biótica
Se colectaron los macroinvertebrados en tramos de 100 m del río, con una red de marco rectangular de 0.5 × 0.3 m y apertura de malla de 500 µm. En cada sitio se realizaron 20 golpes con la red sobre el lecho del río para dispersar y colectar a los organismos, de forma proporcional al tipo de hábitat, de acuerdo con el método propuesto por Barbour et al. (1999). Por ejemplo, si en un sitio el 50 % del hábitat era de cantos rodados, entonces se daban 10 golpes con la red y el resto de los golpes se distribuían en los hábitats restantes (tocones, bancos con vegetación, macrófitas sumergidas y arena). Los organismos se colocaron en una charola de plástico marcada con cuadrículas de 6 × 6 cm, numeradas del 1 al 36. Con una tabla de números aleatorios se seleccionó una cuadrícula para colectar los organismos presentes repitiéndose el procedimiento hasta obtener 250 organismos. Los organismos se identificaron hasta el nivel de familia, utilizando un microscopio estereoscópico y las claves de identificación de Bouchard (2004), Oscoz (2009), Serrano y Zepeda (2010) y Gutiérrez-Fonseca (2010).
Se determinó la tolerancia de los organismos identificados a la perturbación de la calidad del agua del río y se identificó el grupo funcional de alimentación al que pertenecen (raspadores, depredadores, colectores, filtradores o fragmentadores), de acuerdo con los cuadros reportados por Bouchard (2004). Con las familias identificadas y el número de individuos por familia se calcularon los valores para cada sitio: número total de familias; número de familias de los órdenes Ephemeroptera, Plecopera, Trichoptera; porcentaje compuesto de Ephemeroptera, Plecoptera, Trichoptera (% de EPT); porcentaje de Ephemeroptera; porcentaje de organismos tolerantes; Índice Biótico de Hilsenhoff (IBH); porcentaje de la familia presente más dominante (% de familia dominante); porcentaje de raspadores; porcentaje de depredadores; porcentaje de colectores; porcentaje de filtradores y porcentaje de fragmentadores. El IBH en cada sitio se calculó con la fórmula: IBH = ƩXiTi/n, donde Xi es el número de ejemplares de cada familia; Ti corresponde al valor de tolerancia para dicha familia y n es el número total de organismos recolectados en la muestra.
Un criterio de selección de medidas para evaluar la condición biótica fue el análisis de correlación de Spearman entre las medidas calculadas. Un segundo criterio fue la relación de medias entre las medidas de los sitios perturbados y los sitios de referencia. Como tercer criterio se graficaron los índices bióticos con los parámetros de perturbación: conductividad eléctrica, cubrimiento de cantos rodados por sedimentos finos y ancho de la zona riparia (Herbst y Silldorf 2009). Se utilizó el cuadro I (Plafkin et al. 1989) para obtener un puntaje de los índices bióticos. Para el cálculo se utilizó el sitio S1 como referencia.
(a) valor del sitio impactado/valor de referencia x 100, (b) valor del sitio de referencia/valor del sitio impactado × 100, (c) la determinación del grupo funcional de alimentación es independiente de la agrupación taxonómica, (d) porcentaje estimado para cada sitio perturbado, sin comparar con el de referencia, se ubica el valor del índice biótico en los rangos del cuadro. EPT = compuesto de Ephemeroptera, Plecoptera, Trichoptera.
Fuente: Plafkin et al. (1989)
Se obtuvo la suma de puntaje de todos los índices bióticos para cada sitio. Se calculó el porcentaje que representó el puntaje de los sitios perturbados en relación con el puntaje del sitio S1. El valor obtenido (% de referencia) se comparó con los valores del cuadro II para determinar la condición biótica del sitio.
EPT = compuesto de Ephemeroptera, Plecoptera, Trichoptera. El índice EPT es la relación % de EPT en el sitio impactado / % de EPT en el sitio de referencia.
Fuente: Plafkin et al. (1989).
Parámetros de hábitat y condición de hábitat
Respecto a la condición de hábitat, para cada parámetro listado en el cuadro III se asignó, a través de una evaluación visual, un valor de condición de hábitat de acuerdo con los lineamientos descritos por Barbour et al. (1999). Los rangos de condición para cada parámetro se describen en el cuadro III. Para cada sitio se sumaron los valores de todos los parámetros y así se obtuvo una calificación agregada.Esta calificación se dividió entre 200, que es el valor máximo posible y se utilizaron los rangos del Departamento de Recursos Naturales de Vermont (VANR 2004) para evaluar la condición de hábitat de cada sitio. Estos rangos son: óptima (0.85 - 1.0), subóptima (0.65 - 0.84), marginal (0.35 - 0.64) y pobre (0.0 - 0.34).
Se realizó un análisis de conglomerados, utilizando el programa SPSS (Statistical Package for Social Sciences) versión 17, con los valores de condición de hábitat para identificar la agrupación de los sitios y observar si los sitios de referencia se separan de los sitios perturbados. Con base en el dendograma obtenido se identificaron grupos con diferente condición de hábitat.
Medición de parámetros físicos y químicos
Se determinaron los parámetros de pH, total de sólidos disueltos, fósforo total, temperatura del agua, oxígeno disuelto, demanda biológica de oxígeno y amonio. La temperatura del agua y el oxígeno disuelto se midieron en los sitios de muestreo con un oxímetro marca Hanna modelo HI 8043. En laboratorio se determinaron el pH y el total de sólidos disueltos con un multiparámetro de calidad del agua marca Hach modelo Sension 156, el fósforo total por el método de PhosVer(r) 3 con digestión de persulfato ácido, el amonio por el método de salicilato y la demanda biológica de oxígeno (DBO5) por el método respirométrico de acuerdo con el procedimiento Hach GGA (glucosa/ácido glutámico), equipo BODTrack II (HACH 2003).
Se calculó el índice de calidad del agua (ICA) propuesto por el Consejo Ambiental de Ministros de Canadá (CCME 2001):
El valor de 1.732 normaliza los valores resultantes a un rango de 0 a 100, donde 0 representa la peor condición de calidad de agua y 100 la mejor condición. F1 es el porcentaje de parámetros que no cumplen con un valor de referencia (número de parámetros que no cumplen con el valor de referencia / número total de variables × 100). F2 es el porcentaje de valores individuales de los parámetros que no cumplen con el valor de referencia (número de valores de los parámetros que no cumplen con el valor de referencia / total de valores de los parámetros medidos × 100). F3 es la amplitud con la que no se cumplió con el valor de referencia, calculada como una función asintótica que escala la suma normalizada de desviaciones (sne) para dar un rango de 0 a 100, mediante la siguiente ecuación: F3 = [sne / (0.01sne + 0.01)]. Sne se estima como: sne = [∑ei] / número de pruebas.
Una desviación (ei) es el número de veces que el valor individual de un parámetro fue mayor (o menor, si el valor de referencia es mínimo) que el máximo (mínimo) permisible, de acuerdo con el objetivo establecido (la protección de la vida acuática). Para el primer caso e = (Valor del parámetro i que no cumple con el valor de referencia / valor de referencia j para el parámetro i) -1. Para el segundo caso e = (valor de referencia j para el parámetro i / valor del parámetro i que no cumple con el valor de referencia) -1.
Debido a que en México no existen valores de referencia para los parámetros físicos y químicos de calidad de agua en ríos, el fósforo total, la demanda biológica de oxígeno (DBO5) y la temperatura se compararon con los valores de la NOM-001-SEMARNAT-1996, la cual establece los límites máximos permisibles de descargas de aguas residuales en aguas y bienes nacionales. Para el oxígeno disuelto, amonio y pH se consideraron valores de la Norma Canadiense para la Protección de Vida Acuática (CCME 2001, CCME 2007). El valor de referencia del total de sólidos disueltos en ríos se basó en los valores establecidos por el Departamento de Recursos Naturales del Estado de Iowa (IDNR 2009). La asignación de la condición del río se realizó conforme a los rangos de la Norma Canadiense para la Protección de Vida Acuática (CCME 2001), descritos en el cuadro IV.
RESULTADOS Y DISCUSIÓN
Factores de perturbación de la integridad ecológica
Las fuentes puntuales de contaminación identificadas en la subcuenca son sitios de descarga de aguas residuales de origen industrial, público-urbano, servicios y pecuario. Los que tienen mayor incidencia sobre el río, tanto por su distribución como por los volúmenes de agua residual descargada, son los desarrollos industriales y los centros urbanos (Fig. 2, Cuadro V).
La agricultura de riego, como fuente de contaminación no puntual, ocupa una superficie pequeña en la subcuenca (Fig. 3). Se presenta en la parte norte, aguas abajo de la presa Atlangatepec (19º 33'33'' Latitud Norte y 98º 11' 10'' Longitud Oeste), con una superficie de 374.8 ha y en la parte sur, en la salida de la subcuenca, con una superficie de 9969.0 ha.
Condición biótica
Se identificaron trece familias de macroinvertebrados (Cuadro VI). Los sitios de referencia tuvieron cuatro familias con tolerancia baja a la contaminación del agua y tres familias con tolerancia moderada. Dos familias de tolerancia alta estuvieron presentes con muy pocos organismos. En los sitios perturbados se registró el mayor número de familias tolerantes a la contaminación. En los sitios S3 y S9 hubo presencia de oligoquetos, lo que coincide con lo reportado para zonas urbanas con contaminación por descargas de aguas residuales (Moreno y Callisto 2006, Correa-Araneda et al. 2010). En los sitios donde se combinan descargas industriales y urbanas (sitios S3, S5 y S9), sólo los organismos tolerantes están presentes, como la familia de Tubificidae y Glossiphoniidae. La presencia de organismos de la familia Tubificidade en los sitios S3 y S9 (el primero se ubica aguas abajo de la presa de Atlangatepec y el segundo en una zona de confluencia de agricultura de riego y descargas urbanas), indica altas concentraciones de materia orgánica, ya que organismos de esta familia se relacionan con ello (Martins et al. 2008). En los sitios S6, S7, S8, S10 y S11 la descarga industrial eliminó a los organismos totalmente, lo que apunta hacia un alto nivel de contaminación del agua, asociado con las descargas de aguas residuales (industrial, público-urbano, servicios y pecuario), la presencia de zonas agrícolas y la fuerte alteración del hábitat por modificación del cauce del río.
El análisis de correlación de Spearman mostró que el porcentaje de EPT, el porcentaje de organismos tolerantes, el porcentaje de fragmentadores y el número de familias EPT tuvieron el menor número de correlaciones altas (r > 0.8) respecto a otros índices bióticos. El análisis de relación de medias entre los sitios de referencia y los perturbados indicó que los índices bióticos sensibles a la perturbación son: el porcentaje del número total de familias, el porcentaje de colectores, el porcentaje de depredadores, el porcentaje de organismos tolerantes, el porcentaje de la familia dominante y el IBH. Las gráficas que relacionan índices bióticos con los parámetros de perturbación (conductividad, cubrimiento de cantos rodados por sedimentos finos y ancho de la zona ribereña), indicaron que el número total de familias, el porcentaje de organismos tolerantes, el IBH, el porcentaje de familia dominante y el porcentaje de colectores se relacionan con los parámetros analizados. Los índices bióticos relevantes para evaluar la condición biótica del río fueron el IBH, el porcentaje de organismos tolerantes, el porcentaje de colectores, el porcentaje de depredadores, el número total de familias y el porcentaje de la familia dominante. El índice indicó que sólo los sitios de referencia presentaron una condición óptima. Los sitios S3 y S4, en cambio, una condición moderadamente impactada (Cuadro VII).Estos sitios reciben un volumen menor de descargas de aguas residuales. Los sitios restantes mostraron una condición de hábitat pobre.
Condición de hábitat
El valor de los parámetros de condición de hábitat de los dos sitios de referencia fue mayor al de los sitios perturbados (Cuadro VIII). El análisis deconglomerados de los valores de condición de hábitat de los sitios, mostró una distinción clara entre los sitios de referencia y los sitios perturbados (Fig. 4). Dentro de los sitios perturbados, se observó la separación de tres grupos, lo que se asocia a la diferencia en el grado de perturbación del hábitat del río, causada por las actividades humanas.
EBI = estabilidad banco izquierdo, EBD = estabilidad banco derecho, AC = alteración del cauce, EC = estado del cauce, CCR = cubrimiento de cantos rodados por sedimentos finos, RVP = regímenes de velocidad y profundidad, SE/CP = sustrato para la epifauna/cobertura disponible, FR = frecuencia de rápidos, AZRBI = ancho de la zona ribereña banco izquierdo, AZRBD = ancho de la zona ribereña banco derecho, DS = depósito de sedimentos,PVBD = protección vegetal banco derecho, PVBI = protección vegetal banco izquierdo
El valor del índice de condición de hábitat indicó que a pesar de las diferencias en el grado de perturbación a lo largo del río, todos los sitios perturbados tienen una condición de hábitat pobre (Cuadro IX).
Condición de calidad del agua
Los valores obtenidos de los nueve parámetros físicos y químicos tienden a ser mayores en los sitios perturbados que en los de referencia (Cuadro X), excepto el oxígeno disuelto, que es mayor en los sitios de referencia. Los valores de pH y temperatura son similares a los obtenidos por García-Nieto et al. (2011) y Muñoz-Nava et al. (2012) en sitios cercanos a la evaluación, excepto en los valores de temperatura registrados en la parte más alta de la subcuenca. En el presente estudio se obtuvieron valores de temperatura más bajos. Esto se atribuye a que en el día de la medición de temperatura en la zona, se presentó un descenso de la misma.
VR = valor de referencia, (a) VR de la NOM-001-SEMARNAT-1996, (b) VR de la Norma Canadiense para la Protección de Vida Acuática (CCME 2007), (c) VR del estado de Iowa (IDNR 2009)
El índice de calidad del agua indicó que los sitios de referencia (S1 y S2) están en la categoría de excelente u óptima, asociada con una mínima perturbación por las actividades humanas (Cuadro XI). El sitio S4 presenta una condición favorable, debido a la presa Atlangatepec, ubicada aguas arriba del sitio, la cual recibe y retiene las descargas de aguas residuales. Los sitios S5 y S11 presentan una condición marginal, debido a un aumento de las descargas de las zonas urbanas e industriales sobre el río. Los sitios S3, S6, S7, S8, S9 y S10 presentan una condición pobre. En el sitio S3, ubicado en la entrada a la presa Atlangatepec, se concentran las aguas residuales del municipio de Tlaxco y su zona industrial. La condición de los otros sitios refleja la influencia permanente de grandes volúmenes de descargas de agua residual industrial y público-urbano.
Asociación de los índices bióticos de condición biótica, condición de hábitat y calidad del agua
El análisis de correlación de Spearman mostró que el número de familias del orden Trichoptera, el porcentaje de EPT, el porcentaje de Ephemeroptera, el porcentaje de raspadores y el porcentaje de fragmentadores se relacionan positivamente con la puntuación total de hábitat (Fig. 5). El valor de estos índices bióticos disminuye si el hábitat se encuentra perturbado (Barbour et al. 1999, Tetra Tech 2000). Estos resultados son similares a los obtenidos por Ríos y Bailey (2006), Sponseller et al. (2001) y Boosoong et al. (2009), quienes encontraron para sitios forestales y una mínima perturbación, valores más altos de porcentaje de EPT, porcentaje de Ephemeroptera y número de familias del orden Trichoptera.
Los parámetros de calidad del agua que mostraron mayor asociación con los índices de condición biótica del río fueron la demanda biológica de oxígeno para un periodo de cinco días (DBO5), el fósforo total y la temperatura (8, 7 y 6 correlaciones significativas, respectivamente). Las correlaciones significativas fueron con el número total de familias, número de familias de EPT, número de familias del orden Trichoptera, porcentaje de Ephemeroptera y porcentaje de raspadores. Todos estos índices bióticos son indicadoras de buena calidad del sitio.
En términos generales, se confirmó que los índices bióticos de condición biótica que se reportan como indicadores de una mejor calidad del agua, se asociaron en mayor proporción con los parámetros físicos y químicos de calidad del agua. Estos resultados indican que tanto la condición del hábitat como los diferentes tipos de descargas de aguas residuales inciden significativamente en la condición biótica del río Zahuapan.
Condición ecológica global del río Zahuapan
La evaluación global de la condición del río indicó que todos los sitios perturbados tienen un estado pobre (la condición más baja o más impactada), a excepción del sitio S4 (Cuadro XII), localizado aguas abajo de la presa Atlangatepec. Esto se debe a que no recibe grandes volúmenes de aguas residuales y a que la presa retiene contaminantes de las descargas vertidas aguas arriba. Se considera que, en términos generales, hubo concordancia en las tres evaluaciones de la condición del río. Aunque la evaluación de calidad de agua (parámetros físicos y químicos), mostró una ligera diferenciación entre los sitios, seguida de la evaluación biótica.
CONCLUSIÓN
El enfoque multimétrico, basado en el protocolo de evaluación biótica en ecosistemas acuáticos, tomando como bioindicadores a los macroinvertebrados, permitió evaluar la condición ecológica de los sitios donde hay un marcado impacto de las actividades humanas. Se logró diferenciar claramente los sitios tomados como referencia, con mínima perturbación de los sitios perturbados por actividades antrópicas. Comparadas con la evaluación global (que integra la condición de hábitat, biótica y de calidad de agua), las evaluaciones de hábitat y biótica asignaron un grado de condición diferente a un sitio, mientras que la evaluación biótica difirió con la evaluación global en dos sitios. En general, todas las evaluaciones indican que las actividades pecuarias, urbanas e industriales tienen un impacto negativo en la condición ecológica del Río Zahuapan. La condición ecológica global de este río, a excepción de las zonas forestales donde la perturbación por las actividades humanas es mínima, se ubica en la categoría de pobre, que corresponde a la peor condición ecológica de los métodos utilizados en esta evaluación.
Los resultados mostraron que las variables de hábitat, bióticas y de calidad del agua están relacionadas entre sí, por lo que la restauración de la condición ecológica del río debe considerar a las tres variables. Esto requiere la adopción de un enfoque de manejo integral de cuencas. Asimismo, las actividades pecuarias y las descargas de aguas residuales urbanas e industriales son los principales factores de deterioro de la condición del río Zahuapan. La influencia de las actividades agrícolas es más localizada y de menor impacto que las derivadas de los asentamientos humanos y desarrollos industriales.