Las comunidades de macroinvertebrados bentonicos (MIB) empleados como bioindicadores acuaticos vienen aumentando en estos ultimos años en lo que respecta a la protecciòn y restauracion de los ambientes acuaticos (Brown et al. 1997; Cole 1998; Wetzel y Likens 2000; Acosta 2001, Marques et al. 2001; Metzeling et al. 2003). En contraste con los analisis fisicoquimicos y microbiologicos, que sslo proporcionan informacion puntual e indirecta, la evaluacion de las comunidades de macroinvertebrados en los ecosistemas acuaticos, con enfasis en insectos, proporcionan unaexcelente alternaliva en el diagnostico de la calidad del agua (Huryn y Wallace 2000; Baptista et al. 2001; Galdean et al. 2001; Rogers et al. 2002; Figueroa et al. 2003).

Muchos representantes de estos grupos de invertebrados pueden distribuirse en distintos tipos de condiciones ambientales, asi como ser muy susceptibles a la contaminacion, motivo por el cual son recomendados en el establecimicnto de sistemas de vigilancia y control de los ecosistemas hidricos (Leslie et al. 1999; Ogbeibu y Oribhabor 2002; Riss et al. 2002; Iannacone et al. 2003). El Peru es uno de los paises Neotropicales que enfrenta serios problemas de contaminacion de sus ambientes loticos (Paredes et al. 2004).

Se han empleado los indices de diversidad y bioticos usando a los insectos acuaticos para caracterizar biologicamente los cuerpos de agua dulce (Iannacone et al. 2000; Moreno 2001; Fenoglio et al. 2002; Goncalves y Aranha 2004; Paredes et al. 2004). El objetivo del pre-sente trabajo fue analizar la composiciòn faunistica, riqueza de familias y calidad del agua empleando el indice biotico BMWP' modif. (Biological Monitoring Working Party modificado) en el rio Rimac aguas abajo de la Atarjea, Callao-Lima, Peru durante el periodo Octubre 2002 a Seliembre del 2003.

Materiales y Metodos

Zona de estudio: Incluyeron seis estaciones de muestreoestablecidas en la trayectoria del rio Rimac desde las inmediaciones del puente Grana (Estacion 6), cercana a la sede de caplacion o bocatoma de la Atarjea de SEDAPAL en el distrito de El Agustino, Lima, Peru hasta las cercanias del puente Gambeta (Estaciòn 1), en el Callao, proxima a la desembocadura al mar (Fig. 1) Este recorrido tiene aproximadamente 17.5 Km de longitud. con una pendiente aproximada de 1.1 % y disminuyendo desde los 195 a los 0 msnm (Fig. 1). Se determino la ubicacion geografica de cada una de las estaciones de muestreo con un GPS 12 Garmin (Fig. 1). El rio Rimac es considerado de cuarto orden segun la clasificaciòn de Strahler (Wetzel y Likens 2000).

Muestreos: Se realizaron muestreos mensuales entre Octubre 2002 y Setiembre 2003. Los analisis fisicoquimicos incluyeron mediciones de: pH y conductividad electrica (CE: ohm.cm^-2) mediante instrumentales portables marca Hanna®, temperatura en °C mediante un termometro de mercurio. Los datos de caudal (m³.s^-1) y de velocidad superficial promedio (m.s^-1) del rio se obtuvieron siguiendo el metodo propuesto por Lind (1979). El Oxigeno Disuelto (OD) y la Demanda Bioquimica de Oxigeno (DBO_5,20 )en mg L^-1 siguieron el protocolo propuesto por Winkler(APHA 1995).

Para tomar las muestras biologicas en las estaciones de estudio se emplearon: el metodo manual y la tecnica rapida de Hynes modificada, con un muestreador Surber modificado, una red estandar de 30 x 30 cm² con abertura de 1 mm (Paredes et al. 2004). Cada muestra biologica fue separada a nivel de familia y se fijaron con EtOH 70 % en viales de vidrio de 5 ml debidamente rotulados. Para los moluscos se siguio el procedimiento descrito por Thiengo et al. (2001) y para la identificaciòn de los Hexapoda se usaron claves taxondmicas especializadas (Stehr 1987; Borroretui. 1995; Merritty Cummins 1996). Para la determinacion de la composiciòn faumistica y la riqueza de familias de los macroinvertebrados bentdnicos se empleo el nivel taxonomico de familia por su relativa facilidad y rapidez, y por su capacidad de aportar informaciòn ecologica (Metzeling et al. 2003). Material representativo se deposito en el Museo de Entomologia del Programa Nacional de Control Biologico (PNCB), Lima, Peru.

