Repositorio Institucional de Documentos

Trabajo de Fin de Grado de Ciencias Ambientales Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca. . Autora: Marta Michavila López Directora: Rocío López Flores Co – director: José Antonio Cuchi Oterino Escuela Politécnica Superior de Huesca. Universidad de Zaragoza. 8/01/2014 Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 2 Índice Resumen ........................................................................................................................................ 4 Abstract ......................................................................................................................................... 4 1. Introducción, el estado ecológico de las masas de agua .......................................................... 6 1.1 Los ríos como sistemas dinámicos y la importancia de la conectividad de hábitat ............ 6 1.2 Calidad del agua y bioindicadores ....................................................................................... 7 1.3 Acción antrópica: las infraestructuras hidráulicas .............................................................. 9 2. Marco de estudio: el río Gállego ............................................................................................. 11 2.1 El entorno del Gállego a su paso por la provincia de Huesca ........................................... 11 2.2 Impacto humano ............................................................................................................... 15 2.3 Antecedentes de estudios ................................................................................................. 16 3. Metodología ............................................................................................................................ 18 3.1 Área de muestreo .............................................................................................................. 18 3.1.1 Ubicación de los puntos ............................................................................................. 18 3.1.2 Principales impactos que afectan a los puntos .......................................................... 20 3.2 Toma y tratamiento de muestras de muestras ................................................................ 27 3.2.1 Diseño del muestreo .................................................................................................. 27 3.2.2 Análisis físico – químico.............................................................................................. 28 3.2.3 Muestreo y determinación de perifiton .................................................................... 28 3.2.4 Muestreo y determinación de Macroinvertebrados .................................................. 29 3.3 Análisis estadístico ............................................................................................................ 32 3. Resultados ............................................................................................................................... 34 4.1 Evolución espacio – temporal de los parámetros físico - químicos .................................. 34 4.2 Evolución espacio – temporal del estado trófico .............................................................. 38 4.3 Estado ecológico del río Gállego: diversidad taxonómica, IBMWP y riqueza ................... 40 4.3.1. Evolución espacio – temporal de los índices de calidad y diversidad. ...................... 40 4.3.2 Índice IBMWP relativo en función a los ecotipos ....................................................... 42 Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 3 4.4 Relación entre índices biológicos. ..................................................................................... 43 4.5 Estructura de la comunidad de macroinvertebrados en el río Gállego ............................ 45 5. Discusión ................................................................................................................................. 50 5.1 Diagnosis del estado ecológico del río Gállego ................................................................. 50 5.1.1 Parámetros físico – químicos ..................................................................................... 50 5.1.2 Estado trófico ............................................................................................................. 51 5.1.3 Evolución de los índices a lo largo del año ................................................................. 51 5.1.4 Comparación de índices de calidad ............................................................................ 55 5.2 Estructura funcional de la comunidad de macroinvertebrados del río Gállego ............... 56 6. Conclusiones............................................................................................................................ 59 7. Agradecimientos ..................................................................................................................... 61 8. Bibliografía .............................................................................................................................. 62 9. Anexo ....................................................................................................................................... 71 Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 4 Resumen Son evidentes los impactos derivados de las infraestructuras hidráulicas que sufren en la actualidad los ecosistemas fluviales. El río Gállego abarca una amplia diversidad y complejidad de ecosistemas a su paso, debido a las diferencias climáticas en su curso, pero se encuentra regulado desde la cabecera. En este trabajo se abordan los efectos que ocasionan las cadenas de embalses y las presiones antrópicas sobre la calidad ecológica y el estado trófico de un río altamente alterado, el Gállego. Se realizaron muestreos estacionales a lo largo de un año en seis puntos situados aguas abajo de los diferentes embalses existentes. Se efectuaron en cada uno de ellos medidas de parámetros físico – químicos y del estado trófico mediante la concentración de clorofila a presente en el perifiton. Además, se determinó para la comunidad de macroinvertebrados el índice IBMWP, la diversidad de Shannon y la riqueza taxonómica. Los macroinvertebrados han sido utilizados en los últimos años como uno de los principales bioindicadores de la calidad del agua, ya que poseen ciclos de vida lo suficientemente extendidos como para detectar alteraciones producidas tiempo atrás. Todo ello se analizó mediante ANOVAS, utilizando además un análisis regresión lineal de Pearson para establecer la relación existente entre índices biológicos. Además se determinó la estructura funcional de dicha comunidad de acuerdo a los rasgos biológicos: estrategia trófica, grado de eutrofia, tolerancia a la contaminación y microhábitat. Se utilizó un FCA (Análisis de correspondencia difuso) y un análisis de co – inercia con el fin de conocer la distribución de cada familia y rasgo biológico en el curso fluvial y las relaciones entre los tramos. Los resultados mostraron que las cadenas de embalses generan una eutrofización aguas abajo de las mismas, conllevando un aumento de la biodiversidad y la disponibilidad de alimento para los organismos. Además se observó cómo el río Gállego está siendo objeto de un grave deterioro provocado por los vertidos de aguas residuales urbanas e industriales sin tratamiento previo, siendo más intensos en la zona de cabecera. Se pudo comprobar como a lo largo de este curso fluvial, las especies generalistas fueron capaces de colonizar rápidamente espacios perturbados, desplazando de este modo a especies más sensibles e indicadoras de buena calidad. Abstract Impacts of hydraulic infrastructures currently suffering river ecosystems are evident. The Río Gállego covers a wide variety and complexity of ecosystems along the whole river, due to climatic differences in its course, but is regulated from the header. This paper addresses the effects which are caused by chains of dams and anthropogenic pressures on the ecological quality and trophic status of a highly altered River, the Gállego. Sampling was carried out seasonal along one year in six-point located downstream from the different existing reservoirs. They were made in each of these, measurements of physical parameters – chemical and of trophic status through the concentration of chlorophyll a present in the periphyton. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 5 In addition, it was established for the macroinvertebrate community the IBMWP index, the Shannon diversity and taxonomic richness. Macroinvertebrates have been used in recent years as one of the main bio-indicators of water quality, since they have life cycles sufficiently widespread as to detect alterations produced long ago. Everything was analyzed by ANOVAS, using the linear regression of Pearson to establish the relationship between biological index. In addition, the functional structure of the community according to the biological traits was determined: food strategy, degree of eutrophication, tolerance to pollution and microhabitat. We used a FCA (diffuse correspondence analysis) and an analysis of co - inertia in order to know the distribution of each family and biological trait in the river course and relations between the sections. The results showed that the chains of dams generate an eutrophication downstream from them, leads to an increase in biodiversity and the availability of food for the organisms. Moreover, it was observed how the Gállego River is under a serious deterioration caused by discharges of urban and industrial wastewater without treatment, being more intense in the header area. It was verified that along this river course, generalist species were able to quickly colonize disturbed areas, thus displacing more sensitive and indicator species of good quality. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 6 1. Introducción, el estado ecológico de las masas de agua 1.1 Los ríos como sistemas dinámicos y la importancia de la conectividad de hábitat Los ríos presentan una conectividad entre sus tres dimensiones que da cabida a diferentes espacios con alta diversidad de hábitat. Tal hecho los hace idóneos para que diferentes especies encuentren las condiciones apropiadas para su desarrollo (Poff et al., 1997; Batalla et al., 2004). A través de la longitud del canal principal se observan gradientes de sus factores hidrológicos y geomorfológicos y variaciones en las comunidades biológicas. Según la geomorfología de la cuenca, el cauce se ensancha y se contrae surgiendo así una dimensión transversal o lateral referente a la anchura del río y su llanura de inundación. La vertical refleja la conexión existente con las zonas sobre las que discurren las aguas como el medio hiporreico y los acuíferos locales o regionales (Ward & Wiens, 2001; González & García, 2007). Todo ello estando relacionado con una cuarta dimensión según Ward, 1989, la temporal (Figura 1), el hilo conductor de las tres anteriores (Ward, 1989). Figura 1. Conceptualización de las cuatro dimensiones naturales que caracterizan la conectividad de un sistema fluvial lótico: lateral, longitudinal, vertical y la temporal (Ward, 1989). El mantenimiento de los patrones naturales de los ecosistemas fluviales en sus dimensiones resulta esencial para la viabilidad de muchas especies. Muchos organismos acuáticos deben su supervivencia a la capacidad de moverse libremente a través de la corriente. La biodiversidad y heterogeneidad característica de estos hábitats dependerá de un modo directo de la conectividad entre cauce, llanura de inundación y el lecho del río (Benn, 2002). A la hora de tratar un sistema fluvial a escala de cuenca, la dimensión longitudinal adquiere mayor importancia. Representa el eje a través del cual se mantiene la continuidad de flujos, se modifica la energía desde las partes altas a las bajas y se organizan las comunidades en función a diferentes factores (FISRGW, 1998; González & García, 2007). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 7 Surge así el concepto de “río continuo” según el cual todos los cambios ocurridos en un tramo fluvial determinado repercutirán aguas abajo sobre los distintos equilibrios que caracterizan al ecosistema fluvial (Vannote et al., 1980; Minshall et al., 1983; Ward, 1989). Uno de los factores determinantes es el régimen natural de caudales, que controla diferentes aspectos físicos y ecológicos del río, como el transporte de sedimentos y los ciclos de nutrientes que afectan directamente a las comunidades (Poff et al., 1997; Batalla et al., 2004). En los ecosistemas fluviales, las perturbaciones naturales, ya sean a pequeña o gran escala (crecidas, grandes avenidas, sequías, renovación de sedimentos) son muy necesarios (Dale & Turner, 2008) ya que provocan que se genere una amplia escala de procesos y modelos responsables de altos niveles de biodiversidad (Ward, 1998). 1.2 Calidad del agua y bioindicadores Un factor muy importante que determina el estado ecológico de un ecosistema fluvial es la calidad biológica del agua. Este término hace referencia al hecho de poseer unas características naturales apropiadas para el desarrollo de organismos que le son propios en función a su hábitat (Alba – Tercedor, 1996). La cuantificación de los parámetros físico – químicos resulta muy importante al influir significativamente en la abundancia y distribución de los organismos acuáticos (Allan, 1995; Wetzel, 2001; Camargo et al., 2011). Es un método tradicional que proporciona información del estado puntual de la masa de agua en el momento de la toma de la muestra. Los sólidos totales, el pH, la temperatura, la conductividad, la dureza, el color, el olor, el oxígeno disuelto, alcalinidad y los nutrientes son algunos ejemplos de estos parámetros (Jiménez, 2000). Sin embargo, se debe tener en cuenta que las corrientes de agua se encuentran en continua dinámica y es difícil estudiar la situación únicamente mediante estos parámetros, por lo que es necesario analizar la calidad desde otro punto de vista más amplio (Alba – Tercedor, 1996). Los bioindicadores son el complemento ideal para lograr un estudio integral del medio. Se trata de una técnica ecológica que se sustenta en la medición de la diversidad y la ausencia o presencia de determinados organismos o conjuntos de ellos. Estos indican la situación actual de un ecosistema y reflejan situaciones que pudieron darse tiempo atrás. Es decir, “contienen información” de su respuesta ante determinados factores del medio (Li et al, 2010). Para su empleo se requiere un conocimiento previo de sus ciclos de vida, variabilidad estacional y variaciones naturales de modo que sea posible comparar sus condiciones antes y después de una perturbación. Los macroinvertebrados (Figura 2, A) han sido utilizados en los últimos años como uno de los principales bioindicadores de la calidad del agua. Las razones fundamentales radican en su tamaño relativamente grande, su facilidad de muestreo y sobre todo el poseer ciclos de desarrollo lo suficientemente extendidos como para detectar alteraciones producidas tiempo atrás. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 8 Tras una gran perturbación puntual necesitan un plazo de alrededor de un mes para volver a recolonizar el lugar afectado (Alba – Tercedor, 1996). Muestran una respuesta inmediata ante un determinado impacto, dando lugar a un patrón de estímulo – respuesta ante alteraciones físico – químicas (Rosenberg & Resh, 1986; Figueroa et al., 2003). Su diversidad presenta una amplia gama de tolerancia frente a diferentes parámetros de contaminación (Hellawell, 1986; Alba – Tercedor, 1996). Otro aspecto añadido es el económico, pues atendiendo a la relación coste/beneficio del procesado e identificación de muestras resultan muy ventajosos (Ohio EPA, 1987; Purschke, 1993; Alba – Tercedor, 1994; Alba – Tercedor, 1996). Mediante el estudio de la diversidad taxonómica y la riqueza de especies se puede obtener una visión general de la calidad ecológica del ecosistema. Sin embargo es más precisa la utilización del índice IBMWP al valorar y dar mayor puntuación a especies indicadoras de buena calidad (Alba – Tercedor et al., 2002) ya que con la diversidad únicamente se tendría en cuenta el número de individuos y no sus condiciones biológicas. Figura 2. (A) Fotografía de un macroinvertebrado (Familia Hydropsychidae) (CHE, 2009). (B) Detalle de roca cubierta de perifiton del río Gállego del tramo de Formigal en Invierno. Otro bioindicador muy importante es el perifiton. Se trata de la mezcla de algas, bacterias y hongos que crece en las rocas, troncos, macrófitos y estructuras del cauce, que se encuentras unidas al sustrato por una matriz de polisacáridos, que evita el arrastre por la corriente (Figura 2, B). Es una fuente de alimento importante de la cadena trófica del ecosistema fluvial (Rutherford & Cuddy ,2005; Kurthesi et al., 2008). La concentración de clorofila que contenga es proporcional a la biomasa algal y será determinante a la hora de diagnosticar el estado trófico del río (Dodds, 1998). La Directiva Marco del Agua (DMA) dicta que los países miembros deben determinar el estado ecológico de sus masas de agua continentales basándose en criterios biológicos, estableciendo caudales ecológicos y utilizando entre otras técnicas los bioindicadores (D.O.C.E., 2000; Cuellar et al., 2002). Uno de los objetivos fundamentales de esta directiva reside en conseguir masas de agua con buena calidad y buen estado ecológico. B Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 9 1.3 Acción antrópica: las infraestructuras hidráulicas Los ecosistemas fluviales vienen siendo objeto de diversas actividades antrópicas, provocando numerosos impactos y perturbaciones de su equilibrio natural, directas o indirectas. La combinación de un gran aumento de la población unido a un auge de la urbanización y de la industrialización, ha generado un aumento de diferentes presiones sobre los recursos hídricos. A los problemas de contaminación, ocupación y alteración de orillas se les une otro no menos importante, la regulación de los ríos mediante infraestructuras hidráulicas (Sloto, 1987; Weaver & Garman, 1994; Oscoz et al., 2006). Se considera la alteración del régimen natural de caudales como una de las amenazas más serias para la sostenibilidad ecológica de los ríos y sus llanuras de inundación (Naiman et al., 1995; Sparks, 1995; Lundqvist, 1998; Ward et al., 1999; Bunn, 2002). Los embalses alteran dicho régimen, generando cambios en la frecuencia y magnitud de las inundaciones y una disminución del caudal natural (Ward & Stanford, 1979; Petts, 1984; Ligon et al., 1995; Ward & Stanford, 1995; Kondolf, 1997; Batalla et al., 2004). Acompañando a los embalses se encuentran otras infraestructuras como pequeños azudes, canalizaciones y centrales hidroeléctricas. Éstas generan cambios súbitos de caudal que provocan daños letales en muchos organismos al ser arrastrados durante los niveles máximos o por lo contrario, aislados en seco durante los niveles más bajos (Cushman, 1985; Díaz -Hernández, 2005). Se dan de este modo perturbaciones en el gradiente de temperatura, en los hábitats, y en la disponibilidad de recursos tróficos, entre otros. Los organismos deben hacer frente a nuevas condiciones que les son impuestas (Ward, 1976; Armitage, 1984; Oscoz et al., 2006.). En un entorno donde aumentan las zonas lénticas frente a las lóticas, los ciclos de vida se ven directamente alterados (Scheidegger & Bain, 1995; Díaz -Hernández, 2005). La ausencia de crecidas periódicas genera un ciclo cerrado de condiciones de contaminación – eutrofización, al propiciarse entre otros, un desarrollo excesivo de algas (Petts, 1984; Díaz -Hernández, 2005) creándose así un escenario de degradación casi permanente. El sustrato donde habitan queda modificado casi en su totalidad (Reiser et al., 1989; Díaz - Hernández, 2005). Tras un embalse, las especies autóctonas y las más sensibles (Alba – Tercedor, 1996) se extinguen o se desplazan mientras que las generalistas (muchas alóctonas) y más tolerantes predominan (Aguirre & Bikuña, 2003), disminuyendo la diversidad de la comunidad de organismos. Sin embargo no ocurre lo mismo cuando se suceden en un curso fluvial varios embalses. En este caso, los tramos fluviales río abajo irán sumando las especies colonizadoras de los tramos superiores, lo que tendrá como consecuencia el aumento de la diversidad (Margalef, 1983). No obstante, este cambio no se verá reflejado en un aumento del índice de calidad ecológica IBMWP ya que tiene en cuenta los requerimientos ecológicos de cada familia otorgándoles un valor mayor a los taxones indicadores de calidad de agua más alta (CHE, 2009). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 16 Estos embalses de cabecera han introducido importantes modificaciones en volumen, régimen e irregularidad de caudales tanto sólidos como líquidos. Hasta su paso por Sabiñánigo el río Gállego no vuelve a ser represado, en este caso por una pequeña presa. Cerca de los mayos de Riglos se haya el embalse de La Peña, el embalse más antiguo de su curso en funcionamiento desde 1913, con 15 hm3. Cierra esta cadena de embalses el pantano de Ardisa de 1931 y 3hm3. Actualmente la capacidad reguladora de tanto de Lanuza como de La Peña, apenas modifica el régimen fluvial. (García & Moreno, 2003). En lo referente a contaminación ya desde su nacimiento se ve afectado por los vertidos del Portalet y la estación de esquí de Formigal, seguido de pequeños núcleos pirenaicos que carecen de depuración de sus aguas residuales. Cabe destacar como gran foco de contaminación el núcleo urbano de Sabiñánigo debido a su importante desarrollo industrial acontecido a principios y mediados del siglo XX. Durante este período se llevaron a cabo malas prácticas ambientales como vertidos de residuos generados por su población e industrias asociadas. La mayor productora de desechos químicos fue la empresa Inquinosa que fundada en 1970, produjo lindano de 1975 – 1992, un producto altamente nocivo (Fernández et al., 2013). Se trata de un plaguicida organoclorado halogenado que tiende a bioacumularse en tejidos, plantas y animales (Prakash et al., 2004; Robles – González et al., 2006). Se caracteriza por una alta persistencia en el suelo, difícil hidrolización y lenta biodegradación (Hernández et al., 1991). Los desechos de Inquinosa no seguían ninguna operación ni eliminación en particular, por lo que eran dispuestos en rellenos sanitarios sin garantía generándose grandes episodios de contaminación incrementados por la lixiviación e infiltración en los suelos derivada de los residuos. En la actualidad existen cuatro lugares, todos ellos ubicados cerca del cauce fluvial, donde todavía persisten señales evidentes de contaminación: la antigua fábrica de Inquinosa, los vertederos de Sardás y Bailín y la presa de Sabiñánigo. Aguas superficiales y subterráneas (el importante acuífero del aluvial del Gállego) han estado y continúan siendo directamente afectadas (Fernández et al., 2013). 2.3 Antecedentes de estudios Son escasos los estudios y la información publicada en referencia al río Gállego. Los existentes tratan temas muy concretos y en zonas específicas. El organismo de la Confederación Hidrográfica del Ebro, del Ministerio de Agricultura, Alimentación y Medio Ambiente, se encarga de realizar estudios y análisis en continuo. Miden datos y parámetros a tiempo real (pluviometría, temperatura, caudal, presión, etc.) en las estaciones de aforo y meteorológicas distribuidas por la cuenca del Ebro mediante el Servicio Automático de Información Hidrológica (SAIH). El proyecto SAICA o Sistema Automático de Información de Calidad de las Aguas toma datos referentes a parámetros – físicos en estaciones situadas bajo focos importantes de vertido. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 17 En el río Gállego hay dos situadas en Sabiñánigo (al consultarla los datos salen no válidos) y en Villanueva de Gállego. Es posible el acceso a una fuente cartográfica de la cuenca del Ebro donde se encuentran capas referentes a múltiples temas así como informes de bioindicadores como son los macroinvertebrados. Hay que destacar que son pocas las estaciones situadas en la cuenca del Gállego y algunas se encuentran fuera de servicio por lo que la información resulta difícil de homogeneizar. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 18 3. Metodología 3.1 Área de muestreo 3.1.1 Ubicación de los puntos Se eligieron seis puntos a lo largo del curso del Río Gállego (Figura 8), a su paso por la provincia de Huesca en función a diferentes criterios. El principal criterio fue la existencia de un embalse aguas arriba de cada uno de ellos (Figura 7) al tratarse del impacto más relevante valorado en este trabajo. Se tuvo en cuenta su accesibilidad y se ha evitó situarlos inmediatamente aguas debajo de perturbaciones hidromorfológicas o fisicoquímicas (CHE, 2009). Cada punto se subdivide a su vez en otros tres subpuntos a lo largo de un transecto de aproximadamente 50 metros de longitud, de cara a que aparezcan los suficientes microhábitats (almohadillas de musgos, piedras y cantos gruesos, macrófitos) en donde poder realizar un muestreo representativo del tramo (Pujante et al., 1995). Figura 7. Localización del área de estudio en la Cuenca del río Ebro, España. Situación de los puntos de muestreo y de los embalses de la Cuenca del río Gállego (Elaboración propia, GIS). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 19 En la Tabla 2 se muestra la ubicación (coordenadas), altura sobre el nivel del mar y término municipal de los seis puntos de muestreo. Punto Nombre Término municipal Altura (msnm) UTM X UTM Y G1 Formigal Sallent de Gállego 1781 711771,731 4741329,186 G2 El Pueyo Panticosa 1140 720691,4002 4733452,772 G3 Búbal Biescas 793 719429,5693 4726948,396 G4 Sabiñánigo Sabiñánigo 719 715857,8898 4707492,884 G5 La Peña Santa Eulalia 417 685421,2553 4682293,829 G6 Ardisa Gurrea de Gállego 309 684982,7612 4654022,837 Tabla 2: Coordenadas, altitud y término municipal de las 6 estaciones de muestreo. Figura 8. Detalle de cada uno de los seis puntos de muestreo. De arriba abajo y de izquierda a derecha: Formigal (G1) y El Pueyo (G2) en Invierno, Búbal (G3) en primavera, Sabiñánigo (G4), Santa Eulalia (G5) y Gurrea de Gállego (G6) en Otoño. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 20 3.1.2 Principales impactos que afectan a los puntos Se realizó un análisis de los impactos que afectan a cada punto de muestreo, obteniendo los datos de los registros del geoportal SITEbro, perteneciente a la Confederación Hidrográfica del Ebro. Se utilizaron los referentes a vertidos de aguas residuales urbanas e industriales, captaciones, abastecimientos y puntos hidroeléctricos inscritos. Todo ello junto con municipios cercanos, embalses y red hidrográfica, se ubicó en mapas GIS con base de ortofotos del Sistema de Información Territorial de Aragón. Punto 1, Formigal. A pocos kilómetros del nacimiento del río Gállego, se encuentra el primer punto de muestreo (Figura 9), a 1781 metros de altitud. Se trata de una zona clasificada como río de alta montaña dentro del curso alto del río (Ollero et al, 2004). Confluyen en esta zona los vertidos con carácter de agua residual urbana procedentes de la estación de esquí de Formigal, así como de las Ventas del Portalet. Se ve afectado por las diferentes obras realizadas tanto de tomas y captaciones de aguas, puentes, instalaciones de las pistas de esquí y edificios. Se trata de un punto con dificultad de acceso en la estación de invierno, dada la cantidad de precipitación en forma de nieve de esta época y su acumulación en los márgenes del cauce por la acción de las máquinas quitanieves en la carretera situada a poca distancia del mismo. Figura 9. Ubicación del punto 1 (Formigal) y representación de los principales impactos que le afectan. Formigal Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 21 Punto 2, El Pueyo. En el curso alto del río Gállego, en zona de río de alta montaña y a 1140 metros de altura se encuentra el segundo punto (Figura 10), que debido al difícil acceso se sitúa cerca de la cola del embalse de Búbal, tras pasar la central hidroeléctrica del embalse de Lanuza y a poca distancia de la población de El Pueyo de Jaca. Se ve afectado en primer lugar por los vertidos de aguas residuales urbanas tanto de poblaciones por las que transcurre el cauce (Formigal, Sallent de Gállego), como de diferentes actividades como embarcaderos o restaurantes. Importantes también los vertidos de agua residual industrial que recibe de una industria de comercio al por menor de combustible. Recibe a los afluentes Aguas Limpias (en la cola del embalse de Lanuza) y al río Escarra. Las aguas de este último van a parar al embalse directamente y recoge los vertidos e impactos de la población y embalse de Escarra. El Barranco de Portet desemboca en el mismo embalse de Lanuza y cuenta con la central hidroeléctrica de Portet y un pequeño azud, que abastece a la población de San Juan de Lanuza. Atraviesa su paso dos Lugares de Importancia Comunitaria (LIC): Monte Pacino y Foz Escarrilla – Cucuraza. Se suceden en su curso diversas captaciones y tomas de agua destinadas principalmente al abastecimiento de poblaciones. Se detectó una tubería pocos metros aguas abajo del punto de muestreo (Figura 11), que quedaba al descubierto en las estaciones de invierno y primavera. Figura 10. Ubicación del punto 2 (El Pueyo), representación de los principales impactos que le afectan. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 22 Figura 11. Detalle de la tubería al descubierto en El Pueyo durante el muestreo de invierno. Se observó un desarrollo excesivo de algas en torno a su ubicación, así como espumas procedentes del vertido. En el muestreo de verano el punto se encontraba completamente cubierto por las aguas del pantano de Búbal (Figura 12), convirtiéndose así en un sistema léntico. Figura 12. Detalles del punto de El Pueyo en Verano, totalmente anegado por las aguas del pantano de Búbal. Punto 3, Búbal. A 793 metros de altura es el último punto dentro del curso alto del río Gállego y en zona de río de alta montaña. Está situado a aproximadamente 2 kilómetros aguas abajo del embalse de Búbal (Figura 13), antes de llegar a la ermita de Santa Elena. Se encuentra dentro del Lugar de Interés Comunitario (LIC) Tendeñera y de la Zona de Especial Protección para las Aves (ZEPA) Viñamala. Desde la alteración morfológica de zona de la presa debida a la regulación de caudales, hasta la zona de estudio, no existen impactos antrópicos sobre el río. Previamente el Gállego recibe varios afluentes: El Caldarés (y su afluente El Bolatica) y el Aguilero que desembocan en el mismo embalse. Se ve afectado por los vertidos de aguas residuales urbana de las poblaciones que atraviesa (Piedrafita de Jaca, Tramacastilla, Búbal y Hoz de Jaca) así como vertidos procedentes de instalaciones hosteleras. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 23 Figura 13. Ubicación del punto 3 (Búbal) y representación de los principales impactos que le afectan. Punto 4, Sabiñánigo. En el curso medio del río Gállego, en zona de río de montaña húmeda – calcárea y a 719 metros de altitud se encuentra el cuarto punto (Figura 14), aguas abajo de la depuradora de Sabiñánigo. Confluyen principalmente entre el punto tres y cuatro, el río Sia (previo a la población de Orós Alto), el río Oliván (tras Senegüe), el río Aurin (en la cola del embalse de Sabiñánigo) y el río Basa (tras el embalse de Sabiñánigo). A su vez se dan captaciones para abastecimiento de poblaciones como las de Orós Alto y Bajo y para comunidades de regantes (San Pedro del río Gállego, acequia del Molino de Senegüe – Sorripas, etc.). Antes de llegar hasta el lugar de muestreo el río Gállego sufre múltiples afecciones. A su paso se encuentra con la central hidroeléctrica de la presa de Santa Elena y con varios azudes: De Biescas I y II (con sus centrales hidroeléctricas), El Molino de Biescas, y el azud subterráneo de los Molinos de Sabiñánigo. Recoge diversos vertidos de aguas residuales urbanas de poblaciones (Betes de Sobremonte, Aso de Sobremonte, Biescas, Gavín, Escuer, Orós Alto y Bajo, Larrede, Senegüe, Sorripas, Javierre del Arzobispo, Satué, Alaurin, Cartipona) y de diferentes actividades (Piscifactoría, viviendas aisladas, camping). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 24 Sin embargo, son los vertidos de origen urbano e industrial tras pasar el embalse y la población de Sabiñánigo los que le afectan de un modo más directo (Industria de transformación de aluminio, estación de transferencia de residuos sólidos urbanos valorizables, lavado de áridos y hormigonado). El embalse de Sabiñánigo genera impactos directos sobre este punto debido a la regulación de caudales, a las fuentes puntuales (autorización de sustancias peligrosas) y difusas de contaminación (agricultura, pastos intensivos). A todo ello se le une la cercanía de los vertederos de Bailín y Sardás, focos de una gran contaminación por lindano. Estudios del 2009 (Dirección General de Calidad Ambiental, Gobierno de Aragón 2010; Fernández et al., 2013) demostraron que todavía siguen presentes altos niveles de residuos químicos afectando directamente a aguas tanto subterráneas como superficiales (Fernández et al., 2013), en este caso el curso del río Gállego y el tramo muestreado en concreto. En el estudio también se detectó mercurio procedente de los vertidos de la industria ERCROS años atrás. En la estación de verano el agua del río tenía un color marrón – grisáceo (Figura 15) y se podía apreciar un característico. Figura 14. Ubicación del punto 4 (Sabiñánigo) y representación de los principales impactos que le afectan. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 25 Figura 15. Detalle del color del agua en el tramo de Sabiñánigo en la estación de verano. Punto 5, Santa Eulalia. En un tramo de río clasificado como eje mediterráneo continental poco mineralizado en el curso bajo del río Gállego y a 417 metros de altura se encuentra el quinto punto (Figura 16), a 15 kilómetros del embalse de la Peña. Pocos vertidos afectan a este punto, únicamente de aguas residuales urbanas de camping, actividad hostelera pero no se registra actividad industrial importante. Uno de los principales impactos generados por el embalse de la Peña se trata de las fuentes puntuales de contaminación como son los vertidos urbanos no saneados y vertidos biodegradables urbanos. A este aspecto se le suma la regulación de caudales de la presa. Se le unen en este recorrido los principales afluentes: el río Abena (tras el punto cuatro), el río Guarga, el río Val de San Vicente, el río Garona, el río Triste y el Asabón (en la cola del embalse de la Peña) y los barrancos del río Moro y de Artaso y San Julián. Atraviesa el Lugar de Interés Comunitario (LIC) Sierras de Santo Domingo y Caballera y la Zona de Especial Protección para las Aves (ZEPA) Sierra de Santo Domingo y Cabellera y Río Onsella. Se dan diversas captaciones de agua para riego y abastecimiento de poblaciones durante su curso. Se encuentra con varios azudes (Del Molino de Barangua, Molino de Yeste, varios innominados, etc.) y diferentes centrales hidroeléctricas como la de Carcavilla en la confluencia con el río Garona y la de Javierrelatre con el río Guarga. Este punto no sufre afecciones antrópicas directas. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 32 3.3 Análisis estadístico Se utilizó análisis de varianza (ANOVA) para determinar las diferencias de las variables continuas (pH, conductividad, clorofila, oxígeno disuelto, temperatura, IBMWP, diversidad taxonómica, riqueza) respecto a las diversas categorías de los factores (estación y punto de muestreo). - En primer lugar se comprobó la normalidad mediante el test de Shapiro - Wilks. Si las variables seguían distribución normal se calculó mediante el test de Barlett la homogeneidad de varianzas. En el caso de no haber normalidad y/o homogeneidad, se probó diversas transformaciones de las variables, empezando por el logaritmo en base 10 y la raíz cuadrada. - Para variables normales y con homogeneidad de varianzas, se realizó un Análisis de Varianza (ANOVA). Se utilizó el test no paramétrico de Kruskal - Wallis cuando la transformación de las variables no mostró normalidad y/o homogeneidad de varianzas. - Se realizó un test post – hoc tras la ANOVA para determinar las diferentes agrupaciones de las diferentes categorías de cada factor. Para determinar la relación entre los índices biológicos se llevaron a cabo análisis de correlación lineal de Pearson. Se tuvo en cuenta el índice IBMWP, la diversidad de Shannon, la riqueza taxonómica y la clorofila a. Para analizar la estructura de la comunidad de macroinvertebrados en el curso del río Gállego se utilizó un FCA o Análisis de Correspondencia Fuzzy (difuso) además de un test co - inercia. Este tipo de análisis multivariable surge para relacionar conjuntos de datos que son difusos por naturaleza, es decir en los que se encuentra información sobre especies (familias en este caso), que asumen varias categorías de una misma variable con diferente grado afinidad. Permite analizar todos los rasgos biológicos tenidos en cuenta al mismo tiempo y su relación con las familias presentes en los puntos de muestreo (Usegglio et al., 2000). Se consideraron los traits o rasgos: tolerancia a la contaminación, grado de eutrofia, alimentación y microhábitat. - Para comenzar se determinó la afinidad de cada familia presente en la totalidad del muestreo (Tablas 7, 8 y 9, Anexo) con cada categoría perteneciente a cada uno de los rasgos (0=Ninguna 1= Poca 2=Buena 3 o >3=Mucha afinidad). También se tuvo en cuenta la abundancia de cada familia en cada punto de muestreo y estación. - Se realizó un FCA o análisis de correspondencia difuso para cada rasgo biológico o trait analizado. Posteriormente mediante un análisis de componentes principales difuso donde se relacionaron entre sí los traits y las familias, seguido de un análisis Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 33 de componentes principales donde se ubicaron las familias y los puntos de muestreo. - Por último se efectuó un análisis de co –inercia con el fin de relacionar los rasgos biológicos, las familias y los puntos de muestreo en conjunto y concluir la estructura de la comunidad de macroinvertebrados en el río Gállego. Para realizar los análisis estadísticos se utilizó el programa R, las ANOVAS mediante el paquete R Commander (R Development Core Team, Open Source) y para el análisis multivariante los paquetes ad4 y Rcommander (R Development Core Team, Open Source). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 34 4. Resultados 4.1 Evolución espacio – temporal de los parámetros físico - químicos A pesar de que no ha sido posible disponer de información detallada del caudal en cada punto de muestreo, sí que se ha tenido acceso a los datos de caudales procedentes de las estaciones de aforo del SAIH dicho organismo, las cuales coinciden prácticamente con los embalses considerados en este estudio. Los datos siguieron una tendencia de un régimen de caudales regulado (Figura 22). Lanuza y Búbal alcanzaron los datos de caudales más bajos a lo largo del año, coincidiendo con estar ubicados en la zona de cabecera. Zuera es la última de las estaciones de aforo sin embargo, no es la que se caracterizó por mayores caudales de agua, sino que son la de Anzánigo y Santa Eulalia las que registraron los máximos valores. En todas las estaciones los mínimos se registraron en los meses de agosto y septiembre y los máximos en Marzo y Junio. Figura 22. Evolución anual de los caudales medios en el río Gállego medidos en las estaciones de aforo de la Confederación Hidrográfica del Ebro (SAIH). Se ha analizado la evolución espacio-temporal de los diferentes parámetros físico – químicos determinados. En general, se ha comprobado que su evolución respecto a los factores estación y punto de muestreo, sigue un comportamiento dentro de lo esperado, excepto algunos casos en los que se presentan anomalías. A continuación se muestran las gráficas de evolución de cada parámetro. Se muestran también los diagramas de caja, aunque sólo para aquellos análisis de varianza que han resultado significativos. La temperatura en el río Gállego fue menor cuanto mayor era la altitud a la que se situaba el punto de muestreo, sin embargo estas pequeñas variaciones no llegaron a ser significativas (F=1,146; gl=3; p= 0,373) entre los puntos de muestreo (Figura 23). Por el contrario sí que se dieron diferencias significativas entre estaciones (F=9,505; gl=5; p= 4,15*10-4) respondiendo a un gradiente de temperaturas donde se alcanzaron los valores más bajos en invierno y los más altos en verano (Figura 24). 0,00 10,00 20,00 30,00 40,00 50,00 60,00 70,00 80,00 Caudal (m3/s) Lanuza Búbal Anzánigo Santa Eulalia Ardisa Zuera Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 35 En invierno la temperatura en El Pueyo incrementó para luego volver a descender. Este mismo punto en verano alcanzó una temperatura menor a la del punto anterior. También en verano, la temperatura descendió en Sabiñánigo. Figura 23. Evolución anual de la temperatura en los seis puntos de muestreo a lo largo del año 2013. Figura 24. Diagrama de cajas de la temperatura (transformación logaritmo en base 10) para las diferentes estaciones del año. Se ha realizado un test ANOVA (F=9,505; gl=3; p=4,15*10-4). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey). Al no seguir una distribución normal (ni con logaritmo ni raíz cuadrada), la variable oxígeno disuelto se analizó mediante un test no paramétrico de Kruskal – Wallis. En general, conforme disminuye la temperatura aumenta la concentración de oxígeno disuelto de modo que se espera que su valor sea mayor en los tramos altos y en las estaciones más frías como otoño e invierno, alcanzando los valores más bajos en verano. No se detectaron diferencias significativas entre los puntos de muestreo (Χ2= 0,76; gl=5; p=0,9795) (Figura 25), mientras que para las estaciones si existieron diferencias (Χ2=19,75; gl=3; p= 1,9*10-4) (Figura 26), coincidiendo con lo observado para la temperatura. Se da un descenso repentino de este parámetro en primavera desde el punto de El Pueyo hasta Santa Eulalia. 0 5 10 15 20 25 Temperatura ºC Invierno Otoño Primavera Verano Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 36 Figura 25. Evolución anual del oxígeno disuelto (mg/L) en los seis puntos de muestreo a lo largo del año 2013. Figura 26. Diagrama de cajas del oxígeno disuelto para las diferentes estaciones del año. Se ha realizado test de Kruskal – Wallis (Χ2=19,75; gl=3; p= 1,9*10-4). La conductividad no mostró diferencias significativas entre las estaciones (F=0,911; gl=3; p=0,453), (Figura 27). Si lo hizo entre los puntos de muestreo (F=17,49; gl=5; p=2,36*10-6) destacando el punto de Sabiñánigo, en el cual se observó un aumento repentino de la conductividad en todas las estaciones, que se va reduciendo progresivamente río abajo. En este sentido los tres puntos de muestreo río abajo estaban claramente diferenciados de los tres puntos río arriba, con una menor conductividad (Figura 28). 8 9 10 11 12 13 14 O2 disuelto mg/L Invierno Otoño Primavera Verano Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 37 Figura 27. Evolución anual de la conductividad en los seis puntos de muestreo del Río Gállego. Figura 28. Diagrama de cajas de la conductividad para los diferentes puntos de muestreo. Se ha realizado un test ANOVA (F=17,49; gl=3; p=2,36*10-6). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey). Los datos del pH se midieron en las estaciones de invierno, primavera y verano (Figura 29). Se observa como en los últimos tres puntos de muestreo los valores permanecieron constantes en comparación con los tres primeros, donde variaban tanto entre estaciones como entre puntos de muestreo. Sin embargo, estas variaciones no resultaron significativas en ninguno de los dos casos (F=3,041; gl=3; p=0,0779 y F=1,197; gl=5; p=0,367 respectivamente). Resaltar el incremento de pH que tuvo lugar en invierno entre Formigal y El Pueyo, así como el elevado valor que alcanzó Búbal en primavera. 0 100 200 300 400 500 600 Condcutividad µS/cm Invierno Otoño Primavera Verano Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 38 Figura 29. Evolución del pH en las estaciones de Invierno 2012 – Verano 2013 en los seis puntos de muestreo del río Gállego. 4.2 Evolución espacio – temporal del estado trófico Los valores de concentración de clorofila bentónica se han comparado con los valores de referencia de estado trófico propuestos por Dodds et al., (1998). Como puede observarse en la Figura 30, la mayoría de los valores de clorofila se sitúan en el rango característico de un río oligotrófico, por debajo de 20mg/m2. Sin embargo, el tramo de Búbal adquirió concentraciones de clorofila propias de un tramo mesotrófico en las estaciones de otoño y primavera, y eutrófico en invierno. Concretamente alcanza los 163,23 mg/m2 estando el límite mesotrofia – eutrofia en 70 mg/m2. También el punto de Formigal obtuvo valores propios de un sistema mesotrófico en otoño e invierno La variable no presentó diferencias significativas entre los puntos de muestreo ni entre las estaciones (Tabla 6), a causa de la alta variabilidad del punto de Búbal, encontrándose muy cerca del límite de significancia En la Figura 31, se muestra la evolución estacional del logaritmo de la clorofila a. Pueden observarse en invierno dos valores atípicos del conjunto de datos de la clorofila en invierno, por una parte el de Búbal (9) y por otra parte el de Sabiñánigo (10). 7,40 7,60 7,80 8,00 8,20 8,40 8,60 8,80 9,00 9,20 pH Primavera Verano Invierno Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 39 Figura 30. Representación de los valores de clorofila a (mg/m2) en cada punto de muestreo en las distintas estaciones del año, con los límites de oligotrofia, mesotrofia y eutrofia establecidos por Dodds et al., (1998). Clorofila F P gl n Estación 2,7121 0,0716 5 4 Punto de Muestreo 2,458 0,0727 3 6 Tabla 6. Valores del estadístico (F), pvalue, grados de libertad (gl) y n para la variable de concentración de clorofila respecto a los factores estación y punto de muestreo del río Gállego, obtenidos del test ANOVA. Figura 31. Diagrama de cajas de la conductividad para las diferentes estaciones del año. Se ha realizado un test ANOVA (F=2,721; gl=3; p=0,0716). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey). Valores atípicos (caso 9)=Búbal; (caso 10)=Sabiñánigo. 0 20 40 60 80 100 120 140 160 180 Clorofila a (mg/m2) Otoño Invierno Primavera Verano Límite Oligotrofia - Mesotrofia Límite Mesotrofia - Eutrofia Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 40 4.3 Estado ecológico del río Gállego: Diversidad taxonómica, IBMWP y riqueza 4.3.1. Evolución espacio – temporal de los índices de calidad y diversidad. Los datos referentes a las familias presentes en cada punto, sus abundancias y valores de índice IBMWP se muestran en Tablas 2 y 3, Anexo y los valores de los índices biológicos Tablas 4, 5 y 6, Anexo. El índice IBMWP de calidad ecológica mostró una distribución sin grandes variaciones (F=0,183; gl=3; p=0,907) en lo que a estaciones del año se refiere (Figura 32). Las medianas comprendieron valores entre 60 y 100 al tratar los datos de todos los puntos en conjunto, lo que implica una calidad de agua buena (CHE, 2009). No es así tratándose de los puntos de muestreo (Figura 33), donde las diferencias fueron claras y significativas entre ellos (F=42,69; gl=5; p = 2,3*10-9). Destaca Búbal por su muy buena calidad, seguida por los puntos de Santa Eulalia y Formigal. En Gurrea de Gállego se obtuvo un agua con buena calidad, mientras que en Sabiñánigo empeoró siendo aceptable únicamente. Se observó un descenso importante en el índice en el punto de El Pueyo que se mantuvo constante en todas las estaciones y permitió clasificarlo como un tramo con calidad de aguas deficiente. Figura 32. Diagrama de cajas de la índice IBMWP para las diferentes estaciones del año. Se ha realizado un test ANOVA (F=0,183; gl=3; p=0,907). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey). Figura 33. Diagrama de cajas del índice IBMWP respecto a los diferentes puntos de muestreo del río Gállego (F=42,69; gl=5; p = 2,3*10-9) realizado mediante test ANOVA. Las letras representan la comparación de medias dos a dos (Test post – hoc Tukey), situándose en los extremos El Pueyo (a) y Búbal (d). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 41 La riqueza taxonómica siguió una distribución similar a la del índice IBMWP tanto a lo largo de las estaciones (F=0,338; gl=3; p=0,798) como atendiendo al punto de muestreo (F=28,66; gl=5; p= 5,65*10-8). Respecto a la estación se observó el valor máximo (caso 15 Figura 34b) de Búbal en primavera (Figura 34, a), al situarse muy por encima del resto de los valores del conjunto de datos de dicha estación. En la Figura 34b, Búbal sobresale marcando los valores máximos de riqueza taxonómica, seguido de Santa Eulalia. Formigal y Gurrea se sitúan en posiciones intermedias mientras que Sabiñánigo y El Pueyo, al igual que sucedía con el IBMWP, obtuvieron los valores más bajos. Figura 34. (a) Diagrama de cajas de la riqueza para las diferentes estaciones del año. Se ha realizado un test ANOVA (F=0,338; gl=3; p=0,798). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey).Valor atípico (caso 15)= Búbal (b) Diagrama de cajas de la riqueza para los diferentes puntos de muestreo. Se ha realizado un test ANOVA (F=28,66; gl=5; p= 5,65*10-8). Las letras representan las comparaciones de medias dos a dos (Test post – hoc Tukey). La diversidad (Figura 35) alcanzó los mínimos y máximos valores en verano. En este caso tampoco se observaron diferencias significativas (F=0,551; gl=3; p =0,653) entre las estaciones del año, a causa de la amplia variabilidad entre los puntos de muestreo (F=7,751; gl=5; p =4,84*10-4). Cabe destacar el valor extremo (caso 20, Figura 35) perteneciente al Pueyo en la estación de verano, donde sólo se encontró un individuo de la familia Chironomidae. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 48 Por último se llevó a cabo un análisis co – inercia en la Figura 41, donde se representan en conjunto las familias presentes y su distribución en los puntos de muestreo y estaciones. Sobre los ejes del FCA, para análisis de los rasgos biológicos. La varianzas explicada por los dos primeros ejes fue del 34,358%. Figura 41. Análisis co – inercia que muestra la estructura de la distribución de las familias de macroinvertebrados y sus rasgos biológicos (tolerancia a la contaminación, grado de eutrofia, estrategia trófica y microhábitat) en el río Gállego. Se ha modificado la gráfica para facilitar la lectura de modo que en los EJE Y (PCA) y EJE X (FCA) las puntos de muestreo y rasgos biológicos que se encuentran dentro del mismo recuadro se encontraban superpuestos. En la gráfica donde se representan las familias los números corresponden a las familias que se encontraban superpuestas. No se ha utilizado la cursiva en las familias por motivos técnicos. (1= Chloroperlidae Glossosomatidae, Hydropsychidae, Limnephilidae, Perlidae y Perlodidae y; 2= Beraeidae, Empididae; 3= Astacidae, Atyidae, Bithyniidae, Capniidae, Ephemeridae, Ephemerillidae, Gammaridae, Leptoceridae, Limoniidae Lymnaeidae, Nemouridae, Oligoneuriidae, Phsychomyidae, Physidae, Planorbidae, Sericostomatidae, Stratiomyidae y Tipulidae; 4= Baetidae; 5=Caenidae, Hydraenidae, Hydrophilidae) Valores propios (Eje 1=0,2064; Eje2=0,2005). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 49 En el eje 1 (17,411% de varianza) positivo destacan en el extremo las familias Raghionidae, con estrategia perforadora y los puntos de muestreo de Santa Eulalia en otoño y verano y Sabiñánigo en Invierno. También se encuentran coleópteros como Gyrinidae, Hydroptilidae e Haliplidae. Los valores negativos de este eje vienen marcados por predadores como Rhyacophilidae, raspadores como Elmidae y filtradores como Hydropsychidae. En el eje 2 (16,947% de varianza) positivo marcado por el requerimiento de oligotrofia y la poca tolerancia a la contaminación, abundan micrófitos y macrófitos encontramos a los puntos de Búbal en otoño y verano y Sabiñánigo en otoño e invierno. De acuerdo con los rasgos descritos, individuos de la orden Trichoptera y Plecoptera son abundantes en estas muestras. En el segundo eje se observa una clara diferencia entre los puntos de muestreo, las familias y sus rasgos. En el extremo negativo los puntos de Formigal, El Pueyo y el punto de Búbal en invierno, vienen marcados por un ambiente eutrófico, donde la estrategia trófica detritívora y la preferencia por los sedimentos son predominantes y con alta tolerancia a la contaminación. Las familias características de estos puntos son Oligochaeta, Potamanthidae, Heptageniidae y Chironomidae, entre otras. Los puntos de Gurrea y Santa Eulalia comparten la eutrofia – mesotrofia de estos puntos además estar dominados por la estrategia trófica de trituradores, el microhábitat es de gravas, arenas y limos y las familias predominantes Leptophiebiidae, Caenidae, Hydraenidae, Baetidae e Hydrophilidae. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 50 5. Discusión 5.1 Diagnosis del estado ecológico del río Gállego En este apartado se muestra un análisis de los diferentes parámetros e índices bióticos estudiados, relacionándolos con las anomalías y los diferentes impactos antrópicos observados en la cuenca del río Gállego, prestando especial atención al encadenamiento de embalses sobre el estado ecológico. 5.1.1 Parámetros físico – químicos El río Gállego desde Agosto del 2012 al 2013 presenta un claro régimen de caudales de río regulado. Se puede observar especialmente en la evolución anual de Ardisa, pues funciona como un azud de derivación para desviar las aguas del Gállego hacia los Riegos del Alto Aragón (García & Moreno, 2000). Pese a ello, los embalses de Lanuza, Búbal y La Peña presentan una capacidad reguladora moderada (García & Moreno, 2003). Respecto a la temperatura, en El Pueyo en invierno ésta aumenta para luego descender en Búbal, lo que podría deberse al entorno de dichos puntos de muestreo. El Pueyo está rodeado de pastos, lo que genera que la luz incida directamente aumentando así la temperatura del agua (Sponseller et al., 2001). Búbal se encuentra rodeado de vegetación y en un congosto donde en Invierno son pocas las horas de luz que recibe. No hay que obviar el posible efecto en este parámetro de la tubería que vierte a poca distancia del punto de El Pueyo. El descenso de temperatura en verano en Sabiñánigo, coincidiendo con un color del agua marrón – grisáceo y un olor característico, puede deberse al bypass de la depuradora situada aguas arriba, a algún vertido no registrado o al desembalse de aguas del fondo del embalse de Sabiñánigo, ya que en verano el agua del fondo de la presa está a menor temperatura generando un impacto en los hábitats existentes aguas abajo (Elosegui & Sabater, 2009). La Conductividad de un río viene marcada por diversos factores. A lo largo de su curso el cauce va recorriendo diferentes materiales geológicos que le aportarán distintas sales en función a su composición. Así mismo, recogerá sales procedentes de fuentes y también de vertidos urbanos, industriales, etc. Otro foco que puede aumentar la conductividad del agua es la agricultura, mediante infiltración y lixiviado de fertilizantes y abonos. En el presente estudio, la conductividad no varía significativamente a lo largo del año, sin embargo, sí se produce un incremento de la conductividad en Sabiñánigo respecto a los demás puntos de muestreo, que puede deberse a la influencia directa tanto de la población y sus industrias como de la depuradora situada aguas arriba. Directamente relacionado con la temperatura, la altitud, los procesos de producción primaria, la descomposición de la materia orgánica y la densidad de vegetación acuática se encuentra el oxígeno disuelto en agua (mg/L) (Ortiz, 2002). Se ha observado que este parámetro experimenta un notable descenso en primavera desde El Pueyo hasta Santa Eulalia coincidiendo con pH alto respecto al resto del año. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 51 El oxígeno disuelto disminuye al darse un exceso de materia orgánica que puede estar ligada principalmente a vertidos procedentes de núcleos urbanos, aguas procedentes de actividades ganaderas o industriales sin el tratamiento de depuración adecuada (Alonso & Camargo, 2005) o bien a la acumulación de biomasa en los embalses situados aguas arriba. A su vez, el pH de un río refleja el tipo de sustrato geológico y la actividad biológica, aumentando cuando las tasas fotosintéticas y la producción primaria son altas y disminuyendo en condiciones de intensa descomposición (Elosegui & Sabater, 2009). Ambos parámetros en conjunto parecen indicar que en estos puntos de muestreo (entre El Pueyo y Santa Eulalia) durante la primavera, estación de notable afluencia de población en la cabecera del Gállego debido a la estación de esquí de Formigal y el turismo rural, se da un incremento de la eutrofización que repercute aguas abajo. 5.1.2 Estado trófico Los valores de clorofila a en el perifiton registrados a lo largo del año permiten clasificar el río Gállego entre oligotrófico y mesotrófico. La altitud, la temperatura, la velocidad del flujo, la vegetación riparia, el sustrato y la luz son algunos de los factores más importantes que influyen en este parámetro (Young et al., 2008). Sin embargo son la materia orgánica y los nutrientes los que estimulan en mayor grado el crecimiento y desarrollo del perifiton (Lohman et al., 1992; Welch et al., 1992; Van Niewenhuyse & Jones, 1996; Dodds et al., 1997; Dodds et al., 1998). Según las condiciones de luz y temperatura, la primavera es la época más propicia para el desarrollo del perifiton (Riaño et al., 1993) hecho que no se puede observar en el río Gállego, ya que en invierno los valores son más altos respecto a las demás estaciones. Puede asociarse por la alta concentración de clorofila detectada, el punto de Búbal a un elevado grado de eutrofización. Ésta probablemente se deba a la acción de los vertidos de aguas residuales de las poblaciones recogidos desde cabecera, infiltraciones de los pastos donde se practica la ganadería, la alta densidad de vegetación de sus orillas, materia orgánica en descomposición proveniente de la vegetación y material acumulado arrastrado por anteriores crecidas, o desembalse de las aguas del fondo del embalse de Búbal. Mientras los valores permanecen casi constantes en todos los puntos de muestreo, Formigal experimenta variaciones a lo largo del año (aunque en menor grado que Búbal), dándose los máximos en otoño e invierno, cuando la afluencia de público al Portalet y a la estación de esquí de Formigal es máxima, y los vertidos de aguas residuales son mayores. 5.1.3 Evolución de los índices a lo largo del año Los índices estudiados para determinar el estado ecológico (IBMWP, riqueza taxonómica y diversidad) siguen una distribución similar a lo largo del río Gállego, donde se alternan puntos con una muy buena calidad de agua con calidad de aguas deficiente. No se da prácticamente variabilidad estacional en los tres índices. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 52 No ocurre así con los puntos de muestreo, cuyas características se deben al entorno inmediato, ya sea por la peculiaridad natural del punto de muestreo, por las alteraciones hidromorfológicas del embalse situado aguas arriba o por otros impactos antrópicos. Pese a contar con un IBMWP que le otorga un estado de calidad de agua buena, Formigal se ve influido directamente por los vertidos de la estación de esquí y por las ventas del Portalet, sobre todo en invierno donde destaca la alta presencia de Chironomidae con más de 5.800 individuos. Esta familia es una clara indicadora de contaminación orgánica (Alonso & Camargo, 2005). Su abundancia es alta durante todas las estaciones (Figura 42), afectando directamente al valor de la diversidad provocando su disminución, concretamente en invierno. Figura 42. Abundancia (nºindividuos) encontrados de cada familia en Formigal en Otoño y valor en orden decreciente del índice IBMWP de cada una de ellas. (Valor IBMWP Formigal Otoño = 64; calidad de agua buena). El Pueyo permanece con un estado ecológico deficiente a lo largo del todo el año a pesar de estar en un tramo de alta montaña. Se trata de un claro ejemplo de la transformación del río como sistema lótico a léntico por el efecto de los embalses (Scheidegger & Bain, 1995; Díaz - Hernández, 2005). Este punto en verano obtuvo un valor del índice IBMWP de 2 hallándose en la muestra únicamente un individuo de Chironomidae. Ya sea por la cercanía del vertido directo de aguas residuales de la población de El Pueyo, por los vertidos realizados aguas arriba, el efecto del embalse de Lanuza o la central hidroeléctrica cercana este punto de muestro se encuentra en un estado ecológico crítico. Este hecho se agrava si se considera que se sitúa en un tramo de alta montaña, cercano a la cabecera del río por lo que debería caracterizarse por un estado ecológico poco alterado. Búbal presenta una muy buena calidad a lo largo del año, con el valor más alto de riqueza de todos los puntos de muestreo y una alta diversidad. Sin embargo, el predominio de familias como Chironomidae y Elmidae (Figura 43) con altas abundancias indican una gran presencia de materia orgánica (Alonso & Camargo, 2005), lo que coincide con lo observado a partir de los parámetros físico – químicos y el estado trófico medidos en este punto. 0 1 2 3 4 5 6 7 8 9 10 0 100 200 300 400 500 600 700 800 Valor IBMWP Abundancia Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 53 Se remarca por tanto, el hecho de que en este tramo existe un elevado grado de eutrofización ya sea por vertidos acumulados, por efecto de los embalses o por causas naturales. Figura 43. Abundancia (nº individuos) de las familias del punto de Búbal en primavera y valor en orden decreciente del índice IBMWP de cada una de ellas (Valor índice IBMWP Búbal primavera = 166; calidad de agua muy buena). En Sabiñánigo, el índice IBMWP le confiere una calidad ecológica al punto de agua moderada. No obstante, tanto la riqueza como la diversidad presentan valores bajos. Además, la abundancia de individuos de las pocas familias presentes es muy baja. Las depuradoras, generalmente, modifican la calidad biológica del agua tras sus vertidos (Perré, 2010), especialmente si no ha dado tiempo a su dilución (Prat & Munné, 2000; Perré, 2010). Es por ello que la situación espacial concreta de este punto de muestreo podría afectar a su estado ecológico directamente por la cercanía de la EDAR de Sabiñánigo. Sería necesario realizar un muestreo aguas arriba de la depuradora para poder saberlo con exactitud. En la estación de verano (Figura 44), la presencia de más de 1.400 individuos de Oligochaeta y de más de 160 Chironomidae en relación a una abundancia mucho menor del resto de las familias indica algún tipo de contaminación orgánica (Alonso & Camargo, 2005). A esto se le suma el estado del río el día del muestreo con olor característico y color marrón – grisáceo y un descenso de la temperatura. Esta situación pudo ser causada por un bypass de la depuradora situada aguas arriba, que al verse sobrepasada de caudal, tuviera que verter parte del agua residual urbana sin depurar. 0 1 2 3 4 5 6 7 8 9 10 0 200 400 600 800 1000 1200 1400 Caenidae Heptageniidae Leptophiebiidae Perlidae Capniidae Perlodidae Sericostomatidae Glossosomatidae Philopotamidae Psychomyidae Ephemerillidae Nemouridae Rhyacophilidae Elmidae Simuliidae Polymitarcidae Hydroptilidae Hydropsychiidae Empididae Psychodidae Baetidae Gyrinidae Scirtidae Lymnaeidae Planorbidae Chironomidae IBMWP Abundancia Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 54 Figura 44. Abundancia (nº individuos) de las familias del punto de Sabiñánigo en primavera y valor en orden decreciente del índice IBMWP de cada una de ellas (Valor índice IBMWP Sabiñánigo primavera=67; calidad de agua moderada). El punto de Santa Eulalia alcanza los valores más altos de diversidad en todo el curso del río Gállego. Presenta calidades del agua muy buenas y buenas. Sin embargo, existe una gran diferencia de familias entre este punto y Búbal puesto que aquí las predominantes adquieren puntuaciones altas del valor IBMWP como Heptageniidae con un valor de 9 (en otoño más de 1790) y Ephemerillidae con un valor de 7 (Figura 45). Esto puede deberse a que desde el embalse de La Peña hasta el punto de muestreo de Santa Eulalia el río Gállego no sufre grandes afecciones ni impactos humanos, por lo que al ecosistema le da tiempo a recuperarse y mejorar su estado ecológico. Figura 45. Abundancia (nº individuos) de las familias del punto de Santa Eulalia en primavera y valor en orden decreciente del índice IBMWP de cada una de ellas (Valor índice IBMWP Santa Eulalia primavera = 108; calidad de agua muy buena). 0 1 2 3 4 5 6 7 8 9 10 0 200 400 600 800 1000 1200 1400 Athericidae Heptageniidae Capniidae Rhyacophilida e Elmidae Scirticidae Simuliidae Empididae Stratiomydae Baetidae Chironomidae Oligochaeta Valor IBMWP Abundancia 0 1 2 3 4 5 6 7 8 9 10 0 100 200 300 400 500 600 700 Valor IBMWP Abundancia Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 55 Por último, Gurrea se caracteriza por una buena calidad del agua a lo largo del año, con valores intermedios de diversidad y riqueza respecto al resto de los puntos de muestreo (Figura 46). Si bien sus valores no son excepcionalmente altos son considerables para el ecotipo en el que se sitúa. Figura 46. Abundancia (nº individuos) de las familias del punto de Gurrea en primavera y valor en orden decreciente del índice IBMWP de cada una de ellas (Valor índice IBMWP Gurrea primavera = 88; calidad de agua buena). Al analizar el índice de calidad ecológica IBMWP relativo respecto a los diferentes ecotipos de la cuenca del Gállego (IBMWP relativo) se observó cómo puntos que adquirían categorías de agua buena o aceptable se alejan del valor de referencia de IBWMP que deberían en función a su ecotipo. A destacar el punto de El Pueyo que situándose en tramo de alta montaña tendría que caracterizarse por altos valores, así como Sabiñánigo en montaña húmeda. En el otro extremo los puntos de Búbal y especialmente Santa Eulalia, alcanzan valores muy por encima de la categoría muy buena para el ecotipo donde se encuentran. 5.1.3 Comparación de índices de calidad Se obtuvo una correlación positiva entre los tres índices de calidad. Los valores del índice IBMWP y de la riqueza taxonómica en cada punto de muestreo no mostraron variaciones al cambiar de estación, lo que no sucedió con la diversidad. Esto se debe a que la diversidad contempla el número de individuos presente en cada punto y no solamente la presencia de la familia. Se observa por ejemplo en Formigal en Invierno donde la alta presencia de Chironomidae genera un descenso de la diversidad que no se refleja en el índice IBMWP. Según la hipótesis de partida de este trabajo, tras una cadena de embalses, río abajo aumentan, por acumulación, las especies colonizadoras de tramos superiores y la diversidad se incrementa (Margalef, 1983), mientras que el índice IBMWP disminuye al tener en cuenta los requerimientos ecológicos (CHE, 2009). 0 1 2 3 4 5 6 7 8 9 10 0 50 100 150 200 Valor IBMWP Abundancia Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 56 Tras los resultados, se puede decir que la hipótesis en general se cumple. Son los puntos de Búbal y Santa Eulalia, coincidiendo con ser los que mejor estado ecológico poseen, donde se observa claramente este hecho. En Búbal el IBMWP es mayor y va disminuyendo en Santa Eulalia y en Gurrea progresivamente. Por el contrario, la diversidad comienza con bajos valores en Formigal, incrementa en Búbal y alcanza máximos en Santa Eulalia. Además es notable el efecto de la cadena de embalses respecto a la eutrofización que generan (Petts, 1984; Díaz -Hernández, 2005), marcada claramente en el punto de Búbal. Sin embargo, hay que tener en cuenta que no sólo han influido las alteraciones hidromorfológicas de las presas sino también otros impactos antrópicos como la gran contaminación orgánica de Formigal o el estado deficiente de El Pueyo. Así mismo la situación del punto de muestreo de Sabiñánigo tras la depuradora conlleva a que los datos obtenidos no sean del todo objetivos para la comprobación del efecto de los embalses. 5.2 Estructura funcional de la comunidad de macroinvertebrados del río Gállego El análisis de correspondencia difuso (FCA) determina como la presencia de determinadas familias de macroinvertebrados es clara indicadora de la situación del medio. De este modo la presencia de órdenes como Plecoptera y Trichoptera estará ligada a aguas oligotróficas, poco contaminadas y que favorezcan a los trituradores por la presencia de materiales orgánicos procedentes del medio terrestre. Estos grupos han mostrado una alta sensibilidad a la contaminación y degradación de los ecosistemas acuáticos en España (García & González, 1986; García-Criado et al., 1999; Ribera et al., 2002; SánchezFernández et al., 2004; Alonso et al., 2005). En aguas eutróficas, con elevada contaminación, y con una alimentación predominante de sedimentos finos (detritívoros) predominarán Oligochaeta y algunos Ephemeroptera, entre otros. El análisis de componentes principales (PCA) permitió observar cómo se distribuyeron las familias de macroinvertebrados en los puntos de muestreo del río Gállego. La familia más abundante resultó ser Chironomidae, con clara diferencia respecto a las demás, predominando en zona de cabecera durante todas las estaciones del año y destacablemente en invierno. Se trata de una familia generalista y variedad de formas de alimentación dentro de sus especies (Merrit & Cummins, 1978; Hachmoller et al., 1991; Riaño et al., 1993). Si además, se considera la baja abundancia existente de familias indicadoras de buena calidad como las pertenecientes a las órdenes Trichoptera y Plecoptera, se observa un desplazamiento de las especies más sensibles por otras de carácter generalistas en los tramos altos. Esto es habitual aguas abajo de un embalse (Aguirre & Bikuña, 2003) o tras una perturbación del estado natural de un río. Lo mismo ocurrió en la última parte del curso alto y en el medio, donde destacó la alta abundancia de Elmidae. Esta familia, con un carácter bastante generalista tanto en requerimientos tróficos como ambientales (Álvarez, 2009), aumenta con el desarrollo de la producción primaria, al utilizar directamente las algas como alimento (Elosegui & Pozo, 1992; Cummins, 1973; Riaño et al., 1993). Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 57 La parte baja del curso medio y el bajo se caracterizó por las altas abundancias de Heptageniidae, familia indicadora de buena calidad de agua, y Baetidae, que posee ciclos de vida cortos y es común que acompañe a Chironomidae en fases de recolonización tras una perturbación (Robinson et al., 2004; Martínez et al., 2006). En general, el PCA mostró como en el curso alto y medio del río Gállego existe un desequilibrio entre las familias presentes y el número de individuos de cada una. Se da una clara dominancia de especies generalistas (Elmidae, Chironomidae) siendo las menos abundantes pertenecientes a las órdenes indicadoras de muy buena calidad (Trichoptera, Plecoptera). De este modo se recalca la existencia de alteraciones provocadas por la acción antrópica, marcada por una predominancia de taxones tolerantes a la contaminación (Rodríguez et al., 1994; Wood et al, 2000; Morais et al., 2004; Martínez et al., 2006). Por el contrario, el curso bajo del Gállego mostró una distribución de individuos en las familias más equitativa. Por último mediante el análisis co-inercia se pudo analizar la estructura funcional de la comunidad de macroinvertebrados del río Gállego en conjunto, incluyendo sus rasgos biológicos y la composición de la comunidad. En un curso fluvial como el que se trata en este estudio, se puede predecir a grandes rasgos, una estructura trófica común relacionada con las características de sus ecosistemas. En primer lugar, la zona de cabecera debería estar dominada por un alto número de detritívoros y trituradores (Monzón et al., 1991). En la cabecera del río Gállego a lo largo del año predominó la estrategia trófica detritívora (alimentación de sedimentos finos), con la presencia de familias como Chironomidae, Potamanthidae y Oligochaeta, seguida de los trituradores en menor número. Los detritívoros se alimentan en gran medida de las partículas que desprenden los trituradores al procesar el alimento, por lo que el aumento de los primeros indica la existencia de aportes externos de materia orgánica. El resto de las familias de cabecera muestran afinidad por ambientes eutróficos y tolerancia a la contaminación, por tanto serían propias de ecosistemas lóticos de tramos bajos, que reciben muchos aportes alóctonos. Se observa el impacto directo que generan los vertidos de Portalet a lo largo del año sobre el nacimiento del río Gállego, además de la estación de esquí de Formigal en las estaciones de invierno y primavera. En la parte baja del curso alto y medio del río Gállego la estructura funcional fue conformada por familias con afinidades por aguas oligotróficas y poca tolerancia a la contaminación. Estos tramos se caracterizaron en su parte alta por una alta riqueza taxonómica y alto índice IBMWP relacionado con la presencia de un elevado número de familias indicadoras de buena calidad de aguas (como familias del orden Trichoptera y Plecoptera). El FCA, determinó la existencia de una clara dominancia de la familia Elmidae, especialmente en el punto de Búbal. Este coleóptero está ligado a macrófitos y micrófitos como forma mayoritaria de microhábitat y de alimentación (su estrategia trófica es de tipo raspador). Los altos valores de los índices biológicos determinados en este tramo, y la presencia de las familias citadas anteriormente, se relacionan con la gran cantidad de Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 64 Dodds, W. K., Jones, J. R. & Welch, E. B. (1998). Suggested classification of stream trophic state: distributions of temperate stream types by chlorophyll, total nitrogen, and phosphorus. Water Research 32, 1455-1462. Elorza, M. S. (2001). Flora y vegetación arvense y ruderal de la provincia de Huesca. Pp 56 – 81. Elosegui, A. &Sabater, S. (2009). Conceptos y técnicas en ecología fluvial. España. 42 pp. Elosegui, A., & Pozo, J. (1992). Physico-chemical characteristics of the Agiiera River (Spain) during an unusual hydrologically period. Annals. Limnology, 28, 85-96. E. P. A. (1987). Biological criteria for the protection of aquatic life. Ohio, State of Ohio Environmental Protection Agency, Division of Water Quality Planning and Assessment. Extence, C. A. (1981). The effect of drought on benthic invertebrate communities in a lowland river. Hydrobiologia, 83(2), 217-224. Fernández, J., Arjol, M. A., & Cacho, C. (2013). POP-contaminated sites from HCH production in Sabiñánigo, Spain. Environmental science and pollution research international, 20(4), 1937-1950. Figueroa, R., Valdovinos, C., Araya, E., & Parra, O. (2003). Macroinvertebrados bentónicos como indicadores de calidad de agua de ríos del sur de Chile. Revista chilena de historia natural, 76(2), 275-285. García-Criado, F. (1999). Impacto de la minería del carbón sobre Hydraenidae y Elmidae (Coleoptera) en la cuenca del río Sil (León, España). (Doctoral dissertation, Ph. D. Thesis. Universidad de León). García de Jalón, D., & González del Tánago, M. (1986). Métodos biológicos para el estudio de la calidad de las aguas. Aplicación a la Cuenca del Duero. Monografía, 45. García, M. A., & Moreno, M. C. (2000). Los aprovechamientos en la Cuenca del Ebro: Afección en el régimen hidrológico fluvial (No. 2000-PH, pp. 24-1). Internal Technical Report. Moreno, J. I. L., & Ruiz, J. M. G. (2003). Influencia de los embalses sobre el régimen fluvial en los Pirineos centrales. Cuadernos de investigación geográfica, 29, 7-22. Gibson, G. R., Barbour, M. T., Stribling, J. B., Gerritsen, J., & Karr, J. R. (1996). Biological Criteria: Technical guidance for streams and small rivers. Environmental Protection Agency, Washington, DC (United States). González del Tánago, M., & García de Jalón, D. (2007). Restauración de ríos. Guía metodológica para la elaboración de proyectos. Centro de publicaciones. Secretaría General Técnica Ministerio de Medio Ambiente. Pp 83 - 104 Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 65 Hachmoller, B., Matthews, R. A., & Brakke, D. F. (1991). Effects of riparian community structure, sediment size, and water quality on the macroinvertebrate communities in a small, suburban stream. Northwest Science, 65(3). Hellawell, J. M. (1986). Biological indicators of freshwater pollution and environmental management. Elsevier Applied Sciencie. Publ 546. London & New York. Hernández, L. M., Fernández, M. A., & González, M. J. (1991). Lindane pollution near an industrial source in northeast Spain. Bulletin of environmental contamination and toxicology, 46(1), 9-13. Jaime, J. M. M. (1993). Balance hídrico e hidrológico de la cuenca media del río Gállego. Geographicalia, 30, 243-258. Jaime, J. M. M. (1989). Balance hídrico e hidrológico de la cuenca alta del río Gállego. Geographicalia, 26, 175-182. Jeffrey, S. W., Mantoura, R. F. C. &Wright, S. W. (1996). Phytoplankton pigments in oceanography: guidelines to modern methods. Paris: UNESCO Jiménez, A. A. (2000). Determinación de los parámetros físico-químicos de calidad de las aguas. Revista interdisciplinar de gestión ambiental, 2(23), 12-19. Kondolf, G. M. (1997). PROFILE: hungry water: effects of dams and gravel mining on river channels. Environmental management, 21(4), 533-551. Kurthesi, K., Vehapi, I., Gecaj, A., & Millaku, F. (2008) Periphyton bioindicators in the Sintica river (Kosovo) Annales 18(2) 265 – 270. Li, L., Zheng, B., & Liu, L. (2010). Biomonitoring and bioindicators used for river ecosystems: definitions, approaches and trends. Procedia environmental sciences, 2, 15101524. Ligon, F. K., Dietrich, W. E., & Trush, W. J. (1995). Downstream ecological effects of dams. BioScience, 45(3), 183-192. Lohman, K., Jones, J. R., & Perkins, B. D. (1992). Effects of nutrient enrichment and flood frequency on periphyton biomass in northern Ozark streams. Canadian Journal of Fisheries and Aquatic Sciences, 49(6), 1198-1205. Lundqvist, J. (1998). Avert looming hydrocide. Ambio, 27(6), 428-433. Margalef, R. (1983). Limnología. Barcelona: Omega. (p. 1010) Martínez Bastida, J. J., Arauzo, M., & Valladolid, M. (2006). Diagnóstico de la calidad ambiental del río Oja (La Rioja, España) mediante el análisis de la comunidad de macroinvertebrados bentónicos. Limnetica, 25 (3), 733-744. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 66 Meritt, R. W., & Cummins, K. W. (1978). An introduction to the aquatic insects of North America. Kendall. Minshall, G. W., Petersen, R. C., Cummins, K. W., Bott, T. L., Sedell, J. R., Cushing, C. E., & Vannote, R. L. (1983). Interbiome comparison of stream ecosystem dynamics. Ecological Monographs, 53(1), 1-25. Monserrat, P. (1990). Vegetación. En Ruiz de la Torre et al., Memoria del mapa forestal de España escala 1:200.000 hoja 8 – 2 Viella. Monzón, A., Casado, C., Montes, C., & De Jalón, D. G. (1991). Organización funcional de las comunidades de macroinvertebrados acuáticos de un sistema fluvial de montaña (Sistema Central, rıo Manzanares, España). Limnetica, 7, 97-112. Morais, M., Pinto, P., Guilherme, P., Rosado, J., & Antunes, I. (2004). Assessment of temporary streams: the robustness of metric and multimetric indices under different hydrological conditions. En Integrated Assessment of Running Waters in Europe (pp. 229249). Springer Netherlands. Naiman, R. J; Magnuson, J.J., McKnight D.M., & Stanford J.A. (Ed.). (1995). The freshwater imperative: a research agenda. Island Press. 165pp Ollero, A., Sánchez, M., Marín, J. M., Fernández, D., Ballarín, D., Mora, D., Montorio, R., Beguería, S. & Zúñiga, M. (2004). Caracterización hidromorfológica del río Gállego. Universidad de Zaragoza e IPE CSIC. Geografía Física de Aragón. Aspectos generales y temáticos, 117 – 129. Oscoz, J., Campos, F., & Escala, M. C. (2006). Variación de la comunidad de macroinvertebrados bentónicos en relación con la calidad de las aguas A. Limnetica, 25(3), 683-692. Oscoz, J., Galicia, D., & Miranda, R. (2009). Macroinvertebrados de la Cuenca del Ebro: descripción de taxones y guía de identificación. Universidad de Navarra. Departamento de Zoología y Ecología. Pearson, R. G. (1984). Temporal changes in the composition and abundance of the macroinvertebrate communities of the River Hull. Archiv für Hydrobiologie, 100(3), 273-298. Peña, J. L., & Lozano, M. V. (2004). Las unidades del relieve aragonés. Geografía Física de Aragón, Aspectos Generales Y Temáticos, 3-14. Peña, J. L., Sancho, C., Lewis, C., McDonald, E., & Rhodes, E. (2004). Datos cronológicos de las morrenas terminales del glaciar del Gállego y su relación con las terrazas fluvioglaciares (Pirineo de Huesca). Geografía Física de Aragón. Aspectos generales y temáticos, 71 – 84. Pérez, L., & Hurlé, J. B. (2011). Una nota sobre la recuperación de costes de los servicios del agua en la cuenca del Gállego. Economía Agraria y Recursos, 7(13), 49-56. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 67 Perré, I., Rieradevall, M., Prat, N., Martin, J., & Céspedes, R. (2010). Cambios en el estado ecológico de tres ríos producidos por el vertido de depuradoras. Tecnología del Agua, 320, 21-29. Petts, G.E. (1984). Impounded Rivers. John Wiley & Sons, Chichester Poff, N. L., Allan, J. D., Bain, M. B., Karr, J. R., Prestegaard, K. L., Richter, B. D.,Sparks,R. J & Stromberg, J. C. (1997). The natural flow regime. BioScience, 47(11), 769-784. Prakash, O., Suar, M., Raina, V., Dogra, C., Pal, R., & Lal, R. (2004). Residues of hexachlorocyclohexane isomers in soil and water samples from Delhi and adjoining areas. Current Science, 87(1), 73-77. Prat, N., & Munné, A. (2000). Water use and quality and stream flow in a Mediterranean stream. Water Research, 34(15), 3876-3881. Pujante, A., Martínez-López, F., & Tapia, G. (1995). Macroinvertebrados y calidad de las aguas de los rios próximos a la central térmica de Andorra (Teruel, España). Limnetica, 11(2), 1-8. Purschke, G., Westheide, W., Rohde, D., & Brinkhurst, R. O. (1993). Morphological reinvestigation and phylogenetic relationship of Acanthobdella peledina (Annelida, Clitellata). Zoomorphology, 113(2), 91-101. Reiser, D. W., Ramey, M. P., Beck, S., Lambert, T. R., & Geary, R. E. (1989). Flushing flow recommendations for maintenance of salmonid spawning gravels in a steep, regulated stream. Regulated Rivers: Research & Management, 3(1), 267-275. Riaño, P., Basaguren, A., & Pozo, J. (1993). Variaciones espaciales en las comunidades de macroinvertebrados del río Agüera (País Vasco-Cantabria) en dos épocas con diferentes condiciones de régimen hidrológico. Limnetica, 9, 19-28. Ribera, I., Aguilera, P., Hernando, C., & Millán, A. (2002). Los coleópteros acuáticos de la península Ibérica. Quercus, 201, 38-42. Rieradevall, M., Bonada, N., & Prat, N. (1999). Community structure and water quality in the Mediterranean streams of a natural park (St. Llorenç del Munt, NE Spain). Limnetica, 17, 45-56. Robinson, C. T., & Minshall, G. W. (1986). Effects of disturbance frequency on stream benthic community structure in relation to canopy cover and season. Journal of the North American Benthological Society 5, 237-248. Robinson, C. T., Uehlinger, U. R. S., & Monaghan, M. T. (2004). Stream ecosystem response to multiple experimental floods from a reservoir. River Research and Applications, 20(4), 359-377. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 68 Robles-González, I. V., Rios-Leal, E., Galindez-Mayer, J., Caffarel-Mendez, S., BarreraCortes, J., Esparza-García, F., & Poggi-Varaldo, H. M. (2006). Adsorptive-desorptive behaviour of lindane in an agricultural soil. Interciencia, 31(4), 305-308. Rodríguez, S. E., Fernández-Aláez, C., Fernández-Aláez, M., & Bécares, E. (1994). La cuenca alta del río Cares (NO España). Limnetica, 10(1), 93-100. Rosenberg, D. M., & Resh, V. H. (1993). Freshwater biomonitoring and benthic macroinvertebrates. Chapman & Hall. 488p. Rutherford, J. C., & Cuddy, S. (2005). Modelling periphyton biomass, photosynthesis and respiration in streams. CSIRO Land and Water Technical Report 23/05, Canberra, December 2005, 55p. Sánchez-Fernández, D., Abellán, P., Velasco, J., & Millán, A. (2004). Vulnerabilidad de los coleópteros acuáticos de la Región de Murcia. Revista Ecosistemas, 13(1). Scheidegger, K. J., & Bain, M. B. (1995). Larval fish distribution and microhabitat use in free-flowing and regulated rivers. Copeia, 125-135. Sloto, R. A. (1987). Effect of urbanization on the water resources of Eastern Chester County, Pennsylvania (No. WRI-87-4098). United States Geological Survey. Sparks, R. E. (1995). Need for ecosystem management of large rivers and their floodplains. BioScience, 45(3), 168-182. Sponseller, R. A., Benfield, E. F., & Valett, H. M. (2001). Relationships between land use, spatial scale and stream macroinvertebrate communities.Freshwater Biology, 46(10), 1409-1424. Statzner, B., Bis. B, Dolédec, S., & Usseglio-Polatera, P. (2001). Perspectives for biomonitoring at large spatial scales: a unified measure for the functional composition of invertebrate communities in European running waters. Basic and Applied Ecology, 2(1), 7385. Tachet, H., Richoux, P., Bournaud, M., & Usseglio-Polatera, P. (2000).Invertébrés d'eau douce: systématique, biologie, écologie (pp. 1-588). Paris: CNRS éditions. Turner, M. G., & Dale, V. H. (1998). Comparing large, infrequent disturbances: what have we learned? Ecosystems, 1(6), 493-496. Usseglio-Polatera, P., Bournaud, M., Richoux, P., & Tachet, H. (2000). Biomonitoring through biological traits of benthic macroinvertebrates: how to use species trait databases? In Assessing the Ecological Integrity of Running Waters (pp. 153-162). Springer Netherlands. Van Nieuwenhuyse, E. E., & Jones, J. R. (1996). Phosphorus chlorophyll relationship in temperate streams and its variation with stream catchment area. Canadian Journal of Fisheries and Aquatic Sciences, 53(1), 99-105. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 69 Vannote, R. L., Minshall, G. W., Cummins, K. W., Sedell, J. R., & Cushing, C. E. (1980). The river continuum concept. Canadian journal of fisheries and aquatic sciences, 37(1), 130137. Vivas, S., Casas, J., Pardo, I., Robles, S., Bonada, N., Mellado, A. & Moyá, G. (2002). Aproximación multivariante en la exploración de la tolerancia ambiental de las familias de macroinvertebrados de los ríos mediterráneos del proyecto GUADALMED. Limnetica, 21(34), 149-173. Ward, J. V. (1976). Effects of flow patterns below large dams on stream benthos: a review. Instream flow needs 2, 235 - 253. Ward, J. V., & Stanford, J. A. (Eds.). (1979). The ecology of regulated streams. New York: Plenum Press. (p. 398). Ward, J. V. (1989). The four-dimensional nature of lotic ecosystems. Journal of the North American Benthological Society 8, 2-8. Ward, J. V., & Stanford, J. A. (1995). Ecological connectivity in alluvial river ecosystems and its disruption by flow regulation. Regulated Rivers: Research & Management, 11(1), 105119. Ward, J. V. (1998). Riverine landscapes: biodiversity patterns, disturbance regimes, and aquatic conservation. Biological conservation, 83(3), 269-278. Ward, J. V., K. Tockner, & F. Schiemer. (1999). Biodiversity of floodplain ecosystems: Ecotones and connectivity. Regulated Rivers: Research and Management 15, 125–139. Ward, J. V., & Wiens, J. A. (2001). Ecotones of riverine ecosystems: role and typology, spatio-temporal dynamics, and river regulation. Ecohydrology & hydrobiology, 1(1-2), 2536. Weaver, L. A., & Garman, G. C. (1994). Urbanization of a watershed and historical changes in a stream fish assemblage. Transactions of the American Fisheries Society, 123(2), 162172. Welch, E. B., Quinn, J. M., & Hickey, C. W. (1992). Periphyton biomass related to pointsource nutrient enrichment in seven New Zealand streams. Water Research, 26(5), 669675. Wetzel, R. G. (2001). Limnology: lake and river ecosystems. (3rd edition). Academic Press, San Diego. Wood, P. J., Agnew, M. D., & Petts, G. E. (2000). Flow variations and macroinvertebrate community responses in a small groundwater‐dominated stream in south‐east England. Hydrological Processes, 14(16‐17), 3133-3147. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 70 Working, F. I. S. C. R. Group (FISRWG). (1998). Stream corridor restoration: Principles, processes, and practices. GPO Item, (0120-A). Young, R. G., Matthaei, C. D., & Townsend, C. R. (2008). Organic matter breakdown and ecosystem metabolism: functional indicators for assessing river ecosystem health. Journal of the North American Benthological Society, 27(3), 605-625. Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 71 9. Anexo Tabla 1. Rasgos biológicos de las familias de macroinvertebrados, descripción y categorías según Tachet et al., (2000). Rasgo biológico Descripción Categoría Eutrofia Mesotrofia Oligotrofia Sedimentos finos Triturador Raspador Filtrador Predador Perforador Piedras Gravas Arenas Limos Macrófitos Micrófitos Ramas Lecho Pozas X1 - Intolerante X2 - Poco tolerante X3 - Relativamente tolerante X4 - Tolerante X5 - Muy Tolerante Grado de Eutrofia Alimentación Microhábitat Tolerancia a la contaminación Este rasgo se relaciona directamente con la química del agua y con los nutrientes contenidos en ella. Organismos resistentes a altas concentraciones serán eutróficos mientras los más sensibles serán oligotróficos. Los macroinvertebrados detritívoros se alimentan de sedimentos finos, los filtradores captan el alimento particulado presente en el agua, los trituradores arrancan trozos de macrófitos, los raspadores con sus piezas bucales rompen y raspan alimentos como los micrófitos. Los predadores se alimentan de otros macroinvertebrados u organismos vivos mientras que los perforadores poseen un aparato que les permite succionar tanto elementos vegetales como partes de otros macroinvertebrados. Los diferentes tipos de microhábitat se distribuyen en nueve modalidades diferentes, pudiendo tener alta afinidad cada familia por varios de ellos a la vez. Hace referencia a la resistencia de los organismos ante la contaminación, especialmente de carácter orgánico. Los macroinvertebrados podrán presentar desde una gran sensibilidad a la misma siendo X1 - intolerantes, pasando por caracteres intermedios hasta tolerantes a la misma, siendo capaces de desarrollarse en medios contaminados Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 72 Familia Valor IBMWP G1O G2O G3O G4O G5O G6O G1I G2I G3I G4I G5I G6I Elmidae 5 5 0 1461 131 318 222 0279 141 138 5 Gyriniidae 3 0 0 9 0 16 000009 5 Haliplidae 4 0 0 2 0 0 0 0 0 1 0 0 0 Hidrophilidae 3 0 1 1 0 0 0 0 0 0 0 0 0 Hydraenidae 5 0 0 0 0 0 0 0 0 0 0 0 0 Scirtidae 3 0 0 3 0 00001000 Chironomidae 2 778 148 428 27 5486 11 1393 8260 332 Athericidae 10 0 0 0 0 3 1 0 0 3 3 0 0 Ceratopogonidae 4 0 0 0 3 0 1 0 0 0 000 Empididae 4 0 0 0 0 0 0 0 0 4 0 0 0 Limoniidae 4 3 0 0 0 00100000 Pyschodidae 4 0 0 0 0 00100300 Rhagionidae 4 0 0 0 0 00000000 Simuliidae 5 0 0 46 011 015 026 22 65 1 Stratyomidae 4 0 0 0 0 00000000 Tabanidae 4 0 0 0 0 00000100 Tipulidae 5 0 0 0 0 01100000 Baetidae 4 60 0137 187 7610 0382 0153 13 Caenidae 4 0 0 0 0 245 2000035 3 Ephemeridae 10 00200000100 0 Ephemerillidae 7 0 0 72 0 0 0 0 0 42 00 0 Heptageniidae 10 83 0376 01794 13 44 0115 0244 157 Leptophiebiidae 10 0 0 67 000 0 0 14 00 0 Oligoneuriidae 5 0 0 0 0 0000000 0 Polymitarcidae 5 0 0 0 0 00002000 Potamanthidae 10 0 0 0 0 20 12 0 0 0 0 4 1 Chloroperlidae 10 54 0 0 0 0 0 33 01000 Capniidae 10 0 0 3 03 0 24 03000 Leuctridae 10 0 0 1 0 00000000 Nemouridae 7 0 0 2 0 0 0 3 0 40 01 0 Perlidae 10 0 0 53 000 1 0 47 00 0 Perlodidae 10 0 0 0 0 0035 07010 Hydropsychidae 5 19 096 4790 1 7 0 220 354 20 Beraeidae 10 0 0 0 0 0 0 0 0 0 1 0 0 Glossosomatidae 8 0 0 9 0 0 0 0 0 1 00 0 Hydroptilidae 6 0 0 0 0 00000000 Leptoceridae 10 0 0 1 0 00 0 0 12 0 0 0 Limnephilidae 7 3 0 111 000500000 Philopotamidae 8 0 0 1 0 00000000 Polycentropodidae 7 0 0 0 0 00000000 Psychomyidae 8 0 0 0 0 00000000 Rhyacophilidae 7 8 0 22 1 6 0 12 051 0 0 1 Sericostomatidae 10 0 0 6 0 00 0 0 19 0 0 0 Planorbidae 3 0 1 5 0 00000110 Ancylidae 6 0 0 0 0 0 0 0 1 0 0 0 0 Bithyniidae 3 0 1 0 0 0 0 0 1 0 1 0 0 Hydrobiidae 3 0 0 0 0 2 0000040 Lymnaeidae 3 0 0 1 0 00000020 Phsydae 3 0 0 0 0 1 0 0 0 0 01 0 Astacidae 8 0 0 0 0 0 0 0 0 0 0 0 0 Atyidae 6 0 0 0 0 1 5 0 0 0 0 0 1 Gammaridae 6 0 0 0 0 232 67 000075 61 Aeshnidae 8 0 0 0 0 1 0 0 0 0 0 0 0 Calopterygidae 8 0 0 0 020000000 Gomphidae 8 0 0 0 0 6 0 0 0 0 00 0 Oligochaeta 1 0 18 0 0 01100200 Tabla 2. Abundancias y valor IBMWP de las familias presentes en las muestras de otoño e invierno Estado trófico y calidad ecológica del río Gállego a su paso por la provincia de Huesca - 2014 73 Familia Valor IBMWP G1P G2P G3P G4P G5P G6P G1V G2V G3V G4V G5V G6V Elmidae 5 30 11392 89 86 550 0612 28 37 5 Gyriniidae 3 0 0702 1 0 0 2 0 00 Haliplidae 4 0 0000 0 0 0 0 0 00 Hidrophilidae 3 0 0901 0 1 0 0 0 00 Hydraenidae 5 1 0 0 00 0 0 0 0 0 00 Scirtidae 3 1 0 2 1 0 0 0 0 14 0 0 0 Chironomidae 2 382 24 763 169 166 41 128 3134 20 23 16 Athericidae 10 0 0 0 1 0 0 0 0 0 0 0 0 Ceratopogonidae 4 0 00 0 0 0 0 0 00 1 0 Empididae 4 0 0 34 100002000 Limoniidae 4 10 2 0 0 0 0 0 0 0 1 0 0 Pyschodidae 4 0 0 1 0 0 0 0 0 0 1 0 0 Rhagionidae 4 1 0 0 0 0 0 0 0 0 0 0 0 Simuliidae 5 0 0 53 62 055 12 073 423 7 Stratyomidae 4 0 0 0 1 0 0 0 0 0 0 0 0 Tabanidae 4 0 0 0 0 0 0 0 0 0 0 0 0 Tipulidae 5 0 5 0 0 0 0 0 0 0 0 0 0 Baetidae 4 93 098 5326 122 77 0228 385 83 Caenidae 4 0 0 118 081 39 1 0 63 024 0 Ephemeridae 10 0 0 0 0 0 00 0 00 0 0 Ephemerillidae 7 0 0 33 0695 6008018 0 Heptageniidae 10 38 0155 1414 195 29 047 068 251 Leptophiebiidae 10 0 0 16 00 0 0 0 0 0 0 0 Oligoneuriidae 5 00 0 0 8 11 0 0 0 0 24 2 Polymitarcidae 5 0 0 4 0 0 3 0 0 1 0 1 15 Potamanthidae 10 0 0 0 0 41 0000019 0 Chloroperlidae 10 55 00 0 0 0 6 0 00 0 0 Capniidae 10 6095 3 6 7 8 0 222 020 0 Leuctridae 10 0 0 0 0 15 0000015 1 Nemouridae 7 2 0 10 000004000 Perlidae 10 0 0 63 0 0 0 0 0 35 0 0 0 Perlodidae 10 9 0 28 010 0803012 0 Hydropsychidae 5 1 1 25 039 1 1 0 24 012 96 Beraeidae 10 0 0 0 0 0 0 0 0 00 0 0 Glossosomatidae 8 0 0 1 00 0 0 0 0000 Hydroptilidae 6 0 0 0 00 0 0 0 0 0 00 Leptoceridae 10 00 0 00 0 0 0 0 0 00 Limnephilidae 7 10 2 0 0 0 0 4 0 1 0 0 0 Philopotamidae 8 0 0 8 0 0 0 0 0 1 0 00 Polycentropodidae 7 0 0 0 0 4 0 0 0 0 0 0 0 Psychomyidae 8 0 0 6 0 0 0 0 0 2 0 0 0 Rhyacophilidae 7 2 0 14 2 6 3 0 0 24 0 2 0 Sericostomatidae 10 3 0 13 000002001 Planorbidae 3 0 0 1 0 0 0 0 0 2 1 0 0 Ancylidae 6 0 0 0 0 0 0 0 0 0 0 0 0 Bithyniidae 3 0 1 0 0 0 0 0 0 00 0 0 Hydrobiidae 3 0 0000 2 0 0 0 0 12 0 Lymnaeidae 3 0 0 5 0 0 0 0 0 2 0 0 0 Phsydae 3 0 0 0 0 0 0 0 0 0 0 0 0 Astacidae 8 0 0 0 0 0 0 0 0 0 0 0 0 Atyidae 6 0 0 0 0 0 1 0 0 0 0 0 0 Gammaridae 6 0 00 0 279 226 0 0 0 0 46 35 Aeshnidae 8 0 0 0 0 0 0 0 0 0 0 0 0 Calopterygidae 8 0 0 0 0 0 00 0 00 0 0 Gomphidae 8 0 0 0 0 0 0 0 0 0000 Oligochaeta 1 4 58 01444 0 4 0 0 0 3 00 Tabla 3. Abundancias y valor IBMWP de las familias presentes en las muestras de primavera y verano