UNIVERSIDAD DE CANTABRIA ____________________________ ESCUELA TÉCNICA SUPERIOR DE INGENIEROS DE CAMINOS, CANALES Y PUERTOS DEPARTAMENTO DE CIENCIAS Y TÉCNICAS DEL AGUA Y DEL MEDIO AMBIENTE Curso 2008-2010 ____________________________ TESINA DE MÁSTER OFICIAL EN GESTIÓN AMBIENTAL DE SISTEMAS HÍDRICOS “Desarrollo de un modelo hidrológico como herramienta de apoyo para la gestión del agua. Aplicación a la cuenca del Río Laja, Chile” PRESENTADA POR: Enrique Muñoz Ortiz DIRIGIDA POR: César Álvarez Díaz Max Billib Santander, Febrero de 2010
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
UNIVERSIDAD DE CANTABRIA
____________________________
ESCUELA TÉCNICA SUPERIOR DE INGENIEROS DE CAMINOS, CANALES Y PUERTOS
DEPARTAMENTO DE CIENCIAS Y TÉCNICAS DEL AGUA Y DEL MEDIO AMBIENTE
Curso 2008-2010
____________________________
TESINA DE MÁSTER OFICIAL EN GESTIÓN AMBIENTAL DE SISTEMAS HÍDRICOS
“Desarrollo de un modelo hidrológico como herramienta de apoyo para la gestión del agua. Aplicación a la cuenca del
Río Laja, Chile”
PRESENTADA POR: Enrique Muñoz Ortiz
DIRIGIDA POR:
César Álvarez Díaz Max Billib
Santander, Febrero de 2010
UNIVERSIDAD DE CANTABRIA ESCUELA TÉCNICA SUPERIOR DE INGENIEROS DE CAMINOS CANALES Y PUERTOS
Profesores Patrocinantes
D. César Álvarez Dr. Max Billib Dr. José Luis Arumí Dr. Diego Rivera
“DESARROLLO DE UN MODELO HIDROLÓGICO COMO HERRAMIENTA DE APOYO PARA LA GESTIÓN DEL AGUA. APLICACIÓN A LA CUENCA DEL RÍO LAJA, CHILE”
ENRIQUE ALEJANDRO MUÑOZ ORTÍZ
Santander, Febrero de 2010
ÍNDICE
CAPÍTULO 1: INTRODUCCIÓN 1
1.1 Introducción 1
1.2 Objetivo General 2
1.3 Objetivos Específicos 2
CAPÍTULO 2: FUNDAMENTOS Y SELECCIÓN DE UN MODELO HIDROLÓGICO 4
2.1 Fundamentos de los Modelos Hidrológicos 4
2.2 Procesos predominantes a diferentes escalas espaciales 7
2.3 Selección de un Modelo Hidrológico Mensual 8
CAPÍTULO 3: DESCRIPCIÓN Y ADAPTACIÓN DEL MODELO HIDROLÓGICO
MENSUAL (MHM) 12
3.1 Descripción del MHM propuesto por Ferrer et al. (1973) 12
3.2 Incorporación de la Componente Nival al MHM 16
3.3 Incorporación Módulo de Extracciones y Aportes 19
CAPÍTULO 4: APLICACIÓN A LA CUENCA DEL RÍO LAJA 21
4.1 Modelo Físico 21
4.2 Modelo Meteorológico 24
4.3 Registro de Caudales 30
4.4 Calibración y Validación 31
4.4.1 Módulo Pluvial 31
4.4.2 Módulo Nival 34
4.4.3 Parámetros de Calibración y Validación 36
4.5 Simulación (1990-2008) 37
CAPÍTULO 5: DISCUSIÓN Y ANÁLISIS DE RESULTADOS 44
5.1 Ajuste del Modelo 44
5.2 Representatividad de los Caudales Simulados y Utilidad del Modelo para Fines
Predictivos y de Gestión de Recursos 46
CAPÍTULO 6: CONCLUSIONES Y FUTURAS LÍNEAS DE INVESTIGACIÓN 48
6.1 Conclusiones 48
6.2 Futuras Líneas de Investigación 48
BIBLIOGRAFÍA 50
ANEXO A: Script Modelo Hidrológico Mensual.
ANEXO B: Base de Datos de Precipitación, Temperatura y Evapotranspiración Potencial
(Calculada) y Series Ponderadas para cada Subcuenca.
ANEXO C: Extracto Guía Explotación Lago Laja (DOH).
ANEXO D: Base de Datos de Caudales Generación Bruta y Nivel del Lago Laja y de Registro
de Caudales de la Cuenca del Laja.
ANEXO E: Gradiente Térmico y Curvas Hipsométricas de Subcuencas de régimen nivo-pluvial.
ANEXO F: Caudales simulados serie 1990 – 2008.
Capítulo 1: Introducción 1
CAPÍTULO 1: INTRODUCCIÓN
1.1 Introducción El agua dulce es un recurso natural limitado pero renovable, en el planeta corresponde sólo al 3%
del total de agua (Figura 1.1), siendo este un elemento imprescindible para la sobrevivencia y el
desarrollo de seres humanos y ecosistemas. Del total de agua dulce, sólo un 1% se encuentra en
superficie y es de fácil acceso, el resto se encuentra congelada (79%) en glaciares y casquetes de
hielo polar, y subterránea (20%) (Hinrichsen et al., 1998).
Figura 1.1 Distribución de Agua en el Planeta (Sther et al., 2008)
Considerando que el agua es un recurso limitado, que la demanda por este aumenta con el
desarrollo y crecimiento poblacional, que la disponibilidad cambia año a año, a causa de
fenómenos locales, globales, naturales y antropogénicos, que es de importancia fundamental el
impacto sobre el abastecimiento causado por el cambio climático (Arnell et al., 1995, Kaczmarek
et al., 1995) y que gran parte de la población mundial ya experimenta un estrés por la demanda
de agua (Vörösmarty et al., 2000), se hace necesario desarrollar herramientas que permitan una
gestión eficiente del recurso y que sirvan de apoyo para predecir condiciones futuras bajo
diferentes escenarios causados por efectos locales (orográficos), regionales ( “El Niño” o “La
Niña”) o globales (cambio climático y/o tendencias graduales de largo plazo), con el fin de evitar
o reducir el estrés por la demanda de agua sobre diferentes cuencas en el planeta.
Capítulo 1: Introducción 2
La cuenca del Río Laja, Chile (Figura 1.2) es un ejemplo de dicho estrés. Históricamente se ha
visto afectada por alteraciones antropogénicas que buscan el uso de sus caudales con fines
socioeconómicos hasta niveles de sobreoferta, un ejemplo de esto es una de las últimas obras
ejecutadas sobre la cuenca, el canal “Laja-Diguillín”, canal que está operativo desde Abril de
2003, que tiene derechos de agua asignados por 65 [m3 s-1] y que a la fecha (Febrero de 2010) no
ha podido explotar más de 20 [m3s-1]
Figura 1.2: Ubicación cuenca del Río Laja
En el presente trabajo se busca desarrollar un modelo hidrológico adaptable a diferentes
escenarios, tanto físicos como meteorológicos e hidrológicos y que sirva como herramienta para
apoyar desde un punto de vista técnico la administración y gestión de agua de una cuenca,
considerando condiciones pasadas y posibilidades de cambio ante diferentes escenarios futuros.
Lo anterior se resume en los siguientes objetivos
1.2 Objetivo General Desarrollar una herramienta de apoyo para la gestión y administración del agua de una
cuenca hidrográfica con alteraciones hidromorfológicas.
1.3 Objetivos Específicos Desarrollar un modelo hidrológico a escala mensual capaz de simular regímenes de tipo
pluvial, nival y nivo-pluvial, y de incluir alteraciones antropogénicas sobre la cuenca, o
directamente sobre el río (por ejemplo la operación de centrales hidroeléctricas, canales
de trasvase o regadío, deforestación de cuencas, modificación del tipo de vegetación,
entre otros).
Capítulo 1: Introducción 3
Aplicar, calibrar y validar el modelo al caso de la cuenca del Río Laja (Chile) y
representar caudales históricos producidos por la cuenca.
Definir el nivel de confiabilidad del modelo y de los caudales simulados para
aplicaciones predictivas y de gestión.
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 4
CAPÍTULO 2: FUNDAMENTOS Y SELECCIÓN DE UN MODELO HIDROLÓGICO
2.1 Fundamentos de los Modelos Hidrológicos Un modelo intenta reproducir bajo diferentes formas y procesos un fenómeno físico que ocurre
sobre un objeto o territorio. Por lo tanto, en hidrología, un modelo busca representar un territorio
delimitado por una divisoria de aguas (cuenca), y los fenómenos de transferencia de lluvia a
caudal y de agua en el interior de esta. El objeto de reproducir dichos procesos es el poder
simular y predecir condiciones futuras con el fin de actuar desde un punto de vista
socioeconómico y de gestión, administración y optimización de los usos del agua.
Una de las principales ventajas de la simulación mediante modelos es la visión que se gana
recogiendo y organizando los datos requeridos como entrada a los algoritmos matemáticos que
comprenden el sistema completo. Lo anterior puede conducir a la readecuación o ampliación de
las redes de monitoreo existentes de tal manera de mejorar los datos de entrada necesarios para la
modelación y con esto lograr una mejor representación del sistema natural (Sther et al., 2008). El
crecimiento en el uso de los modelos se debe a la necesidad de solucionar problemas complejos,
la mejora en las técnicas de medición de los parámetros de entrada y a la actual disponibilidad de
tecnología. Sin embargo no se puede perder de vista que los modelos solo son buenos en la
medida que los datos de entrada lo son (Rasmussen et al., 1998).
Los modelos pueden ser diferenciados por su arquitectura y estructura, la escala a la que son
típicamente aplicados, la integración o falta de integración de los procesos de superficie y
subsuperficie, el grado de complejidad de los datos de entrada y el numero de parámetros a ser
determinado (Mishra & Singh 2004). En la actualidad existen diferentes tipos de modelos, los
que se pueden clasificar en dos grandes grupos, determinísticos y estocásticos. Los modelos
determinísticos son aquellos que las variables existentes en el algoritmo matemático están
determinadas por leyes físicas, ya sea empíricas, conceptuales o teóricas y por lo tanto su
variabilidad está físicamente justificada. En cambio los modelos probabilísticos rigen una parte o
la totalidad de sus variables por leyes del azar (Estrela, 1992). Estos modelos se pueden
diferenciar a su vez según la representación espacial del modelo (agregados, semidistribuidos o
distribuidos), según la forma de representar los procesos hidrológicos (métricos, conceptuales o
físicamente basados) y según la extensión temporal en la que se puede aplicar el modelo
(modelos de episodios o continuos) (Ponce, 1989, Beck, 1991, Wheater et al., 1993, Sempere,
1996, Francés, 1996, Kokkonen y Jakeman, 2001).
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 5
Los diferentes tipos de modelos se describen a continuación.
Según su distribución espacial.
Modelos agregados: Modelos con distribución espacial uniforme en la cuenca, utiliza
como variable la precipitación media y asume que los parámetros de los diferentes
submodelos de los procesos hidrológicos son uniformes para toda la cuenca y
permanecen constantes a lo largo de un paso de tiempo o de una simulación.
Modelos semidistribuidos: Son aquellos que permiten una cierta variabilidad espacial
de la lluvia y de los parámetros de los submodelos que lo componen mediante la división
de la cuenca en multitud de pequeñas subcuencas con lluvia y parámetros constantes en
cada uno de ellos.
Modelos distribuidos: Son aquellos que permiten la variabilidad espacial de la lluvia y
de los parámetros, mediante la división de la cuenca en celdas o pixeles, en las que se
simulan los diferentes procesos hidrológicos.
Según su forma de representar los procesos hidrológicos.
Modelos métricos: Son aquellos con gran dependencia respecto de los datos observados.
Caracterizan la respuesta del sistema mediante un método de extracción de la
información a partir de datos existentes. Estos modelos presentan una consideración baja
o nula de los procesos físicos que ocurren en el sistema hidrológico.
Los modelos métricos utilizan la representación más simple del comportamiento de una
cuenca hidrológica, ya que la relación entre lluvia y caudal es simulada mediante una
variable de pérdida volumétrica, que engloba los procesos de pérdida por evaporación,
almacenamiento de humedad en el suelo y recarga de aguas subterráneas, y una función
distribución temporal, que simula los diferentes modos dinámicos de respuesta de la
cuenca.
El ejemplo más común son los modelos basados en el concepto de hidrograma unitario,
que simula la respuesta de caudal a un episodio de lluvia mediante una función de
pérdidas no lineal y una función de transferencia lineal (Sherman, 1932).
Una ventaja es que requieren una cantidad mínima de datos, pero como desventaja se
tiene que el rango de aplicación es limitado por la variabilidad de los datos observados,
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 6
los cuales no son capaces de tener en cuenta los cambios producidos en la cuenca (por
ejemplo deforestación, forestación, construcción de represas o embalses, etc.).
Modelos conceptuales: Son modelos que reproducen los procesos hidrológicos más
predominantes mediante una base de conocimiento inicial en forma de representación
conceptual de los mismos (estos modelos comienzan a ser utilizados con el desarrollo y
aumento de la capacidad computacional). La representación hidrológica se realiza
mediante relaciones simplificadas con parámetros que no se pueden medir físicamente en
la realidad. Su aplicación está condicionada a la calibración con datos observados en la
cuenca.
Modelos físicamente basados: Estos modelos incluyen el comportamiento físico de los
procesos hidrológicos, realizando la simulación de una cuenca mediante el uso de las
ecuaciones de continuidad clásicas, resolviendo ecuaciones diferenciales de forma
numérica mediante la aplicación de métodos de diferencias finitas o elementos finitos.
Estos modelos son necesariamente de tipo distribuido, y por lo tanto, permiten la
descripción de los fenómenos en cada una de las celdas del sistema, tanto desde el punto
de vista de los parámetros como de la resolución de los sistemas de ecuaciones
diferenciales, por tanto el sistema físico se representa mediante un sistema de ecuaciones
diferenciales que expresan la masa, el momento y el balance de energías (Mediero,
2007).
La principal ventaja de un modelo físicamente basado es que utiliza parámetros que
pueden ser medidos en terreno y que tienen un significado físico directo, por tanto, si los
valores de estos parámetros se pueden determinar a priori, estos modelos pueden ser
aplicados a cuencas sin datos observados, e incluso se pueden tener en cuenta los cambios
en la misma cuenca. Por otra parte, los modelos físicamente basados permiten una
representación matemáticamente ideal de un fenómeno real, pero la desventaja es que
requieren una cantidad muy grande de datos que son muy difíciles de medir en la realidad
y que tienen un alto costo asociado. Además requieren de una capacidad y tiempo
computacional elevado.
Según la extensión temporal en la cual se puede aplicar el modelo.
Modelos de episodio: Son modelos desarrollados para simulaciones de cortos intervalos
de tiempo, normalmente para un único episodio de lluvia. Estos modelos se centran en la
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 7
simulación de los procesos de infiltración y escorrentía superficial, ya que su principal
objetivo es la evaluación de la escorrentía directa para cálculo de avenidas. Su aplicación
se enfoca principalmente en casos donde la escorrentía directa es la principal fuente de
generación de la escorrentía total. Una desventaja es su limitación para la simulación de
caudales de estiaje y diarios.
Modelos continuos: Permiten la simulación de caudales diarios, mensuales y
estacionales, es decir, permiten la simulación entre episodios de lluvia. Estos modelos
consideran todos los procesos que influyen en la generación de escorrentía, tanto la
escorrentía directa mediante el flujo de superficie y subterráneo. Los modelos continuos
se centran en la simulación de la evapotranspiración y de los procesos que influyen a
largo plazo en la recuperación de la humedad durante los periodos en los que no hay
ocurrencia de precipitación. Por tanto, el objetivo principal de los modelos continuos es
la simulación del balance de humedad total de la cuenca durante períodos largos de
tiempo. Los modelos continuos se aplican fundamentalmente para el pronóstico de
volúmenes de escorrentía a largo plazo y estimación de las reservas de agua.
Para la optimización, protección y gestión de los recursos de agua en una cuenca requieren de
una modelación satisfactoria y por lo tanto de la selección del modelo adecuado en función del
objetivo de la simulación.
2.2 Procesos predominantes a diferentes escalas espaciales La respuesta hidrológica cambia según las distintas escalas espaciales (Bloschl y Sivapalan 1995,
Bloschl 1996).En términos generales se tiene que a nivel de micro cuenca (1 – 10 km2), la
respuesta a las precipitaciones es dominada por procesos de generación de escorrentía en laderas
y áreas cercanas al río (Anderson y Burt, 1990; Montgomery et al., 1997), por lo tanto las
propiedades del suelo y uso de suelo predominan en la respuesta hidrológica, y la distribución de
precipitaciones puede ser representada en términos medios y uniformes. A escala de meso-
cuencas (10 – 103 km2), la variación espacial y temporal toma mayor importancia y esta
dependerá del tipo de precipitación que normalmente se produce sobre la cuenca (frontal,
convectiva o por orografía). Los procesos hidrológicos predominantes son la escorrentía directa
y subterránea. A escala de macro cuencas (103 - 104 km2), la distribución espacial y temporal de
las precipitaciones y el rastreo del caudal domina considerablemente el comportamiento de la
respuesta (Uhlenbrook et al., 2004) y en los procesos hidrológicos influyen variables como la
evapotranspiración y las características del suelo y nivel freático en cuanto a capacidad de campo
y permeabilidad y profundidad o características de los estratos de suelo.
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 8
2.3 Selección de un Modelo Hidrológico Mensual La selección del modelo a utilizar depende directamente del objetivo del estudio. En este caso se
requiere de una herramienta que sirva como apoyo técnico para la gestión y distribución de
caudales medios de una cuenca y para predecir caudales a corto y mediano plazo. Por lo tanto es
necesario un modelo continuo a escala mensual (capaz de simular condiciones pasadas y
futuras). Por otra parte, considerando la poca disponibilidad de información distribuida en cuanto
a tipos de suelo, vegetación, estratos, evaporación y temperatura, se hace más conveniente y
práctico el uso de un modelo agregado, el cual en su aplicación puede ser discretizado en
subcuencas emulando un modelo semidistribuido. Finalmente, según las condicionantes
anteriores, se tiene que el modelo debe ser de tipo conceptual, ya que un modelo físicamente
basado requiere de información distribuida y un modelo métrico no permite considerar cambios
en la cuenca y por lo tanto es muy limitado en cuanto a los escenarios de simulación posibles.
Entre los modelos de tipo conceptual, agregado y continuo disponibles se encuentran los
siguientes:
“T-M Model”: Este modelo fue desarrollado inicialmente por Thornthwaite y Mather
(1955) como un modelo de dos parámetros, la capacidad de almacenamiento en términos
de humedad en el suelo y la fracción de agua excedente o infiltrada en el suelo.
“Palmer Model”: Con la elaboración de un índice de sequía, Palmer (1965) modificó el
“T-M Model” al separar el suelo en dos capas con diferentes capacidades de humedad y
por lo tanto de evaporación e infiltración condicionando la respuesta de la capa inferior a
la evolución de la capa superior. Con esto pasa a ser modelo de tres parámetros.
“abcd Water Balance Model”: Thomas (1981), en base al “T-M Model” y “Palmer
Model”, propuso un modelo con cuatro parámetros donde el suelo se divide en tres zonas,
una primera zona relacionada con el almacenamiento superficial, cuyo caudal es
despreciable respecto de la escorrentía directa y subterránea, una zona media de suelo no
saturado que asimila la recarga que experimenta el acuífero y una tercera zona de suelo
saturado (nivel freático). Los parámetros del modelo representan el nivel de saturación
del suelo (expresado en forma de fracción respecto de la humedad máxima), la humedad
máxima del suelo (y por lo tanto la máxima evaporación real), la fracción del
almacenamiento subterráneo que se transforma en escorrentía y el recíproco del tiempo
de residencia del agua subterránea.
Capítulo 2: Fundamentos y Selección de un Modelo Hidrológico Mensual 9
“Monthly Water Balance Model”: McCabe y Markstrom (2007) realizaron una
adaptación del “T-M Model” donde se incluye la componente nival y se modifican
determinados parámetros, pero sin alterar la estructura ni la base conceptual del modelo
original. Así este modelo pasa a tener siete parámetros (coeficiente de escorrentía,
coeficiente de escorrentía directa, la capacidad de almacenamiento en términos de
humedad en el suelo, la ubicación geográfica de la cuenca, una temperatura nival que
diferencia la precipitación nival de la pluvial, una temperatura pluvial que permite
determinar una proporción entre la precipitación nival y pluvial, y la máxima tasa de
Capítulo 3: Descripción y Adaptación del Modelo Hidrológico Mensual (MHM) 20
Donde Qsalida es el caudal de salida de la cuenca, Qproducido es el caudal producido en conjunto por
la cuenca pluvial y nival, Qaportes es el caudal proveniente de cuencas ubicadas aguas arriba o de
canales de trasvase y Qextracciones, son extracciones desde la cuenca para un uso agrícola (riego) o
consumo humano o industrial.
Con las modificaciones indicadas en los puntos 3.2 y 3.3 se convierte el modelo original de
carácter agregado y aplicable sólo a régimen pluvial, en un modelo semidistribuido aplicable a
un régimen nivo-pluvial (el script del modelo se adjunta en el Anexo A).
Capítulo 4: Aplicación a la Cuenca del Río Laja 21
CAPÍTULO 4: APLICACIÓN A LA CUENCA DEL RÍO LAJA
El Río Laja es un río de régimen nivo-pluvial de uso múltiple donde se produce una compleja
interacción entre los componentes naturales, sistemas económicos y sociales que controlan el uso
y manejo de los recursos (Mardones y Vargas, 2005). Sobre este tipo de cuencas, de carácter
montañoso y de abundante disponibilidad de recursos hídricos y forestales existe un especial
interés ecológico y económico, siendo fundamental en la ordenación territorial y en los planes
energéticos nacionales (García-Ruiz, 1990, Sala y Batalla, 1999).
Con el fin de validar el MHM descrito y de comprobar su capacidad de adaptación a una cuenca
con alteraciones hidrológicas, este se aplica sobre la cuenca del Río Laja, la cual históricamente
ha sido alterada en su régimen por hidroeléctricas, trasvases y extracciones por medio de canales
de regadío.
4.1 Modelo Físico La cuenca del Río Laja se ubica entre los 36° 52’ S, 72° 38’ O y 37° 39’ S, 71° 12’O. Tiene un
área total de 4668 [km2] y se emplaza entre la cordillera de los andes por el Este, entre los
complejos volcánicos de Chillán y Antuco-Sierra Velluda por el Norte y Sur respectivamente y
por la depresión intermedia por el Oeste en su desembocadura en el Río Bíobio. Su altitud va
desde los 50 y 3534 [msnm]. El sector andino del río (600 a 2000 [msnm]) está estructurado por
rocas sedimentarias y plutónicas sobre las cuales se encuentran rocas volcánicas y depósitos
volcano-clásticos, procedentes del complejo volcánico Antuco-Sierra Velluda (Moreno y Varela,
1987, Niemeyer y Muñoz, 1983). En la depresión central de la cuenca se encuentra un gran
depósito aluvial llamado “Cono de Arenas Negras del Laja” (entre 100 a 300 [msnm]), el cual se
encuentra sobre una formación lahárica. Estos depósitos se terminan al entrar en contacto con
cordillera de la costa, la cual se estructura en granito (este límite representa todo el borde
noroccidental de la cuenca del Río Claro).
La combinación de arenas y el depósito lahar consiste en un suelo permeable de base
impermeable que favorece la creación y mantenimiento de reservas aguas subterráneas en la
depresión central.
La Figura 4.1 muestra el modelo físico a utilizar en las simulaciones, y se describen las
interconexiones entre los diferentes elementos y subcuencas.
Capítulo 4: Aplicación a la Cuenca del Río Laja 22
Polcura a/des Central El Toro (CS)Laja en
Tucapel (CS)Laja en P. Perales (CS)
Claro en Yumbel (CS)
Laja en San Rosendo (CS) Lago Laja
Río
Pol
cura
Río Laja
Río Laja
Río
Cla
ro
Cor
dille
ra d
e Lo
s A
ndes
Complejo Antuco-Sierra Velluda
Complejo Chillán
Depresión Intermedia
Sector Andino
Cor
dille
ra d
e la
Cos
ta
Desembocadura en Río Bíobio
SC-1(A=236.5 km2)
SC-2(A=977.3 km2)
SC-3(A=635.4 km2)
SC-4(A=947.4 km2)
SC-5(A=807.7 km2)
SC-6(A=847.8 km2)
SC-7(A=215.6 km2)
N
Polcura en Cuatro Juntas (CS)
Can
al A
lto P
olcu
ra
Desc. C. El Toro
Filtraciones del Laja
Extr
acci
ones
Extr
acci
ones
Figura 4.1: Esquema modelo físico cuenca del Laja
En la figura anterior se presenta esquema en donde se indican las siete subcuencas del modelo
(SC-1 a SC-7), y el área de cada una. En azul se presentan los principales ríos y en rojo los
canales de trasvase o extracción. En cuanto al Lago Laja, cabe mencionar que es un lago de
carácter natural, regulado artificialmente por el uso hidroeléctrico, cuyo caudal medio es del
orden de 34 [m3s-1] y llega al río a través de la descarga de la “Central El Toro”. Por otra parte,
desde el lago se producen filtraciones, las cuales se concentran en “Los Ojos del Laja”, desde
donde nace el río con un caudal medio en torno a los 27 [m3s-1]. A continuación se describen las
características y alteraciones sobre cada subcuenca.
SC-1, “Polcura Alto”: Es una cuenca nivo-pluvial donde el caudal producido por la
subcuenca es dividido en dos salidas, una que va a la SC-3 y que es controlada por la
estación “Polcura en Cuatro Juntas” y otra que va al Lago Laja (SC-2) transportado por el
“Canal Alto Polcura”. La cuenca se ubica entre las cotas 1350 y 2910 [msnm] y se
encuentra sobre la cadena montañosa de la cordillera de los andes, donde predomina un
suelo volcánico de baja permeabilidad.
SC-2, “Lago Laja”: El Lago Laja se ubican en una ladera de la cordillera de los andes,
sobre una depresión de suelo volcánico de baja permeabilidad, este genera dos caudales
de salida que descargan sobre la SC-4. Una salida es producida por las filtraciones del
Laja, cuyo caudal de salida depende del nivel del lago. La otra salida es producida por la
hidroeléctrica “El Toro”, la cual descarga aguas abajo de la estación “Polcura Antes
Descarga Central El Toro”. En este caso los caudales producidos varían según la
Capítulo 4: Aplicación a la Cuenca del Río Laja 23
disponibilidad de agua en el lago y la demanda de electricidad, por otra parte, el cálculo
del caudal se realiza en base a la estadística de energía producida.
SC-3, “Polcura”: Esta cuenca recibe una entrada proveniente de la SC-1, descarga sobre
la SC-4 y es controlada por la estación “Polcura Antes Descarga Central El Toro”. Se
ubica sobre la zona montañosa del complejo volcánico Chillán, entre las cotas 800 y 3089
[msnm], cuyo suelo es poco permeable y de pendientes fuertes ( %6.1≥ ).
SC-4, “Laja Alto”: Esta subcuenca es la más compleja dentro del sistema Laja, debido a
que recibe tres entradas (descarga central El Toro, infiltraciones del Laja y la descarga de
la SC-3), y por otra parte, tiene ocho extracciones (canales Zañartu, Collao, Mirrihue, El
Litre, Bulnes, Ortiz, Laja-Diguillín y Laja Sur). Cinco de estos canales presentan
registros de caudal, y el caudal extraído por los tres restantes se puede estimar como el
22% de las extracciones realizadas por los canales Zañartu, Collao y Mirrihue (DOH,
2007). En cuanto al régimen, es una cuenca nivo-pluvial, ubicada en el sector andino del
río, entre las cotas 289 y 3535 [msnm], presenta laderas de pendientes medias y roca
volcánica, y con una planicie o depresión central de suelo granular con una subcapa
impermeable.
Esta subcuenca se encuentra controlada por la estación “Laja en Tucapel”, la cual además
separa la zona alta y baja del Río Laja y significa un punto de concentración de los
caudales superficiales y subsuperficiales producidos por la cuenca.
SC-5, “Laja Medio”: Esta subcuenca puede ser modelada como una cuenca de régimen
pluvial que recibe los caudales provenientes de la SC-4 y descarga sobre la SC-7. En su
salida se encuentra la estación “Laja en Puente Perales” que cumple la función de control
y permite la calibración del modelo con una serie continua y confiable para el período de
estudio (1990-2008). Además representa el 77% del área total de la cuenca del Laja,
siendo fundamental en la calibración y validación del estudio. La cuenca se ubica en un
depresión que recibe los depósitos aluviales provenientes de aguas arriba llamado “Cono
de Arenas Negras del Laja”, depósitos de una alta permeabilidad, limitados por debajo
por una capa lahárica impermeable que favorece la presencia de un embalse y flujo
subterráneo.
SC-6, “Río Claro”: Esta es una subcuenca independiente, sin entradas de caudal desde
aguas arriba y de régimen pluvial. El caudal producido descarga sobre la SC-7. En cuanto
Capítulo 4: Aplicación a la Cuenca del Río Laja 24
a la geología y geomorfología, es una cuenca muy particular, ya que el área ubicada al
oeste del Río Claro es una zona montañosa, forestada, de suelo de granito de baja
permeabilidad y de pendientes fuertes ( %1≥ ), en cambio la zona ubicada al este del río
es una zona de bajas pendientes ( %8.03.0 − ), con baja densidad de bosques, cuyo suelo
es mayoritariamente usado para agricultura y de mayor permeabilidad que el suelo del
lado oeste.
La calibración de la SC-6 se lleva a cabo con la estación “Río Claro Camino a Yumbel”,
la cual controla el 81% del área de la cuenca del Río Claro. Luego, para calcular el caudal
total producido por la SC-6, se realiza una simulación con la cuenca completa y con los
parámetros obtenidos de la calibración.
SC-7, “Río Laja en Desembocadura”: Esta subcuenca recibe dos entradas, una
proveniente de la descarga del Río Claro y una proveniente de la subcuenca Laja Medio.
Es una cuenca sin control fluviométrico, por lo tanto, se utilizan los mismos parámetros
de calibración de la SC-5, ya que las características geológicas, geomorfológicas e
hidrológicas son prácticamente equivalentes. La descarga de esta subcuenca corresponde
al caudal de salida de toda la cuenca del Río Laja en la desembocadura en el Río Bíobio.
4.2 Modelo Meteorológico Para el desarrollo de la simulación es necesario determinar la precipitación (PM), evaporación
potencial (EM) y la temperatura media mensual (Tm, necesaria para una cuenca con régimen
nival) representativa de cada subcuenca del modelo físico. En la cuenca del Laja y sus cercanías
existen 12 estaciones pluviométricas controladas por la Dirección General de Aguas de Chile
(DGA), con las cuales se puede determinar una precipitación representativa de cada subcuenca,
pero no existe información continua y confiable de EM y Tm. Por otra parte, el Departamento de
Geografía de la Universidad de Delaware ha publicado una base de datos de PM y Tm para el
período 1900 – 20081, la cual se ha obtenido en base a una interpolación espacial de 0.5 x 0.5
grados [latitud/longitud]. Esta interpolación llamada “Climatologically aided interpolation
(CAI)” (Willmott y Robeson, 1995) fue utilizada para determinar los valores mensuales de PM y
Tm del mundo, en base a la siguiente red de estaciones.
1 http://climate.geog.udel.edu/~climate/html_pages/archive.html Universidad de Delaware, Departamento de Geografía. Base de datos meteorológica consultada en Enero de 2010.
Finalmente se hace una corrección del valor calculado, por el número de días de cada
mes (N) y horas de sol que hay en promedio por día en el mes (d).
3012dNEMEM nc ⋅⋅= [4.5]
Donde el número máximo de horas de sol al día es una función de la latitud y la
época del año, la cual se extrae de la Tabla 4.2 (Doreenbos y Pruitt, 1977).
Tabla 4.3: Cantidad de horas de sol al día según latitud Lat. Norte Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov DicLat. Sur Jul Ago Sep Oct Nov Dic Ene Feb Mar Abr May Jun
La tabla anterior opera de la siguiente manera. Por ejemplo, para calcular la PM (DGA) de la
SC-7 de Enero de 2005, se debe multiplicar la PM(Ene, 2005) de las estaciones pluviométricas
Laja y Achiras Fundo por el área representativa de cada estación, en este caso 183.7 y 31.9 [km2]
respectivamente, luego, la suma de ambas se debe dividir por el área total de la cuenca (215.6
[km2]), con lo que se obtiene el valor de la PM (DGA) ponderado para la SC-7.
3 Nota: Si bien, el período de simulación seleccionado en el punto 4.5 es 1990-2008, en el Anexo B se adjunta la base de datos del período 1960-2008 ya que dicha información se espera se utilizada en futuras investigaciones según se indica en el Capítulo 5.
Capítulo 4: Aplicación a la Cuenca del Río Laja 30
4.3 Registro de Caudales En los puntos anteriores se han definido las entradas necesarias para realizar una simulación,
pero, considerando la existencia de conexiones entre cuencas, las extracciones realizadas por
canales de regadío y trasvase y los aportes desde el Lago Laja, es necesario definir aquellos
caudales medidos y necesarios para una simulación que incluya las alteraciones mencionadas.
Por otra parte, para la calibración y validación, y para determinar el nivel de ajuste de los
caudales simulados respecto de los medidos, se presentan los registros de las estaciones
fluviométricas de control.
Tabla 4.5: Resumen de las series de registros de caudales
Estación
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
Canal Alto PolcuraPolcura en Cuatro JuntasPolcura a/d Central El ToroLaja en Tucapel *Laja en Puente Perales * *Claro camino a YumbelLaja en San Rosendo * * * *Caudal Central El Toro(1)
Infiltraciones Lago Laja(2)
Canal ZañartuCanal CollaoCanal MirrihueCanal El Litre, Bulnes y Ortiz(3)
Canal Laja DiguillinCanal Laja Sur
Estación
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
Canal Alto Polcura *Polcura en Cuatro Juntas * *Polcura a/d Central El Toro * * * *Laja en Tucapel * * *Laja en Puente PeralesClaro camino a Yumbel * * *Laja en San RosendoCaudal Central El Toro(1)
Infiltraciones Lago Laja(2)
Canal ZañartuCanal Collao * *Canal Mirrihue * *Canal El Litre, Bulnes y Ortiz(3) * *Canal Laja Diguillin * *Canal Laja Sur * *
Sin InformaciónRegistros de caudal medio mensualRegistros de caudal medio mensual (Estación de Control)
* Año incompleto(1) Caudal calculado a partir de registros de potencia (ver en Anexo C y D)(2) Caudal calculado a partir de registros del nivel del Lago Laja (ver en Anexo C y D)(3) Caudal estimado según manual de operación del Lago Laja (ver en Anexo C y D)
En el Anexo C se presenta un extracto del documento “Guía de Explotación del Sistema Lago
Laja” (DOH, 2007), donde se indica la metodología de cálculo de los caudales descargados por
la central “El Toro”, la metodología de cálculo de las “Infiltraciones del Laja” y la metodología
de estimación del caudal pasante por los canales El Lite, Bulnes y Ortiz. Posteriormente, en el
Capítulo 4: Aplicación a la Cuenca del Río Laja 31
Anexo D se presentan los registros de energía producida por la central “El Toro” (Tabla D.1), los
niveles del Lago Laja (Tabla D.2) y los caudales descritos en la Tabla 4.5 (Tabla D.3 a D.18).
En la tabla anterior se observa que el período común de registros en la mayoría de las estaciones
es entre 2003 y 2007. En este período se tiene información completa entre Agosto de 2003 y
Junio de 2006, por lo que se decide realizar la calibración y validación en dicho período.
La estación “Claro camino a Yumbel” no posee registros en el período mencionado, pero
considerando que esta estación controla la SC-6, la cual es completamente independiente del
resto del sistema Laja y que desemboca directamente en la SC-7, esta se puede calibrar y validar
para el período con información disponible (Enero 1986 – Diciembre 1989 calibración y Enero
1991 – Diciembre 1995 validación).
4.4 Calibración y Validación La calibración se lleva a cabo a través de un proceso iterativo, en donde se definen unos
parámetros de partida considerando características del suelo, vegetación, morfología, geología,
ubicación y tipo de régimen de la cuenca, y luego estos se van modificando según el nivel de
ajuste de los caudales simulados respecto de los medidos.
La metodología y caracterización de los parámetros del proceso de calibración se describe a
continuación.
4.4.1 Módulo Pluvial El módulo pluvial está compuesto por dos series que constituyen las variables de entrada y el
área de la cuenca (PM, EM y AR), tres valores de partida (humedad inicial del suelo, escorrentía
subterránea inicial y percolación profunda inicial), seis parámetros de calibración (Cmax, Hmax, D,
Plim, PORC y Ck) y dos parámetros de modificación de los inputs (A y B). A continuación se
describen y caracterizan los componentes presentes en este módulo.
Variables de Entrada: La precipitación y evaporación potencial se deben determinar
en base a mediciones y deben ser representativas de la cuenca. Definen la longitud
temporal de la simulación y de los resultados. Estos valores se determinan en el
modelo meteorológico (series representativas adjuntas en el Anexo B). Por otra parte,
el modelo requiere como entrada el área de la cuenca, la cual se obtiene a partir de
una base física y permite al término de la simulación transformar las escorrentías
estimadas [mm mes-1] en caudales de salida [m3 s-1].
Capítulo 4: Aplicación a la Cuenca del Río Laja 32
Valores de Partida: Debido a que el modelo realiza los cálculos un paso de tiempo
hacia delante, requiere de valores de partida para la simulación. Estos valores se
pueden estimar inicialmente o definir como nulos, y luego de una simulación se
pueden modificar en función de lo observado en el mismo mes del año siguiente y
subsiguiente. Se debe tener en consideración la variabilidad interanual de la humedad
en el suelo, de la escorrentía subterránea y de percolación profunda, para esto es
recomendable observar los gráficos de variación de dichas variables en el tiempo.
Además es importante considerar la variabilidad del flujo base, el tipo de suelo en
relación a su permeabilidad, el espesor estimado del estrato y condiciones geológicas
que favorezcan la retención de humedad o la mantención de un flujo subterráneo.
Parámetros de Calibración: Estos parámetros emulan condiciones físicas
distribuidas de modo uniforme en la subcuenca. A continuación se describen los
diferentes parámetros del modelo, sus rangos de variación, y los valores asociados a
diferentes características del suelo, subsuelo y de la cuenca, los cuales se estiman en
función de la los parámetros obtenidos en el desarrollo del presente trabajo y de
experiencias previas en el uso de este modelo.
Coeficiente de escorrentía máximo (Cmax): Este coeficiente está relacionado
directamente con el tipo de respuesta de la cuenca ante la precipitación de entrada,
siendo la escorrentía inmediata lo primero que se satisface con el agua caída. El valor
del coeficiente depende de la geomorfología de la cuenca, en particular de la
pendiente de las laderas y del río, y del tipo y espesor de la capa de suelo superficial.
La tabla siguiente presenta algunos rangos recomendados de este parámetro.
Tabla 4.6: Rango valores Cmax según características geomorfológicas Muy baja Baja Media Alta Muy Alta(k < 10-6) (10-5 < k < 10-6) (10-3 < k < 10-5) (10-1 < k < 10-3) (k > 10--1)
Muy Alta (> 0.100 ) 0.80 - 0.90 0.75 - 0.80 0.70 - 0.80 0.70 - 0.75 0.60 - 0.70
Nota: Estos valores pueden ser menores incluso hasta en un 50% al incluir el módulo nival.
Permeabilidad del Suelo (k) [cm s-1]
Pend
. Lad
eras
[m m
-1]
Rango Cmax
Humedad máxima que puede retener el suelo (Hmax): Este parámetro está
relacionado con el espesor del estrato superficial y el tipo de suelo (permeabilidad). A
su vez influye en el cálculo de la evaporación y en el aporte de agua a la capa
subterránea a través de la percolación profunda. Una característica importante a
Capítulo 4: Aplicación a la Cuenca del Río Laja 33
considerar para su estimación es la capacidad de amortiguación de la cuenca y el
desfase de los caudales respecto de los eventos de precipitación. En términos
generales, este parámetro debe fluctuar entre 100 y 500 [mm], siendo el valor
máximo en cuencas con alta capacidad de amortiguación y el valor mínimo en
cuencas con baja capacidad de amortiguación. Un valor recomendado para iniciar la
calibración es 200 [mm], y luego se puede modificar en función del nivel de ajuste y
la respuesta de la cuenca.
D (porcentaje de precipitación sobre un límite que se transforma en
percolación profunda directa) y Plim (límite de precipitación sobre el cual existe
un aporte directo al embalse subterráneo): Estos dos parámetros se encuentran
interrelacionados y tienen influencia sobre la recarga del embalse subterráneo. Sirven
para mejorar el ajuste del flujo base y los valores recomendados para D son entre 0%
cuando no existe recarga y 40% cuando la capa superficial tiene una permeabilidad
media a muy alta y existe un flujo libre hasta el embalse subterráneo. Para Plim, se
estima un valor aceptable entre 3 y 8% de la precipitación media anual.
Porcentaje que define Hcrit (PORC): Este parámetro define un nivel de humedad
(a partir de Hmax), por sobre el cual no existe dificultad para la evapotranspiración. Es
importante destacar que no limita la existencia o no de evapotranspiración, sino que
sólo define determinados límites para la existencia de dificultad para dicho proceso.
Este valor tiende a ser bajo o nulo en cuencas de suelo poco permeable con alta
densidad de vegetación y con cuerpos de agua de gran extensión, ya que en dichos
casos la evapotranspiración real tiende a aproximarse a la potencial. En caso
contrario, en suelos permeables, donde el agua que no se transforma en escorrentía
directa pasa a la capa subsuperficial o subterránea, se tiene una mayor dificultad para
la evapotranspiración, y por lo tanto dicho valor tiende a ser mayor. La siguiente tabla
presenta algunas recomendaciones para el valor del PORC.
Tabla 4.7: Rango de valores de PORC según permeabilidad del suelo Muy baja Baja Media Alta Muy Alta(k < 10-6) (10-5 < k < 10-6) (10-3 < k < 10-5) (10-1 < k < 10-3) (k > 10--1)
PORC [%] 0 - 20 20 - 50 30 - 70 50 - 80 80 - 100
Nota: Estos valores pueden ser menores incluso hasta en un 50% al incluir el módulo nival.
Permeabilidad del Suelo (k) [cm s-1]
La Tabla 4.7 presenta valores recomendados del PORC sólo según el tipo de suelo.
Es importante considerar al momento de estimar este parámetro que cualquier otra
variable puede influir en su valor, como cuerpos de agua y densidad y tipo de
vegetación entre otros.
Capítulo 4: Aplicación a la Cuenca del Río Laja 34
Coeficiente de escorrentía subterránea (Ck): Este parámetro se relaciona con la
forma de descarga del embalse subterráneo. Se recomiendan valores bajos (0 a 0.2) en
cuencas donde el flujo base es estable a lo largo del año, valores medios (0.2 a 0.4) en
cuencas donde el flujo base presenta variabilidad interanual y valores altos (0.4 a 0.7)
en cuencas donde el flujo base presenta una alta variabilidad interanual y sigue de
manera notoria la forma del régimen de precipitaciones.
Parámetros de modificación de las variables de entrada (A y B): Estos parámetros
se deben considerar iguales a uno, a menos que se fundamente el modificar las series
de precipitación o evaporación potencial, en cuyo caso se debe justificar y calcular el
valor de cada parámetro en base fenómenos físicos y meteorológicos. Estos
parámetros, si bien pertenecen al modelo no deben ser utilizados como herramientas
de calibración.
4.4.2 Módulo Nival Este módulo requiere de cinco variables de entrada enfocadas a determinar la cota de la isoterma
cero y el área por sobre y bajo dicho nivel (hcg, gT, Tm y la curva hipsométrica, la cual es una
relación del área y la elevación de la cuenca), un valor de partida (el almacenamiento nival
inicial), cuatro parámetros del modelo destinados a estimar la tasa de derretimiento (M, Tb, DM
y F) y un parámetro de modificación de los inputs (FgT).
Cabe mencionar que con el módulo nival activado la precipitación caída se distribuye de manera
uniforme en dos cuencas de igual área, pero diferente régimen (pluvial y nival), luego, al ser
menor la precipitación pluvial, la tasa de escorrentía disminuye porque un mayor porcentaje de la
precipitación se infiltra y además, al ser menor la cantidad de agua caída sobre la cuenca pluvial,
menor es la saturación del suelo y mayor es la capacidad de retención de agua. Todo esto influye
en el coeficiente de escorrentía del módulo pluvial, el cual tiende a ser menor al considerar este
efecto.
Variables de Entrada: Las variables de entrada tienen como función caracterizar
aspectos de la morfología de la cuenca. Para ello se requiere de la temperatura media
mensual representativa de la cuenca (Tm). Esta temperatura se asume que actúa en el
centro de gravedad (hcg). Por otra parte, con el gradiente térmico (gT) se determinar
la elevación de la isoterma cero y con la curva hipsométrica (elevación versus área)
se determina el área bajo y sobre la isoterma cero y con esto se logra diferenciar entre
precipitación pluvial y nival y calcular la evolución del almacenamiento nival en el
Capítulo 4: Aplicación a la Cuenca del Río Laja 35
tiempo. Los datos de Tm se presentan en el Anexo B, los valores de hcg se presentan
en el punto 4.4.3 del presente capítulo, los datos del gradiente térmico y de la curva
hipsométrica de cada subcuenca se presentan en el Anexo E (Tabas E.1 a E.4).
Valores de Partida: Al igual que en el módulo pluvial, este módulo requiere para el
inicio de los cálculos el valor del almacenamiento nival inicial. La estimación de este
valor dependerá de la época de inicio de la simulación, pudiendo ser nulo si la
simulación comienza en verano.
Este valor se puede estimar o definir como nulo inicialmente, y luego de una primera
simulación se puede estimar en base al almacenamiento calculado el mismo mes del
año siguiente o subsiguiente. En la estimación de esta variable es importante
considerar las temperaturas medias anteriores al mes de inicio de la simulación, ya
que influyen en el derretimiento y por lo tanto en el almacenamiento nival de la
cuenca.
Parámetros de Calibración: Estos parámetros tienen como función determinar la
tasa de derretimiento del almacenamiento nival y la forma en que se transfiere a la
cuenca pluvial. A continuación se describen las funciones de cada parámetro.
Temperatura base (Tb): Corresponde a una temperatura de inicio del
derretimiento. Cualquier valor de Tm mayor que Tb produce un derretimiento
proporcional a la diferencia entre ambas. En general esta temperatura se define como
0 [°C], pero puede ser modificada en casos fundamentados.
Parámetro del modelo de derretimiento (M): Este parámetro define el
derretimiento potencial, el cual es una función de la diferencia entre Tm y Tb.
Define el volumen por unidad de área de nieve derretida por cada grado en que
Tm es mayor que Tb. Los valores este parámetro varían según la ubicación y
temperatura media de la cuenca, siendo mayor en cuencas más frías y su valor
normalmente debería fluctuar entre 1 y 12 [mm °C-1].
Coeficiente de derretimiento mínimo (M): Este parámetro define un porcentaje
del almacenamiento que se transforma en escorrentía y que sólo se produce cuando
Tm es menor que Tb. Su valor debe ser cercano a cero. En general, este parámetro
debe variar entre 0 y 0.1 lo que significa que entre el 0 y 10% del almacenamiento
nival se puede derretir cuando Tm es menor que Tb.
Capítulo 4: Aplicación a la Cuenca del Río Laja 36
Una forma de simplificar el módulo nival es asumir que no existe un
derretimiento cuando Tm es menor que Tb, en dicho caso basta con definir el valor de
M igual a cero.
Factor de transferencia del derretimiento a escorrentía inmediata (F): Este
factor define el porcentaje de nieve derretida que se transforma en EI. El
complemento (1-F) se infiltra aportando humedad al suelo. La escorrentía producida
por el derretimiento de nieve tiende a producir pequeños cauces que llegan
directamente al cauce principal o afluentes, por lo que este factor tiene a ser cercano a
uno (100%).
Parámetros de Modificación del Gradiente Térmico (FgT): Este parámetro
permite modificar el valor del gradiente térmico a la entrada del modelo. Por ejemplo,
en esta simulación es necesario, ya que el valor del gradiente térmico se subestima
debido a que su cálculo se realiza en base a las series sintéticas de temperatura, las
cuales presentan baja variabilidad entre cuadrantes.
4.4.3 Parámetros de Calibración y Validación Las variables, valores de partida y parámetros resultantes de la calibración y los coeficientes de
correlación obtenidos en la calibración y validación de cada subcuenca se resumen en las Tablas
4.8 a 4.10.
Tabla 4.8: Variables de partida y parámetros de calibración módulo pluvial Variables de
Entrada SC-1 SC-2* SC-3 SC-4 SC-5 SC-6 SC-7
PM y EM Anexo D Anexo D Anexo D Anexo D Anexo D Anexo D Anexo DAR [km2] 236.5 977.3 635.4 947.4 807.7 685.8** 215.6
Valores de Partida SC-1 SC-2 SC-3 SC-4 SC-5 SC-6 SC-7
FgT 1.85 1.85 1.10Nota: Las subcuencas restantes no presentan régimen nival.
MODULO NIVAL
Tabla 4.10: Coeficientes de correlación R2 obtenido entre caudales simulados y medidos (período Ago.2003 – Dic.2004 de calibración y Ene.2005 – Jun.2006 de validación).