Analisis de datos: Se empleo el analisis de componentes principales (ACP) como un criterio de reducciòn y ordenaciòn de los siete parametros fisicoquimicos y de los catorce indices de diversidad alfa, obtenidos de cada una de las seis estaciones de muestreo, para producir variables compuestas no relacionadas. Se realizo la rotacion varimax para facilitar la interpretaciòn de los componentes del ACP (Zar 1996; Vivanco 1999). Se utilizo el coeficiente de correlaciòn de Pearson entre los indices bioticos y algunos taxa mas abundantes. Se aplico el indice biotico BMWP' modificado por Zuñiga de Cardoso et al. (1997) citado por Dominguez y Fernandez (1998) propuesto para este estudio, que al igual que el BMWP (Hellawell1978). identifica un nivel de calidad de agua de acuerdo con un puntaje otorgado a ciertas familias de MIB que se encuentren en el cuerpo de agua, con la diferencia que esta modificacion fue adaptada para la Fauna local de la región Neotropical (Tabla 1). Similarmente el puntaje total o por estación corresponde a un valor de calidad previamente establccidu (Tabla 2). También se emplearan los siguientes índices de Marques et al. (2001 ): (Ephemeroptera + Plecoptera + Trichoptera = EPT) / Chironamidae (= Q); EPT / (EPT + Q) y porcentaje de Chironomidae en Diptera. Así como el índice de Goodnight Whittley que se basa en la relacion: (Numero de oligoquetos (Annelida) / Total de organismos benticos) x 100, siendo > 80% como altamente contaminado: 60 -80%, dudoso y < 60% de buena condiciòn (Washington 1984). Con la finalidad de evaluar la diversidad alta en las comunidades maerozoohenlonicas de la zona estudiada, se utilizaron los siguientes indices: Shannon-Wiener (Washington 1984), el indice de diversidad de Margalef (D,_mlt), el indice de equidad de Pielou, los numeros de la serie de Hill, el indice de dominancia de Mcintosh (D) y el indice de Berger-Parker (d) (Moreno 2001). Para la diversidad beta de similitud entre las seis estaciones de la zona de estudio, se uso el indice cualitativo de Sorensen (I_s) y el indice de Sorensen cuantitativo (IS_cuant) (Moreno 200I). Asi mismo. para las familias mas abundantes se determino) la distribucion espacial (al azar, amontonada o uniforme), utilizando el criterio de S²/X (Varianza/ Promedio de individuos por estacion). Para todos los casos se empleo un nivel de significancia de alfa = 0.05 (Zar 1996).

Resultados y discusion

Los promedios de siete parametros fisicoquimicos evaluados mensualmente de Octubre 2002 a Setiembre del 2003 son señalados en la Tabla 3. La DBO_5,20 para las estaciones uno a la cinco, indico perturbacion e impacto en el rio rimac. El valor promedio anual obtenido en este estudio fue mas alto que el valor permisible por la Ley General de Aguas (LGA) peruana vigente. El coeficiente de variacion (CV) se vio incrementado expresamente para este parametro (Tabla 3). La estaciòn dos presento niveles de degradacion ambiental en terminos de DBO_5,20, OD y CL, aparentemente como consecuencia de los efluentes vertidos por las inmediaciones del fondo San Agustin (Callao). Tanto el OD y el pH se encontraron dentro del limite permisible por la LGA peruana vigente. Hay que resaltar que a pesar de la pendiente que presenta el rio Rimac en la zona de estudio (1,1%) y ademas de la presencia de pantallas, que logran incrementar la velocidad superficial del agua y los niveles de oxigeno, estos no son suficientes para oxidar la tremenda carga organica disuelta en sus aguas como lo demuestran los niveles de DBO_5,20 (Tabla 3).

