Tesis Rodrigo Meza L
Tesis Rodrigo Meza L
Tesis Rodrigo Meza L
PROFESOR GUÍA:
XIMENA VARGAS MESA
MIEMBROS DE LA COMISIÓN:
YARKO NIÑO CAMPOS
GINGER MARTÍNEZ ZAROR
SANTIAGO DE CHILE
OCTUBRE 2009
RESUMEN
La zona de estudio seleccionada para esta tesis fue un tramo del río Huequecura,
afluente al río Biobío. Como muchas cuencas pluviales del centro y sur de Chile,
carece de fluviometría, lo cual condujo a resolver la simulación de la calidad del hábitat
de peces nativos empleando el modelo de precipitación-escorrentía GR4J acoplado al
modelo hidráulico HEC-RAS y al modelo de simulación de calidad de hábitat CASiMiR,
el cual es calificado como un sistema experto, al emplear la información cualitativa de
un biólogo en su implementación. La información de peces utilizada corresponde a
aquella levantada por estudios de impacto y seguimiento ambiental en la zona, siendo
ésta la mejor fuente de datos aquí y en muchos ríos del país. Todos estos elementos
representan una innovación en la manera de enfrentar el problema de evaluar el
impacto del caudal sobre los ecosistemas fluviales, ya que se toma un caso real y se
enfrenta desde un enfoque ampliamente validado por la comunidad internacional, pero
con los recursos que Chile posee.
Los resultados indican, en el caso de los bagres, que existe una marcada idoneidad del
hábitat para los individuos juveniles por sobre los adultos, lo cual explica en parte la
mayor presencia de juveniles en las capturas empleadas para analizar el modelo. Sin
embargo, el arte de pesca solo permite acceder a las orillas y el modelo solo permite
diferenciar zonas a través de los mapas, limitando la posibilidad de validar a través de
una sectorización genérica, dado que las capturas analizadas no indican el punto
exacto y solo hacen referencia de su ejecución a la orilla del tramo como un todo.
La investigación indicó que los esfuerzos por simular la calidad de hábitat de peces
nativos podrían estar fuertemente limitados por una política ambiental chilena poco
clara, incluyendo falta de datos en todas las áreas esenciales que un estudio medio
ambiental requiere. El estudio concluye con recomendaciones para mejorar la
integración de políticas ambientales en el manejo de recursos hídricos de Chile, por
ejemplo, modificar el sistema político-legal del país para permitir un manejo flexible y
eficaz del agua, dado el escenario de escasez futura que plantean los estudios más
recientes. Además, se sugiere una metodología alternativa para estimar el caudal
ecológico de los ríos, basada en los índices de calidad de hábitat.
AGRADECIMIENTOS
Nos es fácil para mi comenzar estos agradecimientos, porque siento que cada persona
que está y ha estado en mi vida tiene una cuota de responsabilidad en lo que he hecho
y en el camino que he decidido seguir en mi vida personal y profesional. Siento que de
una u otro forma todo se conjuga en este trabajo, en el cual comprometí gran parte de
lo que soy y de lo que deseo para mi y para los demás. Sé que el camino que debemos
transitar en esta vida es largo, enrevesado e incierto, pero la materialización de este
trabajo significa para mí una suerte de faro que estará siempre ahí para orientarme, de
ahí su significado para mí.
Por último, quiero agradecer a mi mujer, Daniela Domínguez, por todo el amor que me
entrega día a día, el apoyo incondicional en las buenas y las malas y por su fe en mí.
Sin ti, muchos de mis logros no tendrían sentido ni hubiesen sido posibles.
1
Respecto a mi tesis, quiero agradecer a mi comisión de tesis por sus comentarios y
críticas, por su dedicación y esmero para que ésta sea lo mejor posible.
También quiero dar un especial agradecimiento al Dr. Manuel Contreras, del Centro de
Ecología Aplicada, por sus comentarios y consejos, por facilitar sus instalaciones,
personal e información en la realización de esta tesis.
2
Lejos de ser una carga, el desarrollo sostenible supone
una oportunidad excepcional: desde el punto de vista
económico, para crear mercados y empleos; desde el
punto de vista social, para integrar a los marginados; y
desde el punto de vista político, para que todos los
hombres y mujeres tengan voz y voto al decidir su propio
futuro.
3
ÍNDICE DE CONTENIDOS
1. INTRODUCCIÓN ________________________________________________12
1.1 MOTIVACIÓN____________________________________________________ 12
1.2 OBJETIVO GENERAL_____________________________________________ 14
1.2.1 Objetivos Específicos _____________________________________________ 14
1.3 ORGANIZACIÓN DEL INFORME __________________________________ 15
2. REVISIÓN BIBLIOGRÁFICA ______________________________________18
2.1 MODELACIÓN DE HÁBITATS ACUÁTICOS_________________________ 20
2.1.1 PHABSIM ______________________________________________________ 21
2.1.2 CASIMIR _______________________________________________________ 23
2.1.2.1 APLICACIONES ___________________________________________________ 27
2.2 MODELACIÓN HIDROLÓGICA ____________________________________ 31
2.2.1 APLICACIÓN DEL MODELO GR4J EN CHILE CENTRAL_________________ 36
2.2.1.1 PARÁMETROS E INDICES GEOMORFOLÓGICOS ______________________ 36
2.3 CHILE EN EL CONTEXTO DE UN MANEJO INTEGRADO DE RECURSOS
HÍDRICOS ______________________________________________________________ 39
2.3.1 LEGISLACIÓN __________________________________________________ 39
2.3.2 SUFICIENCIA DE LA INFORMACIÓN DISPONIBLE _____________________ 45
2.3.2.1 HIDROMETRÍA EN EL CONTEXTO ACTUAL__________________________ 46
2.3.2.2 INVESTIGACIÓN BIOLÓGICA ______________________________________ 48
2.3.2.3 HERRAMIENTAS DE SELECCIÓN DE LA UNIDAD DE ESTUDIO_________ 49
2.3.3 ESPECIES DULCEACUÍCOLAS CHILENAS ___________________________ 51
2.3.3.1 PERCILIA IRWINI _________________________________________________ 52
2.3.3.2 DIPLOMYSTES NAHUELBUTAENSIS ________________________________ 54
2.3.3.3 TRICHOMYCTERUS AREOLATUS ___________________________________ 55
2.3.3.4 TRICHOMYCTERUS CHILTONI _____________________________________ 57
2.3.3.5 BULLOCKIA MALDONADOI________________________________________ 59
2.3.3.6 BASILICHTHYS AUSTRALIS________________________________________ 60
2.4 ESTADO DEL ARTE EN LA PROBLEMÁTICA DE LA GESTIÓN DE
RECURSOS HÍDRICOS CON ENFOQUE MULTIDISCIPLINARIO ____________ 62
2.4.1 FORMULACIÓN DE MARCOS CONCEPTUALES _______________________ 63
2.4.2 PROYECTOS DE ESTANDARIZACIÓN E INTEGRACIÓN A MACROESCALA 71
2.4.3 SISTEMAS EXPERTOS ___________________________________________ 73
2.4.4 MODELOS ENLAZADOS __________________________________________ 75
2.5 COMENTARIOS __________________________________________________ 78
3. CARACTERIZACIÓN ZONA DE ESTUDIO _________________________83
3.1 UBICACIÓN DE LA ZONA DE ESTUDIO ____________________________ 83
3.1.1 CARACTERÍSTICAS DEL TRAMO ESTUDIADO ________________________ 84
3.2 CLIMA __________________________________________________________ 85
3.2.1 PRECIPITACIONES ______________________________________________ 85
3.2.2 TEMPERATURAS MEDIAS ANUALES _______________________________ 87
4
3.2.3 EVAPORACIÓN MEDIA ANUAL_____________________________________ 87
3.2.4 RED DE ESTACIONES HIDROLÓGICAS _____________________________ 88
3.3 CALIDAD DEL AGUA _____________________________________________ 88
3.4 USO DE SUELOS__________________________________________________ 91
3.5 AGUAS SUBTERRÁNEAS__________________________________________ 93
3.6 COMENTARIOS __________________________________________________ 93
4. MODELACIÓN HIDROLÓGICA ___________________________________96
4.1 PARÁMETROS GEOMORFOLÓGICOS _____________________________ 96
4.2 ÍNDICES GEOMORFOLÓGICOS ___________________________________ 99
4.3 ESTADÍSTICA Y MEDICIONES ___________________________________ 102
4.4 MODELACIÓN Y RESULTADOS __________________________________ 105
4.5 COMENTARIOS _________________________________________________ 111
5. MODELACIÓN HIDRÁULICA ____________________________________114
5.1 CALIBRACIÓN MODELO HEC-RAS _______________________________ 114
5.1.1 CONSIDERACIONES ____________________________________________ 114
5.1.1.1 MÉTODO DE APROXIMACIÓN _____________________________________ 119
5.1.1.2 MÉTODO DE ESTIMACIÓN DE LA RUGOSIDAD______________________ 121
5.2 MODELACIÓN DEL PERÍODO DE ESTUDIO _______________________ 125
5.3 COMENTARIOS _________________________________________________ 127
6. MODELACIÓN DE LA CALIDAD DE HÁBITAT ____________________129
6.1 RECOLECCIÓN Y PROCESAMIENTO DE LOS DATOS ______________ 129
6.1.1 REGLAS DIFUSAS ______________________________________________ 130
6.1.2 CONJUNTOS DIFUSOS __________________________________________ 131
6.1.2.1 AGRUPAMIENTO DIFUSO _________________________________________ 134
6.2 RESULTADOS GENERALES ______________________________________ 142
6.2.1 PERCILIA IRWINI _______________________________________________ 143
6.2.2 DIPLOMYSTES NAHUELBUTAENSIS _______________________________ 147
6.2.3 TRICHOMYCTERUS AREOLATUS _________________________________ 151
6.2.4 TRICHOMYCTERUS CHILTONI____________________________________ 155
6.2.5 BULLOCKIA MALDONADOI _______________________________________ 159
6.2.6 BASILICHTHYS AUSTRALIS ______________________________________ 163
6.3 VALORACIÓN EXPERTA DE LOS RESULTADOS___________________ 167
6.4 COMENTARIOS _________________________________________________ 169
7. ANÁLISIS Y DISCUSIÓN DE RESULTADOS _______________________173
7.1 INFORMACIÓN DE CAPTURAS EN EL RÍO HUEQUECURA _________ 173
7.2 EFECTO DE LA VARIABILIDAD FLUVIAL SOBRE LA ABUNDANCIA DE
LA FAUNA ÍCTICA EN EL RÍO HUEQUECURA ___________________________ 180
7.2.1 T. AREOLATUS ________________________________________________ 181
5
7.2.2 D. NAHUELBUTAENSIS __________________________________________ 186
7.2.3 P. IRWINI _____________________________________________________ 190
7.3 DEFINIENDO UN CAUDAL ECOLÓGICO __________________________ 196
7.3.1 RÉGIMEN SUSTENTABLE DE CAUDALES __________________________ 199
7.4 COMENTARIOS _________________________________________________ 203
8. CONCLUSIONES________________________________________________206
8.1 VISIÓN MULTIDISCIPLINARIA EN EL MANEJO DE LOS RECURSOS
HÍDRICOS _____________________________________________________________ 207
8.2 MODELOS MATEMÁTICOS EN UN CONTEXTO DE ESCASEZ DE
INFORMACIÓN ________________________________________________________ 209
8.3 OBJETO Y MÉTODO: REPLICABILIDAD DEL ESTUDIO ____________ 212
8.4 RÉGIMEN SUSTENTABLE DE CAUDALES _________________________ 214
REFERENCIAS _____________________________________________________217
REFERENCIAS WEB ________________________________________________230
6
ÍNDICE DE CUADROS Y FIGURAS
CUADROS
7
Cuadro 7.4: Muestra de P. irwini capturados entre el 2000 y 2005............................. 178
Cuadro 7.5: Resumen de la calidad de hábitat (HHS), registro de capturas y frecuencia
de caudales (Adaptado de Meza & Vargas, 2009)...................................................... 181
FIGURAS
8
Figura 4.4: Relación entre IG y x1 para Módulo. La proyección de los datos extremos en
el eje X indica el rango definido de cada IG (Fuente: Salvatierra, 2008) .................... 102
Figura 4.5: Contraste del caudal pasante en el tramo de estudio. Arriba el río
Huequecura en marzo de 2007, cuando se realizó aforo exitoso. Abajo el río en
septiembre de 2007, cuando no se pudo aforar.......................................................... 104
Figura 4.6: Ubicación de cuenca respecto a curvas de isoevaporación (en verde).
Estación Quilaco fue ubicada geográficamente, pues código BNA no coincide con lo
señalado por DGA (1987). .......................................................................................... 105
Figura 4.7: Resultado de la modelación preliminar, en el período de ajuste, para los
parámetros de Nash-Sutcliffe y Módulo. ..................................................................... 110
Figura 4.8: Contraste del resultado final con el preliminar, en el período de ajuste, para
los parámetros de Módulo. .......................................................................................... 110
Figura 4.9: Caudales medios diarios en el período de la modelación 2002-2004. Se
compara el resultado preliminar de Módulo con la serie definitiva, producto de los
parámetros ajustados.................................................................................................. 111
Figura 5.1: Configuración del río de acuerdo a los requerimientos de HEC-RAS. Flecha
indica el sentido del flujo. Coordenadas relativas en metros. ..................................... 115
Figura 5.2: Eje hidráulico medido del tramo en estudio. ............................................. 116
Figura 5.3: Calibración modelo hidráulico. Perfil modelado versus el medido. ........... 124
Figura 5.4: Resultado gráfico calibración HEC-RAS. .................................................. 125
Figura 5.5: Curvas de descarga por perfil. Para todo el rango de caudales a modelar.
.................................................................................................................................... 126
Figura 6.1: Medición de variables ............................................................................... 132
Figura 6.2: Resultado del agrupamiento difuso sobre los datos de la campaña
efectuada en septiembre de 2007. .............................................................................. 135
Figura 6.3: Aproximación lineal hecha sobre el resultado del agrupamiento difuso ... 136
Figura 6.4: Conjuntos difusos definitivos de la etapa uno en contraste con valoración
realizado por el experto en la campaña de septiembre de 2007................................. 138
Figura 6.5: Conjuntos difusos modificados de acuerdo a la apreciación del experto
sobre las variables medidas en cada punto. ............................................................... 141
Figura 6.6: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie: Percilia
irwini. Estado de vida: Adulto. ..................................................................................... 144
Figura 6.7: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8
[m3/s]. Especie: Percilia irwini. Estado de vida: Adulto................................................ 144
Figura 6.8: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie: Percilia
irwini. Estado de vida: Juvenil. .................................................................................... 146
Figura 6.9: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8
[m3/s]. Especie: Percilia irwini. Estado de vida: Juvenil............................................... 146
Figura 6.10: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Diplomystes nahuelbutaensis. Estado de vida: Adulto................................................ 148
Figura 6.11: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Diplomystes nahuelbutaensis. Estado de vida: Adulto. ............... 148
Figura 6.12: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Diplomystes nahuelbutaensis. Estado de vida: Juvenil............................................... 150
9
Figura 6.13: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Diplomystes nahuelbutaensis. Estado de vida: Juvenil. .............. 150
Figura 6.14: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Trichomycterus areolatus. Estado de vida: Adulto. ..................................................... 152
Figura 6.15: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Trichomycterus areolatus. Estado de vida: Adulto. ..................... 152
Figura 6.16: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Trichomycterus areolatus. Estado de vida: Juvenil. .................................................... 154
Figura 6.17: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Trichomycterus areolatus. Estado de vida: Juvenil. .................... 154
Figura 6.18: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Trichomycterus chiltoni. Estado de vida: Adulto.......................................................... 156
Figura 6.19: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Trichomycterus chiltoni. Estado de vida: Adulto. ......................... 156
Figura 6.20: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Trichomycterus chiltoni. Estado de vida: Juvenil......................................................... 158
Figura 6.21: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Trichomycterus chiltoni. Estado de vida: Juvenil. ........................ 158
Figura 6.22: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Bullockia maldonadoi. Estado de vida: Adulto............................................................. 160
Figura 6.23: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Bullockia maldonadoi. Estado de vida: Adulto............................. 160
Figura 6.24: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Bullockia maldonadoi. Estado de vida: Juvenil............................................................ 162
Figura 6.25: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Bullockia maldonadoi. Estado de vida: Juvenil............................ 162
Figura 6.26: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Basilichthys australis. Estado de vida: Adulto. ............................................................ 164
Figura 6.27: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Basilichthys australis. Estado de vida: Adulto. ............................ 164
Figura 6.28: Evolución de los índices WUA y HHS con el caudal (izquierda) y la
evolución de la calidad de hábitat disponible con el caudal (derecha). Especie:
Basilichthys australis. Estado de vida: Juvenil. ........................................................... 166
Figura 6.29: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de
0.8 [m3/s]. Especie: Basilichthys australis. Estado de vida: Juvenil. ........................... 166
Figura 7.1: Distribución de frecuencias normalizadas por talla de P. irwini, en el período
2000-2005. .................................................................................................................. 179
Figura 7.2: Valores modales identificados para P. irwini............................................. 179
Figura 7.3: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: T. areolatus. Estado de vida: Adulto.
.................................................................................................................................... 182
10
Figura 7.4: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: T. areolatus. Estado de vida: Juvenil.
.................................................................................................................................... 183
Figura 7.5: Presencia detectada de T. areolatus en relación a la calidad media de
hábitat en el tramo de estudio. .................................................................................... 184
Figura 7.6: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: D. nahuelbutaensis. Estado de vida:
Adulto. ......................................................................................................................... 187
Figura 7.7: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: D. nahuelbutaensis. Estado de vida:
Juvenil. ........................................................................................................................ 188
Figura 7.8: Presencia detectada de D. nahuelbutaensis en relación a la calidad media
de hábitat en el tramo de estudio. ............................................................................... 189
Figura 7.9: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: P. irwini. Estado de vida: Adulto. .... 191
Figura 7.10: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la
zona de estudio. Período 2002 - 2004. Especie: P. irwini. Estado de vida: Juvenil. ... 192
Figura 7.11: Presencia detectada de P. irwini en relación a la calidad media de hábitat
en el tramo de estudio................................................................................................. 193
Figura 7.12: Relación entre la proporción de de juveniles sobre todas las capturas de la
especie y la modelación de CASiMiR, a través del HHS. ........................................... 195
11
“…la naturaleza sólo se ofrece a través de la reja de las denominaciones..."
Michel Foucault
1. INTRODUCCIÓN
1.1 MOTIVACIÓN
Hoy en día las iniciativas relacionadas con el agua son incontables y se extienden por
todo el mundo, ya no solo como un problema en torno al hombre, sino como un desafío
que debe ser abordado en forma multisectorial y multidisciplinaria, considerando al
hombre como parte de un sistema natural mayor.
Los ríos han sido uno de los focos de estudio en la búsqueda por comprender y
manejar los recursos hídricos de agua dulce. Debido a la canalización y construcción
de estructuras de protección de crecidas, así como la construcción de embalses para la
producción de energía o para la irrigación, se ha producido un impacto algunas veces
irreversible en muchos cauces naturales alrededor del mundo (Harby et al., 2004).
También la urbanización, el cambio de uso de los suelos y la explotación de acuíferos
han sido factores importantes en el cambio de los regímenes hidrológicos, todo lo cual
ha derivado en impactos negativos sobre los ecosistemas fluviales (Harby et al., 2004).
En esta línea se han dedicado gran cantidad de recursos y esfuerzos, tanto en ciencia
básica como aplicada y desde la biología de las comunidades acuáticas hasta modelos
de simulación de hábitats. El objetivo, bajo un escenario de constante crecimiento de la
población y cambio climático local y global, es lograr un manejo sustentable de los
recursos hídricos, ya sea desde un enfoque económico o desde el punto de vista
social-ambiental. Hoy en día existe un amplio consenso en estas materias, lo cual ha
llevado a diversos organismos, como el Banco Mundial (Hirji y Panella, 2003), a crear
plataformas de apoyo a países emergentes y a fomentar estudios y políticas dentro de
aquellos países que se esfuercen en el manejo sustentable, integrado y adaptativo de
los recursos hídricos. En este contexto, la Comisión sobre el Desarrollo Sostenible
(CDS) declaró en el año 2002 que “erradicar la pobreza, cambiar los patrones de
12
producción y consumo insostenibles y proteger y administrar los recursos naturales del
desarrollo social y económico constituyen los objetivos primordiales y la exigencia
esencial de un desarrollo sostenible” (WWAP, 2003).
13
1.2 OBJETIVO GENERAL
Desarrollar un modelo que simule la calidad de hábitat disponible para peces nativos
chilenos en un tramo del río Huequecura, en la Región del Bio Bío, y analizar el
desarrollo de la modelación en el contexto chileno, a través de elaborar un análisis
crítico de los resultados y de las condiciones técnicas y operacionales requeridas por el
modelo de simulación de calidad de hábitat CASiMiR y diferenciar éstas con aquellas
disponibles en la zona de estudio y el país.
14
1.3 ORGANIZACIÓN DEL INFORME
15
subterráneas permite asegurar la consistencia del modelo hidrológico mientras que la
red de estaciones hidrológicas asegura que la información empleada es la de mejor
calidad disponible. La calidad del agua asegura que no existe un elemento, además de
las condiciones físicas del escurrimiento, que controlen la disponibilidad de hábitat en
la zona. El uso de suelos permite confirmar los resultados de calidad de agua y entrega
antecedentes sobre el grado de intervención de la cuenca, lo cual también es factor
forzante sobre la calidad de hábitat de sus cauces.
16
El Capítulo 8 presenta las conclusiones finales de esta investigación, agrupadas bajo
cuatro líneas de análisis surgidas del proceso de modelación, obtención de resultados
y análisis de antecedentes bibliográficos. No se busca resumir aquí conclusiones de
cada uno de los Capítulos desarrollados, pues en cada uno se incluyó un apartado con
comentarios críticos que resumen lo más relevante en cada caso. Este Capítulo busca
integrar el desarrollo del informe a los objetivos propuestos, con una revisión
retrospectiva de los procesos, resultados y análisis obtenidos en los Capítulos previos.
En este proceso se busca, esencialmente, dilucidar de qué forma es posible aplicar
este modelo en Chile y cuales son los requerimientos para replicar su uso en otras
cuencas.
17
2. REVISIÓN BIBLIOGRÁFICA
Surgen preguntas básicas, como ¿qué variables son relevantes?, ¿cuánta complejidad
es necesaria?, ¿qué funciones de relación serán las apropiadas? Además, existe la
problemática de la información, pues cualquier modelo que intente concentrar varios
enfoques requerirá información bastante específica. Esto significa costos y dificultades
que tal vez, dada la complejidad del problema, no justifiquen su aplicación.
18
determina el carácter de las instituciones, leyes, políticas y aspectos sociales
relacionados con el recurso. Como se verá, no es trivial la forma en que un país
considera el recurso agua al momento de buscar un manejo integrado y adaptativo de
los recursos. Si se considera el trabajo de Rogers (2006) se podrá visualizar que la
solución del manejo de recursos atraviesa a toda la sociedad, sus instituciones y las
relaciones que entre ellos se establecen, lo cual estará determinado, de acuerdo a
Bauer (2005) en gran parte por los conceptos que definen al recurso en cuestión.
El uso de herramientas científicas integradoras se ha vuelto cada vez más común. Por
ejemplo, la valoración del estado ecológico de aguas superficiales, principalmente ríos,
y el modelamiento de hábitat físicos, como una herramienta de enfoque
interdisciplinario, se ha vuelto un método frecuente, especialmente enfocado a los
peces (Schneider et al., 2004). Estos son utilizados como indicadores para analizar la
problemática de los caudales ecológicos, pues analizar el efecto de la variable caudal
en todo el ecosistema resulta muy difícil e innecesario (Silva et al., 1994). Otras
razones para el uso de peces como indicadores son, primero, que existen especies de
interés (nativas, en peligro de extinción, etc.) y segundo, que los peces son
relativamente fáciles de estudiar y sus preferencias, en cuanto a sus hábitats físicos,
son comparativamente mejor comprendidas y documentadas que otras especies
acuáticas de ríos (Schneider et al., 2004).
El presente capítulo busca hacer una revisión extensa sobre una amplia variedad de
temas relacionados no solo con la modelación de hábitats, sino que también con la
discusión en torno a la gestión de recursos hídricos en la cual se inserta y forma parte
esta clase de herramientas.
19
2.1 MODELACIÓN DE HÁBITATS ACUÁTICOS
20
2.1.1 PHABSIM
A fines de los años setenta, el U.S. Fish and Wildlife Service desarrolló la Metodología
Incremental para la Determinación de Caudales Mínimos Aconsejables IFIM (Instream
Flow Incremental Methodology), la cual incluye un sistema de simulación de hábitat de
tipo modular llamado PHABSIM. Éste está compuesto por una librería de modelos de
simulación interconectados (Espinoza et al., 1999) que permiten simular los
requerimientos físicos de hábitat de una especie en particular.
El enfoque de PHABSIM ha sido utilizado probablemente en más del 90% de todos los
esfuerzos mundiales por modelar hábitat de peces (Jorde, 2006). En Chile se aplicó
concretamente por primera vez el año 2000, en el río Laja (Parra et al., 2000).
Aunque existen adaptaciones de PHABSIM para ser utilizado con modelos hidráulicos
de una a tres dimensiones, el procedimiento de aplicación es esencialmente el mismo
en todos los casos y consta de 4 partes: recolección de los datos de terreno,
simulación hidráulica, construcción y selección de las curvas de preferencia y
simulación del hábitat físico (Diez Hernández, 2005). La Figura 2.1 muestra un
esquema ampliado de las partes señaladas anteriormente.
21
El tipo de información que requiere el modelo de simulación puede ser agrupada en
dos clases: información hidráulica e información biológica (Parra et al., 2000).
Aunque las variables que permiten definir un hábitat acuático son varias, a escala local
en un determinado tramo son clave las variables de velocidad, profundidad y tipo de
sustrato de fondo. Únicamente con ellas basta para determinar los hábitats potenciales
para una especie dada (Parra et al., 2000).
22
recolectada en terreno (se puede sumar el conocimiento experto), pues en cada punto
donde se pesca un ejemplar, en el caso de los peces, se miden las tres variables
mencionadas: velocidad, profundidad y tipo de sustrato, pudiendo así obtener las
curvas de preferencia.
Dicho índice global se denomina área ponderada utilizable (APU), la cual recibe la
denominación genérica “microhábitat” (Parra et al., 2000; Universidad de Chile, 1998).
Un supuesto elemental del índice global de habitabilidad generado mediante esta
metodología es que a cada velocidad, profundidad y tipo de sustrato existe un único
valor de preferencia, lo cual señala el carácter determinista que posee PHABSIM.
Una lista bastante amplia de aplicaciones de este software puede ser encontrada en
Universidad de Chile (1998) y en Diez Hernández (2005).
2.1.2 CASIMIR
23
en el río Huequecura, fue desarrollado por el Instituto de Ingeniería Hidráulica de la
Universidad de Stuttgart (Schneider et al., 2004) y está basado en PHABSIM (Diez
Hernández, 2005), como la mayoría de los modelos de su tipo.
Los requerimientos de información del software son básicamente los mismos que en
los demás modelos de su clase, pero tiene un enfoque nuevo que lo hace más versátil
y que permite aprovechar el conocimiento de los expertos, el cual muchas veces está
poco documentado y es el de mayor grado de certeza disponible. En Chile, esto es
particularmente cierto.
Para la parte biológica, CASiMiR utiliza un enfoque basado en reglas difusas (Fuzzy
Rule-Based), la cual permite el uso de “verdades parciales” o, en otras palabras, que
tienen cierta pertenencia al conjunto “verdad”. El uso de este enfoque está sustentado
por la teoría de conjuntos difusos (Fuzzy Sets Theory), la cual es una extensión de la
teoría clásica de conjuntos (Jorde, 2006; Schneider et al., 2004).
Una de las principales ventajas de los conjuntos difusos es que permiten la definición
de variables lingüísticas, razón por la cual son capaces de integrar el conocimiento
24
semi cuantitativo de los expertos a la modelación numérica (Scheneider et al., 2004).
La conexión entre la definición imprecisa y verbal del experto y los valores numéricos
de las variables, indicadoras de cierta calidad de hábitat, es posible lograrla gracias a
la lógica difusa (Fuzzy Logic) que está involucrada en los conjuntos difusos. Aun más,
recoge el hecho de que no existe en la naturaleza un valor único de velocidad,
profundidad y tipo de sustrato que implique un valor fijo y único de habitabilidad o de
preferencia de hábitat, ya que la complejidad del sistema está siendo simplificada y
dichas variables solo permiten estimar aproximadamente cuáles serán las preferencias
de los peces, dada cierta habitabilidad calculada.
Por estas razones CASiMiR puede ser catalogado como un sistema experto, el cual, de
acuerdo a Puppe (1993) es aquel que permite reconstruir la experiencia y la capacidad
de razonamiento de especialistas calificados dentro de su dominio de aplicación.
Según este mismo autor, los dominios de estos sistemas son considerablemente
difusos, es decir, de aplicación a problemas que no son claradamente entendidos, lo
cual se debe a que el conocimiento y los datos en los dominios de los sistemas
expertos son usualmente inciertos, incompletos y dependientes del tiempo.
Será el experto quien definirá las reglas difusas (Fuzzy Rules) en primera instancia,
pudiendo éstas ser ajustadas posteriormente con información de terreno o de la
literatura disponible. Adicionalmente, en terreno se medirán las variables de velocidad,
profundidad y tipo de sustrato en conjunto con el experto, de tal modo de obtener los
conjuntos difusos (Fuzzy Sets) que relacionarán las variables lingüísticas con valores
numéricos, es decir, será en terreno donde, por ejemplo, la afirmación “velocidad alta”
encontrará la valoración numérica que se necesita para la modelación.
25
puedan ser extraídos desde la literatura disponible, simplificando enormemente la
modelación en su etapa inicial.
Figura 2.2: Esquema de las etapas necesarias en el uso de CASiMiR (Adaptado de Schneider et al.,
2004)
En resumen:
26
CASiMiR es un modelo de simulación de hábitat que utiliza la lógica difusa para
incorporar la incerteza inherente a los datos que se obtienen de monitorear un sistema,
en este caso, las variables ecológicas de un río (Adriaenssens et al., 2006).
Para lo anterior utiliza el conocimiento experto para construir Fuzzy Sets y Fuzzy
Rules, con los cuales es capaz de modelar la disponibilidad de hábitat en un sistema
fluvial (Schneider et al., 2002; Ahmadi-Nedushan et al., 2006; Adriaenssens et al.,
2006)
2.1.2.1 APLICACIONES
Existen diversas aplicaciones en las cuales se puede utilizar CASiMiR, por ejemplo, en
la evaluación de alternativas para caudales ecológicos o escenarios de restauración de
cauces. El despliegue gráfico permite observar qué zonas de los tramos estudiados se
ven afectadas en cada escenario o si las medidas de reparación del cauce darán los
resultados esperados.
Figura 2.3: Despliegue gráfico de las variables física introducidas en CASiMiR (Adaptación de
Schneider et al., 2002)
27
En la Figura 2.3 cada celda contiene un color que representa un valor, para cada
variable, los cuales variarán según el caudal pasante. En un enfoque tradicional la
modelación del hábitat físico se hace utilizando curvas de preferencia univariadas,
donde cada parámetro es considerado independiente de los demás y son unidos
únicamente por el procedimiento de la modelación de hábitat (Schneider et al., 2002).
Figura 2.4: Forma en que CASiMiR determina la consecuencia de cada regla aplicada en las celdas
del tramo modelado (Adaptado de Jorde (2006)).
Los atributos físicos de cada celda son comparados con los conjuntos difusos,
verificando el grado de cumplimiento (DOF: degree of fulfillment) de la(s) regla(s)
28
asociada(s) a dicha celda (Jorde, 2006). En la Figura 2.4 se observa que la velocidad
permite dos reglas en la misma celda, ya que tiene un grado de pertenencia tanto al
conjunto “Alta” como “Media”. Para determinar el DOF, se utiliza un operador “y” (el
equivalente matemático al operador lógico “y”).
En cada situación (celda) existirá un cierto número de reglas que son aplicables.
Luego, se podrá utilizar un variado número de métodos para generar consecuencias
compuestas de las diferentes reglas aplicables. CASiMiR utiliza un método de suma
normalizada y ponderada para combinar los DOF que tiende a minimizar la influencia
de los más pequeños (Jorde, 2006).
El SI se calcula para cada celda y permite obtener los mapas de idoneidad de hábitat
para cada una de las especies y estados de vida estudiados. Estos mapas son
similares a los de la Figura 2.3 y se incluyen en el Capítulo 6 para el tramo de estudio
de esta tesis.
Por otra parte, CASiMiR permite expresar a través de ciertos índices un valor integral
de la calidad de hábitat. Este permite asignar cierta calidad a todo el tramo, asociada a
1
Adaptación de un término (defuzzyfication) que no tiene traducción al español y que se refiere a la transformación de
variables lingüísticas en valores no difusos (crispy numbers).
29
cierto caudal. Estos índices asumen que cada celda puede ser evaluada
independientemente (Jorde, 2006), lo cual debe considerarse solo una aproximación
que no es necesariamente cierta.
El área usable ponderada WUA (Weighted Usable Area) es empleada también por
PHABSIM y se basa en el supuesto de que la sumatoria de las áreas, ponderadas con
la disponibilidad de hábitat en cada celda, es indicador de las condiciones totales del
tramo modelado, para un caudal dado. La siguiente expresión permite calcular éste
índice:
n
WUA = ∑ Ai ⋅ C i = f (Q )
i =1
n
= f (Q )
1
HHS = n ∑ A ⋅C i i
∑A
i =1
i
i =1
30
Este índice corresponde a la calidad promedio del hábitat en el tramo, lo cual permite
expresar los cambios en la calidad global debido a las variaciones de caudal (Jorde,
2006).
31
El modelo, en primer lugar, determina la precipitación neta como la diferencia entre la
precipitación y la evaporación.
Si P ≥ E entonces PN = P – E y EN = 0
Si P < E entonces EN = E – P y PN = 0
Figura 2.6: Esquema general del modelo GR4J descrito por Perrin et al. (2003) (Adaptación)
32
⎛ ⎛ S ⎞2 ⎞ ⎛P ⎞
x1 ⋅ ⎜1 − ⎜⎜ ⎟⎟ ⎟ ⋅ tanh⎜⎜ N ⎟⎟
⎜ ⎝ x1 ⎠ ⎟ ⎝ x1 ⎠
⎝ ⎠
PS =
S ⎛P ⎞
1+ ⋅ tanh⎜⎜ N ⎟⎟
x1 ⎝ x1 ⎠
Cuando la evaporación neta EN es mayor que cero, la tasa de evaporación real también
se calcula en función del nivel S, obteniendo la evaporación ES que se producirá desde
la reserva de la producción (Perrin et al., 2003; Carvajal & Roldán, 2007).
⎛ S⎞ ⎛E ⎞
S ⋅ ⎜⎜ 2 − ⎟⎟ ⋅ tanh⎜⎜ N ⎟⎟
⎝ x1 ⎠ ⎝ x1 ⎠
ES =
⎛ S⎞ ⎛E ⎞
1 + ⎜⎜1 − ⎟⎟ ⋅ tanh⎜⎜ N ⎟⎟
⎝ x1 ⎠ ⎝ x1 ⎠
S = S – ES + PS
⎧ ⎡ 4 −1 / 4 ⎫
⎪ ⎛4 S ⎞ ⎤ ⎪
Perc = S ⋅ ⎨1 − ⎢1 + ⎜⎜ ⋅ ⎟⎟ ⎥ ⎬
⎪⎩ ⎢⎣ ⎝ 9 x1 ⎠ ⎥⎦ ⎪⎭
S = S – Perc
33
Debido a la ley de potencia en la ecuación de la percolación, este término contribuye
principalmente a la simulación de pequeños flujos (Perrin et al., 2003; Carvajal &
Roldán, 2007), como aquel flujo base que persiste luego ocurrido un evento de
precipitación.
El monto total de agua PR que alcanza las funciones de enrutamiento o flujo superficial
es:
7/2
⎛R⎞
F = x 2 ⋅ ⎜⎜ ⎟⎟
⎝ x3 ⎠
34
R = max(0; R + Q9 + F)
⎧ ⎡ 4
⎤
−1 / 4
⎫
⎪ ⎛R ⎞ ⎪
Q R = R ⋅ ⎨1 − ⎢1 + ⎜⎜ ⎟⎟ ⎥ ⎬ donde QR < R para todo instante de tiempo.
⎪⎩ ⎢⎣ ⎝ x 3 ⎠ ⎥⎦ ⎪⎭
R = R - QR
QD = max(0; Q1 + F)
La escorrentía total está dada de esta forma por la suma de los flujos de ambas
componentes:
Q = QR + QD
De acuerdo a Perrin et al. (2003) tanto la división del flujo en 10 – 90% como los
demás coeficientes fueron calibrados en sucesivas pruebas, no demostrándose una
mejora significativa en el desempeño del modelo. Liberar dichos parámetros sería, por
lo tanto, un aumento de la complejidad que no se justificaría de acuerdo a las pruebas
35
efectuadas por estos autores. Las conclusiones de Salvatierra (2008) son las mismas,
tras probar las combinaciones 20 – 80% y 30 – 70% con las funciones objetivo de
Nash-Sutcliffe y Módulo para varias cuencas de Chile central.
El trabajo de Salvatierra (2008) con el modelo GR4J permitió desarrollar una serie de
índices geomorfológicos (IG). Estos permiten determinar los parámetros de calibración
del modelo GR4J en cuencas pluviales no controladas de Chile central. Su estimación
fue realizada mediante la aplicación y calibración con algoritmos genéticos (AG) de
GR4J sobre varias cuencas pluviales de Chile central, razón por la cual se considera a
esta zona como válida para la aplicación de la metodología de Salvatierra (2008).
• Parámetros Geomorfológicos
Los parámetros geomorfológicos (PG) definidos por Salvatierra (2008) para las
cuencas estudiadas en Chile central se resumen en el Cuadro 2.1. Estos son
determinados mediante el uso del software WMS 7.1 (Watershed Modeling System) y
la plataforma SIG de la Comisión Nacional de Riego (CNR).
36
Cuadro 2.1: Parámetros geomorfológicos utilizados por Salvatierra (2008)
Parámetro Unidad Significado
2
A [km ] Área de la cuenca
E [m] Elevación media de la cuenca
i - Pendiente media de la cuenca
Per [m] Perímetro de la cuenca
Lp [m] Largo del cauce principal
Lcg [m] Largo del cauce principal hasta el punto más cercano al centro de gravedad de la cuenca
Dcg [m] Menor distancia entre el cauce principal y el centro de gravedad de la cuenca
SL [m] Suma de los largos de los cauces presentes en la cuenca
AOFD [m] Promedio de distancias de escurrimiento en la superf. de la cuenca, al curso de agua más cercano
Shape - Largo de la cuenca dividido por su ancho
Sin - Óptimo de sinuosidad del cauce principal en la cuenca
2
Aac [km ] Porción de acuífero en superficie localizado dentro de la cuenca
e_ac [m] Espesor del acuífero localizado dentro de la cuenca
• Índices Geomorfológicos
Los IG obtenidos por Salvatierra (2008) se dividen en dos grupos: los determinados
empleando la función objetivo de Nash-Sutcliffe y aquellos derivados mediante la
función objetivo Módulo. Tales funciones objetivo corresponden a la utilizadas en la
calibración de los xi del modelo GR4J, mediante AG.
i. Nash-Sutcliffe
( ) [ ]
n
I N ( x 2 ) = ∑ {Aac j ⋅ SL ⋅ eac j ⋅ i ⋅ AOFD ⋅ sin⋅ shape / A ⋅ Per ⋅ Lp ⋅ Dcg ⋅ E 2 } m −2
j =1
37
( x3 ) N = 91.264 ⋅ I N ( x3 ), 2,1 ⋅ 10 −3 < I N ( x3 ) < 1,4 ⋅ 10 −2 ;
Además, el modelo requiere las condiciones iniciales de los dos estaques o reservas
que lo conforman. Estas condiciones iniciales son determinadas mediante las
siguientes expresiones.
v1 N = 0,75 ⋅ ( x1 ) N + 0,23
ii. Módulo
38
( x 3 )M = 57.488 ⋅ I M ( x3 ) + 309,3, 2,1 ⋅ 10 −3 < I M ( x3 ) < 1,8 ⋅ 10 − 2 ;
Las condiciones iniciales de los estanques son determinadas mediante las siguientes
expresiones.
( )
v1 M = 7 ⋅ 10 −6 ⋅ ( x1 )M − 7 ⋅ 10 −7 ⋅ (E ⋅ i ⋅ Per ⋅ Lcg ) / ( A ⋅ Lp )
v 2 M = 0,246( x 2 )M + 61,52
2.3.1 LEGISLACIÓN
La escasez del agua que se vislumbra para el futuro debido al aumento de su demanda
y que incluso ya se siente en muchos lugares del mundo ha llevado a que el tema sea
debatido en conferencias importantes a nivel internacional, como la Conferencia de
Naciones Unidas sobre el Medio Ambiente y el Desarrollo, desarrollada en Río de
Janeiro en 1992 y el Segundo Foro Mundial del Agua, realizado el 2000 en La Haya
(Bauer, 2005).
39
El consenso general ha sido que las reformas en las leyes relacionadas con el agua
deben enfocarse hacia un manejo integrado de los recursos hídricos (MIRH) a partir de
un enfoque exhaustivo, interdisciplinario y holístico que incluya los aspectos sociales,
políticos, económicos, técnicos y medioambientales (Bauer, 2005). Pahl-Wostl (2008)
va más allá, pues concluye que un manejo integrado y adaptativo de los recursos
hídricos es esencial para garantizar la sustentabilidad de éstos. Para ello señala que es
importante considerar los sistemas de estudios en toda su complejidad, evitando
reducir sus grados de libertad. Además, enfatiza el rol de las plataformas de actores
involucrados y el proceso social de aprendizaje en los múltiples niveles de los sistemas
de gobierno de los recursos (Pahl-Wostl, 2008). No obstante, aunque se reconoce la
necesidad de incorporar la complejidad e incerteza de los sistemas, los proyectos de
manejo adaptativo resultan ser únicamente funcionales y reducen tales elementos, en
parte debido al riesgo de fallar en el estudio de estos sistemas (Allan, 2008). Se trata
de una cultura de manejo adversa al riesgo, lo cual agrega un elemento adicional a la
ya compleja trama de elementos que se deben equilibrar para realizar la tarea del
manejo integrado y adaptativo de los recursos hídricos.
Es más, esto tiene una relevancia mayor en aquellas sociedades con niveles de
desigualdad y pobreza importantes, pues para abordar los nexos entre el
medioambiente, el agua y la pobreza de una comunidad resulta crítico comprender los
vínculos biofísicos, sociales y económicos que existen en dicha composición, porque
son las comunidades pobres las que más dependen, directa o indirectamente, de los
ecosistemas terrestres y acuáticos para su subsistencia y bienestar (Hirji & Panella,
2003).
Sin embargo, no existe un acuerdo completo sobre todos los elementos en juego, pues
aún persiste un fuerte debate acerca del rol de la economía en las reformas perfiladas
en los acuerdos, especialmente sobre las ventajas y desventajas del libre mercado y la
privatización del agua, el cual es tratado como un “bien económico” (Bauer, 2005).
40
neoliberal en el manejo del agua, la cual es tratada como un “bien económico” que
puede ser libremente transado y sin restricciones de un uso específico y/o
“beneficioso”.
Bauer (2005) realiza un extenso análisis del Código de Aguas de 1981 y examina,
entre otras cosas, si su enfoque neoliberal es compatible con el objetivo más general y
de largo plazo del MIRH.
El Código de Aguas fue creado como un marco general que propicia la creación de un
mercado del agua, sin consideraciones de los diversos usos de ésta, sin distinciones
según las diferencias inherentes a la diversidad climática existentes en el país, con una
gran flexibilidad en la entrega de derechos que hoy tiene al sistema dentro de una caja
negra, sin certidumbre adecuada de la cantidad, carácter y ubicación de los derechos
otorgados. La única excepción la representa la definición del uso no consuntivo, el cual
favorece el aprovechamiento del agua como fuente de energía hidroeléctrica pero que,
debido a su pobre definición, ha traído conflictos bien documentados que la
institucionalidad no ha sido capaz de resolver satisfactoria y permanentemente (Bauer,
1998b).
2
Existe un estudio más acabado acerca de las bases de la constitución de 1980 y el rol del poder judicial dentro del
nuevo marco regulatorio neoliberal en Carl Bauer. 1998. Derecho y economía en la Constitución de 1980,
Perspectivas en Política, Economía y Gestión Vol.2, No.1, pg.23-47, (Santiago, Chile).
41
El mismo marco regulatorio y la falta de recursos ha impedido que las atribuciones
dadas por ley a la DGA sean cumplidas a cabalidad, por ejemplo, respecto a la
medición y control de la disponibilidad de recursos hídricos existentes en el país. Esto,
sumado a la incertidumbre acerca de los derechos otorgados, dificulta y encarece
sustancialmente cualquier iniciativa que busque conocer y controlar la distribución de
los recursos hídricos al interior de una cuenca, especialmente si ésta no pertenece al
centro del país.
42
• La protección del medio ambiente y los ecosistemas, incluyendo la mantención
de un caudal ecológico.
• Asistencia pública para la pequeña agricultura en la búsqueda de una mayor
equidad social en temas de derechos y mercados de las aguas.
A lo señalado anteriormente por Bauer (2005) habría que añadir la necesidad de una
asistencia e incentivos para crear verdaderas y fuertes asociaciones de usuarios, con
un importante énfasis en la disposición de recursos humanos y económicos para tales
fines. Sin dichas asociaciones como base, a nivel de cuencas, resulta imposible
generar políticas de manejo de recursos hídricos en forma integrada.
Aun así, ni siquiera la modificación realizada al código el año 2005 ha producido los
efectos deseados. Ésta buscaba, entre otras cosas, otorgar mayor disponibilidad de
recursos hídricos para quienes tuviesen proyectos, favorecer la libre competencia y
eliminar las prácticas monopólicas, la especulación y el acaparamiento del recurso
(Yañez, 2008). Esta reforma requirió un fuerte capital político para ser aprobada, sin
embargo, no logro modificar nada fundamental, como fortalecer la capacidad
regulatoria del Estado o definir qué es un derecho de propiedad de agua (Bauer, 2009).
43
Aunque se han producido importantes remates, como en el río Pescador, aún existen
cerca de 1.400 derechos sin utilizar (Moraga, 2007), pues las empresas han preferido
pagar dichas patentes y conservar sus derechos, por lo cual ya se han comenzado a
gestar nuevas iniciativas que alteren este panorama de eterno acaparamiento.
Por definición, el MIRH concibe a la gestión de los recursos hídricos como un proceso
entre múltiples partes interesadas (multi-stakeholder), en el cual las relaciones entre la
sociedad, la economía y la política son inherentes al desarrollo y gestión de los
recursos hídricos (Iza & Stein, 2009). A través del estudios de casos en diferentes
países, Iza & Stein (2009) logran generar una guía que permite conceptualizar e
implementar mecanismos mejorados y fortalecidos de gestión de los recursos hídricos.
Algunas conclusiones del estudio de Iza & Stein (2009) son las siguientes:
44
formulación de políticas y procesos de negociación en relación al agua debe ser
de responsabilidad pública.
5. El uso de instrumentos económicos es un mecanismo alternativo que apoya y
facilita el cumplimiento de las leyes por parte de los usuarios del agua y otras
partes interesadas.
Las conclusiones de Iza & Stein (2009) respecto a la gestión de recursos hídricos
refuerza el análisis de varios autores consultados sobre el caso chileno (Bauer, 2005;
Latorre, 2006; Moraga, 2007; Yañez, 2008). Por una parte, la política que guió el
Código de Aguas en 1981 afianza los derechos de propiedad sobre el agua, donde
ésta es considerada como un bien económico, con carácter de commodity, susceptible
a la visión y plan estratégico de cualquier empresa o particular que ostente algún
derecho. En segundo lugar, las reformas realizadas al Código de Aguas en los
noventas y en la primera década del siglo XXI han buscado mejorar los incentivos para
un uso provechoso del agua, castigando el acaparamiento. Llegando al final de esta
década, temas como el cambio climático y la desertificación han impulsado a nivel
mundial la visión del recurso hídrico como un bien escaso y primordial para el ser
humano, lo cual a llevado a ciertas autoridades nacionales a plantear la necesidad de
nacionalizar el agua (Bauer, 2009). La contraposición de filosofías en la ley sobre
aguas se hace cada vez más patente y puede seguir creciendo, trabando la posibilidad
de hacer cambios que le permitan a Chile hacer un uso eficiente y sustentable de sus
recursos hídricos.
45
En Chile, haciendo referencia exclusivamente al sistema de estudio de la presente
investigación, existen vacíos importantes de información que ya han sido evaluados al
intentar aplicar ciertas metodologías que han sido utilizadas en otras partes del mundo,
para modelar las necesidades de caudales ecológicos en los ríos.
• Planificar el desarrollo del recurso en las fuentes naturales con el fin de formular
recomendaciones para su aprovechamiento.
• Investigar y medir el recurso.
• Ejercer la policía y vigilancia de las aguas en los cauces naturales de uso
público.
• Supervigilar el funcionamiento de las juntas de vigilancia.
En la práctica, sin embargo, sumado a las trabas descritas en el apartado 2.3.1, la falta
de recursos ha significado que su labor sea parcial e insuficiente para cumplir a
cabalidad lo dispuesto en el Código de Aguas.
Por ejemplo, la red hidrométrica, esencial para una gran cantidad de fines productivos
y sociales, es insuficiente para controlar adecuadamente los recursos hídricos a nivel
de cuencas, como se señala en Universidad de Chile (1998), donde se analizó la
aplicación de la Metodología IFIM, llegando a la conclusión de que la situación en el
país es precaria con respecto a su red de monitoreo, la cual no cuenta con registros
completos en cinco áreas fundamentales: información hidrológica, hidráulica, calidad
físico-química, calidad microbiológica y biótica.
46
En 2005, el Informe País: Estado del Medio Ambiente en Chile (Universidad de Chile),
entregó la disponibilidad de agua por habitante y por región, sin embargo, los datos
empleados son del año 1992 y la fuente es la Dirección General de Aguas. En el
mismo informe, se señala que la demanda de agua del país en el año 1999 es de un
67,8% en el caso del uso hidroeléctrico (uso no consuntivo), proyectando un
crecimiento de esta demanda a un 93% hacia el año 2017. No existen antecedentes
entre 1992 y 2005.
47
subterráneas, parámetros relevantes para la estimación y control correcto de los
recursos hídricos del país.
Desde que Eigenmann publicara en 1928 “The Fresh-water Fishes of Chile” han
existido pequeños cambios en el entendimiento de la diversidad dulceacuícola chilena
(Dyer, 2000). Sin embargo, en los últimos años, principalmente debido al fuerte
aumento de la demanda por nuevas fuentes energéticas y al incremento de las
exigencias ambientales por parte de la comunidad internacional, se han realizado o se
encuentran en curso varios estudios cuantitativos en diversas zonas del país que
intentan determinar de esta forma las preferencias de hábitat de las diversas especies
muestreadas en las zonas estudiadas.
Existe una experiencia con la metodología IFIM, realizada por la Universidad de Chile
(1998) en la que se evaluó la aplicabilidad de dicha metodología en Chile, lo que fue
encargado por la DGA y que significó una serie de campañas de terreno que, sin
embargo, resultaron insuficientes para la magnitud de datos que requiere dicha
metodología. La principal conclusión de dicho estudio fue que la información disponible
para determinar las curvas de preferencia es insuficiente, ya que la información
proveniente de la literatura especializada tiene un carácter eminentemente catastral y
la recopilación en terreno de los datos necesarios debe obligatoriamente cubrir un
período de al menos dos o tres años y cubrir una extensión importante de los ríos
estudiados.
48
legislación existente en su momento haya obligado a concretar dichos estudios,
especialmente en lo que se refiere a la modelación hecha siguiendo la metodología
IFIM en el caso de la central Quilleco.
49
La mayoría de los enfoques biológicos para evaluar el estado o calidad de los ríos y
cauces naturales trabajan sobre tramos individuales de río, sin embargo, no existe una
razón a priori para pensar que lo que ocurre en un tramo es similar en composición a lo
que ocurre en otro (Harby et al., 2004), ni menos en todo el río. A pesar de esto, las
clasificaciones de los ríos a nivel de mesoescala han sido usadas por décadas por la
hidrología y la ecología. Ello ha llevado a buscar métodos de clasificación que faciliten
el escalamiento desde el nivel de modelación al nivel del manejo de recursos, como la
metodología propuesta por Borsányi et al. (2004), que utiliza una aproximación física,
es decir, utiliza parámetros hidromorfológicos que sirvan como herramienta en la
evaluación de mesohábitats en ríos pequeños y medianos. Schwartz & Herricks (2008)
plantean como crítica la necesidad de contar con esquemas de clasificación de hábitats
en ríos que permitan sustentar los programas de restauración de daños, por lo que
plantean en su trabajo un esquema de clasificación de hábitats basado en principios de
la ecohidráulica, lo cual provee un marco conceptual de ecosistema sobre el cual
proponer iniciativas de restauración.
Schwartz & Herricks (2008) señalan que la escala de hábitat más relevante para la
restauración de cauces naturales y la más estudiada ecológica y geomorfológicamente
es la llamada mesoescala, la cual ha sido calificada de variadas formas por diversos
autores, por ejemplo, la escala de poza-rápido (pool-riffle), pozas, caídas, rápidos,
cascadas, vertederos, etc. En general, todas estas clasificaciones están basadas en
condiciones hidrogeomorfológicas, incluso aquellas definidas por los biólogos de
peces. En la actualidad, Schwartz & Herricks (2008) indican el uso de una perspectiva
ecohidráulica que define biotipos físicos diferenciados según biotipos hidráulicos, estos
últimos definidos por parámetros como la velocidad de corte y el número de Froude.
50
Otro camino que cobra cada vez más importancia debido a la capacidad de la
herramienta es el uso del agrupamiento difuso (Fuzzy Clustering). Debido a que los
ríos (referido a su estructura física) son sistemas complejos que necesitan ser
interpretados dentro de un contexto local e histórico, no existe un sistema de
clasificación de tramos de río que pueda sustituir a un observador alerta y bien
entrenado (Montgomery y Buffington, 1998). En este sentido, la experiencia de
Legleiter y Goodchild (2005) al realizar una clasificación mediante el uso de percepción
remota, variables hidráulicas y lógica difusa, demuestra un gran potencial para definir
clasificaciones dinámicas y dependientes de las especies, incorporando la
incertidumbre a través de la medición de otros parámetros y del conocimiento experto.
En Chile, esta clase de estudios no existen, lo cual ha llevado a que no sea posible
distinguir aquellas unidades de estudio más relevantes desde el punto de vista
ecológico, salvo por la experiencia de los biólogos, acostumbrados a encontrar
especimenes en ciertas zonas de los ríos, condiciones de escurrimiento o unidades
geomorfológicas.
La región del Bio Bío exhibe la mayor riqueza de especies, al ser comparada con
cualquier otra cuenca del país. A partir de un estudio realizado por Campos et al.
(1993) se determinó que la hoya hidrográfica del Bio Bío posee 17 especies de peces,
de las cuales 13 son autóctonas o nativas y cuatro son introducidas o exóticas.
Además, 7 especies se clasifican como vulnerables y 7 en peligro de extinción. Una
lista de dichas especies y su categoría de conservación puede ser encontrada en
51
Campos et al. (1993) y en Habit et al. (2006). Además, una lista completa de todas las
especies nativas e introducidas en el país puede ser encontrada en Dyer (2000).
Dado que las especies nativas tienen un nulo interés comercial, su estudio se centra
principalmente en el interés científico que éstas representan. Debido a lo anterior, se
considera que los esfuerzos por proteger los ecosistemas fluviales deben enfocarse en
dichas especies. A continuación se indica la información disponible sobre las especies
nativas que se encuentran en la zona estudiada del río Huequecura, las cuales han
sido capturadas en las campañas de terreno reportadas por el CEA en el año 2003 y
2004.
52
Es descrita como una especie litoral, pues habita dichas zonas, utilizándolas como
áreas de refugio y alimentación, además prefiere los fondos pedregosos y/o con
vegetación acuática sumergida. En las áreas de refugio se reducen los riesgos por
depredación de salmónidos adultos y los costos energéticos de mantención en tramos
con altas velocidades de escurrimiento (CEA, 2003; UdeC, 2004).
Morfológicamente, son peces similares a las percas, pero alcanza tamaños mucho
menores, con un límite aproximado de 10 [cm] de longitud en su estado adulto (UdeC,
2004; CONAMA) y tamaño promedio de 5 a 6 [cm] (CONAMA), siendo clasificados
como adultos a partir de los 5 [cm] (García, 2006). De acuerdo a CONAMA, los
ejemplares adultos de la especie presentan coloración naranja amarillenta brillante. El
dorso es oscuro y el vientre claro. La coloración del dorso es gris oscura, con tonos
verde amarillentos, ligeramente más claro hacia los flancos. Vientre claro, cremoso
amarillento. Por lo general presenta aletas claras con bandas oscuras.
En cuanto a sus hábitos alimenticios, han sido detectados patrones diferentes según su
ubicación. En los tributarios del río Bio Bío su dieta está compuesta principalmente por
Efemerópteros (40.9%), Coleópteros (24.4%) y Dípteros (21.1%), mientras que en el
río Bio Bío los Coleópteros desaparecen de su dieta, concentrándose principalmente
en los Efemerópteros (55.9%) y Dípteros (35.9%) (CEA, 2003).
53
2.3.3.2 DIPLOMYSTES NAHUELBUTAENSIS
Se encuentra presente en las hoyas hidrográficas de los ríos Itata, Bio Bío e Imperial
(CONAMA). Ha sido reportada desde la cuenca del río Maule a la del río Imperial
(Dyer, 2000). Se ha recolectado en los ríos Cautín, Lautaro (ríos Loncomilla y San
Javier), Estero Nongén, Altos del Bio Bío, Malleco y en el río Laja (CONAMA).
Actualmente está clasificada como “en peligro de extinción” (CONAMA; Dyer, 2000;
CEA, 2003; UdeC, 2004).
Antecedentes sobre su biología indican que existe una separación entre los ejemplares
juveniles y adultos en cuanto al hábitat, debido que los adultos viven en aguas más
torrentosas y profundas. En cambio, los juveniles se asocian preferentemente a zonas
menos torrenciales ribereñas, además se han descrito en ambientes de ritrón,
principalmente en zonas de pozones profundos con bentos, sustratos de gravilla a
bolones con bajas temperaturas y alta oxigenación (CONAMA).
54
Esta especie habita las zonas que poseen como sustrato principalmente piedras y
bolones. Los juveniles habitan los sectores ribereños de baja velocidad de corriente y
profundidad, mientras que los adultos habitan sectores más profundos y con mayor
velocidad de corriente (UdeC, 2004), siendo vistos incluso a más de 20 metros de
profundidad en la cuenca del Bio Bío (Manuel Contreras, Comunicación directa, 2007).
Esta especie tiene una tasa baja de fecundidad, por lo que podría alcanzar sólo dos
eventos reproductivos a lo largo de su vida (CONAMA).
55
T. areolatus, si bien es considerada como una especie vulnerable (desde la III a la X
región según CONAMA), en la mayoría de los ríos chilenos que ha sido reportada
presenta altas abundancias en sectores de bolones y de aguas someras, los que
constituyen su hábitat preferencial. Además, T. areolatus presenta cierta plasticidad
trófica, ampliando su nicho trófico en determinados sistemas, tanto lóticos como ritrales
(Habit et al., 2003). Arratia (1983), por su parte, clasifica a T. areolatus como un típico
habitante de ritrón.
3
Esta zona es definida como aquella ubicada en el fondo del río, donde, debido a las macrorugosidades, se produce el
fenómeno de despegue de la capa límite y aparece la resistencia de forma. Debido estas macrorugosidades se originan
zonas de aguas muertas entre las rocas de fondo que la especie aprovecharía como refugio.
56
pedregoso y con una velocidad de agua superficial que varía entre los 0.1 a 0.3 m/s en
diferentes ríos y esteros, en diferentes épocas del año.
Las zonas ribereñas de escasa profundidad y de velocidad de corriente casi nula, son
habitadas por formas juveniles entre octubre y abril, siendo el período más frecuente
entre diciembre y abril (Arratia, 1983).
57
Bío y Andalién; está clasificada como “en peligro de extinción” (CONAMA; Dyer, 2000;
Habit et al., 2006).
Este ejemplar presenta una coloración con manchas negras de variados tamaños y
gran cantidad. Se diferencia de T. areolatus por la posición en que se ubica la aleta
dorsal (CONAMA).
Sobre sus preferencias de hábitat, tanto en estado juvenil como adulto, Arratia (1983)
señala que éstas son relativamente similares a las de T. areolatus, es decir, que sería
una especie preferentemente de ritrón, donde las pendientes son elevadas, con
dominio de sustratos duros y falta de sustratos fangosos y agua de elevada
transparencia. La única diferencia entre T. areolatus y T. chiltoni se produciría en su
estado adulto, con respecto a la variable velocidad, ya que T. chiltoni sería más
sensible a velocidades superficiales altas, prefiriendo un rango de éstas que varía entre
los 0.1 y 0.92 m/s (Arratia, 1983).
Al igual que T. areolatus, las formas juveniles de esta especie ocupan las zonas
ribereñas de escasa profundidad y de velocidad de corriente casi nula, entre octubre y
abril, siendo el período más frecuente entre diciembre y abril (Arratia, 1983).
58
depredador carnívoro, como ocurre en el río Maipo, lo que demuestra una flexibilidad
trófica. El período de reproducción de estos peces está indicado entre los meses de
primavera-verano (CONAMA).
59
noviembre hasta principios de diciembre (CONAMA). No obstante, se registra que a lo
largo de su vida, B. maldonadoi cambia su preferencia de hábitat desde un ambiente
arenoso a otro constituido fundamentalmente de grava y piedras pequeñas (Arratia,
1983; Habit y Victoriano, 2004).
60
Su talla puede ser mayor a 30 cm. de longitud total. Su coloración se distingue por
presentar una banda oscura longitudinal, plateada en su borde inferior, flancos verde
amarillentos, más blancos; el dorso está recorrido por una banda oscura angosta, de
coloración verde azulada (CONAMA). La Figura 2.12 permite apreciar en parte esta
descripción.
4
Nuevos peces que se suman a la población existente cada año.
61
2.4 ESTADO DEL ARTE EN LA PROBLEMÁTICA DE LA GESTIÓN DE
RECURSOS HÍDRICOS CON ENFOQUE MULTIDISCIPLINARIO
• Conceptuales
• Basados en la experiencia, modelos expertos
• Matemáticos o mecanicistas
• Estadísticos o estocásticos
• Basados en lógica difusa
• Modelos a escala o modelos físicos
• Basados en inteligencia artificial
62
El objetivo de la siguiente revisión es señalar el rumbo actual, en discusión a nivel
internacional, respecto al manejo de los recursos hídricos, entendido esto como el
estudio del recurso agua y de todos los elementos que interactúan con él, así como el
análisis respecto a la forma en que estos elementos se relacionan, las metodologías
empleadas y los actores de la sociedad que deben involucrarse en dicho manejo.
63
entender la estructura y dinámica del sistema social – ecológico en el cual se desarrolla
la distribución de costos y beneficios debido al uso de los recursos.
i. Aspectos concretos del trabajo científico en los estudios, entendidos estos como
las herramientas usadas para la construcción conceptual de un fenómeno
observado: la generalización que permite construir el patrón, la explicación
causal que permite cuantificar las relaciones y las pruebas conducentes a validar
lo anterior. El uso sistemático de una sola herramienta en cada disciplina limita el
entendimiento entre éstas.
64
son estudiados a diferentes niveles, dependiendo de las disciplinas, pero en
conjunto permiten extraer los elementos claves que conducen al entendimiento
de la red de causalidades más elementales y que permite a los administradores
manipular esta red.
65
En forma resumida, la Figura 2.14 presenta el marco conceptual definido en el trabajo
de Redman et al. (2004), donde se puede ahondar en las características particulares
de dicho marco.
Figura 2.14: Marco conceptual definido por Redman et al. (2004) (Adaptación)
66
diferentes disciplinas y que poseen diferentes perspectivas, de modo que lo anterior
sea aplicable.
Esta herramienta consta de cuatro módulos (King et al., 2003; King et al., 2006):
i. Un módulo biofísico diseñado para predecir los efectos del cambio en el régimen
de flujos sobre el ecosistema fluvial.
ii. Un módulo sociológico diseñado para identificar a los usuarios vitales en riesgo
debido a la intervención del río, cuantificando las relaciones de las personas con
el recurso.
iii. Un módulo constructor de escenarios que permite reunir las diversas predicciones
de cambios en el río y el impacto sobre las personas.
iv. Un módulo económico que permite resumir los costos y beneficios de cada
escenario, lo que permite considerar mitigaciones y compensaciones.
67
de cualquier proyecto. Un caso representativo de aplicación es el presentado por
Arthington et al. (2003), donde se utiliza a los peces como especie característica de
estudio. Así, se obtiene un protocolo de diez pasos que involucra la componente de
estudio de los peces y se compara a DRIFT con otras metodologías de evaluación
basadas en los paneles de científicos y expertos.
68
d) Un Nuevo Procedimiento de Evaluación para Cuantificar la Compatibilidad
Medio Ambiental de la Producción Hidroeléctrica
Para lograr los objetivos del proyecto se utilizaron tres enfoques metodológicos:
69
ii. La herramienta de evaluación será un modelo de calidad de hábitats (en el
estudio se emplea AQUASIM para modelar la temperatura y la calidad del agua
en conjunto con CASiMiR).
iii. La literatura será aquella que permita cumplir el punto 1 utilizando las
herramientas definidas en el punto 2.
Figura 2.16: Esquema general de matriz de manejo medioambiental (Adaptado de Bratrich et al.,
2004)
70
2.4.2 PROYECTOS DE ESTANDARIZACIÓN E INTEGRACIÓN A MACROESCALA
a) AQEM
5
Como parte del proyecto se desarrollo el programa de evaluación ASTERICS, el cual puede ser descargado, junto
con la lista de taxas de macroinvertebrados, manual y complementos desde su sitio http://www.aqem.de
71
llamados métricas (metrics), con el objetivo de identificar aquellos atributos más
efectivos en la medición de la degradación en el río, lo cual fue estimado mediante
correlaciones de dichas métricas con la degradación de la calidad del río, medida a
través de parámetros químicos o hidromorfológicos. El proceso permitió determinar
más de 18 métricas esenciales para los tipos de ríos definidos, lo cual fue combinado
en diferentes índices multimétricos en cada país (Hering et al., 2004).
b) PAEQANN
72
Formalmente el objetivo de aplicación de este proyecto es proponer un conjunto de
políticas y herramientas para el manejo del agua, facilitando de esta forma la
evaluación de la calidad y de las perturbaciones en los ecosistemas fluviales,
entregando información acerca de la calidad del flujo de agua tanto como de la
estructura de la comunidad que conforma al ecosistema.
a) HSSDex
6
En el sitio web http://aquaeco.ups-tlse.fr/index2.htm es posible descargar algunos artículos relacionados con el
desarrollo de esta herramienta, la gran mayoría relacionados con el uso e implementación de redes neuronales en la
modelación de ecosistemas.
73
de entrevistas a científicos de la industria, de estudios basados en el trabajo hecho en
terreno y de la literatura disponible sobre el tema.
Se utilizó un enfoque conceptual del tipo bottom up7 para desarrollar el modelo, de tal
forma que mediante reglas heurísticas se definieron relaciones entre los factores
ecológicos, biológicos y motivacionales y sus efectos en el comportamiento de los
peces. De esos factores, son considerados como claves la comida, la predación y el
estado reproductivo. La interrelación dinámica de estos atributos, combinada con los
cambios temporales en las prioridades de vida de los peces, dan como resultado una
serie de trade-offs8, en la forma de respuestas de comportamiento que se manifiestan
en la estructura, dinámica y distribución del banco de peces (Mackinson, 1999).
7
En el diseño Bottom-up las partes individuales se diseñan con detalle y luego se enlazan para formar componentes
más grandes, que a su vez se enlazan hasta que se forma el sistema completo. Las estrategias basadas en el flujo de
información "bottom-up" se antojan potencialmente necesarias y suficientes porque se basan en el conocimiento de
todas las variables que pueden afectar los elementos del sistema (Fuente: http://es.wikipedia.org/wiki/Bottom-up).
8
Un trade-off se refiere a perder una cualidad o aspecto de algo a cambio de ganar otra cualidad o aspecto. Esto
implica que la decisión debe ser realizada con una comprensión completa de lo que se gana y lo que se pierde en una
decisión particular (Fuente: http://en.wikipedia.org/wiki/Trade-off)
74
especie). La validación se realizó con los datos físico-químicos y de
macroinvertebrados muestreados.
Se creó un modelo para cada especie, aunque son aplicables a dos especies de
Asellus y varias especies de Gammarus, debido a las pocas diferencias en sus nichos
y en las tolerancias a la polución orgánica, respectivamente. Esto le da mayor
aplicabilidad a esta herramienta, que pretende ayudar en el manejo de la calidad de
ríos y en la evaluación de impactos ambientales mediante el uso de estos
macroinvertebrados como bio-indicadores de la calidad del agua (Adriaenssens et al.,
2006).
El estudio fue realizado en el lago Shasta, el embalse más grande de California, con el
fin de investigar los efectos que tendría sobre la red de alimentos (food web) del
embalse el uso de mecanismo de control de temperatura9 (TCD: temperature control
device) (Saito et al., 2000).
9
Un mecanismo que permite a los operadores de embalse tomar el agua desde diferentes profundidades, incluyendo el
agua más cálida cercana la superficie, regulando no solo la temperatura a la salida, sino que modificando lo que
sucede en el mismo embalse (Fuente: http://www.gcdamp.gov/tcd/index.html).
75
El modelo de transferencia food web-energy fue generado mediante un análisis de
isótopos estables10 el cual permitió descubrir las relaciones tróficas entre los
predadores y múltiples presas, estimando la posición trófica de cada consumidor,
definiendo además la importancia relativa de varias presas en la dieta de cada
consumidor.
Con este modelamiento fue posible estimar los efectos bottom up relativos que nuevas
reglas de operación pueden provocar en varias especies de peces (Saito et al., 2000).
10
Existen análisis isotópicos en estudios metabólicos, en biodegradación, en ecología microbiana, en arqueología, en
estudios paleoambientales. Una búsqueda simple en Internet bastará para conocer más sobre este tema (Fuente:
http://www.eez.csic.es/files/Los%20Isotopos%20Estables.pdf).
76
• Datos de calidad de agua (N, P, coniformes totales, turbiedad, etc.) para el
QUAL2E
• Profundidad del agua, velocidad de la corriente, tamaño del sustrato y datos
sobre uso de hábitat de los peces (mediante buceo con snorkel) para el modelo
HABITAT.
Se midió el impacto para diferentes caudales en cada uno de los modelos, lo cual
permitió a los científicos integrar los resultados en la recomendación final de un caudal
ecológico.
77
de estudio aguas abajo de las represas. Para esto se lleva a cabo el
modelamiento bioenergético y del hábitat de los peces.
• El estudio del comportamiento de los peces y la selección de los tipos de
refugio debido a la rápida variación del nivel del agua, investigando en
laboratorio y subsecuentemente en situaciones naturales empleando telemetría
de alta precisión.
• El monitoreo de la vegetación acuática y del bentos, unido a la modelación
hidráulica de las condiciones de crecida en los ríos. Además, se realiza un
análisis de series de tiempo para mostrar el impacto de la variación de las
condiciones hidráulicas en la composición y distribución de la vegetación
acuática.
• El estudio del consumo de energía, el crecimiento y el estrés en los peces
juveniles para cuantificar los costos energéticos debido a los rápidos y
frecuentes cambios en los niveles del agua.
2.5 COMENTARIOS
Sin embargo, estas definiciones son muy generales y vagas, pues no existe claridad
absoluta acerca de lo que se denomina desarrollo, integrador, flexible, sustentable.
Cada uno de estos términos posee definiciones diferentes y las formas en que las
diferentes sociedades se aproximan a ellos también difieren.
78
Quizás existe un amplio acuerdo acerca de la importancia de la ciencia, ya sea como
medio para alcanzar la verdad o como herramienta de conocimiento, de comprensión,
de gestión.
79
Sin duda, lo más complejo del actual escenario pasa por el modelo de desarrollo que
se ha dado hasta hoy y si éste es compatible con un desarrollo sustentable. En este
punto, los sistemas políticos y económicos serán los más tensionados en el futuro.
En Chile existe un modelo de desarrollo neoliberal que exige cierta rigidez del sistema
político-legal, tal que asegure jurídicamente los derechos de propiedad para dar
confianza a los inversionistas. Esta falta de flexibilidad no permitirá una adaptación
suficiente ante la escasez de agua. Un cambio en la disponibilidad de este recurso
afectará nuestra economía en todos sus campos de desarrollo, tanto en la minería
como la exportación de recursos naturales primarios.
Culturalmente hay retrasos en Chile y con un gran desafío por delante. Muchos autores
han señalado la importancia de ciertos elementos señalados en la presente revisión,
como la libertad de acceso a la información, menor centralización, etc. Sin embargo,
son pocos los que ven la integración de todos estos elementos, siendo propuestos
aisladamente por la mayoría. Por ejemplo, aunque existe acuerdo acerca de la
importancia de la adaptabilidad del sistema, algunos persisten en confiar en el mercado
de derechos de agua para aquello. Otros, incluyendo al presente autor, consideran que
un modelo de esas características no permite una adaptación suficiente. No sólo por la
falta de mercados en la actualidad, bajo el supuesto cuestionable de que agua puede
ser considerado un bien transable, sino por el conjunto de medidas centralizadoras y
totalizadoras necesarias para su sustento y que hoy impiden un modelo de desarrollo
regional y acorde con las necesidades de cada zona.
Otra área básica y paralela a la anterior, tiene relación con el énfasis que se da a la
ciencia en Chile. Aunque las inversiones en investigación continúan aumentando, no
existe un horizonte de desarrollo nacional. Por ejemplo, la falta de inversiones en la
80
medición de los recursos hídricos chilenos produce una base científica permeada por la
incertidumbre. Los resultados son demasiado vagos para permitir que las autoridades
los tomen en cuenta en el desarrollo de las políticas nacionales. No es posible reformar
el actual sistema de derechos de agua si no se puede estudiar adecuadamente como
se verán afectados por un cambio en la cantidad y distribución del mismo. La falta de
información le da poca credibilidad a los resultados y esto impide que se invierta en
mejorarlos. El círculo es claro, pero muy difícil de romper.
Algo similar ocurre en proyectos privados, por ejemplo, aquellos relacionados con
estudios de impacto ambiental. Se estudia la naturaleza que se verá afectada por un
proyecto de inversión, generando una línea base de los recursos naturales, sociales y
culturales que se podrán ver afectados, pero los resultados son expuestos en forma
parcial y gran parte de la información base utilizada para generar los resultados es
reservada. Hoy en día, los miembros de la sociedad (civil, científica, autoridades, etc.)
que deseen acceder a esta información para revisar, comprobar o utilizarla, no pueden,
a menos que formen parte de quienes realizaron o encargaron los estudios.
81
Esta clase de asimetría de la información es interpelada desde diferentes frentes. La
Organización de Cooperación y Desarrollo Económicos (OECD), en el informe de
evaluación ambiental realizada a Chile, enfatizó este punto como uno de los más
importantes a ser corregidos. Algunas de las recomendaciones realizadas por este
organismo son las siguientes (OECD, 2005):
Como señala Pahl-Wostl (2008), no existe una panacea para desarrollar procesos de
aprendizaje y el manejo adaptativo de los recursos. Nuestro entendimiento de las
dinámicas complejas de los sistemas acuáticos, incluyendo componentes humanas,
ambientales y técnicas es aún muy pobre, no obstante, estamos todos invitados a
contribuir desde nuestra propia experiencia para desarrollar aún más este enfoque
integral y adaptativo de desarrollo (Pahl-Wostl, 2008).
82
3. CARACTERIZACIÓN ZONA DE ESTUDIO
El río Bio Bío es un sistema fluvial fuertemente intervenido, con aproximadamente 380
km de largo, y una cuenca hidrográfica de 24260 km², situada entre los 36°45'S y
39°S, la tercera más grande del país, después de las de los ríos Loa y Baker (Parra et
al., 2004). El río Huequecura, ubicado en la comuna de Santa Bárbara, VIII región, es
afluente al río Bio Bío, formando parte de una de las cuencas más importantes y
estudiadas de Chile.
Por su parte, la cuenca del río Huequecura, cuya salida fue definida en la zona de
estudio, en las coordenadas 71.786º de longitud Oeste y 37.706º de latitud Sur, posee
un área total de 338.9 km2. La Figura 3.1 señala la ubicación general y tamaño relativo
de la cuenca (Fuente: Elaboración propia).
Figura 3.1: Ubicación General Zona de Estudio. De izquierda a derecha Chile, río Bio Bío, cuenca
del río Huequecura delimitada en zona de estudio.
83
La comuna de Santa Bárbara ocupa una superficie de 1255 km2 con una población,
según el Instituto Nacional de Estadística (INE), de 14271 habitantes, dando una
densidad poblacional de 11.4 habitantes por km2, revelando así el carácter rural de la
zona. Las campañas de terreno indican que no existen poblados o grupos de viviendas
importantes. Los asentamientos están dispersos, sin una agrupación clara que permita
definir poblados o localidades.
84
Figura 3.2: Diversidad de tipos geomorfológicos en el río Huequecura (Fuente: registro y
elaboración propia)
3.2 CLIMA
3.2.1 PRECIPITACIONES
85
El régimen de precipitaciones en la zona de estudio fue determinado utilizando la
estación pluviométrica Cerro El Padre, ubicada en las coordenadas 71.867º de longitud
Oeste y 37.783º de latitud Sur, a una altitud de 400 m.s.n.m., a poco kilómetros de la
salida de la cuenca del río Huequecura. La Figura 3.3 muestra el régimen de
precipitaciones que domina a la zona.
850
800
750
700
Precipitación Total Mensual [mm]
650
600
550
500
450
400
350
300
250
200
150
100
50
0
Gumbel Pearson Pearson Pearson Gumbel Normal Gumbel Pearson Gumbel Pearson Pearson Gumbel
Enero Febrero Marzo Abril Mayo Junio Julio Agosto Septiembre Octubre Noviembre Diciembre
Pnoexc 5% 0.0 0.0 6.8 30.05 29.70 177.65 99.92 91.17 35.79 26.37 7.55 0.0
Pnoexc 25% 12.33 8.00 33.77 90.45 175.42 317.32 205.17 171.86 96.52 82.38 42.78 26.35
Pnoexc 50% 56.52 30.84 59.06 149.19 306.50 414.30 299.85 234.51 151.15 136.26 81.40 70.12
Pnoexc 75% 112.49 70.69 90.87 225.04 472.81 511.27 419.97 302.88 220.46 205.30 134.14 128.40
Pnoexc 85% 153.81 104.26 108.21 267.10 580.80 563.38 497.97 342.06 265.46 243.39 169.08 166.08
Pnoexc 90% 183.74 129.23 121.92 300.67 662.76 598.67 557.17 369.59 299.62 273.72 195.58 191.31
Pnoexc 95% 264.20 171.49 143.61 354.24 798.89 650.95 655.49 411.90 356.35 321.99 239.35 237.96
Del área total de la cuenca, el 84% se encuentra bajo los 1500 m.s.n.m. y 98% bajo los
1800 m.s.n.m. lo cual indica indirectamente la escasa influencia de las precipitaciones
sólidas en el régimen de caudales del río Huequecura, pudiendo clasificar a éste como
de régimen pluvial. La cota mínima de la cuenca es 308 m.s.n.m. y la máxima es 2187
m.s.n.m.
86
3.2.2 TEMPERATURAS MEDIAS ANUALES
El área del piedmont y del valle del Bio Bío se caracteriza por una disminución suave
de las temperaturas medias anuales, hasta los 8 ºC. La Alta Cordillera se extiende
sobre los 2000 m.s.n.m. y se encuentra enmarcada por temperaturas inferiores a los 4
ºC (EIC, 1996).
Debido la inexistencia de una red de estaciones más densa se dará un valor referencial
de la tasa de evaporación media anual, correspondiente a lo registrado por EIC (1996)
en la zona de Pangue y Troyo, donde se observa un comportamiento homogéneo,
fluctuando entre 1133 y 1221 mm anuales respectivamente.
87
3.2.4 RED DE ESTACIONES HIDROLÓGICAS
88
No obstante, los estudios llevados a cabo por el Centro de Ecología Aplicada (CEA)
con motivo de la construcción de la central hidroeléctrica Pangue revelan que, en la
mayoría de los parámetros relevantes, el río Huequecura mostró las menores
concentraciones de entre todos los tributarios al río Bio Bío, lo que refleja un río con
escasos o poco importantes impactos antropogénicos.
Esta razón, sumada a los muestreos de peces realizados por el CEA y al poco grado
de intervención de la cuenca, fueron los principales criterios que llevaron a seleccionar
al río Huequecura para realizar el presente estudio.
Los parámetros medidos por el CEA (ver Cuadro 3.1 a Cuadro 3.4) (CEA, 2004b;
2004c; 2004d; 2004e), incluyen pH, temperatura (Tº), oxígeno disuelto (OD),
conductividad eléctrica (Cond. Elec.), alcalinidad (Alc-total), sólidos suspendidos totales
(S.T.S.) y disueltos (S.T.D.), amonio (NH4), nitritos (NO2), nitratos (NO3), nitrógeno
orgánico total (N. Org. Total), ortofosfato (P-PO4), fósforo total (P Total), sodio (Na+),
potasio (K+), magnesio (Mg+2), calcio (Ca+2), cloruro (Cl-), bicarbonato (HCO3-),
carbonato (CO3-2), sulfato (SO4-2), sílice (SiO3-2), clorofila a, demanda bioquímica de
oxígeno (D.B.O.), coliformes fecales (C.F.) y coliformes totales (C.T.). Además, se
tienen parámetros no detectados (ND) y detectados (D) pero con una concentración
insuficiente para ser cuantificados.
El Cuadro 3.5 contiene los límites máximos establecidos por norma chilena para la
descarga a cuerpos de agua fluviales, indicando solo aquellos parámetros que
coinciden con los medidos por el CEA, salvo la DBO5, que corresponde a la DBO
degradable en 5 días y que, en el caso de descargas domiciliarias, es equivalente al
70% de la DBO total carbonácea (Sancha, 2004).
89
Cuadro 3.2: Nutrientes en río Huequecura, año 2004. ND: No Detectado. D: Detectado.
Mes Año NH4 NO2 NO3 N Org.-Total P-PO4 P-Total
(µg / l) (µg / l) (µg / l) (µg / l) (µg / l) (µg / l)
Febrero 2004 ND ND ND 148.80 ND ND
Mayo 2004 ND D ND D ND ND
Agosto 2004 D 0.5 ND ND ND ND
Noviembre 2004 ND ND ND 90.00 ND ND
Cuadro 3.5: Límites máximos permitidos para la descarga de residuos líquidos a cuerpos de agua
fluviales según el Decreto Supremo Nº 90.
90
En el río Huequecura la DBO total es menor a aquella de más rápida degradación
(DBO5), indicando un muy bajo contenido de materia orgánica biodegradable. En este
sentido, de la carga de nutrientes detectada (Cuadro 3.2) solo destaca la presencia de
nitrógeno orgánico total en la época de estiaje, al mismo tiempo que no se detectó
ninguna otra forma de compuestos orgánicos nitrogenados. La elevada concentración
de oxígeno detectada en las mismas campañas indica que se la presencia únicamente
de nitrógeno orgánico se debe a una contaminación reciente y/o a escasa distancia del
punto de muestreo. Su origen pude ser explicado por la mayor presencia humana,
debido a los usos recreativos del río que se producen principalmente en verano. Sin
embargo, la ausencia de compuestos fosforados, la nula presencia de coliformes y la
baja concentración de DBO (Sancha, 2004), apuntan a que el origen del nitrógeno es
natural, producto del metabolismo de la vida biológica presente en el río.
Este resumen permite verificar lo señalado por CEA (2000), en relación a que no se
observa un decaimiento de la calidad entre el año 2000 y 2004. Algunos valores fuera
del rango habitual no son justificados por los demás parámetros, no pudiéndose
concluir ninguna tendencia definitiva en la zona.
91
La Figura 3.5 muestra mediante colores los diversos usos que existen en los registros
de la CNR. El Cuadro 3.6 resume cuantitativamente dichas áreas.
92
3.5 AGUAS SUBTERRÁNEAS
La Figura 3.6 indica el acuífero que es reportado por la CNR en su plataforma SIG. El
área del acuífero es 86.6 km2 y la profundidad media es de 56 metros.
A modo de descripción, se tiene que el volumen estimado del acuífero es de 4.5 km3 y
el coeficiente de almacenamiento es de 0.9.
3.6 COMENTARIOS
Las condiciones de la zona de estudio resultan ser idóneas para aplicar un modelo de
calidad de hábitat desacoplado de otros estudios, debido a la poca intervención
93
humana del lugar. Además, se simplifican los balances y es posible plantear cualquier
escenario futuro como posible, reflejando directamente dichas posibilidades en los
resultados de la modelación. Por una parte, la calidad del agua permite cumplir con el
principal supuesto de CASiMiR, el cual se refiere a que son las condiciones físicas del
río (velocidad del agua, profundidad y sustrato del lecho) las que determinan la
habitabilidad de una determinada especie. Por otro lado, el bajo uso de los recursos
hídricos permite que los peces estudiados en la zona deban sus poblaciones actuales
casi exclusivamente a las condiciones climáticas propias del lugar.
Con respecto a la selección del tramo estudiado (ver Punto 3.1), no existe una
justificación con respecto a su validez como unidad representativa del río desde el
94
punto de vista ecológico. No obstante, la experiencia empírica señala que la zona
permite condiciones apropiadas para varias especies nativas, lo cual sustenta la
selección, sin que esto sea lo deseable ni lo más adecuado para llevar adelante esta
clase de estudios en el futuro.
95
4. MODELACIÓN HIDROLÓGICA
96
Cuadro 4.1: Parámetros geomorfológicos río Huequecura.
Parámetro Unidad Valor Significado
2
A [km ] 338.85 Área de la cuenca
E [m] 1026.2 Elevación media de la cuenca
i - 0.395 Pendiente media de la cuenca
Per [m] 122040 Perímetro de la cuenca
Lp [m] 30304.7 Largo del cauce principal
Lcg [m] 16605.3 Largo del cauce principal hasta el punto más cercano al centro de gravedad de la cuenca
Dcg [m] 4196.8 Menor distancia entre el cauce principal y el centro de gravedad de la cuenca
SL [m] 75836.2 Suma de los largos de los cauces presentes en la cuenca
AOFD [m] 1855 Promedio de distancias de escurrimiento en la superf. de la cuenca, al curso de agua más cercano
Shape - 1.88 Largo de la cuenca dividido por su ancho
Sin - 1.2 Óptimo de sinuosidad del cauce principal en la cuenca
2
Aac [km ] 86.56 Porción de acuífero en superficie localizado dentro de la cuenca
e_ac [m] 52 Espesor del acuífero localizado dentro de la cuenca
e_ac En ArcGIS utilizar herramientas para delimitar porción contenida dentro de cuenca
La definición empleada corresponde a cuencas pluviales (en Chile central) con cursos
de agua permanentes a lo largo del año, es decir, no se consideran quebradas o
cursos eventuales surgidos durante un evento de precipitación.
No se ha definido un criterio que señale hasta qué período de retorno, por ejemplo,
debe considerarse dicho flujo ni se ha verificado en terreno si los cursos definidos en la
Figura 4.1 son efectivamente permanentes. La definición está sujeta principalmente al
97
software WMS 7.1 y al valor umbral (threshold value) que se debe definir. Un criterio
más amplio es el propuesto por Salvatierra (2008b) referido al flujo permanente, pero
es necesario definir algún criterio, tal como cierto período de retorno o probabilidad de
excedencia. Debido a que se trata de una cuenca no controlada por una estación
fluviométrica, esta definición debe estar ligada a la estadística pluviométrica existente.
Figura 4.1: Cuenca definida en WMS 7.1. Valor Umbral = 5. Cauces definidos se destacan en negro.
Cauces en azul corresponden al total definido dada la resolución de la topografía digital (escala 30
m) (Fuente: Elaboración propia con topografía USGS).
Otra variable influenciada por el threshold value es el AOFD (ver Cuadro 4.1), debido a
que el agua recorrerá más o menos distancia, desde donde cae la gota de lluvia hasta
el cauce más cercano, dependiendo de los cauces efectivamente definidos. También
se verá afectado el valor SL (ver Cuadro 4.1), debido a las mismas razones. No
obstante, es claro que la definición del número de cauces relevantes afecta a todas las
variables por igual, siendo parte integral de la definición de los IG determinados por
Salvatierra (2008).
98
Los parámetros Aac y e_ac corresponden a la información contenida en la plataforma
SIG de la CNR, por lo que su incertidumbre está contenida en la calidad de la
información utilizada para obtener los mapas empleados. Es indudable que un espesor
único del acuífero no es real, pero debido a que GR4J es un modelo concentrado, no
es posible aumentar la resolución a menos que se divida la modelación en
subcuencas, con los balances hídricos correspondientes. La información, sin embargo,
no es suficiente para realizar tal aumento de la resolución.
Los IG obtenidos se presentan en el Cuadro 4.3, junto con la evaluación del rango o
límite de aplicación definido por Salvatierra (2008). Dicho límite fue definido a partir de
los resultados obtenidos por Salvatierra (2008), sin embargo, no significa que sean
absolutos (César Salvatierra, Comunicación directa, 2008). En el Cuadro 4.4 se
presentan los parámetros de GR4J obtenido a partir de los IG.
99
El índice IN(x2) (Nash-Sutcliffe) se excede por el límite superior definido, el cual está
asociado al parámetro de mayor sensibilidad del modelo (Salvatierra, 2008). No
obstante, la sensibilidad definida se refiere a cómo varía la respuesta del modelo GR4J
ante la dispersión de la correlación del que se deriva IN(x2) y no a como varía x2 al
variar dicho índice.
La Figura 4.2 fue extraída del trabajo de Salvatierra (2008) y da cuenta de cuál es la
incertidumbre asociada a cada parámetro. Como se señaló, la incertidumbre se refiere
a la dispersión de la correlación.
Esto permite aceptar que IN(x2) esté fuera del rango definido, pues la variación de x2 (y
el caudal resultante) es muy pequeña, en contraste con la incertidumbre propia de la
correlación que permite su estimación.
Figura 4.2: Esquema del procedimiento de obtención de parámetros estadísticos para análisis de
incertidumbre (Fuente: Salvatierra, 2008)
100
El resultado asociado a la función objetivo Módulo entrega el índice IM(x1) fuera del
rango. Debido a que la variación de x1 es más importante que en el caso de Nash-
Sutcliffe, se realizó un análisis de sensibilidad idéntico al anterior. Este demostró que,
ajustado IM(x1) para cumplir la restricción, el error absoluto medio de la modelación es
de 2.2 [m3/s]. Este error, sin embargo, se explica casi por completo en las crecidas.
Para caudales posteriores a los eventos de precipitación, el error tiende a cero.
La validez de esto puede ser cuestionable, pero las variaciones en los resultados caen
dentro de la incertidumbre propia de esta metodología, como se puede ver en el
análisis de sensibilidad e incertidumbre realizado por Salvatierra (2008).
Figura 4.3: Relación entre IG y x2 para Nash-Sutcliffe. La proyección de los datos extremos en el eje
X indica el rango definido de cada IG (Fuente: Salvatierra, 2008)
101
Figura 4.4: Relación entre IG y x1 para Módulo. La proyección de los datos extremos en el eje X
indica el rango definido de cada IG (Fuente: Salvatierra, 2008)
La estadística de precipitación presentó algunos días con vacíos, los cuales fueron
rellenados utilizando la estación pluviométrica Cerro El Padre (BNA 08317003-4), que
cuenta con datos diarios completos en el período concurrente. El relleno se realizó
correlacionando los montos mensuales y distribuyendo, posteriormente, según la
proporción observada en Cerro El Padre. El Cuadro 4.5 resume las coordenadas,
altitud y distancia entre las estaciones utilizadas.
102
La evaporación presentó unos pocos días con vacíos, los cuales fueron rellenados con
el valor de evaporación promedio mensual de julio, obtenido de la estadística
disponible. Este nivel de exactitud es suficiente, de acuerdo a lo señalado por Perrin et
al. (2003).
Cuadro 4.5: Estaciones utilizadas para la modelación hidrológica del río Huequecura.
Distancia
Altitud
Estación Tipo UTM Norte UTM Este a Cerro El
[m.s.n.m.]
Padre [km]
Cerro El Padre Pluviométrica 5814639 248065 0 400
Quilaco Evaporímetro 5825671 235734 17 225
Quillaileo Pluviométrica 5833333 264610 25 500
Adicionalmente, se estudió el resultado del modelo para los años 2006 y 2007, período
definido entre el 1 de abril de 2006 y el 12 de diciembre de 2007, con el propósito de
contrastar los dos aforos realizados el 17 y 18 marzo de 2007 y hacer el ajuste de los
parámetros de GR4J. En septiembre de 2007 se intentó realizar un tercer aforo, pero
las condiciones de escurrimiento impidieron dicha medición, debido a la mayor
velocidad y profundidad del río en ese momento. Sin embargo, la inspección visual, en
contraste con el aforo correctamente efectuado, permite estimar un caudal entre dos y
tres veces el medido en marzo de 2007. La Figura 4.5 contrasta ambos flujos, dando
cuenta de lo observado.
103
variación de menos del 1% en los resultados de caudales, lo cual indica que la
sensibilidad del modelo respecto a este parámetro es muy baja.
Figura 4.5: Contraste del caudal pasante en el tramo de estudio. Arriba el río Huequecura en marzo
de 2007, cuando se realizó aforo exitoso. Abajo el río en septiembre de 2007, cuando no se pudo
aforar.
104
adecuadamente la distribución de la precipitación. Una ponderación y distribución más
acuciosa de la lluvia está fuera del alcance de la información disponible.
Figura 4.6: Ubicación de cuenca respecto a curvas de isoevaporación (en verde). Estación Quilaco
fue ubicada geográficamente, pues código BNA no coincide con lo señalado por DGA (1987).
Como primer paso, se modeló con los parámetros determinados según ambas
funciones objetivo: Nash-Sutcliffe y Módulo. La modelación se realizó en ambos
períodos definidos previamente, antes y después de aplicar la corrección por las curvas
de isoevaporación. En adelante se hará referencia al “modelo preliminar”, en alusión a
los resultados después de la corrección de la evaporación, pero sin el ajuste de los
parámetros en base a las observaciones de caudales en terreno.
105
mensuales cercanos a la realidad. Sin embargo, se recomienda que el criterio de
elección se defina, además, por el tipo de caudales que se desea modelar. Para definir
correctamente los caudales bajos se recomienda Módulo.
El caudal máximo estimado para el año 2006, que afectó a toda la cuenca del Bio Bío,
fue de 456 [m3/s] según Nash-Sutcliffe y de 148 [m3/s] según Módulo (modelo
preliminar). La inspección visual de la caja del río, en la zona de estudio, y la cota
indicada por los lugareños desestima el primer monto. Los aforos y la inspección visual
realizada en septiembre de 2007 indican que los mínimos son 10 veces más bajos que
los determinados por Nash-Sutcliffe, contra las 5 veces que sobrestima Módulo, según
el modelo preliminar.
Por otra parte, el caudal dominante durante gran parte de 2003, 2004 y 2007, predicho
con el modelo preliminar, es marcadamente bajo, en relación a las crecidas
importantes señaladas previamente. Tal circunstancia y evidencia conducen a utilizar el
modelo definido por la función objetivo Módulo, debido a su capacidad superior para
determinar correctamente los caudales más bajos (Salvatierra, 2008).
106
Además, cabe señalar que, mientras no se tenga el resultado de la calidad de hábitat,
no es posible conocer a priori la sensibilidad de los peces a las fluctuaciones del
caudal. Sin embargo, se estima que los mayores cambios de la calidad del hábitat
ocurrirán para aquellos caudales más bajos, pues será entonces cuando sea posible
observar variaciones relevantes de esta calidad (simuladas por CASiMiR). Si bien las
crecidas representan eventos extremos, de alto estrés para las especies analizadas,
que producen los cambios geomorfológicos más importantes en los ríos, estos
significaran posiblemente estados de baja calidad de hábitat para todas las especies.
Es decir, a partir de cierta magnitud, que puede ser llamada caudal umbral, la calidad
de hábitat no variará significativamente y será marcadamente baja. Se espera que las
crecidas sean bien representadas en volumen, sin que se pueda señalar nada respecto
a las pendientes de las curvas de concentración y recesión de los hidrogramas de
crecida.
Los parámetros del modelo preliminar GR4J determinados para el río Huequecura, en
el Cuadro 4.4, corresponden a los asignados por las ecuaciones de correlación
indicadas en el Punto 2.2.1.1. El Cuadro 4.6 indica los valores extremos de x1, x2 y x3
definidos por el intervalo asignado por Salvatierra (2008) para cada uno de los IG (ver
Punto 2.2.1.1). Dentro de dicho intervalo se enmarca la incertidumbre asociada a la
estimación de los parámetros de GR4J. No se incluye el parámetro x4, correspondiente
al tiempo base del hidrograma unitario, pues posee una incertidumbre muy pequeña,
de acuerdo a lo estimado por Salvatierra (2008).
107
Cuadro 4.6: Valores extremos de parámetros analizados de acuerdo al intervalo de incertidumbre
definido por Salvatierra (2008).
Módulo
Parámetros GR4J Unidad Ec. Salvatierra Banda Sup. Banda Inf.
X1 mm 126.18 144.08 106.35
X2 mm 6.49 23.51 -8.74
X3 mm 437.03 942.95 129.88
Los valores definidos en el Cuadro 4.6 corresponden al intervalo dentro del cual los
parámetros de GR4J tienen un 95% de probabilidad de ser los “correctos”. En el caso
de una cuenca no controlada, existirá todo un espectro de caudales posibles, producto
de la combinación de estos parámetros.
Para determinar cual es la combinación correcta solo cabe tener caudales medidos in
situ con los cuales contrastar los resultados del modelo, dadas ciertas variaciones
dentro de los intervalos definidos.
Los caudales estimados mediante aforo e inspección visual (ver Figura 4.5) permiten
determinar los parámetros definidos en el Cuadro 4.7.
Cuadro 4.7: Valores calibrados de GR4J para río Huequecura. Modelo final.
Parámetros GR4J Módulo
X1 144.08
X2 -8.70
X3 140
X4 1.45
v1 70.56
v2 63.12
108
El Cuadro 4.8 contrasta los valores medidos e inspeccionados con lo modelado por
GR4J. Además, se incluye la observación de los días 11, 12 y 13 de febrero de 2007,
correspondientes a la primera campaña de terreno. En esta visita de inspección se
recorrió el río desde el tramo estudiado hasta el puente ubicado cerca de la confluencia
con el río Bio Bío.
Cuadro 4.8: Valores modelados por GR4J, en contraste a las observaciones realizadas en las tres
campañas de terreno efectuadas en 2007.
Modelado
Fecha Nº Terreno Observado
[m3/s]
11/02/2007 0.856
Caudal muy similar a lo
12/02/2007 0.819 1
observado en 2º terreno
13/02/2007 0.785
17/03/2007 0.876 0.826 [m3/s]
2
18/03/2007 0.837 0.842 [m3/s]
08/09/2007 1.353 Entre 2 y 3 veces lo
3
09/09/2007 1.326 observado en 2º terreno
La Figura 4.8 contrasta el resultado del modelo final (parámetros de Cuadro 4.7) con el
resultado del modelo preliminar, en el período de ajuste, para los parámetros de
Módulo en ambos casos. Como se observa, el ajuste permitió disminuir los caudales
mínimos sin afectar las mayores crecidas, lo cual es importante para el modelo de
hábitat. Así, será posible contrastar claramente las diferencias, en la calidad disponible
de hábitat, entre el período de lluvias y el de estiaje.
109
Finalmente, la Figura 4.9 contrasta el resultado del modelo final, obtenido a partir del
ajuste, con el modelo preliminar, en el período de modelación 2002-2004.
500
450
400
350
Caudal [m3/s]
300
250
200
150
100
50
0
1-4-06
1-5-06
1-6-06
1-7-06
1-8-06
1-9-06
1-10-06
1-11-06
1-12-06
1-1-07
1-2-07
1-3-07
1-4-07
1-5-07
1-6-07
1-7-07
1-8-07
1-9-07
1-10-07
1-11-07
1-12-07
Fecha [día]
Módulo Nash-Sutcliffe
Figura 4.7: Resultado de la modelación preliminar, en el período de ajuste, para los parámetros de
Nash-Sutcliffe y Módulo.
160
140
120
Caudal [m3/s]
100
80
60
40
20
0
1-4-06
1-5-06
1-6-06
1-7-06
1-8-06
1-9-06
1-10-06
1-11-06
1-12-06
1-1-07
1-2-07
1-3-07
1-4-07
1-5-07
1-6-07
1-7-07
1-8-07
1-9-07
1-10-07
1-11-07
1-12-07
Fecha [día]
Figura 4.8: Contraste del resultado final con el preliminar, en el período de ajuste, para los
parámetros de Módulo.
110
100
90
80
C au d al [m 3/s]
70
60
50
40
30
20
10
0
1-1-02
1-3-02
1-5-02
1-7-02
1-9-02
1-11-02
1-1-03
1-3-03
1-5-03
1-7-03
1-9-03
1-11-03
1-1-04
1-3-04
1-5-04
1-7-04
1-9-04
1-11-04
Fecha [día]
Módulo Serie Definitiva
4.5 COMENTARIOS
111
comunicación fluida de las mediciones, además, para registrar la evolución de
los caudales.
• Utilizar un sistema de aforo más elaborado.
Se consideró, en un comienzo, utilizar la cuenca alta del río Chillán para calibrar GR4J
y, posteriormente, ajustar los parámetros en base a las características geomorfológicas
comparadas de dicho río y el río Huequecura. Sin embargo, la incertidumbre en tal
caso es de la misma medida y características que la existente con la metodología de
Salvatierra (2008). No es posible resolver el problema de forma satisfactoria sin las
mediciones suficientes en terreno. Esto representa el mayor problema en proyectos de
esta naturaleza, pues la inversión debe ser proporcional a la precisión deseada.
112
Aproximar un análisis de calidad de hábitat a través de un modelo hidrológico
representa, no obstante, una innovación importante en el campo de la ecohidráulica.
Los estudios revisados y señalados en las referencias de este trabajo, corresponden a
aplicaciones en condiciones de amplia colecta y disponibilidad de información, por lo
que la necesidad de una aproximación de esta clase no ha sido estudiada.
113
5. MODELACIÓN HIDRÁULICA
Uno de los ejes principales sobre los que se sustenta la modelación de hábitat
propuesta en esta investigación corresponde a la modelación hidráulica del tramo
analizado en el río Huequecura.
5.1.1 CONSIDERACIONES
114
para su uso en HEC-RAS. Los perfiles medidos fueron enumerados de uno a ocho en
el sentido del flujo, tal como se indica en la Figura.
2550
8
2540
2530
7
2520
6
2510
Eje Y
5
2500
2490
4
3
2480
Water Level SECCIÓN
2
DERECHA
2470 Water Level SECCIÓN
IZQUIERDA
1
Eje Río
2460
4940 4950 4960 4970 4980 4990 5000 5010 5020 5030
Eje X
Figura 5.1: Configuración del río de acuerdo a los requerimientos de HEC-RAS. Flecha indica el
sentido del flujo. Coordenadas relativas en metros.
La Figura 5.2 muestra las cotas de fondo del río y del eje del río, definido en la Figura
5.1, respecto al nivel del espejo de agua medido en terreno. Como se puede ver, la
sección derecha indica un nivel del agua menor a la salida. La observación en terreno
conduce a plantear que aquí se produce crisis, sin embargo, ésta no ocupa todo el
perfil. La zona izquierda tiene un banco de sedimento importante que desvía el flujo
hacia la salida, en la sección derecha del perfil.
115
98.80
98.60
98.40
98.20
Cota [m]
98.00
97.80
Cota de Fondo
97.60
Eje Río
Water Level SECCIÓN IZQUIERDA
Water Level SECCIÓN DERECHA
97.40
1 2 3 4 5 6 7 8
Sección
La gran irregularidad del lecho del río, sumada a la escasa longitud del tramo
modelado, introduce incertidumbre respecto a la pendiente del río. Esto toma
relevancia cuando se debe definir la condición de borde de la modelación hidráulica del
tramo estudiado.
116
Como se observa, es difícil estimar una pendiente específica, dada las condiciones de
irregularidad del fondo y la escasa longitud del tramo medido.
d S1/ 6
n =
26.4
LOB (Left Overbank) y ROB (Right Overbank) corresponden a las orillas izquierda y
derecha del cauce definido a través de la topografía realizada y fueron determinados
arbitrariamente durante la modelación con HEC-RAS, manteniendo en el canal
principal el espejo de agua medido durante la campaña de marzo de 2007.
117
tanto, plantear un conjunto de coeficientes base inicial a partir del cual poder calibrar el
modelo hidráulico mediante estos coeficientes de rugosidad.
118
5.1.1.1 MÉTODO DE APROXIMACIÓN
• Crisis
119
El error promedio resultante es de un 30%, el cual es significativo respecto a las
profundidades totales registradas.
• Altura Conocida
Debido a que la altura a la salida (hS) es variable, pues cambia de izquierda a derecha
y con el caudal, se probó con diferentes valores entre ambos extremos de altura
medida durante la campaña de topografía (WL izquierdo = 98.46 y WL derecho = 98.3).
El comportamiento general fue de disminución del error a mayor nivel a la salida.
Considerando esta condición de borde solo para el caudal aforado, se corrió el modelo
para una serie de caudales, entre 0.2 y 95 [m3/s], rango en el cual está contenido el
máximo y el mínimo caudal del período de modelación. Se impuso crisis como
condición de borde para los demás caudales.
120
promedio a un 14%. El Cuadro 5.3 resume los principales valores derivados de esta
modelación.
Este método consiste en el empleo de una rutina de HEC-RAS que permite estimar
cualquiera de los parámetros contenidos en la ecuación de Manning, introduciendo los
valores de los demás. Con respecto al coeficiente de Manning, introduce varios
métodos para su estimación, pudiéndose seleccionar cualquiera de ellos, por ejemplo,
Manning (a partir de la propia ecuación), Keulegan, Strickler, Limerinos, Brownlie, etc.
(USACE, 2002).
121
La información necesaria varía según el método, aunque los parámetros básicos son
pendiente, caudal, y cota del nivel del agua (WL) en el perfil. Para utilizar Keulegan se
debe indicar la temperatura, para Limerinos se requiere el D84, para Browlie se requiere
el D84, D50 y D16.
Como se señaló en las consideraciones del Punto 5.1.1, existen al menos seis
pendientes que representan lo que sucede en el tramo. Debido a lo accidentado y corto
del tramo estudiado, se desestimó emplear pendientes entre perfiles, principalmente
por obtenerse subtramos en contrapendiente. Como primeras dos pruebas, se
presentan el Cuadro 5.4 y Cuadro 5.5. Notar la nula variación de la altura en el perfil 8
respecto a las pendientes usadas, para ambos casos. En este sentido, dicho perfil
presentará en todos los casos un elevado error a causa de su fuerte diferencia con la
observación en terreno. Por tal razón, no se considerará el perfil 8 para la estimación
del error, debido a las condiciones irrepetibles de escurrimiento que se producen en
este punto.
Las pendientes utilizadas buscan definir el intervalo [0.001; 0.005], donde se espera
encontrar el mejor resultado, de tal forma que las seis pendientes definidas en el Punto
5.1.1 pertenecen a dicho intervalo. Los resultados indican que mientras en la cota
superior (ver Cuadro 5.4) se sobrestimó la altura de escurrimiento, en la cota inferior
(ver Cuadro 5.5) el resultado fue de subestimación.
122
Cuadro 5.4: Cálculo de coeficiente de Manning. Pendiente i = 0.005. Error medio total de la
modelación: 7%.
WL HEC-RAS Coef. Manning Altura [m]
Perfil [m] Main Channel Modelada Medido Error %
1 98.66 0.167 0.86 0.77 0.12
2 98.65 0.153 0.81 0.73 0.12
3 98.63 0.661 1.14 1.07 0.06
4 98.61 0.318 0.88 0.82 0.07
5 98.55 0.325 0.74 0.72 0.03
6 98.51 0.332 0.8 0.79 0.01
7 98.45 0.158 0.5 0.55 0.08
8 98.14 0.183 0.18 0.42 0.57
Cuadro 5.5: Cálculo de coeficiente de Manning. Pendiente i = 0.001. Error medio total de la
modelación: 14%.
WL HEC-RAS Coef. Manning Altura [m]
Perfil [m] Main Channel Modelada Medido Error %
1 98.49 0.075 0.69 0.77 0.10
2 98.48 0.068 0.64 0.73 0.12
3 98.47 0.296 0.98 1.07 0.09
4 98.46 0.142 0.73 0.82 0.11
5 98.42 0.146 0.61 0.72 0.15
6 98.39 0.147 0.68 0.79 0.14
7 98.36 0.071 0.41 0.55 0.25
8 98.14 0.081 0.18 0.42 0.57
Aunque fue posible encontrar varias pendientes representativas del tramo, la mayor
incertidumbre del método es, justamente, determinar aquella más adecuada. No existe
una sola que defina lo que sucede, como tampoco existe una sola cota del nivel de
agua que sea referencia completa de lo que ocurre en el tramo. La irregularidad de
fondo plantea la imposibilidad de definir un eje de escurrimiento particular, en especial
si se considera que existe un flujo importante que se produce entre los granos más
gruesos, por lo cual la topografía del fondo del río se puede considerar como
referencial, pues no logra definir apropiadamente las pendientes de escurrimiento que
enfrenta el agua.
Una vez definido el intervalo de pendientes posibles, se realizó una búsqueda manual
de la mejor combinación de coeficientes de rugosidad que minimizara el error medio
123
total. El resultado final se presenta en el Cuadro 5.6 y la Figura 5.3. La Figura 5.4
muestra la salida gráfica de HEC-RAS para el caudal de calibración.
Cuadro 5.6: Cálculo de los coeficientes de Manning. Error medio total = 4%.
WL HEC-RAS Coef. Manning Altura [m]
Perfil [m] Main Channel Modelada Medido Error %
1 98.6 0.071 0.8 0.77 0.04
2 98.6 0.061 0.76 0.73 0.05
3 98.59 0.418 1.1 1.07 0.03
4 98.58 0.170 0.85 0.82 0.03
5 98.54 0.295 0.73 0.72 0.01
6 98.51 0.302 0.8 0.79 0.01
7 98.45 0.160 0.5 0.55 0.08
8 98.14 0.163 0.18 0.42 0.57
1.2
Modelada
1 Medido
Profundidad [m]
0.8
0.6
0.4
0.2
0
1 2 3 4 5 6 7 8
Perfil
124
en el modelo difuso de CASiMiR. Esto permite utilizar este resultado para modelar lo
que ocurre en el río durante el período de modelación definido.
3.75*
4.2* 4
4.6*
5
5.33333*
5.66666*
6.14285*
6.42857*
6.71428*
7.33333* 7
7.66666*
8
125
Cuadro 5.7: Cotas del espejo de agua para cada perfil y caudal señalado.
Caudal P E R F I L E S [m]
3
[m /s] 1 2 3 4 5 6 7 8
0.2 98.36 98.36 98.36 98.35 98.33 98.32 98.29 98.09
0.4 98.47 98.46 98.46 98.45 98.42 98.4 98.36 98.11
0.6 98.54 98.54 98.53 98.53 98.49 98.46 98.41 98.12
0.8 98.6 98.6 98.59 98.58 98.54 98.51 98.45 98.14
1 98.66 98.65 98.64 98.63 98.59 98.55 98.48 98.15
2 98.84 98.83 98.82 98.8 98.75 98.71 98.6 98.2
3 98.97 98.96 98.95 98.92 98.86 98.81 98.69 98.23
4 99.06 99.06 99.04 99.02 98.95 98.89 98.76 98.27
5 99.15 99.14 99.13 99.1 99.02 98.96 98.82 98.29
10 99.45 99.44 99.42 99.38 99.27 99.18 99.03 98.41
15 99.65 99.64 99.61 99.57 99.44 99.33 99.17 98.5
20 99.81 99.8 99.76 99.71 99.57 99.44 99.29 98.58
25 99.93 99.92 99.87 99.82 99.65 99.54 99.38 98.65
30 100.04 100.03 99.97 99.92 99.74 99.62 99.46 98.71
35 100.14 100.13 100.06 100.01 99.82 99.69 99.52 98.77
40 100.23 100.22 100.14 100.09 99.89 99.75 99.58 98.83
50 100.4 100.38 100.28 100.23 100.01 99.87 99.69 98.93
65 100.6 100.58 100.47 100.41 100.17 100.01 99.82 99.12
75 100.73 100.7 100.58 100.53 100.27 100.1 99.9 99.21
95 100.95 100.92 100.79 100.73 100.43 100.26 100.05 99.34
101.5
101 Perfiles
1
Cota Nivel del Agua [m]
100.5
2
100 3
4
99.5
5
99 6
7
98.5 8
98
97.5
0 20 40 60 80 100
Caudal [m3/s]
Figura 5.5: Curvas de descarga por perfil. Para todo el rango de caudales a modelar.
126
5.3 COMENTARIOS
El objetivo fue, en todos los casos, plantear una línea base de rugosidad y estudiar la
sensibilidad ante variaciones de ésta. El uso de una aproximación teórica como la de
Strickler parece ser una buena solución, pero es necesario estudiar más
profundamente su alcance real en esta clase de aplicaciones. Para ello se requiere,
además una cantidad y calidad de información mayor. La información granulométrica
es una limitante en este sentido, siendo la estimación a priori de la rugosidad poco
confiable. Esto significa que es necesario realizar pruebas que permitan ajustar estos
coeficientes mediante las mediciones realizadas al espejo de agua.
127
expresiones más acabadas, como Limerinos; Parker y Peterson o Ayala y Oyarce, en
cuyo caso es necesario conocer el D84 y el D90 (Niño, 2005).
El elemento de mayor incerteza resultó ser la pendiente, debido a lo irregular del fondo
y escasa longitud del tramo modelado. Además, la influencia de aguas arriba es
desconocida. Esto es independiente de que la condición de borde esté aguas abajo, ya
que se debe considerar al tramo en su contexto. En el caso estudiado, el río se
encuentra en una curva y el sector modelado se encuentra al final de un tramo recto en
el cual dicha curva desemboca. Además, el perfil de salida no evacua el agua a lo largo
de todo su ancho, debido a un banco de sedimentos ubicado en su costado izquierdo.
Esto último implica que no sea posible reproducir con un menor error lo que sucede
con la altura de escurrimiento en ese perfil (el error en ese punto es de 57%). Significa
también que es cuestionable suponer un flujo uniforme que permita emplear como
condición de borde (CB) una altura normal, lo cual condujo al uso de la crisis como CB
definitiva.
Es indudable que existen otras fuentes de incertidumbre, como la calidad de los puntos
topográficos tomados en terreno, la orientación de los perfiles transversales con
respecto al flujo, las correcciones realizadas, los puntos eliminados, etc. Cada uno de
estos elementos forma parte de las limitantes propias del modelo, además de su
aproximación 1D mediante HEC-RAS.
Sin embargo, la evaluación del modelo planteado será dentro del contexto de su
aplicación, de los resultados que se obtengan al final de la investigación, pues no es
posible hacer un juicio sin contrastar su respuesta al problema que intenta resolver.
Así, aunque las limitantes se conocen, al igual que la incertidumbre de las variables
involucradas, es posible que todo esto y no mucho más, sea suficiente para dar
solución al problema estudiado.
128
6. MODELACIÓN DE LA CALIDAD DE HÁBITAT
Para aprovechar y manejar adecuadamente los recursos naturales que dan sustento a
la vida y que cada vez son más demandados, es primordial entender la variabilidad
inherente a los procesos físicos, biológicos y químicos que los originan (Mackinson,
1999; Rogers, 2006). Sin embargo, con el tiempo se ha descubierto que los sistemas
naturales no se comportan bajo leyes simples y fáciles de deducir. Es más, se ha
reconocido su comportamiento complejo, no lineal, no estacionario y caótico
(Mackinson, 1999).
Dentro de este marco, la lógica difusa (Fuzzy Logic) se ha convertido en una técnica
innovadora para modelar procesos ecológicos debido a su capacidad para incorporar la
incerteza inherente a las variables ecológicas, además de las relaciones no lineales
que se establecen entre éstas (Adriaenssens et al., 2006; Ahmadi-Nedushan et al.,
2006).
129
realizó con el biólogo de la Universidad de Chile Leopoldo Fuentes, experto en
especies dulceacuícolas del Centro de Ecología Aplicada (CEA), con el fin de dar
forma a la información biológica que servirá de entrada al modelo CASiMiR.
El número total de reglas, para cualquier especie, está definido por todas las
combinaciones de las variables de entrada y por el número de categorías lingüísticas
que el experto emplea para calificar dichas variables, fijadas previamente. De esta
forma cada variable será calificada mediante los calificativos: HIGH (alto o grande),
MEDIUM (medio), LOW (baja o pequeño).
El experto recibe una serie de Tablas con todas las combinaciones de calificativos por
variable (27 reglas por especie y estado de vida), las cuales completa de acuerdo a su
apreciación de lo que dichas combinaciones de variables significan para la idoneidad
del hábitat. No existe un rango previo de cada variable que permita calificarlas, pues la
calificación y medición se hace en paralelo en la campaña de terreno correspondiente.
El Cuadro 6.1 muestra las reglas difusas establecidas para la especie nativa B.
australis, en la cual se aprecia el cambio de sus requerimientos de hábitat,
especialmente para velocidades bajas, para dos estados de vida diferente: adulto y
juvenil. Además se establecieron las reglas difusas para T. aerolatus, T. chiltoni, D.
130
nahuelbutaensis, P. irwini y B. maldonadoi, tanto para su estado adulto como juvenil,
las cuales se encuentran en el anexo digital que acompaña a la presente investigación.
El procedimiento para obtener los conjuntos difusos consta de una etapa en terreno y
otra donde la información es procesada para derivar dichos conjuntos difusos (ver
Fuzzy Clustering en el Anexo A).
131
Durante la campaña de terreno efectuada los días 8 y 9 de Septiembre de 2007 se
midieron las tres variables involucradas en la modelación (velocidad del río,
profundidad y tamaño representativo del sustrato en el punto de medición). Dichas
mediciones fueron acompañadas por el biólogo de peces, quien, previo a cada
medición, realizó una valoración verbal de cada una de dichas variables, sin conocer el
valor de la variable medida. El proceso llevó a seleccionar un total de 36 puntos que
abarcaran la diversidad del hábitat estudiado, siendo la única limitación la fuerte
corriente en la zona media del eje transversal del río. La Figura 6.1 muestra el proceso
de medición que se llevó a cabo y el problema asociado a la fuerte corriente en
algunas zonas.
132
concentrando en el lado izquierdo la valoración verbal de los puntos registrados en el
lado derecho.
133
6.1.2.1 AGRUPAMIENTO DIFUSO
Los resultados obtenidos a partir del agrupamiento difuso (fuzzy clustering), resumidos
en la Figura 6.2, son tres conjuntos difusos tipo 2 (ver Anexo A) correspondientes a las
tres variables utilizadas para delimitar la habitabilidad del río.
Esta aproximación se muestra en la Figura 6.3, donde es posible apreciar que cada
uno de los puntos originales es modificado de tal manera que se conserva la forma
difusa original, salvo para aquellos valores de las variables cercanos e iguales cero y
para los valores más altos definidos por el algoritmo. Esto se hace para dar sentido a
los resultados y a las valoraciones del experto. Así, si α > x > 0 tiene un grado de
pertenencia igual a uno en el conjunto difuso “bajo”, entonces ∀x ∈ X se tendrá que si
x < α el grado de pertenencia será completo al conjunto “bajo”, eliminando los mínimos
grados de pertenencia asignados a los conjuntos “medio” y “alto”. Esto es análogo para
los valores altos.
134
Resultado Fuzzy Clustering
1.2
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
Velocidad [m/s]
Low Medium High
1
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Profundidad [m]
Low Medium High
1
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
Tamaño del Grano [m]
Low Medium High
Figura 6.2: Resultado del agrupamiento difuso sobre los datos de la campaña efectuada en
septiembre de 2007.
135
1.2
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2
Velocidad [m/s]
1.2
1
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Profundidad [m]
1.2
1
Grado de Pertenencia
0.8
0.6
0.4
0.2
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
Tamaño del Grano [m]
Figura 6.3: Aproximación lineal hecha sobre el resultado del agrupamiento difuso
136
Esta aproximación de los datos es, sin embargo, la etapa uno de la fuzificación11, en la
cual se asignan las funciones de pertenencia a cada grupo. Esto permite contrastar la
información dada por el experto, al valorar las mismas variables que se emplearon para
obtener los conjuntos difusos, con lo obtenido mediante el algoritmo. La etapa dos es
ajustar los conjuntos difusos a la imprecisión del experto.
La Figura 6.4 contiene los conjuntos difusos definitivos, derivados de los conjuntos
indicados en la Figura 6.3. Los conjuntos definitivos son idénticos a los de la Figura
6.3, salvo leves diferencias, y no hacen más que simplificar su formulación matemática
al definir unos pocos puntos para delimitar cada cambio lineal del grado de
pertenencia. En esta Figura, además, se indica el intervalo de valores en los que el
experto realizó su valoración de las variables correspondientes. Estos intervalos son
obtenidos de los datos resumidos en el Cuadro 6.3.
Cuadro 6.3: Definición de índices empleados en CASiMiR para caracterizar el sustrato (Adaptado de
Schneider et al., 2004).
Índice Sustrato Asignación
0 Arcilla, Limo, Franco
1 Materia Orgánica, Detritus
2 Fango
3 Arena 0.063 - 2 mm
4 Grava Fina 2 - 20 mm
5 Grava 2 - 5 cm
6 Grava Gruesa 5 - 20 cm
7 Bolones 20 - 40 cm
8 Bolones > 40 cm
9 Roca
11
Este término viene del término inglés fuzzyfication, el cual no tiene una traducción limpia al español y que significa
la determinación de la forma y los parámetros de la función de pertenencia.
137
1
0.9
0.8
Grado de Pertenencia
0.7
Low
0.6 Medium
0.5 High
0.4 LOW-Exp
MED-Exp
0.3
HI-Exp
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
Velocidad [m/s]
0.9
0.8
Grado de Pertenencia
0.7
Low
0.6
Medium
0.5 High
0.4 LOW-Exp
MED-Exp
0.3
HI-Exp
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Profundidad [m]
1
0.9
0.8 Low
Grado de Pertenencia
0.7 Medium
High
0.6
LOW-Exp
0.5
MED-Exp
0.4
HI-Exp
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9
Sustrato
Figura 6.4: Conjuntos difusos definitivos de la etapa uno en contraste con valoración realizado por
el experto en la campaña de septiembre de 2007.
138
Es destacable la definición de “sustrato medio” del experto, ya que ésta, sin excepción,
corresponde a los puntos contenidos dentro del índice “grava gruesa”. Lo difuso
resulta, en este caso, de haber definido datos tipo “grava gruesa” una vez como
“sustrato pequeño” y en dos oportunidades como “sustrato grande”. Esta inconsistencia
(si no ocurriese sería inusual) tiene su origen en la subjetividad del experto para
realizar la definición, lo cual se traduce en cierta dispersión de la clasificación hecha.
Por otra parte, el conjunto difuso “sustrato pequeño” fue modificado considerando la
definición del experto del sustrato denominado “Roca” en el Cuadro 6.3. De acuerdo al
experto, los afloramientos de roca en el río deben ser considerados como “sustrato
pequeño”, debido a la nula protección que ofrecen a los peces. Schneider
(Comunicación directa, 2008) señala una experiencia similar respecto a la definición de
roca como “sustrato pequeño” por parte de los expertos y recomienda como solución la
modificación introducida en la Figura 6.4.
A partir de la Figura 6.4, la cercanía entre las definiciones realizadas por el experto y
las agrupaciones que hizo el algoritmo, basado únicamente en la estructura de los
datos, permite cumplir con la etapa dos de la fuzificación: ajustar los grados de
pertenencia de acuerdo a la apreciación del experto.
El conjunto “velocidad” debe ser corregido en el límite inferior y superior del conjunto
“velocidad media”, pero de manera leve. El conjunto “velocidad alta” debe ser ajustado
ligeramente, pues es más difuso para el experto que para el algoritmo, lo cual es
incuestionable, ya que la valoración humana se basa en la apreciación visual de dicha
variable.
139
• El día en que se realizó la medición coincidió con un alto caudal del río, lo cual
impidió realizar mediciones en zonas de profundidades relativamente altas.
• La turbulencia y los sedimentos arrastrados por el río dificultan la visión e
impiden discriminar certeramente las profundidades.
• La dificultad aumenta debido a lo heterogéneo del tamaño de los granos. No es
lo mismo pararse sobre una roca que sobre otra. No existe una combinación de
profundidad – tamaño de grano correcta. Esta incertidumbre es manejada por el
experto, pero no es 100% efectiva su reducción.
Los ajustes que se realizan al conjunto “profundidad” son más importantes, lo cual
aumenta la posibilidad de que los resultados no sean los esperados. Sin embargo, aún
queda una última etapa que permitirá corregir los conjuntos difusos hasta su forma
definitiva.
∑μ
i =1
ik =1
La Figura 6.5 contiene los conjuntos difusos modificados de acuerdo a la valoración del
experto. En una entrevista posterior, con algunos resultados como parte de la
conversación, Leopoldo Fuentes (Comunicación directa, 2008) no realizó
modificaciones a los conjuntos y reglas difusas. Sus comentarios son revisados a
continuación de los resultados generales obtenidos.
140
1
0.9
0.8
Grado de Pertenencia
0.7
Low
0.6 Medium
0.5 High
0.4 LOW-Exp
MED-Exp
0.3
HI-Exp
0.2
0.1
0
0 0.2 0.4 0.6 0.8 1
Velocidad [m/s]
0.9
0.8
Grado de Pertenencia
0.7
Low
0.6
Medium
0.5 High
0.4 LOW-Exp
MED-Exp
0.3
HI-Exp
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Profundidad [m]
1
0.9
0.8 Low
Grado de Pertenencia
0.7 Medium
High
0.6
LOW-Exp
0.5
MED-Exp
0.4
HI-Exp
0.3
0.2
0.1
0
0 1 2 3 4 5 6 7 8 9
Sustrato
Figura 6.5: Conjuntos difusos modificados de acuerdo a la apreciación del experto sobre las
variables medidas en cada punto.
141
6.2 RESULTADOS GENERALES
Los resultados del modelo corresponden a una variedad de índices y salidas gráficas,
explicadas en el Punto 2.1.2, que permiten estudiar y analizar la calidad de hábitat
disponible para cada especie y estado de vida. A continuación se presentan los
resultados básicos de CASiMiR para cada una de las especies y estados de vida
modelados. Además, se explica brevemente cómo se interpretan estos resultados,
dejando para los capítulos sucesivos el análisis más acabado de estos.
El formato de los resultados que se entregan es el mismo para todas las especies y
sus estados de vida. Las figuras presentadas son un resumen, cuyo contenido general
es el siguiente:
ii. Segunda Figura: indica el mapa del tramo modelado, donde los colores indican el
SI correspondiente. De esta forma es posible ver la distribución espacial de las
diversas calidades de hábitat en el tramo de estudio. Se utiliza como ejemplo un
caudal cercano al aforado (0.8 m3/s) en la campaña de marzo de 2007 (ver
Cuadro 4.8), aunque el modelo permite estudiar cualquier caudal dentro de cierto
rango previamente definido.
142
6.2.1 PERCILIA IRWINI
• Adulto
Para el adulto en la Figura 6.6 derecha se indica cómo el aumento del caudal
disminuye el porcentaje de área con condiciones medias de idoneidad (0.3 < SI ≤ 0.6)
al mismo tiempo que aumenta el porcentaje de área con condiciones bajas de
idoneidad (0 ≤ SI ≤ 0.3). También es posible notar que existe un aumento del
porcentaje de área con condiciones de idoneidad alta (0.6 < SI ≤ 1.0) para condiciones
de flujo similares a las aforadas, lo cual explica la convexidad de las curvas de WUA y
HHS en la Figura 6.6 izquierda.
En crecidas importantes, el tramo contará con escasas áreas de refugio, lo cual revela
que la zona no es apta para la especie en el período de lluvias. En el período estival
(por ejemplo, Figura 6.7) se cuenta con una amplia zona de buenas condiciones para
el adulto, lo cual permitiría su reproducción, dado que las orillas, en las mismas
condiciones, resultan aptas para los juveniles, como se ve en la Figura 6.9.
143
400 0.4 90
0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
HHS [-] 70
10 [m³/s]
HHS [-]
50
200 0.2 15 [m³/s]
40
30
100 0.1 20
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.6: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Percilia irwini. Estado de vida:
Adulto.
Simbología SI
Figura 6.7: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Percilia irwini. Estado de vida: Adulto.
144
• Juvenil
Para el juvenil en la Figura 6.8 derecha se indica cómo el aumento del caudal
disminuye las condiciones altas de idoneidad (0.6 < SI ≤ 1.0) al mismo tiempo que
aumenta el porcentaje de área con condiciones bajas de idoneidad (0 ≤ SI ≤ 0.3).
También es posible notar que existe un aumento del porcentaje de área con
condiciones de idoneidad media (0.3 < SI ≤ 0.6) para condiciones de flujo similares a
las aforadas.
La Figura 6.8 izquierda indica que, a pesar de que las condiciones medias de
habitabilidad (HHS) son altas en un comienzo, el individuo juvenil sufre un cambio de
condiciones importante debido al aumento del caudal, lo cual se ve reflejado en el
fuerte gradiente de descenso de HHS.
La Figura 6.9 muestra que las zonas de preferencia del pez en su estado juvenil se
ubicarán en las orillas, donde las velocidades son menores. Una simulación de
caudales superiores muestra que, en general, el individuo contará con algunas zonas
de refugio durante el período de lluvias. A pesar del estrés que implica una crecida, el
individuo juvenil contará con más áreas de refugio que el adulto, aunque estas estarán
acotadas a zonas muy específicas del tramo estudiado.
Existe alta complementariedad entre individuos adultos y juveniles que justificaría los
montos de capturas reportados en el sector (Leopoldo Fuentes, Comunicación directa,
2008). No obstante, la bibliografía consultada no señala los períodos de reproducción
de la especie. Como se ve en el siguiente capítulo, esta especie cuenta con
información parcial acerca de sus estados de vida, lo que dificulta el análisis de los
resultados.
145
90
500 0.5 0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
60 10 [m³/s]
HHS [-]
300 0.3 50 15 [m³/s]
200 0.2 40
30
100 0.1 20
10
0 0
0
0 10 20 30 40 50 60 70 80 90 100
0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s]
Idoneidad [-]
Figura 6.8: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Percilia irwini. Estado de vida:
Juvenil.
Simbología SI
Figura 6.9: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Percilia irwini. Estado de vida: Juvenil.
146
6.2.2 DIPLOMYSTES NAHUELBUTAENSIS
• Adulto
Los resultados indicados en la Figura 6.10 muestran que el individuo adulto se siente
considerablemente más cómodo con caudales relativamente grandes. No obstante, los
resultados obtenidos en el Capítulo 4 indican que las condiciones de flujo base no son
adversas para esta especie, tal como se confirma en la Figura 6.11. A pesar de que
caudales bajos como los de la temporada de estiaje no representan el óptimo, las
zonas habitables poseen un SI considerable, lo cual indica que en la poca área
disponible, las condiciones son buenas.
Esto significa que el cambio en los patrones de velocidad y profundidad para pequeñas
crecidas no afectaría mayormente al individuo adulto. Un aumento mayor, sin embargo,
reduce las zonas óptimas, debido a que, en las zonas más profundas, las velocidades
aumentan. La profundidad, por su parte, no es limitante, pues ha sido reportada a
varios metros de profundidad (Manuel Contreras, Comunicación directa, 2007). En las
orillas, el sustrato es inadecuado para la protección del individuo, pues el grano es
demasiado pequeño para el tamaño del pez.
147
80
600 0.4 0.2 [m³/s]
WUA [m²] 70 0.8 [m³/s]
500
400 10 [m³/s]
50
HHS [-]
15 [m³/s]
300 0.2 40
200 30
0.1 20
100
10
0 0
0
0 10 20 30 40 50 60 70 80 90 100
0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.10: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Diplomystes nahuelbutaensis.
Estado de vida: Adulto.
Simbología SI
Figura 6.11: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Diplomystes nahuelbutaensis. Estado de vida: Adulto.
148
• Juvenil
Para el juvenil en la Figura 6.12 derecha se señala que las condiciones óptimas se
encuentran en mayor proporción de área con el caudal más bajo. Un aumento del
caudal lleva las condiciones hacia un nivel medio.
Se observa que el individuo juvenil puede afrontar pequeñas crecidas sin que esto
produzca desmejoras importantes en la habitabilidad del tramo estudiado. La curva de
WUA indica que la relación entre calidad y cantidad de hábitat será máxima dentro de
un rango que varía entre los 4 y los 10 [m3/s]. En este rango, de acuerdo a la Figura
6.12 derecha, un aumento de caudal producirá una disminución de las condiciones
promedio debido al aumento de área inundada con SI bajas, sin embargo, también
aumentarán las zonas de SI medias, lo cual permitirá que esta especie se adapte bien
a las condiciones de estrés en el período de lluvias.
La Figura 6.13 señala que las zonas de mayor calidad de hábitat (SI) para el individuo
juvenil se encuentran en las riberas, evidenciando complementariedad con el individuo
adulto, tal como ocurre en el caso de la P. irwini. Esto concuerda con los reportes de
captura (Leopoldo Fuentes, Comunicación directa, 2008), que indican una mayor
frecuencia de captura de individuos juveniles. Debido a que las pescas en el tramo de
estudio se realizaron solo con instrumentos de pesca eléctrica, las zonas de captura se
reducen únicamente a las orillas (Leopoldo Fuentes, Comunicación directa, 2008).
El cambio de gradiente de descenso de HHS alrededor de los 25 [m3/s] indica que las
zonas con SI > 0 se mantienen con valores altos de este índice. La disminución de la
WUA es menos pronunciada entonces, pues aunque el aumento del área inundada
conlleva a un aumento de zonas de bajo SI, persisten aún áreas de resguardo,
principalmente a causa del sustrato y bajas profundidades, que evitan el aumento
excesivo de la velocidad.
149
80
500 0.5 0.2 [m³/s]
WUA [m²] 70 0.8 [m³/s]
400 0.4
10 [m³/s]
50
300 0.3
HHS [-]
15 [m³/s]
40
200 0.2 30
100 0.1 20
10
0 0
0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.12: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Diplomystes nahuelbutaensis.
Estado de vida: Juvenil.
Simbología SI
Figura 6.13: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Diplomystes nahuelbutaensis. Estado de vida: Juvenil.
150
6.2.3 TRICHOMYCTERUS AREOLATUS
• Adulto
Para el adulto en la Figura 6.14 se indica que la condición general del tramo estudiado
es media y para caudales muy bajos la zona es poco apta para la especie. Las mejores
condiciones de habitabilidad se encuentran entre los 5 y los 15 [m3/s], con un máximo
claro a los 10 [m3/s]. Esto indica que la especie se encontrará con la mejor calidad
disponible durante la temporada de lluvias, mientras que en la época estival la
habitabilidad se verá disminuida. La habitabilidad también será mínima en crecidas de
un gran período de retorno.
La Figura 6.15 señala las zonas de mejor SI disponibles, los cuales se ubican en las
zonas de mayor profundidad. Las reglas difusas permiten entender que las condiciones
óptimas de la especie se encuentran en valores medios, siendo poco adecuado la
combinación de velocidades altas y medias en conjunto con profundidades altas. Los
mapas de habitabilidad indican que la calidad diminuirá en el centro del cauce a
medida que aumente el caudal, desplazando las áreas aptas hacia los costados, pero
nunca en las zonas próximas a las riberas. De esta forma, se observa que el individuo
adulto tiene requerimientos mucho más específicos que los demás casos analizados,
razón por la cual, la habitabilidad media del tramo, representada por el HHS, es la más
baja de las especies estudiadas.
151
90
400 0.3 0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
70
0.2 60 10 [m³/s]
HHS [-]
50 15 [m³/s]
200
40
0.1 30
100
20
10
0 0
0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.14: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Trichomycterus areolatus. Estado
de vida: Adulto.
Simbología SI
Figura 6.15: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Trichomycterus areolatus. Estado de vida: Adulto.
152
• Juvenil
La Figura 6.17 señala que las áreas de SI > 0 cubren casi toda el área inundada,
siendo en su mayoría de valor medio. Esto es indicativo de la diversidad de patrones
que el individuo juvenil es capaz de tolerar en cierto grado. Las condiciones, si bien son
medias, poseen una cobertura importante que puede ser justificado por los
antecedentes de captura (Leopoldo Fuentes, Comunicación directa, 2008), donde la
mayoría de los individuos capturados corresponden a juveniles.
153
90
400 0.4 0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
10 [m³/s]
HHS [-]
50 15 [m³/s]
200 0.2
40
30
100 0.1
20
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.16: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Trichomycterus areolatus. Estado
de vida: Juvenil.
Simbología SI
Figura 6.17: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Trichomycterus areolatus. Estado de vida: Juvenil.
154
6.2.4 TRICHOMYCTERUS CHILTONI
• Adulto
Para el adulto en la Figura 6.18 se indica que la condición general del tramo estudiado
es media y para caudales muy bajos la zona es poco apta para la especie. Las mejores
condiciones de habitabilidad se encuentran entre los 5 y los 15 [m3/s], con un máximo
claro a los 10 [m3/s], según la WUA. Esto indica que la especie se encontrará con la
mejor calidad disponible durante la temporada de lluvias, mientras que en la época
estival la habitabilidad se verá disminuida. La habitabilidad también será mínima en
crecidas de un elevado período de retorno.
Al comparar los resultados de T. chiltoni (Figura 6.18 y Figura 6.19) con los resultados
de T. areolatus (Figura 6.14 y Figura 6.15), en conjunto con las reglas difusas
correspondientes, es posible apreciar una gran similitud. Ambas especies poseen
preferencias muy similares, sin embargo, la bibliografía sobre T. chiltoni es
relativamente escasa y su presencia en Chile es también mucho más limitada. Los
resultados de capturas indican que es posible encontrar a la especie en diferentes
épocas del año, sin embargo, la bibliografía y la cantidad de capturas es insuficiente
para diferenciar entre individuos adultos y juveniles. Esto limita en forma importante el
análisis de los resultados.
La Figura 6.19 señala las zonas de mejor SI disponibles, los cuales se ubican en las
zonas de mayor profundidad. Las reglas difusas permiten entender que las condiciones
óptimas de la especie se encuentran en valores medios, siendo poco adecuado la
combinación de velocidades altas y medias en conjunto con profundidades altas. Los
mapas de habitabilidad indican que la calidad disminuirá en el centro del cauce a
medida que aumente el caudal, desplazando las áreas aptas hacia los costados, pero
nunca en las zonas próximas a las riberas.
155
Respecto a T. areolatus, la especie T. chiltoni demuestra tener menor sensibilidad a las
fluctuaciones de caudal, lo cual eleva los valores de HHS y WUA respecto a la primera
especie. Aunque ambas parecen alcanzar un óptimo en el período de lluvias, T. chiltoni
aparece como mejor adaptado a las condiciones del tramo analizado.
90
500 0.3 0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
400
0.2 60 10 [m³/s]
300
HHS [-]
50 15 [m³/s]
200 40
0.1 30
100 20
10
0 0
0
0 10 20 30 40 50 60 70 80 90 100
0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s]
Idoneidad [-]
Figura 6.18: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Trichomycterus chiltoni. Estado de
vida: Adulto.
Simbología SI
Figura 6.19: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Trichomycterus chiltoni. Estado de vida: Adulto.
156
• Juvenil
La Figura 6.21 señala que las áreas de SI > 0 cubren casi toda el área inundada,
siendo en su mayoría de valor medio. Esto es indicativo de la diversidad de patrones
que el individuo juvenil es capaz de tolerar dentro de cierto rango. Debido a que los
antecedentes de captura no permiten diferenciar entre individuos adultos y juveniles,
no es posible analizar con mayor profundidad los resultados. No obstante, es esperable
que los individuos capturados sean en su mayoría juveniles.
157
400 0.4 90
0.2 [m³/s]
WUA [m²] 80
0.8 [m³/s]
10 [m³/s]
HHS [-]
50
200 0.2 15 [m³/s]
40
30
100 0.1 20
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.20: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Trichomycterus chiltoni. Estado de
vida: Juvenil.
Simbología SI
Figura 6.21: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Trichomycterus chiltoni. Estado de vida: Juvenil.
158
6.2.5 BULLOCKIA MALDONADOI
• Adulto
La Figura 6.22 derecha indica que los individuos adultos tienen una preferencia
elevada por condiciones de poco escurrimiento. La época con mejores condiciones de
habitabilidad será en el período estival. No obstante, existe una presencia permanente
de condiciones medias de habitabilidad, lo cual, sumado a las condiciones altas,
permite una zona de habitabilidad aceptable durante la época de lluvias. Alrededor de
los 15 [m3/s] las zonas de preferencia se restringen casi únicamente a las orillas, sin
embargo, aquí las condiciones continúan siendo altas, por lo cual la WUA mantiene
valores altos respecto a las demás especies.
159
80.0
600 0.6 0.2 [m³/s]
WUA [m²] 70.0 0.8 [m³/s]
500 0.5
HHS [-]
15 [m³/s]
300 0.3 40.0
Figura 6.22: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Bullockia maldonadoi. Estado de
vida: Adulto.
Simbología SI
Figura 6.23: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Bullockia maldonadoi. Estado de vida: Adulto.
160
• Juvenil
161
80
600 0.6 0.2 [m³/s]
WUA [m²] 70
500 0.5 0.8 [m³/s]
HHS [-]
300 0.3 40 15 [m³/s]
200 0.2 30
20
100 0.1
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.24: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Bullockia maldonadoi. Estado de
vida: Juvenil.
Simbología SI
Figura 6.25: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Bullockia maldonadoi. Estado de vida: Juvenil.
162
6.2.6 BASILICHTHYS AUSTRALIS
• Adulto
De los casos analizados, esta especie es la que posee menor capacidad de soportar
los eventos extremos del río Huequecura. Como indican los caudales modelados, el río
presenta fuertes crecidas, al menos una vez al año y, en el caso de un término
temprano de las lluvias, caudales muy pequeños en verano.
La época de desove comenzaría en agosto, lo cual coincidiría con las últimas lluvias
del año y montos menores de precipitación respecto al período mayo – julio. De
acuerdo a los resultados, esta especie requiere una zona cuyo flujo base sea elevado,
lo cual no ocurre en el río Huequecura.
163
400 0.3 100
90 0.2 [m³/s]
WUA [m²]
0.8 [m³/s]
0.2
60 10 [m³/s]
HHS [-]
200 50 15 [m³/s]
40
0.1
100 30
20
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.26: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Basilichthys australis. Estado de
vida: Adulto.
Simbología SI
Figura 6.27: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Basilichthys australis. Estado de vida: Adulto.
164
• Juvenil
La Figura 6.28 señala que el individuo juvenil posee muy buena adaptación a
condiciones de escurrimiento muy bajas e incluso acepta cierto rango de crecidas.
Alrededor de los 5 [m3/s] las condiciones de habitabilidad decrecen rápidamente, sin
embargo, la WUA decrece a una menor tasa que la HHS debido a que las pocas áreas
disponibles con caudales elevados, poseen una idoneidad alta.
Los datos de captura no permiten comprobar los resultados obtenidos, pues la especie
no fue detectada en el río Huequecura. Los antecedentes revisados ubican a la
especie en la zona alta del río Queuco y en el río Bio Bío, alrededor de la
desembocadura del río Queuco (Leopoldo Fuentes, Comunicación directa, 2008) y
aguas abajo de la desembocadura del río Huequecura. Se trata de una especie
vulnerable y poco abundante en la cuenca del Bio Bío, como se concluye ante la
escasez de capturas reportadas.
165
100
600 0.6 0.2 [m³/s]
WUA [m²] 90
0.8 [m³/s]
500 0.5 80
HHS [-]
15 [m³/s]
300 0.3 50
40
200 0.2
30
100 0.1 20
10
0 0 0
0 10 20 30 40 50 60 70 80 90 100 0<=SI<=0.3 0.3<SI<=0.6 0.6<SI<=1.0
Caudal [m3/s] Idoneidad [-]
Figura 6.28: Evolución de los índices WUA y HHS con el caudal (izquierda) y la evolución de la
calidad de hábitat disponible con el caudal (derecha). Especie: Basilichthys australis. Estado de
vida: Juvenil.
Simbología SI
Figura 6.29: Mapa de la calidad del hábitat en el tramo de estudio para un caudal de 0.8 [m3/s].
Especie: Basilichthys australis. Estado de vida: Juvenil.
166
6.3 VALORACIÓN EXPERTA DE LOS RESULTADOS
1. Desde el punto de vista del experto, las valoraciones plasmadas en las reglas
difusas se relacionan con su experiencia en varios ríos, por lo cual intentan ser
dependientes de condiciones específicas y de zonas particulares. No es
posible, por lo tanto, definir las condiciones generales del tramo estudiado como
correctas o incorrectas. Se busca únicamente evaluar una zona o punto de
muestreo por las condiciones físicas locales observadas en ella.
167
en otros ríos. Actualmente se utilizan, en Chile, curvas de preferencia de peces
generadas en el río Laja (Parra et al., 2000) para otros ríos en los cuales se
encuentran las mismas especies, sin embargo, se desconoce o no se considera
hasta qué punto la especie está adaptada a las condiciones de variabilidad
hidrológica del río en cuestión12. En este sentido, Victoriano et al. (2003)
presentan un estudio de variabilidad genética donde dos de las tres especies
analizadas poseen curvas de preferencia en el estudio de Parra et al. (2000) y,
a su vez, la mayor variabilidad genética, incluso en la misma cuenca, a medida
que se desciende por ella. Es significa que, dependiendo de la especie, es
posible encontrar adaptaciones de una misma especie a condiciones locales de
un río, lo cual obligaría al menos a considerar este factor antes de emplear las
curvas de preferencia generadas en otra cuenca o tramo del río.
12
Existe evidencia de que la variabilidad de los caudales juega un rol fundamental en la estructura de las comunidades
acuáticas en sistemas lóticos. Para más información ver Poff y Ward (1989).
168
tramo o tipo geomorfológico analizado. Debido a la conectividad natural de los
sistemas lóticos, los individuos que se encuentren una zona determinada
deberán su presencia no solo a las condiciones del área de permanencia, sino
de las condiciones en la zona en la cual se encontraban previamente. Por lo
tanto, se debe considerar seriamente ampliar la complejidad espacial y
temporal de esta clase de estudios. Es posible que existan desfases temporales
de información durante su traslado espacial. En otras palabras, el individuo
acumula una historia que lo lleva a un sitio en un momento determinado. Con
los estudios actuales, se ignora el grado de importancia que tiene esta historia
en los resultados de la modelación del hábitat realizada. Al menos se requiere
validar fuertemente una aproximación simplificada de la ecología del río, como
la que representa la presente investigación.
6.4 COMENTARIOS
169
es ajustar la valoración experta en los conjuntos definitivos y que estos sean lineales.
En la presente investigación se empleó un método de agrupamiento difuso que no está
justificado en la bibliografía estudiada, pues no existen antecedentes acerca de las
metodologías empleadas en cada parte de la modelación con CASiMiR. De este modo,
el método empleado representa una alternativa no validada de agrupamiento, que es
relativamente sencilla y rápida, una vez estudiada la teoría.
Respecto a los resultados obtenidos hay varios aspectos importantes que requieren
mayor profundización en estudios posteriores. Por ejemplo, la convexidad conjunta de
HHS y WUA, observada en algunas especies (gráficos no incluidos), refleja cierto
grado de histéresis, pues la idoneidad no es la misma, para un mismo caudal, entre
una crecida y la recesión de ésta. Sin embargo, con los datos existentes no es posible
entender si representa lo que sucede o es una particularidad del modelo. Es posible
que el uso de CASiMiR pudiese ayudar a dilucidar esta cuestión, pero esto es solo una
suposición que requiere mayor análisis.
170
(Tello, 1999). La historia de tales bifurcaciones es una dimensión imprescindible para
comprender el sistema y sus cambios de estado. Estudiarlos supone identificar
aquellos parámetros cuya acumulación o combinación pueden desbordar los sucesivos
umbrales de adaptación y resiliencia (Tello, 1999).
Por lo tanto, las capturas de un individuo no deben verse justificadas únicamente por lo
que sucede en la zona de estudio. Debido a la conectividad natural del río, los
individuos que se encuentren una zona específica deberán su presencia tanto a las
condiciones del área de permanencia, como a las condiciones en la zona en la cual se
encontraban previamente. La historia del sistema puede ser importante, por lo cual se
requiere ahondar en este aspecto de la investigación.
171
Los datos de captura se consideran como un antecedente que permite profundizar el
análisis, pero los resultados y comentarios de las condiciones de habitabilidad,
contrastadas con las capturas, no son consideradas definitivas.
172
7. ANÁLISIS Y DISCUSIÓN DE RESULTADOS
En este capítulo se busca profundizar en el análisis de los resultados por medio del
procesamiento de la información disponible. Se intenta verificar si la información es
suficiente para encontrar la relación entre la variabilidad hidrológica – fluvial y la
abundancia de una especie. También, se indica brevemente la relación encontrada
entre la biodiversidad de cada campaña de captura y el caudal medio mensual
modelado. Finalmente, y en base al trabajo realizado, se intenta una forma de definir
caudal “ecológico”, tanto en la forma como en el fondo. Para esto se entrega una breve
referencia de casos y antecedentes existentes en la literatura, como medio de respaldo
de las indicaciones entregadas.
El Cuadro 7.1 representa un ejemplo de los datos obtenidos por el CEA, donde la
Estación 7 corresponde al tramo en estudio. Como se observa, la información no
detalla el estado de vida del individuo, solo sus dimensiones, peso y estado gonadal.
Sobre esto último, el procedimiento es abrir a los peces, por lo cual solo es permitido
hacerlo en 5 individuos cada vez, de los cuales, algunos darán como resultado
“indeterminado” (Leopoldo Fuentes, Comunicación directa, 2008). Esto significa que es
173
escasa la información que se puede obtener de cada campaña respecto al estado de
vida de los individuos capturados, haciendo en la práctica muy difícil una separación
entre adultos y juveniles.
Cuadro 7.1: Datos de captura, mayo de 2002. L.T.: Largo total. P.T.: Peso total. G: Desarrollo
gonadal. IND: Indeterminado. M(n): Macho (estado desarrollo gonadal).
Especie L.T. P.T G Estación
(cm) (g)
Diplomystes nahuelbutaensis 14.5 48.0 M4 3
16.0 46.1 M2 4
5.3 1.9 IND 7
Trichomycterus areolatus 3.7 0.4 IND 7
3.7 0.4 IND 7
3.7 0.4 IND 7
4.3 0.5 IND 7
3.3 0.3 IND 8
4.6 0.7 IND 8
4.8 0.8 IND 8
Trichomycterus chiltoni 3.2 0.2 IND 7
4.6 0.7 IND 7
Percilia irwini 3.5 0.6 IND 7
3.5 0.6 IND 7
3.6 0.6 IND 7
3.7 0.6 IND 7
3.8 0.7 IND 7
3.8 0.8 IND 7
4.0 0.9 IND 7
4.4 1.2 M 7
174
Cuadro 7.2: Información recopilada sobre relación entre tallas y estados de vida de las especies dulceacuícolas nativas estudiadas
(elaboración propia).
Tamaño Máximo [cm]
Nombre Científico Nombre Común Bibliografía
Juvenil Adulto
En el río Laja, la longitud total promedio de los individuos adultos fue de 149.1±24.6mm, con
valores límite de 110.8 y 217.0mm. En los juveniles, el promedio fue de 54.4±27.3mm y fluctuó
entre 12.9 a 101.9mm (Habit, 2005). Existe, además, varianza de tamaño, para un mismo
Diplomystes nahuelbutaensis Tollo 10.2 30
estado de vida según época del año, estación de muestreo y cuenca (Habit, 2005). La longitud
se correlaciona significativa y positivamente con la velocidad de la corriente, profundidad
máxima y temperatura (Habit, 2005).
Bullockia maldonadoi Bullockia 6 - -
Especímenes juveniles de 7.6 mm a 25 mm de longitud total prefieren zonas ribereñas, de 0.5
a 10 mm de profundidad, con fondo barroso, con pequeñas piedras, escasa o nula vegetación y
aguas quietas o lentas (Arratia, 1983). Por otra parte, especímenes juveniles de 25 a 60 mm de
longitud total se distribuyen de preferencia en zonas próximas a las riberas, entre la vegetación
acuática. Parte de los individuos se ubican entre las rocas de los torrentes o zona de aguas
Trichomycterus aerolatus Pez gato ó Bagre 6 15
muertas, según Arratia (1983), cuya profundidad varía entre los 5 y 10 cm. Finalmente, los
individuos de 70 a 130 mm de longitud total prefieren la zona de aguas muertas de los torrentes
en profundidades de hasta 40 cm, de fondo pedregoso y con una velocidad de agua superficial
que varía entre los 0.1 a 0.3 m/s en diferentes ríos y esteros, en diferentes épocas del año
(Arratia, 1983)
Trichomycterus chiltoni Pez gato ó Bagre - - -
Son pequeños peces, de alrededor de 7 cm de longitud total, no llegan a sobrepasar los 10 cm,
Percilia irwini Carmelita 5 10
siendo el promedio 5-6 cm (CONAMA)
Basilichthys australis Pejerrey Chileno - - Su talla puede ser mayor a 30 cm. de longitud total (CONAMA).
175
Cuadro 7.3: Información sobre épocas de desove y presencia de juveniles de las especies dulceacuícolas nativas estudiadas (elaboración
propia).
Nombre Científico Nombre Común Época de Desove Hábitat de los Juveniles
Con los caudales más bajos del año, es decir, entre fines del
Diplomystes nahuelbutaensis Tollo
verano y comienzos de otoño (Habit, 2004).
Ocupan las zonas ribereñas de
Ocurre entre fines de noviembre hasta principios de diciembre escasa profundidad y de velocidad de
Bullockia maldonadoi Bullockia (CONAMA). De acuerdo a Habit y Victoriano (2004), su corriente casi nula, entre octubre y
reproducción se produce en primavera-verano. abril, siendo el período más frecuente
entre diciembre y abril (Arratia, 1983).
176
Como se observa en el Cuadro 7.1, solo es posible diferenciar los estados de vida
(adulto y juvenil) de cuatro de las seis especies estudiadas. No existen antecedentes,
dentro de la literatura revisada, que permitan reconocer los estados de vida de las
demás especies, lo cual imposibilita comparar sus resultados de captura con lo
obtenido mediante CASiMiR. El Cuadro 7.3 ordena la información conocida sobre
momentos importantes para la sobrevivencia de las especies, lo que facilita su
contraste con los resultados obtenidos.
Existen métodos directos e indirectos para determinar la edad de una especie. Dentro
de los primeros generalmente se recurre a muestras de escamas, otolitos, espinas, etc.
Dentro de los segundos el más corriente es la interpretación de la frecuencia de tallas o
método de Petersen, que se basa en el análisis de datos de frecuencia por clases de
longitud, y asume que los sucesivos valores modales de la distribución representan
distintos grupos de edades (UBA, 2004). Existen métodos más complejos y certeros,
177
como el método de Bhattacharya, pero su uso requiere de herramientas
computacionales que no están disponibles en forma libre (FAO, 1997).
178
artificialmente el número de grupos. Sin embargo, es evidente que existe al menos un
gran grupo, cuyo peak se ubica en el intervalo (4; 4.3]. Se encontraron cinco grupos o
valores modales, pero esto no significa necesariamente que dichos grupos existan,
pues la aparición de estos depende del tamaño de los intervalos fijados previamente y
en forma arbitraria, sin embargo, lo importante es que se puedan diferenciar
claramente, lo cual indica que la división está correctamente hecha.
Capturas 2000-2005
0.25
Fracción de la Captura Total
0.20
0.15
0.10
0.05
0.00
1 1.6 2.2 2.8 3.4 4 4.6 5.2 5.8 6.4 7 7.6 8.2 8.8
Largo Total [cm ]
Figura 7.1: Distribución de frecuencias normalizadas por talla de P. irwini, en el período 2000-2005.
Capturas 2000-2005
0.25
Fracción de la Captura Total
0.20
0.15
0.10
0.05
0.00
1
8
1.
2.
2.
3.
4.
5.
5.
6.
7.
8.
8.
179
De acuerdo a Ojeda (2004), se han encontrado antecedentes de peces cuya edad de
madurez ha variado, ya sea por causas naturales o por intervención humana, pero la
talla se ha mantenido. Por ejemplo, en casos donde la pesca excesiva aumenta la
mortalidad de los especímenes, los sobrevivientes alcanzan una talla mínima de
madurez más deprisa por un crecimiento más rápido, madurando a una edad más
temprana (Ojeda, 2004). Esto significa que el resultado de la Figura 7.2 representa
más un estado que una edad, pues no se conocen antecedentes del período que toma
a P. irwini alcanzar la madurez. De acuerdo a Habit y Belk (2007), esta especie alcanza
un máximo de cuatro años de edad, con el más rápido crecimiento durante el primer
año. En definitiva, resulta claro que existe un importante estado de vida alrededor de
los 4 [cm] de talla, lo cual podría esta ocurriendo entre el primer y segundo año de vida.
Estos antecedentes llevan a fijar la talla de madurez en 4 [cm], la mínima talla a la cual
se ha detectado desarrollo gonadal (ver Cuadro 7.4).
Por ende, se cuenta con cuatro especies con las cuales hacer un análisis más
profundo de los resultados, de las cuales solo se estudiarán tres: T. areolatus, D.
nahuelbutaensis y P. irwini. El caso de B. maldonadoi se desestimó en esta tesis,
debido a sus preferencias de hábitat similares con T. areolatus (ver Puntos 6.2.3 y
6.2.5), pero puede ser encontrado en el trabajo de Meza & Vargas (2009), que se
incluye en el anexo digital de esta tesis.
180
El Cuadro 7.5 resume los principales resultados de las tres especies que se
presentarán a continuación.
7.2.1 T. AREOLATUS
Estas figuras permiten tener una idea general de lo que sucede en las temporadas de
lluvias y estiaje, pero son muy generales. Para resumir los resultados y contrastarlos
con los datos de captura disponibles, se presenta la Figura 7.5.
181
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
-01
0
0.1
0.2
0.3
-01
0
0.1
0.2
0.3
-01
-04 -03 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
HHS
HHS
05 06 06
HHS
Día
Día
05 06 06
-08 -08 -08
-04 -03 -02
Caudal
Caudal
Caudal
05 06 06
-09 -09 -09
-04 -03 -02
06 07 07
-10 -10 -10
-04 -03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
10
10
10
0.1
0.1
0.1
100
100
Figura 7.3: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de estudio.
182
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
0.4
0
0.1
0.2
0.3
0.4
-01 -01
0
0.1
0.2
0.3
0.4
-01
-04 -03 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
HHS
HHS
HHS
05 06 06
-07 -07 -07
-04 -03 -02
Día
Día
Día
05 06 06
-08 -08 -08
-04 -03 -02
Caudal
Caudal
Caudal
05 06 06
-09 -09 -09
-04 -03 -02
06 07 07
-10 -10 -10
-04 -03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
10
10
10
0.1
0.1
0.1
100
100
Figura 7.4: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de estudio.
183
Trichomycterus Areolatus % Presencia Juvenil
% Presencia Adulto
0.35 120
Habitat Juvenil
0.3 100 Habitat Adulto
0.25
80
% Captura
HHS [-]
0.2
60
0.15
40
0.1
0.05 20
0 0
02
03
04
02
03
04
02
03
04
2
4
-0
-0
-0
v-
v-
v-
b-
b-
b-
o-
o-
o-
ay
ay
ay
no
no
no
fe
fe
fe
ag
ag
ag
m
m
Fecha
Figura 7.5: Presencia detectada de T. areolatus en relación a la calidad media de hábitat en el tramo
de estudio.
Empleando los tamaños incluidos en el Cuadro 7.2, se separó entre adultos y juveniles
el total de individuos capturado en la zona de estudio, lo cual permitió generar las
barras de la Figura 7.5, las que indican, en términos de porcentaje, la presencia de
individuos capturados en el mes correspondiente, respecto al total de adultos y
juveniles capturados en esa fecha. Tal como se señaló anteriormente, no se cuenta
con la fecha exacta, por tanto, se incluye el parámetro HHS, como promedio mensual,
en forma de líneas.
184
principalmente en la orilla. De acuerdo al Cuadro 7.3, esta es la zona principal
que habitan los juveniles, entre octubre y abril. Las capturas, sin embargo,
señalan individuos juveniles fuera de este período, el año 2002. Los años 2003 y
2004, por otra parte, corroboran lo indicado en la literatura.
iii. No se observa un patrón entre la calidad promedio modelada y la presencia
detectada a través de las capturas. El año 2002 fue muy lluvioso, sin embargo,
solo se capturaron individuos juveniles, fuera del período esperado según la
bibliografía consultada. Los dos años posteriores entregan resultados que son
esperables, pero debido a la baja cantidad y frecuencia de capturas durante los
tres años, una conclusión sería injustificada.
185
7.2.2 D. NAHUELBUTAENSIS
186
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
0.4
-01
0
0.1
0.2
0.3
0.4
0
0.1
0.2
0.3
0.4
-01 -01
-04 -03 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
05 06 06
-07
HHS
-07 -07
HHS
HHS
Día
Día
Día
05 06 06
-08 -08 -08
-04 -03 -02
05 06 06
Caudal
-09
Caudal
-09
Caudal
-09
-04 -03 -02
06 07 07
-10 -10 -10
-04 -03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
10
10
0.1
0.1
0.1
100
100
Figura 7.6: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de estudio.
187
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
0.4
0.5
0
0.1
0.2
0.3
0.4
0.5
-01 -01
0
0.1
0.2
0.3
0.4
0.5
-01
-04 -03 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
05 06 06
HHS
-07
HHS
-07 -07
-04 -03 -02
Día
Día
Día
05 06 06
HHS
-08 -08 -08
-04 -03 -02
05 06 06
Caudal
Caudal
Caudal
-10 -10 07
-04 -10
-03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
1
1
10
0.1
0.1
0.1
100
100
Figura 7.7: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de estudio.
188
Diplomystes nahuelbutaensis % Presencia Juvenil
% Presencia Adulto
0.45 120
Habitat Juvenil
0.40
100 Habitat Adulto
0.35
0.30 80
% Captura
HHS [-]
0.25
60
0.20
0.15 40
0.10
20
0.05
0.00 0
02
03
04
02
03
04
02
03
04
2
4
-0
-0
-0
v-
v-
v-
b-
b-
b-
o-
o-
o-
ay
ay
ay
no
no
no
fe
fe
fe
ag
ag
ag
m
m
Fecha
i. Existe un dominio claro de presencia juvenil en el tramo, sin embargo, esto puede
ser atribuido casi por completo al arte de pesca utilizado. Aun así, la presencia se
registra en períodos de una habitabilidad media – alta, siendo la presencia del
adulto en febrero de 2004 la única inconsistencia o variabilidad detectada.
ii. Al igual que en el caso de T. areolatus, existe una marcada idoneidad del hábitat
para los individuos juveniles por sobre los adultos. No obstante, en períodos muy
lluviosos esta tendencia se puede revertir.
iii. Aunque las condiciones de habitabilidad del tramo son levemente mejores para
D. nahuelbutaensis que para T. areolatus, el número total de capturas de la
primera, en el mismo período, es menor en un 25%.
Aunque los resultados de capturas parecen mostrar cierta consistencia en este caso,
12 campañas en 3 años no es estadísticamente representativo, sumado a la baja
abundancia detectada de la especie. Por otra parte, la escasez de información impide
responder a través del modelo si la abundancia encontrada se debe a las fluctuaciones
del caudal.
189
7.2.3 P. IRWINI
Para resumir los resultados y contrastarlos con los datos de captura disponibles, se
presenta la Figura 7.11. En ella es posible observar una variedad de porcentajes de
captura mucho mayor, pues la abundancia de esta especie supera a todas las
demás13.
13
Es el punto de muestreo del CEA que posee mayor abundancia de esta especie.
190
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
0.4
-01
0
0.1
0.2
0.3
0.4
-01
0
0.1
0.2
0.3
0.4
-03 -01
-04 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
HHS
HHS
HHS
05 06 06
-07 -07 -07
-04 -03 -02
Día
Día
Día
05 06 06
-08 -08 -08
-04 -03 -02
Caudal
Caudal
Caudal
05 06 06
-09 -09 -09
-04 -03 -02
06 07 07
-10 -10 -10
-04 -03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
1
10
0.1
10
10
100
0.1
0.1
100
100
Figura 7.9: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de estudio.
191
Idoneidad [-] Idoneidad [-] Idoneidad [-]
01 01 01
0
0.1
0.2
0.3
0.4
0.5
0
0.1
0.2
0.3
0.4
0.5
-01 -01
0
0.1
0.2
0.3
0.4
0.5
-01
-04 -03 -02
01 01 01
-02 -02 -02
-04 -03 -02
03 04 04
-03 -03 -03
-04 -03 -02
03 04 04
-04 -04 -04
-04 -03 -02
04 05 05
-05 -05 -05
-04 -03 -02
04 05 05
-06 -06 -06
-04 -03 -02
HHS
HHS
05 06 HHS
-07 -07 06
-07
-04 -03 -02
Día
Día
Día
05 06 06
-08 -08 -08
-04 -03 -02
Caudal
Caudal
Caudal
05 06 06
-09 -09 -09
-04 -03 -02
06 07 07
-10 -10 -10
-04 -03 -02
06 07 07
-11 -11 -11
-04 -03 -02
07 08 08
-12 -12 -12
-04 -03 -02
1
1
10
10
10
0.1
0.1
0.1
100
100
100
Figura 7.10: Variación diaria del caudal y la calidad promedio del hábitat (HHS) en la zona de
192
Percilia irwini % Presencia Juvenil
0.45 120 % Presencia Adulto
0.40 Habitat Juvenil
100 Habitat Adulto
0.35
0.30 80
% Captura
HHS [-]
0.25
60
0.20
0.15 40
0.10
20
0.05
0.00 0
02
03
04
02
03
04
02
03
04
2
4
-0
-0
-0
v-
v-
v-
b-
b-
b-
o-
o-
o-
ay
ay
ay
no
no
no
fe
fe
fe
ag
ag
ag
m
m
Fecha
Figura 7.11: Presencia detectada de P. irwini en relación a la calidad media de hábitat en el tramo
de estudio.
193
aspecto debe ser parte de la discusión en la aplicación de modelos de calidad de
hábitat como CASiMiR.
iii. Aunque las condiciones generales de calidad se encuentran en la misma fase, no
sucede lo mismo con la presencia, donde se observa un desfase claro entre
adultos y juveniles. La abundancia y la frecuencia de captura en el total de
campañas permite determinar que los juveniles se encuentran en mayor
proporción en el período de estiaje, mientras que los adultos dominan en el
período de lluvias. En esta temporada es posible que los individuos adultos se
desplacen hacia las orillas debido al aumento de la escorrentía, razón por la cual
son capturados a través de la pesca eléctrica. Esta aparente contradicción revela
que el ciclo natural de la especie está adecuado a las mejores condiciones
disponibles para la sobrevivencia de los juveniles, pues los peaks de habitabilidad
coinciden con los peaks de presencia de éstos. El Cuadro 7.2 puede
complementar esta observación, pues la mayor población de P. irwini se
encuentra en torno a la edad de madurez sexual, mientras que los demás
estados cuentan con muy pocos individuos detectados.
iv. Dentro del período de modelación, el año 2002 presenta las peores condiciones
de hábitat, lo cual se ve confirmado por la menor presencia de individuos
capturados. Los años 2003 y 2004 cuentan con condiciones medias de
habitabilidad relativamente altas durante todo el año, lo cual explicaría la
presencia constante de la especie en casi todas las campañas de captura.
Los resultados asociados a P. irwini resultan ser los más consistentes, en gran medida
gracias a la abundancia y frecuencia de capturas. En este caso, las 12 campañas
resultan ser suficientes para obtener un patrón que, sin embargo, no es explicado del
todo por el modelo. En este sentido, es necesario comprender más acerca de la
biología del pez y la función evaluativa del modelo en un dominio poco comprendido
como el de la ecología de los peces nativos en el río Huequecura. No obstante, se
deben analizar las zonas de captura respecto a los mapas de habitabilidad. En este
ámbito, los resultados han sido consistentes para todas las especies revisadas.
194
No resulta útil interpretar los resultados de CASiMiR como predictivos, pues no es su
objetivo ni su función. Sin embargo, la presente aplicación permite conocer
antecedentes que pueden complementar o guiar los estudios acerca de la biología de
P. irwini y de las demás especies, de ahí la llamada función evaluativa del modelo.
Además, permite identificar las zonas del río en que su presencia tiene mayor
posibilidad de ocurrencia.
0.45
HHS Promedio del Tramo en el Mes
0.40
0.35
0.30
de Captura [-]
0.25
0.20
0.15
0.10
0.05
0.00
0 20 40 60 80 100
% de Presencia de Juveniles sobre el Total de las Capturas
Figura 7.12: Relación entre la proporción de de juveniles sobre todas las capturas de la especie y la
modelación de CASiMiR, a través del HHS.
195
Si la Figura 7.12 fuese hecha en base a suficientes capturas, estaría indicando que el
modelo es capaz de modelar adecuadamente la calidad del hábitat para la P. irwini
juvenil, ya que a mejores condiciones modeladas (HHS) mayor es la presencia
detectada.
Por otra parte, el concepto de caudal “ecológico” resulta ambiguo, pues no existe una
definición clara de parte de los diferentes organismos, en Chile, que tienen
competencia en esta materia (CONAMA, DGA, DOH). Ormazábal (2004) presenta
estas definiciones y señala como ejemplo el caso del embalse Ralco, donde Endesa
presentó una propuesta de caudal ecológico, para la aprobación ambiental del proyecto
a través del Sistema de Evaluación de Impacto Ambiental (SEIA), basada en criterios
seleccionados ad hoc y con datos de preferencias de peces nativos distorsionados14. El
caudal mínimo aconsejable propuesto por Endesa fue de 6.1 [m3/s], el cual fue
rechazado por la autoridad, estableciéndose el mínimo histórico registrado en la
zona15: 39.3 [m3/s] (Ormazábal, 2004).
14
En referencia a D. nahuelbutaensis, T. areolatus y otras especies Endesa señaló en su EIA (Estudio de Impacto
Ambiental) que estas “especies son de tamaño pequeño, con características morfológicas de desarrollo longitudinal, lo
cual las hace especialmente adaptables a aguas someras” (Ormazábal, 2004).
15
El criterio de la autoridad se basó en la ausencia de información biológica que sustente una metodología diferente a
la utilizada por la DGA, sin que esto impida, en base a nuevos antecedentes, modificar el caudal en el futuro
(Ormazábal, 2004).
196
En Chile, el caudal “ecológico”, también llamado reservado, recomendado, regulado,
medioambiental, mínimo aconsejable, etc. (Palau, 1994; DGA, 2000), se aplica sólo en
el tramo comprendido entre el muro del embalse y la zona de restitución de las aguas.
Dicha zona, que en el caso de Ralco alcanza los 9 [km] (Ormazábal, 2004), busca ser
protegida de ser secada por completo, lo cual es conseguido por criterios de base
hidrológica, en ausencia de más información. Sin embargo, no se considera dentro del
área de influencia del proyecto lo que sucede luego de la restitución, donde el régimen
hidrológico se ve fuertemente afectado. Aquí se producen descargas importantes en el
período de estiaje, cuando los embalse de generación hidroeléctrica tienen un mayor
uso.
Con las exigencias de la Ley de Bases del Medio Ambiente 19.300, promulgada en
1994, el Estado chileno busca que los proyectos de inversión importantes, por ejemplo,
los embalses, se sometan al SEIA de tal modo que incorporen los impactos
ambientales del proyecto dentro de sus costos, antes de la ejecución del proyecto, bajo
el criterio del principio precautorio: siempre será más barato prevenir los impactos
ambientales que reparar el daño ya hecho. No se desea desincentivar la inversión
privada través de exigencias mayores, sino que se consideren dentro de tales
inversiones aquellas externalidades que afectan a terceros y que, por falta de
información, son ignoradas o desestimadas.
197
Los elementos esenciales relacionados con los ecosistemas fluviales (y otros) que
deben ser considerados en todo proyecto de importancia han sido descritos y
abordados desde diferentes enfoques en distintas partes del mundo. Un ejemplo
atingente es la metodología IFIM. Por lo tanto, sólo cabe mencionar aquellos
elementos que, desde el punto de vista de la presente investigación, deben ser
reforzados, tomando en cuenta el punto 1 de los elementos a considerar: proyecto =
oportunidad.
Estos tres elementos básicos deben ser parte de las medidas de mitigación de
proyectos de inversión importantes que afecten a los ríos, como ocurre con los
embalses de generación hidroeléctrica. Si se considera una zona como prístina o poco
intervenida, la necesidad de estudios tendrá un carácter puramente científico. No
obstante, al intervenir o afectar este lugar para fines productivos o de sobrevivencia,
tales estudios se convierten en una necesidad social básica. Para asegurar la
sustentabilidad y rentabilidad social que el país necesita, se requiere que el costo del
proyecto considere tal necesidad. No es excesivo, entonces, considerar, como parte de
la mitigación al impacto del proyecto, los recursos necesarios para realizar las
investigaciones en las tres áreas básicas mencionadas.
198
Atendiendo lo anterior, es posible cambiar el paradigma del caudal “ecológico”, para lo
cual se propone el concepto “régimen sustentable de caudales”.
199
su operación. El tramo intermedio, entre la barrera y la restitución, representa un
desafío mayor, si se da énfasis a la zona posterior a la restitución, pues el agua
turbinada en las centrales de embalse de Chile solo siguen la orden de producción del
Centro de Despacho Económico de Carga (CDEC) (McPhee, 2009), ignorando
variables como época del año y el régimen hidrológico del río. Sin duda, a falta de más
información, se pueden continuar utilizando metodologías de base hidrológica para
estimar un cierto caudal ecológico, pero sin omitir la existencia de cierta variabilidad de
caudales, aunque sea amortiguada respecto a lo que es habitual en el río analizado.
En este sentido, el proyecto hidroeléctrico del río San Pedro (EULA, 2007a) representa
un ejemplo de régimen sustentable de caudales, al igual que varias iniciativas de
hidroelectricidad en desarrollo actualmente, donde las empresas solo poseen derechos
no consuntivos sobre los río, lo cual les impide embalsar las aguas, pero no darles
altura de caída mediante el uso de muros que atraviesen el cauce de los ríos.
Se debe permitir replicar las condiciones naturales del sistema fluvial, en lo referente a
la variabilidad temporal y de volumen. Respecto a las cantidades, la presente
investigación demuestra, en forma preliminar, que a partir de cierto monto de caudal,
todas las especies han alcanzado un grado de calidad de hábitat máximo.
200
la medida que los estudios financiados por el proyecto entreguen nuevos
antecedentes. En tal sentido, la validación requiere de un monitoreo constante durante
un período de al menos 10 años posterior a la entrada en vigencia del proyecto (Hardy,
2009). Idealmente se requiere información de las mismas características en el período
previo al embalsamiento, aunque esto es más difícil de obtener (Hardy, 2009).
Esto significa que el embalse hipotético en el río Huequecura podrá aprovechar toda el
agua por sobre este máximo, que fluctúa entre los 10 y los 15 [m3/s], para ser
embalsada y, además, podrá liberar hasta dicho monto máximo en caso de generar
electricidad durante el período de lluvias (abril – agosto), si es un año seco. Si las
condiciones naturales lo permiten, se puede optar entre almacenar agua o generar más
electricidad en ese momento. En este sentido es esencial generar un modelo
hidrológico predictivo que permita tomar decisiones año a año, a través de un manejo
adaptativo de los recursos embalsados.
Como caudal mínimo se debe permitir el monto definido por la especie más exigente,
que en el caso de este estudio corresponde a los bagres D. nahuelbutaensis, T.
areolatus y T. chiltoni. El monto mínimo es de 0.5 [m3/s], pero en caso de que las
condiciones sean de mayor estrés, se puede liberar menos, igualando el caudal
afluente al embalse. El objetivo es permitir el almacenamiento aún en condiciones de
bajo caudal, aumentando ligeramente la tasa de descenso del flujo y manteniendo
estable el caudal en 0.5 [m3/s], mientras el río llega naturalmente a este caudal. Si el
caudal natural desciende bajo este límite, no se podrá almacenar nada.
Para acumular agua, el embalse puede alterar, bajo ciertas condiciones, la variabilidad
natural introduciendo una nueva, similar en montos a la natural. Si se requiere
almacenar agua en condiciones de ausencia de lluvias, se podrá hacer durante un
201
período definido por la especie menos tolerante a períodos de estrés. El aumento
posterior deberá ser gradual y acorde a la variabilidad afluente al embalse. La limitante
relacionada con los derechos de agua existentes puede ser resuelta mediante el
arriendo de derechos, por el período en que se desee almacenar el agua.
Aunque se pueden aprovechar las crecidas por sobre los 10 o 15 [m3/s] para
almacenar, se debe permitir que dicho monto sea excedido cada cierto tiempo, si las
condiciones naturales lo permiten, con un tope máximo de 40 [m3/s], si aún queda
capacidad de almacenamiento ociosa. Se podrá liberar más de 40 [m3/s] si no se desea
almacenar más agua, siempre que las condiciones naturales impongan esta condición.
Por lo tanto, se requiere profundizar en los estudios de tolerancia al estrés por períodos
prolongados, tanto de altos caudales como bajos. Para ello se propone estudiar la serie
histórica de caudales e identificar los períodos de estrés y frecuencias asociadas a los
caudales menores a 0.5 [m3/s] y superiores a 15 [m3/s]. La presencia de las especies
en el río Huequecura es indicadora de que dichas condiciones son tolerables, mientras
que los argumentos de Poff y Ward (1989) incluso justifican esto.
202
Finalmente, es esencial comprender de qué forma se verá alterado el equilibrio
sedimentológico de la cuenca producto de una barrera hipotética en el cauce del río
Huequecura, debido al muro del embalse, para buscar las medidas de mitigación
necesarias para no alterar el hábitat acuático aguas abajo de éste. En este sentido, el
proyecto de restauración del río Trinity (TRRP) y su programa de restauración de la
dinámica aluvial (Gaeuman, 2007), es un ejemplo de la relevancia que este elemento
tiene para los ecosistemas acuáticos en ríos.
Estas medidas permitirán que se pueda establecer una operación clara y específica
para un embalse hipotético ubicado en el río Huequecura. Será flexible de acuerdo al
año hidrológico en curso, pero única para el río en cuestión.
Para extender los resultados a otras cuencas, es necesario validar las necesidades de
las especies a través de estudios biológicos, en conjunto con la modelación del hábitat.
Es posible que la variabilidad genética de éstas no sea relevante y que las condiciones
requeridas en una zona sean muy similares en otras cuencas cercanas.
7.4 COMENTARIOS
Una de las limitantes más importantes que existen hoy en día para estudiar la ecología
de los ecosistemas fluviales es la relacionada con la ausencia de información biológica
acabada sobre un número importante de especies nativas dulceacuícolas.
Estos datos corresponden a la clase de estudios que Endesa financia para mantener
un monitoreo de las especies dulceacuícolas, mientras que los estudios cuantitativos
se basan principalmente en las curvas de preferencia generadas por el centro EULA
para el proyecto hidroeléctrico Quilleco (2000).
203
Sin embargo, la autoridad continúa validando las metodologías de base hidrológica,
pues está conciente que la información es insuficiente para recurrir a otras
metodologías (Ormazábal, 2004). En este sentido, CASiMiR busca ser una alternativa,
considerando que la mayor cantidad de conocimientos, en Chile, se encuentra en
manos de los biólogos expertos en esta área, y no en la literatura. Sin embargo, los
resultados demuestran que la necesidad de validarlos continúa existiendo, en igual
medida que la requerida por PHABSIM y por cualquier otro modelo basado en el
mismo concepto.
De acuerdo a los resultados, las abundancias de las especies son insuficientes, lo cual
no se podrá corregir aumentando el esfuerzo de muestreo durante las campañas de
terreno. Existe una baja densidad de varias especies, lo cual debe ser compensado,
por ejemplo, aumentando la frecuencia de las campañas.
El caso del río San Pedro (EULA, 2007b) es un ejemplo de manejo sustentable de los
recursos hídricos, pues evita la regulación del régimen hidrológico del río al actuar
204
como una central de pasada. Sin embargo, no todos los embalses podrán cumplir con
una restricción de tal magnitud, en especial si poseen derechos de aprovechamiento
que les permitan retener el agua un tiempo determinado.
205
8. CONCLUSIONES
206
8.1 VISIÓN MULTIDISCIPLINARIA EN EL MANEJO DE LOS RECURSOS
HÍDRICOS
Uno de los ejes que motivan esta investigación se orienta a abordar la gestión de los
recursos hídricos desde un enfoque multidisciplinario. Aunque no se explicita en los
objetivos, la aplicación y análisis de CASiMiR para evaluar la calidad del hábitat de
peces nativos busca, a través de la práctica y el estudio de su estructura y origen,
introducir al lector en el contexto mundial de la investigación desarrollada durante las
últimas dos décadas en torno al recurso agua.
Como todo campo de investigación, los nuevos desarrollos tecnológicos han permitido
la entrada de nuevas visiones en lugares antes impensados. Lo que antes era una
disciplina de carácter técnico, ligado casi totalmente al uso de herramientas
matemáticas, hoy ha cambiado para siempre. Nunca más existirá un solo campo de
acción para una problemática cualquiera. En el tema de los recursos hídricos existen
economistas, abogados, ecólogos, geólogos, ingenieros, geógrafos, sociólogos,
antropólogos, políticos y una larga lista de profesionales cuyas especializaciones se
cruzan entre si, eliminado la posibilidad de enmarcar el problema dentro de un solo
sector.
Durante la revisión bibliográfica se buscó abordar una serie de aspectos que, desde
afuera, pueden parecer inconexos. Desde modelos matemáticos de hábitat, modelos
conceptuales de investigación ecológica, pasando por legislación, investigación
biológica sobre especies nativas, modelación hidrológica e institucionalidad. De cada
elemento se busca una sola cosa: definir el estado del arte correspondiente.
207
La investigación de otros trabajos y autores en el área de la modelación de hábitats
acuáticos es la que definió los temas revisados. De cada uno se puede extraer un tema
de investigación nuevo y poco explorado, pero como mosaico presentan los elementos
que conllevan al manejo integral y adaptativo de los recursos hídricos. Cualquier
ausencia de uno de estos elementos debilita y limita la posibilidad de hacer un uso
responsable y sustentable del agua.
i. En Chile existe un modelo de desarrollo neoliberal que exige cierta rigidez del
sistema político-legal, tal que asegure jurídicamente los derechos de propiedad
para dar confianza a los inversionistas. Esta falta de flexibilidad no permitirá una
adaptación rápida y eficaz en un contexto de escasez de agua futura, por lo cual
todo se resolverá por la fuerza de la crisis.
ii. La falta o reserva de información, en general, acerca de la ecología de las
especies afectadas en proyectos de inversión, de la disponibilidad de recursos
hídricos y de los procesos geomorfológicos de las cuencas intervenidas le da
poca credibilidad a los resultados de investigaciones en dichos campos. Esto,
asociado a una cultura chilena de escasa innovación, impide que se invierta en
mejorarlos.
iii. No existe, hoy en día, un modelo de desarrollo sustentable basado únicamente
en un actor que fija las directrices y las ejecuta. Los autores revisados
concuerdan con esta conclusión. Es más, el modelo chileno de mercados del
agua (como forma de gestión de los recursos hídricos) es visto hoy en día como
algo que hay que evitar, más que seguir (Bauer, 2009). Es imprescindible generar
plataformas de actores que se relacionen entre si, vertical y transversalmente. La
sociedad civil organizada puede gestar cambios que conlleven a un cambio de
políticas en las autoridades gobernantes. Los científicos pueden educar a la
sociedad, dar respuestas a necesidades de la comunidad. Es por esto que la
labor de la ciencia debe estar completamente ligada a su contexto social.
208
Cómo concluye la OECD (2005), se debe fortalecer la educación y conciencia
ambientales. Se necesita una estrategia de aprendizaje ambiental de largo plazo y un
plan nacional de educación ambiental. Esto debe incluir la integración de los temas
ambientales en los planes de estudios de los colegios y el desarrollo del conocimiento
sobre el medio ambiente, tanto en las asociaciones de profesionales, como en los
sistemas de gestión ambiental de las empresas públicas y privadas (e.g., ISO 14.001).
209
los parámetros de ajuste del modelo GR4J en cuencas pluviales no controladas de
Chile central. Estos IG se estiman a través de ciertas relaciones basadas en una lista
de parámetros geomorfológicos definidos por Salvatierra (2008). A partir de estos
índices es posible obtener los parámetros del modelo GR4J para la cuenca no
controlada, en este caso, el río Huequecura en Tío Lolo.
Salvatierra (2008) estima que sólo dos funciones (Módulo y Nash-Sutcliffe) son
consistentes con las observaciones, sin embargo, cada una da énfasis a un aspecto
diferente en el ajuste de las regresiones. De esta forma, mientras el Módulo tiene una
buena estimación de los caudales más bajos, Nash-Sutcliffe permite realizar un buen
ajuste de los montos medios mensuales, si lo que se busca es tener una buena
estimación de los recursos hídricos totales.
En el presente estudio se busca modelar lo que sucede a nivel diario, con un mayor
énfasis en caudales medios y bajos, pues son los de mayor frecuencia y relevancia
para los peces. De esta forma se selecciona el criterio de Módulo.
Si los resultados de Salvatierra fuesen otros, el problema solo llegaría hasta este
punto, sin embargo, como en toda aproximación, existe cierta incertidumbre.
Salvatierra (2008) pudo establecer un límite de incertidumbre para cada uno de los
parámetros de GR4J, dentro del cual la variación de los caudales modelados no es
despreciable. En este caso la única solución es contrastar el modelo con mediciones
en terreno.
Una situación similar ocurre con la caracterización del eje hidráulico del río, para una
serie de caudales diferentes, necesaria para el modelo de calidad de hábitat. Además
210
de la topografía del río y el caudal, se requieren datos del eje hidráulico o nivel del
espejo de agua, para distintos caudales. En general se recomienda que estas
mediciones no sean menos de tres, midiendo al menos un caudal bajo y otro alto
(Hardy, 2009). Otras alternativas están ligadas a la caracterización de la granulometría
del río, permitiendo el uso de relaciones empíricas para determinar el coeficiente de
rugosidad de Manning.
La campaña de terreno realizada en marzo de 2007 entregó los datos para ajustar el
modelo hidrológico, pues en base a ésta se puede estimar el caudal pasante en las
otras dos campañas realizadas. Esta misma campaña permite obtener los
antecedentes de caudal, rugosidad y altura de escurrimiento (topografía) empleados en
el ajuste del modelo hidráulico.
En este punto, sin embargo, los datos de captura están incompletos. En varios casos
se desconoce el día exacto de la campaña de terreno, lo cual restringe la evaluación
de los resultados únicamente a nivel mensual. Surgen interrogantes que solo pueden
ser examinadas desde la experiencia y la teoría, a falta de más información de terreno.
211
8.3 OBJETO Y MÉTODO: REPLICABILIDAD DEL ESTUDIO
Las artes de pesca u observación utilizadas deben ser evaluadas, pues cada una
introducirá un sesgo que debe ser considerado al evaluar los resultados. En la
presente investigación se cuenta con peces capturados en las orillas, debido al uso de
pesca eléctrica. Existen antecedentes, entre otros, del uso de redes de arrastre,
espineles y observación directa con el uso de cámaras y personas calificadas.
El período en que se realizan las colectas es otro aspecto importante. ¿Se desea
estudiar la presencia a lo largo del año?, ¿la respuesta ante períodos de lluvias o
sequía?, ¿las zonas de desove? Cada pregunta tendrá una metodología que se adapte
mejor a sus requerimientos, cambiando el método de observación, el período y
frecuencia de ésta. Si no se conoce mucho acerca de la ecología de la especie,
entonces estudios exploratorios, como la presente investigación, son adecuados, pues
son relativamente económicos y permiten individualizar los vacíos que se deben
rellenar.
En este sentido, las capturas de P. irwini indican que el modelo puede responder o
guiar hacia nuevas preguntas. Por ejemplo, el desfase de la presencia detectada en
relación a la compatibilidad de hábitats de adultos y juveniles no es una respuesta
esperada, pero puede ser explicada una vez que se analiza junto con los antecedentes
biológicos y las preferencias. Es decir, lo que puede ser evidente, desde el punto de
vista de la ecología, no es necesariamente bien entendido por los biólogos, a menos
que se pueda tener un enfoque diferente sobre la misma materia, como el que ofrece
CASiMiR a través de sus mapas de habitabilidad e índices de calidad.
212
Aunque no existen antecedentes sobre el período de reproducción de P. irwini, si se
tiene sobre Percilia gillissi, especie que pondría sus huevos entre noviembre y enero,
mientras que las larvas aparecerían entre enero y marzo (EULA, 2007b). Durante este
período se observó, en el tramo de estudio, un descenso de la presencia de individuos
adultos y un aumento de la presencia de juveniles. Sin embargo, no se tienen
antecedentes de preguntas como ¿ocupan los juveniles los mismos nichos que huevos
y larvas? o ¿se retiran los adultos de las zonas de desove una vez efectuado éste? El
modelo no podrá responder preguntas de esta clase, pues no pretende dar respuestas
ontológicas de las especies estudiadas.
Desde este punto de vista, es posible replicar el uso del modelo en otros tramos y ríos,
pero su validez debe ser corroborada en base a las observaciones. Estas tienen la
facultad única de darle sentido a la convención asumida respecto a los índices y mapas
que entrega como respuesta el modelo.
Las observaciones necesarias son, en este caso, aquellas que permitan describir los
ciclos de vida de los peces estudiados, las áreas que ocupan en cada etapa, etc. Sin
embargo, esta información debe ser adquirida desde el enfoque metodológico a través
del cual se busca describir las relaciones entre tales ciclos y áreas. Tal cosa no existe
en la cantidad suficiente y adecuada, por lo cual se puede decir que la replicabilidad
está limitada en cuanto a nuevas observaciones.
213
8.4 RÉGIMEN SUSTENTABLE DE CAUDALES
Es desde esta perspectiva que se emplean los resultados para elaborar una propuesta
de caudal “ecológico”. En la forma que se construye, el término “régimen sustentable
de caudales” surge del conjunto de acciones directas e indirectas que se ejecutan
antes y durante la ejecución de un proyecto para reducir al mínimo el impacto de éste
sobre los ecosistemas. En este caso se puede hablar de un “régimen sustentable de
operación”. El objetivo es quitar el acento de la mera consecuencia y ponerlo también
en las acciones necesarias para que el resultado sea el mejor posible.
214
El caso del proyecto hidroeléctrico en el río San Pedro (EULA, 2007a) representa un
ejemplo de operación en base a las reglas de la misma naturaleza, al igual que en el
caso de las centrales de pasada. Será la variabilidad la que indique el caudal pasante y
la capacidad de producción hidroeléctrica. Sin embargo, muchas centrales son
construidas con el objetivo de cubrir necesidades de demanda en horas punta, por lo
cual deben regular el río. ¿Pueden los ecosistemas resistir una regulación parcial de
las aguas? ¿Hasta qué punto?
Los resultados indican, en forma preliminar, que es posible permitir una regulación
parcial. Pero también se requieren medidas de mitigación, por ejemplo, trasportando
sedimentos gruesos aguas abajo del muro (Gaeuman, 2007), liberando crecidas
importantes que movilicen el lecho, tal como ocurre naturalmente, a través de la
translocación de algunas especies para suplir la falta de libre desplazamiento (Habit et
al., 2002), etc.
215
asegurar la sustentabilidad y el beneficio social que el país necesita y persigue en el
desarrollo de nuevos proyectos, se requiere que sus costos consideren tal necesidad.
La propuesta de mitigación al impacto del proyecto, señala como aspecto esencial los
recursos necesarios para realizar las investigaciones en las tres áreas básicas
mencionadas. La propuesta hecha en este estudio es un medio para lograr el manejo
integral, adaptativo y sustentable de los recursos hídricos del país.
216
REFERENCIAS
Abonyi, J.; Babuska, R.; Szeifert, F. (2001). “Fuzzy Modeling with Multivariate
Membership Functions: Gray-Box Identification and Control Design”. IEEE Transactions
on Systems, Man, and Cybernetics – Part B: Cybernetics, Vol. 31, No. 5, pp. 755-767.
Arthington, A.; Rall, J.; Kennard, M.; Pusey, B. (2003). “Environmental Flow
Requirements of Fish in Lesotho Rivers using the DRIFT Methodology”. River Research
and Applications, Vol. 19, No. 5-9, pp. 641-666.
217
Babuska, R. (1998). “Fuzzy Modelling for Control”. Boston, Kluwer Academic
Publishers. 257 pp.
Bauer, C. (2005). “In the Image of the Market: the Chilean Model of Water Resources
Management”. International Journal of Water, Vol. 3, No. 2, pp. 146-165.
Bauer, C. (1998b). “Slippery property rights: Multiple water uses and the neoliberal
model in Chile, 1981-1995”. Natural Resources Journal Vol. 38, No. 1, pp. 109-155.
Bezdek, J.; Keller, J.; Krisnapuram, R.; Pal, N. (1999). “Fuzzy Models and Algorithms
for Pattern Recognition and Image Processing”. Editorial Springer. 776 pp.
Borsányi, P; Alfredsen, K.; Harby, A.; Ugendal, O.; Kraxner, C. (2004). “A Meso-
scale Habitat Classification Method for Production Modelling of Atlantic Salmon in
Norway”. Hydroécologie Appliquée, Tomo 14, pp. 119.
Bratrich, C.; Truffer, B.; Jorde, K.; Markard, J.; Meier, W.; Peter, A.; Schneider, M.;
Bernhard, W. (2004). “Green Hydropower: A New Assessment Procedure for River
Management”. River Research and Applications Vol. 20, No. 7, pp. 865-882.
Campos, H., J. F. Gavilán, F. Alay & V. Ruiz (1993). Comunidad íctica de la hoya
hidrográfica del río Bío Bío. En: Evaluación de la calidad del agua y ecología del
sistema limnético y fluvial del río Bío Bío. Monografía EULA 12:249-278
218
Centro de Ecología Aplicada (CEA) (2001). “Estudio de Migraciones de Peces en el
Río Biobío, Central Ralco”. Estudio encargado por Endesa. Informe entregado por
Manuel Contreras L.
Díaz, B.; Moniche, L.; Morillas, A. (2005). “Identificación de Sectores Clave Mediante
Análisis de Conglomerados Difuso. El Caso de España”. En: I Jornadas de Análisis
Input-Output, Oviedo, España. Grupo de Investigación en Análisis Input-Output y la
Universidad de Oviedo.
219
Diez Hernández, J. M. (2005). “Análisis comparativo de los métodos actuales de
simulación hidráulica en PHASIM-IFIM y su influencia en la evaluación del hábitat
fluvial: propuesta de nuevos modelos hidráulicos 1D”. ProQuest-Information and
Learning España. ISBN: 0-496-97826-8. 428 pp.
Emery, J.; Gurnell, A.; Clifford, N.; Petts, G.; Morrissey, I.; Soar, P. (2003).
“Classifying the Hydraulic Performance of Riffle-Pool Bedforms for Habitat Assessment
and River Rehabilitation Design”. River Research and Applications, Vol. 19, pp. 533-
549.
220
Fredes, M. (2000). “El Derecho de Acceso a la Información”. Ambiente y Desarrollo,
Vol. XVI, No. 4, pp. 16-24.
Galvão, C. O. (1999). Introdução à teoria dos conjuntos difusos. En: Carlos de Oliveira
Galvão; Mêuser Jorge Silva Valença. (Org.). Sistemas inteligentes; aplicações a
recursos hídricos e ciências ambientais. 1 ed. Porto Alegre: Editora da Universidade
Federal do Rio Grande do Sul e Associação Brasileira de Recursos Hídricos, Vol. 1, pp.
167-191.
García, A. (2006). “Aquatic Habitat Modeling of Chilean Native Fish in a Reach of the
Biobío River”. Thesis presented in partial fulfillment for the Requirements for the Degree
of Master of Science, University of Idaho.
Habit, E. & Belk, M. (2007). “Threatened fishes of the world: Percilia irwini Eigenmann,
1927 (Perciliidae)”. Environmental Biology of Fishes, Vol. 78, pp. 213-214.
Habit, E.; Belk, M.; Tuckfield, C.; Parra, O. (2006). “Response of the Fish Community
to Human-induced Changes in the Biobío River in Chile”. Freshwater Biology, Vol. 51,
pp. 1-11.
Habit, E.; Dyer, B.; Vila, I. (2006). “Estado de Conocimiento de los Peces
Dulceacuícolas de Chile”. En Gayana, Vol. 70, No. 1, pp. 100-112.
221
Habit, E.; Parra, O.; Valdovinos, C. (2005). “Ictiofauna de un Sistema Fluvial Receptor
de Aguas Servidas: Respuesta a una Nueva Planta de Tratamiento (Río Quilque, Chile
Central)”. Gayana, Vol. 69, No.1, pp. 94-103.
Harby, A.; Alfredsen, K.T.; Fjeldstad, H.P.; Halleraker, J.H.; Arnekleiv, J.V.;
Borsányi, P.; Flodmark, L.E.W.; Saltveit, S.J.; Johansen, S.W.; Vehanen, T.;
Huusko, A.; Clarke, K.;. Scruton, D.A (2001). “Ecological Impacts of Hydro Peaking in
Rivers”. En: hydropower in the new millennium: Proceedings of the 4th International
Conference on Hydropower Development, Hydropower '01, Bergen, Norway, 20-22
June 2001, Taylor & Francis.
Hardy, T. (2009). “Master Course: Field Techniques, Hydraulic Modeling, and Habitat
Assessments”. Apuntes. Universidad de Concepción.
Hering, D.; Moog, O.; Sandin, L.; Verdonschot, P.F.M. (2004). “Overview and
Application of the AQEM Assessment System”. Hydrobiologia, Vol. 516, pp. 1-20.
Hirji, R. & Panella, T. (2003). “Envolving Policy Reforms and Experiences for
Addressing Downstream Impacts in World Bank Water Resources Projects”. River
Research and Applications, Vol. 19, pp. 667-681.
Ibáñez, C. (2001). “El Impacto Ambiental de los Trasvases: El Caso del Ebro”. En: El
Plan Hidrológico Nacional a Debate. Colección Nueva Cultura del Agua. pp. 245-253.
222
Instituto Geográfico Militar (IGM) [200-]. Cartas electrónicas en formatos SIG escala
1:50.000 VIII Región.
Iza, A. & Stein, R. (Eds) (2009). RULE – Reforming Water Governance. Gland,
Switzerland: IUCN.
Jorde, K. (2006). Aquatic Habitat Modeling: Lecture Notes CE 526. University of Idaho.
College of Engineering-Boise. Ecohydraulics Research Group.
Jorde, K. & Bratrich, C. (1998). “River Bed Morphology and Flow Regulations in
Diverted Streams: Effects on Bottom Shear Stress Patterns and Hydraulic Habitat”. En:
Bretschko G. & Helešic J. (Eds.), Advances in River Bottom Ecology IV, Backhuys
Publishers, pp. 47-63.
Jorde, K.; Schneider, M.; Peter, A. and Zoellner, F. (2001). “Fuzzy Based Models for
the Evaluation of Fish Habitat Quality and Instream Flow Assessment”. Proceeding of
the 2001International Symposium on EnvironmentalHydraulics (ISEH).
223
King, J. & Brown, C. (2006). “Enviromental Flows: Striking the Balance between
Development and Resource Protection”. Ecology and Society Vol. 11, No. 2, Art. 26.
224
Meier, C. (2004). “Aspectos Hidrológicos e Hidráulicos de las Inundaciones por
Crecidas en Ríos: El Caso del BioBio en Hualqui”. En: Curso Intensivo Técnico en
Geomorfología Fluvial y Modelación de Hábitat Acuáticos. Universidad de Concepción.
Meza, R.I. & Vargas, H.X. (2009). “Fuzzy logic-based expert system for native fish
habitat assessment in a scarcity information context”. Proc. of JS.1 at the Joint IAHS &
IAH Convention, Hyderabad, India, September 2009. IAHS Publ. 328, 2009.
Ojeda, D. (2004). “Estudio del Ciclo de Vida y Parámetros Reproductivos del Puye
(Galaxias maculatus Jenyns 1842), en Condiciones de Cultivo Experimental”. Tesis
para el Grado de Licenciado en Ciencias de la Acuicultura. Temuco, Chile. Escuela de
Acuicultura, Universidad Católica de Temuco. 71 pp.
225
Olivares, M. & Lund, J. (2008). “Re-Regulation of Hourly Hydropower Operations”.
XXXI Congreso Interamericano AIDIS, Santiago, Chile.
Parra, O.; Valdovinos, C.; Habit, E. (2000). “Determinación del Caudal Mínimo
Ecológico del Proyecto Hidroeléctrico Quilleco. Unidad de Sistemas Acuáticos, Centro
EULA-Chile, Universidad de Concepción.
Parra, O.; Valdovinos, C.; Habit, E.; Figueroa, R. (2004). “Programa de Monitoreo de
la Calidad del Agua del Sistema Río Biobío”. Centro de Ciencias Ambientales EULA-
Chile, Universidad de Concepción.
226
Puppe, F. (1993). “Systematic Introduction to Expert Systems: Knowledge
Representations and Problem-Solving Methods”. EditorialSpringer-Verlag Berlin
Heidelberg, Cap. 1, pp. 3-6.
Redman, C.; Morgan, J.; Kuby, L. (2004). “Integrating Social Science into the Long-
Term Ecological Research (LTER) Network: Social Dimensions of Ecological Change
and Ecological Dimensions of Social Change”. Ecosystems, Vol. 7, pp. 161-171.
Roger, K.H. (2006). “The Real River Management Challenge: Integrating Scientists,
Stakeholders and Service Agencies”. River Research and Applications, Wiley
InterScience, Vol. 22, No. 2, pp. 269-280.
Saito, L.; Johnson, B.; Bartholow, J.; Hanna, B. (2001). “Assessing Ecosystem
Effects of Reservoir Operations Using Food Web-Energy Transfer and Water Quality
Models”. Ecosystems, Vol. 4, pp. 105-125.
Salazar, C.; Cuevas, T.; Cabrera, G.; Orphanópoulos, D.; Silva, A. (1994). Ecología
de Ríos. En: Terceras Jornadas de Hidráulica Profesor Francisco Javier Domínguez.
Sancha, A.M. (2004). “CI51I Calidad de Aguas”. Apuntes del curso Calidad del Agua.
Departamento de Ingeniería Civil, Universidad de Chile.
227
Schneider, M.; Jorde, K.; Kerle, F. (2004). “A Fuzzy Rule-Based Model for de
Assessment of Fish Habitats in Running Waters”. Submitted of Journal of Hydrology, in
review.
Schneider, M.; Jorde, K.; Zöllner, F.; Kerle, F. and Eisner, A. (2002). “Use of Habitat
Models for Decision Suport in Water Resources Management”. Proceeding of the 3rd
International Conferenceon Water Resources and Environment Research (ICWRER),
Dresden 2002.
Tello, E. (1999). “La Formación Histórica de los Paisajes Agrarios Mediterráneos: Una
Aproximación Coevolutiva”. Historia Agraria, No. 19, pp. 195-212.
Tickner, D; Armitage, P.; Bickerton, M.; Hall, K. (2000). “Assessing Stream Quality
Using Information on Messohabitat Distribution and Character”. Aquatic Conservation:
Marine and Freshwater Ecosystems, Vol. 10, pp. 179-196.
228
Universidad de Chile (1998). “Análisis de Criterios Hidroambientales en el Manejo de
Recursos Hídricos: Diseño de Plan de Monitoreo para la determinación de Caudales
Ecológicos”. División de Recursos y Medio Ambiente Hídrico, Departamento de
Ingeniería Civil, Universidad de Chile.
U.S. Army Corps of Engineers (USACE) (2002). “HEC-RAS River Analysis System:
Hydraulic Reference Manual, version 3.1”. Institute for Water Resources, Hydrologic
Engineering Center, California.
Vargas, J. & Mellado, M. (1994). Modelos de Calidad de Agua Fluvial. En: Terceras
Jornadas de Hidráulica Profesor Francisco Javier Domínguez.
Zadeh, L. (1965). “Fuzzy Sets”. Information and Control, Vol. 8, pp. 338-353.
229
REFERENCIAS WEB
230
Centro de Ecología Aplicada (CEA) (2004c). “Monitoreo de la Calidad de Agua y
Biota en el río Bíobío y Embalse Pangue” [en línea] <http://www.endesa.cl/Endesa_
Chile/monitoreo_biobio/Monitoreo_Mayo_2004.pdf > [consulta: 1 Diciembre 2008]
EULA (2007a). “Estudio de Impacto Ambiental Central Hidroeléctrica San Pedro”. [en
línea] <http://www.e-seia.cl/archivos/20071106.211445.pdf> [consulta: 13 enero 2009]
231
<http://www.eula.cl/evelyn-habit/pdf/Resultados%20San%20Pedro.pdf> [consulta: 13
enero 2008]
Harby, A., Baptist, M., Dunbar, M. J. and Schmutz, S. (editors) (2004). “State-of-the-
art in Data Sampling, Modelling Analysis and Applications of River Habitat Modelling”.
Cost Action 626, European Aquatic Modelling Network. ” [en línea]
<http://www.eamn.org/documents/COST%20626-State-of-the-art_new.pdf> [consulta:
19 enero 2009]
232
Kerle, F.; Zöllner, F.; Schneider, M.; Böhmer, J.; Kappus, B.; Baptist, M.J. (2002).
“Modelling of Long-term Fish Habitat Changes in Restored Secondary Floodplain
Channels of the River Rhine”. Conference Proceedings of the fourth Ecohydraulics
Symposium, 3-8 March 2002, Cape Town, South Africa. [en línea]
<http://www.sjeweb.de/pdf/Kerle_etal2002_1.pdf> [consulta: 05 septiembre 2008]
Manual for AQEM European Stream Assessment Program, version 2.3 (2004). [en
línea] <http://www.fliessgewaesser-bewertung.de/downloads/aqem_assessment_
software_manual_english_2_3.pdf> [consulta: 05 septiembre 2008]
Moraga, J. (2007). “El Agua: Derechos y Deberes”. Revista Capital, No. 200. [en línea]
<http://www.capital.cl/reportajes-y-entrevistas/el-agua-derechos-y-deberes.html>
[consulta: 28 noviembre 2008]
233
Ormazábal, C. (2004). “Caudal Ecológico, Metodologías y Casos Aplicados: El Caso
del Proyecto Ralco”. Seminario Ingendesa. [en línea] < http://www.aic.cl/pdf/
1015%20Caudal%20ecologico%20Ralco%20(Cesar%20Hormazabal).pdf> [consulta:
08 enero 2009].
Palau, A. (1994). “Los Mal Llamados Caudales ‘Ecológicos’. Bases para una Propuesta
de Cálculo”. Revista del Colegio de Ingenieros de Caminos, Canales y Puertos, nº 28.
[en línea] <http://hispagua.cedex.es/documentacion/revistas/op/28/op28_8.htm>
[consulta: 08 enero 2009]
Universidad de Chile (2005). “Informe País: Estado del Medio Ambiente en Chile”. [en
línea] <http://dspace.utalca.cl/retrieve/11060/informepais2006.pdf> [consulta: 30 Julio
2009]
U.S. Geological Service (USGS). Topografía digital escala 30 [m] de VIII región. [en
línea] <http://seamless.usgs.gov> [consulta: 11 octubre 2008]
234
Yañez, N. (2008). Taller Derechos de Agua Indígenas y de Usuarios Tradicionales los
Paradigmas para una Reforma de la Legislación de Aguas en Chile. [en línea]
<https://www.u-cursos.cl/derecho/2008/2/D121T0774/8/material_alumnos/objeto/
24242> [consulta: 28 noviembre 2008]
235