El ACP produjo tres componentes con Eigenvalues > 1 para los parametros fisicoquimicos (Tabla 4). Los parametros fisicoquimicos fueron ordenados en un espacio de dos dimensiones, de acuerdo con el ACP (Fig. 2). La velocidad, DBO, CE y OD del agua estuvieron correlacionados con el primer componente que explico el 64.804 % de la variabilidad. El caudal y pH estuvieron correlacionados con el segundo componente que explico el 18,245 % de la variabilidad (Tabla 4). El tercer componente describio la variacion en la tempcratura del agua (Tabla 4). El CP1 va en relacion con el factor oxigeno. Muy correlacionado al OD se encuentra la DBO, esta ultima a su vez esta relacionada con la carga organica que afecta a la CE. Finalmente. la velocidad del cuerpo de agua permite el ingreso de dosis imponantes de oxigeno superficial al cuerpo de agua lotico. EL CP2 correlaciona el caudal y el pH. El caudal se incrementa o se reduce, por la incorporacion de aguas residuales de distintos tipos que caracterizan las seis estaciones, ya sea industrial, domesitica o agricola, y que indudablemente afectan los niveles de pH (Bispo et al, 2002).

Se obtuvo un puntaje biotico de 29 (Tabla 6), que segun el BMWP' modif. caracteriza al rio Rimac como de calidad critica o de aguas muy contaminadas (Tabla 2). Segun el puntaje obtenido para cada estacion, se puede observar que las estaciones uno y dos presentan aguas fuertemente contaminadas o de calidad muy critica, las estaciones tres y cuatro presentan aguas muy contaminadas o de calidad critica, y las estaciones cinco y seis presentan aguas contaminadas o de calidad dudosa (Tabla 2). Es importante mencionar que solo las estaciones cinco y seis cuentan con cercos perimetricos que de alguna manera impiden una mayor influencia antropica. La Tabla 6 muestra los 14 indices bioticos empleados en el estudio de ecosistemas acuaticos. Podemos notar comparando los indices de dominancia evaluados que existe una notable variaciòn. El valor promedio anual del indice de Simpson (C) y el de Berger-Parker (d) podrian estar diagnosticando buena calidad a las aguas evaluadas, en contraste, para el indice de Melntosh (D) estas son de pobre calidad. Tal incongruencia de resultados se observa tambien para el indice de Goodnight y Whittley, el cual valora a las aguas del rio Rimac como de buena condicion, Finalmente, evaluando el indice de Shannon-Wiener (H) (1,85), el rio Rimac presentaria aguas moderadamente contaminadas (Tabla 6). Segun Oliveira y Bispo (2001) las variaciones en la abundancia de los MIB pueden atribuirse a la disponibilidad de recursos, competencia, depredaciòn y el ciclo de vida. Sin embargo, en ambientes donde impera la degradacion ambiental, es importante considerar el efecto antrdpico del mismo en la comunidad biotica, la cual evidentemente determina la distribucion y estructura de los MIB (Riss et al. 2002; Figueroa et al. 2003; Goncalves y Aranha 2004).

De los cinco taxa mas abundantes como lo valora el N₂ de la serie de Hill, solo la familia Chironomidae obtuvo una correlaciòn positiva significativa con la BMWP' (r = 0,87; p = 0,02). Los otros cuatro taxa abundantes no presentaron correlaciones con la BMWP'. La DBO_5,20 encontro una correlacion negativa con BMWP' (r = -0,82; p - 0,04) y con H (r = -0.81; p = 0,04); y solo con el taxa Psychodidae una correlacion positiva (r = 0,98; p « 0.05). Los otros cuatro taxa abundantes no presentaron correlacion significativa con la DBO._5,20, Por otro lado, Dixidae y Oligochaeta presentaron una correlacion positiva (r = 0,88; p = 0,02). El aumento de la abundancia de Chironomidae resulto en un aumento en la BMWP'. Hecho que no coincide con lo referido por algunos autores (Figueroa et al. 2003) y en la afirmacion de que el aumento del grado de contaminacion organica produciria un incremento de Chironomidae. Este resultado incorporaria nuevas variables para nuevos estudios en ecosistemas de calidad critica a muy critica donde parece que la familia Chironomidae no podria ser empleada individualmente como un indicador adecuado de aguas degradadas (Figueroa et al. 2003).

En adicion, se nota la presencia de la familia Dixidae en un ecosistema netamente lotico, cuando esta familia se encuentra relacionada con ambientes lenticos, como lo afirman Merritt y Cummins (1996). Posiblemente han sido traidos por arrastre al rio Rimac o quiza estariamos tratando con una excepcion del comportamiento de la familia, lo cual tambien quedaria para ser dilucidado en proximos estudios.

Los dendrogramas cualitativos y cuantitativos construidos con base en los indices de diversidad beta evaluados para comparar las estaciones de muestreo se aprecian en las Figuras. 3 y 4. El analisis cualitativo agrupo dos sectores, el primero conformado por las estaciones cuatro, cinco y seis, y el segundo conformado por las estaciones uno, dos y tres; distinguiendose a si un sector alto y uno bajo, este ultimo mas comprometido con la actividad antropica. El analisis cuantitativo agrupo las estaciones en tres secrores diferenciables; el primer sector conformado por las estaciones tres y cuatro, de bastante similitud; el segundo sector conformado por las estaciones uno y cinco, de regular similitud; finalmente el sector conformado por las estaciones dos y seis, de poca similitud. La estacion dos presenta el mayor grado de degradacion ambiental y en contraste la estacion seis presentd el menor (Figs. 3-4).

El ACP produjo dos components con Eigenvalues > 1 para los indices bioticos (Tabla 7). Los indices bioticos de diversidad alfa fueron ordenados en un espacio de dos dimensiones, de acuerdo al ACP (Tabla 7; Fig. 5). Ocho indices bioticos (Numeros de Hill-N_1, McIntosh, Numeros de Hill-N₂, Simpson, Pielou, Berger-Parker, Shannon-Wiener, Margalef, Goodnight & Whittley, Numeros de Hill-₀) estuvieron correlacionados con el primer componente, que explico el 71,084 % de la variabilidad. Otros cuatro indices bioticos (ETP/ ETP + Chironomidae, ETP/Chironomidae/ Chironomidae en Diptera, BMWP' modif.) estuvieron correlacionados con el segundo componente, que explico el 20,323 % de la variabilidad (Tabla 7).

El estudio de los MIB como organismos indicadores de calidad de agua permite estimar si han sido afectados por cambios fisicos o quimicos de su habitat, pero no indica cual es la causa especifica que los afeta, lo cual debe ser estimado mediante asociacion con informacion del habitat fisico (e.g. substrato de fondo, velocidad de corrientes, profundidad) y con potenciales fuentes de estres (e.g. contaminantes acualicos) (Figueroa et al. 2003). Los resultados sugieren que el BMWP' es un buen indicador de la calidad de las aguas del rio Rimac. Sin embargo, la aparicion del BMWP' en otras cuencas debiera realizarse con precaucion a que pueden existir factores diferentes que pudieian afectar a las comunidades acuaticas. Al respecto se sugiere complementar el estudio de los MIB con adecuadas caracterizaciones de la columna de agua, ictiofauna, actividades productivas desarrolladas en la cuenca y evaluaciones experimentales de taxa locales frente a diferentes tipos de contaminantes (Figueroa et al. 2003: Paredes et al. 2004).

Conclusiones

- Se registraron 35 familias de macroinvertebrados bentonicos (MIB): Hexapoda (27), Annelidae (2), Mollusca (2), Arachnidae (2). Plathyhelmintes (I) y Chilopoda (1). De los 2.166 especimenes coleccionados Oligochaeta obtuvo la mayor abundancia absoluta; scguido por Psychodidae, Physidac, Chironomidae y Dixidac.

- El indice Biological Monitoring Working Party modificado (BMWP' modif.). valoro el agua del Rio Rimac, aguas abajo de la Atarjea, como de calidad critica o de aguas muy contaminadas (29 puntos). - La Demanda Bioquimica de Oxigeno (DBO_5,20), indico efectos de perturbacion en el ecosistema acuatico.

Agradecimientos

A Alvaro Torres, Maria Isabel La Torre. Rafael La Rosa, Alfonso Lizarraga, Carlos Santa Cruz, Ana Gutierrez. Margarita Robles, Augusto Mendoza. Edith Aguilar, Ana Carolina Paucar. Vicenta Morales. Jorge de la Cruz, Eric Carreño. Maria del Carmen Leon, Francisco Valdez. denia ,Salinas, Pedro Jose Paredes. Fabiola Isabel Correa, Gloria Asuncion Espinal, Victor Miguel Espinal, Adrian Ernesto Paredes y Oscar Ivan Garcia por su ayuda y apoyo en esta investigacion.

Literatura citada

ACOSTA, C.R.R. 2001. Patrones de diversidad espacio-temporal de insectos acuaticos y caracterizacion limnologica en la quebrada Canton. Tributaria del rio Rimac (Huarochiri, Lima). Resumenes. XLIII Convencion Nacional de Entomologia. Peru, p 75. Huancayo.        [ Links ]

APHA. 1995. Standard Methods for the examination of water and wastewater. Public Health Assoc. Washington. 19^th Ed. 8:1-21 p. 10:106-127p. 10:153-173p.        [ Links ]

BAPTISTA. D.F.. BUSS. D.F., DORVILLE. L.F., NESSIMIAN. J.L. 2001. Diversity and habitat preference of aquatic insects along longitudinal gradient of the Macae River basin. Rio de Janeiro. Brazil. Brazilian Journal of Biology 61: 249-258.        [ Links ]

BISPO. PC. FROEHLICH. C.G. OLIVEIRA. L.G. 2002. Stonefly (Plecoptera) fauna of streams in a mountainous area of Central Brazil: abiotic factors and nymph density. Revista BrasileiradeZoologia 19 (Supl. I): 325-334.        [ Links ]

BORROR. D., DE LONG. D.. THRIPLE-HORN, C. 1995. An introduction to the study of insects. Saunders College Publishing. USA. 6^th Ed. 827 p.        [ Links ]

BROWN. S.C., SMITH. K.. BATZER. D. 1997. Macroinvertebrate response to wetland restoration in northern New York. Environmental Entomology 26:1016-1024.        [ Links ]

COLE. G.A. 1998. Manual de Limnologia. Ed. Hemisferio Sur S.A. Argentina. 405 pp.        [ Links ]

DOMINGUEZ. E.. FERNANDEZ. H.R. 1998. Cal idad de los Ribs de la Cuenca del Sali (Tueuman. Argentina) medidoporun indice biotico. Serie Conservacion de la Naturaleza. Fundacion Miguel Lillo, Tueuman. Argentina. 38 p.        [ Links ]

FENOGLIO. S.. BADINO. G.. BONA. F. 2002. Benihic macroinvertebratc communities as indicators of river environment quality: an experience in Nicaragua. Revista biologia tropical 50: 1125-1131.        [ Links ]

FIGUEROA. R.. VALDOVINOS. C. ARAYA. E.. PARRA. O. 2003. Macroinvertebrados bentonicos como indicadores de calidad de agua de rios del sur de Chile. Revista Chilena de Historia Natural 76: 275-285.        [ Links ]

GALDEAN. N.. CALLISTO. M.. BARBOSA. F.A. 2001. Biodiversity assessment of benihic macroinvertebrates in altitudinal ecosystems of SerradoCipo (MG. Brazil). Brazilian Journal of Biology 61: 239-248.        [ Links ]

GONCALVES. F.B.. ARANHA. J.M.R. 2004. Spatial occupation and spatio-temporal by benthonic macroinvetebrates in Riberao river. Paranagua. PR (Brazil). Acta Biologica Parana. Curitiba 33: 181-191.        [ Links ]

HELLAWELL.J.M. 1978. Biological surveillance of rivers- (A biological monitoring handbook). Natural Environment Research Council. Water Research Center-Regional Water Authorities. England, p. 135-262.        [ Links ]

HURYN. A.D.. WALLACE. J.B. 2000. Life history and production of stream insects. Annual Review of Entomology45:83-1 10.        [ Links ]

IANNACONE. J.. ALVARINO. L.. MORENO. R.. REYES. M.. CHAUCA. J. 2000. Culicids (Dipslera) of the Chillon River and adjacent areas from Constitutional Province of Callao. Peru, during "El Nino" event 1997-1998. Acta Entomologica Chilena 24: 51-60.        [ Links ]

IANNACONE. J., MANSILLA J.. VENTURA. J. 2003. Macroinvertebrados en las lagunas de Puerto Viejo. Lima -Peru. EcologiaAplicada2: 116-124.        [ Links ]

LESLIE. H.A.. PAVLUK. T. DE VAATE. A.B., KRAAK. M.H. 1999. Triad assessment of the impact of chromium contamination on benthic macroinvertebrates in the Chusovaya river (Urals. Russia). Archives of Environmental Contamination and Toxicology 37: 182-189.        [ Links ]

LIND. O.T. 1979. Handbook of common methods in Limnology. C.V. Mosby Company. U.K. 199 p.        [ Links ]

MARQUES. M.J., MARTINEZ-CONDE, E.. ROVIRA. J.V. 2001. Los macroinvertebrados como indices de evaluacion rapida de ecosistemas acuaticos contaminados por metales pesados. Ecotoxicology and Environmental Restoration 4: 25-31.        [ Links ]

METZELING. L.. CHESSMAN. B.. HARD-WICK. R.. WONG. V. 2003. Rapid assessment of rivers using macroinvertebrates: the role of experience, and comparison with quantitative methods. Hydrobiologia 510:39-52. MERRITT. R.W.. CUMMINS. K.M. 1996. An introduction to the aquatic insects of North America. Kendall Hunt Publishing Co. Dubuque. IA. U.S.A. 3^th. Ed. XIII + 862 p.        [ Links ]         [ Links ]

MORENO. C. E. 2001. Metodos para medir la biodiversidad. Manuales y Tesis SEA. Sociedad Entomologica Aragonesa Ed. Madrid. España. 80 p.        [ Links ]

OGBEIBU. A.E.. ORIBHABOR. B.J. 2002. Ecological impact of river impoundment using benthic macroinvertebrates as indicators. Water Research 36: 2427-2436.        [ Links ]

OLIVEIRA. L.G.. BISPO. PC. 2001. Ecologia de comunidadaes das larvas de Trichoptera Kirby (Insecla) em dois corregos de primeira orden da Serra dos Pireneus. Pirenopolis: Goias. Brazil. Revista BrasileiraZoologials: 1245-1252.        [ Links ]

PAREDES. C. IANNACONE, J.. ALVARINO. L. 2004. Macroinvertebrados bentonicos como indicadores biologicos de la calidad de agua en dos rios de Cajamarca y Amazonas. Peru. Revista Peruana de Entomologia 44: 107-1 18.        [ Links ]

RISS. W.. OSPINA. R.. GUTIERREZ. J.D. 2002. Establecimientos de valores de bioindicacitin para macroinvertebrados acuaticos de la sabana de Bogota. Caldasia 24: 135-156.        [ Links ]

ROGERS. C.E.. BRABANDER. D.J.. BAR-BOUR. M.T.. HEMOND. H.F. 2002. Use of physical, chemical, and biological indices to assess impacts of contaminants and physical habitat alteration in urban streams. Environmental Toxicology and Chemistry 21: 1156-1167.        [ Links ]

STEHR.W. F. 1987. Immature insects. Kendal/ Hunt Publishing Company. Dubuque. Iowa. US. 754 p.        [ Links ]

THIENGO. S.C.. FERNANDEZ, M.A.. BOA-VENTURA. M.F.. GRAULT C.E.. SILVA. H.F.R.. MATTOS. A.C. SANTOS S. B. 2001. Freshwater snails and schistosomiasis mansoni in the state of Rio de Janeiro. Brazil: l-Metropolitan Mesoregion. Memo-rias do Instituto Oswaldo Cruz. Rio de Janeiro 96: 177-184.        [ Links ]

VIVANCO. M. 1999. Analisis estadistico multivariable. Teoria y practica. Editorial Universitaria. Santiago de Chile. 234 p.        [ Links ]

WASHINGTON, H.G. 1984. Diversity, Biotic and similarity indices: A review with special relevance to aquatic ecosystems. Water Research 18: 653-694.        [ Links ]

WETZEL. R.G.. LIKENS. G.E. 2000. Limno-logical analyses. Springer - Verlag. New York. 3,hEd. pp. 424-429.        [ Links ]

ZAR J. H. 1996. Bioestatistical analysis. 3^th Ed. Prentice-Hall. Inc. Upper Saddle River. New Jersey. 662 p.        [ Links ]

Recibido: 01 -mar-05 - Aceptado: 11 -jul-05